Search

Mountainous watersheds are characterized by variability in functional traits, including vegetation, topography, geology, and geomorphology, which determine nitrogen (N) retention, and release. Coal Creek and East River are two contrasting catchments within the Upper Colorado River Basin that differ markedly in total nitrate (NO3−) export. The East River has a diverse vegetation cover, and sinuous floodplains, and is underlain by N-rich marine shale. At 0.21 ± 0.14 kg ha−1 yr−1, the East River exports ∼3.5 times more NO3− relative to the conifer-dominated Coal Creek (0.06 ± 0.02 kg ha−1 yr−1). While this can partly be explained by the larger size of the East River, the distinct watershed traits of these two catchments imply different mechanisms controlling the aggregate N-export signal. A causality analysis shows physical and biogenic processes were critical in determining NO3− export from the East River catchment. Stable isotope ratios of NO3− (δ15NNO3 and δ18ONO3) show the East River catchment is a strong hotspot for biogeochemical processing of NO3− at the hillslope soil-saprolite. By contrast, the conifer-dominated Coal Creek retained nearly all atmospherically deposited NO3−, and its export was controlled by catchment hydrological traits (i.e., snowmelt periods and water table depth). The conservative N-cycle within Coal Creek is likely due to the abundance of conifer trees, and smaller riparian regions, retaining more NO3− overall and reduced processing prior to export. This study highlights the value of integrating isotope systematics to link watershed functional traits to mechanisms of watershed element retention and release.

In complex terrain, non-parallel surfaces receive emitted radiation from adjacent surfaces. Qualitatively, where surface skin temperatures and lower tropospheric temperature and humidity are not uniform, the downwelling longwave radiation (DLR) will be determined not just by radiation from the atmosphere above a given location, but also by adjacent surface temperatures. We quantify this three-dimensional longwave radiative effect over the Upper Colorado River Basin in clear-sky conditions by calculating surface DLR with observed land-surface temperatures from ECOSTRESS. We find that this effect is due to terrain-subtended sky-view and represents ∼22% of the surface longwave flux, rising to ∼28% and ∼24% in the East and Southeast of the Basin, respectively, and can be >50% in extreme cases. The common omission of this effect in atmospheric radiation models leads to an underestimation of DLR in complex terrain, especially at higher elevations, which has significant implications for mountainous ecohydrology simulations.

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.

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 50 years 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.

In recent years, the availability of airborne imaging spectroscopy (hyperspectral) data has expanded dramatically. The high spatial and spectral resolution of these data uniquely enable spatially explicit ecological studies including species mapping, assessment of drought mortality and foliar trait distributions. However, we have barely begun to unlock the potential of these data to use direct mapping of vegetation characteristics to infer subsurface properties of the critical zone. To assess their utility for Earth systems research, imaging spectroscopy data acquisitions require integration with large, coincident ground-based datasets collected by experts in ecology and environmental and Earth science. Without coordinated, well-planned field campaigns, potential knowledge leveraged from advanced airborne data collections could be lost. Despite the growing importance of this field, documented methods to couple such a wide variety of disciplines remain sparse. We coordinated the first National Ecological Observatory Network Airborne Observation Platform (AOP) survey performed outside of their core sites, which took place in the Upper East River watershed, Colorado. Extensive planning for sample tracking and organization allowed field and flight teams to update the ground-based sampling strategy daily. This enabled collection of an extensive set of physical samples to support a wide range of ecological, microbiological, biogeochemical and hydrological studies. We present a framework for integrating airborne and field campaigns to obtain high-quality data for foliar trait prediction and document an archive of coincident physical samples collected to support a systems approach to ecological research in the critical zone. This detailed methodological account provides an example of how a multi-disciplinary and multi-institutional team can coordinate to maximize knowledge gained from an airborne survey, an approach that could be extended to other studies. The coordination of imaging spectroscopy surveys with appropriately timed and extensive field surveys, along with high-quality processing of these data, presents a unique opportunity to reveal new insights into the structure and dynamics of the critical zone. To our knowledge, this level of co-aligned sampling has never been undertaken in tandem with AOP surveys and subsequent studies utilizing this archive will shed considerable light on the breadth of applications for which imaging spectroscopy data can be leveraged.

In the headwater catchments of the Rocky Mountains, plant productivity and its dynamics are largely dependent upon water availability, which is influenced by changing snowmelt dynamics associated with climate change. Understanding and quantifying the interactions between snow, plants and soil moisture is challenging, since these interactions are highly heterogeneous in mountainous terrain, particularly as they are influenced by microtopography within a hillslope. Recent advances in satellite remote sensing have created an opportunity for monitoring snow and plant dynamics at high spatiotemporal resolutions that can capture microtopographic effects. In this study, we investigate the relationships among topography, snowmelt, soil moisture and plant dynamics in the East River watershed, Crested Butte, Colorado, based on a time series of 3-meter resolution PlanetScope normalized difference vegetation index (NDVI) images. To make use of a large volume of high-resolution time-lapse images (17 images total), we use unsupervised machine learning methods to reduce the dimensionality of the time lapse images by identifying spatial zones that have characteristic NDVI time series. We hypothesize that each zone represents a set of similar snowmelt and plant dynamics that differ from other identified zones and that these zones are associated with key topographic features, plant species and soil moisture. We compare different distance measures (Ward and complete linkage) to understand the effects of their influence on the zonation map. Results show that the identified zones are associated with particular microtopographic features; highly productive zones are associated with low slopes and high topographic wetness index, in contrast with zones of low productivity, which are associated with high slopes and low topographic wetness index. The zones also correspond to particular plant species distributions; higher forb coverage is associated with zones characterized by higher peak productivity combined with rapid senescence in low moisture conditions, while higher sagebrush coverage is associated with low productivity and similar senescence patterns between high and low moisture conditions. In addition, soil moisture probe and sensor data confirm that each zone has a unique soil moisture distribution. This cluster-based analysis can tractably analyze high-resolution time-lapse images to examine plant-soil-snow interactions, guide sampling and sensor placements and identify areas likely vulnerable to ecological change in the future.

This study aims to investigate the microtopographic controls that dictate the heterogeneity of plant communities in a mountainous floodplain-hillslope system, using remote sensing and surface geophysical techniques. Working within a lower montane floodplain-hillslope study site (750 m × 750 m) in the Upper Colorado River Basin, we developed a new data fusion framework, based on machine learning and feature engineering, that exploits remote sensing optical and light detection and ranging (LiDAR) data to estimate the distribution of key plant meadow communities at submeter resolution. We collected surface electrical resistivity tomography data to explore the variability in soil properties along a floodplain-hillslope transect at 0.50-m resolution and extracted LiDAR-derived metrics to model the rapid change in microtopography. We then investigated the covariability among the estimated plant community distributions, soil information, and topographic metrics. Results show that our framework estimated the distribution of nine plant communities with higher accuracy (87% versus 80% overall; 85% versus 60% for shrubs) compared to conventional classification approaches. Analysis of the covariabilities reveals a strong correlation between plant community distribution, soil electric conductivity, and slope, indicating that soil moisture is a primary control on heterogeneous spatial distribution. At the same time, microtopography plays an important role in creating particular ecosystem niches for some of the communities. Such relationships could be exploited to provide information about the spatial variability of soil properties. This highly transferable framework can be employed within long-term monitoring to capture community-specific physiological responses to perturbations, offering the possibility of bridging local plot-scale observations with large landscape monitoring.

Surface fluxes and states can recur and remain consistent across various spatial and temporal scales, forming space‐time patterns. Quantifying and understanding the observed patterns is desirable, as they provide information about the dynamics of the processes involved. This study introduces the empirical spatio‐temporal covariance function and a corresponding parametric covariance function as tools to identify and characterize spatio‐temporal patterns in remotely sensed fields. The method is demonstrated using 2 km hourly GOES‐16/17 land surface temperature (LST) data over the Contiguous United States by splitting the area into 1.0° × 1.0° domains. The summer day‐time LST ESTCFs for 2018 to 2022 are derived for each domain, and a parametric covariance model is fitted. Clustering analysis is applied to detect areas with similar spatio‐temporal LST patterns. Six main zones within CONUS are identified and characterized based on their variance and temporal and spatial characteristic length scales (i.e., scales for which the temperature variations are temporally and spatially related), respectively: (a) Eastern plains with 3 K2, ∼6 hr, and 0.15°, (b) Gulf of California with 60 K2, ∼8 hr, and 0.34°, (c) mountains and coasts transition 1 with 16 K2, ∼11 hr, and 0.25°, (d) central US, Midwest, and South cities with 5.5 K2, ∼8 hr, and ∼0.2°, (e) mountains and coasts transition 2 with ∼10 K2, ∼8 hr, and 0.2°, and (f) largest mountains and coastlines with ∼19 K2, ∼13 hr, and 0.3°. The tools introduced provide a pathway to formally identify and summarize the spatio‐temporal patterns observed in remotely sensed fields and relate those to more complex processes within the Soil‐Vegetation‐Atmosphere System. Plain Language Summary: Specific processes on Earth's surface, like heat and water fluxes and air movement, often exhibit coherent patterns across space and time. Figuring out where these patterns occur can be helpful as they can aid to explain how the underlying physical processes work. This research introduces a new tool to find and describe these patterns and applies it to temperature data derived from satellites. The United States is divided into smaller areas, and the developed tool is used to analyze the land temperature data within those areas; then, regions with similar temperature patterns are identified. These results help to determine that certain landscape features, like coastlines, mountains, and cities, influence how temperature patterns behave in space, time and both. These tools help us to recognize and describe the patterns we see in satellite observations and ultimately, connect them to more complicated processes happening in the Soil‐Vegetation‐Atmosphere System. Key Points: Space‐time covariance can summarize the space‐time structure of time‐varying spatial fields of surface fluxes and statesThe observed space‐time covariance of land surface temperature (LST) reveals variability in its characteristic length and time scales over the United StatesClustering of the fitted LST space‐time covariance parameters reveals the role of environmental characteristics in the observed patterns [ABSTRACT FROM AUTHOR]

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 hydrologydriven biogeochemical watershed dynamics remains high. With an aim to reduce this uncertainty, an approximately 300-km2 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.

Science data gateways are effective in providing complex science data collections to the world-wide user communities. In this paper we describe a gateway for the Advanced Simulation Capability for Environmental Management (ASCEM) framework. Built on top of established web service technologies, the ASCEM data gateway is specifically designed for environmental modeling applications. Its key distinguishing features include (1) handling of complex spatiotemporal data, (2) offering a variety of selective data access mechanisms, (3) providing state-of-the-art plotting and visualization of spatiotemporal data records, and (4) integrating seamlessly with a distributed workflow system using a RESTful interface. ASCEM project scientists have been using this data gateway since 2011.

Shallow permafrost distribution and characteristics are important for predicting ecosystem feedbacks to a changing climate over decadal to century timescales because they can drive active layer deepening and land surface deformation, which in turn can significantly affect hydrologic and biogeochemical responses, including greenhouse gas dynamics. As part of the U.S. Department of Energy Next-Generation Ecosystem Experiments- Arctic, we have investigated shallow Arctic permafrost characteristics at a site in Barrow, Alaska, with the objective of improving our understanding of the spatial distribution of shallow permafrost, its associated properties, and its links with landscape microtopography. To meet this objective, we have acquired and integrated a variety of information, including electric resistance tomography data, frequency-domain electromagnetic induction data, laboratory core analysis, petrophysical studies, high-resolution digital surface models, and color mosaics inferred from kite-based landscape imaging. The results of our study provide a comprehensive and high-resolution examination of the distribution and nature of shallow permafrost in the Arctic tundra, including the estimation of ice content, porosity, and salinity. Among other results, porosity in the top 2 m varied between 85% (besides ice wedges) and 40%, and was negatively correlated with fluid salinity. Salinity directly influenced ice and unfrozen water content and indirectly influenced the soil organic matter content. A relatively continuous but depth-variable increase in salinity led to a partially unfrozen saline layer (cryopeg) located below the top of the permafrost. The cryopeg environment could lead to year-round microbial production of greenhouse gases. Results also indicated a covariability between topography and permafrost characteristics including icewedge and salinity distribution. In addition to providing insight about the Arctic ecosystem, through integration of lab-based petrophysical results with field data, this study also quantified the key controls on electric resistivity at this Arctic permafrost site, including salinity, porosity, water content, ice content, soil organic matter content, and lithologic properties.

Abstract: Microtopographic variation that develops among features (troughs, rims, and centers) within polygonal landforms of coastal arctic tundra strongly affects movement of surface water and snow and thereby affects soil water contents (θ) and active layer depth (ALD). Spatial variation in ALD among these features may exceed interannual variation in ALD caused by changes in climate and so needs to be represented in projections of changes in arctic ALD. In this study, increases in near‐surface θ with decreasing surface elevation among polygon features at the Barrow Experimental Observatory (BEO) were modeled from topographic effects on redistribution of surface water and snow and from lateral water exchange with a subsurface water table during a model run from 1981 to 2015. These increases in θ caused increases in thermal conductivity that in turn caused increases in soil heat fluxes and hence in ALD of up to 15 cm with lower versus higher surface elevation which were consistent with increases measured at BEO. The modeled effects of θ caused interannual variation in maximum ALD that compared well with measurements from 1985 to 2015 at the Barrow Circumpolar Active Layer Monitoring (CALM) site (R2 = 0.61, RMSE = 0.03 m). For higher polygon features, interannual variation in ALD was more closely associated with annual precipitation than mean annual temperature, indicating that soil wetting from increases in precipitation may hasten permafrost degradation beyond that caused by soil warming from increases in air temperature. This degradation may be more rapid if increases in precipitation cause sustained wetting in higher features. [ABSTRACT FROM AUTHOR]