SEARCH · Geology Reports

Results for “Geophysical Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30Linked to original sources

USGS2016

A synthesis of the basal thermal state of the Greenland Ice Sheet

The basal thermal state of an ice sheet (frozen or thawed) is an important control upon its evolution, dynamics and response to external forcings. However, this state can only be observed directly within sparse boreholes or inferred conclusively from the presence of subglacial lakes. Here we synthesize spatially extensive inferences of the basal thermal state of the Greenland Ice Sheet to better constrain this state. Existing inferences include outputs from the eight thermomechanical ice-flow models included in the SeaRISE effort. New remote-sensing inferences of the basal thermal state are derived from Holocene radiostratigraphy, modern surface velocity and MODIS imagery. Both thermomechanical modeling and remote inferences generally agree that the Northeast Greenland Ice Stream and large portions of the southwestern ice-drainage systems are thawed at the bed, whereas the bed beneath the central ice divides, particularly their west-facing slopes, is frozen. Elsewhere, there is poor agreement regarding the basal thermal state. Both models and remote inferences rarely represent the borehole-observed basal thermal state accurately near NorthGRIP and DYE-3. This synthesis identifies a large portion of the Greenland Ice Sheet (about one third by area) where additional observations would most improve knowledge of its overall basal thermal state.

Journal of Geophysical Research F: Earth Surface↗

USGS2015

Soil surface organic layers in Arctic Alaska: spatial distribution, rates of formation, and microclimatic effects

Organic layers of living and dead vegetation cover the ground surface in many permafrost landscapes and play important roles in ecosystem processes. These soil surface organic layers (SSOLs) store large amounts of carbon and buffer the underlying permafrost and its contained carbon from changes in aboveground climate. Understanding the dynamics of SSOLs is a prerequisite for predicting how permafrost and carbon stocks will respond to warming climate. Here we ask three questions about SSOLs in a representative area of the Arctic Foothills region of northern Alaska: (1) What environmental factors control the thickness of SSOLs and the carbon they store? (2) How long do SSOLs take to develop on newly stabilized point bars? (3) How do SSOLs affect temperature in the underlying ground? Results show that SSOL thickness and distribution correlate with elevation, drainage area, vegetation productivity, and incoming solar radiation. A multiple regression model based on these correlations can simulate spatial distribution of SSOLs and estimate the organic carbon stored there. SSOLs develop within a few decades after a new, sandy, geomorphic surface stabilizes but require 500–700 years to reach steady state thickness. Mature SSOLs lower the growing season temperature and mean annual temperature of the underlying mineral soil by 8 and 3°C, respectively. We suggest that the proximate effects of warming climate on permafrost landscapes now covered by SSOLs will occur indirectly via climate's effects on the frequency, extent, and severity of disturbances like fires and landslides that disrupt the SSOLs and interfere with their protection of the underlying permafrost.

Alaska↗

USGS2015

Riparian vegetation, Colorado River, and climate: five decades of spatiotemporal dynamics in the Grand Canyon with river regulation

Documentation of the interacting effects of river regulation and climate on riparian vegetation has typically been limited to small segments of rivers or focused on individual plant species. We examine spatiotemporal variability in riparian vegetation for the Colorado River in Grand Canyon relative to river regulation and climate, over the five decades since completion of the upstream Glen Canyon Dam in 1963. Long-term changes along this highly modified, large segment of the river provide insights for management of similar riparian ecosystems around the world. We analyze vegetation extent based on maps and imagery from eight dates between 1965 and 2009, coupled with the instantaneous hydrograph for the entire period. Analysis confirms a net increase in vegetated area since completion of the dam. Magnitude and timing of such vegetation changes are river stage-dependent. Vegetation expansion is coincident with inundation frequency changes and is unlikely to occur for time periods when inundation frequency exceeds approximately 5%. Vegetation expansion at lower zones of the riparian area is greater during the periods with lower peak and higher base flows, while vegetation at higher zones couples with precipitation patterns and decreases during drought. Short pulses of high flow, such as the controlled floods of the Colorado River in 1996, 2004, and 2008, do not keep vegetation from expanding onto bare sand habitat. Management intended to promote resilience of riparian vegetation must contend with communities that are sensitive to the interacting effects of altered flood regimes and water availability from river and precipitation.

Arizona↗

USGS2016

Effects of permafrost aggradation on peat properties as determined from a pan-Arctic synthesis of plant macrofossils

Permafrost dynamics play an important role in high-latitude peatland carbon balance and are key to understanding the future response of soil carbon stocks. Permafrost aggradation can control the magnitude of the carbon feedback in peatlands through effects on peat properties. We compiled peatland plant macrofossil records for the northern permafrost zone (515 cores from 280 sites) and classified samples by vegetation type and environmental class (fen, bog, tundra and boreal permafrost, and thawed permafrost). We examined differences in peat properties (bulk density, carbon (C), nitrogen (N) and organic matter content, and C/N ratio) and C accumulation rates among vegetation types and environmental classes. Consequences of permafrost aggradation differed between boreal and tundra biomes, including differences in vegetation composition, C/N ratios, and N content. The vegetation composition of tundra permafrost peatlands was similar to permafrost-free fens, while boreal permafrost peatlands more closely resembled permafrost-free bogs. Nitrogen content in boreal permafrost and thawed permafrost peatlands was significantly lower than in permafrost-free bogs despite similar vegetation types (0.9% versus 1.5% N). Median long-term C accumulation rates were higher in fens (23 g C m −2   yr −1 ) than in permafrost-free bogs (18 g C m −2  yr −1 ) and were lowest in boreal permafrost peatlands (14 g C m −2  yr −1 ). The plant macrofossil record demonstrated transitions from fens to bogs to permafrost peatlands, bogs to fens, permafrost aggradation within fens, and permafrost thaw and reaggradation. Using data synthesis, we have identified predominant peatland successional pathways, changes in vegetation type, peat properties, and C accumulation rates associated with permafrost aggradation.

Journal of Geophysical Research: Biogeosciences↗

USGS2016

Variation of energy and carbon fluxes from a restored temperate freshwater wetland and implications for carbon market verification protocols

Temperate freshwater wetlands are among the most productive terrestrial ecosystems, stimulating interest in using restored wetlands as biological carbon sequestration projects for greenhouse gas reduction programs. In this study, we used the eddy covariance technique to measure surface energy carbon fluxes from a constructed, impounded freshwater wetland during two annual periods that were 8 years apart: 2002–2003 and 2010–2011. During 2010–2011, we measured methane (CH 4 ) fluxes to quantify the annual atmospheric carbon mass balance and its concomitant influence on global warming potential (GWP). Peak growing season fluxes of latent heat and carbon dioxide (CO 2 ) were greater in 2002–2003 compared to 2010–2011. In 2002, the daily net ecosystem exchange reached as low as −10.6 g C m −2  d −1 , which was greater than 3 times the magnitude observed in 2010 (−2.9 g C m −2  d −1 ). CH 4 fluxes during 2010–2011 were positive throughout the year and followed a strong seasonal pattern, ranging from 38.1 mg C m −2  d −1 in the winter to 375.9 mg C m −2  d −1 during the summer. The results of this study suggest that the wetland had reduced gross ecosystem productivity in 2010–2011, likely due to the increase in dead plant biomass (standing litter) that inhibited the generation of new vegetation growth. In 2010–2011, there was a net positive GWP (675.3 g C m −2  yr −1 ), and when these values are evaluated as a sustained flux, the wetland will not reach radiative balance even after 500 years.

California↗

USGS2015

Incorporating microbial dormancy dynamics into soil decomposition models to improve quantification of soil carbon dynamics of northern temperate forests

Soil carbon dynamics of terrestrial ecosystems play a significant role in the global carbon cycle. Microbial-based decomposition models have seen much growth recently for quantifying this role, yet dormancy as a common strategy used by microorganisms has not usually been represented and tested in these models against field observations. Here we developed an explicit microbial-enzyme decomposition model and examined model performance with and without representation of microbial dormancy at six temperate forest sites of different forest types. We then extrapolated the model to global temperate forest ecosystems to investigate biogeochemical controls on soil heterotrophic respiration and microbial dormancy dynamics at different temporal-spatial scales. The dormancy model consistently produced better match with field-observed heterotrophic soil CO 2 efflux ( R H ) than the no dormancy model. Our regional modeling results further indicated that models with dormancy were able to produce more realistic magnitude of microbial biomass (<2% of soil organic carbon) and soil R H (7.5 ± 2.4 Pg C yr −1 ). Spatial correlation analysis showed that soil organic carbon content was the dominating factor (correlation coefficient = 0.4–0.6) in the simulated spatial pattern of soil R H with both models. In contrast to strong temporal and local controls of soil temperature and moisture on microbial dormancy, our modeling results showed that soil carbon-to-nitrogen ratio (C:N) was a major regulating factor at regional scales (correlation coefficient = −0.43 to −0.58), indicating scale-dependent biogeochemical controls on microbial dynamics. Our findings suggest that incorporating microbial dormancy could improve the realism of microbial-based decomposition models and enhance the integration of soil experiments and mechanistically based modeling.

Journal of Geophysical Research G: Biogeosciences↗

USGS2015

Abundance, stable isotopic composition, and export fluxes of DOC, POC, and DIC from the Lower Mississippi River during 2006–2008

Sources, abundance, isotopic compositions, and export fluxes of dissolved inorganic carbon (DIC), dissolved and colloidal organic carbon (DOC and COC), and particulate organic carbon (POC), and their response to hydrologic regimes were examined through monthly sampling from the Lower Mississippi River during 2006–2008. DIC was the most abundant carbon species, followed by POC and DOC. Concentration and δ 13 C of DIC decreased with increasing river discharge, while those of DOC remained fairly stable. COC comprised 61 ± 3% of the bulk DOC with similar δ 13 C abundances but higher percentages of hydrophobic organic acids than DOC, suggesting its aromatic and diagenetically younger status. POC showed peak concentrations during medium flooding events and at the rising limb of large flooding events. While δ 13 C-POC increased, δ 15 N of particulate nitrogen decreased with increasing discharge. Overall, the differences in δ 13 C between DOC or DIC and POC show an inverse correlation with river discharge. The higher input of soil organic matter and respired CO 2 during wet seasons was likely the main driver for the convergence of δ 13 C between DIC and DOC or POC, whereas enhanced in situ primary production and respiration during dry seasons might be responsible for their isotopic divergence. Carbon export fluxes from the Mississippi River were estimated to be 13.6 Tg C yr −1 for DIC, 1.88 Tg C yr −1 for DOC, and 2.30 Tg C yr −1 for POC during 2006–2008. The discharge-normalized DIC yield decreased during wet seasons, while those of POC and DOC increased and remained constant, respectively, implying variable responses in carbon export to the increasing discharge.

Louisiana↗

USGS2016

Wetland tree transpiration modified by river-floodplain connectivity

Hydrologic connectivity provisions water and nutrient subsidies to floodplain wetlands and may be particularly important in floodplains with seasonal water deficits through its effects on soil moisture. In this study, we measured sapflow in 26 trees of two dominant floodplain forest species ( Celtis laevigata and Quercus lyrata ) at two hydrologically distinct sites in the lower White River floodplain in Arkansas, USA. Our objective was to investigate how connectivity-driven water table variations affected water use, an indicator of tree function. Meteorological variables (photosynthetically active radiation and vapor pressure deficit) were the dominant controls over water use at both sites; however, water table variations explained some site differences. At the wetter site, highest sapflow rates were during a late-season overbank flooding event, and no flood stress was apparent. At the drier site, sapflow decreased as the water table receded. The late-season flood pulse that resulted in flooding at the wetter site did not affect the water table at the drier site; accordingly, higher water use was not observed at the drier site. The species generally associated with wetter conditions ( Q. lyrata ) was more positively responsive to the flood pulse. Flood water subsidy lengthened the effective growing season, demonstrating ecological implications of hydrologic connectivity for alleviating water deficits that otherwise reduce function in this humid floodplain wetland.

Arkansas↗

USGS2016

Component greenhouse gas fluxes and radiative balance from two deltaic marshes in Louisiana: Pairing chamber techniques and eddy covariance

Coastal marshes take up atmospheric CO 2 while emitting CO 2 , CH 4 , and N 2 O. This ability to sequester carbon (C) is much greater for wetlands on a per-area basis than from most ecosystems, facilitating scientific, political, and economic interest in their value as greenhouse gas sinks. However, the greenhouse gas balance of Gulf of Mexico wetlands is particularly understudied. We describe the net ecosystem exchange (NEE c ) of CO 2 and CH 4 using eddy covariance (EC) in comparison with fluxes of CO 2 , CH 4 , and N 2 O using chambers from brackish and freshwater marshes in Louisiana, USA. From EC, we found that 182 g C m -2 y -1 was lost through NEE c from the brackish marsh. Of this, 11 g C m -2 y -1 resulted from net CH 4 emissions and the remaining 171 g C m -2 y -1 resulted from net CO 2 emissions. In contrast, -290 g C m 2 y -1 was taken up through NEE c by the freshwater marsh, with 47 g C m -2 y -1 emitted as CH 4 and -337 g C m -2 y -1 taken up as CO 2 . From chambers, we discovered that neither site had large fluxes of N 2 O. Sustained-flux greenhouse gas accounting metrics indicated that both marshes had a positive (warming) radiative balance, with the brackish marsh having a substantially greater warming effect than the freshwater marsh. That net respiratory emissions of CO 2 and CH 4 as estimated through chamber techniques were 2-4 times different from emissions estimated through EC requires additional understanding of the artifacts created by different spatial and temporal sampling footprints between techniques.

Journal of Geophysical Research: Biogeosciences↗

USGS2016

The differing biogeochemical and microbial signatures of glaciers and rock glaciers

Glaciers and rock glaciers supply water and bioavailable nutrients to headwater mountain lakes and streams across all regions of the American West. Here we present a comparative study of the metal, nutrient, and microbial characteristics of glacial and rock glacial influence on headwater ecosystems in three mountain ranges of the contiguous U.S.: The Cascade Mountains, Rocky Mountains, and Sierra Nevada. Several meltwater characteristics (water temperature, conductivity, pH, heavy metals, nutrients, complexity of dissolved organic matter (DOM), and bacterial richness and diversity) differed significantly between glacier and rock glacier meltwaters, while other characteristics (Ca 2+ , Fe 3+ , SiO 2 concentrations, reactive nitrogen, and microbial processing of DOM) showed distinct trends between mountain ranges regardless of meltwater source. Some characteristics were affected both by glacier type and mountain range (e.g. temperature, ammonium (NH 4 + ) and nitrate (NO 3 - ) concentrations, bacterial diversity). Due to the ubiquity of rock glaciers and the accelerating loss of the low latitude glaciers our results point to the important and changing influence that these frozen features place on headwater ecosystems.

California, Colorado, Oregon, Washington, Wyoming↗

USGS2016

Climate change and dissolved organic carbon export to the Gulf of Maine

Ongoing climate change is affecting the concentration, export (flux), and timing of dissolved organic carbon (DOC) exported to the Gulf of Maine (GoM) through changes in hydrologic regime. DOC export was calculated for water years 1950 through 2013 for 20 rivers and for water years 1930 through 2013 for 14 rivers draining to the GoM. DOC export was also estimated for the 21st century based on climate and hydrologic modeling in a previously published study. DOC export was calculated by using the regression model LOADEST to fit seasonally adjusted concentration discharge (C-Q) relations. Our results are an analysis of the sensitivity of DOC export to changes in hydrologic conditions over time since land cover and vegetation were held constant over time. Despite large interannual variability, all rivers had increasing DOC export during winter and these trends were significant ( p < 0.05) in 10 out of 20 rivers for 1950 to 2013 and in 13 out of 14 rivers for 1930 to 2013. All rivers also had increasing annual export of DOC although fewer trends were statistically significant than for winter export. Projections for DOC export during the 21st century were variable depending on the climate model and greenhouse gas emission scenario that affected future river discharge through effects on precipitation and evapotranspiration. The most consistent result was a significant increase in DOC export in winter in all model-by-emission scenarios. DOC export was projected to decrease during the summer in all model-by-emission scenarios, with statistically significant decreases in half of the scenarios.

Maine↗

USGS2016

Induced earthquake magnitudes are as large as (statistically) expected

A major question for the hazard posed by injection-induced seismicity is how large induced earthquakes can be. Are their maximum magnitudes determined by injection parameters or by tectonics? Deterministic limits on induced earthquake magnitudes have been proposed based on the size of the reservoir or the volume of fluid injected. However, if induced earthquakes occur on tectonic faults oriented favorably with respect to the tectonic stress field, then they may be limited only by the regional tectonics and connectivity of the fault network. In this study, we show that the largest magnitudes observed at fluid injection sites are consistent with the sampling statistics of the Gutenberg-Richter distribution for tectonic earthquakes, assuming no upper magnitude bound. The data pass three specific tests: (1) the largest observed earthquake at each site scales with the log of the total number of induced earthquakes, (2) the order of occurrence of the largest event is random within the induced sequence, and (3) the injected volume controls the total number of earthquakes rather than the total seismic moment. All three tests point to an injection control on earthquake nucleation but a tectonic control on earthquake magnitude. Given that the largest observed earthquakes are exactly as large as expected from the sampling statistics, we should not conclude that these are the largest earthquakes possible. Instead, the results imply that induced earthquake magnitudes should be treated with the same maximum magnitude bound that is currently used to treat seismic hazard from tectonic earthquakes.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Subsidence rates at the southern Salton Sea consistent with reservoir depletion

Space geodetic measurements from the Envisat satellite between 2003 and 2010 show that subsidence rates near the southeastern shoreline of the Salton Sea in Southern California are up to 52mmyr −1 greater than the far-field background rate. By comparing these measurements with model predictions, we find that this subsidence appears to be dominated by poroelastic contraction associated with ongoing geothermal fluid production, rather than the purely fault-related subsidence proposed previously. Using a simple point source model, we suggest that the source of this proposed volumetric strain is at depths between 1.0 km and 2.4 km (95% confidence interval), comparable to generalized boundaries of the Salton Sea geothermal reservoir. We find that fault slip on two previously imaged tectonic structures, which are part of a larger system of faults in the Brawley Seismic Zone, is not an adequate predictor of surface velocity fields because the magnitudes of the best fitting slip rates are often greater than the full plate boundary rate and at least 2 times greater than characteristic sedimentation rates in this region. Large-scale residual velocity anomalies indicate that spatial patterns predicted by fault slip are incompatible with the observations.

California↗

USGS2016

Alternative source models of very low frequency events

We present alternative source models for very low frequency (VLF) events, previously inferred to be radiation from individual slow earthquakes that partly fill the period range between slow slip events lasting thousands of seconds and low-frequency earthquakes (LFE) with durations of tenths of a second. We show that VLF events may emerge from bandpass filtering a sum of clustered, shorter duration, LFE signals, believed to be the components of tectonic tremor. Most published studies show VLF events occurring concurrently with tremor bursts and LFE signals. Our analysis of continuous data from Costa Rica detected VLF events only when tremor was also occurring, which was only 7% of the total time examined. Using analytic and synthetic models, we show that a cluster of LFE signals produces the distinguishing characteristics of VLF events, which may be determined by the cluster envelope. The envelope may be diagnostic of a single, dynamic, slowly slipping event that propagates coherently over kilometers or represents a narrowly band-passed version of nearly simultaneous arrivals of radiation from slip on multiple higher stress drop and/or faster propagating slip patches with dimensions of tens of meters (i.e., LFE sources). Temporally clustered LFE sources may be triggered by single or multiple distinct aseismic slip events or represent the nearly simultaneous chance occurrence of background LFEs. Given the nonuniqueness in possible source durations, we suggest it is premature to draw conclusions about VLF event sources or how they scale.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

A large refined catalog of earthquake relocations and focal mechanisms for the Island of Hawai'i and its seismotectonic implications

We present high-quality focal mechanisms based on a refined earthquake location catalog for the Island of Hawai'i, focusing on Mauna Loa and Kīlauea volcanoes. The relocation catalog is based on first-arrival times and waveform data of both compressional and shear waves for about 180,000 events on and near the Island of Hawai'i between 1986 and 2009 recorded by the seismic stations at the Hawaiian Volcano Observatory. We relocate all the earthquakes by applying ray tracing through an existing three-dimensional velocity model, similar event cluster analysis, and a differential-time relocation method. The resulting location catalog represents an expansion of previous relocation studies, covering a longer time period and consisting of more events with well-constrained absolute locations. The focal mechanisms are obtained based on the compressional-wave first-motion polarities and compressional-to-shear wave amplitude ratios by applying the HASH program to the waveform cross correlation relocated earthquakes. Overall, the good-quality (defined by the HASH parameters) focal solutions are dominated by normal faulting in our study area, especially in the active Ka'ōiki and Hīlea seismic zones. Kīlauea caldera is characterized by a mixture of approximately equal numbers of normal, strike-slip, and reverse faults, whereas its south flank has slightly fewer strike-slip events. Our relocation and focal mechanism results will be useful for mapping the seismic stress and strain fields and for understanding the seismic-volcanic-tectonic relationships within the magmatic systems.

Hawai'i↗

USGS2016

Volcano dome dynamics at Mount St. Helens: Deformation and intermittent subsidence monitored by seismicity and camera imagery pixel offsets

The surface deformation field measured at volcanic domes provides insights into the effects of magmatic processes, gravity- and gas-driven processes, and the development and distribution of internal dome structures. Here we study short-term dome deformation associated with earthquakes at Mount St. Helens, recorded by a permanent optical camera and seismic monitoring network. We use Digital Image Correlation (DIC) to compute the displacement field between successive images and compare the results to the occurrence and characteristics of seismic events during a 6 week period of dome growth in 2006. The results reveal that dome growth at Mount St. Helens was repeatedly interrupted by short-term meter-scale downward displacements at the dome surface, which were associated in time with low-frequency, large-magnitude seismic events followed by a tremor-like signal. The tremor was only recorded by the seismic stations closest to the dome. We find a correlation between the magnitudes of the camera-derived displacements and the spectral amplitudes of the associated tremor. We use the DIC results from two cameras and a high-resolution topographic model to derive full 3-D displacement maps, which reveals internal dome structures and the effect of the seismic activity on daily surface velocities. We postulate that the tremor is recording the gravity-driven response of the upper dome due to mechanical collapse or depressurization and fault-controlled slumping. Our results highlight the different scales and structural expressions during growth and disintegration of lava domes and the relationships between seismic and deformation signals.

Washington↗

USGS2016

Insights into shallow magmatic processes at Kīlauea Volcano, Hawaiʻi, from a multiyear continuous gravity time series

Continuous gravity data collected near the summit eruptive vent at Kīlauea Volcano, Hawaiʻi, during 2011–2015 show a strong correlation with summit-area surface deformation and the level of the lava lake within the vent over periods of days to weeks, suggesting that changes in gravity reflect variations in volcanic activity. Joint analysis of gravity and lava level time series data indicates that over the entire time period studied, the average density of the lava within the upper tens to hundreds of meters of the summit eruptive vent remained low—approximately 1000–1500 kg/m 3 . The ratio of gravity change (adjusted for Earth tides and instrumental drift) to lava level change measured over 15 day windows rose gradually over the course of 2011–2015, probably reflecting either (1) a small increase in the density of lava within the eruptive vent or (2) an increase in the volume of lava within the vent due to gradual vent enlargement. Superimposed on the overall time series were transient spikes of mass change associated with inflation and deflation of Kīlauea's summit and coincident changes in lava level. The unexpectedly strong mass variations during these episodes suggest magma flux to and from the shallow magmatic system without commensurate deformation, perhaps indicating magma accumulation within, and withdrawal from, void space—a process that might not otherwise be apparent from lava level and deformation data alone. Continuous gravity data thus provide unique insights into magmatic processes, arguing for continued application of the method at other frequently active volcanoes.

Hawaii↗

USGS2016

Coseismic slip and early afterslip of the 2015 Illapel, Chile, earthquake: Implications for frictional heterogeneity and coastal uplift

Great subduction earthquakes are thought to rupture portions of the megathrust, where interseismic coupling is high and velocity-weakening frictional behavior is dominant, releasing elastic deformation accrued over a seismic cycle. Conversely, postseismic afterslip is assumed to occur primarily in regions of velocity-strengthening frictional characteristics that may correlate with lower interseismic coupling. However, it remains unclear if fixed frictional properties of the subduction interface, coseismic or aftershock-induced stress redistribution, or other factors control the spatial distribution of afterslip. Here we use interferometric synthetic aperture radar and Global Position System observations to map the distribution of coseismic slip of the 2015 M w 8.3 Illapel, Chile, earthquake and afterslip within the first 38 days following the earthquake. We find that afterslip overlaps the coseismic slip area and propagates along-strike into regions of both high and moderate interseismic coupling. The significance of these observations, however, is tempered by the limited resolution of geodetic inversions for both slip and coupling. Additional afterslip imaged deeper on the fault surface bounds a discrete region of deep coseismic slip, and both contribute to net uplift of the Chilean Coastal Cordillera. A simple partitioning of the subduction interface into regions of fixed frictional properties cannot reconcile our geodetic observations. Instead, stress heterogeneities, either preexisting or induced by the earthquake, likely provide the primary control on the afterslip distribution for this subduction zone earthquake. We also explore the occurrence of coseismic and postseismic coastal uplift in this sequence and its implications for recent hypotheses concerning the source of permanent coastal uplift along subduction zones.

Journal of Geophysical Research B: Solid Earth↗