Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Snow”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Turbulent Pressure Spectrum Within the Roughness Sublayer of a Subarctic Forest Canopy

The turbulent static pressure spectrum E pp (k x ) as a function of longitudinal wavenumber k x in the roughness sublayer of forested canopies is of interest to a plethora of problems such as pressure transport in the turbulent kinetic energy budget, pressure pumping from snow or forest floor, and coupling between flow within and above canopies. Long term static pressure measurements above a sub-arctic forested canopy for near-neutral conditions during the winter and spring were collected and analyzed for three snow cover conditions: trees and ground covered with snow, trees are snow free but the ground is covered with snow, and snow free cover. In all three cases, it is shown that E pp (k x ) obeys the attached eddy hypothesis at low wavenumbers (1/δ < k x < 1/z) —with E pp (k x ) ∝ u$^4_*$k$^{-1}_x$ and Kolmogorov scaling in the inertial subrange at higher wavenumbers—with E pp (k x ) ∝ ϵ 4/3 k x -7/3 , where u * is the friction velocity at the canopy top, is the mean turbulent kinetic energy dissipation rate, z is the distance from the snow top, and is the boundary layer depth. The implications of these two scaling laws to the normalized root-mean squared pressure C p = o p /u$^2_*$ and its newly proposed logarithmic scaling with normalized wall-normal distance z/δ are discussed for snow covered and snow free vegetation conditions. The work here also shows that k$^{-1}_x$ the in the E pp (k x ) appears more extensive and robust than its longitudinal velocity counterpart.

Aslan, Toprak [Finnish Meteorological Inst., Helsi↗
DOE OSTI2022

Forest Canopy Density Effects on Snowpack Across the Climate Gradients of the Western United States Mountain Ranges

Limitations inherent to field studies such as limited spatial and temporal coverage has made it very challenging, if not impossible, to explicitly disentangle the effect of canopy density from climate on snowpack dynamics. Hence, our understanding of how variations in canopy density can affect snow processes under diverse climate conditions has been limited. To address this knowledge gap, this study uses a physics-based modeling approach to represent most key processes that drive snow accumulation and melt in the open and forest with varied forest density, across the climate gradients of the Western U.S. mountain ranges, as represented by 228 Snow Telemetry (SNOTEL) locations. Simulations suggest that, under most winter climate conditions, canopy density exerts a strong influence on under-canopy snowpack dynamics and their controlling processes. Changing canopy density can shift the direction of the modeled differential peak snow accumulation and snow disappearance date (SDD) between the open and forest. For example, while simulated snowpack lasts longer in the open than the high-density forest at most wet/warm and dry/cold locations by up to 55 days, about 90% of all locations show a longer snow duration in the low-density forest by up to 27 days as compared to the open. Interannual variability in climate can also change the magnitude and direction of canopy impact on snow processes. By comparing model simulations to published field observations, and we also identify conditions that require additional field observations to improve process-scale understanding of forest-snow interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assessing Prior Emergent Constraints on Surface Albedo Feedback in CMIP6

An emergent constraint (EC) is a popular model evaluation technique, which offers the potential to reduce intermodel variability in projections of climate change. Two examples have previously been laid out for future surface albedo feedbacks (SAF) stemming from loss of Northern Hemisphere (NH) snow cover (SAF snow ) and sea ice (SAF ice ). These processes also have a modern-day analog that occurs each year as snow and sea ice retreat from their seasonal maxima, which is strongly correlated with future SAF across an ensemble of climate models. The newly released CMIP6 ensemble offers the chance to test prior constraints through out-of-sample verification, an important examination of EC robustness. Here, we show that the SAF snow EC is equally strong in CMIP6 as it was in past generations, while the SAF ice EC is also shown to exist in CMIP6, but with different, slightly weaker characteristics. We find that the CMIP6 mean NH SAF exhibits a global feedback of 0.25 ± 0.05 W m -2 K -1 , or ~61% of the total global albedo feedback, largely in line with prior generations despite its increased climate sensitivity. Additionally, the NH SAF can be broken down into similar contributions from snow and sea ice over the twenty-first century in CMIP6. Crucially, intermodel variability in seasonal SAF snow and SAF ice is largely unchanged from CMIP5 because of poor outlier simulations of snow cover, surface albedo, and sea ice thickness. These outliers act to mask the noted improvement from many models when it comes to SAF ice , and to a lesser extent SAF snow .

58 GEOSCIENCES↗
DOE OSTI2024

Bridging gaps in permafrost-shrub understanding

Permafrost, permanently frozen ground which underlies much of the Arctic and sub-Arctic, is influenced by deciduous shrubs, yet our understanding of permafrost-shrub interactions is limited. This is largely due to a lack of widespread, long-term, high-quality observations across Arctic and sub-Arctic systems, which are difficult to study due to their remote locations. Shrubs are rapidly expanding in many areas of Arctic tundra and can either amplify or inhibit permafrost thaw, making it crucial that we can understand and predict permafrost-shrub interactions. As we have limited time and resources to make observations, we must design our field campaigns to be as impactful as possible. Below we outline the current state of knowledge and give suggestions about how we can maximize our fieldwork’s impact. From field studies conducted around the Arctic, we know that shrubs generally have a cooling effect on permafrost in the summer, largely caused by the shrub canopy shading the soil. However, summertime cooling can be counteracted by a warming effect in winter, caused by the trapping of blowing snow by shrubs, which leads to increased snow depths. In addition to these commonly observed effects, both observations and models show that other local factors such as climate, soil, snow, and disturbances can interact to cause contrasting dominant effects on permafrost temperatures. For example, shrubs that protrude through the snowpack have been observed to cool permafrost in winter via thermal bridging between the air above the snow and the soil below. However, at other research sites, protruding shrubs lower the surface albedo and cause earlier snowmelt in spring, amplifying permafrost thaw by exposing the soil to thawing earlier. Protruding shrubs have also been observed to induce snow melt in autumn, creating ice layers at the snow surface that prevent further snow accumulation and eliminate the ability of shrubs to trap blowing snow, thereby reducing the winter warming effect. While we understand the physical processes that lead to these effects, we cannot explain why we observe different effects between sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hydrologic Model Files (PRMS) for Historical Conditions in the East River Watershed, Colorado between 1987-2019

The data package contains model input and output file for the East River watershed (85 km2) located in the headwaters of the Colorado River. The model is the U.S. Geological Survey (USGS) Precipitation-Modeling Runoff System (PRMS) for the historical, or baseline, condition (1987-2019) to assess the efficiency of summer monsoon rains to generate streamflow. The finite difference grid resolution is 100 m with elevations resampled from the USGS National Elevation Dataset. LANDFIRE raster data is used to derive parameters of dominant cover type, summer and winter cover density, canopy interception characteristics for snow and rain and transmission coefficient for shortwave radiation. Climate forcing uses minimum and maximum daily temperature lapse rates defined by the two Snow Telemetry (SNOTEL) stations located in the vicinity of the East River adjusted for aspect. Schofield SNOTEL snowfall is spatially distributed using LiDAR derived snow depth observations from the Airborne Snow Observatory (ASO) flown 4 April 2016. Snow depths are converted to snow water equivalent based on ground surveys and density modeling. Rainfall is spatially distributed using the monthly Parameter-elevation Relationships on Independent Slopes Model 30-year monthly averages. Simulated solar radiation is calibrated to match weather station observations. Model verification of snow water equivalent (SWE) accumulation and ablation relies on the 2018 and 2019 ASO maps and weather station snow depth. The data package contains model input files in text format, model documentation in pfd and a readme.txt file to help users understand the file structure to run the program.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Xanthos-Lake Model Source Code

This repository contains the source code for Xanthos-Lake, a lake-modeling extension of the Xanthos framework that introduces a coupled lake component comprising the Xanthos-Lake Snow and Ice Model (xLSIM) and the Xanthos-Lake Water Balance Model (xLWBM). xLSIM is a basin-aware machine-learning model for lake snow, ice, and thermal conditions. It predicts monthly lake ice thickness, snow depth, snow-cover fraction, mixing-layer temperature, and lake ice fraction from meteorological forcing and lake surface-area information. It uses sequence-based deep-learning architectures, including Transformer and hybrid Long Short-Term Memory–Transformer (LSTM–Transformer) models, together with seasonal encoding, multi-lake learning, physical masking, and basin-level cryospheric and non-cryospheric classification. The training workflow uses Ray for scalable execution and includes optional Ray Tune hyperparameter optimization. Model predictions, observations, diagnostics, and feature-importance outputs are written in NetCDF. xLWBM is the water-balance component of the new lake framework. It simulates monthly lake storage, surface area, evaporation, inflow, outflow, and lake–groundwater exchange. It combines physical water-balance equations with calibrated bathymetric relationships, weir-based outlet flow, modified Penman open-water evaporation, groundwater head relaxation, Penman–Monteith snow and ice sublimation, and snow, ice, and thermal conditions supplied by xLSIM. The model calibrates lake parameters against satellite-derived surface-area data, using evaporation-based calibration where surface-area data are unavailable, and supports small, medium, and large lake classes. For large lakes, xLWBM is integrated with the managed-routing workflow so that lake storage and outflow interact directly with downstream river routing and reservoir operations. Together, xLSIM and xLWBM provide Xanthos with a coupled lake-modeling capability. xLSIM supplies the snow, ice, and thermal conditions that affect lake evaporation and snow- and ice-related water exchanges, while xLWBM translates those conditions into dynamic lake storage, surface area, evaporation, and discharge. In return, xLWBM supplies evolving lake surface area to xLSIM. This coupling enables Xanthos to represent lakes as active hydrologic components within basin-scale water-availability and routing simulations.

Machine Learning↗
DOE OSTI2020

Light-absorbing impurities accelerating glacial melting in southeastern Tibetan Plateau

Deposition of light-absorbing particles on glacier surfaces poses a series of adverse impacts on the cryospheric environment, climate and human health. Broad attention of the scientific community has been paid on insoluble light-absorbing impurities (ILAIs) in snow and ice on glaciers over the Tibetan Plateau (TP). However, systematic investigation of ILAIs in snowpack of glaciers on the TP is scarce. In this study, the properties and darkening effect of ILAIs in snowpack on glaciers are extensively investigated in the southeast of TP. Results show that ILAIs concentrations in multiple types of snow and ice samples were significantly different. Snowpit depths varied substantially from one profile to another during May and June 2016. The average concentrations of ILAIs in snowpits increase as snow melting progresses. Black carbon (BC) and dust cause snow albedo reduction more in snow with larger grain size Re. Based on a radiative transfer model calculation, the average albedo reduction induced by BC in the snowpack was 0.141 ± 0.02, and associated daily maximum radiative forcing (RF) was 72.97 ± 12.7 W m-2. BC is a controlling light-absorbing factor in snowpack and causes substantial albedo reduction and thus the associated daily maximum RF. The maximum reduction of snow cover duration was 4.56 ± 0.71 days caused by BC and dust in snowpack in southeastern TP. The average mass absorption cross-section (MAC) of BC from multiple snowpits was 3.26 ± 0.46 m2 g-1, which represents a typical value of MAC in snow on glaciers, but it is type-dependent of snow/ice samples. Tropospheric aerosols vertically extended up to 8 km over the TP and its surrounding areas, which indicates the transport of aerosols from remote sources through elevated pathways. A large amount of carbon stored in the brittle glaciers can be potentially released with meltwater runoff under a warming climate. This study provides a new insight for investigating carbonaceous and light-absorbing particles in glacierization areas.

Niu, Hewen↗
DOE OSTI2020

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Theory for the Balance between Warm Rain and Ice Crystal Processes of Precipitation in Mixed-Phase Clouds

Abstract Mixed-phase clouds contain both supercooled cloud liquid and ice crystals. In principle, precipitation may be initiated either by the liquid phase or by the ice phase. Ice crystals may grow by vapor diffusion to become snow (“ice crystal process”), forming “cold” precipitation. Equally, cloud droplets, when large enough, coalesce to form “warm” precipitation by the “warm rain process.” Warm rain could be supercooled and freeze as “warm” graupel. In the present paper, a new simplified theoretical analysis is provided to examine the microphysical system consisting of three species of hydrometeor, namely, cloud liquid, “cold ice” (crystals, snow), and “warm rain” (frozen or supercooled). This is obtained by nondimensionalizing and simplifying the evolution equations for the mass of each species. Analytical formulas are given for equilibria. Feedback analysis shows that the sign of the feedback is linked to the abundance of precipitation, with a neutral surface in the 3D phase space. The system’s precipitation amount explodes while in the initial unstable regime, crossing the neutral surface and approaching the equilibrium point that is a stable attractor. Positive and negative feedbacks are elucidated. In a standard case, the cold ice mass is about 1000 times larger than the warm rain mass. To illustrate the physical behavior of the theory, sensitivity tests are performed with respect to environmental conditions (e.g., aerosol, updraft speed) and microphysical parameters (e.g., riming and sedimentation rates for cold ice). Cold ice prevails, especially in fast ascent, due to its low bulk density, favoring slow sedimentation and a wide cross-sectional area for riming. Significance Statement The theory elucidates how the ice phase can prevail in the precipitation from any mixed-phase clouds with supercooled cloud liquid and crystals. The ice phase radically suppresses cloud liquid by riming when active and “wins” the competition against coalescence. This prevalence of ice is shown to arise from the low bulk density of snow. The cloud is viewed as a system of negative and positive feedbacks that prevail in realms of stability and instability in a 3D phase space.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Sea ice heat and mass balance measurements from four autonomous buoys during the MOSAiC drift campaign

As part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC), four autonomous seasonal ice mass balance buoys were deployed in first- and second-year ice. These buoys measured position, barometric pressure, snow depth, ice thickness, ice growth, surface melt, bottom melt, and vertical profiles of temperature from the air, through the snow and ice, and into the upper ocean. Observed air temperatures were similar at all four sites; however, snow–ice interface temperatures varied by as much as 10°C, primarily due to differences in snow depth. Observed winter ice growth rates (November to May) were <1 cm day -1 , with summer melt rates (June to July) as large as 5 cm day -1 . Air temperatures changed as much as 2°C hour -1 but were dampened to <0.3°C hour -1 at the snow–ice interface. Initial October ice thicknesses ranged from 0.3 m in first-year ice to 1.2 m in second-year ice. By February, this range was only 1.20–1.46 m, due in part to differences in the onset of basal freezing. In second-year ice, this delay was due to large brine-filled voids in the ice; propagating the cold front through this ice required freezing the brine. Mass balance results were similar to those measured by autonomous buoys deployed at the North Pole from 2000 to 2013. Winter average estimates of the ocean heat flux ranged from 0 to 3 W m -2 , with a large increase in June 2020 as the floe moved into warmer water. Estimates of average snow thermal conductivity measured at two buoys during periods of linear temperature profiles were 0.41 and 0.42 W m -1 °C -1 , higher than previously published estimates. Results from these ice mass balance buoys can contribute to efforts to close the MOSAiC heat budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A parameterization of sub-grid topographical effects on solar radiation in the E3SM Land Model (version 1.0): implementation and evaluation over the Tibetan Plateau

Abstract. Topography exerts significant influences on the incoming solar radiation at the land surface. A few stand-alone regional and global atmospheric models have included parameterizations for sub-grid topographic effects on solar radiation. However, nearly all Earth system models (ESMs) that participated in the Coupled Model Intercomparison Project (CMIP6) use a plane-parallel (PP) radiative transfer scheme that assumes that the terrain is flat. In this study, we incorporated a well-validated sub-grid topographic (TOP) parameterization in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) version 1.0 to quantify the effects of sub-grid topography on solar radiation flux, including the shadow effects and multi-scattering between adjacent terrain. We studied the role of sub-grid topography by performing ELM simulations with the PP and TOP schemes over the Tibetan Plateau (TP). Additional ELM simulations were performed at multiple spatial resolutions to investigate the role of spatial scale on sub-grid topographic effects on solar radiation. The Moderate Resolution Imaging Spectroradiometer (MODIS) data was used to compare with the ELM simulations. The results show that topography has non-negligible effects on surface energy budget, snow cover, snow depth, and surface temperature over the TP. The absolute differences in surface energy fluxes for net solar radiation, latent heat flux, and sensible heat flux between TOP and PP exceed 20, 10, and 5 W m−2, respectively. The differences in land surface albedo, snow cover fraction, snow depth, and surface temperature between TOP and PP exceed 0.1, 0.1, 10 cm, and 1 K, respectively. The magnitude of the sub-grid topographic effects is dependent on seasons and elevations and is also sensitive to the spatial scales. Although the sub-grid topographic effects on solar radiation are larger with more spatial details at finer spatial scales, they cannot be simply neglected at coarse spatial scales. When compared to MODIS data, incorporating the sub-grid topographic effects overall reduces the biases of ELM in simulating surface energy balance, snow cover, and surface temperature, especially in the high-elevation and snow-covered regions over the TP. The inclusion of sub-grid topographic effects on solar radiation parameterization in ELM will contribute to advancing our understanding of the role of the surface topography on terrestrial processes over complex terrain.

Hao, Dalei (ORCID:0000000334979774)↗
DOE OSTI2025

Seasonality and Albedo Dependence of Cloud Radiative Forcing in the Upper Colorado River Basin

Mountains create and enhance their own clouds, which both scatter and absorb shortwave radiation from the sun and absorb and re-emit land surface and atmospheric longwave radiation. However, the impacts of clouds on the surface radiation balance in high elevation snowy mountain terrain are poorly explored. In this study, we use data collected by the SAIL field campaign and partner organizations in the upper elevations (2,880 m.a.s.l) of the Upper Colorado River Basin (UCRB) over a 21-month period from September 2021 to June 2023 to estimate Cloud Radiative Forcing (CRF) in the shortwave, longwave, and the net effect. Longwave warming effects dominate during the winter when snow albedos are high (0.8–0.9) and the background atmospheric precipitable water vapor is low (<0.5 cm), yielding a maximum monthly average net CRF of +34.7 W·m -2 , meaning that clouds increase the net radiation relative to clear skies during this time period. The sign of net CRF switches in the warm season as snow recedes, sun-angles increase, and the North American monsoon arrives, yielding a minimum monthly average net CRF of -47.6 W·m -2 with hourly minima of -600 W·m -2 . The sign of net CRF is typically positive, even at solar noon, when the surface is snow covered, except for a brief period over melting, low-albedo snow (0.5–0.6) impacted by dust impurities. Sensitivity tests elucidate the role of the surface albedo on the net CRF. The results suggest that net CRF will increase in magnitude and lead to a more persistent cooling effect on the surface net radiation budget as the snow cover declines.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Daily, 30 m Resolution NDSI Data for the East River Watershed, CO for 2000-2020

This dataset contains daily Normalized Difference Snow Index (NDSI) values at 30 m spatial resolution for the East River watershed in Colorado, USA. The temporal range of these data includes water years 2001-2020. These data were created using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). This model fuses low spatial and high temporal resolution data from MODIS (500 m, daily) with high spatial and low temporal resolution data from Landsat (30 m, 16 days) to create a 30m synthetic daily snow product. This product allows for the analysis of historical snow covered area trends in the East River Watershed at fine spatiotemporal resolutions where it was not available previously. This research was performed as a part of the Department of Energy’s Subsurface Biogeochemical Research Program with the primary intent of better understanding the timing and spatial patterns of water delivery to the Critical Zone in mountain watersheds. Each .zip file contains one "water year" of data (October 1 - September 30; i.e., water year 2010 starts October 1, 2010 and ends September 30, 2011). Each zip file contains the following: STARFM daily Normalized Difference Snow Index (NDSI) fusion data files in GeoTiff format with one layer for each day between Landsat data acquisition dates (i.e., for dates of Landsat acquisition, the Landsat image is included for that date). The study area is located in an area of Landsat path overlap, so Landsat dates acquisitions are every 7-9 days. Landsat NDSI files containing the high spatial (30m), low temporal (7-9 days due to Landsat path overlap) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Dates for which no Landsat data were obtained are included as NoData layers. MODIS NDSI files containing the high temporal (daily), low spatial (500m) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Please note the MODIS data were resampled to 30m pixels for input into the STARFM model. The data have a scale factor of 10,000 and a no data value of -32767. The projection of all datasets is WGS 84 (EPSG: 4326), which has a latitude/longitude based degree resolution of 0.0002694946 X 0.0002694946, and approximates to the 30 m spatial resolution mentioned above. The Layer Index files in .csv format. They contain information for each layer in the above GeoTiff files regarding the corresponding date for each layer, the fraction of pixels in the image that contain valid data (missing data is due to either cloud cover or poor data quality; these values are not percent snow cover). Dates of Landsat overpass are indicated in these files. If no Landsat data were able to be obtained due to cloud cover or lack of Landsat Tier 1 data available on Google Earth Engine, this is also noted.

EARTH SCIENCE > CRYOSPHERE > SNOW/ICE↗
DOE OSTI2023

Midwinter Dry Spells Amplify Post‐Fire Snowpack Decline

Abstract Increasing wildfire and declining snowpacks in mountain regions threaten water availability. We combine satellite‐based fire detections with snow seasonality classifications to examine fire activity in California's seasonal and ephemeral snow zones. We find a nearly tenfold increase in fire activity during 2020–2021 versus 2001–2019. Accumulation season broadband snow albedo declined 25%–71% at two burned sites (2021 and 2022) according to in‐situ data relative to un‐burned conditions, with greater declines associated with increased burn severity. By enhancing snowpack susceptibility to melt, both decreased snow albedo and canopy drove midwinter melt during a multi‐week dry spell in 2022. Despite similar meteorological conditions in December–February 2013 and 2022–linked to persistent high pressure weather regimes–minimal melt occurred in 2013. Post‐fire snowpack differences are confirmed with satellite measurements. With growing geographical overlap between wildfire and snow, our findings suggest California's snowpack is increasingly vulnerable to the compounding effects of dry spells and wildfire.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The understudied winter: Evidence of how precipitation differences affect stream metabolism in a headwater

Climate change is causing pronounced shifts during winter in the US, including shortening the snow season, reducing snowpack, and altering the timing and volume of snowmelt-related runoff. These changes in winter precipitation patterns affect in-stream freeze-thaw cycles, including ice and snow cover, and can trigger direct and indirect effects on in-stream physical, chemical, and biological processes in ~60% of river basins in the Northern Hemisphere. We used high-resolution, multi-parameter data collected in a headwater stream and its local environment (climate and soil) to determine interannual variability in physical, chemical, and biological signals in a montane stream during the winter of an El Niño and a La Niña year. We observed ~77% greater snow accumulation during the El Niño year, which caused the formation of an ice dam that shifted the system from a primarily lotic to a lentic environment. Water chemistry and stream metabolism parameters varied widely between years. They featured anoxic conditions lasting over a month, with no observable gross primary production (GPP) occurring under the ice and snow cover in the El Niño year. In contrast, dissolved oxygen and GPP remained relatively high during the winter months of the La Niña year. These redox and metabolic changes driven by changes in winter precipitation have significant implications for water chemistry and biological functioning beyond the winter. Our study suggests that as snow accumulation and hydrologic conditions shift during the winter due to climate change, hot-spots and hot-moments for biogeochemical processing may be reduced, with implications for the downstream movement of nutrients and transported materials.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impacts of climate and disturbance on nutrient fluxes and stoichiometry in mixed-conifer forests

Elucidating climatic impacts on stream nutrient export and stoichiometry will improve the understanding of forest carbon (C) storage in a warmer world. We constructed C, nitrogen (N), and phosphorus (P) budgets in four watersheds within a rain-snow transition site and another four within a higher-elevation, snow-dominated site, in California’s mixed-conifer zone. We used these two sites in a space-for-time substitution to assess the potential effects of warming on nutrient cycles in currently snow-dominated areas that will become more rain-dominated. During a non-drought period (water year (WY) 2004 – 2011), mean annual stream exports of C and N in particulate forms at the transition site were double that at the snow-dominated site, suggesting that sediment-associated nutrient losses may increase with warming. The transition site had 12% lower N but double the P content in mineral horizons, lower N:P mass ratios in organic horizons, and lower stream export of dissolved inorganic N than the snow-dominated site. Furthermore, these differences suggest that montane forests may receive lower supply of N relative to P with warming. In addition, given strong interests in forest thinning to increase drought resiliency, we examined changes in stream nutrient export after thinning and during a major drought period (WY 2013 – 2015). Stream exports of C, N, and P were similar between unthinned and thinned watersheds during drought, suggesting negligible thinning impacts on stream nutrient export during excessively dry periods. Taken together, our results suggest that as the climate warms, California’s montane forests may lose more nutrients through erosion and receive less N input than P.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improved Convective Ice Microphysics Parameterization in the NCAR CAM Model

Partitioning deep convective cloud condensates into components that sediment and detrain, known to be a challenge for global climate models, is important for cloud vertical distribution and anvil cloud formation. In this study, we address this issue by improving the convective microphysics scheme in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 (CAM5.3). The improvements include: (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. Results from the default and modified CAM5.3 models were evaluated against observations from the U.S. Department of Energy Tropical Warm Pool-International Cloud Experiment (TWP-ICE) field campaign. The default model overestimates ice amount, which is largely attributed to the underestimation of convective ice particle sedimentation. By considering cloud ice sedimentation and rimed ice particles and applying a new convective snow terminal velocity parameterization, the vertical distribution of ice amount is much improved in the midtroposphere and upper troposphere when compared to observations. The vertical distribution of ice condensate also agrees well with observational best estimates upon considering snow detrainment. Further, comparison with observed convective updrafts reveals that current bulk model fails to reproduce the observed updraft magnitude and occurrence frequency, suggesting spectral distributions be required to simulate the subgrid updraft heterogeneity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Photolytic Degradation of Water‐Soluble Organic Carbon in Snowmelts: Changes in Molecular Characteristics, Brown Carbon Chromophores, and Radiative Effects

Abstract Water‐soluble organic carbon (WSOC) deposited in ambient snowpack play key roles in regional carbon cycle and surface energy budget, but the impacts of photo‐induced processes on its optical and chemical properties are poorly understood yet. In this study, melted samples of the seasonal snow collected from northern Xinjiang, northwestern China, were exposed to ultraviolet (UV) radiation to investigate the photolytic transformations of WSOC. Molecular characteristics and chemical composition of WSOC and its brown carbon (BrC) constituents were investigated using high‐performance liquid chromatography interfaced with a photodiode array detector and a high‐resolution mass spectrometer. Upon illumination, formation of nitrogen‐ and sulfur‐containing species with high molecular weight was observed in snow samples influenced by soil‐ and plant‐derived organics. In contrast, the representative sample collected from remote region showed the lowest molecular diversity and photolytic reactivity among all samples, in which no identified BrC chromophores decomposed upon illumination. Approximately 65% of chromophores in urban samples endured UV irradiation. However, most of BrC composed of phenolic/lignin‐derived compounds and flavonoids disappeared in the illuminated samples containing WSOC from soil‐ and plant‐related sources. Effects of the photochemical degradation of WSOC on the potential modulation of snow albedo were estimated. Apparent half‐lives of WSOC estimated as albedo reduction in 300–400 nm indicated 0.1–0.4 atmospheric equivalent days, which are shorter than typical photolysis half‐lives of ambient biomass smoke aerosol. This study provides new insights into the roles of WSOC in snow photochemistry and snow surface energy balance.

Meteorology & Atmospheric Sciences↗