Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

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 937 records · Page 52Linked to original sources
USGS2022

A model of the spatiotemporal dynamics of soil carbon following coastal wetland loss applied to a Louisiana salt marsh in the Mississippi River Deltaic Plain

The potential for carbon sequestration in coastal wetlands is high due to protection of carbon (C) in flooded soils. However, excessive flooding can result in the conversion of the vegetated wetland to open water. This transition results in the loss of wetland habitat in addition to the potential loss of soil carbon. Thus, in areas experiencing rapid wetland submergence, such as the Mississippi River Delta, coastal wetlands could become a significant source of carbon emissions if land loss is not mitigated. To accurately assess the capacity of wetlands to store (or emit) carbon in dynamic environments, it is critical to understand the fate of soil carbon following the transition from vegetated wetland to open water. We developed a simple soil carbon model representing soil depths to 1 m using the data collected from a Louisiana coastal salt marsh in the Mississippi River Deltaic Plain to predict soil carbon density and stock following the transition from a vegetated salt marsh to an open water pond. While immediate effects of ponding on the distribution of carbon within the 1-m soil profile were apparent, there were no effects of ponding on the overall, integrated, carbon stocks 14 years, following wetland submergence. Rather, the model predicts that soil carbon losses in the first meter will be realized over long periods of time (∼200 years) due to changes in the source of carbon (biomass vs. mineral sediment) with minimal losses through mineralization.

Louisiana↗
USGS2022

A process-model perspective on recent changes in the carbon cycle of North America

Continental North America has been found to be a carbon (C) sink over recent decades by multiple studies employing a variety of estimation approaches. However, several key questions and uncertainties remain with these assessments. Here we used results from an ensemble of 19 state-of-the-art dynamic global vegetation models from the TRENDYv9 project to improve these estimates and study the drivers of its interannual variability. Our results show that North America has been a C sink with a magnitude of 0.37 ± 0.38 (mean and one standard deviation) PgC year −1 for the period 2000–2019 (0.31 and 0.44 PgC year −1 in each decade); split into 0.18 ± 0.12 PgC year −1 in Canada (0.15 and 0.20), 0.16 ± 0.17 in the United States (0.14 and 0.17), 0.02 ± 0.05 PgC year −1 in Mexico (0.02 and 0.02) and 0.01 ± 0.02 in Central America and the Caribbean (0.01 and 0.01). About 57% of the new C assimilated by terrestrial ecosystems is allocated into vegetation, 30% into soils, and 13% into litter. Losses of C due to fire account for 41% of the interannual variability of the mean net biome productivity for all North America in the model ensemble. Finally, we show that drought years (e.g., 2002) have the potential to shift the region to a small net C source in the simulations (−0.02 ± 0.46 PgC year −1 ). Our results highlight the importance of identifying the major drivers of the interannual variability of the continental-scale land C cycle along with the spatial distribution of local sink-source dynamics.

Journal of Geophysical Research Biogeosciences↗
USGS2023

Variation in leaf reflectance spectra across the California flora partitioned by evolutionary history, geographic origin, and deep time

Evolutionary relatedness underlies patterns of functional diversity in the natural world. Hyperspectral remote sensing has the potential to detect these patterns in plants through inherited patterns of leaf reflectance spectra. We collected leaf reflectance data from across the California flora from plants grown in a common garden. Regions of the reflectance spectra vary in the depth and strength of phylogenetic signal. We also show that these differences are much greater than variation due to the geographic origin of the plant. At the phylogenetic extent of the California flora, spectral variation explained by the combination of ecotypic variation (divergent evolution) and convergent evolution of disparate lineages was minimal (3 to 7 %) but statistically significant. Interestingly, at the extent of a single genus ( Arctostaphylos ) no unique variation could be attributed to geographic origin. However, up to 18% of the spectral variation among Arctostaphylos individuals was shared between phylogeny and intraspecific variation stemming from ecotypic differences (i.e., geographic origin). Future studies could conduct more structured experiments (e.g., transplants or observations along environmental gradients) to disentangle these sources of variation and include other intraspecific variation (e.g., plasticity). We constrain broad scale spectral variability due to ecotypic sources, providing further support for the idea that phylogenetic clusters of species might be detectable through remote sensing. Phylogenetic clusters could represent a valuable dimension of biodiversity monitoring and detection.

California↗
USGS2024

Reply to comment on "Five decades of observed daily precipitation reveal longer and more variable drought events across much of the western United States"

Paciorek and Wehner raise important questions around our use of the Mann-Kendall nonparametric trend test on smoothed data for analyzing long-term hydrometeorological trends in Zhang et al. (2021, https://doi.org/10.1029/2020gl092293 ). We thank them for initiating this important conversation and their gracious cooperation in exploring the issues addressed in their comment. In this reply we confirm the inflation of significant p -values by our choice to smooth, illustrate the relatively minor impacts on the main conclusions of our paper, and add our voices to those of Paciorek and Wehner in highlighting the lack of methodology for hypothesis testing across multiple stations that have spatial structure (i.e., testing for regionally consistent trends).

western United States↗
USGS2023

Accuracy of finite fault slip estimates in subduction zone regions with topographic Green's functions and seafloor geodesy

Until recently, the lack of seafloor geodetic instrumentation and the use of unrealistically simple, half-space based forward models have resulted in poor resolution of near-trench slip in subduction zone settings. Here, we use a synthetic framework to investigate the impact of topography and geodetic data distribution on coseismic slip estimates in various subduction zone settings. We calculate surface displacements in two synthetic topographic domains that have topography similar to that of Chile and Japan, respectively. We then attempt to image target slip distributions by using a Bayesian approach to solve for slip with two sets of Green's functions—one that accounts for topography and one that does not—and five sets of 50 or more observation points selected from the synthetic surface displacements. Three of these sets of observation points are entirely onland, and two include 5–10 seafloor geodetic sites. We find that the use of topographic Green's functions always improves inferred slip models, and with seafloor geodetic data, it enables an almost perfect recovery of a target slip model, even in the near-trench region. Critically, our results demonstrate that it would be impossible for non-topographic Green's functions to properly recover the true slip distribution, particularly in the near-trench region. We also perform a parameter study with approximately 4,000 slip models estimated using a least-square approach, and find that topographic Green's functions yield significantly more accurate slip models in cases where good data (well distributed and reasonably dense) are available, even in the absence of seafloor geodetic sites.

Journal of Geophysical Research: Solid Earth↗
USGS2023

Pre-eruptive outgassing and pressurization, and post-fragmentation bubble nucleation, recorded by vesicles in breadcrust bombs from vulcanian activity at Guagua Pichincha Volcano, Ecuador

Breadcrust bombs formed during Vulcanian eruptions are assumed to originate from the shallow plug or dome. Their rim to core texture reflects the competition between cooling and degassing timescales, which results in a dense crust with isolated vesicles contrasting with a highly vesicular vesicle network in the interior. Due to relatively fast quenching, the crust can shed light on pre- and syn-eruptive conditions prior to or during fragmentation, whereas the interior allows us to explore post-fragmentation vesiculation. Investigation of pre- to post-fragmentation processes in breadcrust bombs from the 1999 Vulcanian activity at Guagua Pichincha, Ecuador, via 2D and 3D textural analysis reveals a complex vesiculation history, with multiple, spatially localized nucleation and growth events. Large vesicles (Type 1), present in low number density in the crust, are interpreted as pre-eruptive bubbles formed by outgassing and collapse of a permeable bubble network during ascent or stalling in the plug. Haloes of small, syn-fragmentation vesicles (Type 2), distributed about large vesicles, are formed by pressurization and enrichment of volatiles in these haloes. The nature of the pressurization process in the plug is discussed in light of seismicity and ground deformation signals, and previous textural and chemical studies. A third population (Type 3) of post-fragmentation small vesicles appears in the interior of the bomb, and growth and coalescence of Type 2 and 3 vesicles causes the transition from isolated to interconnected bubble network in the interior. We model the evolution of viscosity, bubble growth rate, diffusion timescales, bubble radius and porosity during fragmentation and cooling. These models reveal that thermal quenching dominates in the crust whereas the interior undergoes a viscosity quench caused by degassing, and that the transition from crust to interior corresponds to the onset of percolation and development of permeability in the bubble network.

Guagua Pichincha Volcano↗
USGS2023

Quantifying site effects and their influence on earthquake source parameter estimations using a dense array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

Oklahoma↗
USGS2024

GRACE and GRACE Follow-On gravity observations of intermediate-depth earthquakes contrasted with those of shallow events

Earthquakes involve mass redistribution within the solid Earth and the ocean, and as a result, perturb the Earth's gravitational field. For most of the shallow (<60 km) earthquakes with M w > 8.0, the GRACE satellite gravity measurements suggest considerable volumetric disturbance of rocks. At a spatial scale of hundreds of km, the effect of volumetric change exceeds gravity change by vertical deformation; for example, negative gravity anomalies associated with volumetric expansion are characteristic patterns after shallow thrust events. In this study, however, we report contrasting observations of gravity change from two intermediate-depth (100–150 km) earthquakes of 2016 & 2017 Mw 8.0 (two combined) Papua New Guinea thrust faulting events and 2019 Mw 8.0 Peru normal faulting and highlight the importance of compressibility in earthquake deformation. The combined 2016/17 thrust events resulted in a positive gravity anomaly of 5–6 microGal around the epicenter, while the 2019 normal faulting produced a negative gravity anomaly of 3–4 microGal. Our modeling found that these gravity changes are manifestation of vertical deformation with limited volumetric change, distinct from gravity changes after the shallow earthquakes. The stronger resistance of rocks to volume change at intermediate-depth results in largely incompressible deformation and thus in a gravity change dominated by vertical deformation. In addition, malleable rocks under high pressure and temperature at depth facilitated substantial afterslip and/or fast viscoelastic relaxation causing additional vertical deformation and gravity change equivalent to the coseismic change. For the Papua New Guinea events, this means that postseismic relaxation enhanced coseismic uplift and relative sea level decrease.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Physics-based forecasts of eruptive vent locations at calderas

Constraining stresses in the Earth's crust in volcanic regions is critical for understanding many mechanical processes related to eruptive activity. Dike pathways, in particular, are shaped by the orientation of principal stress axes. Therefore, accurate models of dike trajectories and future vent locations rely on accurate estimates of stresses in the subsurface. This work presents a framework for probabilistic constraint of the stress state of calderas by combining three-dimensional physics-based dike pathway models with observed past vent locations using a Monte Carlo approach. The retrieved stress state is then used to produce probability maps of future vent opening across a caldera. We test our stress inversion and vent forecast approach on synthetic scenarios, and find it successful depending on the distribution of the available vents and the complexity of the volcano's structural history. We explore the potential and limitations of the approach, show how its performance is sensitive to the assumptions in the models and available prior information, and discuss how it may be applied to real calderas.

Journal of Geophysical Research: Solid Earth↗
USGS2024

Modeled coastal-ocean pathways of land-sourced contaminants in the aftermath of Hurricane Florence

Extreme precipitation during Hurricane Florence, which made landfall in North Carolina in September 2018, led to breaches of hog waste lagoons, coal ash pits, and wastewater facilities. In the weeks following the storm, freshwater discharge carried pollutants, sediment, organic matter, and debris to the coastal ocean, contributing to beach closures, algae blooms, hypoxia, and other ecosystem impacts. Here, the ocean pathways of land-sourced contaminants following Hurricane Florence are investigated using the Regional Ocean Modeling System (ROMS) with a river point source with fixed water properties from a hydrologic model (WRF-Hydro) of the Cape Fear River Basin, North Carolina's largest watershed. Patterns of contaminant transport in the coastal ocean are quantified with a finite duration tracer release based on observed flooding of agricultural and industrial facilities. A suite of synthetic events also was simulated to investigate the sensitivity of the river plume transport pathways to river discharge and wind direction. The simulated Hurricane Florence discharge event led to westward (downcoast) transport of contaminants in a coastal current, along with intermittent storage and release of material in an offshore (bulge) or eastward (upcoast) region near the river mouth, modulated by alternating upwelling and downwelling winds. The river plume patterns led to a delayed onset and long duration of contaminants affecting beaches 100 km to the west, days to weeks after the storm. Maps of the onset and duration of hypothetical water quality hazards for a range of weather conditions may provide guidance to managers on the timing of swimming/shellfishing advisories and water quality sampling.

Journal of Geophysical Research: Oceans↗
USGS2023

Bayesian hierarchical modeling for probabilistic estimation of tsunami amplitude from far-field earthquake sources

Evaluation of tsunami disaster risk for a coastal region requires reliable estimation of tsunami hazard, for example, wave amplitude close to the shore. Observed tsunami data are scarce and have poor spatial coverage, and for this reason probabilistic tsunami hazard analysis (PTHA) traditionally relies on numerical simulation of “synthetic” tsunami generation and propagation toward the coast. Such an approach has been extensively studied in the past and it is widely recognized as an important disaster-risk mitigation tool. PTHA can not only provide less uncertain and spatially coherent hazard estimates in comparison with classical empirical data analysis which is restricted at the tide gauge stations, but also local inundation information. In this paper, we explore a purely statistical alternative to traditional PTHA for evaluation of tsunami amplitude hazard. Here, we use tide gauge measurements of tsunami amplitude along the western United States, specifically California and Oregon, and develop a spatial Bayesian hierarchical model (BHM) to assess tsunami hazard from far-field earthquake sources at various recurrence intervals. The configuration of our model incorporates latent Gaussian fields that utilize information on the distance between tide gauges as well as on the continental shelf width, that is, a covariate linked to potential dissipative effects on wave energy as the tsunami travels over shallow water. Through our BHM, we produce spatially continuous probabilistic maps of far-field tsunami hazard which can aid comprehensive tsunami disaster risk reduction and management.

California, Oregon, Washington↗
USGS2023

VogCast: A framework for modeling volcanic air pollution and its application to the 2022 eruption of Mauna Loa Volcano, Hawai'i

Volcanic activity and the associated gas emissions into the atmosphere often result in adverse air quality conditions and present a hazard to human health and the environment. Building on a decade-long effort to provide operational surface sulfur dioxide and sulfate aerosol forecasts for the State of Hawai'i, we present an air quality modeling framework called VogCast. VogCast is designed to simplify ensemble air quality prediction on a regional scale by linking together multiple state-of-the-art models of meteorology, emissions, and dispersion. The framework is open-source and introduces a new dynamic plume-rise algorithm for distributing pollutants vertically. Using radar and satellite data, we demonstrate that VogCast reasonably captured the mean injection height, the location, and the general envelope of the vog plume during Mauna Loa's 2022 eruption. The results suggest that during the 12-day eruption period model performance varied between days with trade and non-trade wind conditions. Our findings also highlight the importance of sulfur dioxide emission rate and vent parameter inputs for improving forecast accuracy. The broad goal of this work is to better our understanding of vog dispersion and improve air quality prediction for impacted communities.

Hawaii↗
USGS2024

Measuring erosional and depositional patterns across Comet 67P's imhotep region

Comet 67P/Churyumov-Gerasimenko displays a pronounced hemispherical dichotomy in surface morphology, where the southern hemisphere exhibits more erosional features than the northern hemisphere due to receiving much greater solar radiation. Consequently, it is generally assumed that particles are ejected from the southern hemisphere through sublimation and a significant fraction eventually descends as airfall, covering the northern terrains. To investigate this south-to-north material transfer during the comet's perihelion passage, we used photoclinometry to measure material redistribution within its most extensive smooth terrain deposit around the Imhotep region. However, our findings do not align with this expected trend. Instead, we show that local-scale processes substantially impact the erosion and accumulation of material, with one area experiencing net erosion while another nearby region, just a few dozen meters away, sees sediment buildup. Our analysis underscores the complex interplay of processes shaping Comet 67P's surface and likely comets more generally.

Journal of Geophysical Research↗
USGS2023

Kinematic evolution of a large paraglacial landslide in the Barry Arm fjord of Alaska

Our warming climate is adversely affecting cryospheric landscapes via glacial retreat, permafrost degradation, and associated slope destabilization. In Prince William Sound, Alaska, the rapid retreat of Barry Glacier has destabilized the slopes flanking the glacier, resulting in numerous landslides. The largest of these landslides (∼500 Mm 3 in volume) is more than 2 km wide and has the potential to generate a tsunami that could affect nearby recreationists, marine traffic, infrastructure, natural and cultural resources, and the community of Whittier, located 60 km from the landslide. Here, we combine landslide structural and kinematic element mapping with data acquired from bi-yearly airborne lidar, multi-week satellite-based synthetic aperture radar (SAR), sub-hourly ground-based SAR, and seismic monitoring from 2020 to 2022 to characterize this landslide and examine its evolution. While some methods serve as a snapshot in time that is a culmination of events, others emphasize the ever-evolving nature of the landslide and associated hazards. Four major kinematic elements define the overall structure of the landslide, which vary in deformation type and rate, from creep (5 mm per day over several months) to episodic movement (2 m in 30 days) and landslide-wide to localized events. In some areas of the landslide, short-term deformation deviates from structures formed by cumulative movement, implying structural and kinematic evolution associated with glacier retreat. These insights are important for assessing landslide hazards and hazard evolution for large, slow-moving bedrock landslides in actively deglaciating environments.

Alaska↗
USGS2023

Tree-ring derived avalanche frequency and climate associations in a high-latitude, maritime climate

Snow avalanches are a natural hazard in mountainous areas worldwide with severe impacts that include fatalities, damage to infrastructure, disruption to commerce, and landscape disturbance. Understanding long-term avalanche frequency patterns, and associated climate and weather influences, improves our understanding of how climate change may affect avalanche activity. We used dendrochronological techniques to evaluate the historical frequency of large magnitude avalanches (LMAs) in the high-latitude climate of southeast Alaska, United States. We collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between 1720 and 2018, which allowed us to reconstruct 82 years with LMA activity across three sub-regions. By combining this tree-ring-derived avalanche data set with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found that February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of LMA activity in the study area. Specifically, LMA activity occurred during winters with substantial February and March precipitation and neutral or negative (cold) ONI values, while years not characterized by LMAs occur more frequently during warm winters (positive ONI values). Our examination of the climate-avalanche relationship in southeast Alaska sheds light on important climate variables and physical processes associated with LMA years. These results can be used to inform long-term infrastructure planning and avalanche mitigation operations in an urban area, such as Juneau, where critical infrastructure is subject to substantial avalanche hazard.

Alaska↗
USGS2024

Pockmarks offshore Big Sur, California provide evidence for recurrent, regional, and unconfined sediment gravity flows

Recent surface ship multibeam surveys of the Sur Pockmark Field, offshore Central California, reveal >5,000 pockmarks in an area that is slated to host a wind farm, between 500- and 1,500-m water depth. Extensive fieldwork was conducted to characterize the seafloor environment and its recent geologic history, including visual observations with remotely operated vehicles, sediment core sampling, and high-resolution, near-bottom Chirp and multibeam surveys collected with autonomous underwater vehicles to capture the morphology and stratigraphy of the pockmarks. No evidence of high methane concentrations in sediments, chemosynthetic biological communities, or methane-derived diagenetic byproducts was found. Chirp data and sediment cores showed alternating layers of slowly accumulating hemipelagic drapes interrupted by more reflective turbidite horizons that extend throughout the pockmark field and beyond. Chirp data showed multiple episodes of lateral migration over time in some of the pockmarks in association with erosion and infilling events. Laterally continuous turbidite horizons that overlay erosional surfaces indicated that pockmark migration occurred synchronously in multiple pockmarks separated by tens of kilometers. These shifts are presumed to be the result of asymmetrical erosion of the pockmark flanks caused by passing sediment gravity flows. While some pockmarks occur in chains, most are not clustered or randomly spaced but are regularly dispersed within the pockmark field. We hypothesize that intermittent, unconfined sediment gravity flows occurring over at least the last 280,000 years are the source of the regionally continuous turbidite deposits and the mechanism that maintained the regularly dispersed pockmarks.

California↗
USGS2024

Predicting characteristic length scales of barrier island segmentation in microtidal environments

Segmented barrier islands can be found in regions with small tidal ranges. In contrast to tidally dominated barriers, where inlet dynamics are thought to control island length scales, the controls on barrier island length scales in wave-dominated environments have not been quantified. These microtidal barriers typically have a curved shoreline, suggesting the influence of wave-driven alongshore sediment transport. Microtidal barriers are also typically hydrodynamically isolated from one another, as weak tidal flows limit interactions between adjoining barriers. To better understand the controls on and scales of barrier segmentation in the relative absence of tides, here we develop a theoretical framework to estimate the alongshore length scales at which a barrier will either breach or heal following a disturbance in the barrier morphology. The non-dimensional framework compares the timescales of overwash (advective) and alongshore sediment transport (diffusive) processes along barrier island chains. We then apply this framework to modern barrier islands in the microtidal Gulf of Mexico using wave hindcast data and the lengths, widths, heights, and lagoon depths measured from remotely sensed geospatial data and topobathymetric data. We find that most of these barriers are currently longer than their critical length scale, often as a result of coastal restoration efforts. Our critical length scale analysis suggests that most of the Gulf of Mexico barriers are vulnerable to segmentation despite coastal restoration efforts intended to protect fisheries and the mainland coasts.

Florida, Louisiana, Mississippi↗
USGS2024

Source, migration pathways, and atmospheric release of geologic methane associated with the complex permafrost regimes of the outer Mackenzie River Delta, Arctic, Canada

Sources and fluxes of methane to the atmosphere from permafrost are significant but poorly constrained in global climate models. We present data collected from the variable permafrost setting of the outer Mackenzie River Delta, including observations of aquatic methane seepage, core determinations of in situ methane occurrence and seep gas isotope geochemistry. The sources and locations of in situ geologic methane occurrence and aquatic and atmospheric gas release appear to be controlled by the regional geology and permafrost conditions. Where permafrost is >250 m thick, thermogenic gas deposits at depth are isolated by laterally continuous, low permeability ice-bearing sediments with few through-going thawed taliks. Thus, the observed in situ methane and aquatic gas seepage appears to be dominated by microbial methane. In contrast, where permafrost is <80 m thick, taliks are more likely to be through-going, providing permeable conduits from depth and migration pathways for both thermogenic and biogenic gas. Continuous annual fluid sampling of two lakes and a river channel documents aquatic methane flux from microbial sources, more deeply buried thermogenic sources, and mixtures of both. Using estimates of in situ methane concentration from deep core samples and observations of in situ free gas occurrences, we conclude that the reservoir of in situ geologic methane within ice bonded permafrost is substantial and that this methane is presently migrating with ongoing atmospheric release. It is our assessment that the permafrost setting, and processes described are sensitive to future climate change as the permafrost warms.

Mackenzie River Delta↗