Search USGS⌕ Search

SEARCH · Search USGS

Results for “JGR - Atmospheres”

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 109 records · Page 6Linked to original sources
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗
USGS2026

Structural evolution and slip rate variations through time of the Puente-Hills blind-thrust fault beneath Los Angeles: Implications for seismic hazard and folding kinematics

Using seismic reflection profiles, historical well logging data, and luminescence and radiocarbon ages, we determine a Pleistocene-Holocene slip history for the central, Santa Fe Springs segment of the Puente Hills blind-thrust fault (PHT), a major seismogenic fault situated beneath the urbanized Los Angeles metropolitan region. We analyze the geometry of correlative stratigraphic units in the forelimb and backlimb of the overlying growth-fold, the Santa Fe Springs anticline, and determine the uplift of seven age-correlative markers. Uplift measurements are converted to thrust displacements on the underlying PHT using a structural method laid out by Don et al. (2022, https://doi.org/10.1785/0120220048 ) that accounts for the geometry of the fault. These data indicate that deep thrust displacement on the PHT is partially consumed updip in the creation of a hanging-wall fault-bend fold, with forelimb growth strata recording <80% of the slip documented within the backlimb. Chronological data from growth strata yield age constraints for folding and faulting on the underlying PHT, providing a detailed incremental slip history derived from both the forelimb and backlimb folding for seven discrete growth horizons, spanning the past 1.4 million years. The resulting six incremental slip rates demonstrate that fault slip has varied through time from the middle Pleistocene to the Holocene. Moreover, these results reveal synchronous acceleration of both the central, Santa Fe Springs and western, Los Angeles segments of the fault system since late Pleistocene time (after 200 ka) and slip rates of greater than 2 mm/yr on the downdip, backlimb, fault-ramp below the anticline.

California↗
USGS2025

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth↗
USGS2026

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗
USGS2026

Compact seismicity bursts have different characteristics from regional seismicity

Earthquakes tend to cluster, developing into sequences driven by stress perturbations and transient fault-zone processes. Depending on the driving process, earthquake sequences show differing behaviors. This variability challenges our ability to observe or distinguish these driving processes in high resolution. Here we systematically identify seismicity bursts throughout southern California using new statistical methods and examine their causes with other independent geophysical observations. Seismicity bursts are defined as statistically significant seismicity rate anomalies. Our approach has the advantage of being data driven, depending on few earthquake occurrence assumptions. We find abundant seismicity bursts across southern California, most frequently occurring along the San Jacinto Fault and in the Salton Sea and Coso geothermal fields. These seismicity bursts are highly compact in space and time, often encompassed by a 5 km radius and have durations less than 10 hr. Many of the seismicity bursts have their largest earthquake near the beginning of the sequence, but the largest earthquake is usually not an obvious mainshock. We utilize a variety of independent geophysical data sets to analyze the characteristics of the seismicity bursts, finding that the seismicity bursts have low b-values, low spectral stress drops, and varied stress ratios compared to regional seismicity. These differences suggest that seismicity bursts are driven by transient processes acting frequently across fault networks.

California↗
USGS2026

Cajon Pass and the southern San Andreas Fault System: Earthquake cycle stress accumulation and present-day loading

With over a century since the last major rupture affecting the wider Los Angeles region, tectonic stress has steadily built along the southern San Andreas and San Jacinto fault systems, raising concerns of an imminent large earthquake. Cajon Pass, located at the junction of these faults, represents a critical site for potential through-going ruptures in Southern California. We constructed new 4D earthquake cycle simulations using a 1000-year paleoseismic rupture history of the San Andreas Fault System (SAFS) to assess spatial and temporal variations in stress. A semi-analytic Fourier transform model was used to compute stress from 3D dislocations in an elastic plate overlying a Maxwell viscoelastic half-space, assuming a complete coseismic reset of resolved shear stress on ruptured elements. Results show highest stress accumulation north of Cajon Pass (∼1.8 MPa/100 years) due to greater slip rates, and lower rates south of Cajon Pass (∼1.0–1.5 MPa/100 years). By 2025, Coulomb stress is estimated at 2.8 MPa on the Mojave South (MOS) segment, 1.8 MPa on the North San Bernardino (NSB1) segment and 3.6 MPa on the San Jacinto Bernardino (SJB) segment. Segments accumulate stresses with characteristic ranges of pre-event stress interpreted as failure thresholds: 1.2–2.7 MPa for MOS, 0.4–1.6 MPa for NSB1, and 1.2–2.9 MPa for SJB. When the stress disparity between segments SJB and MOS narrows, the faults appear to rupture jointly, suggesting that stress levels may control how Cajon Pass acts as an earthquake gate. These results may inform seismic hazard assessments by linking stress evolution to fault interactions.

California↗
USGS2026

Tectonic controls on volcanism and associated hydrothermal activity in a sediment-dominated mid-ocean ridge; Escanaba Trough

Mid-ocean ridges, the Earth's most extensive volcanic system, exhibit unique characteristics in sediment-dominant environments. Thick sediment cover insulates the crust and channels fluid along pathways that can lead to the formation of distinct crustal alteration patterns, exceptionally large mineral deposits, and specialized chemosynthetic ecosystems. This study presents an interdisciplinary investigation into the tectonics of the Escanaba Trough, a heavily sedimented axial valley at the southern Gorda Ridge in the Northeast Pacific Ocean. A primary challenge in such environments is overcoming the masking effect of thick sediments on basement structures that control magmatic and hydrothermal activity. We address this by employing three-dimensional (3D) magnetic modeling of high-resolution near-seafloor magnetic data collected by an autonomous underwater vehicle (AUV). The 2022 surveys with AUV Sentry provided data for 3D magnetic susceptibility models, refining our understanding of the geometry of sub-sediment laccoliths/saucer-shaped sills and hydrothermal alteration. In conjunction with a new 1:100,000 scale lithostratigraphic map, we outline the tectonic controls on the emplacement of Escanaba Trough's three main volcanic centers, characterize the geometry of its spreading segments, and provide volumetric data on the distribution of sub-sediment volcanism in the southern Gorda Ridge.

Escanaba Trough, Pacific Ocean↗
USGS2026

Fully contained laboratory earthquakes: The effect of asperity aspect ratio and free surfaces

Corner frequency ( f c ) and seismic moment ( M 0 ) are key parameters derived from seismic signals that are used to characterize earthquake stress drop, rupture area, and slip. These parameters are also affected by fault geometry and boundary conditions. However, the systematic study of these effects in laboratory settings has been challenging. This study presents laboratory earthquake experiments that examine how rupture dynamics are influenced by (a) the aspect ratio of rectangular PMMA velocity-weakening (VW) asperities surrounded by the Teflon velocity-strengthening (VS) patches, and (b) whether the sides of a VW asperity are confined with VS patches or are free surfaces. We found that increasing confinement by reducing free surfaces or increasing the VW asperity aspect ratio stabilizes fault slip, so that higher normal stress is required to transition from aseismic to seismic slip. Increased confinement and high aspect ratios also reduced M 0 and increased f c , both of which were determined from the radiated seismic waves. M 0 and f c were primarily controlled by the shorter dimension of the VW asperity. Analysis of high-frequency acoustic emission signals revealed that ruptures on high-aspect-ratio VW asperities propagated more unidirectionally, whereas ruptures on square VW asperities were more complex. Further, the high-aspect-ratio asperities were more likely to be eroded by surrounding VS regions while low-aspect-ratio asperities were more likely to rupture into the VS surroundings. These results demonstrate that both the confinement from surrounding stable areas and the geometry of the seismogenic patch can affect rupture nucleation, propagation, and seismic source characteristics.

JGR Solid Earth↗
USGS2025

Wind pumping dominates landward salt transport in a weakly tidal estuary

In tidally energetic estuaries, salinity dynamics vary with tidal forcing as well as factors such as river discharge and bathymetry. However, in weakly tidal estuaries, mechanisms governing salt transport remain poorly understood. Here, we investigate salt transport processes in a weakly tidal estuary, Albemarle Sound on the U.S. East Coast, using a deterministic numerical model. In 2022, the volume averaged salinity of Albemarle Sound ranged from ∼2 to 10 psu. Although seasonal river discharge typically induces gradual salinity changes of less than 1 psu per month, four wind events caused sharp salinity increases of 1–2 psu within just a few days. We find that wind pumping, defined as the net salt transport due to temporal correlations between velocity and salinity associated with wind events, is the dominant mechanism contributing to these episodic increases in volume averaged salinity. In Albemarle Sound, wind pumping operates on timescales of several days, which is longer than the tidal timescale, and occurs intermittently throughout the year rather than following a regular repetitive pattern. Our analysis suggests that salt transport by wind pumping is typically an order of magnitude larger than by gravitational circulation and serves as the primary mechanism balancing salt loss from river flow. Isolating wind pumping requires averaging over variability of multiple days, longer than typical tidal averaging. Wind pumping likely also controls material exchange with coastal ocean, influencing biogeochemical processes in the estuary. This mechanism may be important in other systems with similar characteristics, particularly in coastal lagoons with multiple subembayments.

North Carolina↗
USGS2026

Short-term estuarine phytoplankton dynamics in response to hurricanes along the Gulf Coast of America: A Variational Autoencoder (VAE) approach with satellite and bio-optical observations

Hurricanes drive diverse estuarine phytoplankton responses and can trigger cascading ecological and physicochemical impacts. Capturing these short-term dynamics requires high spatiotemporal resolution. Here, we applied a globally-applicable coastal ocean color algorithm, Variational Autoencoder (VAE), to Sentinel-2 MSI imagery for chlorophyll- a (Chl- a ) estimation and validated its strong performance across the northern Gulf coast of America (GoA) estuaries, including Galveston Bay (TX), Barataria-Terrebonne Estuary (LA), Apalachicola Estuary (FL) and Tampa Bay (FL). The test set showed strong performance (MAE: 1.44 mg m −3 ; RMSE: 17.7 mg m −3 ; slope: 0.86; median symmetric accuracy: 30.33%). The validated VAE was then applied to 76 Sentinel-2 MSI images to assess phytoplankton biomass responses to hurricanes Harvey (2017), Michael (2018), Ida (2021), Francine (2024), Helene (2024), and Milton (2024) in the GoA estuaries. Results showed that hurricane disturbances on Chl- a typically lasted 3–5 weeks. Estuarine waters west (left) of hurricane tracks showed a rapid decline in Chl- a (∼5 mg m −3 ) due to elevated turbidity from heavy rainfall, and wind-driven flushing in the estuary, followed by a rebound over about two weeks, with Chl- a increasing approximately 10–15 mg m −3 above pre-storm levels. In contrast, right-side waters showed a slower response, likely from oligotrophic seawater intrusion driven by the hurricane's counterclockwise rotation. Post-storm observations showed increased freshwater phytoplankton like chlorophytes and cyanobacteria dominating estuaries, while shelf-waters exhibited elevated dinoflagellates (e.g., Karenia brevis bloom after Hurricane Milton). These results highlight the spatial heterogeneity of hurricane impacts on estuarine phytoplankton dynamics, which may trigger cascading effects on biogeochemical cycling and food webs, potentially prolonging ecosystem recovery.

Alabama, Florida, Louisiana, Mississippi, Texas↗
USGS2025

Numerical simulation of sound-side barrier-island inundation and breaching during Hurricane Dorian (2019)

Hurricane-induced morphological changes and associated community hazards along sandy, barrier-island coastlines have been studied primarily from the perspective of ocean-side attack by storm-driven ocean surge and large waves. Thus, our understanding of long-term barrier island morphological change focuses on beach erosion, overwash, and inlet formation. In contrast, outwash events with inundation from the sound side, such as one that occurred in Cape Lookout National Seashore, North Carolina, USA during Hurricane Dorian (September 2019), are understudied. Studying such events can improve understanding of barrier island response and stability for a broader range of conditions. Here, we model the hydrodynamics and morphological evolution of a barrier island using a coupled wave-current-sediment transport modeling system. Wind-driven surge in Pamlico Sound led to overtopping from the sound side, which eroded outwash channels and transported sediment seaward into the nearshore. Simulations reproduce the channel features observed with aerial imagery and provide information not available from the remote-sensing observations, including channel depths (>2 m) and the fate of the eroded sand. We found that >99% of the eroded sand was deposited in the nearshore, within 1,000 m of the shoreline in depths <10 m, suggesting that the deposited sediment remains available for littoral transport and beach recovery. Simulations with combinations of coarse or fine sediment and vegetated or unvegetated landcover indicate that channel position did not vary with grain size or vegetation, while volume of erosion and channel morphology were more responsive to variations in grain size and less responsive to presence of vegetation.

North Carolina↗
USGS2025

Modeling the influence of upper and lower shoreface dynamics on barrier island evolution

Barrier island resilience to climate impacts depends on sediment redistribution between the subaqueous shoreface and subaerial barrier during sea-level rise and storms. However, autogenic interactions between the upper and lower shoreface and their influence on the subaerial barrier are poorly characterized. Here, we explore the influences of various shoreface components on barrier morphology using a model of barrier and shoreface evolution under sea-level rise, the Articulated Barrier Shoreface (ABSF) Model. This reduced-complexity model divides the shoreface into upper and lower shoreface panels that respond independently to sea-level rise and deviations from the equilibrium slope. We couple the ABSF with the Lorenzo-Trueba & Ashton, 2014, https://doi.org/10.1002/2013jf002941 model (LTA), a barrier island evolution model driven by overwash and sea-level rise. Through this coupled framework, we examine the influences of upper and lower shoreface slopes, their respective depths, and sensitivity to wave climate on long-term barrier evolution. Results show that the relative depths of the upper and lower shoreface toes influence barrier response to rising seas, alongside overwash flux and closure depth. Notably, the lower shoreface response to sea-level change lags that of the upper shoreface over decades, diminishing the resilience of the barrier over centennial timescales by slowing the overall barrier response. In fact, the ABSF model predicts barriers will drown faster and more than predicted with a linear shoreface. Results highlight the shoreface as an important sediment reservoir for barrier islands and that differences in upper and lower shoreface responses can reduce barrier resilience to sea-level rise due to limited lower shoreface sediment accessibility.

JGR Earth Surface↗
USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗
USGS2025

Are equilibrium shoreline models just convolutions?

Yes. Equilibrium shoreline models, which simulate wave-driven cross-shore erosion and accretion, are mathematically equivalent to a discrete convolution (i.e., a weighted, moving average) of a time series of wave-forcing conditions with a parameterized memory-decay kernel function. The direct equivalence between equilibrium shoreline models and convolutions reveals key theoretical aspects of equilibrium behavior. Convolutions (representing quasi-low-pass filter operations) provide an intuitive theoretical description of shoreline erosion and accretion behavior in response to waves: that is, shoreline position often mirrors the weighted moving average of wave time series. Model-convolution equivalence also provides a conceptual basis to interpret, evaluate, and construct data-driven Machine-Learning/Deep-Learning (ML/DL) models that use convolutions to extract features from data and then apply them for prediction (e.g., Convolutional Neural Networks (CNNs)). Finally, our findings provide a methodological pathway (based on Fourier transforms) for future understanding of wave-driven shoreline change, which can be used to interpret the coherence between the frequency spectrum of the processes of waves and shoreline change and construct more computationally efficient and effective shoreline-modeling approaches.

JGR Earth Surface↗
USGS2026

Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico

Landscapes in humid climates often experience rapid chemical and mechanical weathering that can change groundwater fluxes, drainage pathways, and frictional strength via alterations in subsurface mineral composition and matrix porosity. These changes to hillslope critical zone structure and strength should control landslide susceptibility, size, and failure style. Here we test the impact of differential weathering on landslide failure processes by exploiting a natural experiment along the dioritic Utuado batholith in Puerto Rico's Cordillera Central, where a sharp weathering gradient follows an upstream-propagating knickzone into an uplifted shoreline platform that was affected by Hurricane Maria in 2017. Through a combination of landslide mapping, field geophysical and geotechnical surveys, topographic analysis, calibrated stream profile inversion, 3D back-analysis of hillslope strength, and debris-flow inundation assessments, we show that differential weathering induces changes to critical zone structure and consequently landslide susceptibility and hazard. Our mapping reveals that the hurricane-induced landslides preferentially cluster along less weathered hillslopes despite a marked increase in frictional strength. Jointly inverted seismic refraction and electrical resistivity profiles show that the frictionally weak, heavily weathered hillslopes are highly porous and dry, which we posit acts to stabilize slopes through limiting deeper hydrologic connectivity that can produce shallow water tables with destabilizing lateral seepage forces during storms. Debris-flow inundation mapping and analysis reveals that major contributors to debris-flow growth, channel steepness and landslide coalescence, scale linearly with hillslope friction angle. This leads to higher-mobility debris flows along knickzone reaches, which disproportionately affects communities at the base of the plateau.

Puerto Rico↗
USGS2025

Evidence of mineral alteration in a salt marsh subterranean estuary: Implications for carbon and trace element cycling

Subterranean estuaries (STE) in salt marshes are biogeochemically active zones where interactions between terrestrial groundwater and seawater drive complex cycling of carbon and trace elements, influenced by mineral dissolution. These systems, characterized by fine-grained organic-rich peat overlying permeable coastal aquifers, play a crucial role as a blue carbon sink, yet their geochemical dynamics remain poorly understood. We investigated dissolved trace elements, carbon, silica, and radium isotopes in a salt marsh STE (Sage Lot Pond, Waquoit Bay, MA) over seasonal and annual cycles. Our results reveal that groundwater and estuarine water circulation through marsh peat and aquifer sediments leads to enrichments of dissolved organic and inorganic carbon (DOC and DIC), Si, Ba, and Mn, with variable source/sink behavior of Fe and net removal of U. Submarine groundwater discharge dominated Ba fluxes, whereas pore water drainage from marsh peat acted as the main sink for U and source of Si. Fe cycling was variable, with terrestrial Fe largely removed as groundwater passed through the STE, consistent with Fe-sulfide and amorphous phase formation. Radium isotope ratios identified two distinct subsurface flow pathways, influenced by metal-oxide cycling and organic matter breakdown. Si production was decoupled from DIC, suggesting Si originates from mineral alteration, whereas DIC results from both mineral weathering and microbial respiration. Silicate mineral alteration, coupled with marsh pore water drainage, accounts for up to 16% of annual DIC exports (66 g C m −2 y −1 ), highlighting the importance of STEs in coastal carbon and trace element cycling, especially as marshes face environmental change.

Massachusetts↗
USGS2026

Biophysical controls on sediment erodibility in shallow estuarine embayments

The erodibility of cohesive sediment is known to vary both spatially and temporally but the factors governing its variation are not well understood. We conducted a field investigation of the influence of hydrodynamic forcing, sediment properties, and benthic infauna on erodibility in the muddy shallows of San Pablo and Grizzly Bays in northern San Francisco Bay in summer 2019 and winter 2020. An erosion rate parameter M c was determined from regressions between near-bed vertical turbulent sediment flux, as a proxy for erosion, and bed shear stress due to currents. During each 2-month study period, we measured benthic infauna abundance and dry bulk density, particle size distribution, percent organic carbon, chlorophyll a , pheophytin a , and carbohydrates carbon concentrations of surficial bed sediments five or six times. M c increased with bed shear stress due to waves in both embayments. In San Pablo Bay, erodibility was approximately 50% lower during the winter than the summer deployment, whereas in Grizzly Bay, there was no significant difference. The factor most strongly related to the decrease in M c in San Pablo Bay was increased abundance of the amphipod Ampelisca abdita . The observed reduction in erodibility may occur in many muddy estuaries because A. abdita is broadly distributed in the coastal waters of North America. Erodibility was also directly related to biomass of the invasive clam Potamocorbula amurensis . Erodibility did not depend on dry bulk density: bulk density did not vary seasonally in San Pablo Bay and was lower in winter than summer in Grizzly Bay.

California↗
USGS2025

Wetlands, groundwater and seasonality influence the spatial distribution of stream chemistry in a low-relief catchment

Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well-established for high-relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low-relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low-relief wetland-rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low-relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment-scale studies for low-relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low-relief wetland-rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes.

Michigan↗