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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

A Composite Analysis of Snowfall Modes from Four Winter Seasons in Marquette, Michigan

Presented are four winter seasons of data from an enhanced precipitation instrument suite based at the National Weather Service (NWS) Office in Marquette (MQT), Michigan (250–500 cm of annual snow accumulation). In 2014 the site was augmented with a Micro Rain Radar (MRR) and a Precipitation Imaging Package (PIP). MRR observations are utilized to partition large-scale synoptically driven (deep)and surface-forced (shallow) snow events. Coincident PIP and NWS MQT meteorological surface observations illustrate different characteristics with respect to snow event category. Shallow snow events are often extremely shallow, with MRR-indicated precipitation heights of less than 1500 m above ground level. Large vertical reflectivity gradients indicate efficient particle growth, and increased boundary layer turbulence inferred from observations of spectral width implies increased aggregation in shallow snow events. Shallow snow events occur 2 times as often as deep events; however, both categories contribute approximately equally to estimated annual accumulation. PIP measurements reveal distinct regime-dependent snow microphysical differences, with shallow snow events having broader particle size distributions and comparatively fewer small particles and deep snow events having narrower particle size distributions and comparatively more small particles. In addition, coincident surface meteorological measurements indicate that most shallow snow events are associated with surface winds originating from the northwest (over Lake Superior), cold temperatures, and relatively high surface pressures, which are characteristics that are consistent with cold-air outbreaks. Deep snow events have meteorologically distinct conditions that are accordant with midlatitude cyclones and frontal structures, with mostly southwest surface winds, warmer temperatures approaching freezing, and lower surface pressures.

Claire Pettersen↗
DOE OSTI2019

Intercomparison and improvement of two-stream shortwave radiative transfer schemes in Earth system models for a unified treatment of cryospheric surfaces

Abstract. Snow is an important climate regulator because it greatly increases the surface albedo of middle and high latitudes of the Earth. Earth system models (ESMs) often adopt two-stream approximations with different radiative transfer techniques, the same snow therefore has different solar radiative properties depending whether it is on land or on sea ice. Here we intercompare three two-stream algorithms widely used in snow models, improve their predictions at large zenith angles, and introduce a hybrid model suitable for all cryospheric surfaces in ESMs. The algorithms are those employed by the SNow ICe and Aerosol Radiative (SNICAR) module used in land models, dEdd–AD used in Icepack, the column physics used in the Los Alamos sea ice model CICE and MPAS-Seaice, and a two-stream discrete-ordinate (2SD) model. Compared with a 16-stream benchmark model, the errors in snow visible albedo for a direct-incident beam from all three two-stream models are small (<±0.005) and increase as snow shallows, especially for aged snow. The errors in direct near-infrared (near-IR) albedo are small (<±0.005) for solar zenith angles θ<75∘, and increase as θ increases. For diffuse incidence under cloudy skies, dEdd–AD produces the most accurate snow albedo for both visible and near-IR (<±0.0002) with the lowest underestimate (−0.01) for melting thin snow. SNICAR performs similarly to dEdd–AD for visible albedos, with a slightly larger underestimate (−0.02), while it overestimates the near-IR albedo by an order of magnitude more (up to 0.04). 2SD overestimates both visible and near-IR albedo by up to 0.03. We develop a new parameterization that adjusts the underestimated direct near-IR albedo and overestimated direct near-IR heating persistent across all two-stream models for θ>75∘. These results are incorporated in a hybrid model SNICAR-AD, which can now serve as a unified solar radiative transfer model for snow in ESM land, land ice, and sea ice components.

58 GEOSCIENCES↗
DOE OSTI2026

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗
NASA NTRS2000

Coupled Model Simulation of Snowfall Events Over the Black Hills

Although many long-term simulations of snow accumulation and oblation have been made using stand-alone land surface models and surface models coupled with GCMs, less research has focused on short-term event simulations. Actually, accurate event simulations of snow-related processes are the basis for successful long-term simulation. Three advantages of event simulations of snowfall and snow melting are availability of: (1) intensive observation data from field experiments for validation; (2) more physically-realistic precipitation schemes for use in atmospheric models to simulate snowfall; and (3) a more detailed analysis of the snow melting processes. In addition to the complexities of snow related processes themselves, terrain-induced effects on snowfall/snow melting make simulations of snow events more difficult. Climatological observations indicate that terrain features such as the Black Hills of South Dakota and Wyoming can exert important effects on snow accumulation and snow oblation processes. One of the primary effects is that the orography causes forced uplift of airflow and causes atmospheric waves to form both upwind and downwind of it. Airflow often splits around the obstacle, converging on the lee side. This convergence may lead to precipitation enhancement. It also provides an elevated heat and moisture source that enhances atmospheric instability. During the period of April 5-May 5, 1999, the Upper Missouri River Basin Pilot Project (UMRBPP) made intensive observations on precipitation events occurring in the Black Hills. Two moderate snowfall events were captured during the period. The resulting high temporal and spatial resolution data provides opportunities to investigate terrain effects on snowfall amount, distribution, and melting. Successful simulation of snowfall amount, distribution, and evolution using atmospheric models is important to subsequent modeling of snow melting using snow sub-models in land surface schemes. In this paper, a coupled model system, consisting of an atmosphere model (ARPS) and a land-surface model (revised NCAR LSM), is used to simulate one of these UMRBPP cases.

Wang, Jianzhong↗
NASA NTRS2004

Clear-Sky Narrowband Albedo Variations Derived from VIRS and MODIS Data

A critical parameter for detecting clouds and aerosols and for retrieving their microphysical properties is the clear-sky radiance. The Clouds and the Earth's Radiant Energy System (CERES) Project uses the visible (VIS; 0.63 m) and near-infrared (NIR; 1.6 or 2.13 m) channels available on same satellites as the CERES scanners. Another channel often used for cloud and aerosol, and vegetation cover retrievals is the vegetation (VEG; 0.86- m) channel that has been available on the Advanced Very High Resolution Radiometer (AVHRR) for many years. Generally, clear-sky albedo for a given surface type is determined for conditions when the vegetation is either thriving or dormant and free of snow. Snow albedo is typically estimated without considering the underlying surface type. The albedo for a surface blanketed by snow, however, should vary with surface type because the vegetation often emerges from the snow to varying degrees depending on the vertical dimensions of the vegetation. For example, a snowcovered prairie will probably be brighter than a snowcovered forest because the snow typically falls off the trees exposing the darker surfaces while the snow on a grassland at the same temperatures will likely be continuous and, therefore, more reflective. Accounting for the vegetation-induced differences should improve the capabilities for distinguishing snow and clouds over different surface types and facilitate improvements in the accuracy of radiative transfer calculations between the snow-covered surface and the atmosphere, eventually leading to improvements in models of the energy budgets over land. This paper presents a more complete analysis of the CERES spectral clear-sky reflectances to determine the variations in clear-sky top-of-atmosphere (TOA) albedos for both snow-free and snow-covered surfaces for four spectral channels using data from Terra and Aqua.. The results should be valuable for improved cloud retrievals and for modeling radiation fields.

Sun-Mack, Sunny↗
NASA NTRS2013

Evaluation of Model Microphysics Within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is ~0.25 meters per second too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were ~0.25 meters per second too slow, while the SBU-YLIN was 0.25 to 0.5 meters per second too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗
NASA NTRS2013

Evaluation of Model Microphysics within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is approx 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were approx 0.25 m/s too slow, while the SBU-YLIN was 0.25 to 0.5 m/s too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗
NASA NTRS2019

Boreal canopy surfaces from spaceborne stereogrammetry

Surface elevation estimates from high resolution spaceborne image (HRSI) stereogrammetry are used to examine fine-scaled structure of boreal forest canopies. These data can depict detailed spatial patterns of vertical forest structure at remote sites across the circumpolar domain where these estimates would otherwise be unavailable. This work examines where these estimates are most effective at describing vertical forest structure to explain which canopy surfaces they represent. We evaluated the variation in canopy surface estimates captured from four general types of HRSI digital surface models (DSMs) across the full range of boreal canopy cover. These DSMs, classified into 4 types by grouping them according to the acquisition's (1) sun elevation angle (low or high) and (2) seasonality-driven ground surface condition (snow presence/absence), vary with acquisition characteristics and the details of this variation continues to be studied. We explored some of this variation by comparing the distributions of differences in boreal canopy percentile heights derived from reference small footprint lidar in Tanana Valley, Alaska with canopy surface elevations derived from these 4 types of HRSI DSMs. We examined how canopy surface estimates from HRSI DSMs differ according to acquisition characteristics and canopy cover, and ultimately which canopy surfaces are represented in these DSMs. Our results help clarify which boreal canopy surfaces are representative of those captured with HRSI DSMs. They show that in the Tanana Valley (1) DSMs grouped by sun elevation angle and ground surface condition provide different surface estimates of boreal canopies; (2) the two DSM types that appear to most differently capture boreal forest canopy surfaces are DSMs from snow-free images acquired at sun elevation angles <30° (Low sun elev. & snow-free) and those with snow-cover at sun elevation angles ≥30° (High sun elev. & snow-free); (3) DSMs with snow most often do not capture upper canopy surfaces; (4) the “Low sun elev. & snow-free” DSMs resolve surfaces that are most representative of upper canopy surfaces (dense forests >60% cover, 70th–80th percentile heights); and (5) in the most dense forests (>80% cover) where canopy gaps are least likely to bias downward the average surface estimates, the snow-free DSM types are representative of 70th - 80th percentile heights (“Low sun elev. & snow-free”) and 60th–70th percentile heights (“High sun elev. & snow-free”). The combination of horizontal structure (canopy cover) and acquisition characteristics affect the boreal vertical structure (canopy surface height) estimates from spaceborne stereogrammetry. These effects should be considered when analyzing products derived from HRSI DSMs, and as part of a comprehensive approach to spaceborne remote sensing of circumpolar boreal forests.

forest structure↗
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↗
NASA NTRS1980

Microwave remote sensing of snowpacks

The interaction mechanisms responsible for the microwave backscattering and emission behavior of snow were investigated, and models were developed relating the backscattering coefficient (sigma) and apparent temperature (T) to the physical parameters of the snowpack. The microwave responses to snow wetness, snow water equivalent, snow surface roughness, and to diurnal variations were investigated. Snow wetness was shown to have an increasing effect with increasing frequency and angle of incidence for both active and passive cases. Increasing snow wetness was observed to decrease the magnitude sigma and increase T. Snow water equivalent was also observed to exhibit a significant influence sigma and T. Snow surface configuration (roughness) was observed to be significant only for wet snow surface conditions. Diurnal variations were as large as 15 dB for sigma at 35 GHz and 120 K for T at 37 GHz. Simple models for sigma and T of a snowpack scene were developed in terms of the most significant ground-truth parameters. The coefficients for these models were then evaluated; the fits to the sigma and T measurements were generally good. Finally, areas of needed additional observations were outlined and experiments were specified to further the understanding of the microwave-snowpack interaction mechanisms.

Stiles, W. H.↗
NASA NTRS1990

Spectral albedo and emissivity of CO2 in Martian polar caps - Model results

In this paper, a snow albedo model previously developed for terrestrial snow is extended to the case of CO2 snow on Mars. Pure CO2 snow is calculated to have high albedo at visible wavelengths but not as high as that of water snow. At any given wavelength, the primary variable controlling albedo and emissivity is the snow grain size, with albedo decreasing and emissivity increasing as grain size increases. Observations that red albedo is much higher than blue albedo in the Martian south polar cap indicates that the snow or the atmosphere is contaminated with red dust. The absorption coefficient of CO2 ice in the thermal infrared is two to three orders of magnitude smaller than that measured for H2O ice. CO2 snow emissivity is therefore much lower than H2O snow, varying substantially with wavelength and quite sensitive to grain size and emission angle. Factors tending to increase emissivity are large grain size, small emission angle, and large concentrations of dust or water.

Warren, Stephen G.↗
NASA NTRS1999

Airborne Spectral Measurements of Surface-Atmosphere Anisotropy for Arctic Sea Ice and Tundra

Angular distributions of spectral reflectance for four common arctic surfaces: snow-covered sea ice, melt-season sea ice, snow-covered tundra, and tundra shortly after snowmelt were measured using an aircraft based, high angular resolution (1-degree) multispectral radiometer. Results indicate bidirectional reflectance is higher for snow-covered sea ice than melt-season sea ice at all wavelengths between 0.47 and 2.3 pm, with the difference increasing with wavelength. Bidirectional reflectance of snow-covered tundra is higher than for snow-free tundra for measurements less than 1.64 pm, with the difference decreasing with wavelength. Bidirectional reflectance patterns of all measured surfaces show maximum reflectance in the forward scattering direction of the principal plane, with identifiable specular reflection for the melt-season sea ice and snow-free tundra cases. The snow-free tundra had the most significant backscatter, and the melt-season sea ice the least. For sea ice, bidirectional reflectance changes due to snowmelt were more significant than differences among the different types of melt-season sea ice. Also the spectral-hemispherical (plane) albedo of each measured arctic surface was computed. Comparing measured nadir reflectance to albedo for sea ice and snow-covered tundra shows albedo underestimated 5-40%, with the largest bias at wavelengths beyond 1 pm. For snow-free tundra, nadir reflectance underestimates plane albedo by about 30-50%.

Arnold, G. Thomas↗
NASA NTRS2014

Evaluation of Mixed-Phase Microphysics Within Winter Storms Using Field Data and In Situ Observations

Snow prediction within models is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33-km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Cold-season Precipitation Experiment (GCPEx) experiment over southern Ontario, as well as a few years (12 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. The GCPex presentation will focus on verification using aircraft spirals through warm frontal snow band event on 18 February 2012. All the BMPs realistically simulated the structure of the band and the vertical distribution of snow/ice aloft, except the SBU-YLIN overpredicted slightly and Thompson (THOM) underpredicted somewhat. The Morrison (MORR) scheme produced the best slope size distribution for snow, while the Stony Brook (SBU) underpredicted and the THOM slightly overpredicted. Those schemes that have the slope intercept a function of temperature (SBU and WSM6) tended to perform better for that parameter than others, especially the fixed intercept in Goddard. Overall, the spread among BMPs was smaller than in other studies, likely because there was limited riming with the band. For the 15 cases at SBNY, which include moderate and heavy riming events, the non-spherical snow assumption (THOM and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM is approx. 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM, and MORR schemes were approx. 0.25 m/s too slow, while the SBU was 0.25 to 0.5 m/s too fast because of some excessive cloud water issues.

Colle, Brian A.↗
NASA NTRS2021

Maximum entropy modeling to identify physical drivers of shallow snowpack heterogeneity using unpiloted aerial system (UAS) lidar

Understanding the spatial variability of the snowpack is valuable for hydrologists and ecologists seeking to predict hydrological processes in a cold region. Snow distribution is a function of interactions among static variables, such as terrain, vegetation, and soil properties, and dynamic meteorological variables, such as solar radiation, wind speed and direction, and soil moisture. However, identifying the dominant physical drivers responsible for spatial patterns of the snowpack, particularly for ephemeral, shallow snowpacks, has been challenging due to the lack of the high-resolution snowpack and physical variables with high vertical accuracy as well as inherent limitations in traditional approaches. This study uses an Unpiloted Aerial System (UAS) lidar-based snow depth and static variables (1-m spatial resolution) to analyze field-scale spatial structures of snow depth and apply the Maximum Entropy (MaxEnt) model to identify primary controls over open terrain and forests at the University of New Hampshire Thompson Farm Research Observatory, New Hampshire, United States. We found that, among nine topographic and soil variables, plant functional type and terrain roughness contribute up to 80% and 76% of relative importance in the MaxEnt framework to predict locations of deeper or shallower snowpacks, respectively, across a mixed temperate forested and field landscape. Soil variables, such as organic matter and saturated hydraulic conductivity, were also important controls (up to 70% and 81%) on snow depth spatial variations for both open and forested landscapes suggesting spatial variations in soil variables under snow can control thermal transfer among soil, snowpack, and surface-atmosphere. This work contributes to improving land surface and snow models by informing parameterization of the sub-grid scale snow depths, down-scaling remotely sensed snow products, and understanding field scale snow states.

Eunsang Cho↗
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 OSTI2019

LLNL/Rescal-snow

When wind blows over snow, it self-organizes. This forms surface features, such as ripples and dunes, that alter the reflectivity and thermal conductivity of the snow. These features have just begun to be studied by the snow and climate science communities (see 1, 2, 3 for recent work). We created rescal-snow to provide a highly capable snow dune modelling toolkit, to enable snow scientists to study snow features in controlled numerical experiments, and produce high-quality quantitative output. We hope that this model will be useful to researchers in snow science, geomorphology, and polar climate.

Kochanski, KellyA↗
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↗