Core Ideas: Development of a 300‐km ² mountainous headwater testbed began in 2016 in the East River. The testbed can be used to explore how watershed changes impact downgradient water availability and quality. System‐of‐system, scale‐adaptive approaches can potentially improve watershed dynamics simulation. We have new approaches to monitor and simulate water partitioning and system responses. The East River watershed has been developed as a “community” testbed. Extreme weather, fires, and land use and climate change are significantly reshaping interactions within watersheds throughout the world. Although hydrological–biogeochemical interactions within watersheds can impact many services valued by society, uncertainty associated with predicting hydrology‐driven biogeochemical watershed dynamics remains high. With an aim to reduce this uncertainty, an approximately 300‐km ² mountainous headwater observatory has been developed at the East River, CO, watershed of the Upper Colorado River Basin. The site is being used as a testbed for the Department of Energy supported Watershed Function Project and collaborative efforts. Building on insights gained from research at the “sister” Rifle, CO, site, coordinated studies are underway at the East River site to gain a predictive understanding of how the mountainous watershed retains and releases water, nutrients, carbon, and metals. In particular, the project is exploring how early snowmelt, drought, and other disturbances influence hydrological–biogeochemical watershed dynamics at seasonal to decadal timescales. A system‐of‐systems perspective and a scale‐adaptive simulation approach, involving the combined use of archetypal watershed subsystem “intensive sites” are being tested at the site to inform aggregated watershed predictions of downgradient exports. Complementing intensive site hydrological, geochemical, geophysical, microbiological, geological, and vegetation datasets are long‐term, distributed measurement stations and specialized experimental and observational campaigns. Several recent research advances provide insights about the intensive sites as well as aggregated watershed behavior. The East River “community testbed” is currently hosting scientists from more than 30 institutions to advance mountainous watershed methods and understanding.

We use 3‐D high‐resolution reactive transport modeling to investigate whether the spatial distribution of organic‐carbon‐rich and chemically reduced sediments located in the riparian zone and temporal variability in groundwater flow direction impact the formation and distribution of nitrogen hot spots (regions that exhibit higher reaction rates when compared to other locations nearby) and hot moments (times that exhibit high reaction rates as compared to longer intervening time periods) within the Rifle floodplain in Colorado. Groundwater flows primarily toward the Colorado River from the floodplain but changes direction at times of high river stage. The result is that oxic river water infiltrates the Rifle floodplain during these relatively short‐term events. Simulation results indicate that episodic rainfall in the summer season leads to the formation of nitrogen hot moments associated with Colorado River rise and resulting river infiltration into the floodplain. The results further demonstrate that the naturally reduced zones (NRZs) present in sediments of the Rifle floodplain have a higher potential for nitrate removal, approximately 70% greater than non‐NRZs for typical hydrological conditions. During river water infiltration, nitrate reduction capacity remains the same within the NRZs, however, these conditions impact non‐NRZs to a greater extent (approximately 95% less nitrate removal). Model simulations indicate chemolithoautotrophs are primarily responsible for the removal of nitrate in the Rifle floodplain. These nitrogen hot spots and hot moments are sustained by microbial respiration and the chemolithoautotrophic oxidation of reduced minerals in the riparian zone.

We develop a multiscale zonation approach to characterize the spatial variability of Arctic polygonal ground geomorphology and to assess the relative controls of these elements on land surface and subsurface properties and carbon fluxes. Working within an ice wedge polygonal region near Barrow, Alaska, we consider two scales of zonation: polygon features (troughs, centers, and rims of polygons) that are nested within different polygon types (high, flat, and low centered). In this study, we first delineated polygons using a digital elevation map and clustered the polygons into four types along two transects, using geophysical and kite‐based landscape‐imaging data sets. We extrapolated those data‐defined polygon types to all the polygons over the study site, using the polygon statistics extracted from the digital elevation map. Based on the point measurements, we characterized the distribution of vegetation, hydrological, thermal, and geochemical properties, as well as carbon fluxes, all as a function of polygon types and polygon features. Results show that nested polygon geomorphic zonation—polygon types and polygon features—can be used to represent distinct distributions of carbon fluxes and associated properties, as well as covariability among those properties. Importantly, the results indicate that polygon types have more power to explain the variations in those properties than polygon features. The approach is expected to be useful for improved system understanding, site characterization, and parameterization of numerical models aimed at predicting ecosystem feedbacks to the climate.