USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

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
USGS2023

Mapping landslide susceptibility over large regions with limited data

Landslide susceptibility maps indicate the spatial distribution of landslide likelihood. Modeling susceptibility over large or diverse terrains remains a challenge due to the sparsity of landslide data (mapped extent of known landslides) and the variability in triggering conditions. Several different data sampling strategies of landslide locations used to train a susceptibility model are used to mitigate this challenge. However, to our knowledge, no study has systematically evaluated how different sampling strategies alter a model's predictor effects (i.e., how a predictor value influences the susceptibility output) critical to explaining differences in model outputs. Here, we introduce a statistical framework that examines the variation in predictor effects and the model accuracy (measured using receiver operator characteristics) to highlight why certain sampling strategies are more effective than others. Specifically, we apply our framework to an array of logistic regression models trained on landslide inventories collected at sub-regional scales over four terrains across the United States. Results show significant variations in predictor effects depending on the inventory used to train the models. The inconsistent predictor effects cause low accuracies when testing models on inventories outside the domain of the training data. Grouping test and training sets according to physiographic and ecological characteristics, which are thought to share similar triggering mechanisms, does not improve model accuracy. We also show that using limited landslide data distributed uniformly over the entire modeling domain is better than using dense but spatially isolated data to train a model for applications over large regions.

California, Kentucky, North Carolina, West Virgini↗
USGS2022

The influence of drying on the aeolian transport of river-sourced sand

Transgression and regression of water levels (stages) have impacted the evolution of aeolian landforms and sedimentary deposits throughout geologic history. We studied this phenomenon over a five-day period of reduced flow on the Colorado River in Grand Canyon National Park, AZ, USA, in March 2021. These transient low flows exposed river-channel sand deposits to the air, causing progressive desiccation (drying) and thereby making these deposits susceptible to aeolian transport. We measured aeolian threshold friction velocities ( u *t ) for sand saltation and PM10 dust emissions, as well as other characteristics, on a subaerially exposed sandbar and downwind aeolian dunefield during each day of the low river flow. The sandbar transitioned from supply-limited to transport-limited aeolian sediment transport conditions during the regression in river water stage. A possible tipping point between the two transport conditions occurred approximately 48 hours after the drop in river flow. The empirically measured u *t decreased as the sandbar sediment dried with increased subaerial exposure time. Theoretical estimates and empirical measurements of u *t corresponded closely on the aeolian dunefield and on the sandbar when it was drier during the third and fourth day of the experiment. Eighty-seven percent of the variability in u *t was explained by empirical models that provide practical estimates of aeolian transport potential of subaerial river sediment deposits using monitoring data that are commonly available in this and other river systems. The work provides theoretical insight into the response of aeolian processes to sediment supply changes driven by periods of anthropogenic activity, drought, and climate change.

Arizona↗
USGS2023

Impacts of spontaneous waterfall development on bedrock river longitudinal profile morphology

River profiles are shaped by climatic and tectonic history, lithology, and internal feedbacks between flow hydraulics, sediment transport and erosion. In steep channels, waterfalls may self-form without changes in external forcing (i.e., autogenic formation) and erode at rates faster or slower than an equivalent channel without waterfalls. We use a 1-D numerical model to investigate how self-formed waterfalls alter the morphology of bedrock river longitudinal profiles. We modify the standard stream power model to include a slope threshold above which waterfalls spontaneously form and a rate constant allowing waterfalls to erode faster or slower than other fluvial processes. Using this model, we explore how waterfall formation alters both steady state and transient longitudinal profile forms. Our model predicts that fast waterfalls create km-scale reaches in a dynamic equilibrium with channel slope held approximately constant at the threshold slope for waterfall formation, while slow waterfalls can create local channel slope maxima at the location of slow waterfall development. Furthermore, slow waterfall profiles integrate past base level histories, leading to multiple possible profile forms, even at steady-state. Consistency between our model predictions and field observations of waterfall-rich rivers in the Kings and Kaweah drainages in the southern Sierra Nevada, California, supports the hypothesis that waterfall formation can modulate river profiles in nature. Our findings may help identify how bedrock channels are influenced by waterfall erosion and aid in distinguishing between signatures of external and internal perturbations, thereby strengthening our ability to interpret past climate and tectonic changes from river longitudinal profiles.

Journal of Geophysical Research - Earth Surface↗
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

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↗
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↗
USGS2023

Viscous relaxation of Oort and Edgeworth craters on Pluto: Possible indicators of an epoch of early high heat flow

Impact craters, with their well-defined initial shapes, have proven useful as heat flow probes of a number of icy bodies, provided characteristics of viscous relaxation can be identified. For Pluto's numerous craters, such identifications are hampered/complicated by infilling and erosion by mobile volatile ices, but not in every case. Large craters offer relatively deep probes of rheological structure, and on Pluto two large old craters in a major dark (volatile-ice free) region are probably the best examples for possible viscous relaxation: Oort (115-km diameter) and Edgeworth (140-km diameter). They are similar in size, location, and apparent age (morphological preservation), but may or may not be coeval. Edgeworth is particularly shallow and its floor appears bowed up above the original ground plane, a classic hallmark of viscous relaxation in which viscosity decreases sharply with depth. Oort is less relaxed, but may be somewhat younger and less affected by an early epoch of high heat flow. Finite element calculations show that this heat flow would have to have been substantial to explain Edgeworth's upbowed floor by viscous relaxation, several times steady-state radiogenic values for present-day surface temperatures. We expect Pluto's brittle ice lithosphere to be fractured and porous, however, markedly reducing thermal conductivity and increasing temperatures at depth and relaxation for a given heat flow. We find that most relaxation occurs well within 100 Myr after impact for Edgeworth and Oort, and focus attention on a temporal (and/or regional) epoch of elevated heat flow, possibly tied to the serpentinization of Pluto's rocky core.

Journal of Geophysical Research: Planets↗
USGS2023

Relationship between explosive and effusive volcanism in the Montes Apenninus region of the Moon

Lunar Pyroclastic Deposits (LPDs) are sites of explosive volcanism and often occur in areas of effusive volcanism on the Moon. On Earth, it has been observed that most volcanism has both effusive and explosive phases, whereas on the Moon, these two types of volcanism have typically been considered separately. We hypothesize that the relationship between explosive and effusive volcanism on the Moon is similar to what is observed on the Earth, where individual eruptions can experience multiple phases rather than one type of volcanism always preceding another or occurring separately. We present observations from the Moon Mineralogy Mapper detailing compositional relationships between volcanic features in the lunar Montes Apenninus region. We evaluated whether co-located LPDs and effusive features (e.g., rilles, mare) could have erupted from the same volcanic vent or even at the same time based on their compositional similarities and stratigraphic relationships. We found that the LPDs have varied stratigraphic relationships with co-located effusive features. We identified LPDs near sinuous rilles that may be related to the formation of the rille, where explosive and effusive volcanism occurred at the same vent (e.g., Mozart Rille), and LPDs that may be unrelated to the rille (e.g., Rimae Bode and Rima Bode LPD). Our results suggest that lunar volcanism can mirror terrestrial volcanism, with explosive and effusive eruptions demonstrating more complex dynamics and relationships than previously thought. This variability suggests that the relationship between LPDs and nearby volcanic features cannot be generalized for studies on their resource potential, eruption styles, or deposit volume.

Journal of Geophysical Research: Planets↗
USGS2024

Evolution of a lake margin recorded in the Sutton Island member of the Murray formation, Gale crater, Mars

This study uses data from the Mars Science Laboratory Curiosity rover to document the facies of the Sutton Island member of the Murray formation, interpret paleoenvironments, and establish key stratigraphic transitions at Gale crater. Two facies associations were identified: Facies Association 1 (FA1) and Facies Association 2 (FA2). Individual facies in FA1 include planar-laminated mudstone with minor intervals of planar sandstone, ripple cross-laminated sandstone, cross-stratified sandstone, and alternating laminated sandstone and mudstone. Meter-thick packages of planar-laminated mudstone in FA1 are interpreted to represent deposition in low-energy ponded environments along the lake margin. Straight- and curve-crested ripple cross-laminated facies are interpreted to represent current-influenced deposition. Cross-stratified sandstone facies consist of dm-thick sets that represent deposition in distal channels. Intercalated mm-scale mudstone and sandstone laminae represent waning flow conditions and possible channel abandonment. Facies in FA1 collectively represent deposition in a distal delta plain. FA2 is comprised of planar-laminated mudstone with minor sandstone and is interpreted to represent deposition in a lacustrine-basin setting by suspension settling linked to density flows. FA1 transitions upward into FA2, defining a rapid transgression substantial enough to facilitate the deposition of distal lake facies above delta plain facies. The abrupt transition from FA2 back to FA1 deltaic deposits is suggestive of forced regression. Facies in FA1 and FA2 are consistent with the prevalence of aqueous environments recorded in other Murray formation members and extend our understanding of the dynamic sedimentary processes that characterized ancient lacustrine systems at Gale crater.

Journal of Geophysical Research - Planets↗
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↗
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↗