SEARCH · Engineering Papers

Results for “snow contribution”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Integrated Instrument Simulator Suites for Earth Science

The NASA Earth Observing System Simulators Suite (NEOS3) is a modular framework of forward simulations tools for remote sensing of Earth's Atmosphere from space. It was initiated as the Instrument Simulator Suite for Atmospheric Remote Sensing (ISSARS) under the NASA Advanced Information Systems Technology (AIST) program of the Earth Science Technology Office (ESTO) to enable science users to perform simulations based on advanced atmospheric and simple land surface models, and to rapidly integrate in a broad framework any experimental or innovative tools that they may have developed in this context. The name was changed to NEOS3 when the project was expanded to include more advanced modeling tools for the surface contributions, accounting for scattering and emission properties of layered surface (e.g., soil moisture, vegetation, snow and ice, subsurface layers). NEOS3 relies on a web-based graphic user interface, and a three-stage processing strategy to generate simulated measurements. The user has full control over a wide range of customizations both in terms of a priori assumptions and in terms of specific solvers or models used to calculate the measured signals.This presentation will demonstrate the general architecture, the configuration procedures and illustrate some sample products and the fundamental interface requirements for modules candidate for integration.

emission↗

NASA NTRS2012

Integrated Instrument Simulator Suites for Earth Science

radiative transfer↗

NASA NTRS1986

Variations in brightness temperature over cold first-year sea ice near Tuktoyaktuk, Northwest Territories

Microwave radiometric temperature T(B) profiles of first-year sea ice were obtained along 70- to 100-m traverses, with sled-mounted radiometers at 10, 18.7, 33.6, and 37 GHz and an effective spot size of 30 cm. Measurements of T(B) as a function of nadir angle were obtained at selected sites along the traverses. Snow and ice properties were recorded and correlated with the T(B) measurements to infer the effect of snow cover and ice conditions on the radiometric temperature. T(B) correlated positively with the brine volume profile in the ice at several sites, suggesting that brine volume has a strong effect on T(B) under these conditions. An overall statistical comparison of snow thickness with T(B), when compared with previously published models, suggests that the effect of snow cover on the microwave transmission coefficient of the snow/ice interface may be an important contribution to the radiometric temperature at these frequencies. A model is proposed to explain the data.

Lohanick, A. W.↗

DOE OSTI2024

Stream water sourcing from high-elevation snowpack inferred from stable isotopes of water: a novel application of d-excess values

About 80 % of the precipitation at the Colorado River's headwaters is snow, and the resulting snowmelt-driven hydrograph is a crucial water source for about 40 million people. Snowmelt from alpine and subalpine snowpack contributes substantially to groundwater recharge and river flow. However, the dynamics of snowmelt progression are not well understood because observations of the high-elevation snowpack are difficult due to challenging access in complex mountainous terrain as well as the cost and labor intensity of currently available methods. We present a novel approach to infer the processes and dynamics of high-elevation snowmelt contributions predicated upon stable hydrogen and oxygen isotope ratios observed in streamflow. We show that deuterium-excess (d-excess) values of stream water could serve as a comparatively cost-effective proxy for a catchment-integrated signal of high-elevation snowmelt contributions to catchment runoff. We sampled stable hydrogen and oxygen isotope ratios of the precipitation, snowpack, and stream water in the East River, a headwater catchment of the Colorado River, and the stream water of larger catchments at sites on the Gunnison River and Colorado River. The d-excess of snowpack increased with elevation; the upper subalpine and alpine snowpack (> 3200 m) had substantially higher d-excess compared to lower elevations (< 3200 m) in the study area. The d-excess values of stream water reflected this because d-excess values increased as the higher-elevation snowpack contributed more to stream water generation later in the snowmelt/runoff season. End-member mixing analyses based on the d-excess data showed that the share of high-elevation snowmelt contributions within the snowmelt hydrograph was on average 44 % and generally increased during melt period progression, up to 70 %. The observed pattern was consistent during 6 years for the East River, and a similar relation was found for the larger catchments on the Gunnison and Colorado rivers. High-elevation snowpack contributions were found to be higher for years with lower snowpack and warmer spring temperatures. Thus, we conclude that the d-excess of stream water is a viable proxy to observe changes in high-elevation snowmelt contributions in catchments at various scales. Inter-catchment comparisons and temporal trends of the d-excess of stream water could therefore serve as a catchment-integrated measure to monitor if mountain systems rely on high-elevation water inputs more during snow drought compared to years of average snowpack depths.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Variability in observed stable water isotopes in snowpack across a mountainous watershed in Colorado

In this study, isotopic information from 81 snowpits was collected over a 5-year period in a large, Colorado watershed. Data spans gradients in elevation, aspect, vegetation, and seasonal climate. They are combined with overlapping campaigns for water isotopes in precipitation and snowmelt, and a land-surface model for detailed estimates of snowfall and climate at sample locations. Snowfall isotopic inputs, describe the majority of δ18O snowpack variability. Aspect is a secondary control, with slightly more enriched conditions on east and north facing slopes. This is attributed to preservation of seasonally enriched snowfall and vapour loss in the early winter. Sublimation, expressed by decreases in snowpack d-excess in comparison to snowfall contributions, increases at low elevation and when seasonal temperature and solar radiation are high. At peak snow accumulation, post-depositional fractionation appears to occur in the top 25 ± 14% of the snowpack due to melt-freeze redistribution of lighter isotopes deeper into the snowpack and vapour loss to the atmosphere during intermittent periods of low relative humidity and high windspeed. Relative depth of fractionation increases when winter daytime temperatures are high and winter precipitation is low. Once isothermal, snowpack isotopic homogenization and enrichment was observed with initial snowmelt isotopically depleted in comparison to snowpack and enriching over time. The rate of δ18O increase (d-excess decrease) in snowmelt was 0.02‰ per day per 100-m elevation loss. Isotopic data suggests elevation dictates snowpack and snowmelt evolution by controlling early snow persistence (or absence), isotopic lapse rates in precipitation and the ratio of energy to snow availability. Hydrologic tracer studies using stable water isotopes in basins of large topographic relief will require adjustment for these elevational controls to properly constrain stream water sourcing from snowmelt.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Sensitivity of Simulated Polar Climate to Improved Partitioning of Spectral Solar Irradiance Between Visible and Near‐Infrared Bands

The solar radiative processes that contribute to Earth's surface and atmospheric energy budgets are strongly dependent on wavelength. For example, snow and water vapor become more absorptive as incident radiation shifts from visible (VIS) to near-infrared (NIR) wavelengths. Some earth system models (ESMs) aggregate solar radiation into just two bands (VIS and NIR) partitioned at 0.7 μm for transmission between the atmospheric and surface components. In the widely used radiative transfer model RRTMG_SW this partition is near the center of the overlap spectral band spanning 0.625–0.778 μm, whose flux is often approximated as being evenly divided between the VIS and NIR surface bands. Using a hyperspectral radiative transfer model, we show that the fractional downwelling surface flux within the overlap band is usually distributed about 55.5%:44.5% VIS:NIR. This improved approximation shifts as much as 4.27 W m −2 from the NIR to the VIS band, leading to an instantaneous decrease in surface absorption of up to 0.91 W m −2 over snow-covered surfaces. Century-long fully coupled ESM simulations show surface absorption over snow decreases by over 1.6 W m −2 . The coupled response in sea ice regions amplifies the initial forcing ten-fold, and increases seasonal sea ice area by up to 4.9%. These results highlight the importance of accurately representing the spectral distribution of solar radiation in the cryosphere.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Deep Learning Advances Arctic River Water Temperature Predictions

The accelerated warming in the Arctic poses serious risks to freshwater ecosystems by altering streamflow and river thermal regimes. However, limited research on Arctic River water temperatures exists due to data scarcity and the absence of robust methodologies, which often focus on large, major river basins. To address this, we leveraged the newly released, extensive AKTEMP data set and advanced machine learning techniques to develop a Long Short-Term Memory (LSTM) model. By incorporating ERA5-Land reanalysis data and integrating physical understanding into data-driven processes, our model advanced river water temperature predictions in ungauged, snow- and permafrost-affected basins in Alaska. Our model outperformed existing approaches in high-latitude regions, achieving a median Nash-Sutcliffe Efficiency of 0.95 and root mean squared error of 1.0°C. The LSTM model learned air temperature, soil temperature, solar radiation, and thermal radiation—factors associated with energy balance—were the most important drivers of river temperature dynamics. Soil moisture and snow water equivalent were highlighted as critical factors representing key processes such as thawing, melting, and groundwater contributions. Glaciers and permafrost were also identified as important covariates, particularly in seasonal river water temperature predictions. Our LSTM model successfully captured the complex relationships between hydrometeorological factors and river water temperatures across varying timescales and hydrological conditions. This scalable and transferable approach can be potentially applied across the Arctic, offering valuable insights for future conservation and management efforts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Measurement report: Molecular composition, optical properties, and radiative effects of water-soluble organic carbon in snowpack samples from northern Xinjiang, China

Abstract. Water-soluble organic carbon (WSOC) in the cryosphere has an important impact on the biogeochemistry cycling and snow–ice surface energy balance through changes in the surface albedo. This work reports on the chemical characterization of WSOC in 28 representative snowpack samples collected across a regional area of northern Xinjiang, northwestern China. We employed multimodal analytical chemistry techniques to investigate both bulk and molecular-level composition of WSOC and its optical properties, informing the follow-up radiative forcing (RF) modeling estimates. Based on the geographic differences and proximity of emission sources, the snowpack collection sites were grouped as urban/industrial (U), rural/remote (R), and soil-influenced (S) sites, for which average WSOC total mass loadings were measured as 1968 ± 953 ng g−1 (U), 885 ± 328 ng g−1 (R), and 2082 ± 1438 ng g−1 (S), respectively. The S sites showed the higher mass absorption coefficients at 365 nm (MAC365) of 0.94 ± 0.31 m2 g−1 compared to those of U and R sites (0.39 ± 0.11 m2 g−1 and 0.38 ± 0.12 m2 g−1, respectively). Bulk composition of WSOC in the snowpack samples and its basic source apportionment was inferred from the excitation–emission matrices and the parallel factor analysis featuring relative contributions of one protein-like (PRLIS) and two humic-like (HULIS-1 and HULIS-2) components with ratios specific to each of the S, U, and R sites. Additionally, a sample from site 120 showed unique pollutant concentrations and spectroscopic features remarkably different from all other U, R, and S samples. Molecular-level characterization of WSOC using high-resolution mass spectrometry (HRMS) provided further insights into chemical differences among four types of samples (U, R, S, and 120). Specifically, many reduced-sulfur-containing species with high degrees of unsaturation and aromaticity were uniquely identified in U samples, suggesting an anthropogenic source. Aliphatic/protein-like species showed the highest contribution in R samples, indicating their biogenic origin. The WSOC components from S samples showed high oxygenation and saturation levels. A few unique CHON and CHONS compounds with high unsaturation degree and molecular weight were detected in the 120 sample, which might be anthraquinone derivatives from plant debris. Modeling of the WSOC-induced RF values showed warming effects of 0.04 to 0.59 W m−2 among different groups of sites, which contribute up to 16 % of that caused by black carbon (BC), demonstrating the important influences of WSOC on the snow energy budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Capturing Sub‐Kilometer Flood Inundation Dynamics During the California Rain‐on‐Snow Events of 2017

The severe impacts of rain‐on‐snow (ROS) extreme events have been widely recognized and studied. However, unlike hydrological processes, flood inundation dynamics and the relative contribution of rainfall and snowmelt during ROS events remain under‐investigated. We diagnosed and documented sub‐kilometer spatio‐temporal dynamics of flood inundation during the 2017 California ROS events, simulated by a 2‐dimensional hydrodynamic model, River Dynamical Core (RDycore). RDycore shows good performance in capturing fine‐scale flood inundation dynamics against gauge measurements (with a median correlation coefficient of 0.81) and satellite observations. On top of rainfall, snowmelt not only increases the mean maximum inundation depth (12.0%–25.1%), but also expands the total flooded area (19.9%–31.9%) and prolongs the mean flood duration (3.4%–7.1%) across the events. Our findings offer an explicit and accurate picture of when and where ROS flooding could occur and how snowmelt increases flood hazard, valuable for risk assessment and infrastructure planning.

Flood inundation↗

NASA NTRS1995

C-Band Backscatter Measurements of Winter Sea-Ice in the Weddell Sea, Antarctica

During the 1992 Winter Weddell Gyre Study, a C-band scatterometer was used from the German ice-breaker R/V Polarstern to obtain detailed shipborne measurement scans of Antarctic sea-ice. The frequency-modulated continuous-wave (FM-CW) radar operated at 4-3 GHz and acquired like- (VV) and cross polarization (HV) data at a variety of incidence angles (10-75 deg). Calibrated backscatter data were recorded for several ice types as the icebreaker crossed the Weddell Sea and detailed measurements were made of corresponding snow and sea-ice characteristics at each measurement site, together with meteorological information, radiation budget and oceanographic data. The primary scattering contributions under cold winter conditions arise from the air/snow and snow/ice interfaces. Observations indicate so e similarities with Arctic sea-ice scattering signatures, although the main difference is generally lower mean backscattering coefficients in the Weddell Sea. This is due to the younger mean ice age and thickness, and correspondingly higher mean salinities. In particular, smooth white ice found in 1992 in divergent areas within the Weddell Gyre ice pack was generally extremely smooth and undeformed. Comparisons of field scatterometer data with calibrated 20-26 deg incidence ERS-1 radar image data show close correspondence, and indicate that rough Antarctic first-year and older second-year ice forms do not produce as distinctively different scattering signatures as observed in the Arctic. Thick deformed first-year and second-year ice on the other hand are clearly discriminated from younger undeformed ice. thereby allowing successful separation of thick and thin ice. Time-series data also indicate that C-band is sensitive to changes in snow and ice conditions resulting from atmospheric and oceanographic forcing and the local heat flux environment. Variations of several dB in 45 deg incidence backscatter occur in response to a combination of thermally-regulated parameters including sea-ice brine volume, snow and ice complex dielectric properties, and snow physical properties.

Drinkwater, M. R.↗

DOE OSTI2022

Contemporary sources dominate carbonaceous aerosol on the North Slope of Alaska

As the Arctic continues to change and warm rapidly, it is increasingly important to understand the organic carbon (OC) contribution to Arctic aerosol. Biogenic sources of primary and secondary OC in the Arctic will be impacted by climate change, including warming temperatures and earlier snow and ice melt. Here, this study focuses on identifying potential sources and regional influences on the seasonal concentration of organic aerosol through analysis of chemical and isotopic composition. Aerosol samples were collected at two sites on the North Slope of Alaska ($Utqia\dot{g}vik$, UQK, and Oliktok Point, OLK, which is in an Arctic oilfield) over three summers from 2015 to 2017. The elemental carbon (EC) trends at each site were used to understand local combustion influences. Local sources drove EC concentrations at Oliktok Point, where high EC was attributed to oil and gas extraction activity, including diesel combustion emissions. $Utqia\dot{g}vik$ had very low EC in the summer. OC was more similar in concentration and well correlated between the two sites with high contributions of contemporary carbon by radiocarbon apportionment (UQK = 74%, OLK = 63%), which could include both marine and terrestrial sources of contemporary carbon (e.g. primary and secondary biogenic, biomass burning and/or associated SOA, and bioaerosols). OC concentrations are strongly correlated to maximum ambient temperatures on the NSA during the summer, which may have implications for predicting future OC aerosol concentrations in a warming Arctic. Biomass burning was determined to be an episodic influence at both sites, based on interpretation of combined aerosol composition, air mass trajectories, and remote sensing of smoke plumes. The results from this study overall strongly suggests contribution from regional sources of contemporary organic aerosol on the NSA, but additional analysis is needed to better constrain contributions from both biogenic sources (terrestrial and/or marine) and bioaerosol to better understand temperature-related aerosol processes in the Arctic.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1988

Search for the Tunguska event in the Antarctic snow

The Tunguska explosion in 1908 is supposed to have been produced by the impact of a small celestial body. The absence of any identifiable crater together with the huge energy released by the event suggest that the impactor exploded in midair and that its material was widely spread over the Earth. The short term contribution of such exceptional events to the total accretion rate of extraterrestrial material by the Earth could be significant. Samples were chosen in a core electromechanically drilled in 1984 near South Pole Station. There, the low temperatures, preventing melting all year long, and the nearly regular snow fall rate provide good conditions for a reliable continuous record of any infalling material. In many samples Ir was below the detection limit of the instrumentation. The iridium infall averaged over 45 samples is given. In a few samples the iridium content is significantly higher than the average: the frequency and amplitude of such fluctuations can be explained by the presence on some filters of finite size cosmic particles. No significant systematic increase above the average level is observed in the part of the core corresponding to the Tunguska event. The two major results of this study are: (1) The presence of Tunguska explosion debris in the Antarctic snow is not confirmed; and (2) The estimate of the average iridium infall, is an order of magnitude lower than the Ganapathy's background but is close to the values measured in Antarctic snow and atmospheric samples by Takahashi et al. The results are also consistent with the flux of micrometeoroids deduced from optical and radar observations or derived from the study of Greenland cosmic dust collection but are lower than the flux at mid-latitude measured in paleocene-oligocene sediments from the central part of the Pacific Ocean.

Rocchia, R.↗

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

NASA NTRSDate not supplied

Applications of MERRA-2 data for avian migration, biomass burning, and dusty atmospheric rivers

Three different applications of MERRA-2 data are presented. 1) Using radar data, we introduced a new concept for spatial patterns of bird migration across the contiguous U.S. This approach allowed us to use MERRA-2 data and learn that remote forcing in the tropical Pacific—through a chain of processes including atmospheric Rossby wave trains— controls the climatic conditions, associated with bird migration in North America. 2) We showed that emissions from biomass burning in the Congo Basin are partly controlled by the low-level winds, which are in turn associated with the intensity of the subtropical high in the Indian Ocean. Using back-trajectory analysis, we found that these emissions combined with their transport mechanism explain the interannual variability of black carbon in West Africa. 3) Our analysis showed that atmospheric rivers in the Middle East contribute to both heavy flood and dust transport within their corridor. We also found that warm advection and rain-on-snow effect of dusty atmospheric rivers further enhance the chance of flood through rapid snowmelt processes.

Amin Dezfuli↗

NASA NTRSDate not supplied

Applications of MERRA-2 Data for Avian Migration, Biomass Burning, and Dusty Atmospheric Rivers

Amin Dezfuli↗

NASA NTRS2012

Iterative Atmospheric Correction Scheme and the Polarization Color of Alpine Snow

Characterization of the Earth's surface is crucial to remote sensing, both to map geomorphological features and because subtracting this signal is essential during retrievals of the atmospheric constituents located between the surface and the sensor. Current operational algorithms model the surface total reflectance through a weighted linear combination of a few geometry-dependent kernels, each devised to describe a particular scattering mechanism. The information content of these measurements is overwhelmed by that of instruments with polarization capabilities: proposed models in this case are based on the Fresnel reflectance of an isotropic distribution of facets. Because of its remarkable lack of spectral contrast, the polarized reflectance of land surfaces in the shortwave infrared spectral region, where atmospheric scattering is minimal, can be used to model the surface also at shorter wavelengths, where aerosol retrievals are attempted based on well-established scattering theories. In radiative transfer simulations, straightforward separation of the surface and atmospheric contributions is not possible without approximations because of the coupling introduced by multiple reflections. Within a general inversion framework, the problem can be eliminated by linearizing the radiative transfer calculation, and making the Jacobian (i.e., the derivative expressing the sensitivity of the reflectance with respect to model parameters) available at output. We present a general methodology based on a Gauss-Newton iterative search, which automates this procedure and eliminates de facto the need of an ad hoc atmospheric correction. In this case study we analyze the color variations in the polarized reflectance measured by the NASA Goddard Institute of Space Studies Research Scanning Polarimeter during a survey of late-season snowfields in the High Sierra. This insofar unique dataset presents challenges linked to the rugged topography associated with the alpine environment and a likely high water content due to melting. The analysis benefits from ancillary information provided by the NASA Langley High Spectral Resolution Lidar deployed on the same aircraft. The results obtained from the iterative scheme are contrasted against the surface polarized reflectance obtained ignoring multiple reflections, via the simplistic subtraction of the atmospheric scattering contribution. Finally, the retrieved reflectance is modeled after the scattering properties of a dense collection of ice crystals at the surface. Confirming that the polarized reflectance of snow is spectrally flat would allow to extend the techniques already in use for polarimetric retrievals of aerosol properties over land to the large portion of snow-covered pixels plaguing orbital and suborbital observations.

Polarized BRDF↗

DOE OSTI2020

Dust dominates high-altitude snow darkening and melt over high-mountain Asia

Westerly-driven, long-range transportation of dust particles in elevated aerosol layers (EALs) is a persistent phenomenon during spring and summer over the Indian subcontinent. During the snow accumulation season, EALs transport substantial dust to the snow-covered slopes of high-mountain Asia (HMA). Here we use multiple satellite-based estimates to demonstrate a robust physical association between the EALs and dust-induced snow darkening over HMA. Results from a fully-coupled atmosphere-chemistry-snow model support these observations, revealing a signature of increasing dust-induced snow darkening with surface elevation over HMA that peaks near 4500 m. Moreover, the influence of dust on snow darkening is greater than that of black carbon above 4000 m. Our findings suggest a discernable role of dust in the observed spatial heterogeneity of snowmelt and snowline trends over HMA and highlight an increasing contribution of dust to snowmelt as the snowline rises with warming.

Elevated aerosol layers (EAL), Light absorbing par↗

NASA NTRS2021

Physics-Based SNOWPACK Model Improves Representation of Near-Surface Antarctic Snow and Firn Density

Estimates of snow and firn density are required for satellite-altimetry-based retrievals of ice sheet mass balance that rely on volume-to-mass conversions. Therefore, biases and errors in presently used density models confound assessments of ice sheet mass balance and by extension ice sheet contribution to sea level rise. Despite this importance, most contemporary firn densification models rely on simplified semi-empirical methods, which are partially reflected by significant modeled density errors when compared to observations. In this study, we present a new drifting-snow compaction scheme that we have implemented into SNOWPACK, a physics-based land surface snow model. We show that our new scheme improves existing versions of SNOWPACK by increasing simulated near-surface (defined as the top 10 m) density to be more in line with observations (near-surface bias reduction from −44.9 to −5.4 kg m−3). Furthermore, we demonstrate high-quality simulation of near-surface Antarctic snow and firn density at 122 observed density profiles across the Antarctic ice sheet, as indicated by reduced model biases throughout most of the near-surface firn column when compared to two semi-empirical firn densification models (SNOWPACK mean bias=−9.7 kg m−3, IMAU-FDM mean bias=−32.5 kg m−3, GSFC-FDM mean bias=15.5 kg m−3). Notably, our analysis is restricted to the near surface where firn density is most variable due to accumulation and compaction variability driven by synoptic weather and seasonal climate variability. Additionally, the GSFC-FDM exhibits lower mean density bias from 7–10 m (SNOWPACK bias=−22.5 kg m−3, GSFC-FDM bias=10.6 kg m−3) and throughout the entire near surface at high-accumulation sites (SNOWPACK bias=−31.4 kg m−3, GSFC-FDM bias=−4.7 kg m−3). However, we found that the performance of SNOWPACK did not degrade when applied to sites that were not included in the calibration of semi-empirical models. This suggests that SNOWPACK may possibly better represent firn properties in locations without extensive observations and under future climate scenarios, when firn properties are expected to diverge from their present state.

Eric Keenan↗