USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Journal”

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

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

At least 1,333 records · Page 74Linked to original sources
USGS2018

Environmental controls, emergent scaling, and predictions of greenhouse gas (GHG) fluxes in coastal salt marshes

Coastal salt marshes play an important role in mitigating global warming by removing atmospheric carbon at a high rate. We investigated the environmental controls and emergent scaling of major greenhouse gas (GHG) fluxes such as carbon dioxide (CO 2 ) and methane (CH 4 ) in coastal salt marshes by conducting data analytics and empirical modeling. The underlying hypothesis is that the salt marsh GHG fluxes follow emergent scaling relationships with their environmental drivers, leading to parsimonious predictive models. CO 2 and CH 4 fluxes, photosynthetically active radiation (PAR), air and soil temperatures, well water level, soil moisture, and porewater pH and salinity were measured during May–October 2013 from four marshes in Waquoit Bay and adjacent estuaries, MA, USA. The salt marshes exhibited high CO 2 uptake and low CH 4 emission, which did not significantly vary with the nitrogen loading gradient (5–126 kg · ha −1 · year −1 ) among the salt marshes. Soil temperature was the strongest driver of both fluxes, representing 2 and 4–5 times higher influence than PAR and salinity, respectively. Well water level, soil moisture, and pH did not have a predictive control on the GHG fluxes, although both fluxes were significantly higher during high tides than low tides. The results were leveraged to develop emergent power law‐based parsimonious scaling models to accurately predict the salt marsh GHG fluxes from PAR, soil temperature, and salinity (Nash‐Sutcliffe Efficiency = 0.80–0.91). The scaling models are available as a user‐friendly Excel spreadsheet named Coastal Wetland GHG Model to explore scenarios of GHG fluxes in tidal marshes under a changing climate and environment.

Massachusetts↗
USGS2018

Estimating soil respiration in a subalpine landscape using point, terrain, climate and greenness data

Landscape carbon (C) flux estimates are necessary for assessing the ability of terrestrial ecosystems to buffer further increases in anthropogenic carbon dioxide (CO2) emissions. Advances in remote sensing have allowed for coarse-scale estimates of gross primary productivity (GPP) (e.g., MODIS 17), yet efforts to assess spatial patterns in respiration lag behind those of GPP. Here, we demonstrate a method to predict growing season soil respiration at a regional scale in a forested ecosystem. We related field measurements (n=144) of growing season soil respiration across subalpine forests in the Southern Rocky Mountains ecoregion to a suite of biophysical predictors with a Random Forest model (30 m pixel size). We found that Landsat Enhanced Vegetation Index (EVI), growing season AI, temperature, precipitation, elevation, and slope aspect explained spatiotemporal variability in soil respiration. Our model had a psuedo-r2 of 0.45 and root mean squared error (RMSE) of roughly one-quarter of the mean value of respiration. Predicted growing season soil respiration across the region was remarkably consistent across 2004, 2005 and 2006 (150-d averages of 542.8, 544.3, and 536.5 g C m-2, respectively). Yet, we observed substantial variability in spatial patterns of soil respiration predictions that varied between years, suggesting that our method is sensitive to changes in respiration drivers. We compared our estimates to MODIS GPP and nocturnal net ecosystem exchange (NEE) derived from eddy covariance towers as a proxy for ecosystem respiration. Averaged across the predictive region, mean predicted growing season soil respiration was 73% of MODIS GPP, while predicted soil respiration was generally within 20% of nocturnal NEE from eddy covariance towers. This study demonstrated that geospatial and remotely-sensed datasets can be used in a statistical modeling framework to estimate soil respiration at landscape scales.

Colorado, Wyoming↗
USGS2019

Seasonal and spatial variation in the location and reactivity of a nitrate-contaminated groundwater discharge zone in a lakebed

Groundwater discharge delivering anthropogenic N from surrounding watersheds can impact lake nutrient budgets. However, upgradient groundwater processes and changing dynamics in N biogeochemistry at the groundwater-lake interface are complex and difficult to resolve. In this study, hydrograph variations in a groundwater flow-through lake altered discharge patterns of a wastewater-derived, groundwater contaminant plume, thereby affecting biogeochemical processes controlling N transport. Groundwater geochemistry 15 cm under the lakebed along transects perpendicular to shore varied from oxic to anoxic with increasing nitrate concentrations (10-75 M) and corresponding gradients in nitrite and nitrous oxide. Porewater depth profiles of nitrate concentrations and stable isotope compositions largely reflected upgradient groundwater N sources and N-cycle processes, with minor additional nitrate reduction in the shallowest lakebed sediments. Potential denitrification rates determined in laboratory microcosms were 10-100 fold higher in near-surface sediments (0-5 cm) than in deeper sediments (5-30 cm) and were correlated with sediment carbon content and abundance of denitrification genes (nirS, nosZI, and nosZII). Potential anammox-driven N2 production was highest in deeper anoxic sediments. Injection of bromide and nitrite in the lake sediments indicated a vertical porewater velocity of 4-5 cm hr-1, with highest nitrite consumption rates above 10 cm. However, short residence times in the shallow sediments allowed only a small fraction of the contaminant nitrate to be removed before discharging into the lake. Results demonstrate the importance of resolving local versus upgradient biogeochemical processes affecting contaminant distribution in discharge areas, and transient migration of local gradients and processes in response to changing lake levels and groundwater flow paths.

Journal of Geophysical Research: Biogeosciences↗
USGS2018

El Niño increases high‐tide flooding in tidal wetlands along the U.S. Pacific coast.

Periodic oscillations between El Niño and La Niña conditions in the Pacific Basin affect oceanographic and meteorological phenomena globally, with impacts on the abundance and distribution of marine species. However, El Niño effects on estuarine hydrology and tidal wetland processes have seldom been examined rigorously. We used detailed wetland elevation and local inundation data from 10 tidal wetlands located along the Pacific coast of the United States to assess changes in flooding during the 2015–2016 El Niño and to determine decadal‐scale relationships between estuarine sea‐level anomalies and Pacific Basin climate indices for this region. During the 2015–2016 El Niño all sites experienced significant increases in high‐tide water levels exceeding those predicted by astronomical tides, and increased flooding frequency during at least one of the El Niño subperiods relative to pre‐El Niño conditions. The magnitude of positive sea‐level anomalies varied by site (4–15 cm), with local hot spots of high water in southern Oregon, northern California, and Pt. Mugu lagoon in the Southern California Bight. Furthermore, over the last three decades of historic tide records, there were positive relationships between high‐tide sea‐level anomalies and equatorial Pacific Basin sea surface temperature anomalies across the region, and negative relationships with the Northern Oscillation Index. Increases of 1 °C in equatorial sea surface temperature were associated with 3–5 cm of increased high‐tide flooding at the sites. Elevated estuarine flooding associated with future El Niños could impact important tidal wetland processes and could be an additive stressor for wetlands facing accelerating sea‐level rise.

California, Oregon, Washington↗
USGS2020

Persistence and plasticity in conifer water-use strategies

The selective use of seasonal precipitation by vegetation is critical to understanding the residence time and flow path of water in watersheds, yet there are limited datasets to test how climate alters these dynamics. Here, we use measurements of the seasonal cycle of tree ring 18O for two widespread conifer species in the Rocky Mountains of North America to provide a multi-decadal depiction of the seasonal origins of forest water use. The results show that while the conifer tree stands had a dominant preference for use of snowmelt, there were multi-annual periods over the last four decades when use of summer precipitation was preferential. Utilization of summer rain emerged during years with increased snowfall and tree growth, suggesting that summer rain enhanced the transpiration stream only during the periods of highest water use. We hypothesize this could be explained through shallowing of the root profile during wetter periods and/or through the influence of changing water table depths on the residence time of summer precipitation in the root zone. We suggest the tree ring proxy approach used here could be applied in other watersheds to provide critical insight into the temporal dynamics of plant water use that could not be inferred from short measurement campaigns. These data on the seasonal origins of forest water are critical for understanding forest vulnerability to drought, the processes that affect precipitation pathways and residence time in watersheds and the interpretation of tree ring proxy data.

Colorado↗
USGS2019

Assessing the chemistry and bioavailability of dissolved organic matter from glaciers and rock glaciers

As glaciers thaw in response to warming, they release dissolved organic matter (DOM) to alpine lakes and streams. The United States contains an abundance of both alpine glaciers and rock glaciers. Differences in DOM composition and bioavailability between glacier types, like rock and ice glaciers, remain undefined. To assess differences in glacier and rock glacier DOM we evaluated bioavailability and molecular composition of DOM from four alpine catchments each with a glacier and a rock glacier at their headwaters. We assessed bioavailability of DOM by incubating each DOM source with a common microbial community and evaluated chemical characteristics of DOM before and after incubation using untargeted gas chromatography–mass spectrometry‐based metabolomics. Prior to incubations, ice glacier and rock glacier DOM had similar C:N ratios and chemical diversity, but differences in DOM composition. Incubations with a common microbial community showed that DOM from ice glacier meltwaters contained a higher proportion of bioavailable DOM and resulted in greater bacterial growth efficiency. After incubation, DOM composition from each source was statistically indistinguishable. This study provides an example of how MS‐based metabolomics can be used to assess effects of DOM composition on differences in bioavailability of DOM. Furthermore, it illustrates the importance of microbial metabolism in structuring composition of DOM. Even though rock glaciers had significantly less bioavailable DOM than ice glaciers, both glacial types still have potential to be important sources of bioavailable DOM to alpine headwaters over the coming decades.

Colorado↗
USGS2019

Identifying the molecular signatures of agricultural expansion in Amazonian headwater streams

Agricultural impacts on aquatic ecosystems are well-studied, however, most research has focused on temperate regions, whereas the forefront of agricultural expansion is currently in the tropics. At the vanguard of this growth is the boundary between the Amazon and Cerrado biomes in Brazil, driven primarily by expansion of soybean and corn croplands. Here we examine the impacts of cropland expansion on receiving lowland Amazon Basin headwater streams in terms of dissolved organic carbon (DOC) concentration and dissolved organic matter (DOM) composition via ultra-high resolution mass spectrometry. Streams draining croplands had lower DOC concentrations and DOM molecular signatures enriched in N- and S-containing formula in comparison to forested streams. Cropland streams were also enriched in aliphatic, peptide-like and highly unsaturated and phenolic (low O/C) compound categories in comparison to forest streams (enriched in polyphenolics, condensed aromatics and highly unsaturated and phenolic (high O/C) compound categories) indicative of the shifting of sources from organic-rich surface soils and litter layers to autochthonous and more microbial biomass. Distinct molecular assemblages were strongly correlated with cropland and forest catchments, highlighting headwater streams as sentinels for detecting change. On investigation of unique molecular formulae present in only cropland sites, four cropland markers provided the ability to track agricultural impacts in the region. Overall, these patterns indicate reduced organic matter inputs in croplands and greater microbial degradation at these sites leading to declining DOC concentrations, and DOM of more microbial character in receiving streams that is more biolabile, with clear ramifications for downstream ecology and biogeochemical cycles.

Amazon River basin, Tanguro Ranch, Xinga River↗
USGS2020

Modeling soil porewater salinity response to drought in tidal freshwater forested wetlands

There is a growing concern about the adverse effects of saltwater intrusion via tidal rivers, streams and creeks into tidal freshwater forested wetlands (TFFW) due to sea‐level rise (SLR) and intense and extended drought events. However, the magnitude and duration of porewater salinity in exceedance of plant salinity stress threshold (2 practical salinity units, psu) and the controlling factors remain unclear. In this study, we developed a TFFW soil porewater salinity model, in which the feedback mechanisms between soil salinity and evapotranspiration and hydraulic conductivity were incorporated. We selected sites (upper, middle, lower tidal freshwater forest sites and oligohaline marsh site) along the coastal floodplains of two rivers, the Waccamaw River (SC, USA) and the Savannah River (GA and SC, USA), that represent landscape salinity gradients from tidal influence of the Atlantic Ocean. The model results agreed well with field measurements and revealed that with drought‐induced saltwater intrusion, the mean annual soil porewater salinity and duration of elevated soil porewater salinity (> 2 psu) increased significantly compared to the normal (non‐drought) condition, posing a threat to the health and ecosystem services of TFFW even in the absence of SLR. Model results also showed more severe salinity stress under drought for the lower forest sites along the two rivers, where soil salinity values have already been at or in exceedance of the 2 psu threshold.

Georgia, South Carolina↗
USGS2019

High rates of inflation during a noneruptive episode of seismic unrest at Semisopochnoi Volcano, Alaska in 2014–2015

Magma intrusion rate is a key parameter in eruption triggering but is poorly quantified in existing geodetic studies. Here we examine two episodes of rapid inflation in this context. Two noneruptive microseismic swarms were recorded at Semisopochnoi Volcano, Alaska in 2014–2015. We use differential SAR techniques and TerraSAR‐X images to document surface deformation from 2011 to 2015, which comprises island‐wide radial inflation totaling ~25 cm (+/−1 cm) line of sight displacement in 2014–2015. Multiple source geometries are tested in an inversion of the deformation data, and InSAR data are best fit by a spheroid trending to the northeast and plunging to the southeast, with a major axis of ~4 km and minor axes of ~1 km, directly under the central caldera of Semisopochnoi. In 2014, a modeled influx of 0.043 km 3 of magma caused line of sight displacement of ~17 cm. This magma was stored at a depth of ~8 km, until 2015 when 0.029 km 3 was added. Along with the definition of inflation source parameters, the recorded seismic events are relocated using differential travel times. These relocated events outline a linear aseismic area within a larger zone of shallow (<10 km) seismicity. This aseismic region aligns with the centroid of the deformation model. Based on these geodetic and seismic models, the plumbing system at Semisopochnoi is interpreted as a spheroidal magma storage zone at a depth of ˜8 km below a linear feature of partial melt. The observed deformation and seismicity appear to result from rapid injection into this main storage region.

Alaska↗
USGS2019

Leakage and increasing fluid pressure detected in Oklahoma's wastewater disposal reservoir

The Arbuckle Group is the principal reservoir used for wastewater disposal in Oklahoma. In Osage County—a seismically quiet part of the state—continuous measurements of fluid pressure reveal that pressure in the reservoir is increasing by at least 5 kPa annually and sometimes at a much higher rate. Tidal analysis reveals that fluid level changes lead the local strain tides, with no apparent influence from transient permeability changes; this indicates a response that is inconsistent with flow in a radially extensive, confined reservoir. We investigate whether this is due to vertical flow to the water table, vertical flow within the Arbuckle, or local distortions from fractures. While none of these alternative models can fully explain both the observed tidal phases and amplitude ratios, the observed response to teleseismic waves supports a mechanism related to leakage rather than fracture effects. At this location fluid influx associated with wastewater disposal is offset by migration into surrounding layers, which include the Precambrian basement below. Thus, our findings suggest the need to monitor for changes in the induced seismicity hazard, while pore pressures increase in a leaky disposal reservoir.

Oklahoma↗
USGS2019

Interseismic quiescence and triggered slip of active normal faults of Kīlauea Volcano’s south flank during 2001-2018

The mobile south flank of Kīlauea Volcano hosts two normal fault systems, the Koa'e fault system (KFS) and the Hilina fault system (HFS). In historical time, at least three M>6.5 earthquakes have occurred on the basal detachment of the Kīlauea Volcano's south flank, with the most recent being the 4 May 2018 M6.9 earthquake. Here we analyze kinematic Global Positioning System data collected from 2001 to 2017 and interferometric synthetic aperture radar data before, during, and after the 2018 M6.9 earthquake to determine the crustal motion across the HFS and KFS faults. Our results indicate that the HFS faults did not significantly slip during the interseismic period from 2007 to 2011. Despite its substantial magnitude, interferometric synthetic aperture radar (InSAR) data show that the 2018 M6.9 earthquake triggered subcentimeter level slip along sections of the previously mapped HFS branches. Up to 20 cm of offset occurred on what appears to be a newly formed (or previously unknown) fault near the eastern end of the HFS. During the 3 months following the M6.9 earthquake, up to ~30 cm of slip occurred along the KFS, which helps accommodate rapid large‐scale subsidence of Kīlauea's summit region as large volumes of summit reservoir magma fed the lower East Rift Zone eruption. The HFS appears to activate only in concert with large earthquakes on the basal detachment. The KFS, on the other hand, moves both seismically during small local earthquakes and aseismically in response to nearby earthquakes and caldera subsidence.

Hawaii↗
USGS2019

Insights into pāhoehoe lava emplacement using visible and thermal structure-from-motion photogrammetry

We present the evolution over 3 months of a 2016–2017 pāhoehoe flow at Kīlauea as it changed from a narrow sheet flow into a compound lava field fed by a stable system of tubes. The portion of the flow located on Kīlauea's coastal plain was characterized using helicopter‐based visible and thermal structure‐from‐motion photogrammetry to construct a series of georeferenced digital surface models and thermal maps on eight different days. Results reveal key influences on the emplacement and evolution of such long‐lived pāhoehoe flows. This region of the flow grew by ~12 × 10 6 m 3 with a near‐constant time‐average discharge rate of 1.2–2.7 m 3 /s. The development of two tube systems is captured and shows an initial nascent tube enhanced by a narrow topographic confinement, which later inflated and created a topographic inversion that modulated the emplacement of a second flow lobe with its own tube system. The analysis of breakouts at various stages of the field's life suggests that the evolution of the thermal and morphological properties of the flow surface reflect its maturity. Thermal properties of breakouts were used to expand the empirical relationship of breakout cooling to longer timescales. This study contributes to the long‐term development and validation of more accurate predictive models for pāhoehoe, required during the management of long‐lasting lava flow crises in Hawai'i and elsewhere.

Hawaii↗
USGS2019

Slow-growing and extended-duration seismicity swarms: Reactivating joints or foliations in the Cahuilla Valley Pluton, Central Peninsular Ranges, Southern California

Three prolific earthquake swarms and numerous smaller ones have occurred since 1980 in the Mesozoic igneous plutonic rocks of the Perris block of the Peninsular Ranges, Southern California. The major swarms occurred in 1980–1981, 1983–1984, and 2016–2018, with the latest swarm still ongoing. These swarms have no clear mainshock, with the largest events of M L 3.6, M L 3.7, and M w 4.4. Each successive swarm had larger cumulative seismic moment release with about 314 and 411 events of M ≥ 1.5, while the third swarm has produced about 451 events of M ≥ 1.5 (as of September 2018). The concurrent strike‐slip faulting occurred on north to northwest striking planes but with no orthogonal northeast trending seismicity alignments. These shallow swarms are probably driven by intrablock Pacific‐North America plate boundary stress loading of the two bounding major late Quaternary strike‐slip faults, the Elsinore and San Jacinto faults. The state of stress within the Cahuilla Valley pluton has a ~40° angle between the maximum principal stress and the average trend of the swarms, suggesting that migrating pore fluid pressures aid in the formation and growth of zones of weakness. These swarms, which last more than 600 days each, exhibit clear bilateral spatial migration for distances of up to ~7–8 km and reach their full length in about 20 months. The slow spatial‐temporal development of the swarms corresponds to a fluid diffusivity of 0.006 to 0.01 m 2 /s, consistent with very low permeability rocks as expected for this block. There is no geodetic or other evidence for a slow slip event driving the swarms.

California↗
USGS2019

Heat flow in the Western Arctic Ocean (Amerasian Basin)

From 1963 to 1973 the U.S. Geological Survey (USGS) measured heat flow at 356 sites in the Amerasian Basin (Western Arctic Ocean) from a drifting ice island (T-3). The resulting measurements, which are unevenly distributed on Alpha-Mendeleev Ridge (AMR) and in Canada and Nautilus basins, greatly expand available heat flow data for the Arctic Ocean. Average T-3 heat flow is ~54.7 ± 11.3 mW m-2, and Nautilus Basin, including Mendeleev Plain, is the only well-surveyed area (~13% of data) with significantly higher average heat flow (63.8 mW m-2). Heat flow and bathymetry are not correlated at a large scale, and turbiditic surficial sediments (Canada and Nautilus basins) have higher heat flow than the sediments that blanket the AMR. Thermal gradients are mostly near-linear, implying that conductive heat transport dominates and that near-seafloor sediments are in thermal equilibrium with overlying bottom waters. Combining the heat flow data with modern seismic imagery suggests that some of the observed heat flow variability may be explained by local changes in sediment thickness or lithology or the presence of basement faults that channel circulating seawater. A thermal model that incorporates thermal conductivity variations along a profile from Canada Basin (thick sediment on mostly oceanic crust) to Alpha Ridge (thin sediment over thick magmatic units associated with the High Arctic Large Igneous Province) predicts heat flow lower than that observed on Alpha Ridge. This, along with other observations, implies that circulating fluids modulate conductive heat flow and contribute to high variability in the T-3 dataset. .

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research↗
USGS2020

How often can Earthquake Early Warning systems alert sites with high intensity ground motion?

Although numerous Earthquake Early Warning (EEW) algorithms have been developed we still lack a detailed understanding of how often and under what circumstances useful ground motion alerts can be provided to end-users. Here we analyze the alerting performance of the PLUM, EPIC and FinDer algorithms by running them retrospectively on the seismic strong motion data of the 219 earthquakes in Japan since 1996 that exceeded Modified Mercalli Intensity (MMI) of 4.5 on at least 10 sites (Mw 4.5-9.1). Our analysis suggests that, irrespective of the algorithm, EEW end-users should be prepared that EEW can often but not always provide useful ground motion alerts. A majority of sites with moderate-strong ground motion (MMI 5-6) can generally get at least a few seconds of warning time from all algorithms. If such shaking is caused by a shallow crustal event, around 50% of such sites receive alerts with warning times >5 s. Many sites with severe-extreme ground motion (MMI >=8) can be alerted successfully in the case of very large offshore earthquakes, but less than 20% can be alerted ahead of time if such shaking is caused by a shallow crustal event. Our results provide detailed quantitative insight into the expected alerting performance for EEW algorithms under realistic conditions. The main caveat is that the largest shallow crustal event in our data set has Mw7.0, i.e. the data set does not contain very large strike slip events.

Journal of Geophysical Research↗
USGS2019

Shear velocity structure from ambient noise and teleseismic surface wave tomography in the Cascades around Mount St. Helens

Mount St. Helens (MSH) lies in the forearc of the Cascades where conditions should be too cold for volcanism. To better understand thermal conditions and magma pathways beneath MSH, data from a dense broadband array are used to produce high‐resolution tomographic images of the crust and upper mantle. Rayleigh‐wave phase‐velocity maps and three‐dimensional images of shear velocity ( Vs ), generated from ambient noise and earthquake surface waves, show that west of MSH the middle‐lower crust is anomalously fast (3.95 ± 0.1 km/s), overlying an anomalously slow uppermost mantle (4.0–4.2 km/s). This combination renders the forearc Moho weak to invisible, with crustal velocity variations being a primary cause; fast crust is necessary to explain the absent Moho. Comparison with predicted rock velocities indicates that the fast crust likely consists of gabbros and basalts of the Siletzia terrane, an accreted oceanic plateau. East of MSH where magmatism is abundant, middle‐lower crust Vs is low (3.45–3.6 km/s), consistent with hot and potentially partly molten crust of more intermediate to felsic composition. This crust overlies mantle with more typical wave speeds, producing a strong Moho. The sharp boundary in crust and mantle Vs within a few kilometers of the MSH edifice correlates with a sharp boundary from low heat flow in the forearc to high arc heat flow and demonstrates that the crustal terrane boundary here couples with thermal structure to focus lateral melt transport from the lower crust westward to arc volcanoes.

Washington↗
USGS2020

Structural control on megathrust rupture and slip behavior: Insights from the 2016 Mw 7.8 Pedernales Ecuador earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw 7.7 earthquake.

Journal of Geophysical Research: Earth Surface↗