Search NASA⌕ Search

SEARCH · Search NASA

Results for “Snow model”

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

General Report of the Researches of Snowpack Properties, Snowmelt Runoff and Evapotranspiration in Japan

A method was developed for estimating the distribution of snow and the snow water equivalent in Japan by combining LANDSAT data with the degree day method. A snow runoff model was improved and applied to the Okutadami River basin. The Martinec Rango model from the U.S. was applied to Japanese river basins to verify its applicability. This model was then compared with the Japanese model. Analysis of microwave measurements obtained by a radiometer on a tower over dry snow in Hokkaido indicate a certain correlation between brightness temperature and snowpack properties. A correlation between brightness temperature and depth of dry snow in an inland plain area was revealed in NIMBUS SMMR data obtained from the U.S. Calculation of evaporation using airborne remote sensing data and a Priestley-Taylor type of equation shows that the differentiation of evaporation with vegetation type is not remarkable because of little evapotransportation in winter.

Takeda, K.↗

Interactions Between Thresholds and Spatial Discretizations of Snow: Insights From Estimates of Wolverine Denning Habitat in the Colorado Rocky Mountains

Thresholds can be used to interpret environmental data in a way that is easily communicated and useful for decision making purposes. However, thresholds are often developed for specific data products and time periods, changing findings when the same threshold is applied to datasets or periods with different characteristics. Here, we test the impact of different spatial discretizations of snow on annual estimates of wolverine denning opportunities in the Colorado Rocky Mountains, defined using a snow water equivalent (SWE) threshold (0.20 m) and threshold date (15 May) from previous habitat assessments. Annual potential wolverine denning area (PWDA) was thresholded from a 36-year (1985 – 2020) snow reanalysis model with three different spatial discretizations: 1) 480 m grid cells (D480), 2) 90 m grid cells (D90), and 3) 480 m grid cells with implicit representations of subgrid snow spatial heterogeneity (S480). Relative to the D480 and S480 discretizations, D90 resolved shallower snow deposits on slopes between 3050 and 3350 m elevation, decreasing PWDA by 10%, on average. In years with warmer and/or drier winters, S480 discretizations with subgrid representations of snow heterogeneity increased PWDA, even within grid cells where mean 15 May SWE was less than the SWE threshold. These simulations increased PWDA by upwards of 30% in low snow years, as compared to the D480 and D90 simulations without subgrid snow heterogeneity. Despite PWDA sensitivity to different snow spatial discretizations, PWDA was controlled more by annual variations in winter precipitation and temperature. However, small changes to the SWE threshold (± 0.07 m) and threshold date (± 2 weeks) also affected PWDA by as much as 82%. Across these threshold ranges, PWDA was approximately 18% more sensitive to the SWE threshold than the threshold date. However, the sensitivity to the threshold date was larger in years with late spring snowfall, when PWDA depended on whether modeled SWE was thresholded before, during, or after spring snow accumulation. Our results demonstrate that snow thresholds are useful but may not always provide a complete picture of the annual variability in snow-adapted wildlife denning opportunities. Studies thresholding spatiotemporal datasets could be improved by including 1) information about the fidelity of thresholds across multiple spatial discretizations, and 2) uncertainties related to ranges of realistic thresholds.

Justin M. Pflug↗

How do tradeoffs in satellite spatial and temporal resolution impact snow water equivalent reconstruction?

Given the tradeoffs between spatial and temporal resolution, questions about resolution optimality are fundamental to the study of global snow. Answers to these questions will inform future scientific priorities and mission specifications. Heterogeneity of mountain snowpacks drives a need for daily snow cover mapping at the slope scale (≤30 m) that is unmet for a variety of scientific users, ranging from hydrologists to the military to wildlife biologists. But finer spatial resolution usually requires coarser temporal or spectral resolution. Thus, no single sensor can meet all these needs. Recently, constellations of satellites and fusion techniques have made noteworthy progress. The efficacy of two such recent advances is examined: (1) a fused MODIS–Landsat product with daily 30 m spatial resolution and (2) a harmonized Landsat 8 and Sentinel 2A and B (HLS) product with 3–4 d temporal and 30 m spatial resolution. State-of-the-art spectral unmixing techniques are applied to surface reflectance products from 1 and 2 to create snow cover and albedo maps. Then an energy balance model was run to reconstruct snow water equivalent (SWE). For validation, lidar-based Airborne Snow Observatory SWE estimates were used. Results show that reconstructed SWE forced with 30 m resolution snow cover has lower bias, a measure of basin-wide accuracy, than the baseline case using MODIS (463 m cell size) but greater mean absolute error, a measure of per-pixel accuracy. However, the differences in errors may be within uncertainties from scaling artifacts, e.g., basin boundary delineation. Other explanations are (1) the importance of daily acquisitions and (2) the limitations of downscaled forcings for reconstruction. Conclusions are as follows: (1) spectrally unmixed snow cover and snow albedo from MODIS continue to provide accurate forcings for snow models and (2) finer spatial and temporal resolution through sensor design, fusion techniques, and satellite constellations are the future for Earth observations, but existing moderate-resolution sensors still offer value.

Edward H. Bair↗

Scattering from a random layer with applications to snow, vegetation and sea ice

Two approaches for computing scattering from a random layer with an irregular interface are shown using the radiative-transfer principle. One approach is applied to a random layer to develop a scattering model for snow and sea ice, while the other is used to generate a scattering model for vegetation. It is noted that to model the scattering characteristics of a special medium, it is necessary to relate the electromagnetic parameters to the measurable parameters of the scatterers in the medium. Such relations are given and used to calculate theoretical estimates to compare with measurements acquired by microwave scatterometers.

Fung, A. K.↗

A Wind-Driven Snow Redistribution Module for Alpine3d V3.3.0: Adaptations Designed for Downscaling Ice Sheet Surface Mass Balance

Ice sheets gain mass via snow accumulation at the ice sheet surface, which is the primary component of surface mass balance. On the Antarctic ice sheet, winds redistribute snow resulting in surface mass balance that is variable in both space and time. Representing wind-driven snow redistribution processes in models is critical for local assessments of surface mass balance, repeat altimetry studies, and interpretation of ice core accumulation records. To this end, we have adapted Alpine3D, an existing distributed snow modeling framework, to downscale Antarctic surface mass balance to horizontal resolutions up to 1 km. In particular, we have introduced a new two-dimensional advection-based wind-driven snow redistribution module that is driven by an offline coupling between WindNinja, a wind downscaling model, and Alpine3D. We then show that large accumulation variability can be at least partially explained by terrain-induced wind speed variations which subsequently redistribute snow around rolling topography. By comparing Alpine3D to airborne-derived snow accumulation measurements within a testing domain over Pine Island Glacier in West Antarctica, we demonstrate that our Alpine3D downscaling approach improves surface mass balance estimates when compared to MERRA-2, a global atmospheric reanalysis which we use as atmospheric forcing. In particular, when compared to MERRA-2, Alpine3D reduces simulated surface mass balance root mean squared error by 23.4 mm w.e.yr−1 (13%) and increases variance explained by 24%. Despite these improvements, Alpine3D still underestimates observed accumulation variability, thus providing an opportunity for future model improvement.

Eric Keenan↗

Estimating Alpine Snow Depth by Combining Multi-Frequency Passive Radiance Observations with Ensemble Snowpack Modeling

This paper presents a physically-based snow depth retrieval algorithm adapted for deep mountainous snowpack and airborne multifrequency (10.7, 18.7, 37.0 and 89.0 GHz) passive microwave (PM) radiance observations from a single flight. The algorithm employs a single forecast-analysis cycle of a traditional sequential assimilation scheme. It uses an ensemble of multi-layer snowpack model runs to resolve snow microstructure and melt-refreeze crusts, andmicrowave radiative transfer models to relate snow properties to microwave measurements.Snow depth was retrieved at a 120 m spatial resolution over three 1 km2 Intensive Study Area (ISA) within the Rabbit Ears Meso-Cell Study Area (MSA) from the NASA Cold Land Processes Experiment (CLPX) in Colorado (United States) for one date in February 2003. When evaluated against in situ observations, root mean square error (RMSE) of the snow depth from the assimilation was 13.3 cm for areas with low (<5 %) forest cover, which was a reduction of 48 % in the RMSE compared with the modeled snow depth when the PM observations were not assimilated; indicating a ~5% relative error of the posterior snow depth with respect to the average snow depth (200 cm) measured at these pixels. For pixels with forest cover ranging from 5-15% and 15-30%, results were improved (R2 increased from 0.65 to 0.71 and from 0 to 0.38, respectively) by introducing a forest radiative transfer model during the assimilation.

Airborne Multifrequency↗

Evaluation of Brightness Temperature Sensitivity to Snowpack Physical Properties Using Coupled Snow Physics and Microwave Radiative Transfer Models

There are multiple existing microwave radiative transfer models (RTMs) to simulate the brightness temperature (Tb) of snowpacks. It is still challenging to have consistent Tb responses from RTMs due to individual physical formulations of the snowpack scattering process. This article examines three of the widely-used multi-layer RTMs: 1) the microwave emission model of layered snowpacks (MEMLS); 2) the dense media radiative transfer based on the quasi-crystalline approximation (QCA) of Mie scattering of densely packed sticky spheres (DMRTQMS); and 3) the Helsinki University of Technology (HUT) model. Interestingly, these models yield slightly different Tb responses when driven by the same physical snowpack properties. Tb variations, dependent on the choice of RTMs, are then evaluated to improve the understanding of model differences in microwave emission from a snowpack. We first perform a sensitivity study of the Tb predictions from the three RTMs as a function of snow grain sizes, densities, and depths. While Tb from all three RTMs decreases with increasing snow grain sizes, it is found that a scaling factor is required to have the same amount of Tb attenuation for small grain sizes within the Rayleigh scattering regime. For larger grain sizes, however, a scaling coefficient is not enough to match the model outputs due to the different scattering assumptions of the RTMs. For a single snow layer with increasing snow depths and densities, all three RTMs exhibit Tb attenuations arising from the increase in path lengths and optical depths. Further evaluations are conducted by feeding the three RTMs with the output of a snow physics model driven by in situ weather forcing in a coupled simulation. Outputs of this coupled model include snowpack physical properties and Tbs. By using snow stratigraphy observations, another set of Tb simulation is also conducted with RTMs driven by in situ snowpit observations. The snow physics outputs from the coupled case are compared against in situ snow stratigraphy observations.

Do Hyuk Kang↗

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.↗

Spectral reflectance of thin snow

A radiative transfer model was used to calculate the spectral reflectance of thin snow overlying dark soil. Model results show that directional hemispherical reflectance depends on density, grain size, and solar and viewing geometries. Measurements of thin snow spectral Bidirectional Reflectance Distribution Function (BRDF) show that the reflectance in the visible wavelengths is reduced when the sensor is near nadir yet, at a viewing zenith of 75 degrees the same snowpack will appear to be optically thick. If the snow is sufficiently thin and snow grains are large, the spectral signature of the substrate may influence snow reflectance. This phenomenon was also detected in reflectance data collected using the ground based PIDAS (Portable Instantaneous Display and Analysis Spectrometer). Concurrent airborne data from the AVIRIS instrument (Advanced Visible/Infrared Imaging Spectrometer) show effects of thin snow and mixed pixels for areas of thin and patchy spring snow.

Nolin, Anne W.↗

Top-of-Atmosphere Shortwave Broadband Observed Radiance and Estimated Irradiance over Polar Regions from Clouds and the Earth's Radiant Energy System (CERES) Instruments on Terra

Empirical angular distribution models for estimating top-of-atmosphere shortwave irradiances from radiance measurements over permanent snow, fresh snow and sea ice are developed using CERES measurements on Terra. Permanent snow angular distribution models depend on cloud fraction, cloud optical thickness, and snow brightness. Fresh snow and sea ice angular distribution models depend on snow and sea ice fraction, cloud fraction, cloud optical thickness, and snow and ice brightness. These classifications lead to 10 scene types for permanent snow and 25 scene types for fresh snow and sea ice. The average radiance over clear-sky permanent snow is more isotropic with satellite viewing geometry than that over overcast permanent snow. On average, the albedo of clear-sky permanent snow varies from 0.65 to 0.68 for solar zenith angles between 60$logical and\circ$ and 80 deg, while the corresponding albedo of overcast scenes varies from 0.70 to 0.73. Clear-sky permanent snow albedos over Antarctica estimated from two independent angular distribution models are consistent to within 0.6%, on average. Despite significant variability in sea ice optical properties with season, the estimated mean relative albedo error is -1 % for very dark sea ice and 0.1% for very bright sea ice when albedos derived from different viewing angles are averaged. The estimated regional root-mean-square (RMS) relative albedo error is 5.6% and 2.6% when the sea ice angular distribution models are applied to a region that contains very dark and very bright sea ice, respectively. Similarly, the estimated relative albedo bias error for fresh snow is -0.1% for very dark snow.

Kato, S.↗

Diagnosing coupled jet-streak circulations for a northern plains snow band from the operational nested-grid model

On 17 March 1989, moderate to heavy snow developed in a 100- to 200-km-wide band extending from South Dakota to northern Michigan. The 4- to 8-inch snowfall within this band was not associated with major cyclogenesis, and developed 500 to 600 km north of a stationary surface front. A diagnostic analysis based on an application of the General Meteorological Package (GEMPAK 5.0) to a numerical simulation from the operational nested-grid model (NGM) is utilized to show that the development of this snow band is related to the interaction of two upper-tropospheric jet streaks and their associated transverse circulation patterns. The eastward propagation of a jet streak from the West Coast toward the middle United States and to the south of a slower-moving jet along the U.S.-Canadian border led to a merger of the ascent maxima associated with the direct and indirect circulations of the northern and southern jets, respectively. The snow band developed as the ascending branches of the jet-streak circulation patterns merged, with the eastward propagation of the heaviest snow linked to the motion of the coupled circulation pattern. The study also demonstrates the usefulness of the operational NGM for providing the higher-resolution datasets required to relate the evolution of jet-streak circulation patterns to the development of mesoscale precipitation bands.

Hakim, Gregory J.↗

Assimilation of Multi-Frequency, Multi-Polarization Passive Microwave Brightness Temperature Observations in North America over Snow-Covered Regions Using Support Vector Machines

Accurately estimating the mass of water within a snowpack (a.k.a. snow water equivalent, or SWE) across regional or continental scales is a challenge. In order to overcome some of the limitations in traditional SWE retrieval algorithms or radiative transfer-based snow emission models, this study explores the use of a support vector machine (SVM) to merge an advanced land surface model within a radiance emission (i.e., brightness temperature) assimilation framework. The goal of direct radiance assimilation is preferable as it avoids inconsistencies in the use of ancillary data between the assimilation system and the independently-generated geophysical retrieval. The impact of assimilating multiple observations simultaneously at different frequency and polarization combinations is then evaluated via comparisons to state-of-the-art SWE and snow depth products as well as available ground-based measurements across North America for the years 2002 through 2011. It is found that assimilation-derived estimates (relative to estimates without assimilation) tend to better agree with state-of-the-art snow products. In addition, an overall improvement in goodness-of-fit statistics for snow estimates is achieved via assimilation when compared against ground-based snow measurements. In addition, these improvements in snow are shown to translate into improvements in streamflow predictions. Specifically, 11 out of the 13 major snow-dominated basins investigated have improved cumulative runoff estimates versus ground-based discharge measurements compared to the no-assimilation scenario. It is proven that a SVM can serve as an efficient and effective observation operator for a snow mass analysis within a radiance assimilation system.

Xue, Yuan↗

Arctic ice surface temperature retrieval from AVHRR thermal channels

The relationship between AVHRR thermal radiances and the surface (skin) temperature of Arctic snow-covered sea ice is examined through forward calculations of the radiative transfer equation, providing an ice/snow surface temperature retrieval algorithm for the central Arctic Basin. Temperature and humidity profiles with cloud observations collected on an ice island during 1986-1987 are used. Coefficients that correct for atmospheric attenuation are given for three Arctic clear sky 'seasons', as defined through statistical analysis of the daily profiles, for the NOAA 7, 9, and 11 satellites. Modeled directional snow emissivities, different in the two split-window (11 and 12 micron) channels, are used. While the sensor scan angle is included explicitly in the correction equation, its effect in the dry Arctic atmosphere is small, generally less than 0.1 K. Using the split-window channels and scan angle, the rms error in the estimated ice surface temperature is less than 0.1 K in all seasons. Inclusion channel 3(3.7 microns) during the winter decreases the rms error by less than 0.003

Key, J.↗

Comparison of SeaWinds scatterometer data with a hydrologic process model for the assessment of snow melt dynamics

The use of remote sensing observations for hydrological purposes is of particular interest in high latitude environments where in situ observations are sparse. Spacebornescatterometer data have the potential for monitoring freeze/thaw transitions and associated processes, and have various attractive attributes including frequent overpasses and all-weather capability.

seawinds Quikscat hydrologic process model↗

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↗

Code Description for "Brief Communication: Monitoring snow depth using small, cheap, and easy-to-deploy ground surface temperature sensors"

Temporally continuous snow depth estimates are vital for understanding changing snow patterns and impacts on permafrost in the Arctic. We train a random forest machine learning model to predict snow depth from variability in ground surface temperature. To our knowledge, this is the first time that small ground surface temperature sensors have been used to estimate snow depth. The model performs well at sites where the model was trained and at pan-arctic evaluation sites (RMSE <= 0.15 m). Small temperature sensors are cheap and easy-to-deploy, so this technique enables spatially distributed and temporally continuous snowpack monitoring to an extent previously infeasible. The model is flexible and can be applied to datasets retroactively to retrieve snow depth estimates at additional sites. This code package includes a *.joblib file of the trained random forest model and a *.ipynb file showing how to clean input data, train the random forest model, and apply the model.

Bachand, Claire↗

Diagnosis of Antarctic Blowing Snow Properties Using MERRA-2 Reanalysis with a Machine Learning Model

This paper presents the work on using a machine learning model to diagnose Antarctic blowing snow (BLSN) properties with the Modern Era Retrospective analysis for Research and Applications v2 (MERRA-2) data. We adopt the random forest classifier for BLSN identification and the random forest regressor for BLSN optical depth and height diagnosis. BLSN properties observed from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) are used as the truth for training the model. Using MERRA-2 fields such as snow age, surface elevation and pressure, temperature, specific humidity, and temperature gradient at the 2m level, and wind speed at the 10m level as input, reasonable results are achieved. Hourly blowing snow property diagnostics are generated with the trained model. Using the year 2010 as an example, it is shown that the Antarctic BLSN frequency is much higher over East than West Antarctica. High frequency months are from April to September, during which BLSN frequency exceeds 20% over East Antarctica. For May 2010, the BLSN snow frequency in the region is as high as 37%. Due to the suppression by strong surface-based inversions, larger values of BLSN height and optical depth are usually limited to the coastal regions, wherein the strength of surface-based inversions is weaker.

Antarctic↗