Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “snowfall”

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

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Turning Satellite Data into Global Precipitation Maps

It's really obvious how to measure rain and snowfall just put a container out in the yard. But, it turns out the physical processes that create precipitation, and the practicalities of taking measurements everywhere, make satellites key tools in mapping rain and snowfall around the globe. Employing satellites in this meteorological detective story involves some science, some technology, a lot of computing, and, yes, rocket science. The payoff is that we know more about how much rain and snow fell in, say, 2010, than we knew at the time that it fell, and we're ready to make the same kind of estimate for tomorrow's storms. There is still a lot of work to be done to make these estimates more reliable and accurate.Dr. George J. Huffman is a Research Meteorologist at NASA Goddard Space Flight Center (GSFC), Greenbelt, Maryland. After completing a B.S. in Physics at The Ohio State University (1976) and a Ph.D. in Meteorology at Massachusetts Institute of Technology (1982), he was an Assistant Professor at University of Maryland, College Park before moving to GSFC in 1988, where he consulted for 24 years before joining NASA in 2012. Dr. Huffman focuses on the design, implementation, and extension of combined estimates of global precipitation, seeking to incorporate all available precipitation-relevant satellite and global surface data. He is the Deputy Project Scientist for the Global Precipitation Measurement (GPM) mission, the Chair of the GPM Multi-satellite Algorithm Working Group, and Chief of the Mesoscale Atmospheric Processes Lab.

Huffman, George J.↗
NASA NTRS2020

Unprecedented Atmospheric Conditions (1948–2019) Drive the 2019 Exceptional Melting Season Over the Greenland Ice Sheet

Understanding the role of atmospheric circulation anomalies on the surface mass balance of the Greenland ice sheet (GrIS) is fundamental for improving estimates of its current and future contributions to sea level rise. Here, we show, using a combination of remote sensing observations, regional climate model outputs, reanalysis data, and artificial neural networks, that unprecedented atmospheric conditions (1948–2019) occurring in the summer of 2019 over Greenland promoted new record or close-to-record values of surface mass balance (SMB), runoff, and snowfall. Specifically, runoff in 2019 ranked second within the 1948–2019 period (after 2012) and first in terms of surface mass balance negative anomaly for the hydrological year 1 September 2018–31 August 2019. The summer of 2019 was characterized by an exceptional persistence of anticyclonic conditions that, in conjunction with low albedo associated with reduced snowfall in summer, enhanced the melt–albedo feedback by promoting the absorption of solar radiation and favored advection of warm, moist air along the western portion of the ice sheet towards the north, where the surface melt has been the highest since 1948. The analysis of the frequency of daily 500 hPa geopotential heights obtained from artificial neural networks shows that the total number of days with the five most frequent atmospheric patterns that characterized the summer of 2019 was 5 standard deviations above the 1981–2010 mean, confirming the exceptional nature of the 2019 season over Greenland.

Greenland ice sheet↗
NASA NTRS2020

Arctic Snow Depth and Sea Ice Thickness From ICESat‐2 and CryoSat‐2 Freeboards: A First Examination

We present a first examination of Arctic sea ice snow depth estimates from differencing satellite lidar (ICESat‐2) and radar (CryoSat‐2) freeboards. These estimates cover the period between 14 October 2018 and the end of April 2019. Snow depth is related to freeboard differences by the refractive index/bulk density of the snow layer—the only free parameter in the approach. Area‐averaged snow depth ranges from 9 cm (on first‐year ice: 5 cm, multiyear ice: 14 cm) in late October to 19 cm (first‐year ice: 17 cm, multiyear ice: 27 cm) in April; on average, this snow is thinner over FYI. Spatial patterns and gradients of snow depth estimates compare well with reconstructions using snowfall from ERA‐Interim and ERA5, although snowfall from ERA5 is systematically higher. For all months, the results suggest that ~50% of the total freeboard is comprised of snow. Retrievals are within a few centimeters of snow depth data acquired by Operation IceBridge in April 2019. Sources of uncertainties associated with this freeboard‐differencing approach are discussed. Further, sea ice thicknesses calculated using the retrieved snow depth and a modified climatology are contrasted. Comparatively, the snow depth and calculated ice thickness using a modified climatology are higher by ~5 cm and 0.33 m, although these differences are not uniform throughout the season. Snow accumulation was slower between October and December but increased between December and January, unlike the modified climatology, which exhibited a monotonic accumulation for all months. Future opportunities for assessment and improvement of these estimates are discussed.

Arctic Ocean↗
NASA NTRS2020

The Precipitation Imaging Package: Assessment of Microphysical and Bulk Characteristics of Snow

Remote-sensing observations are needed to estimate the regional and global impacts of snow. However, to retrieve accurate estimates of snow mass and rate, these observations require augmentation through additional information and assumptions about hydrometeor properties. The Precipitation Imaging Package (PIP) provides information about precipitation characteristics and can be utilized to improve estimates of snowfall rate and accumulation. Here, the goal is to demonstrate the quality and utility of two higher-order PIP-derived products: liquid water equivalent snow rate and an approximation of volume-weighted density called equivalent density. Accuracy of the PIP snow rate and equivalent density is obtained through intercomparison with established retrieval methods and through evaluation with colocated ground-based observations. The results confirm the ability of the PIP-derived products to quantify properties of snow rate and equivalent density, and demonstrate that the PIP produces physically realistic snow characteristics. When compared to the National Weather Service (NWS) snow field measurements of six-hourly accumulation, the PIP-derived accumulations were biased only +2.48% higher. Additionally, this work illustrates fundamentally different microphysical and bulk features of low and high snow-to-liquid ratio events, through assessment of observed particle size distributions, retrieved mass coefficients, and bulk properties. Importantly, this research establishes the role that PIP observations and higher-order products can serve for constraining microphysical assumptions in ground-based and spaceborne remotely sensed snowfall retrievals.

Claire Pettersen↗
NASA NTRS2022

Tracking the Impacts of Precipitation Phase Changes Through the Hydrologic Cycle in Snowy Regions: From Precipitation to Reservoir Storage

Cool season precipitation plays a critical role in regional water resource management in the western United States. Throughout the twenty-first century, regional precipitation will be impacted by rising temperatures and changing circulation patterns. Changes to precipitation magnitude remain challenging to project; however, precipitation phase is largely dependent on temperature, and temperature predictions from global climate models are generally in agreement. To understand the implications of this dependence, we investigate projected patterns in changing precipitation phase for mountain areas of the western United States over the twenty-first century and how shifts from snow to rain may impact runoff. We downscale two bias-corrected global climate models for historical and end-century decades with the Weather Research and Forecasting (WRF) regional climate model to estimate precipitation phase and spatial patterns at high spatial resolution (9 km). For future decades, we use the RCP 8.5 scenario, which may be considered a very high baseline emissions scenario to quantify snow season differences over major mountain chains in the western U.S. Under this scenario, the average annual snowfall fraction over the Sierra Nevada decreases by >45% by the end of the century. In contrast, for the colder Rocky Mountains, the snowfall fraction decreases by 29%. Streamflow peaks in basins draining the Sierra Nevada are projected to arrive nearly a month earlier by the end of the century. By coupling WRF with a water resources model, we estimate that California reservoirs will shift towards earlier maximum storage by 1–2 months, suggesting that water management strategies will need to adapt to changes in streamflow magnitude and timing.

Melissa L. Wrzesien↗
NASA NTRSDate not supplied

A LaGrangian Perspective of the Surface Energy Budget Contributions to Sea Ice Evolution in the ‘New Arctic’

Recent Arctic warming outpaces global warming and is characterized by dramatic changes in sea ice (record sea ice loss, thinner and more seasonal ice cover), warming SSTs, and a lengthening melt season. These changes impact the surface energy budget and the sensitivity of sea ice to radiative forcing, namely an increased sensitivity of thinner ice to local radiative forcings and a weaker buffering of surface-atmosphere exchanges of heat and moisture. However, these processes are a large source of uncertainty in global climate models with prior studies finding that poor representations of downwelling longwave radiation and sea ice albedo lead to model spread in sea ice conditions. A limitation of prior work is examining sea ice-atmosphere interactions through an Eulerian framework. This study utilizes a sea ice parcel dataset in a LaGrangian framework, allowing for tracking of surface energy budget-sea ice interactions over the lifespan of individual sea ice parcels; this framework also allows for the quantification of the sea ice response to episodic weather events along with seasonal drivers, including the drivers of rapid growth and melt events. First, lag-lead relationships between relevant surface energy budget terms and sea ice thickness changes are explored to understand the atmospheric conditions preceding rapid growth and melt. The impact of episodic events- moisture intrusion, warm air outbreaks, cyclones- on sea ice growth are also investigated by stratifying first-year parcels by sea ice thickness at the end of the growth season and analyzing the effect of atmospheric variability (number and intensity of episodic events) on parcel thickness change. Lastly, parcel growth and melt and the associated radiative anomalies are compared to snowfall events/snow thickness to determine how snowfall modifies the interaction between the local surface energy budget and sea ice mass balance.

Patrick Taylor↗
NASA NTRS2021

Chapter 6: Satellite Rainfall Estimation

The measurement of precipitation (rainfall and snowfall) across the globe is important for a variety of scientific and social applications. Knowledge of the occurrence and amount of precipitation, together with its distribution and changes are crucial for improving our understanding of the global energy and water cycle, as well as for monitoring water resource availability and for hydrological modelling to help alleviate flood and drought impacts. The use of conventional instruments (gauge or radar) to map global precipitation is essentially limited to land areas and thus satellite observations must be used to provide estimates of global precipitation. Many satellite sensors operating over the last 50 years have provided data for a range of techniques, algorithms and schemes developed to generate quantitative precipitation estimates. Current satellite-based precipitation products can provide estimates at up to 4 km every 30 minutes. This chapter outlines the basis of satellite precipitation estimation, the satellites and sensors used, and the range of techniques and schemes used to generate the precipitation products.

Earth observation↗
NASA NTRS2021

Chapter 10: Development of a Precipitation Retrieval Scheme for Cross-track Passive Microwave Sounding Instruments

The measurement of precipitation (snowfall and rainfall) across the globe is of great importance to both science and society. Despite the importance, measurements from surface-based instruments are limited in extent and are generally absent from sparsely populated or inhospitable regions, including over the world’s oceans. Satellite-based observations have the potential to provide the necessary information to generate estimates of precipitation across the whole of the Earth’s surface. However, adequate temporal and spatial sampling is required to capture vagaries of precipitation and it is necessary to exploit all possible satellite (and, where appropriate surface) information. While much research has been carried out into the retrieval of precipitation from observations gathered by conically-scanning passive microwave imagers in previous decades, relatively little has been done to fully utilise the observations from the cross-track scanning passive microwave sounding instruments. This paper will review a number of cross-track precipitation retrieval schemes, together an evaluation of their performance, and future directions.

Precipitation retrievals↗
NASA NTRS2021

The Climates of Earth's Next Supercontinent: Effects of Tectonics, Rotation Rate, and Insolation

We explore two possible Earth climate scenarios, 200 and 250 million years into the future, using projections of the evolution of plate tectonics, solar luminosity, and rotation rate. In one scenario, a supercontinent forms at low latitudes, whereas in the other it forms at high northern latitudes with an Antarctic subcontinent remaining at the south pole. The climates between these two end points are quite stark, with differences in mean surface temperatures approaching several degrees. The main factor in these differences is related to the topographic height of the high latitude supercontinents where higher elevations promote snowfall and subsequent higher planetary albedos. These results demonstrate the need to consider multiple boundary conditions when simulating Earth-like exoplanetary climates.

insolation↗
NASA NTRS2023

Sensitivity of the MAR Regional Climate Model Snowpack to the Parameterization of the Assimilation of Satellite-Derived Wet-Snow Masks on the Antarctic Peninsula

Both regional climate models (RCMs) and remote sensing (RS) data are essential tools in understanding the response of polar regions to climate change. RCMs can simulate how certain climate variables, such as surface melt, runoff and snowfall, are likely to change in response to different climate scenarios but are subject to biases and errors. RS data can assist in reducing and quantifying model uncertainties by providing indirect observations of the modeled variables on the present climate. In this work, we improve on an existing scheme to assimilate RS wet snow occurrence data with the “Modèle Atmosphérique Régional” (MAR) RCM and investigate the sensitivity of the RCM to the parameters of the scheme. The assimilation is performed by nudging the MAR snowpack temperature to match the presence of liquid water observed by satellites. The sensitivity of the assimilation method is tested by modifying parameters such as the depth to which the MAR snowpack is warmed or cooled, the quantity of water required to qualify a MAR pixel as “wet” (0.1 % or 0.2 % of the snowpack mass being water), and assimilating different RS datasets. Data assimilation is carried out on the Antarctic Peninsula for the 2019–2021 period. The results show an increase in meltwater production (+66.7 % on average, or +95 Gt), along with a small decrease in surface mass balance (SMB) (−4.5 % on average, or −20 Gt) for the 2019–2020 melt season after assimilation. The model is sensitive to the tested parameters, albeit with varying orders of magnitude. The prescribed warming depth has a larger impact on the resulting surface melt production than the liquid water content (LWC) threshold due to strong refreezing occurring within the top layers of the snowpack. The values tested for the LWC threshold are lower than the LWC for typical melt days (approximately 1.2 %) and impact results mainly at the beginning and end of the melting period. The assimilation method will allow for the estimation of uncertainty in MAR meltwater production and will enable the identification of potential issues in modeling near-surface snowpack processes, paving the way for more accurate simulations of snow processes in model projections.

regional climate models↗
DOE OSTI2025

Long-Term Photovoltaic System Performance in Cold, Snowy Climates

As countries around the world transition towards renewable energy, there is increasing interest in using photovoltaic (PV) technologies to help decarbonize northern and alpine communities due to their scalability and affordability. However, a barrier to large-scale adoption of PV in cold climates is long-term performance uncertainty under snowfall, freeze-thaw cycles, low temperatures, and high winds. In this work, we provide a comprehensive review of published silicon degradation rates in cold Koppen-Geiger climate classifications of Dfb (humid continental), Dfc (subarctic), and ET (tundra). We first analyze the system degradation rates of three subarctic ground-mounted photovoltaic sites in North America using the RdTools year-on-year method: an Al-BSF double-axis tracking site in Fairbanks, Alaska (65degrees N); a PERC and silicon heterojunction bifacial vertical and south-tilted site in Fairbanks, Alaska; and a PERC south-facing fixed-tilt site in Fort Simpson, Northwest Territories (62degrees N). Degradation rates of these newly analyzed sites vary between -0.4%/year and -1.5%/year. Combining these data with previously reported cold climate degradation rates, we show that the distribution of cold climate degradation peaks at -0.1%/year to -0.2%/year but has a large tail with rates above -0.5%/year. The average reported cold climate degradation rate is -0.45%/year, whereas the median value is -0.33%/year. These results suggest that despite frequent freeze-thaw cycles and potential exposure to high wind and snow loads, PV systems in cold climates tend to degrade slower than PV systems in warmer climates. The limited sample size of reported degradation rates in cold climates (27) motivates the need for further data acquisition and monitoring efforts as new technologies are deployed.

14 SOLAR ENERGY↗
DOE OSTI2024

Anticipating how rain-on-snow events will change through the 21st century: lessons from the 1997 new year’s flood event

The California-Nevada 1997 New Year’s flood was an atmospheric river (AR)-driven rain-on-snow (RoS) event and remains the costliest in their history. The joint occurrence of saturated soils, rainfall, and snowmelt generated inundation throughout northern California-Nevada. Although AR RoS events are projected to occur more frequently with climate change, the warming sensitivity of their flood drivers across scales remains understudied. We leverage the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) to recreate the 1997 New Year’s flood with horizontal grid spacings of 3.5 km across California, with forecast lead times of up to 4 days, and across six warming levels ranging from pre-industrial conditions to +3.5° C. We describe the sensitivity of the flood drivers to warming including AR duration and intensity, precipitation phase, intensity and efficiency, snowpack mass and energy changes, and runoff efficiency. Our findings indicate current levels of climate change negligibly influence the flood drivers. At warming levels ≥ 1.7° C, AR hazard potential increases, snowpack nonlinearly decreases, antecedent soil moisture decreases (except where the snowline retreats), and runoff decreases (except in the southern Sierra Nevada where antecedent snowpack persists). Storm total precipitation increases, but at rates below warming-induced increases in saturation-specific humidity. Warming intensifies short-duration, high-intensity rainfall, particularly where snowfall-to-rainfall transitions occur. This study highlights the nonlinear tradeoffs in 21st-century RoS flood hazards with warming and provides water management and infrastructure investment adaptation considerations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impacts of Topography-Based Subgrid Scheme and Downscaling of Atmospheric Forcing on Modeling Land Surface Processes in the Conterminous US

The effects of small-scale topography-induced land surface heterogeneity are not well represented in current Earth System Models (ESMs). In this study, a new topography-based subgrid structure referred to as topographic units (TGU) designed to better capture subgrid topographic effects, and methods to downscale atmospheric forcing to the land TGUs have been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Effects of the subgrid scheme and downscaling methods on ELM simulated land surface processes are evaluated over the conterminous United States (CONUS). For this purpose, ELM simulations are performed using two configurations without (NoD ELM) and with (D ELM) downscaling, both using TGUs derived for the 0.5-degree grids and the same land surface parameters. Simulations using the two ELM configurations are compared over the CONUS domain, regional levels, and at observational sites (e.g., SNOTEL). The CONUS-level results suggest that D ELM simulates more snowfall and snow water equivalent (SWE), higher runoff, and less ET during spring and summer. Regional-level results suggest more pronounced impacts of downscaling over regions dominated by higher elevation TGUs and regions with maximum precipitation occurring during cool seasons. Results at the SNOTEL sites suggest that D ELM has superior capability of reproducing the observed SWE at 83% of the sites, with more pronounced performance over topographically heterogeneous TGUs with their maximum precipitation occurring during cool seasons. The results highlight the importance of improving representation of small-scale surface heterogeneity in ESMs and motivate future research to understand their effects on land-atmosphere interactions, streamflow, and water resources management over mountainous regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Unprecedented Shifts in Hydrology Are Emerging Across California's Critical Basins: An Evaluation From 0.5 to 3.5°C

With advances in climate models and downscaling techniques, stakeholders anticipate high-resolution analysis to inform regional to local changes in water management. Here, we produce hydrologic projections from an ensemble of Earth System Models (ESMs) that were selected and downscaled to support California's 5th Climate Assessment. An ensemble of 19 ESMs was downscaled to a 3-km resolution across California using a statistical-dynamical downscaling approach and subsequently run through two calibrated hydrology models. Although California has been extensively studied in the context of climate change, we provide the first evaluation of the warming thresholds at which hydroclimate metrics demonstrate statistically significant shifts. We show that present-day to near-term warming levels in Klamath and Northern Sierra Nevada basins, which serve as a critical source of water for California, show statistically significant decreases in snowfall and peak snowpack and associated decreases in summer snowmelt and runoff. More generally, shifts in these hydroclimate metrics occur for intermediate elevation basins ranging from 1,315 to 1,455 m (4,314.3–4,773.6 ft), while the warming level of emergence is delayed for lower and higher elevation basins. We also find that several basins already demonstrate significant increases in 5- to 100-year runoff intensities, primarily due to the increasing influence of precipitation falling as rain. Hydroclimate metrics with trends that demonstrate near-term warming levels of emergence are reflected in reanalysis-based observations, suggesting California is entering a fundamentally different hydroclimate regime. While this will likely stress California's water management, the research provided can support when to implement adaptation efforts.

Bass, B. [Univ. of California, Los Angeles, CA (Un↗
DOE OSTI2024

Climate change-resilient snowpack estimation in the Western United States

Abstract In the 21st century, warmer temperatures and changing atmospheric circulation will likely produce unprecedented changes in Western United States snowfall 1–3 , with impacts on the timing, amount, and spatial patterns of snowpack 4–7 . The ~900 snow pillow stations are indispensable to water resource management by measuring snow-water equivalent (SWE) 8,9 in strategic but fixed locations 10,11 . However, this network may not be impacted by climate change in the same way as the surrounding area 12 and thus fail to accurately represent unmeasured locations; climate change thereby threatens our ability to measure the effects of climate change on snow. In this work, we show that maintaining the current peak SWE estimation skill is nonetheless possible. We find that explicitly including spatial correlations—either from gridded observations or learned by the model—improves skill at predicting distributed snowpack from sparse observations by 184%. Existing artificial intelligence methods can be useful tools to harness the many available sources of snowpack information to estimate snowpack in a nonstationary climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

X-Band Radar and Surface-Based Observations of Cold-Season Precipitation in Western Colorado’s Complex Terrain

Abstract Hydrologic processes associated with intermountain cold-season precipitation in the Upper Colorado River basin have important impacts on avalanche forecasting and water resource management. However, traditional weather radar networks struggle with observations in this complex terrain. Data collected during the Study of Precipitation, the Lower Atmosphere, and the Surface for Hydrometeorology (SPLASH) and its sister campaign, Surface Atmosphere Integrated Field Laboratory (SAIL) in the East River watershed of western Colorado, are used to examine a multistorm period from 23 December 2021 to 1 January 2022 that contributed 35% of the total winter precipitation in this watershed. Dual-polarization X-band radar and disdrometer measurements show ∼30-mm differences in precipitation amount at two sites in proximity over four distinct storm events within the period. Wind patterns, synoptic forcings, microphysical characteristics of precipitation, and surface meteorology are analyzed to explain the observed spatial variability of cold-season precipitation in complex mountainous terrain. Analysis shows that differences over time within this event are mainly accounted for by synoptic forcings, such as frontal passages; differences between sites are accounted for by the impact of variations in local wind patterns on precipitation microphysics. Patterns of surface precipitation intensity are compared and found to be correlated with X-band radar signatures; a relationship between a strong dendritic growth stage and intense low-density surface precipitation is reinforced by this study. This relationship demonstrates the importance of particle growth mechanisms on surface snowfall patterns in high-altitude complex terrain, underscoring the importance of realistic microphysical parameterizations. Significance Statement The amount and density of snowpack from western Colorado winter storms have significant impacts on water resources in the Upper Colorado River basin. Snowpack characteristics are affected by small-scale differences in how snow forms in the atmosphere. These differences are hard to study in the complex terrain of the Rockies, but data from the SPLASH and SAIL field campaigns allows us to investigate how snow crystal formation and mountain-driven wind patterns affect snow near the surface. Our study finds that snow crystal growth varies over small space and time scales and is likely controlled by the terrain beneath a given location and resultant local wind patterns. These results imply that predicting snowpack in the Rockies requires properly representing local wind patterns and crystal growth processes in models.

Heflin, Stella↗
DOE OSTI2025

Characteristics and effects of aerosols during blowing snow events in the central Arctic

Sea salt aerosol (SSaer) significantly impacts aerosol-radiation and aerosol-cloud interactions, and sublimated blowing snow is hypothesized to be an important SSaer source in polar regions. Understanding blowing snow and other wind-sourced aerosols’ climate relevant properties is needed, especially during winter when Arctic amplification is greatest. However, most of our understanding of blowing snow SSaer comes from modeling studies, and direct observations are sparse. Additionally, SSaer can originate from multiple sources, making it difficult to disentangle emission processes. Here, we present comprehensive observations of wind-sourced aerosol during blowing snow events from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic. High wind speed strongly enhances total aerosol number, submicron sodium chloride mass, cloud condensation nuclei concentrations, and scattering coefficients. Generally, the relative response of aerosol properties to wind speed enhancement is strongest in fall when Arctic aerosol concentrations are lowest. Blowing snow events showed similar aerosol and environmental properties across events, apart from occasions with high snow age (>6 days since last snowfall). Coarse-mode number concentrations (>1 μm) are better explained by variability in wind speed averaged over 12-h air mass back trajectories arriving at the MOSAiC site compared to local, instantaneous wind speed, suggesting the importance of regional transport and consideration of air mass history for wind-driven aerosol production. These MOSAiC observations provide new insights into wind-driven aerosol in the central Arctic and may help validate modeling studies and improve model parameterizations particularly for aerosol direct and indirect radiative forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗