Search

Spatial variation in snow depth is a main driver of heterogeneity in discontinuous permafrost landscapes, exerting a strong control on thermal and hydrological processes, vegetation dynamics, and carbon cycling. Topography and vegetation are understood to play an important role in driving variation in snow depth, but complex morphology often impedes efforts to disentangle these drivers. Maps of ground, vegetation and snow surface elevation were collected using an Unmanned Aerial Vehicle (UAV) over multiple years across a watershed on the Seward Peninsula in Alaska. Here, we quantify drivers of snow depth variation using the inferred maps of snow depth during peak snow accumulation in 2019 and 2022 and collocated ground surface elevation and vegetation height. A novel approach to extract microtopographic information from complex landscape morphologies is used to classify different features (e.g. drainage paths, risers and terraces, thermokarst patterned ground) and characterize their relationships with snow depth variation. A simple model developed using topographic information alone is shown to correlate strongly with local snow depth variation where vegetation height is low. We build a machine learning model to quantify snow trapping by shrub canopies in the watershed and show that snow trapping can be characterized by an exponential function of canopy height above snow (RMSE = 0.12 m, R2 = 0.5). Finally, we demonstrate that relationships between microtopography, vegetation height, and snow depth hold in years of deep and shallow snowpack. These results can be applied to improve representation of heterogeneity and vegetation-snow feedbacks in Earth System Models and to increase the spatial resolution of pan-arctic estimates of snow depth.

Measuring soil and snow temperature with high vertical and lateral resolution is critical for advancing the predictive understanding of thermal and hydro-biogeochemical processes that govern the behavior of environmental systems. Vertically resolved soil temperature measurements enable the estimation of soil thermal regimes, frozen-/thawed-layer thickness, thermal parameters, and heat and/or water fluxes. Similarly, they can be used to capture the snow depth and the snowpack thermal parameters and fluxes. However, these measurements are challenging to acquire using conventional approaches due to their total cost, their limited vertical resolution, and their large installation footprint. This study presents the development and validation of a novel distributed temperature profiling (DTP) system that addresses these challenges. The system leverages digital temperature sensors to provide unprecedented, finely resolved depth profiles of temperature measurements with flexibility in system geometry and vertical resolution. The integrated miniaturized logger enables automated data acquisition, management, and wireless transfer. A novel calibration approach adapted to the DTP system confirms the factory-assured sensor accuracy of ±0.1 ∘C and enables improving it to ±0.015 ∘C. Numerical experiments indicate that, under normal environmental conditions, an additional error of 0.01 % in amplitude and 70 s time delay in amplitude for a diurnal period can be expected, owing to the DTP housing. We demonstrate the DTP systems capability at two field sites, one focused on understanding how snow dynamics influence mountainous water resources and the other focused on understanding how soil properties influence carbon cycling. Results indicate that the DTP system reliably captures the dynamics in snow depth and soil freezing and thawing depth, enabling advances in understanding the intensity and timing in surface processes and their impact on subsurface thermohydrological regimes. Overall, the DTP system fulfills the needs for data accuracy, minimal power consumption, and low total cost, enabling advances in the multiscale understanding of various cryospheric and hydro-biogeochemical processes.

Water resources, including groundwater and prominent rivers worldwide, are under duress because of excessive contaminant and nutrient loads. To help mitigate this problem, the United States Department of Energy (DOE) has supported research since the late 1980s to improve our fundamental knowledge of processes that could be used to help clean up challenging subsurface problems. Problems of interest have included subsurface radioactive waste, heavy metals, and metalloids (e.g. uranium, mercury, arsenic). Research efforts have provided insights into detailed groundwater biogeochemical process coupling and the resulting geochemical exports of metals and nutrients to surrounding environments. Recently, an increased focus has been placed on constraining the exchanges and fates of carbon and nitrogen within and across bedrock to canopy compartments of a watershed and in river–floodplain settings, because of their important role in driving biogeochemical interactions with contaminants and the potential of increased fluxes under changing precipitation regimes, including extreme events. While reviewing the extensive research that has been conducted at DOE’s representative sites and testbeds (such as the Oyster Site in Virginia, Savannah River Site in South Carolina, Oak Ridge Reservation in Tennessee, Hanford in Washington, Nevada National Security Site in Nevada, Riverton in Wyoming, and Rifle and East River in Colorado), this review paper explores the nature and distribution of contaminants in the surface and shallow subsurface (i.e. the critical zone) and their interactions with carbon and nitrogen dynamics. We also describe state-of-the-art, scale-aware characterization approaches and models developed to predict contaminant fate and transport. The models take advantage of DOE leadership-class high-performance computers and are beginning to incorporate artificial intelligence approaches to tackle the extreme diversity of hydro-biogeochemical processes and measurements. Recognizing that the insights and capability developments are potentially transferable to many other sites, we also explore the scientific implications of these advances and recommend future research directions.

Soil thickness plays a central role in the interactions between vegetation, soils, and topography, where it controls the retention and release of water, carbon, nitrogen, and metals. However, mapping soil thickness, here defined as the mobile regolith layer, at high spatial resolution remains challenging. Here, we develop a hybrid model that combines a process-based model and empirical relationships to estimate the spatial heterogeneity of soil thickness with fine spatial resolution (0.5 m). We apply this model to two aspects of hillslopes (southwest- and northeast-facing, respectively) in the East River watershed in Colorado. Two independent measurement methods – auger and cone penetrometer – are used to sample soil thickness at 78 locations to calibrate the local value of unconstrained parameters within the hybrid model. Sensitivity analysis using the hybrid model reveals that the diffusion coefficient used in hillslope diffusion modeling has the largest sensitivity among all input parameters. In addition, our results from both sampling and modeling show that, in general, the northeast-facing hillslope has a deeper soil layer than the southwest-facing hillslope. By comparing the soil thickness estimated between a machine-learning approach and this hybrid model, the hybrid model provides higher accuracy and requires less sampling data. Modeling results further reveal that the southwest-facing hillslope has a slightly faster surface soil erosion rate and soil production rate than the northeast-facing hillslope, which suggests that the relatively less dense vegetation cover and drier surface soils on the southwest-facing slopes influence soil properties. With seven parameters in total for calibration, this hybrid model can provide a realistic soil thickness map with a relatively small amount of sampling dataset comparing to machine-learning approach. Integrating process-based modeling and statistical analysis not only provides a thorough understanding of the fundamental mechanisms for soil thickness prediction but also integrates the strengths of both statistical approaches and process-based modeling approaches.

Isotopic information from 81 snowpits was collected over a 5-year period in a large, Colorado watershed. Data spans gradients in elevation, aspect, vegetation, and seasonal climate. They are combined with overlapping campaigns for water isotopes in precipitation and snowmelt, and a land-surface model for detailed estimates of snowfall and climate at sample locations. Snowfall isotopic inputs, describe the majority of δ18O snowpack variability. Aspect is a secondary control, with slightly more enriched conditions on east and north facing slopes. This is attributed to preservation of seasonally enriched snowfall and vapour loss in the early winter. Sublimation, expressed by decreases in snowpack d-excess in comparison to snowfall contributions, increases at low elevation and when seasonal temperature and solar radiation are high. At peak snow accumulation, post-depositional fractionation appears to occur in the top 25 ± 14% of the snowpack due to melt-freeze redistribution of lighter isotopes deeper into the snowpack and vapour loss to the atmosphere during intermittent periods of low relative humidity and high windspeed. Relative depth of fractionation increases when winter daytime temperatures are high and winter precipitation is low. Once isothermal, snowpack isotopic homogenization and enrichment was observed with initial snowmelt isotopically depleted in comparison to snowpack and enriching over time. The rate of δ18O increase (d-excess decrease) in snowmelt was 0.02‰ per day per 100-m elevation loss. Isotopic data suggests elevation dictates snowpack and snowmelt evolution by controlling early snow persistence (or absence), isotopic lapse rates in precipitation and the ratio of energy to snow availability. Hydrologic tracer studies using stable water isotopes in basins of large topographic relief will require adjustment for these elevational controls to properly constrain stream water sourcing from snowmelt.

Understanding the partitioning of snow and rain contributing to either catchment streamflow or evapotranspiration (ET) is of critical relevance for water management in response to climate change. To investigate this partitioning, we use endmember splitting and mixing analyses based on stable isotope (18O) data from nine headwater catchments in the East River, Colorado. Our results show that one third of the snow partitions to ET and 13% of the snowmelt sustains summer streamflow. Only 8% of the rainfall contributes to the summer streamflow, because most of the rain (67%) partitions to ET. The spatial variability of precipitation partitioning is mainly driven by aspect and tree cover across the sub-catchments. Catchments with higher tree cover have a higher share of snow becoming ET, resulting in less snow in summer streamflow. Summer streamflow did not contain more rain with higher rainfall sums, but more rain was taken up in ET.

Atmospheric nitrous oxide contributes directly to global warming, yet models of the nitrogen cycle do not account for bedrock, the largest pool of terrestrial nitrogen, as a source of nitrous oxide. Although it is known that release rates of nitrogen from bedrock are large, there is an incomplete understanding of the connection between bedrock-hosted nitrogen and atmospheric nitrous oxide. Here, we quantify nitrogen fluxes and mass balances at a hillslope underlain by marine shale. We found that, at this site, bedrock weathering contributes 78% of the subsurface reactive nitrogen, while atmospheric sources (commonly regarded as the sole sources of reactive nitrogen in pristine environments) account for only the remaining 22%. About 56% of the total subsurface reactive nitrogen denitrifies, including 14% emitted as nitrous oxide. The remaining reactive nitrogen discharges in porewaters to a floodplain where additional denitrification probably occurs. We also found that the release of bedrock nitrogen occurs primarily within the zone of the seasonally fluctuating water table and suggest that the accumulation of nitrate in the vadose zone, often attributed to fertilization and soil leaching, may also include contributions from weathered nitrogen-rich bedrock. Our hillslope study suggests that, under oxygenated and moisture-rich conditions, weathering of deep, nitrogen-rich bedrock makes an important contribution to the nitrogen cycle. Weathering of deep bedrock releases reactive nitrogen into the subsurface, which contributes to the flux of nitrous oxide to the atmosphere, according to a field study that combines soil, rock and groundwater data within a river catchment.

Understanding the interactions among agricultural processes, soil, and plants is necessary for optimizing crop yield and productivity. This study focuses on developing effective monitoring and analysis methodologies that estimate key soil and plant properties. These methodologies include data acquisition and processing approaches that use unmanned aerial vehicles (UAVs) and surface geophysical techniques. In particular, we applied these approaches to a soybean farm in Arkansas to characterize the soil–plant coupled spatial and temporal heterogeneity, as well as to identify key environmental factors that influence plant growth and yield. UAV-based multitemporal acquisition of high-resolution RGB (red–green–blue) imagery and direct measurements were used to monitor plant height and photosynthetic activity. We present an algorithm that efficiently exploits the high-resolution UAV images to estimate plant spatial abundance and plant vigor throughout the growing season. Such plant characterization is extremely important for the identification of anomalous areas, providing easily interpretable information that can be used to guide near-real-time farming decisions. Additionally, high-resolution multitemporal surface geophysical measurements of apparent soil electrical conductivity were used to estimate the spatial heterogeneity of soil texture. By integrating the multiscale multitype soil and plant datasets, we identified the spatiotemporal co-variance between soil properties and plant development and yield. Our novel approach for early season monitoring of plant spatial abundance identified areas of low productivity controlled by soil clay content, while temporal analysis of geophysical data showed the impact of soil moisture and irrigation practice (controlled by topography) on plant dynamics. Our study demonstrates the effective coupling of UAV data products with geophysical data to extract critical information for farm management.

Working within a vineyard in the Pessac Léognan Appellation of Bordeaux, France, this study documents the potential of using simple statistical methods with spatially-resolved and increasingly available electromagnetic induction (EMI) geophysical and normalized difference vegetation index (NDVI) datasets to accurately estimate Bordeaux vineyard soil classes and to quantitatively explore the relationship between vineyard soil types and grapevine vigor. First, co-located electrical tomographic tomography (ERT) and EMI datasets were compared to gain confidence about how the EMI method averaged soil properties over the grapevine rooting depth. Then, EMI data were used with core soil texture and soil-pit based interpretations of Bordeaux soil types (Brunisol, Redoxisol, Colluviosol and Calcosol) to estimate the spatial distribution of geophysically-identified Bordeaux soil classes. A strong relationship (r = 0.75, p

Diverse observational and simulation datasets are needed to understand and predict complex ecosystem behavior over seasonal to decadal and century time-scales. Integration of these datasets poses a major barrier towards advancing environmental science, particularly due to differences in the structure and formats of data provided by various sources. Here, we describe BASIN-3D (Broker for Assimilation, Synthesis and Integration of eNvironmental Diverse, Distributed Datasets), a data integration framework designed to dynamically retrieve and transform heterogeneous data from different sources into a common format to provide an integrated view. BASIN-3D enables users to adopt a standardized approach for data retrieval and avoid customizations for the data type or source. We demonstrate the value of BASIN-3D with two use cases that require integration of data from regional to watershed spatial scales. The first application uses the BASIN-3D Python library to integrate time-series hydrological and meteorological data to provide standardized inputs to analytical and machine learning codes in order to predict the impacts of hydrological disturbances on large river corridors of the United States. The second application uses the BASIN-3D Django framework to integrate diverse time-series data in a mountainous watershed in East River, Colorado, United States to enable scientific researchers to explore and download data through an interactive web portal. Thus, BASIN-3D can be used to support data integration for both web-based tools, as well as data analytics using Python scripting and extensions like Jupyter notebooks. The framework is expected to be transferable to and useful for many other field and modeling studies.

Land-atmosphere carbon exchange is known to be extremely heterogeneous in arctic ice-wedge polygonal tundra regions. In this study, a Kalman filter-based method was developed to estimate the spatio-temporal dynamics of daytime average net ecosystem exchange (NEEday) at 0.5-m resolution over a 550 m by 700 m study site. We integrated multi-scale, multi-type datasets, including normalized difference vegetation indices (NDVIs) obtained from a novel automated mobile sensor system (or tram system) and a greenness index map obtained from airborne imagery. We took advantage of the significant correlations between NDVI and NEEday identified based on flux chamber measurements. The weighted average of the estimated NEEday within the flux-tower footprint agreed with the flux tower data in term of its seasonal dynamics. We then evaluated the spatial variability of the growing season average NEEday, as a function of polygon geomorphic classes, i.e., the combination of polygon types—which are known to present different degradation stages associated with permafrost thaw—and microtopographic features (i.e., troughs, centers and rims). Our study suggests the importance of considering microtopographic features and their spatial coverage in computing spatially aggregated carbon exchange.

Author(s): Kakalia, Z; Varadharajan, C; Alper, E; Brodie, EL; Burrus, M; Carroll, RWH; Christianson, DS; Dong, W; Hendrix, VC; Henderson, M; Hubbard, SS; Johnson, D; Versteeg, R; Williams, KH; Agarwal, DA | Abstract: The U.S. Department of Energy's (DOE) Colorado East River Community Observatory (ER) in the Upper Colorado River Basin was established in 2015 as a representative mountainous, snow-dominated watershed to study hydrobiogeochemical responses to hydrological perturbations in headwater systems. The ER is characterized by steep elevation, geologic, hydrologic and vegetation gradients along floodplain, montane, subalpine, and alpine life zones, which makes it an ideal location for researchers to understand how different mountain subsystems contribute to overall watershed behaviour. The ER has both long-term and spatially-extensive observations and experimental campaigns carried out by the Watershed Function Scientific Focus Area (SFA), led by Lawrence Berkeley National Laboratory, and researchers from over 30 organizations who conduct cross-disciplinary process-based investigations and modelling of watershed behaviour. The heterogeneous data generated at the ER include hydrological, genomic, biogeochemical, climate, vegetation, geological, and remote sensing data, which combined with model inputs and outputs comprise a collection of datasets and value-added products within a mountainous watershed that span multiple spatiotemporal scales, compartments, and life zones. Within 5 years of collection, these datasets have revealed insights into numerous aspects of watershed function such as factors influencing snow accumulation and melt timing, water balance partitioning, and impacts of floodplain biogeochemistry and hillslope ecohydrology on riverine geochemical exports. Data generated by the SFA are managed and curated through its Data Management Framework. The SFA has an open data policy, and over 70 ER datasets are publicly available through relevant data repositories. A public interactive map of data collection sites run by the SFA is available to inform the broader community about SFA field activities. Here, we describe the ER and the SFA measurement network, present the public data collection generated by the SFA and partner institutions, and highlight the value of collecting multidisciplinary multiscale measurements in representative catchment observatories.

Soil temperature has been recognized as a property that strongly influences a myriad of hydro-biogeochemical processes and reflects how various properties modulate the soil thermal flux. In spite of its importance, our ability to acquire soil temperature data with high spatial and temporal resolution and coverage is limited because of the high cost of equipment, the difficulties of deployment, and the complexities of data management. Here we propose a new strategy that we call distributed temperature profiling (DTP) for improving the characterization and monitoring near-surface thermal properties through the use of an unprecedented number of laterally and vertically distributed temperature measurements. We developed a prototype DTP system, which consists of inexpensive, low-impact, low-power, and vertically resolved temperature probes that independently and autonomously record soil temperature. The DTP system concept was tested by moving sequentially the system across the landscape to identify near-surface permafrost distribution in a discontinuous permafrost environment near Nome, Alaska, during the summertime. Results show that the DTP system enabled successful acquisition of vertically resolved profiles of summer soil temperature over the top 0.8 m at numerous locations. DTP also enabled high-resolution identification and lateral delineation of near-surface permafrost locations from surrounding zones with no permafrost or deep permafrost table locations overlain by a perennially thawed layer. The DTP strategy overcomes some of the limitations associated with – and complements the strengths of – borehole-based soil temperature sensing as well as fiber-optic distributed temperature sensing (FO-DTS) approaches. Combining DTP data with co-located topographic and vegetation maps obtained using unmanned aerial vehicle (UAV) and electrical resistivity tomography (ERT) data allowed us to identify correspondences between surface and subsurface property distribution and in particular between topography, vegetation, shallow soil properties, and near-surface permafrost. Finally, the results highlight the considerable value of the newly developed DTP strategy for investigating the significant variability in and complexity of subsurface thermal and hydrological regimes in discontinuous permafrost regions.

Author(s): Tokunaga, TK; Wan, J; Williams, KH; Brown, W; Henderson, A; Kim, Y; Tran, AP; Conrad, ME; Bill, M; Carroll, RWH; Dong, W; Xu, Z; Lavy, A; Gilbert, B; Carrero, S; Christensen, JN; Faybishenko, B; Arora, B; Siirila-Woodburn, ER; Versteeg, R; Raberg, JH; Peterson, JE; Hubbard, SS | Abstract: Major components of hydrologic and elemental cycles reside underground, where their complex dynamics and linkages to surface waters are obscure. We delineated seasonal subsurface flow and transport dynamics along a hillslope in the Rocky Mountains (USA), where precipitation occurs primarily as winter snow and drainage discharges into the East River, a tributary of the Gunnison River. Hydraulic and geochemical measurements down to 10 m below ground surface supported application of transmissivity feedback of snowmelt to describe subsurface flow and transport through three zones: soil, weathering shale, and saturated fractured shale. Groundwater flow is predicted to depths of at least 176 m, although a shallower limit exists if hillslope-scale hydraulic conductivities are higher than our local measurements. Snowmelt during the high snowpack water year 2017 sustained flow along the weathering zone and downslope within the soil, while negligible downslope flow occurred along the soil during the low snowpack water year 2018. We introduce subsurface concentration-discharge (C-Q) relations for explaining hillslope contributions to C-Q observed in rivers and demonstrate their calculations based on transmissivity fluxes and measured pore water specific conductance and dissolved organic carbon. The specific conductance data show that major ions in the hillslope pore waters, primarily from the weathering and fractured shale, are about six times more concentrated than in the river, indicating hillslope solute loads are disproportionately high, while flow from this site and similar regions are relatively smaller. This methodology is applicable in different representative environments within snow-dominated watersheds for linking their subsurface exports to surface waters.

Water is essential to human health and economic development due to its utilization in sanitation, agriculture, and energy. Supplying water to an expanding world population requires simultaneous consideration of multiple societal sectors competing for limited resources. Water conservation, supply augmentation, distribution, and treatment of contaminants must work in concert to ensure water sustainability. Water is linked to other sectors, and the quantity and quality of water resources are changing. The efficient use of water in agriculture, the largest user of water worldwide, via drip irrigation is described as is the use of energy-intensive reverse osmosis to supplement freshwater supplies. Efforts to manage watersheds and model their responses to severe weather events are discussed along with efforts to improve the predictability of their function. The regional competition for water resources impacts both energy and water supply reliability, which requires that nations balance both for sustainable econ...

Background Biogeochemical exports from watersheds are modulated by the activity of microorganisms that function over micron scales. Here, we tested the hypothesis that meander-bound regions share a core microbiome and exhibit patterns of metabolic potential that broadly predict biogeochemical processes in floodplain soils along a river corridor. Results We intensively sampled the microbiomes of floodplain soils located in the upper, middle, and lower reaches of the East River, Colorado. Despite the very high microbial diversity and complexity of the soils, we reconstructed 248 quality draft genomes representative of subspecies. Approximately one third of these bacterial subspecies was detected across all three locations at similar abundance levels, and ~ 15% of species were detected in two consecutive years. Within the meander-bound floodplains, we did not detect systematic patterns of gene abundance based on sampling position relative to the river. However, across meanders, we identified a core floodplain microbiome that is enriched in capacities for aerobic respiration, aerobic CO oxidation, and thiosulfate oxidation with the formation of elemental sulfur. Given this, we conducted a transcriptomic analysis of the middle floodplain. In contrast to predictions made based on the prominence of gene inventories, the most highly transcribed genes were relatively rare amoCAB and nxrAB (for nitrification) genes, followed by genes involved in methanol and formate oxidation, and nitrogen and CO2 fixation. Within all three meanders, low soil organic carbon correlated with high activity of genes involved in methanol, formate, sulfide, hydrogen, and ammonia oxidation, nitrite oxidoreduction, and nitrate and nitrite reduction. Overall, the results emphasize the importance of sulfur, one-carbon and nitrogen compound metabolism in soils of the riparian corridor. Conclusions The disparity between the scale of a microbial cell and the scale of a watershed currently limits the development of genomically informed predictive models describing watershed biogeochemical function. Meander-bound floodplains appear to serve as scaling motifs that predict aggregate capacities for biogeochemical transformations, providing a foundation for incorporating riparian soil microbiomes in watershed models. Widely represented genetic capacities did not predict in situ activity at one time point, but rather they define a reservoir of biogeochemical potential available as conditions change.

In this study, we develop a watershed zonation approach for characterizing watershed organization and function in a tractable manner by integrating multiple spatial data layers. Recognizing the coupled ecohydrogeological-biogeochemical interactions that occur across bedrock through canopy compartments of a watershed, we hypothesize that (1) suites of above/belowground properties co-varying with each other, (2) hillslopes are representative units for capturing watershed-scale heterogeneity, (3) remote sensing data layers and clustering methods can be used to identify watershed hillslope zones having unique distributions of bedrock-through-canopy properties relative to neighboring parcels, and (4) property suites associated with the identified zones can be used to understand zone-based functions, such as response to early snowmelt or drought, and associated solute exports to the river. We demonstrate this concept using unsupervised clustering methods that synthesizes airborne remote sensing data (LiDAR, hyperspectral, and electromagnetic surveys) along with satellite and streamflow data collected in the East River Watershed, Crested Butte, Colorado, USA. Results show that, (1) hillslope-average elevation and slope are significantly correlated with near-surface bedrock electrical resistivity (top 20 m), (2) elevation and aspect are independent controls on plant and snow signatures, (3) the correlation between hillslope-averaged above- and below- ground properties are significantly higher than pixel-by-pixel correlation and (4) K-means, hierarchical clustering, and Gaussian mixture clustering methods generate similar zonation patterns across the watershed. Using independently collected data, it is shown that the identified zones provide information about zone-based watershed functions, including foresummer drought sensitivity and river nitrogen exports. The approach is expected to be extensible to other sites and generally useful for guiding the selection of hillslope experiment locations and informing model parameterization.

Patterns of watershed nitrogen (N) retention and loss are shaped by how watershed biogeochemical processes retain, biogeochemically transform, and lose incoming atmospheric deposition of N. Loss patterns represented by concentration, discharge, and their associated stream exports are important indicators of integrated watershed N retention behaviors. We examined continental United States (CONUS) scale N deposition (e.g., wet and dry atmospheric deposition), vegetation trends, and stream trends as potential indicators of watershed N-saturation and retention conditions, and how watershed N retention and losses vary over space and time. By synthesizing changes and modalities in watershed nitrogen loss patterns based on stream data from 2200 U.S. watersheds over a 50years record, our work revealed two patterns of watershed N-retention and loss. One was a hysteresis pattern that reflects the integrated influence of hydrology, atmospheric inputs, land-use, stream temperature, elevation, and vegetation. The other pattern was a one-way transition to a new state. We found that regions with increasing atmospheric deposition and increasing vegetation health/biomass patterns have the highest N-retention capacity, become increasingly N-saturated over time, and are associated with the strongest declines in stream N exports—a pattern, that is, consistent across all land cover categories. We provide a conceptual model, validated at an unprecedented scale across the CONUS that links instream nitrogen signals to upstream mechanistic landscape processes. Our work can aid in the future interpretation of in-stream concentrations of DOC and DIN as indicators of watershed N-retention status and integrators of watershed hydrobiogeochemical processes.

Recent studies in snowmelt-dominated catchments have documented changes in nitrogen (N) retention over time, such as declines in watershed exports of N, though there is a limited understanding of the controlling processes driving these trends. Working in the mountainous headwater East River Colorado watershed, we explored the effects of riparian hollows as N-cycling hotspots and as important small-scale controls on observed watershed trends. Using a modeling-based approach informed by remote sensing and in situ observations, we simulated the N-retention capacity of riparian hollows with seasonal and yearly hydrobiogeochemical perturbations imposed as drivers. We then implemented a scaling approach to quantify the relative contribution of riparian hollows to the total river corridor N budget. We found that riparian hollows primarily serve as N sinks, with N-transformation rates significantly limited by periods of enhanced groundwater upwelling and promoted at the onset of rainfall events. Given these observed hydrologic controls, we expect that the nitrate (NO3-) sink capacity of riparian hollows will increase in magnitude with future climatic perturbations, specifically the shift to more frequent rainfall events and fewer snowmelt events, as projected for many mountainous headwater catchments. Our current estimates suggest that while riparian hollows provision ~5–20% of NO3- to the river network, they functionally act as inhibitors to upland NO3- reaching the stream. Our work linking transient hydrological conditions to numerical biogeochemical simulations is an important step in assessing N-retaining features relative to the watershed N budget and better understanding the role of small-scale features within watersheds.

Soil thickness plays a central role in the interactions between vegetation, soils, and topography where it controls the retention and release of water, carbon, nitrogen, and metals. However, mapping soil thickness, here defined as the mobile regolith layer, at high spatial resolution remains challenging. Here, we develop a hybrid model that combines a process-based model and empirical relationships to estimate the spatial heterogeneity of soil thickness with fine spatial resolution (0.5 m). We apply this model to two examples of hillslopes (south-facing and north-facing, respectively) in the East River Watershed in Colorado that validates the effectiveness of the model. Two independent measurement methods – auger and cone penetrometer – are used to sample soil thickness at 78 locations to calibrate the local value of unconstrained parameters within the hybrid model. Sensitivity analysis using the hybrid model reveals that the diffusion coefficient used in hillslope diffusion modelling has the largest sensitivity among all input parameters. In addition, our results from both sampling and modeling show that, in general, the north-facing hillslope has a deeper soil layer than the south-facing hillslope. By comparing the soil thickness estimated between a machine learning approach and this hybrid model, the hybrid model provides higher accuracy and requires less sampling data. Modeling results further reveal that the south-facing hillslope has a slightly faster surface soil erosion rate and soil production rate than the north-facing hillslope, which suggests that the relatively less dense vegetation cover and drier surface soils on the south-facing slopes may influence soil characteristics. With only seven parameters for calibration, this hybrid model can provide a realistic soil thickness map at other study sites by with a relatively small amount of sampling dataset. Integrating process-based modeling and statistical analysis not only provides a thorough understanding of the fundamental mechanisms for soil thickness prediction, but integrates the strengths of both statistical approaches and process-based modeling approaches.

Summary Remote sensing capabilities have revolutionized several fields including geology and vegetation ecology. With these new capabilities to observe important interactions between subsurface hydrobiogeochemical properties and surface biological response, there is vast potential to rigorously test theories of ecosystem function and evolution. For example, the functional trait distributions of plants can be quantified at increasingly broader scales allowing ecological theories to be tested across plots to watersheds and beyond. Geophysical methods now allow us to interrogate how regolith properties vary across the same scales and provide an insight into the factors regulating plant species composition and functional trait distributions. Given the fundamental constraints of water, energy and nutrients are shared across biology, and that plants and microorganisms share a long history of coevolution, it is likely that microbial functional traits distributions follow similarly interpretable and predictable patterns at ecosystem scales. In this presentation I will cover our work combining theory, observation and modeling to advance how microbial traits are deciphered from genomic information, as well as approaches to ground truth these hypotheses. I will also discuss how we combine various sensing approaches, including vegetation, geophysical, and microbial to evaluate patterns of co-evolution across a watershed with consequences for ecosystem hydrobiogeochemistry.

Climate change is causing rapid changes of Arctic ecosystems. Yet, data needed to unravel complex subsurface processes are very rare. Using geophysical and in situ sensing, this study closes an observational gap associated with thermohydrological dynamics in discontinuous permafrost systems. It highlights the impact of vegetation and snow thickness distribution on subsurface thermohydrological properties and processes. Large snow accumulation near tall shrubs insulates the ground and allows for rapid and downward heat flow. Thinner snow pack above graminoid results in surficial freezing and prevents water from infiltrating into the subsurface. Analyzing short‐term disturbances, we found that lateral flow could be a driving factor in talik formation. Interannual measurements show that deep permafrost temperatures increased by about 0.2°C over 2 years. The results, which suggest that snow‐vegetation‐subsurface processes are tightly coupled, will be useful for improving predictions of Arctic feedback to climate change, including how subsurface thermohydrology influences CO2 and CH4 fluxes., Geophysical Research Letters, 48 (6), ISSN:0094-8276, ISSN:1944-8007

Summary In this work, a novel method for obtaining maps about the physiochemical nature of soils from electromagnetic induction (EMI) data is investigated for a French vineyard. Specifically, K means clustering is used to produce maps from both the quadrature and in-phase components of EMI measurements. The commonly discarded in-phase component of EMI measurements exhibits a higher correlation with plant vigor, as measured by NDVI. In addition to the non-intrusive characterization, extensive soil samples are used to explore the correlations and causations that link soil physiochemical properties with EMI and NDVI data. It is demonstrated that the in-phase component of EMI is related to metal concentrations at the site, these properties co-vary with crucial properties related to plant vigor, e.g. cation exchange capacity and organic matter content.

The U.S. Department of Energy’s (DOE) East River community observatory (ER) in the Upper Colorado River Basin was established in 2015 as a representative mountainous, snow-dominated watershed to study hydrobiogeochemical responses to hydrological perturbations in headwater systems. Led by the Watershed Function Science Focus Area (SFA), the ER has both long-term and spatially-extensive observations paired with experimental campaigns. The Watershed Function SFA, led by Berkeley Laboratory, includes researchers from over 30 organizations who conduct cross-disciplinary process-based investigations and mechanistic modeling of watershed behavior in the ER. The data generated at the ER are extremely heterogeneous, and include hydrological, biogeochemical, climate, vegetation, geological, remote sensing, and model data that together comprise an unprecedented collection of data and value-added products within a mountainous watershed, across multiple spatiotemporal scales, compartments, and life zones. Within 5 years of data collection, these datasets have already revealed insights into numerous aspects of watershed function such as factors influencing snow accumulation and melt timing, water balance partitioning, and impacts of floodplain biogeochemistry and hillslope ecohydrology on riverine geochemical exports. Data generated by the SFA are managed and curated through its Data Management Framework. The SFA has an open data policy, and over sixty ER datasets are publicly available through relevant data repositories. A public interactive map of data collection sites run by the SFA is available to inform the broader community about SFA field activities. Here, we describe the ER and the SFA measurement network, present the public data collection generated by the SFA and partner institutions, and highlight the value of collecting multidisciplinary multiscale measurements in representative catchment observatories.

Concentration-discharge (C-Q) relationships have been widely used as “hydrochemical tracers” to determine the variability in riverine solute exports across event, seasonal, annual, and decadal time scales. However, these C-Q relationships are limited to investigating solute transport dynamics at individual sampling stations, such that they create an incomplete understanding of the solute behavior upstream or downstream of the sampling station. Therefore, the objective of this study is to develop, apply and assess a differential C-Q approach that can characterize spatial variability in solute behavior across stations, as well as investigate their controls, by following a different spatial scheme and organizing the river into multiple sections. The differential C-Q approach captures the difference in concentration in a river segment over the difference in discharge, thereby accounting for gains, losses or fractional solute turnover between sampling stations. Using water quality data collected over four water years (2015–2018) in a mountainous headwater catchment of the East River, Colorado, this study compares traditional and differential C-Q relationships in predicting solute behavior between three sampling stations distributed throughout the river. Results from the differential C-Q analysis demonstrate significant differences in solute behavior within upstream vs. downstream reaches of the East River watershed. In particular, the meandering downstream section is marked by significant gains in both groundwater and solute concentrations as opposed to the dilution and the declining trends observed in the high-relief, steep terrain upstream reach. Shale mineralogy was determined to have a major influence on in-stream concentrations pertaining to Ca, DIC, DOC, Mg, Mo, NO3, and SO4. The analyses further revealed that total P concentration in the downstream reach exceeded the U.S. Environmental Protection Agency's desired goal for control of eutrophication (110 ppb). Overall, differential C-Q metrics yield a better understanding of the lateral storage and interactions within catchments than traditional analyses, and holds potential for aiding water quality managers in the identification of critical stream reaches that assimilate harmful chemicals.