USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ocean Modeling”

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 487 records · Page 27Linked to original sources
USGS2013

Insights and issues with simulating terrestrial DOC loading of Arctic river networks

Terrestrial carbon dynamics influence the contribution of dissolved organic carbon (DOC) to river networks in addition to hydrology. In this study, we use a biogeochemical process model to simulate the lateral transfer of DOC from land to the Arctic Ocean via riverine transport. We estimate that, over the 20th century, the pan-Arctic watershed has contributed, on average, 32 Tg C/yr of DOC to river networks emptying into the Arctic Ocean with most of the DOC coming from the extensive area of boreal deciduous needle-leaved forests and forested wetlands in Eurasian watersheds. We also estimate that the rate of terrestrial DOC loading has been increasing by 0.037 Tg C/yr2 over the 20th century primarily as a result of climate-induced increases in water yield. These increases have been offset by decreases in terrestrial DOC loading caused by wildfires. Other environmental factors (CO2 fertilization, ozone pollution, atmospheric nitrogen deposition, timber harvest, agriculture) are estimated to have relatively small effects on terrestrial DOC loading to Arctic rivers. The effects of the various environmental factors on terrestrial carbon dynamics have both offset and enhanced concurrent effects on hydrology to influence terrestrial DOC loading and may be changing the relative importance of terrestrial carbon dynamics on this carbon flux. Improvements in simulating terrestrial DOC loading to pan-Arctic rivers in the future will require better information on the production and consumption of DOC within the soil profile, the transfer of DOC from land to headwater streams, the spatial distribution of precipitation and its temporal trends, carbon dynamics of larch-dominated ecosystems in eastern Siberia, and the role of industrial organic effluents on carbon budgets of rivers in western Russia.

Ecological Applications↗
USGS2010

Paleoclimates: Understanding climate change past and present

The field of paleoclimatology relies on physical, chemical, and biological proxies of past climate changes that have been preserved in natural archives such as glacial ice, tree rings, sediments, corals, and speleothems. Paleoclimate archives obtained through field investigations, ocean sediment coring expeditions, ice sheet coring programs, and other projects allow scientists to reconstruct climate change over much of earth's history. When combined with computer model simulations, paleoclimatic reconstructions are used to test hypotheses about the causes of climatic change, such as greenhouse gases, solar variability, earth's orbital variations, and hydrological, oceanic, and tectonic processes. This book is a comprehensive, state-of-the art synthesis of paleoclimate research covering all geological timescales, emphasizing topics that shed light on modern trends in the earth's climate. Thomas M. Cronin discusses recent discoveries about past periods of global warmth, changes in atmospheric greenhouse gas concentrations, abrupt climate and sea-level change, natural temperature variability, and other topics directly relevant to controversies over the causes and impacts of climate change. This text is geared toward advanced undergraduate and graduate students and researchers in geology, geography, biology, glaciology, oceanography, atmospheric sciences, and climate modeling, fields that contribute to paleoclimatology. This volume can also serve as a reference for those requiring a general background on natural climate variability.

Book↗
USGS2016

Processes of multibathyal aragonite undersaturation in the Arctic Ocean

During 3 years of study (2010–2012), the western Arctic Ocean was found to have unique aragonite saturation profiles with up to three distinct aragonite undersaturation zones. This complexity is produced as inflow of Atlantic-derived and Pacific-derived water masses mix with Arctic-derived waters, which are further modified by physiochemical and biological processes. The shallowest aragonite undersaturation zone, from the surface to ∼30 m depth is characterized by relatively low alkalinity and other dissolved ions. Besides local influence of biological processes on aragonite undersaturation of shallow coastal waters, the nature of this zone is consistent with dilution by sea-ice melt and invasion of anthropogenic CO 2 from the atmosphere. A second undersaturated zone at ∼90–220 m depth (salinity ∼31.8–35.4) occurs within the Arctic Halocline and is characterized by elevated p CO 2 and nutrients. The nature of this horizon is consistent with remineralization of organic matter on shallow continental shelves bordering the Canada Basin and the input of the nutrients and CO 2 entrained by currents from the Pacific Inlet. Finally, the deepest aragonite undersaturation zone is at greater than 2000 m depth and is controlled by similar processes as deep aragonite saturation horizons in the Atlantic and Pacific Oceans. The comparatively shallow depth of this deepest aragonite saturation horizon in the Arctic is maintained by relatively low temperatures, and stable chemical composition. Understanding the mechanisms controlling the distribution of these aragonite undersaturation zones, and the time scales over which they operate will be crucial to refine predictive models.

Journal of Geophysical Research: Oceans↗
USGS1987

Subsidence, crustal structure, and thermal evolution of Georges Bank basin

A geophysical study of Georges Bank basin defines a deep crustal structure that is interpreted in terms of the basin's tectonic and thermal history. Gravity models along three basin cross sections delineate two zones of crustal thinning at the basement hinge zone and oceanic crustal margins. These two zones bound rift-stage crust (about 25 km thick) which underlies the central portion of the basin. Subsidence analysis of the basin, using data from multichannel seismic reflection lines and two COST wells, suggests a rifting and (uniform) extensional origin. Two-dimensional finite difference modeling of the basin defines a crustal structure that concurs with the gravity and subsidence studies. The resulting isotherms show no major changes in the thermal structure since the ate Jurassic. In some areas of the basin, temperatures sufficient for oil generation are determined from maturation studies of Jurassic sediments. Hydrocarbon generation is questionable, however, because of the probable lack of proper and sufficient kerogen in the Jurassic deposits.

American Association of Petroleum Geologists Bulle↗
USGS2008

Important observations and parameters for a salt water intrusion model

Sensitivity analysis with a density-dependent ground water flow simulator can provide insight and understanding of salt water intrusion calibration problems far beyond what is possible through intuitive analysis alone. Five simple experimental simulations presented here demonstrate this point. Results show that dispersivity is a very important parameter for reproducing a steady-state distribution of hydraulic head, salinity, and flow in the transition zone between fresh water and salt water in a coastal aquifer system. When estimating dispersivity, the following conclusions can be drawn about the data types and locations considered. (1) The "toe" of the transition zone is the most effective location for hydraulic head and salinity observations. (2) Areas near the coastline where submarine ground water discharge occurs are the most effective locations for flow observations. (3) Salinity observations are more effective than hydraulic head observations. (4) The importance of flow observations aligned perpendicular to the shoreline varies dramatically depending on distance seaward from the shoreline. Extreme parameter correlation can prohibit unique estimation of permeability parameters such as hydraulic conductivity and flow parameters such as recharge in a density-dependent ground water flow model when using hydraulic head and salinity observations. Adding flow observations perpendicular to the shoreline in areas where ground water is exchanged with the ocean body can reduce the correlation, potentially resulting in unique estimates of these parameter values. Results are expected to be directly applicable to many complex situations, and have implications for model development whether or not formal optimization methods are used in model calibration.

Ground Water↗
USGS2010

The PRISM3D paleoenvironmental reconstruction

The Pliocene Research, Interpretation and Synoptic Mapping (PRISM) paleoenvironmental reconstruction is an internally consistent and comprehensive global synthesis of a past interval of relatively warm and stable climate. It is regularly used in model studies that aim to better understand Pliocene climate, to improve model performance in future climate scenarios, and to distinguish model-dependent climate effects. The PRISM reconstruction is constantly evolving in order to incorporate additional geographic sites and environmental parameters, and is continuously refined by independent research findings. The new PRISM three dimensional (3D) reconstruction differs from previous PRISM reconstructions in that it includes a subsurface ocean temperature reconstruction, integrates geochemical sea surface temperature proxies to supplement the faunal-based temperature estimates, and uses numerical models for the first time to augment fossil data. Here we describe the components of PRISM3D and describe new findings specific to the new reconstruction. Highlights of the new PRISM3D reconstruction include removal of Hudson Bay and the Great Lakes and creation of open waterways in locations where the current bedrock elevation is less than 25m above modern sea level, due to the removal of the West Antarctic Ice Sheet and the reduction of the East Antarctic Ice Sheet. The mid-Piacenzian oceans were characterized by a reduced east-west temperature gradient in the equatorial Pacific, but PRISM3D data do not imply permanent El Niño conditions. The reduced equator-to-pole temperature gradient that characterized previous PRISM reconstructions is supported by significant displacement of vegetation belts toward the poles, is extended into the Arctic Ocean, and is confirmed by multiple proxies in PRISM3D. Arctic warmth coupled with increased dryness suggests the formation of warm and salty paleo North Atlantic Deep Water (NADW) and a more vigorous thermohaline circulation system that may have provided the enhanced ocean heat transport necessary to move warm surface water to the Arctic. New deep ocean temperature data also suggests greater warmth and further southward penetration of paleo NADW.

Stratigraphy↗
USGS2012

A circulation modeling approach for evaluating the conditions for shoreline instabilities

Analytical models predict the growth (instability) of shoreline salients when deep-water waves approach the coast from highly oblique angles, contrary to classical shoreline change models in which shoreline salients can only dissipate. Using the process-based wave, circulation, and sediment transport model Delft3D, we test this prediction for simulated bathymetric and wave characteristics approximating the open-ocean conditions at Duck, North Carolina. We consider two cases: a uniform coast with a varying wave approach angle, and a bathymetry with coastal salients and a single high-angle boundary wave condition. Incident wave conditions include a swell case with no wind and a wind-wave case with active local wave regeneration by wind. The uniform-coast tests predict transport maxima at oblique wave angles for both wave cases, indicating the potential for shoreline instabilities, similar to the analytical models. However, the critical angle for instability is much higher in the wind-wave case. Our tests with coastal salients agree with previous findings that a minimum salient length scale may be required for the instability effect to be active. Here, a salient with a longshore scale of 4 km results in transport divergence (erosion; no instability) at the salient crest while an 8 km salient results in transport convergence (accretion; instability) at the crest.

North Carolina↗
USGS2022

Data-driven modeling of wind waves in upper Delaware Bay with living shorelines

Living shoreline projects have been built to preserve coastal ecosystems under future climate change and sea level rise. To quantify the wave power variation across living shorelines, the wave characteristics around the constructed oyster reefs (CORs) in upper Delaware Bay were investigated in this study. Wave parameters seaward and shoreward of CORs were recorded by wave gauges in early 2018. Four winter storms happened in this period and induced strong winds and coastal flooding at the study site. To estimate the wind wave characteristics across the CORs on a yearly basis, soft computing-based models combining fully connected neural networks and long short-term memory were developed to extend the two-month energetic wave measurements. The results show that when CORs were emergent or slightly submerged, the averaged wave height attenuation was about 39.8% from the offshore gauge to the nearshore gauge (behind CORs) during 2018–2020, owing to the combined effect of nearshore bathymetric changes and CORs. Furthermore, it was found that the annually averaged wave power reduction from offshore to nearshore at the study site was about 30.0% in 2018, 2019, and 2020. This study provides a novel framework to predict long-term wave characteristics based on short-term wave measurements using soft computing-based models.

New Jersey↗
USGS2024

Physics-based satellite-derived bathymetry (SDB) using Landsat OLI images

The estimation of depth in optically shallow waters using satellite imagery can be efficient and cost-effective. Active sensors measure the distance traveled by an emitted laser pulse propagating through the water with high precision and accuracy if the bottom peak intensity of the waveform is greater than the noise level. However, passive optical imaging of optically shallow water involves measuring the radiance after the sunlight undergoes downward attenuation on the way to the sea floor, and the reflected light is then attenuated while moving back upward to the water surface. The difficulty of satellite-derived bathymetry (SDB) arises from the fact that the measured radiance is a result of a complex association of physical elements, mainly the optical properties of the water, bottom reflectance, and depth. In this research, we attempt to apply physics-based algorithms to solve this complex problem as accurately as possible to overcome the limitation of having only a few known values from a multispectral sensor. Major analysis components are atmospheric correction, the estimation of water optical properties from optically deep water, and the optimization of bottom reflectance as well as the water depth. Specular reflection of the sky radiance from the water surface is modeled in addition to the typical atmospheric correction. The physical modeling of optically dominant components such as dissolved organic matter, phytoplankton, and suspended particulates allows the inversion of water attenuation coefficients from optically deep pixels. The atmospheric correction and water attenuation results are used in the ocean optical reflectance equation to solve for the bottom reflectance and water depth. At each stage of the solution, physics-based models and a physically valid, constrained Levenberg–Marquardt numerical optimization technique are used. The physics-based algorithm is applied to Landsat Operational Land Imager (OLI) imagery over the shallow coastal zone of Guam, Key West, and Puerto Rico. The SDB depths are compared to airborne lidar depths, and the root mean squared error (RMSE) is mostly less than 2 m over water as deep as 30 m. As the initial choice of bottom reflectance is critical, along with the bottom reflectance library, we describe a pure bottom unmixing method based on eigenvector analysis to estimate unknown site-specific bottom reflectance.

Remote Sensing↗
USGS2024

Spatio-temporal distribution of adult Pacific lamprey Entosphenus tridentatus relative to habitat fragmentation

Pacific lamprey ( Entosphenus tridentatus ), a fish species native to the Pacific Northwest (USA), have distinctive cultural and ecological value but determining their spatial and temporal distribution is challenging due to a general lack systematic monitoring. In this study, we used counts of Pacific lamprey redds to model the probability of occurrence and abundance of Pacific lamprey based on environmental covariates including artificial barriers, assuming higher predicted lamprey redds translates to more suitable spawning habitats. Using generalized linear mixed zero-inflated models, results suggest that Pacific lamprey abundance was generally lower in high gradient streams, further from the ocean. Stream reaches with warmer spring water temperatures and greater historical median spring flows supported higher abundances. Lamprey occurrence was primarily influenced by spring water temperatures and distance from the ocean. We further observed that when streams warm beyond 18°C, confidence intervals around the abundance estimates widen and zero-inflation increases, indicating a decrease in occurrence. One objective of the study was to recommend where barrier removal or restoration should be prioritized to increase passage and thus access to upstream habitats. We considered artificial barriers to primarily influence the probability of occurrence through access. The barrier variable in this model had a negative effect on the probability of lamprey occurrence, but it was not a strong predictor in the model. While we are not able to suggest specific locations that would most benefit barrier removal or improvement based on these model results, we can identify the watersheds with a higher probability to support Pacific lamprey and provide potential additional habitats by improving habitat connectivity. Focusing restoration and/ or removal of barriers on watersheds in the Mid-South region of the Oregon Coast (i.e., Alsea, Siuslaw, Coos, Coquille, and Sixes rivers) with higher habitat suitability could prioritize use of limited funds, increase the probability of benefiting Pacific lamprey, and potentially other native lampreys and migratory (e.g., salmon, steelhead; Oncorhynchus ) species. Although this manuscript focuses on the Oregon Coast region, the methods are transferrable to other regions where Pacific lamprey are present.

Oregon↗
USGS2021

Sea state from single optical images: A methodology to derive wind-generated ocean waves from cameras, drones and satellites

Sea state is a key variable in ocean and coastal dynamics. The sea state is either sparsely measured by wave buoys and satellites or modelled over large scales. Only a few attempts have been devoted to sea state measurements covering a large domain; in particular its estimation from optical images. With optical technologies becoming omnipresent, optical images offer incomparable spatial resolution from diverse sensors such as shore-based cameras, airborne drones (unmanned aerial vehicles/UAVs), or satellites. Here, we present a standalone methodology to derive the water surface elevation anomaly induced by wind-generated ocean waves from optical imagery. The methodology was tested on drone and satellite images and compared against ground truth. The results show a clear dependence on the relative azimuth view angle in relation to the wave crest. A simple correction is proposed to overcome this bias. Overall, the presented methodology offers a practical way of estimating ocean waves for a wide range of applications.

Remote Sensing↗
USGS2011

Modeling sulfate reduction in methane hydrate-bearing continental margin sediments: Does a sulfate-methane transition require anaerobic oxidation of methane?

The sulfate‐methane transition (SMT), a biogeochemical zone where sulfate and methane are metabolized, is commonly observed at shallow depths (1–30 mbsf) in methane‐bearing marine sediments. Two processes consume sulfate at and above the SMT, anaerobic oxidation of methane (AOM) and organoclastic sulfate reduction (OSR). Differentiating the relative contribution of each process is critical to estimate methane flux into the SMT, which, in turn, is necessary to predict deeper occurrences of gas hydrates in continental margin sediments. To evaluate the relative importance of these two sulfate reduction pathways, we developed a diagenetic model to compute the pore water concentrations of sulfate, methane, and dissolved inorganic carbon (DIC). By separately tracking DIC containing 12 C and 13 C, the model also computes δ 13 C‐DIC values. The model reproduces common observations from methane‐rich sediments: a well‐defined SMT with no methane above and no sulfate below and a δ 13 C‐DIC minimum at the SMT. The model also highlights the role of upward diffusing 13 C‐enriched DIC in contributing to the carbon isotope mass balance of DIC. A combination of OSR and AOM, each consuming similar amounts of sulfate, matches observations from Site U1325 (Integrated Ocean Drilling Program Expedition 311, northern Cascadia margin). Without AOM, methane diffuses above the SMT, which contradicts existing field data. The modeling results are generalized with a dimensional analysis to the range of SMT depths and sedimentation rates typical of continental margins. The modeling shows that AOM must be active to establish an SMT wherein methane is quantitatively consumed and the δ 13 C‐DIC minimum occurs. The presence of an SMT generally requires active AOM.

Geochemistry, Geophysics, Geosystems↗
USGS2006

Bottom instrumented tripods: History, applications, and impacts

Instrumented bottom tripods have provided important data on sediment transport processes on continental shelves and in estuaries for four decades. Since the initial deployment in a tidal channel in Puget Sound, WA, in 1965 numerous tripods have been constructed to investigate bottom boundary layer and sediment dynamics worldwide. Tripod data have led to new understanding of near-bottom wave and current flows in the coastal ocean, and have been crucial to the development of shelf circulation and sediment transport models. Calculations of bottom stress, bottom roughness, and sediment flux that resulted directly from tripod data have been compared to bottom boundary layer model results. Where these have differed, new or revised model components have been developed to improve the skill of the models. The many discoveries that have been made from tripod experiments include dense, near-bottom fluid mud layers that transport large quantities of suspended sediment offshore into deeper regions of the continental shelf. This process has been linked to the seaward progradation of subaqueous deltas and to the boundaries of mid-shelf mud deposits off rivers with high fine-sediment discharge.

Continental Shelf Research↗
USGS2026

Tracing mercury from land to river: Global sources, retention, and implications for sustainability

Mercury (Hg) pollution in river systems is a global sustainability challenge. Yet the transport, transformation, and retention of Hg within global rivers remain poorly quantified, particularly in regions with sparse observations such as Southeast Asia and Africa, hindering effective pollution mitigation and reinforcing geographic inequities in scientific knowledge and environmental governance. Here, we present the first global, high-resolution simulation of riverine Hg dynamics using a process-based model that traces Hg from land-based sources through river networks to the ocean. Under a realistic scenario, we estimate that ∼1900 megagrams per year (Mg/yr) of Hg enters global rivers, including 1500 Mg/yr from human-induced sources and 400 Mg/yr from soil erosion. Nearly half of this flux (∼1000 Mg/yr) is retained in reservoirs and dams, which act as major sinks. While such retention limits downstream delivery to the oceans, it also heightens in-reservoir Hg methylation risks. By bridging the gap between Hg releases and observed riverine exports, our framework offers a scalable tool for data-limited regions, promotes data access, and supports global freshwater and pollution-management strategies.

Environmental Science & Technology↗
USGS2015

Precise interpolar phasing of abrupt climate change during the last ice age

The last glacial period exhibited abrupt Dansgaard–Oeschger climatic oscillations, evidence of which is preserved in a variety of Northern Hemisphere palaeoclimate archives 1 . Ice cores show that Antarctica cooled during the warm phases of the Greenland Dansgaard–Oeschger cycle and vice versa 2 , 3 , suggesting an interhemispheric redistribution of heat through a mechanism called the bipolar seesaw 4 , 5 , 6 . Variations in the Atlantic meridional overturning circulation (AMOC) strength are thought to have been important, but much uncertainty remains regarding the dynamics and trigger of these abrupt events 7 , 8 , 9 . Key information is contained in the relative phasing of hemispheric climate variations, yet the large, poorly constrained difference between gas age and ice age and the relatively low resolution of methane records from Antarctic ice cores have so far precluded methane-based synchronization at the required sub-centennial precision 2 , 3 , 10 . Here we use a recently drilled high-accumulation Antarctic ice core to show that, on average, abrupt Greenland warming leads the corresponding Antarctic cooling onset by 218 ± 92 years (2 σ ) for Dansgaard–Oeschger events, including the Bølling event; Greenland cooling leads the corresponding onset of Antarctic warming by 208 ± 96 years. Our results demonstrate a north-to-south directionality of the abrupt climatic signal, which is propagated to the Southern Hemisphere high latitudes by oceanic rather than atmospheric processes. The similar interpolar phasing of warming and cooling transitions suggests that the transfer time of the climatic signal is independent of the AMOC background state. Our findings confirm a central role for ocean circulation in the bipolar seesaw and provide clear criteria for assessing hypotheses and model simulations of Dansgaard–Oeschger dynamics.

Nature↗
USGS2010

Analysis of the Arctic system for freshwater cycle intensification: Observations and expectations

Hydrologic cycle intensification is an expected manifestation of a warming climate. Although positive trends in several global average quantities have been reported, no previous studies have documented broad intensification across elements of the Arctic freshwater cycle (FWC). In this study, the authors examine the character and quantitative significance of changes in annual precipitation, evapotranspiration, and river discharge across the terrestrial pan-Arctic over the past several decades from observations and a suite of coupled general circulation models (GCMs). Trends in freshwater flux and storage derived from observations across the Arctic Ocean and surrounding seas are also described. With few exceptions, precipitation, evapotranspiration, and river discharge fluxes from observations and the GCMs exhibit positive trends. Significant positive trends above the 90% confidence level, however, are not present for all of the observations. Greater confidence in the GCM trends arises through lower interannual variability relative to trend magnitude. Put another way, intrinsic variability in the observations tends to limit confidence in trend robustness. Ocean fluxes are less certain, primarily because of the lack of long-term observations. Where available, salinity and volume flux data suggest some decrease in saltwater inflow to the Barents Sea (i.e., a decrease in freshwater outflow) in recent decades. A decline in freshwater storage across the central Arctic Ocean and suggestions that large-scale circulation plays a dominant role in freshwater trends raise questions as to whether Arctic Ocean freshwater flows are intensifying. Although oceanic fluxes of freshwater are highly variable and consistent trends are difficult to verify, the other components of the Arctic FWC do show consistent positive trends over recent decades. The broad-scale increases provide evidence that the Arctic FWC is experiencing intensification. Efforts that aim to develop an adequate observation system are needed to reduce uncertainties and to detect and document ongoing changes in all system components for further evidence of Arctic FWC intensification.

Journal of Climate↗
USGS2022

Waveform signatures of earthquakes located close to the subducted Gorda Plate interface

Complex seismic velocity structure near the earthquake source can affect rupture dynamics and strongly modify the seismic waveforms recorded near the fault. Fault‐zone waves are commonly observed in continental crustal settings but are less clear in subduction zones due to the spatial separation between seismic stations and the plate boundary fault. We observed anomalously long duration S waves from earthquake clusters located near the interface of the subducted Gorda plate north of the Mendocino triple junction. In contrast, earthquakes located just a few kilometers below each cluster show impulsive S waves. A nodal array experiment was conducted around the Northern California Seismic Network station KCT for two months to investigate the origin of the complex S waves. Beamforming analysis shows that the S waves contain three arrivals that have different horizontal slownesses, which we term S 1, S 2, and S coda. Similar analysis on P waves also show two arrivals with different horizontal slownesses, which we term P 1 and P 2. P 1 and S 1 have larger horizontal slowness than P 2 and S 2, respectively, indicating that the phase pairs are body waves with different ray paths. Building upon a seismic refraction profile, we construct 1D velocity models and test different thicknesses and V P / V S ratios for the subducted oceanic crust. The arrival times and relative slownesses of P 1/ P 2 and S 1/ S 2 phases indicate that they are the direct and the Moho reflected phases, respectively. Their properties are consistent with a crustal thickness of ∼6 km and a moderate V P / V S ratio (∼1.8). The S coda is more difficult to characterize but has a clear dominant frequency that likely reflects the near‐source velocity and attenuation structure. Our study indicates that waveforms from earthquakes near the interface of the subducted slab can be used to infer detailed structural information about the plate‐boundary zone at seismogenic depths.

California↗
USGS2025

Tracing mercury from land to river: Global sources, retention, and implications for sustainability

Mercury (Hg) pollution in river systems is a global sustainability challenge. Yet the transport, transformation, and retention of Hg within global rivers remain poorly quantified, particularly in regions with sparse observations such as Southeast Asia and Africa, hindering effective pollution mitigation and reinforcing geographic inequities in scientific knowledge and environmental governance. Here, we present the first global, high-resolution simulation of riverine Hg dynamics using a process-based model that traces Hg from land-based sources through river networks to the ocean. Under a realistic scenario, we estimate that ~1,900 megagrams per year (Mg/yr) of Hg enters global rivers, including 1,500 Mg/yr from human-induced sources and 400 Mg/yr from soil erosion. Nearly half of this flux (~1,000 Mg/yr) is retained in reservoirs and dams, which act as major sinks. While such retention limits downstream delivery to the oceans, it also heightens in-reservoir Hg methylation risks. By bridging the gap between Hg releases and observed riverine exports, our framework offers a scalable tool for data-limited regions, promotes data access, and supports global freshwater and pollution-management strategies.

EarthArXiv↗