Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research”

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 991 records · Page 55Linked to original sources
USGS2019

Carbon chemistry of intact versus chronically drained peatlands in the southeastern USA

The Great Dismal Swamp (GDS) is a large temperate swamp in Virginia/North Carolina with peat soils historically resistant to microbial decomposition. However, this peatland has been subject to ~200 years of disturbance during which extensive drainage, fire suppression, and wide-spread logging have increased decomposition and dramatically decreased the distribution of Atlantic white cedar (AWC). The purpose of this study was to determine the impact of long-term drainage and AWC loss on the carbon chemistry of GDS peats. Peat cores were collected from three drained GDS vegetation communities (pocosin, AWC, and red maple-black gum) and compared to cores collected from an intact, undrained AWC peatland at the Alligator River National Wildlife Refuge (AR) in North Carolina, USA. The AR peats had higher lignin content in the deeper peat intervals, and lignin content and % organic carbon were largely invariant with depth compared to the GDS peats. The concentrations of syringyl group phenols were greater in the surface layers of GDS peats, likely reflecting the selective removal of AWC and transition from gymnosperms to angiosperms. Acid to aldehyde ratios for vanillyl and syringyl group phenols indicated that the GDS peats were more decomposed, particularly at depth, and that this occurred under aerobic conditions. Moreover, solid state 13C NMR confirmed a coincident loss of carbohydrates and increase in recalcitrant byproducts of carbohydrate degradation with depth. These data indicate that long-term drainage has accelerated the decomposition of peat at the GDS, reducing the capacity and stability of the carbon sink.

Virginia, North Carolina↗
USGS2019

Complex response of sediment phosphorus to land use and management within a river network

Rivers affected by anthropogenic nutrient inputs can retain some of the phosphorus (P) load through sediment retention and burial. Determining the influence of land use and management on sediment P concentrations and P retention in fluvial ecosystems is challenging because of different stressors operating at multiple spatial and temporal scales. In this study, we sought to determine how land use and management influenced sediment P concentrations and P retention within a river network draining a watershed modified by mixed land use activities; the Fox River, Wisconsin, USA. Results show current land use had no detectable effect on sediment P concentrations, and only a small potential effect on P retention capacity. However, sites draining predominantly forested areas were associated with riverbed sediments less saturated in P, whereas sites draining mainly agricultural areas were more likely to release sediment-bound P. Current management actions, including the implementation of best management practices, had a small positive effect on P retention capacity but had no observed effect on sediment P concentrations. Our results suggest that fine riverbed sediment throughout the Fox River network is saturated with P and has the capability to release P when water column P concentrations are low. These P-saturated riverbed sediments represent a potential legacy effect of past land use disturbances on P dynamics in the Fox River network.

Wisconsin↗
USGS2019

Heterogeneity in hyporheic flow, pore water chemistry, and microbial community composition in an alpine streambed

The hyporheic zone, where surface water and groundwater mix, is an important microbial habitat where biogeochemical reactions influence water quality. We show that spatial variability in hyporheic flow in the East River near Crested Butte, Colorado, drives heterogeneity in streambed geochemical conditions and microbial community assemblages, but the diversity of assemblages remains nearly constant throughout the reach. In July 2018, we collected approximately 100 pore water samples at 20 cm depth and analyzed them for anions, cations, dissolved organic carbon, dissolved organic matter (DOM) quality, and basic water quality parameters. Vertical hydraulic head gradients were also measured to assess the potential for upward or downward flow, and heat tracing was used to quantify vertical flux rates at a subset of locations. We found that regions of the streambed that are more groundwater-dominated contain less dissolved oxygen, higher concentrations of reduced metals, and more microbially-processed, recalcitrant DOM, while more surface water-dominated locations contain higher dissolved oxygen concentrations and terrestrially-derived, labile DOM. 16S rRNA gene sequencing of extracted DNA revealed that microbial community composition varies with geochemical gradients related to hyporheic flow. These findings provide a better understanding of hyporheic controls on streambed biogeochemistry during the baseflow season, which is expected to lengthen with climate change in alpine watersheds due to earlier snowmelt onset and reduced snowpack.

Colorado↗
USGS2020

Increasing rates of carbon burial in southwest Florida coastal wetlands

Rates of organic carbon (OC) burial in some coastal wetlands appear to be greater in recent years than they were in the past. Possible explanations include ongoing mineralization of older OC or the influence of an unaccounted‐for artefact of the methods used to measure burial rates. Alternatively, the trend may represent real acceleration in OC burial. We quantified OC burial rates of mangrove and coastal freshwater marshes in southwest Florida through a comparison of rates derived from 210 Pb, 137 Cs, and surface marker horizons (MHs). Age/depth profiles of lignin: OC were used to assess whether down‐core remineralization had depleted the OC pool relative to lignin, and lignin phenols were used to quantify the variability of lignin degradation. Over the past 120 years, OC burial rates at seven sites increased by factors ranging from 1.4 to 6.2. We propose that these increases represent net acceleration. Change in relative sea‐level rise is the most likely large‐scale driver of acceleration, and sediment deposition from large storms can contribute to periodic increases. Mangrove sites had higher OC and lignin burial rates than marsh sites, indicating inherent differences in OC burial factors between the two habitat types. The higher OC burial rates in mangrove soils mean that their encroachment into coastal freshwater marshes has the potential to increase burial rates in those locations even more than might be expected from the acceleration trends. Regionally, these findings suggest that burial represents a substantially growing proportion of the coastal wetland carbon budget.

Florida↗
USGS2020

Impacts of seagrass dynamics on the coupled long-term evolution of barrier-marsh-bay systems

Seagrass provides a wide range of economically and ecologically valuable ecosystem services, with shoreline erosion control often listed as a key service, but can also alter the sediment dynamics and waves within back-barrier bays. Here we incorporate seagrass dynamics into an existing barrier-marsh exploratory model, GEOMBEST++, to examine the coupled interactions of the back-barrier bay with both adjacent (marsh) and non-adjacent (barrier island) subsystems. While seagrass reduces marsh edge erosion rates and increases progradation rates in many of our 288 model simulations, seagrass surprisingly increases marsh edge erosion rates when sediment export from the back-barrier basin is negligible because the ability of seagrass to reduce the volume of marsh sediment eroded matters little for back-barrier basins in which all sediment is conserved. Our model simulations also suggest that adding seagrass to the bay subsystem leads to increased deposition in the bay, reduced sediment available to the marsh, and enhanced marsh edge erosion until the bay reaches a new, shallower equilibrium depth. In contrast, removing seagrass liberates previously-sequestered sediment that is then delivered to the marsh, leading to enhanced marsh progradation. Lastly, we find that seagrass reduces barrier island migration rates in the absence of back-barrier marsh by filling accommodation space in the bay. These model observations suggest that seagrass meadows operate as dynamic sources and sinks of sediment that can influence the evolution of coupled marsh and barrier island landforms in unanticipated ways.

Virginia↗
USGS2020

Key components and contrasts in the nitrogen budget across a US-Canadian transboundary watershed

Watershed nitrogen (N) budgets provide insights into drivers and solutions for groundwater and surface water N contamination. We constructed a comprehensive N budget for the transboundary Nooksack River Watershed (British Columbia, Canada, and Washington, USA) using locally derived data, national statistics, and standard parameters. Feed imports for dairy (mainly in the United States) and poultry (mainly in Canada) accounted for 30% and 29% of the total N input to the watershed, respectively. Synthetic fertilizer was the next largest source contributing 21% of inputs. Food imports for humans and pets together accounted for 9% of total inputs, lower than atmospheric deposition (10%). N imported by returning salmon representing marine‐derived nutrients accounted for <0.06% of total N input. Quantified N export was 80% of total N input, driven by ammonia emission (32% of exports). Animal product export was the second largest output of N (31%) as milk and cattle in the United States and poultry products in Canada. Riverine export of N was estimated at 28% of total N export. The commonly used crop nitrogen use efficiency (NUE) metric alone did not provide sufficient information on farming activities but in combination with other criteria such as farm‐gate NUE may better represent management efficiency. Agriculture was the primary driver of N inputs to the environment as a result of its regional importance; the N budget information can inform management to minimize N losses. The N budget provides key information for stakeholders across sectors and borders to create environmentally and economically viable and effective solutions.

Washington, British Columbia↗
USGS2021

Estimating the impact of seep methane oxidation on ocean pH and dissolved inorganic radiocarbon along the U.S. mid‐Atlantic Bight

Ongoing ocean warming can release methane (CH 4 ) currently stored in ocean sediments as free gas and gas hydrates. Once dissolved in ocean waters, this CH 4 can be oxidized to carbon dioxide (CO 2 ). While it has been hypothesized that the CO 2 produced from aerobic CH 4 oxidation could enhance ocean acidification, a previous study conducted in Hudson Canyon shows that CH 4 oxidation has a small short‐term influence on ocean pH and dissolved inorganic radiocarbon. Here we expand upon that investigation to assess the impact of widespread CH 4 seepage on CO 2 chemistry and possible accumulation of this carbon injection along 234 km of the U.S. Mid‐Atlantic Bight. Consistent with the estimates from Hudson Canyon, we demonstrate that a small fraction of ancient CH 4 ‐derived carbon is being assimilated into the dissolved inorganic radiocarbon (mean fraction of 0.5 ± 0.4 %). The areas with the highest fractions of ancient carbon coincide with elevated CH 4 concentration and active gas seepage. This suggests that aerobic CH 4 oxidation has a greater influence on the dissolved inorganic pool in areas where CH 4 concentrations are locally elevated, instead of displaying a cumulative effect downcurrent from widespread groupings of CH 4 seeps. An upper limit approximation of the input rate of ancient‐derived DIC into the waters overlying the northern U. S Mid‐Atlantic Bight further suggests that oxidation of ancient CH 4 ‐derived carbon is not negligible on the global scale and could contribute to deep‐water acidification over longer time scales.

Journal of Geophysical Research- Biogeosciences↗
USGS2020

Combined seismic and geodetic analysis before, during and after the 2018 Mt. Etna eruption

In December 2018, Etna volcano experienced one of the largest episodes of unrest since the installation of geophysical monitoring networks in 1970. The unrest culminated in a short eruption with a small volume of lava erupted, a significant seismic crisis and deformation of the entire volcanic edifice of magnitude never recorded before at Mount Etna. Here we describe the evolution of the 2018 eruptive cycle from the analysis of seismic and geodetic data collected in the months preceding, during, and following the intrusion. We model the space-time evolution of high-rate deformation data starting from the active source previously identified from deformation data and the propagation of seismicity in a 3-D velocity model. The intrusion model suggests emplacement of two dikes: a smaller dike located beneath the eruptive fissure and a second, deeper dike between 1 and 5 km below sea level that opened ~2 m. The rise and eruption of magma from the shallower dike did not interrupt the pressurization of a long-lasting deeper reservoir (~6 km) that induced continuous inflation and intense deformation of the eastern flank. Shortly after the intrusion, on 26 December 2018, a M L 4.8 earthquake occurred near Pisano, destroying buildings and roads in two villages. We propose a time-dependent intrusion model that supports the hypothesis of the inflation inducing flank deformation and that this process has been active since September 2018.

Mt. Etna↗
USGS2020

Hydrologically induced deformation in Long Valley Caldera and adjacent Sierra Nevada

Vertical and horizontal components of GNSS displacements in the Long Valley Caldera and adjacent Sierra Nevada range show a clear correlation with hydrological trends at both multiyear and seasonal time scales. We observe a clear vertical and horizontal seasonal deformation pattern primarily attributable to the solid earth response to hydrological surface loading at large-to-regional (Sierra Nevada range) scales. Several GNSS sites, located at the eastern edge of the Sierra Nevada along the southwestern rim of Long Valley Caldera, also show significant horizontal deformation that cannot be explained by elastic deformation from surface loading. Due to the location of these sites and the strong correlation between their horizontal displacements and spring discharge, we hypothesize that this deformation reflects poroelastic processes related to snowmelt runoff water infiltrating into the Sierra Nevada slopes around Long Valley Caldera. Interestingly, this is also an area where water infiltrates to feed the local hydrothermal system, and where snowmelt-induced earthquake swarms have been recently detected.

California, Nevada↗
USGS2020

The role of pre-magmatic rifting in shaping a volcanic continental margin: An example from the Eastern North American Margin

Both magmatic and tectonic processes contribute to the formation of volcanic continental margins. Such margins are thought to undergo extension across a narrow zone of lithospheric thinning (~100 km). New observations based on existing and reprocessed data from the Eastern North American Margin contradict this hypothesis. With ~64,000 km of 2‐D seismic data tied to 40 wells combined with published refraction, deep reflection, receiver function, and onshore drilling efforts, we quantified along‐strike variations in the distribution of rift structures, magmatism, crustal thickness, and early post‐rift sedimentation under the shelf of Baltimore Canyon Trough (BCT), Long Island Platform, and Georges Bank Basin (GBB). Results indicate that BCT is narrow (80–120 km) with a sharp basement hinge and few rift basins. The seaward dipping reflectors (SDR) there extend ~50 km seaward of the hinge line. In contrast, the GBB is wide (~200 km), has many syn‐rift structures, and the SDR there extend ~ 200 km seaward of the hinge line. Early post‐rift depocenters at the GBB coincide with thinner crust suggesting “uniform” thinning of the entire lithosphere. Models for the formation of volcanic margins do not explain the wide structure of the GBB. We argue that crustal thinning of the BCT was closely associated with late syn‐rift magmatism, whereas the broad thinning of the GBB segment predated magmatism. Correlation of these variations to crustal terranes of different compositions suggests that the inherited rheology determined the premagmatic response of the lithosphere to extension.

Journal of Geophysical Research-- Solid Earth↗
USGS2020

Shallow basin structure and attenuation are key to predicting long shaking duration in Los Angeles Basin

Ground motions in the Los Angeles Basin during large earthquakes are modulated by earthquake ruptures, path effects into the basin, basin effects, and local site response. We analyzed the direct effect of shallow basin structures on shaking duration at a period of 2–10 s in the Los Angeles region through modeling small magnitude, shallow, and deep earthquake pairs. The source depth modulates the basin response, particularly the shaking duration, and these features are a function of path effect and not site condition. Three-dimensional simulations using the CVM-S4.26.M01 velocity model show good fitting to the initial portion of the waveforms at periods of 5 s and longer but fail to predict the long shaking duration during shallow events, especially at periods less than 5 s. Simulations using CVM-H do not match the timing of the initial arrivals as well as CVM-S4.26.M01, and the strong late arrivals in the CVM-H simulation travel with an apparent velocity slower than observed. A higher-quality factor than traditionally assumed may produce synthetics with longer durations but is unable to accurately match the amplitude and phase. Beamforming analysis using dense array data further reveals the long duration surface waves have the same back azimuth as the direct arrivals and are generated at the basin edges, while the later coda waves are scattered from off-azimuth directions, potentially due to strong, sharp boundaries offshore. Improving the description of these shallow basin structures and attenuation model will enhance our capability to predict long-period ground motions in basins.

California↗
USGS2020

Origin and properties of hydrothermal tremor at Lone Star Geyser, Yellowstone National Park, USA

Geysers are rare geologic features that intermittently discharge liquid water and steam driven by heating and decompression boiling. The cause of variability in eruptive styles and the associated seismic signals are not well understood. Data collected from five broadband seismometers at Lone Star Geyser, Yellowstone National Park are used to determine the properties, location, and temporal patterns of hydrothermal tremor. The tremor is harmonic at some stages of the eruption cycle and is caused by near‐periodic repetition of discrete seismic events. Using the polarization of ground motion, we identify the location of tremor sources throughout several eruption cycles. During preplay episodes (smaller eruptions preceding the more vigorous major eruption), tremor occurs at depths of 7–10 m and is laterally offset from the geyser's cone by ~5 m. At the onset of the main eruption, tremor sources migrate laterally and become shallower. As the eruption progresses, tremor sources migrate along the same path but in the opposite direction, ending where preplay tremor originates. The upward and then downward migration of tremor sources during eruptions are consistent with warming of the conduit followed by evacuation of water during the main eruption. We identify systematic relations among the two types of preplays, discharge, and the main eruption. A point‐source moment tensor fit to low‐frequency waveforms of an individual tremor event using half‐space velocity models indicates average V S ≳ 0.8 km/s, source depths ~4–20 m, and moment tensors with primarily positive isotropic and compensated linear vector dipole moments.

Wyoming↗
USGS2020

Along-margin variations in breakup volcanism at the Eastern North American Margin

We model the magnetic signature of rift-related volcanism to understand the distribution and volumeofmagmatic activity that occurred during the breakup of Pangaea and early Atlantic opening at the Eastern North American Margin (ENAM).Along-strike variations in the amplitude and character of the prominent East Coast Magnetic Anomaly (ECMA) suggest that the emplacement of the volcanic layers producing this anomaly similarly varied along the margin. We use three-dimensional magnetic forward modeling constrained by seismic interpretationsto identify along-margin variations in volcanic thickness and width that can explain the observed amplitude and character of the ECMA. Our model results suggest that the ECMA is produced by a combination of both first-order (~600-1000 km)and second-order (~50-31100 km) magmatic segmentation. The first-order magmatic segmentation could have resulted from preexisting variations in crustal thickness and rheology developed during the tectonic amalgamation of Pangaea. The second-order magmatic segmentation developed during continental breakup and likely influenced the segmentation and transform fault spacing of the initial, and modern, Mid-Atlantic Ridge. These variations in magmatism showhow extension and thermal weakening was distributed at the ENAM during continental breakup and how this breakup magmatism was related to both previous and subsequent Wilson Cycle stages.

East Coast of United States, Atlantic Ocean↗
USGS2020

The induced Mw 5.0 March 2020 west Texas seismic sequence

On March 26, 2020, a M 5.0 earthquake occurred in the Delaware Basin, Texas, near the border between Reeves and Culberson Counties. This was the third largest earthquake recorded in Texas and the largest earthquake in the Central and Eastern United States since the three M 5.0–5.8 induced events in Oklahoma during 2016. Using multistation waveform template matching, we detect 3,940 earthquakes in the sequence with the first event in the area occurring in May 2018. The M 5.0 earthquake sequence occurred on a ENE (∼082°) normal fault dipping ∼37° toward the south. The earthquake caused 6 mm of oblique surface deformation, and geodetic slip inversion suggests slip was isolated above 6 km depth. We find that the sequence was most likely induced by nearby wastewater disposal operations, and seismicity rates in the region surrounding the M 5.0 will likely continue to increase in the future if disposal operations continue unaltered.

Texas↗
USGS2021

Airborne geophysical imaging of weak zones on Iliamna Volcano, Alaska: Implications for slope stability

Water‐saturated, hydrothermally altered rocks reduce the strength of volcanic edifices and increase the potential for sector collapses and far‐traveled mass flows of unconsolidated debris. Iliamna Volcano is an andesitic stratovolcano located on the western side of the Cook Inlet, ∼225 km southwest of Anchorage and is a source of repeated avalanches. The widespread snow and ice cover on Iliamna Volcano make surface alteration difficult to identify. However, intense hydrothermal alteration significantly reduces both the electrical resistivity and magnetization of volcanic rock and can therefore be identified with airborne geophysical measurements. We use airborne electromagnetic and magnetic data to map snow and ice thickness and identify underlying alteration zones at Iliamna Volcano, Alaska. Resistivities were calculated to an average depth of >300 m, and a 3‐D susceptibility model extends from the surface to the base of the volcano, about 3,000 m below the summit. Geophysical models image low resistivity (<30 ohm‐m) and low susceptibilities near the summit of Iliamna and below its older vent complex, with the low susceptibilities indicating alteration up to ∼800 m in thickness. Thin conductors (∼50–100 m thick) on the edifice slopes coincide with recorded locations of repeated debris avalanches over the past ∼60 years and are attributed to saturated zones at high elevation. Three‐dimensional slope stability models based upon the geophysically constrained alteration distribution suggest the edifice of Iliamna is unstable and could lead to collapse scars ∼400 m deep near the current and former vent complexes.

Alaska↗
USGS2021

Characteristics of frequent dynamic triggering of microearthquakes in Southern California

Dynamic triggering of earthquakes has been reported at various fault systems. The triggered earthquakes are thought to be caused either directly by dynamic stress changes due to the passing seismic waves, or indirectly by other nonlinear processes that are initiated by the passing waves. Distinguishing these physical mechanisms is difficult because of the general lack of high‐resolution earthquake catalogs and robust means to quantitatively evaluate triggering responses, particularly, delayed responses. Here we use the high‐resolution Quake Template Matching catalog in Southern California to systematically evaluate teleseismic dynamic triggering patterns in the San Jacinto Fault Zone and the Salton Sea Geothermal Field from 2008 to 2017. We develop a new statistical approach to identify triggered cases, finding that approximately 1 out of every 5 global M w ≥ 6 earthquakes dynamically trigger microearthquakes in Southern California. The triggering responses include both instantaneous and delayed triggering, showing a highly heterogeneous pattern and indicating possible evolving triggering thresholds. We do not observe a clear peak ground velocity triggering threshold that can differentiate triggering earthquakes from non‐triggering events, but there are subtle differences in the frequency content that may possibly differentiate the earthquakes. In contrast to the depth distribution of background seismicity, the identified triggered earthquakes tend to concentrate at the edges of the seismogenic zones. Although instantaneously triggered earthquakes are likely a result of dynamic Coulomb stress changes, the cases of delayed dynamic triggering are best explained by nonlinear triggering processes, including cyclic material fatigue, accelerated transient creep, and stochastic frictional heterogeneities.

California↗
USGS2021

Comparative rhyolite systems: Inferences from vent patterns and eruptive episodicities: Eastern California and Laguna del Maule

Distilling my experience in having field mapped in detail the volcanic fields at Laguna del Maule and Long Valley and having worked out their time-volume-composition magmatic histories, I compare and contrast the postglacial rhyolites of the former with six multi-vent eruptive sequences of rhyolite in California. Compilations and discussions are made of volcanic-field areas and longevities, their compositions, vent distributions, individual batch and total volumes, eruptive episodicities, and tectonic influences. Growth of long-lived pluton-scale reservoirs of granitic crystal mush, from which the rhyolite melts separated, are interpreted in terms of conceptual models I published previously—(1) fundamentally basaltic transcrustal magmatism, 1981; (2) the deep-crustal MASH zone model, 1988; and (3) the rhyolite-melt crystal-mush model, 2001. Inferences and speculations are advanced concerning processes and timescales of rhyolite-melt separation from granitic mush and of prompt or long-delayed subsequent eruption.

California↗
USGS2021

The birth of a Hawaiian fissure eruption

Most basaltic explosive eruptions intensify abruptly, allowing little time to document processes at the start of eruption. One opportunity came with the initiation of activity from fissure 8 (F8) during the 2018 eruption on the lower East Rift Zone of Kīlauea, Hawaii. F8 erupted in four episodes. We recorded 28 min of high‐definition video during a 51‐min period, capturing the onset of the second episode on 5 May. From the videos, we were able to analyze the following in‐flight parameters: frequency and duration of explosions; ejecta heights; pyroclast exit velocities; in‐flight total mass and estimated mass eruption rates; and the in‐flight total grain size distributions. The videos record a transition from initial pulsating outgassing, via spaced, but increasingly rapid, discrete explosions, to quasisustained, unsteady fountaining. This transition accompanied waxing intensity (mass flux) of the F8 eruption. We infer that all activity was driven by a combination of the ascent of a coupled mixture of small bubbles and melt, and the buoyant rise of decoupled gas slugs and/or pockets. The balance between these two types of concurrent flow determined the exact form of the eruptive activity at any point in time, and changes to their relative contributions drove the transition we observed at early F8. Qualitative observations of other Hawaiian fountains at Kīlauea suggest that this physical model may apply more generally. This study demonstrates the value of in‐flight parameters derived from high‐resolution videos, which offer a rapid and highly time‐sensitive alternative to measurements based on sampling of deposits posteruption.

Island of Hawai'i↗