Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hillslope”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Vegetation heterogeneity reflects soil thermal state and surface soil displacement in a thawing permafrost landscape

Thawing permafrost has the potential to dramatically alter the physical and ecological structure of northern landscapes. Warming of the Arctic and subsequent degradation of permafrost have created a need to assess the stability and movement of soils on hillslopes and the potential impacts on ecosystem structure. In this work, we explore the relationships among vegetation heterogeneity, soil temperature, and soil surface displacements observed from 2019 to 2022 in a watershed in the discontinuous permafrost region on the Seward Peninsula of Alaska. Vegetation heterogeneity was measured as the standard deviation (SD) of the normalized difference vegetation index (NDVI) from 3 m PlanetScope satellite imagery around each soil temperature and active layer thickness observation. Locations of observations were clustered into three soil thermal groups, warm, intermediate, and cold, based on soil temperature and active layer thickness. Average annual horizontal surface displacements were significantly lower for soils within the warm thermal group (median = 0.033 m yr −1 ) compared to soils within the cold thermal group (median = 0.090 m yr −1 ; p < 0.001). Conversely, vegetation heterogeneity was significantly higher in the warm (median = 0.014 SD NDVI; p = 0.002) and intermediate (median = 0.015 SD NDVI; p = 0.002) groups compared with the cold thermal group (median = 0.012 SD NDVI), suggesting a warming-induced shift in vegetation community complexity. Because of the observed associations of ground surface displacement rates and vegetation heterogeneity with soil thermal state, we hypothesize that warming soil conditions induce changes in the rates and patterns of hillslope erosion due to an increase in surface movement as near-surface permafrost thaws, followed by a decrease as the permafrost table deepens and excess ice content diminishes. The transition to warm soils promotes surface ecosystem transformation, shifting the dominant vegetation at the site, given the warming climatic conditions of the region. We integrated our observations of soil temperature, vegetation heterogeneity, and soil surface displacements into a conceptual model that describes the co-evolution of hillslopes and vegetation in warming permafrost environments, which is currently unrepresented in earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Model and remote-sensing-guided experimental design and hypothesis generation for monitoring snow-soil–plant interactions

In this study, we develop a machine-learning (ML)-enabled strategy for selecting hillslope-scale ecohydrological monitoring sites within snow-dominated mountainous watersheds, with a particular focus on snow-soil–plant interactions. Data layers rely on spatial data layers from both remote sensing and hydrological model simulations. Specifically, a Landsat-based foresummer drought sensitivity index is used to define the dependency of the annual peak plant productivity on the Palmer drought severity index in the early growing season. Hydrological simulations provide the spatiotemporal dynamics of near-surface soil moisture and snow depth. In this framework, a regression analysis identifies the key hydrological variables relevant to the spatial heterogeneity of drought sensitivity. We then apply unsupervised clustering to these key variables, using the Gaussian mixture model, to group hillslopes into several zones that have divergent relationships regarding soil moisture, snow dynamics, and drought sensitivity. Using the datasets collected in the East River Watershed (Crested Butte, Colorado, United States), results show that drought sensitivity is significantly correlated with model-derived soil moisture and snow-free timing over space and time. The relationship is, however, non-linear, such that the correlation decreases above a threshold elevation and in a heavy snow year due to large snowpacks, lateral flow, and soil storage limitations. Clustering is then able to define the zones that have high or low sensitivity to drought, as well as the mid-elevation regions where sensitivity is associated with the topographic aspect and net potential radiation. In addition, the algorithm identifies the most representative hillslopes with road/trail access within each zone for installing monitoring sites. Our method also aims to significantly increase the use of ML and model-simulation results to guide critical zone and watershed monitoring activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Importance of Accounting for Landscape Position When Investigating Grasslands: A Multidisciplinary Characterisation of a California Coastal Grassland

Abstract Grasslands are one of the most common land‐cover types, providing important ecosystem services globally, yet few studies have examined grassland critical‐zone functioning throughout hillslopes. This study characterised a coastal grassland over a small hillslope at Point Reyes National Seashore, California, using multidisciplinary techniques, combining remotely‐sensed, geophysical, plant, and soil measurements. Clustering techniques delineated the study area into four landscape zones, up‐, mid‐, and down‐slope, and a bordering riparian ecotone, which had distinct environmental properties that varied spatially across the site, with depth, and time. Soil moisture increased with depth and down slope towards a bordering riparian zone, and co‐varied with soil CO 2 flux rates both spatially and temporally. This highlighted three distinct controls of soil moisture on soil respiration: CO 2 fluxes were inhibited by high moisture content in the down‐slope during the wet winter months, and converged across landscape positions in the dry summer months, while also displaying post‐rain pulses. The normalised difference vegetation index (NDVI) ranged from 0.32 (September)–0.80 (April) and correlated positively with soil moisture and aboveground biomass, moving down slope. Yet, NDVI, aboveground biomass, and soil moisture were not correlated to soil organic carbon (SOC) content (0.4%–4.5%), which was highest in the mid‐slope. The SOC content may instead be linked to shifts in dominant grassland species and their rhizosphere properties with landscape position. This multidisciplinary characterisation highlighted significant heterogeneity in grassland properties with landscape position, and demonstrated an approach that could be used to characterise other critical‐zone environments on hillslopes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

Abstract Lateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Water Level Data from Wells PLM1 and PLM6 for the East River Watershed, Colorado

This dataset (Williams et al., 2020) contains the original un-QA/QC-ed water level data for PLM1 and PLM6 and has been obsoleted. The data contained within this dataset is not to be used. Refer to Faybishenko et al., 2022 (DOI: 10.15485/1866836) for the latest QA/QC-ed data available via ESS-DIVE.This data set contains water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location. These wells used to monitor subsurface water and carbon inventories and fluxes at the East River Watershed, Colorado, USA. Complete metadata information on the PLM1 and PLM6 wells are available in the related data package reference Varadharajan C, et al (2020). https://doi.org/10.15485/1660962.Data are reported in .csv files per well. The latitude and longitude of each location are given in a file called locations.csv. These data are used for determining the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry can be used to estimate rates of solute export from the hillslope to the floodplain and river.These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

QA/QC-ed Groundwater Level Time Series in PLM-1 and PLM-6 Monitoring Wells, East River, Colorado (2016-2022)

This data set contains QA/QC-ed (Quality Assurance and Quality Control) water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location at the East River Watershed, Colorado, USA. These wells are used to monitor subsurface water and carbon inventories and fluxes, and to determine the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry (see related references) can be used to estimate rates of solute export from the hillslope to the floodplain and river. QA/QC analysis of measured groundwater levels in monitoring wells PLM-1 and PLM-6 included identification and flagging of duplicated values of timestamps, gap filling of missing timestamps and water levels, removal of abnormal/bad and outliers of measured water levels. The QA/QC analysis also tested the application of different QA/QC methods and the development of regular (5-minute, 1-hour, and 1-day) time series datasets, which can serve as a benchmark for testing other QA/QC techniques, and will be applicable for ecohydrological modeling. The package includes a Readme file, one R code file used to perform QA/QC, a series of 8 data csv files (six QA/QC-ed regular time series datasets of varying intervals (5-min, 1-hr, 1-day) and two files with QA/QC flagging of original data), and three files for the reporting format adoption of this dataset (InstallationMethods, file level metadata (flmd), and data dictionary (dd) files).QA/QC-ed data herein were derived from the original/raw data publication available at Williams et al., 2020 (DOI: 10.15485/1818367). For more information about running R code file (10.15485_1866836_QAQC_PLM1_PLM6.R) to reproduce QA/QC output files, see README (QAQC_PLM_readme.docx). This dataset replaces the previously published raw data time series, and is the final groundwater data product for the PLM wells in the East River. Complete metadata information on the PLM1 and PLM6 wells are available in a related dataset on ESS-DIVE: Varadharajan C, et al (2022). https://doi.org/10.15485/1660962. These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales. 2022/09/09 Update: Converted data files using ESS-DIVE’s Hydrological Monitoring Reporting Format. With the adoption of this reporting format, the addition of three new files (v1_20220909_flmd.csv, V1_20220909_dd.csv, and InstallationMethods.csv) were added. The file-level metadata file (v1_20220909_flmd.csv) contains information specific to the files contained within the dataset. The data dictionary file (v1_20220909_dd.csv) contains definitions of column headers and other terms across the dataset. The installation methods file (InstallationMethods.csv) contains a description of methods associated with installation and deployment at PLM1 and PLM6 wells. Additionally, eight data files were re-formatted to follow the reporting format guidance (er_plm1_waterlevel_2016-2020.csv, er_plm1_waterlevel_1-hour_2016-2020.csv, er_plm1_waterlevel_daily_2016-2020.csv, QA_PLM1_Flagging.csv, er_plm6_waterlevel_2016-2020.csv, er_plm6_waterlevel_1-hour_2016-2020.csv, er_plm6_waterlevel_daily_2016-2020.csv, QA_PLM6_Flagging.csv). The major changes to the data files include the addition of header_rows above the data containing metadata about the particular well, units, and sensor description. 2023/01/18 Update: Dataset updated to include additional QA/QC-ed water level data up until 2022-10-12 for ER-PLM1 and 2022-10-13 for ER-PLM6. Reporting format specific files (v2_20230118_flmd.csv, v2_20230118_dd.csv, v2_20230118_InstallationMethods.csv) were updated to reflect the additional data. R code file (QAQC_PLM1_PLM6.R) was added to replace the previously uploaded HTML files to enable execution of the associated code. R code file (QAQC_PLM1_PLM6.R) and ReadMe file (QAQC_PLM_readme.docx) were revised to clarify where original data was retrieved from and to remove local file paths.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Runoff evaluation in an Earth System Land Model for permafrost regions in Alaska

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce. In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-based integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM's parameterized representation of total runoff. Results show that (1) ELM's total runoff was the same order of magnitude as the ATS simulations, and both models were similarly variable over time; (2) minor adjustments to coefficients in ELM's runoff parameterization improved the match between the ATS simulation and ELM's parameterized representation of annual and seasonal total runoff; (3) overall, runoff responses in ATS and ELM are more similar in flat hillslope environments compared to steep hillslopes; and (4) shallower active layer thicknesses and higher precipitation simulations resulted in lower correlations between the two models due to greater total runoff. By incorporating the optimized runoff coefficients from the Sagavanirktok River basin into ELM, the simulated total runoff better matched the streamflow observations at a small watershed located on the Seward Peninsula of Alaska. Our findings revealed important insights into the effectiveness of runoff parameterizations in land surface models and pathways for improving runoff coefficients in typical Arctic regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluating simplifications of subsurface process representations for field-scale permafrost hydrology models

Abstract. Permafrost degradation within a warming climate poses a significant environmental threat through both the permafrost carbon feedback and damage to human communities and infrastructure. Understanding this threat relies on better understanding and numerical representation of thermo-hydrological permafrost processes and the subsequent accurate prediction of permafrost dynamics. All models include simplified assumptions, implying a tradeoff between model complexity and prediction accuracy. The main purpose of this work is to investigate this tradeoff when applying the following commonly made assumptions: (1) assuming equal density of ice and liquid water in frozen soil, (2) neglecting the effect of cryosuction in unsaturated freezing soil, and (3) neglecting advective heat transport during soil freezing and thaw. This study designed a set of 62 numerical experiments using the Advanced Terrestrial Simulator (ATS v1.2) to evaluate the effects of these choices on permafrost hydrological outputs, including both integrated and pointwise quantities. Simulations were conducted under different climate conditions and soil properties from three different sites in both column- and hillslope-scale configurations. Results showed that amongst the three physical assumptions, soil cryosuction is the most crucial yet commonly ignored process. Neglecting cryosuction, on average, can cause 10 %–20 % error in predicting evaporation, 50 %–60 % error in discharge, 10 %–30 % error in thaw depth, and 10 %–30 % error in soil temperature at 1 m beneath the surface. The prediction error for subsurface temperature and water saturation is more obvious at hillslope scales due to the presence of lateral flux. By comparison, using equal ice–liquid density has a minor impact on most hydrological metrics of interest but significantly affects soil water saturation with an averaged 5 %–15 % error. Neglecting advective heat transport presents the least error, 5 % or even much lower, in most metrics of interest for a large-scale Arctic tundra system without apparent influence caused by localized groundwater flow, and it can decrease the simulation time at hillslope scales by 40 %–80 %. By challenging these commonly made assumptions, this work provides permafrost hydrology scientists an important context for understanding the underlying physical processes, including allowing modelers to better choose the appropriate process representation for a given modeling experiment.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Correspondence of surface temperatures and terrain variables over a tallgrass prairie

The time-dependent correspondence between maps of surface brightness temperature derived from Landsat Thematic Mapper data and mapped terrain variables over a tallgrass prairie in northeastern Kansas is examined. Individual terrain variables including burning treatment, vegetation cover type (agriculture, prairie, woody vegetation), hillslope position, and greenness exhibit varying degrees of association with surface temperature. Burning treatment is most strongly associated with mid-morning surface temperature. Examination of terrain strata based on combinations of terrain variables, notably burning treatment and hillslope position, suggest that terrain variables interact in affecting surface temperature. Interaction between hillslope position, burning treatment, and surface temperature is more important in August than in May.

Friedl, Mark A.↗
DOE OSTI2021

Amount and reactivity of dissolved organic matter export are affected by land cover change from old‐growth to second‐growth forests in headwater ecosystems

Abstract Headwater forest ecosystems of the western USA generate a large portion of the dissolved organic matter (DOM) transported in streams across North America. Land cover changes that alter forest structure and species composition affect the quantity and composition of DOM transferred to aquatic ecosystems. Clear‐cut harvesting affects ~1% of the forest area of North America annually, leaving most western forests in varying stages of regrowth and the total area of old‐growth forest is decreasing. The consequences of this widespread management practice on watershed carbon cycling remain unknown. We investigated the role of land cover change, because of clear‐cut harvesting, from mixed‐species old‐growth to lodgepole pine‐dominated second‐growth forest on the character and reactivity of hillslope DOM exports. We evaluated inputs of DOM from litter leachates and export of DOM collected at the base of trenched hillslopes during a 3‐year period (2016–2018) at the Fraser Experimental Forest in north‐central Colorado, USA. Dissolved organic carbon and total dissolved nitrogen were higher in lateral subsurface flow draining old‐ versus second‐growth forest. Fluorescence spectroscopy showed that the DOM exported from the old‐growth forest was more heterogeneous and aromatic and that proteinaceous, microbially processed DOM components were more prevalent in the second‐growth forest. Biological oxygen demand assays revealed much lower microbial metabolism of DOM in litter leachate and subsurface exports from the old‐growth forest relative to second growth. Old‐growth and second‐growth forests are co‐mingled in managed ecosystems, and our findings demonstrate that land cover change from a mixture of conifer species to lodgepole pine dominance influences DOM inputs that can increase the reactivity of DOM transferred from terrestrial to aquatic ecosystems.

Fegel, Timothy S.↗
DOE OSTI2023

Deep root activity overprints weathering of petrogenic organic carbon in shale

The oxidation of organic carbon in sedimentary bedrock (petrogenic OC, OC petro ) is increasingly recognized as a potential source of CO 2 to the atmosphere. Recent studies provide evidence for the mobilization and oxidation of OC petro in sedimentary bedrock during rock weathering. However, the mechanisms and rates remain uncertain, particularly where overlying soils and vegetation drive contemporaneous oxidation of recently fixed organic carbon. Here, in this study, we quantify OC petro weathering across a 16 m shale depth profile in a steep, rapidly eroding forested hillslope in the Northern California Coast Ranges. We report solid and gas phase radiocarbon and stable isotope analyses of samples extracted from specialized in-situ samplers, and a supporting laboratory incubation experiment of the shale regolith. OC petro is removed from the weathered bedrock at a rate of approximately 0.12 gC/m 3 yr, which is orders of magnitude lower than the rate of OC petro oxidation we achieved in the laboratory with crushed samples (557.1 gC/m 3 yr). This disparity occurs despite high O 2(g) content across the depth profile, indicating that physical accessibility of OC petro can regulate oxidative weathering. There is no direct radiocarbon evidence of OC petro oxidation in CO 2(g) across the upper 13 m of the weathering profile during both wet and dry seasons. Instead, vadose zone CO 2(g) production at the site is dominated by respiration of recently fixed carbon associated with deep rooting. OC petro is clearly mobilized across the vadose zone during weathering in this rapidly eroding, oxygen-rich, biologically dynamic hillslope, but at rates far below what can be measured given the contribution of root-derived CO 2(g) .

58 GEOSCIENCES↗
DOE OSTI2022

Dynamic coevolution of baseflow and multiscale groundwater flow system during prolonged droughts

Field and numerical studies suggest that baseflow is composed of waters from a spectrum of groundwater flow paths termed the Groundwater Flow System (GWFS) – from shallow hillslope contributions to watershed-scale deep circulation originating in headwaters and discharging into lowland rivers. Here, we explore the evolution of the GWFS under prolonged droughts to understand its dynamics and multiscale nature, and to elucidate its role in baseflow generation and recession at the watershed scale. In this work, we consider three drought scenarios of varying severity and simulate groundwater flow in a 2-D cross-section of an idealized watershed with deep permeable bedrock, tracking the evolution of flow paths, baseflow, and residence times during the recession process. We find that baseflow generation at different drainage stages, and within different subwatersheds, is influenced distinctly by flow paths of different scales, depending on the relative strength of the flow paths and the position of the subwatersheds relative to the recharge/discharge zones of the deeper watershed-scale groundwater circulation. Despite having the same local relief, geology, and climate, baseflow from each subwatershed has a distinct recession behavior and time-dependent residence time distribution. Also, the hydraulic and transport characteristics of baseflow generation co-evolve and are strongly affected by the connection state of the water table to subwatersheds. These findings suggest that asynchrony and dissimilarity of baseflow generation from hillslopes under the impact of the watershed-scale groundwater flow, and interactions with local-scale and intermediate-scale groundwater flow, must be taken into account when interpreting baseflow recession data and building conceptual baseflow models at the watershed scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dynamics of streamflow permanence in a headwater network: Insights from catchment-scale model simulations

The hillslope and channel dynamics that govern streamflow permanence in headwater systems have important implications for ecosystem functioning and downstream water quality. Recent advancements in process-based, semi-distributed hydrologic models that build upon empirical studies of streamflow permanence in well-monitored headwater catchments show promise for characterizing the dynamics of streamflow permanence in headwater systems. However, few process-based models consider the continuum of hillslope-stream network connectivity as a control on streamflow permanence in headwater systems. The objective of this study was to expand a process-based, catchment-scale hydrologic model to better understand the spatiotemporal dynamics of headwater streamflow permanence and to identify controls of streamflow expansion and contraction in a headwater network. Further, we aimed to develop an approach that enhanced the fidelity of model simulations, yet required little additional data, with the intent that the model might be later transferred to catchments with limited long-term and spatially explicit measurements. This approach facilitated network-scale estimates of the controls of streamflow expansion and contraction, albeit with higher degrees of uncertainty in individual reaches due to data constraints. Our model simulated that streamflow permanence was highly dynamic in first-order reaches with steep slopes and variable contributing areas. The simulated stream network length ranged from nearly 98±2% of the geomorphic channel extent during wet periods to nearly 50±10% during dry periods. The model identified a discharge threshold of approximately 1 mm d –1 , above which the rate of streamflow expansion decreases by nearly an order of magnitude, indicating a lack of sensitivity of streamflow expansion to hydrologic forcing during high-flow periods. Overall, we demonstrate that process-based, catchment-scale models offer important insights on the controls of streamflow permanence, despite uncertainties and limitations of the model. We encourage researchers to increase data collection efforts and develop benchmarks to better evaluate such models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sulfur Biogeochemical Cycling and Redox Dynamics in a Shale-Dominated Mountainous Watershed

Sulfur (S) is an essential macronutrient and important component of the earth’s crust, and its cycling has critical impacts on trace metal mobility, water quality, and human health. Pyrite weathering is the primary pathway by which sulfur enters surface waters. However, biogeochemical cycling of sulfur in soils and the river corridor mediates sulfate exports. In this study, we identified the major forms of sulfur across multiple compartments and scales in a pristine mountainous watershed, including shale bedrock weathering profiles, hillslope soils, and alluvial floodplain sediments, in order to provide insight into biogeochemical sulfur cycling in a hydrologically variable alpine system. X-ray absorption near-edge spectroscopy (XANES) analysis of shale weathering profiles showed clear evidence of pyrite oxidation to sulfate, with large accumulations of intermediate S(0) (20%–53%). Micro-scale XANES showed evidence of reprecipitation of pyrite at fracture surfaces within the permanently saturated zone. Organic sulfur dominated S speciation in shallow hillslope soil and floodplain sediment, with little evidence of reduced inorganic S. However, mackinawite formation, representing active sulfate reduction, was observed in saturated oxbow sediments and saturated weathered shale underlying floodplain sediments. Further evidence of sulfate reduction from aqueous sulfur isotopic analysis was observed in shallow groundwater transects across an Fe-reducing meander, whereas increases in pore water sulfate concentrations implied sulfur oxidation at other locations. Furthermore, the data present an integrated picture of sulfur cycling in a shale-dominated watershed, where riverine sulfate exports are mediated by biological cycling, particularly in redox-stratified and temporally dynamic hyporheic zone sediments.

58 GEOSCIENCES↗
DOE OSTI2023

Fracture Intensity Impacts on Reaction Front Propagation and Mineral Weathering in Three‐Dimensional Fractured Media

Abstract Studying reaction front propagation in heterogeneous natural settings is challenging, but numerical simulations can provide insight into the varying spatial and temporal scales of reaction front propagation. Here, the impact of increasing fracture intensity on mineral dissolution rates, and the extent of reaction front propagation is investigated using reactive transport simulations in upscaled discrete fracture network domains with varied fracture intensity. Domain‐averaged dissolution rates vary less than 0.5 log units regardless of the fracture intensity, but the spatial distribution of reactions is controlled by the location and number of dead‐end fractures and the number of connected flowpaths through the domain. Higher fracture intensities lead to more weathering in the domain because of more available mineral for water‐rock interactions. We find that reaction fronts propagate through the primary flowpaths in the first 10,000 years of the simulation for a 10‐m length domain, then propagate into secondary flowpaths and dead‐end fractures between 10,000 and 100,000 years, and finally into the matrix over timescales of hundreds of thousands of years. The domain‐averaged reaction rates decrease through time corresponding to a transition from dissolution in advection‐dominated, fast‐flowing pathways, to dissolution in transport‐limited zones of disconnected fractures and matrix. Matrix dissolution, or dissolution under transport‐limited conditions, is the dominant process at late times in these simulations. The results of these simulations recreate the observed paradox found in nature where highly fractured hillslopes tend to be more weathered but have slower weathering rates, while hillslopes with fewer fractures, are less weathered but have higher dissolution rates.

58 GEOSCIENCES↗
DOE OSTI2023

Constraining Bedrock Groundwater Residence Times in a Mountain System With Environmental Tracer Observations and Bayesian Uncertainty Quantification

Groundwater residence time distributions provide fundamental insights on the hydrological processes within watersheds. Yet, observations that can constrain groundwater residence times over broad timescales remain scarce in mountain catchment studies. We use environmental tracers (CFC-12, SF 6 , 3 H, and 4 He) to investigate groundwater residence times along a hillslope in the East River Watershed, Colorado, USA. We develop a Bayesian inference framework that applies a Markov-chain Monte Carlo (MCMC) approach to estimate noble gas recharge temperature, elevation, and excess-air parameters and the resulting environmental tracer concentrations. MCMC is then used to propagate the environmental tracer uncertainties to estimates of groundwater mean residence times inferred with lumped parameter models. All samples contain 3 H, CFC-12, and SF 6 in addition to terrigenic 4 He, suggesting a mixture of water characterized by modern and premodern residence times. 4He exponential mean residence times range from hundreds of years at the upslope well to thousands of years at the toe-slope well assuming average crustal production rates. We find that binary mixing residence time distributions with separate young and old mixing fractions are needed to predict the 4 He, CFC-12, SF 6 , and 3 H observations, supporting the importance of flow path mixing in this bedrock system. Our findings that the fractured bedrock hosts groundwater with a mixture of residence times ranging from decades to millennia suggest variable recharge dynamics and flow path mixing along the hillslope and highlight the importance of characterizing groundwater systems with observations that are sensitive to transport over a broad range of residence times.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hydrobiogeochemical Controls on the Delivery of Dissolved Organic Matter to Boreal Headwater Streams

Understanding controls on the delivery of dissolved organic matter (DOM) from terrestrial to aquatic systems is key for constraining carbon balances and fluxes in boreal headwater catchments. The largest export of DOM generally occurs during intense periods of water delivery to the landscape (e.g., rain events); however, the timing and magnitude of DOM is complicated by geomorphological, hydrometeorological, and biogeochemical variability. To better assess mechanisms controlling the delivery of DOM to boreal streams, DOM was investigated in two morphologically distinct catchments of an experimental forest watershed in western Newfoundland, one dominated by low relief wetlands and the other by steep hillslopes. The DOM was compared during two fall storms of similar magnitude but contrasting moisture conditions during the autumn transition. Variable concentrations and optical character highlighted antecedent conditions are important for the delivery of DOM to boreal streams. Concentration-discharge relationships with stream hydrograph separation analysis suggest preferential flowpaths through shallow mineral horizons as a key pathway for delivery of DOM in the hillslope dominated catchment compared to rapid input from near-stream wetlands in the low relief catchment. Loss of DOM via preferential pathways suggests an additional mechanism for surface soil carbon loss not directly informed by lateral flow and soil hydrology, which has important implications for our conceptualization and modeling of carbon fluxes in boreal systems. These results suggest watershed scale responses of DOM in boreal headwater catchments are complex, and better linkages between catchment scale hydrology and landscape hydrobiogeochemistry are required to constrain terrestrial to aquatic carbon fluxes in boreal landscapes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Constraints and Drivers of Dissolved Fluxes of Pyrogenic Carbon in Soil and Freshwater Systems: A Global Review and Meta‐Analysis

Abstract Pyrogenic carbon (PyC) is a significant component of the global soil carbon pool due to its longer environmental persistence than other soil organic matter components. Despite PyC's persistence in soil, recent work has indicated that it is susceptible to loss processes such as mineralization and leaching, with the significance and magnitude of these largely unknown at the hillslope and watershed scales. We present a review of the work concerning dissolved PyC transport in soil and freshwater. Our analysis found that the primary environmental controls on dissolved PyC (dPyC) transport are the formation conditions and quality of the PyC itself, with longer and higher temperature charring conditions leading to less transport of dPyC. While correlations between dPyC and dissolved organic carbon in rivers and other pools are frequently reported, the slope of these correlations was pool‐dependent (i.e., soil‐water, precipitation, lakes, streams, rivers), suggesting site‐specific environmental controls. However, the lack of consistency in analytical techniques and sample preparation remains a major challenge to quantifying environmental controls on dPyC fluxes. We propose that future research should focus on the following: (a) consistency in methodological approaches, (b) more quantitative measures of dPyC in pools and fluxes from soils to streams, (c) turnover times of dPyC in soils and aquatic systems, and (d) improved understanding of how mechanisms controlling the fate of dPyC in dynamic post‐fire landscapes interact. With more refined quantitative information about the controls on dPyC transport at the hillslope and landscape scale, we can increase the accuracy and utility of global carbon models.

58 GEOSCIENCES↗