USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology Journal”

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 1,261 records · Page 70Linked to original sources

Geochemistry and hydrology of perched groundwater springs: assessing elevated uranium concentrations at Pigeon Spring relative to nearby Pigeon Mine, Arizona (USA)

The processes that affect water chemistry as the water flows from recharge areas through breccia-pipe uranium deposits in the Grand Canyon region of the southwestern United States are not well understood. Pigeon Spring had elevated uranium in 1982 (44 μg/L), compared to other perched springs (2.7–18 μg/L), prior to mining operations at the nearby Pigeon Mine. Perched groundwater springs in an area around the Pigeon Mine were sampled between 2009 and 2015 and compared with material from the Pigeon Mine to better understand the geochemistry and hydrology of the area. Two general groups of perched groundwater springs were identified from this study; one group is characterized by calcium sulfate type water, low uranium activity ratio 234 U/ 238 U (UAR) values, and a mixture of water with some component of modern water, and the other group by calcium-magnesium sulfate type water, higher UAR values, and radiocarbon ages indicating recharge on the order of several thousand years ago. Multivariate statistical principal components analysis of Pigeon Mine and spring samples indicate Cu, Pb, As, Mn, and Cd concentrations distinguished mining-related leachates from perched groundwater springs. The groundwater potentiometric surface indicates that perched groundwater at Pigeon Mine would likely flow toward the northwest away from Pigeon Spring. The geochemical analysis of the water, sediment and rock samples collected from the Snake Gulch area indicate that the elevated uranium at Pigeon Spring is likely related to a natural source of uranium upgradient from the spring and not likely related to the Pigeon Mine.

Arizona↗

Linking channel hydrology with riparian wetland accretion in tidal rivers

The hydrologic processes by which tide affects river channel and riparian morphology within the tidal freshwater zone are poorly understood, yet are fundamental to predicting the fate of coastal rivers and wetlands as sea level rises. We investigated patterns of sediment accretion in riparian wetlands along the non-tidal through oligohaline portion of two coastal plain rivers in Maryland, U.S.A., and how flow velocity, water level, and suspended sediment concentration (SSC) in the channel may have contributed to those patterns. Sediment accretion was measured over a one year period using artificial marker horizons, channel hydrology was measured over a one month period using acoustic Doppler current profilers, and SSC was predicted from acoustic backscatter. Riparian sediment accretion was lowest at the non-tidal sites (mean and standard deviation = 8 ± 8 mm yr -1 ), highest at the upstream tidal freshwater forested wetlands (TFFW) (33 ± 28 mm yr -1 ), low at the midstream TFFW (12 ± 9 mm yr -1 ), and high at the oligohaline (fresh-to-brackish) marshes (19 ± 8 mm yr -1 ). Channel maximum flood and ebb velocity was 2-fold faster at the oligohaline than tidal freshwater zone on both tidal rivers, corresponding with the differences in in-channel SSC: the oligohaline zone's SSC was more than double the tidal freshwater zone's, and was greater than historical SSC at the non-tidal gages. The tidal wave characteristics differed between rivers, leading to significantly greater in-channel SSC during floodplain inundation in the weakly convergent than the strongly convergent tidal river. Overall sediment accretion was higher in the embayed river likely due to a single storm discharge and associated sedimentation.

Maryl↗

Predictions of hydrothermal alteration within near-ridge oceanic crust from coordinated geochemical and fluid flow models

Coordinated geochemical and hydrological calculations guide our understanding of the composition, fluid flow patterns, and thermal structure of near-ridge oceanic crust. The case study presented here illustrates geochemical and thermal changes taking place as oceanic crust ages from 0.2 to 1.0 Myr. Using a finite element code, we model fluid flow and heat transport through the upper few hundred meters of an abyssal hill created at an intermediate spreading rate. We use a reaction path model with a customized database to calculate equilibrium fluid compositions and mineral assemblages of basalt and seawater at 500 bars and temperatures ranging from 150 to 400??C. In one scenario, reaction path calculations suggest that volume increases on the order of 10% may occur within portions of the basaltic basement. If this change in volume occurred, it would be sufficient to fill all primary porosity in some locations, effectively sealing off portions of the oceanic crust. Thermal profiles resulting from fluid flow simulations indicate that volume changes along this possible reaction path occur primarily within the first 0.4 Myr of crustal aging. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Detection, attribution, and sensitivity of trends toward earlier streamflow in the Sierra Nevada

Observed changes in the timing of snowmelt dominated streamflow in the western United States are often linked to anthropogenic or other external causes. We assess whether observed streamflow timing changes can be statistically attributed to external forcing, or whether they still lie within the bounds of natural (internal) variability for four large Sierra Nevada (CA) basins, at inflow points to major reservoirs. Streamflow timing is measured by "center timing" (CT), the day when half the annual flow has passed a given point. We use a physically based hydrology model driven by meteorological input from a global climate model to quantify the natural variability in CT trends. Estimated 50-year trends in CT due to natural climate variability often exceed estimated actual CT trends from 1950 to 1999. Thus, although observed trends in CT to date may be statistically significant, they cannot yet be statistically attributed to external influences on climate. We estimate that projected CT changes at the four major reservoir inflows will, with 90% confidence, exceed those from natural variability within 1-4 decades or 4-8 decades, depending on rates of future greenhouse gas emissions. To identify areas most likely to exhibit CT changes in response to rising temperatures, we calculate changes in CT under temperature increases from 1 to 5??. We find that areas with average winter temperatures between -2??C and -4??C are most likely to respond with significant CT shifts. Correspondingly, elevations from 2000 to 2800 in are most sensitive to temperature increases, with CT changes exceeding 45 days (earlier) relative to 1961-1990. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Concentrations, loads, and yields of organic carbon in streams of agricultural watersheds

Carbon is cycled to and from large reservoirs in the atmosphere, on land, and in the ocean. Movement of organic carbon from the terrestrial reservoir to the ocean plays an important role in the global cycling of carbon. The transition from natural to agricultural vegetation can change the storage and movement of organic carbon in and from a watershed. Samples were collected from 13 streams located in hydrologically and agriculturally diverse watersheds, to better understand the variability in the concentrations and loads of dissolved organic carbon (DOC) and particulate organic carbon (POC) in the streams, and the variability in watershed yields. The overall annual median concentrations of DOC and POC were 4.9 (range: 2.1–6.8) and 1.1 (range: 0.4–3.8) mg C L −1 , respectively. The mean DOC watershed yield (± SE) was 25 ± 6.8 kg C ha −1 yr −1 . The yields of DOC from these agricultural watersheds were not substantially different than the DOC yield from naturally vegetated watersheds in equivalent biomes, but were at the low end of the range for most biomes. Total organic carbon (DOC + POC) annually exported from the agricultural watersheds was found to average 0.03% of the organic carbon that is contained in the labile plant matter and top 1 m of soil in the watershed. Since the total organic carbon exported from agricultural watersheds is a relatively small portion of the sequestered carbon within the watershed, there is the great potential to store additional carbon in plants and soils of the watershed, offsetting some anthropogenic CO2 emissions.

California;Indiana;Iowa;Maryl;Mississippi;Nebraska↗

Correlated factors in amphibian decline: Exotic species and habitat change in western Washington

Amphibian declines may frequently be associated with multiple, correlated factors. In western North America, exotic species and hydrological changes are often correlated and are considered 2 of the greatest threats to freshwater systems. Bullfrog ( Rana catesbeiana ) introductions are frequently cited as a threat to lentic-breeding anurans native to western North America and are a suspected factor in the decline of red-legged frogs ( Rana aurora ) in California. Introduced fish and habitat change are cited less frequently but are equally viable hypotheses. I examined the relation among introduced species, habitat, and the distribution and abundance of red-legged frogs in western Washington. Red-legged frog occurrence in the Puget Lowlands was more closely associated with habitat structure and the presence of exotic fish than with the presence of bull-frogs. The spread of exotics is correlated with a shift toward greater permanence in wetland habitats regionally. Conservation of more ephemeral wetland habitats may have direct benefits for some native amphibians and may also reduce the threat of exotic fish and bullfrogs, both of which were associated with permanent wetlands. Research and conservation efforts for lowland anurans in the West should emphasize the complexities of multiple contributing factors to amphibian losses.

Washington↗

Assessment of basin-scale hydrologic impacts of CO2 sequestration, Illinois basin

Idealized, basin-scale sharp-interface models of CO2 injection were constructed for the Illinois basin. Porosity and permeability were decreased with depth within the Mount Simon Formation. Eau Claire confining unit porosity and permeability were kept fixed. We used 726 injection wells located near 42 power plants to deliver 80 million metric tons of CO2/year. After 100 years of continuous injection, deviatoric fluid pressures varied between 5.6 and 18 MPa across central and southern part of the Illinois basin. Maximum deviatoric pressure reached about 50% of lithostatic levels to the south. The pressure disturbance (>0.03 MPa) propagated 10-25 km away from the injection wells resulting in significant well-well pressure interference. These findings are consistent with single-phase analytical solutions of injection. The radial footprint of the CO2 plume at each well was only 0.5-2 km after 100 years of injection. Net lateral brine displacement was insignificant due to increasing radial distance from injection well and leakage across the Eau Claire confining unit. On geologic time scales CO2 would migrate northward at a rate of about 6 m/1000 years. Because of paleo-seismic events in this region (M5.5-M7.5), care should be taken to avoid high pore pressures in the southern Illinois basin. ?? 2010 Elsevier Ltd.

International Journal of Greenhouse Gas Control↗

Monitoring changes in seismic velocity related to an ongoing rapid inflation event at Okmok volcano, Alaska

Okmok is one of the most active volcanoes in the Aleutian Arc. In an effort to improve our ability to detect precursory activity leading to eruption at Okmok, we monitor a recent, and possibly ongoing, GPS-inferred rapid inflation event at the volcano using ambient noise interferometry (ANI). Applying this method, we identify changes in seismic velocity outside of Okmok’s caldera, which are related to the hydrologic cycle. Within the caldera, we observe decreases in seismic velocity that are associated with the GPS-inferred rapid inflation event. We also determine temporal changes in waveform decorrelation and show a continual increase in decorrelation rate over the time associated with the rapid inflation event. Themagnitude of relative velocity decreases and decorrelation rate increases are comparable to previous studies at Piton de la Fournaise that associate such changes with increased production of volatiles and/ormagmatic intrusion within the magma reservoir and associated opening of fractures and/or fissures. Notably, the largest decrease in relative velocity occurs along the intrastation path passing nearest to the center of the caldera. This observation, along with equal amplitude relative velocity decreases revealed via analysis of intracaldera autocorrelations, suggests that the inflation sourcemay be located approximately within the center of the caldera and represent recharge of shallow magma storage in this location. Importantly, there is a relative absence of seismicity associated with this and previous rapid inflation events at Okmok. Thus, these ANI results are the first seismic evidence of such rapid inflation at the volcano.

Alaska↗

Evaluation of potential impacts on Great Lakes water resources based on climate scenarios of two GCMs

The results of general circulation model predictions of the effects of climate change from the Canadian Centre for Climate Modeling and Analysis (model CGCM1) and the United Kingdom Meteorological Office's Hadley Centre (model HadCM2) have been used to derive potential impacts on the water resources of the Great Lakes basin. These impacts can influence the levels of the Great Lakes and the volumes of channel flow among them, thus affecting their value for interests such as riparians, shippers, recreational boaters, and natural ecosystems. On one hand, a hydrological modeling suite using input data from the CGCM1 predicts large drops in lake levels, up to a maximum of 1.38 m on Lakes Michigan and Huron by 2090. This is due to a combination of a decrease in precipitation and an increase in air temperature that leads to an increase in evaporation. On the other hand, using input from HadCM2, rises in lake levels are predicted, up to a maximum of 0.35 m on Lakes Michigan and Huron by 2090, due to increased precipitation and a reduced increase in air temperature. An interest satisfaction model shows sharp decreases in the satisfaction of the interests of commercial navigation, recreational boating, riparians, and hydropower due to lake level decreases. Most interest satisfaction scores are also reduced by lake level increases. Drastic reductions in ice cover also result from the temperature increases such that under the CGCM1 predictions, most of Lake Erie has 96% of its winters ice-free by 2090. Assessment is also made of impacts on the groundwater-dependent region of Lansing, Michigan.

Great Lakes basin↗

A database of natural monthly streamflow estimates from 1950 to 2015 for the conterminous United States

Quantifying and understanding the natural streamflow regime, defined as expected streamflow that would occur in the absence of anthropogenic modification to the hydrologic system, is critically important for the development of management strategies aimed at protecting aquatic ecosystems. Water balance models have been applied frequently to estimate natural flows, but are limited in the number of predictor variables that can be included. Here, a statistical machine learning technique — random forest modeling — was applied to estimate natural flows at a monthly time‐step from 1950 to 2015 for >2.5 million stream reaches in the conterminous United States (U.S.) using 200 potential predictor variables. We describe the development and documentation of this dataset and assess model performance. Model fit statistics (mean Nash–Sutcliffe efficiency = 0.85; observed/expected ratio = 0.94) indicate good correspondence between predicted and observed flows at nearly 2,000 streamgages. As an example application of the dataset, the observed streamflow record at a site prior to and after the construction of an upstream reservoir was compared with estimated natural flows to demonstrate the magnitude of seasonal depletions in streamflow due to the reservoir. This dataset can be applied to quantify natural and anthropogenic processes contributing to streamflow depletion or augmentation, and assess associated ecological effects.

Journal of the American Water Resources Associatio↗

Late Quaternary paleoenvironments of an ephemeral wetland in North Dakota, USA: Relative interactions of ground-water hydrology and climate change

This study of fossils (pollen, plant macrofossils, stomata and fish) and sediments (lithostratigraphy and geochemistry) from the Wendel site in North Dakota, USA, emphasizes the importance of considering ground-water hydrology when deciphering paleoclimate signals from lakes in postglacial landscapes. The Wendel site was a paleolake from about 11,500 14C yr BP to 11,100 14C yr BP. Afterwards, the lake-level lowered until it became a prairie marsh by 9,300 14C yr BP and finally, at 8,500 14C yr BP, an ephemeral wetland as it is today. Meanwhile, the vegetation changed from a white spruce parkland (11,500 to 10,500 14C yr BP) to deciduous parkland, followed by grassland at 9,300 14C yr BP. The pattern and timing of these aquatic and terrestrial changes are similar to coeval kettle lake records from adjacent uplands, providing a regional aridity signal. However, two local sources of ground water were identified from the fossil and geochemical data, which mediated atmospheric inputs to the Wendel basin. First, the paleolake received water from the melting of stagnant ice buried under local till for about 900 years after glacier recession. Later, Holocene droughts probably caused the lower-elevation Wendel site to capture the ground water of up-gradient lakes. ?? 2007 Springer Science+Business Media, Inc.

Journal of Paleolimnology↗

Responses of biological soil crusts to rehabilitation strategies

Biological soil crusts (biocrusts) are common to dryland ecosystems and can influence a broad suite of soil ecological functions including stability and surface hydrology. Due to long recovery times following disturbance, there is a clear need for rehabilitation strategies to enhance the recovery of biocrust communities. Essential to biocrust recovery are exopolysaccharides (EPS): secretions comprised mainly of high molecular weight polymers that protect cyanobacteria from harsh environmental conditions. We examined whether biocrust rehabilitation strategies (inoculation in combination with surface shading and artificial soil stabilization) promote EPS production. To test if responses varied by soil texture, we measured biocrust recovery on two fine-textured soil types (clay and sandy clay loam) in a cool desert ecosystem. Shade coupled with inoculum addition resulted in the highest biocrust recovery, especially on clay soils. Independent of rehabilitation strategies, natural recovery of biocrusts occurred more rapidly on clay soils, reflected by greater increases in chlorophyll a (chl a). Chl a, a proxy for cyanobacterial biomass, was correlated to EPS amounts, suggesting that cyanobacteria are significant contributors to EPS production in biocrust development. Despite the role of EPS in biocrust establishment, EPS amounts had negligible effects on soil stability due inherent properties of fine soil texture.

Utah↗

Spatiotemporal variability of modeled watershed scale surface-depression storage and runoff for the conterminous United States

This study uses the explores the viability of a proxy model calibration strategy through assessment of the spatiotemporal variability of surface-depression storage and runoff generated with the U.S. Geological Survey’s National Hydrologic Model (NHM) infrastructure for hydrologic response units (HRUs; n=109,951) across the conterminous United States (CONUS). Simulated values for each HRU of daily surface-depression storage (treated as a decimal fraction of total possible volume) and monthly normalized runoff (0 to 1) values were calculated using Spearman’s rho at monthly and annual aggregations. Locations where values are correlated show where previously-developed proxy calibration strategies are likely to be effective. In addition, differences in the correlation for monthly and annual time scale aggregations show which time scale drives surface-depression storage processes in the NHM. Results show overall long-term (annual) correlation is more common than short-term (monthly) correlation over the CONUS; however, summary statistics for eighty-six ecoregions show five with higher ranges of monthly relative to annual Spearman’s rank coefficient values. This landscape-scale analysis shows simulations aggregated to an annual time scale are generally more dominant for the CONUS; however, simulations aggregated to monthly, short-term time scales are more dominant in focused areas where surface-depression storage processes are investigated.

Journal of the American Water Resources Associatio↗

Geochemistry and Soils of the Big Smoky Valley Fens, Nevada

Fens are groundwater-fed wetlands that can provide habitat for plants and animals. Due to anthropogenic activities and climate change, many fens around the world are at risk. This paper presents the results of a study of the hydrology and geochemistry of fens in Big Smoky Valley, central Nevada to support the Bureau of Land Management’s activities in the area. A water sample from the largest fen in the study area was analyzed for its water chemistry and compared to a nearby alluvial aquifer and hot spring. The high SiO 2 concentration of the fen sample implies that the fen water may originate from geothermal water. A soil core was taken to analyze radiocarbon age and soil type. A majority of the core was composed of silt and clay interlayered with water-filled voids. Changes in the character of the clay with depth suggest that there may have been changes in the depositional environment over time. Radiocarbon dating of Ruppia seeds showed longevity of the fen, with the minimum 14 C age of the core as 4,375±40 years. This paper provides reconnaissance-level information on the Big Smoky Valley fens, but further information would be needed to better understand the source of water to the fens or how the fen environment has changed over time with climate.

Nevada↗

Adapting standardized trout monitoring to a changing climate for the upper Yellowstone River, Montana, USA

Objective Long‐term standardized monitoring programs are fundamental to assessing how fish populations respond to anthropogenic stressors. Standardized monitoring programs may need to adopt new methods to adapt to rapid environmental changes that are associated with a changing climate. In the upper Yellowstone River, Montana, biologists have used a standardized, mark–recapture monitoring protocol to annually estimate the abundance of trout since 1978 to assess population status and trends. However, within the past two decades, climate change has caused changes in discharge timing that have prevented standardized monitoring from occurring annually. Methods We investigated the feasibility of using two analytical methods, N‐mixture models and mean capture probability, for estimating the abundance of three trout species in the upper Yellowstone River using the historical long‐term data set; these methods allow abundance to be estimated when a mark–recapture estimate cannot be obtained due to hydrologic conditions. Result When compared with abundance estimates from mark–recapture methods, N‐mixture models most often resulted in negatively biased abundance estimates, whereas mean capture probability analyses resulted in positively biased abundance estimates. Additionally, N‐mixture models produced negatively biased estimates when tested against true abundance values from simulated data sets. The bias in the N‐mixture model estimates was caused by poor model fit and variation in capture probability. The bias in the mean capture probability estimates was caused by heterogeneity in capture probability, likely caused by variable environmental conditions, which were not accounted for in the models. Conclusion N‐mixture models and mean capture probability are not viable alternatives for estimating abundance in the upper Yellowstone River. Thus, exploring additional adaptations to sampling methodologies and analytical approaches, including models that require individually marked fish, will be valuable for this system. Climate change will undoubtedly necessitate changes to standardized sampling methods throughout the world; thus, developing alternative sampling and analytical methods will be important for maintaining the utility of long‐term data sets.

Montana↗

An empirical assessment of which inland floods can be managed

Riverine flooding is a significant global issue. Although it is well documented that the influence of landscape structure on floods decreases as flood size increases, studies that define a threshold flood-return period, above which landscape features such as topography, land cover and impoundments can curtail floods, are lacking. Further, the relative influences of natural versus built features on floods is poorly understood. Assumptions about the types of floods that can be managed have considerable implications for the cost-effectiveness of decisions to invest in transforming land cover (e.g., reforestation) and in constructing structures (e.g., storm-water ponds) to control floods. This study defines parameters of floods for which changes in landscape structure can have an impact. We compare nine flood-return periods across 31 watersheds with widely varying topography and land cover in the southeastern United States, using long-term hydrologic records (≥20 years). We also assess the effects of built flow-regulating features (best management practices and artificial water bodies) on selected flood metrics across urban watersheds. We show that landscape features affect magnitude and duration of only those floods with return periods ≤10 years, which suggests that larger floods cannot be managed effectively by manipulating landscape structure. Overall, urban watersheds exhibited larger (270 m 3 /s) but quicker (0.41 days) floods than non-urban watersheds (50 m 3 /s and 1.5 days). However, urban watersheds with more flow-regulating features had lower flood magnitudes (154 m 3 /s), but similar flood durations (0.55 days), compared to urban watersheds with fewer flow-regulating features (360 m 3 /s and 0.23 days). Our analysis provides insight into the magnitude, duration and count of floods that can be curtailed by landscape structure and its management. Our findings are relevant to other areas with similar climate, topography, and land use, and can help ensure that investments in flood management are made wisely after considering the limitations of landscape features to regulate floods.

Journal of Environmental Management↗

The Metzger marsh restoration: A vegetation-centric look after 27 years

We investigated wetland vegetation before, during, and after dike construction at the Metzger Marsh project in western Lake Erie, which was designed to restore a 300-ha wetland that had been degraded following the loss of a protective barrier beach. A dike was constructed in 1995 to replace the function of the eroded barrier beach, but it contained a water-control structure to allow managed hydrologic connection to the lake. The control structure contained a fish passageway to allow movement of fish across the dike, while restricting entry of large common carp. Color-infrared aerial photos from project start in 1994 through 2010 (and 2022) were analyzed to track vegetation changes, and major vegetation types were sampled quantitatively. Drawdown of water levels in 1996 after dike construction elicited a response of mudflat species from the seed bank, as well as tree seedlings. Over half of the marsh was vegetated then and in subsequent years. The water-control structure was opened in 1998, and by 2000, invasive Phragmites australis had gained dominance. Most trees were eventually eliminated by herbicide treatment and flooding, and extent of Phragmites was reduced by management actions. Typha spp. and emergents Sagittaria latifolia and Schoenoplectus tabernaemontani became dominant by 2022. This restoration project increased habitat values for fish and wildlife; it also provided lessons for future projects on lands managed by multiple agencies with differing missions. More importantly, it showed that long-term monitoring data are critical for assessing wetland restoration projects and guiding management decisions.

Ohio↗

Nutrient dynamics in partially drained arctic thaw lakes

Thaw lakes are ubiquitous on arctic coastal plains (ACPs). While many thaw lakes have steep banks, stable water levels, and static surface areas, others only partially fill their basins and vary in area over the summer. These partially drained lakes (PDLs) are hydrologically connected to the wetlands immediately surrounding them. Heat and nutrient availability limit aquatic productivity on ACPs, and we hypothesized that shallow shorelines and greater hydrologic connectivity with the landscape should result in greater nutrient concentrations and biogeochemical cycling in PDLs. We tested this by monitoring water chemistry in lakes with varying levels of seasonal drainage in sandy and silty peaty lowland sites on the ACP of Alaska. One highly drained lake (N1) was significantly warmer than minimally drained lakes (minDLs) related to earlier ice off, reaching temperatures as high as 16 °C in June when minDLs still contained ice. Ammonia, total dissolved phosphorus, and dissolved organic carbon and nitrogen concentrations were higher in lakes with greater drainage, and concentrations in N1 rivaled those in the small, biologically productive ponds. Many PDLs displayed a midsummer decrease in nutrients consistent with assimilation by the aquatic ecosystem, and a late‐summer increase most likely related to runoff from drained lake margins following precipitation. N1 exported kilograms of ammonium and total dissolved phosphorus to the stream network over the summer. Given increased warming and drying in the arctic, the proportion of PDLs may be changing, which in turn may affect nutrient and organic matter availability in arctic lakes and export to downstream environments.

Journal of Geophysical Research: Biogeosciences↗