Search NASA⌕ Search

SEARCH · Search NASA

Results for “Desert Climate”

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.

89 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Trends and meteorological drivers of extreme daily reservoir evaporation events in the western United States

Extreme daily evaporation from reservoir surfaces can lead to significant short-term water losses, affecting water quality, water supply, and reservoir operation strategies. Historical trends in daily reservoir evaporation events have eluded the scientific and operational communities, largely due to a lack of long-term, consistent data record. This study quantifies trends in extreme daily reservoir evaporation events at 165 major reservoirs located in the western U.S. Here, we use the place-based energy balance and aerodynamic Daily Lake Evaporation Model (DLEM) driven by multiple meteorological data products (RTMA, gridMET, Daymet) to estimate daily evaporation rates at these reservoirs from 1981 to 2022. The results—while are based on different meteorological forcing datasets—consistently indicate that the California, Lower Colorado, and Rio Grande hydrologic regions are more prone to higher daily evaporation extremes. Compared to the rest of western U.S, these three regions also experience a more pronounced increasing trend in the annual maximum daily evaporation rate, at about 0.3 mm day -1 decade -1 during 1981-2022. The results show that heat and dryness are the main drivers to the increasing trend of extreme evaporation, while extreme wind speed is the primary contributor to exceptionally high daily evaporation events across all regions. This phenomenon is particularly prominent in the arid Lower Colorado region, but less significant in the cold and humid Pacific Northwest region. We also find that the correlation between extreme wind speed and extreme evaporation degrades as the time scale increases from daily, to monthly and seasonal. Our findings have strong implications for the pattern and distribution of extreme evaporation events at the western U.S. reservoirs, and illustrate how various drivers influence extreme evaporation across different time scales.

13 HYDRO ENERGY↗
DOE OSTI2022

Ten Years of VIIRS Land Surface Temperature Product Validation

The Visible Infrared Imaging Radiometer Suite (VIIRS) Land Surface Temperature (LST) has been operationally produced for a decade since the Suomi National Polar-orbiting Partnership (SNPP) launched in October 2011. A comprehensive evaluation of its accuracy and precision will be helpful for product users in climate studies and atmospheric models. In this study, the VIIRS LST is validated with ground observations from multiple high-quality radiation networks, including six stations from the Surface Radiation budget (SURFRAD) network, two stations from the Baseline Surface Radiation Network (BSRN), and 13 stations from the Atmospheric Radiation Measurement (ARM) network, to evaluate its performance over various land-cover types. The VNP21A1 LST was validated against the same ground observations as a reference. The results yield a close agreement between the SNPP VIIRS LST and ground LSTs with a bias of -0.4 K and a RMSE of 1.96 K over six SURFRAD sites; a bias of -0.2 K and a RMSE of 1.93 K over two BSRN sites; and a bias of -0.1 K and a RMSE of 1.7 K over the 13 ARM sites. The time series of the LST errors over individual sites indicate seasonal cycles. The data anomaly over the BSRN site in Cabauw and the SURFRAD site in Desert Rock is revealed and discussed in this study. In addition, a method using Landsat-8 data is applied to quantify the heterogeneity level of each ground station and the results provide promising insights. The validation results demonstrate the maturity of the JPSS VIIRS LST products and their readiness for various application studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluation and improvement of the parameterization of aerosol hygroscopicity in global climate models using in-situ surface measurements (Final Report)

Aerosols are tiny particles suspended on the atmosphere that can interact with incoming solar radiation and affect the Earth radiative budget. They do so by scattering and absorbing solar radiation, and these properties vary depending on the aerosol size and chemical composition. Moreover, by taking up water from the surrounding air, hygroscopic aerosol particles will grow in size and change their chemical composition, thus modifying their optical properties (scattering and absorption) and their final impact on radiative forcing calculations. An accurate knowledge of aerosol hygroscopicity is crucial for estimating the net radiative impact of aerosols. We took a three-pronged approach to improve our understanding of aerosol hygroscopicity and how it is implemented in Earth system models. In the first part of our project, we developed a benchmark dataset from existing aerosol hygroscopic growth measurements made by tandem nephelometer humidogram systems. We analyzed, using a standardized methodology, observations from 26 in-situ stations around the globe. Measurement data was collected from multiple data providers, reviewed and harmonized to create a consistent dataset of the scattering enhancement factor due to aerosol water uptake. This dataset is archived in several publicly available databases for use by interested researchers. In the second part of the project, we used the benchmark hygroscopicity dataset to perform a global study on aerosol hygroscopicity and aerosol optical properties. Measurements show a global picture of scattering enhancement with larger values for Arctic and marine sites and lower for urban and desert sites. We assessed the RH dependence of aerosol radiative forcing and showed that the overall effect of aerosol hygroscopicity on DARF is an increase in the absolute forcing effect (negative sign) by a factor of up to 4 compared to dry conditions (RH<40%). Finally, we explored using aerosol single scattering albedo (SSA) and scattering Angstrom exponent (SAE) as possible proxies for aerosol hygroscopicity. SSA showed more promise as a surrogate for the scattering enhancement factor than SAE, but neither was ideal. In the third part of the project we evaluated the output of ten Earth system models (ESMs) against the benchmark hygroscopicity dataset. ESMs utilize various schemes for aerosol hygroscopicity which had not been previously tested against observations on a global scale. We found that ESMs currently overestimate scattering enhancement due to hygroscopic growth. Model parameterizations of hygroscopicity and model chemistry are two main factors driving the observed diversity in hygroscopicity simulations among the models. In addition, our study makes several suggestions for modelers, including improving the parameterizations of organic and sea salt aerosol hygroscopicity. Future hygroscopicity model evaluation experiments should include the model data related to particles size which was not available for our study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

EvoNet: A phylogenomic and systems biology approach to identify genes underlying plant survival in marginal, low‐N soils

The DOE‐BER “EvoNet” project investigates the genetic and molecular basis of plant resilience in extreme environments. We do this by identifying key genes that enable “extreme survivor” species to thrive in the nitrogen-poor soils of Chile’s hyper-arid Atacama Desert. Our collections focus on 32 Atacama extremophile species, including seven grass species with potential biofuel applications. To identify genes-of-importance to survival we compared genomic and transcriptomic profiles of extremophile species that thrive in the Atacama to those of closely related “sister” species from nitrogen-rich arid and mesic regions of California. Deep RNA sequencing and de novo transcriptome assembly across these triplet species sets supported a phylogenomic framework for identifying positively selected genes associated with adaptive divergence. Our integrative analysis combined ecological and environmental data, metagenomics, evolutionary and systems biology, and metabolomics. This enabled us to create an unprecedented framework for systematically understanding how non-model plants have adapted to survive in extreme conditions. Our resulting database of positively selected ortholog groups in the extremophile plants offers promising targets for engineering crop and biofuel species with enhanced resilience to drought and extreme weather. Additionally, our newest dataset explores and exploits a complementary metabolomic approach. This new aspect provides innovative strategies to manipulate plant cell metabolism, further supporting efforts to improve agricultural productivity in the face of extreme climates. Importantly, our combined evolutionary- and metabolomic-based strategies focused on convergent patterns of adaptation, providing a genetic and metabolomic toolkit for improving crop and biofuel resilience across diverse plant species. Finally, our novel exploration of ecological and evolutionary dynamics delivered to the community a phylogenomic computational pipeline called “PhyloGeneious.” Our continued adaptations of this pipeline are publicly available to expedite evolutionary genomic research for future scientific discoveries. In total, our DOE-BER has provided genomic, metabolomic, and computational strategies to understand how extremophile plants provide evolutionary and physiological targets for improving agricultural and biofuel production.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

U.S. Department of Energy Free-Air CO 2 Enrichment Experiments: FACE Results, Lessons, and Legacy

In the 1970s, growing concern about potential long-range consequences of carbon dioxide (CO 2 ) emissions from the combustion of fossil fuels prompted the U.S. Department of Energy (DOE) to begin developing a CO 2 research program. The goal was to obtain the scientific information needed to reduce uncertainty and develop a knowledgebase for decision-making. At its start in 1980, this program had three components: the global carbon cycle, climate effects, and vegetation effects. The vegetation effects component recognized that increased atmospheric CO 2 could be a potential benefit to vegetation and crops because CO 2 is an essential resource for plant growth. In addition to increasing crop yield, program developers and researchers recognized that this fertilizer effect could also result in greater storage of fossil fuel CO 2 by forests. Initial experiments were small in scale and scope and focused on the photosynthesis, physiology, and water use of crop plants, but the focus shifted toward responses of plants in natural systems, as needed for connections to the global carbon cycle. Program research investigated CO 2 interactions with other environmental variables, emphasizing field research whenever possible and recognizing that exposure to elevated CO 2 concentrations under realistic field conditions should lead to better predictions of plant response. Whole-ecosystem field experiments were instituted in low-stature ecosystems, but experiments to address forest responses necessarily relied on tree seedlings or saplings in open-top chambers (OTCs) and not within intact forest ecosystems. The need for larger-scale, longer-duration experiments in each of the world’s six major biomes (tundra, boreal forest, temperate forest, tropical forest, grassland, and desert) became clear, and for these experiments, free-air CO 2 exposure facilities were necessary. This report summarizes the objectives and main results of DOE-sponsored experiments employing such facilities, along with the outcomes from research on ecosystem responses using field chambers (see Fig. ES.1, p. x). Additionally, this report highlights the many successes of DOE CO 2 enrichment studies and considers the lessons learned that can be used to guide future initiatives.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗
DOE OSTI2025

Remote-sensing detectability of airborne Arctic dust

Remote-sensing (RS)-based estimates of Arctic dust are oftentimes overestimated due to a failure in separating out the dust contribution from that of spatially homogeneous clouds or low-altitude cloud-like plumes. A variety of illustrations are given with a particular emphasis on questionable claims of using brightness temperature differences (BTDs) as a signature indicator of Arctic dust transported from mid-latitude deserts or generated by local Arctic sources. While there is little dispute about the presence of both Asian and local dust across the Arctic, the direct RS detectability of airborne dust, as ascribed to satellite (MODIS and AVHRR) measurements of significantly negative brightness temperature differences at 11 and 12 µm (BTD 11–12 ), has been misrepresented in certain cases. While it is difficult to account for all examples of strongly negative BTD 11–12 values in the Arctic, it is unlikely that airborne dust plays a significant role. One much more likely contributor would be water clouds in the Arctic inversion layer. The RS detectability of the impact of Arctic dust (notably due to Arctic dust from local sources) can, however, be of significance. Sustained dust deposition can substantially decrease (visible to shortwave IR) snow and ice reflectance albedo (pan-chromatic reflectance) and the signal measured by satellite sensors. Significantly negative BTD 11–12 values would, however, only represent a limited area near the drainage basin sources according to our event-level case studies. The enhanced ice-nucleating particle (INP) role of local Arctic dust can, for example, induce significant changes in the properties of low-level mixed-phase clouds (cloud optical depth changes <~ 1) that can readily be detected by active and passive RS instruments. It is critical that the distinction between the RS detectability of airborne Arctic dust versus the RS detectability of the impacts of that dust be understood if we are to appropriately parameterize, for example, the radiative forcing influence of dust in this climate-sensitive region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Assessing the climate impact of an improved volcanic sulfate aerosol representation in E3SM

Accurately simulating historical surface temperature variations is essential for evaluating climate models, yet many struggle to reproduce the mid-20th-century temperature trends associated with significant volcanic eruptions. This study examines the impact of volcanic sulfate aerosol representation on these biases using the Energy Exascale Earth System Model (E3SM). The standard CMIP6 protocol prescribes volcanic forcing through radiative perturbations, omitting volcanic aerosol–cloud interactions (VACIs). Here, we implement an emission-based approach with an updated volcanic eruption inventory that directly incorporates volcanic sulfur dioxide (SO 2 ) emissions, enabling a more process-based representation of volcanic forcing. This approach leads to improved surface temperature variability and a modest reduction in cold biases between 1940 and 1980 compared to the CMIP6 setup. Additionally, we assess cloud property responses to a more realistic volcanic sulfate aerosol representation, which weakens cloud-induced cooling during periods of lower volcanic activity. However, despite these refinements, a significant temperature cold bias remains, indicating that further improvements in atmospheric chemistry, aerosol microphysics, cloud processes, and model parameterizations are needed to fully resolve this issue in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Multi-ignition fire complexes drive extreme fire years and impacts

Climate change is intensifying fire behavior, with the largest and fastest-spreading fires causing the greatest impacts on people and ecosystems. Yet the mechanisms driving variability and trends in large fires remain poorly understood. Using 12-hour satellite-derived fire tracking data from 2012 to 2023, we show that the merging of separate ignitions into multi-ignition complexes is a key process amplifying fire size and destructive potential across temperate and boreal ecoregions. Multi-ignition fires account for 31% of the burned area in California and 59% in the Arctic-boreal domain, spread faster and persist longer than single-ignition fires, and disproportionately contribute to extreme fire years in California, Canada, and Siberia. They also generate stronger atmospheric feedbacks, produce more pyrocumulonimbus events, and strain firefighting capacity by dispersing suppression resources. Recognizing and accounting for fire-merging dynamics are critical for improving wildfire prediction, risk assessment, and management.

Environmental sciences↗
DOE OSTI2021

Climatic Controls on Soil Carbon Accumulation and Loss in a Dryland Ecosystems

Arid and semiarid ecosystems drive year-to-year variability in the strength of the terrestrial carbon (C) sink, yet there is uncertainty about how soil C gains and losses contribute to this variation. For this work, to address this knowledge gap, we embedded C-depleted soil mesocosms, containing litter or biocrust C inputs, within an in situ dryland ecosystem warming experiment. Over the course of one year, changes in microbial biomass and total soil organic C pools were monitored alongside hourly measurements of soil CO 2 flux. We also developed a biogeochemical model to explore the mechanisms that gave rise to observed soil C dynamics. Field data and model simulations demonstrated that water exerted much stronger control on soil biogeochemistry than temperature, with precipitation events triggering large CO 2 pulses and transport of litter- and biocrust-derived C into the soil profile. We expected leaching of organic matter would result in steady accumulation of C within the mineral soil over time. Instead, the size of the total organic C pool fluctuated throughout the year, largely in response to microbial growth: increases in the size of microbial biomass were negatively correlated with the quantity of C residing in the top 2 cm, where most biogeochemical changes were observed. Our data and models suggest that microbial responses to precipitation events trigger rapid metabolism of dissolved organic C inputs, which strongly limit accumulation of autotroph-derived C belowground. Accordingly, changes in the magnitude and/or frequency of precipitation events in this dryland ecosystem could have profound impacts on the strength of the soil C sink.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Energy Exascale Earth System Model Version 3: 1. Overview of the Atmospheric Component

This paper describes the atmospheric component of the US Department of Energy's Energy Exascale Earth System Model (E3SM) version 3. Significant updates have been made to the atmospheric physics compared to earlier versions. Specifically, interactive gas chemistry has been implemented, along with improved representations of aerosols and dust emissions. A new stratiform cloud microphysics scheme more physically treats ice processes and aerosol‐cloud interactions. The deep convection parameterization has been largely improved with sophisticated microphysics for convective clouds, making model convection sensitive to large‐scale dynamics, and incorporating the dynamical and physical effects of organized mesoscale convection. Improvements in aerosol wet removal processes and parameter re‐tuning of key aerosol and cloud processes have improved model aerosol radiative forcing. The model's vertical resolution has increased from 72 to 80 layers with the extra eight layers added in the lower stratosphere to better simulate the Quasi‐Biennial Oscillation. These improvements have enhanced E3SM's capability to couple aerosol, chemistry, and biogeochemistry and reduced some long‐standing biases in simulating tropical variability. Compared to its predecessors, the model shows a much stronger signal for the Madden‐Julian Oscillation, Kelvin waves, mixed Rossby‐gravity waves, and eastward inertia‐gravity waves. Aerosol radiative forcing has been considerably reduced and is now better aligned with community best estimates, leading to significantly improved skill in simulating historical temperature records. Its simulated mean‐state climate is largely comparable to E3SMv2, but with some notable degradation in shortwave cloud radiative effect, precipitable water, and surface wind stress, which will be addressed in future updates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Strategies to mitigate urban heat: Effects on overheating and cloud formation

This study evaluates the effectiveness of various urban heat island (UHI) overheating mitigation strategies and their associated impacts on urban cloud dynamics and thermal processes. This study shows cloud-resolving and urban-resolving modeling results estimating the impact of Houston-Galveston heat mitigation scenarios and other resilient strategies contemplated in the Climate Adaptation Plan and Resilient Houston reports. The simulated scenarios include high intensity green rooftops, rooftop photovoltaic solar panels, enhanced urban irrigation, white/cool roofs and roads, and street trees. We contrast the adaptation scenarios against a present baseline case, a no city scenario and a larger and denser city as projected by the Houston-Galveston Area Council for 2045. During the daytime, cooling strategies such as cool roofs, cool roads, and green roofs exhibit superior performance in mitigating urban overheating. At night, enhanced urban irrigation emerges as the most effective cooling intervention. Cooling strategies significantly reduce sensible heat flux partitioning during the day, a process that reduces the uplift of air, suppressing the formation of urban shallow cumulus clouds. The extent of urban cloud mixing ratio is reduced in proportion to the decrease in sensible heating. Under the BEP-Tree scenario, which includes wind effects and evapotranspiration driven by a stomatal conductance model, urban trees demonstrated negligible environmental cooling effects and minimal urban cloud modifications. In contrast, the scenarios with more urban cooling and higher latent heat fluxes led to suppressed urban clouds. The net cooling effect achieved by the heat mitigation strategies is influenced by a combination of indirect processes, including the reduction of downwelling longwave radiation flux, due to reduced cloud presence, while some warming is attributed to a modest increase in shortwave radiation that offsets the cooling benefit. Additionally, reduced heat dissipation, weakened thermal gradients, and diminished vertical mixing over urban areas further moderate the cooling potential. These findings highlight the pivotal role of clouds and moist atmospheric processes in shaping the UHI effect and offer insights for designing more effective urban cooling strategies.

albedo changes↗
DOE OSTI2025

Mesoscale Convective Systems Tracking Method Intercomparison (MCSMIP): Application to DYAMOND Global km‐Scale Simulations

Abstract Global kilometer‐scale models represent the future of Earth system modeling, enabling explicit simulation of organized convective storms and their associated extreme weather. Here, we comprehensively evaluate tropical mesoscale convective system (MCS) characteristics in the DYAMOND (DYnamics of the atmospheric general circulation modeled on non‐hydrostatic domains) simulations for both summer and winter phases. Using 10 different feature trackers applied to simulations and satellite observations, we assess MCS frequency, precipitation, and other key characteristics. Substantial differences (a factor of 2–3) arise among trackers in observed MCS frequency and their precipitation contribution, but model‐observation differences in MCS statistics are more consistent across trackers. DYAMOND models are generally skillful in simulating tropical mean MCS frequency, with multi‐model mean biases ranging from −2%–8% over land and −8%–8% over ocean (summer vs. winter). However, most DYAMOND models underestimate MCS precipitation amount (23%) and their contribution to total precipitation (17%). Biases in precipitation contributions are generally smaller over land (13%) than over ocean (21%), with moderate inter‐model variability. While models better simulate MCS diurnal cycles and cloud shield characteristics, they overestimate MCS precipitation intensity and underestimate stratiform rain contributions (up to a factor of 2), particularly over land, albeit observational uncertainties exist. Additionally, models exhibit a wide range of precipitable water in the tropics compared to reanalysis and satellite observations, with many models showing exaggerated sensitivity of MCS precipitation intensity to precipitable water. The MCS metrics developed here provide process‐oriented diagnostics to guide future model development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Survival and behavior of Mojave desert tortoises head-started with and without outdoor rearing

Mojave desert tortoise (Gopherus agassizii) populations in some regions have declined by >50% since 2004, prompting the need for more research on ways to recover populations. One possible recovery tool is head-starting (i.e., the act of protecting and raising juvenile tortoises to sizes that increase survival upon release); however, head-starting can have high start-up and maintenance costs that can limit its feasibility. Strategies that reduce cost and rearing duration may foster broader and more effective use. We released and radio-tracked 60 juvenile tortoises in the Mojave National Preserve in California, USA, that had been reared under 2 treatments: those reared 1 year indoors after hatching, then 1 year outdoors (combo) and those reared just 1 year indoors (indoor-only). We tested whether indoor-only rearing alone could be a more efficient means of producing robust head-started tortoises. We examined the behavior, movement, and survival of tortoises after release into the wild from 2020 to 2021 to determine whether these outcomes differed between husbandry treatments. Combo tortoises tended to perform settling behaviors (mean ± 1 SE days to building first burrow = 6.7 ± 0.8, entering dormancy = 23.3 ± 2.1, and emerging from dormancy = 189.6 ± 4.4) earlier than indoor-only tortoises (7.4 ± 0.9, 31.5 ± 2.6, and 193.9 ± 5.9, respectively), but this difference was not significant, suggesting the rearing method did not greatly alter settling behavior. Indoor-only tortoises dispersed at least twice as far from their release site (156.2 ± 26.3 m compared with 77.3 ± 20.6 m for combo tortoises), had larger mean use areas (3.7 ± 0.1 ha compared with 2.8 ± 0.3 ha for combo tortoises for 95% Brownian bridge movement model estimates), and greater variability in their movements than combo tortoises (daily average step length post-emergence: 4.3 ± 0.2 m compared with 2.8 ± 0.1 m for combo tortoises). Despite differences in their movements, indoor-only and combo tortoises had similar survival rates over the study, 51% versus 42%, respectively, during a period of extreme drought in 2021. The similarity in survival between groups gives head-starting practitioners freedom in their rearing methods. The indoor-only group had lower site fidelity, which should be considered when this is an undesirable trait for released tortoises.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"

This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Location Identifiers, Metadata, and Map for Field Measurements at the East-Taylor Watershed Community Observatory, Colorado, USA (Version 3.3)

This dataset contains identifiers, metadata, and a map of the locations where field measurements have been conducted at the East-Taylor Watershed Community Observatory located in the Upper Colorado River Basin, United States. This is version 3.3 of the dataset and replaces the prior version 3.2 (see below for details on changes between the versions). Dataset description: The East River-Taylor Watershed is the primary field site of the Watershed Function Scientific Focus Area (WFSFA) and the Rocky Mountain Biological Laboratory. Researchers from several institutions generate highly diverse hydrological, biogeochemical, climate, vegetation, geological, remote sensing, and model data at the East-Taylor Watershed in collaboration with the WFSFA. Thus, the purpose of this dataset is to maintain an inventory of the field locations and instrumentation to provide information on the field activities in the East-Taylor Watershed and coordinate data collected across different locations, researchers, and institutions. The dataset contains (1) a README file with information on the various files, (2) three csv files describing the metadata collected for each surface point location, plot and region registered with the WFSFA, (3) csv files with metadata and contact information for each surface point location registered with the WFSFA, (4) a csv file with with metadata and contact information for plots, (5) a csv file with metadata for geographic regions and sub-regions within the watershed, (6) a compiled xlsx file with all the data and metadata which can be opened in Microsoft Excel, (7) a kml map of the locations plotted in the watershed which can be opened in Google Earth, (8) a jpg image of the kml map which can be viewed in any photo viewer, and (9) a zipped file with the registration templates used by the SFA team to collect location metadata. The zipped template file contains two csv files with the blank templates (point and plot), two csv files with instructions for filling out the location templates, and one compiled xlsx file with the instructions and blank templates together. Additionally, the templates in the xlsx include drop down validation for any controlled metadata fields. Persistent location identifiers (Location_ID) are determined by the WFSFA data management team and are used to track data and samples across locations. Dataset uses: This location metadata is used to update the Watershed SFA’s publicly accessible Field Information Portal (an interactive field sampling metadata exploration tool; https://wfsfa-data.lbl.gov/watershed/), the kml map file included in this dataset, and other data management tools internal to the Watershed SFA team. Version Information: The latest version of this dataset publication is version 3.3. This version contains 167 new point locations, 1 new plot, and 2 new geographic regions. Overall, there are a total of 1439 point locations, 75 plots, and 54 geographic regions. Additionally, the kml map of locations and image now includes two boundaries (Upper Ohio Creek (UO) and Carbon Creek (CA)) outside of the East River watershed (USGS HUC-10) and accompanying stream network that represents areas of focus. Refer to methods for further details on the version history. This dataset will be updated on a periodic basis with new measurement location information. Researchers interested in having their East-Taylor Watershed measurement locations added to this list should reach out to the WFSFA data management team at wfsfa-data@googlegroups.com. Acknowledgments: Please cite this dataset if using any of the location metadata in other publications or derived products. If using the location metadata for the 2018 NEON hyperspectral campaign, additionally cite Chadwick et al. (2020). doi:10.15485/1618130. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. Part of this work was performed at SLAC Accelerator Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-76SF00515.

2018 NEON and 2025 CHESS Campaigns↗