skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Isotope Ratio – Discharge Relationships of Solutes Derived From Weathering Reactions
To date, the vast majority of studies seeking to link discharge to solute concentrations have been based on representations of fluid age distributions in watersheds that are time-invariant. As increasingly detailed spatial and temporal datasets become available for weathering-derived riverine solute concentrations, the capacity to link this mass flux to transient routing of reactive fluids through Critical Zone environments is vital to quantitative interpretation. Relationships between fluid age distributions and the stable isotope ratios of these geogenic solutes are even less developed, yet these signatures are vital to parsing the suite of water-rock-life interactions that create concentration-discharge relationships. Here we offer the first merging of a hydrological model featuring time-variant fluid age distributions with a geochemical model for isotopically fractionating weathering reactions. Using SiO2(aq)and the corresponding silicon isotope ratio δ 30Si as an example, we show that the stable isotope signatures of riverine solutes produced by weathering reactions reflect a component of the fluid age distribution that is unique to the corresponding solute concentrations. This distinct sensitivity is the result of a stronger link between isotope ratios and the age distribution parameters describing a given watershed. This novel modeling framework is used to provide a quantitative basis for the interpretation of SiO2(aq)and δ 30Si in six low-order streams spread across a diversity of climates, geologies, and ecosystems. To our knowledge, this is the first forward and process-based model to describe the isotopic signatures of solutes derived from weathering reactions in watersheds subject to time-varying discharge.  more » « less
Award ID(s):
2135405 2134453
PAR ID:
10482747
Author(s) / Creator(s):
;
Publisher / Repository:
American Journal of Science
Date Published:
Journal Name:
American Journal of Science
Volume:
323
ISSN:
0002-9599
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Basaltic watersheds such as those found in Iceland are thought to be important sites of CO2sequestration via silicate weathering. However, determining the magnitude of CO2uptake depends on accurately interpreting river chemistry. Here, we compile geochemical data from Iceland and use them to constrain weathering processes. Specifically, we use a newly developed inverse model to quantify solute supply from rain and hydrothermal fluids as well as allow for variable silicate end‐member compositions, solutes to be removed via secondary phase formation, and some Ca to be sourced from carbonate dissolution. While some of these processes have been considered previously, they have not been considered together allowing us to newly determine their relative contributions. We find that weathering in Iceland is incongruent in two ways. First, solute release from primary silicates is characterized by a higher proportion of Na than would be expected from bulk basalts, which may reflect preferential weathering or some contribution from rhyolites. This Na enrichment is further enhanced by preferential Mg and K uptake by secondary phases. No samples in our data set (n = 537) require carbonate dissolution even if isotopic data (δ26Mg,δ30Si,δ44Ca, and/or87Sr/86Sr) are included. While some carbonate weathering is allowable, silicate weathering likely dominates. The complexity we observe in Iceland underscores the need for inverse models to account for a wide range of processes and end‐members. Given that riverine fluxes from Iceland are more Na‐rich than expected for congruent basalt weathering, the characteristic timescale of CO2drawdown is likely affected. 
    more » « less
  2. Stream channel burial drastically alters watershed flowpaths by routing surface waters underground and increasing the potential for interactions between stream water and urban infrastructure such as storm and sanitary sewers. While numerous studies have investigated storm event solute loads from urban watersheds, the influences of stream channel burial and sewer overflows are often overlooked. This study uses grab samples and natural abundance stable isotope tracers to quantify the event dynamics of solute concentration-discharge relationships as well as cumulative loads in a buried urban stream. Our results demonstrate that different solutes, as well as different sources of the same solute (atmospheric NO3and sewer-derived NO3differentiated by the Δ17O tracer), are delivered via separate watershed flowpaths and thus have different timings within the event and contrasting relationships to flow. This inter-event variability reveals dynamics that result from temporal and spatial heterogeneity in infiltration, exfiltration, and pipe overflows. These results can help guide system-wide infrastructure maintenance as cities seek to meet challenges in sustaining and improving water quality as infrastructural systems age. 
    more » « less
  3. Abstract Blood Falls is a hypersaline, iron‐rich discharge at the terminus of the Taylor Glacier in the McMurdo Dry Valleys, Antarctica. In November 2014, brine in a conduit within the glacier was penetrated and sampled using clean‐entry techniques and a thermoelectric melting probe called the IceMole. We analyzed the englacial brine sample for filterable iron (fFe), total Fe, major cations and anions, nutrients, organic carbon, and perchlorate. In addition, aliquots were analyzed for minor and trace elements and isotopes including δD and δ18O of water, δ34S and δ18O of sulfate,234U,238U, δ11B,87Sr/86Sr, and δ81Br. These measurements were made in order to (1) determine the source and geochemical evolution of the brine and (2) compare the chemistry of the brine to that of nearby hypersaline lake waters and previous supraglacially sampled collections of Blood Falls outflow that were interpreted asend‐memberbrines. The englacial brine had higher Clconcentrations than the Blood Falls end‐member outflow; however, other constituents were similar. The isotope data indicate that the water in the brine is derived from glacier melt. The H4SiO4concentrations and U and Sr isotope suggest a high degree of chemical weathering products. The brine has a low N:P ratio of ~7.2 with most of the dissolved inorganic nitrogen in the form of NH4+. Dissolved organic carbon concentrations are similar to end‐member outflow values. Our results provide strong evidence that the original source of solutes in the brine was ancient seawater, which has been modified with the addition of chemical weathering products. 
    more » « less
  4. Abstract Rock weathering impacts atmospheric CO2levels with silicate rock dissolution removing CO2,and carbonate dissolution, pyrite oxidation, and organic rock carbon oxidation producing CO2. Glacierization impacts the hydrology and geomorphology of catchments and glacier retreat due to warming can increase runoff and initiate landscape succession. To investigate the impact of these changes on catchment scale weathering CO2balances, we report monthly samples of solute chemistry and continuous discharge records for a sequence of glacierized watersheds draining into Kachemak Bay, Alaska. We partition solute and acid sources and estimate inorganic weathering CO2balances using an inverse geochemical mixing model. Furthermore, we investigated how solutes vary with discharge conditions utilizing a concentration‐runoff framework. We develop an analogous fraction‐runoff framework which allows us to investigate changes in weathering contributions at different flows. Fraction‐runoff relationships suggest kinetic limitations on all reactions in glacierized catchments, and only silicate weathering in less glacierized catchments. Using forest cover as a proxy for landscape age and stability, multiple linear regression shows that faster reactions (pyrite oxidation) contribute less to the solute load with increasing forest cover, whereas silicate weathering (slow reaction kinetics) contributes more. Overall, in glacierized catchments, we find elevated weathering fluxes at high runoff despite significant dilution effects. This makes flux estimates that account for dilution more important in glacierized catchments. Our findings quantify how glaciers modify the inorganic weathering CO2balance of catchments through hydrologic and geomorphic forcings, and support the previous hypothesis that deglaciation will be accompanied by a shift in inorganic weathering CO2balances. 
    more » « less
  5. Abstract Lithium isotope ratios (δ7Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ7Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ7Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water‐rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope‐enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long‐term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ7Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near‐surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust. 
    more » « less