Spatial Distribution Of Snow Research Articles

Three Centuries of Snowpack Decline at an Alpine Pass Revealed by Cosmogenic Paleothermometry and Luminescence Photochronometry

AbstractThe spatial and temporal distribution of Alpine snow is a sensitive gauge of environmental change. While understanding past snow dynamics is essential for reconstructing past climate and forecasting future trends, reliable snowpack data prior to the instrumental record are scarce. We present a novel pairing of cosmogenic paleothermometry and luminescence photochronometry which constrain the temperature and insolation history of bedrock outcrops at the Gotthard Pass, Switzerland, over the last ∼15,000 years. By coupling these results with cosmogenic 14C‐10Be chronology and modern in situ rock thermometry, we infer a ∼70‐day reduction of snowpack at the topographic mid‐slope. Our data indicate stable environmental conditions throughout the Holocene, followed by a 6.6 ± 2.9°C increase of ground surface temperature, coeval with an order‐of‐magnitude or more increase in ground surface insolation. Bracketing the onset of these changes between 1504 and 1807 CE, our findings tie the snowpack decline with the onset of human industrialization.

Regional Impact of Snow‐Darkening on Snow Pack and the Atmosphere During a Severe Saharan Dust Deposition Event in Eurasia

AbstractLight‐absorbing impurities such as mineral dust can play a major role in reducing the albedo of snow surfaces. Particularly in spring, deposited dust particles lead to increased snow melt and trigger further feedbacks at the land surface and in the atmosphere. Quantifying the extent of dust‐induced variations is difficult due to high variability in the spatial distribution of mineral dust and snow. We present an extension of a fully coupled atmospheric and land surface model system to address the impact of mineral dust on the snow albedo across Eurasia. We evaluated the short‐term effects of Saharan dust in a case study. To obtain robust results, we performed an ensemble simulation followed by statistical analysis. Mountainous regions showed a strong impact of dust deposition on snow depth. We found a mean significant reduction of −1.4 cm in the Caucasus Mountains after 1 week. However, areas with flat terrain near the snow line also showed strong effects despite lower dust concentrations. Here, the feedback to dust deposition was more pronounced as increase in surface temperature and air temperature. In the region surrounding the snow line, we found an average significant surface warming of 0.9 K after 1 week. This study shows that the impact of mineral dust deposition depends on several factors. Primarily, these are altitude, slope, snow depth, and snow cover fraction. Especially in complex terrain, it is therefore necessary to use fully coupled models to investigate the effects of mineral dust on snow pack and the atmosphere.

Snowfall Variation in Eastern Mediterranean Catchments

This study aims to present and analyze the time series of the snow parameters focusing on representative geographical areas of the Eastern Mediterranean (i.e., Greece and Italy) and to examine their seasonal variability, in terms of region and geography. The satellite retrievals were firstly validated against in-situ retrievals for 67 common days, with a mean bias equal to −0.018 cm, with a near-Gaussian distribution, showing the good performance of the satellite snow detection. The satellite-based analysis resulted in increasing trends of snow water equivalent, attributed to the enhanced values between 2000 and 2009; however, decreasing trends are found starting from 2010 until now of −1.79 × 10−17 and −2.31 × 10−18 over the two representative areas of Greece (e.g., Thessaloniki and Kozani). A similar pattern is found for the snow water equivalent in the Italian study area, with a decreasing trend of −4.45 × 10−18. The presented results contribute to a better understanding of the spatial snow distribution and the snow coverage seasonality that could be crucial for the long-term groundwater management, by combining snow data trends from in-situ data and satellite statistics.

Impact of distributed meteorological forcing on simulated snow cover and hydrological fluxes over a mid-elevation alpine micro-scale catchment

Abstract. From the micro- to the mesoscale, water and energy budgets of mountainous catchments are largely driven by topographic features such as terrain orientation, slope, steepness, and elevation, together with associated meteorological forcings such as precipitation, solar radiation, and wind speed. Those topographic features govern the snow deposition, melting, and transport, which further impacts the overall water cycle. However, this microscale variability is not well represented in Earth system models due to coarse resolutions. This study explores the impact of precipitation, shortwave radiation, and wind speed on the water budget distribution over a 15.28 ha small, mid-elevation (2000–2200 m) alpine catchment at Col du Lautaret (France). The grass-dominated catchment remains covered with snow for 5 to 6 months per year. The surface–subsurface coupled distributed hydrological model ParFlow-CLM is used at a very high resolution (10 m) to simulate the impacts on the water cycle of meteorological variability at very small spatial and temporal scales. These include 3D simulations of hydrological fluxes with spatially distributed forcing of precipitation, shortwave radiation, and wind speed compared to 3D simulations of hydrological fluxes with non-distributed forcing. Our precipitation distribution method encapsulates the spatial snow distribution along with snow transport. The model simulates the dynamics and spatial variability of snow cover using the Common Land Model (CLM) energy balance module and under different combinations of distributed forcing. The resulting subsurface and surface water transfers are computed by the ParFlow module. Distributed forcing leads to spatially heterogeneous snow cover simulation, which becomes patchy at the end of the melt season and shows a good agreement with the remote sensing images (mean bias error (MBE) = 0.22). This asynchronous melting results in a longer melting period compared to the non-distributed forcing, which does not generate any patchiness. Among the distributed meteorological forcings tested, precipitation distribution, including snow transport, has the greatest impact on spatial snow cover (MBE = 0.06) and runoff. Shortwave radiation distribution has an important impact, reducing evapotranspiration as a function of the slope orientation (decreasing the slope between observed and simulated evapotranspiration from 1.55 to 1.18). For the primarily east-facing catchment studied here, distributing shortwave radiation helps generate realistic timing and spatial heterogeneity in the snowmelt at the expense of an increase in the mean bias error (from 0.06 to 0.22) for all distributed forcing simulations compared to the simulation with only distributed precipitation. Distributing wind speed in the energy balance calculation has a more complex impact on our catchment, as it accelerates snowmelt when meteorological conditions are favorable but does not generate snow patches at the end of our test case. This shows that slope- and aspect-based meteorological distribution can improve the spatio-temporal representation of snow cover and evapotranspiration in complex mountain terrain.

Emulating the Adaptation of Wind Fields to Complex Terrain with Deep Learning

Abstract Estimating the impact of wind-driven snow transport requires modeling wind fields with a lower grid spacing than the spacing on the order of 1 or a few kilometers used in the current numerical weather prediction (NWP) systems. In this context, we introduce a new strategy to downscale wind fields from NWP systems to decametric scales, using high-resolution (30 m) topographic information. Our method (named “DEVINE”) is leveraged on a convolutional neural network (CNN), trained to replicate the behavior of the complex atmospheric model ARPS, and was previously run on a large number (7279) of synthetic Gaussian topographies under controlled weather conditions. A 10-fold cross validation reveals that our CNN is able to accurately emulate the behavior of ARPS (mean absolute error for wind speed = 0.16 m s−1). We then apply DEVINE to real cases in the Alps, that is, downscaling wind fields forecast by the AROME NWP system using information from real alpine topographies. DEVINE proved able to reproduce main features of wind fields in complex terrain (acceleration on ridges, leeward deceleration, and deviations around obstacles). Furthermore, an evaluation on quality-checked observations acquired at 61 sites in the French Alps reveals improved behavior of the downscaled winds (AROME wind speed mean bias is reduced by 27% with DEVINE), especially at the most elevated and exposed stations. Wind direction is, however, only slightly modified. Hence, despite some current limitations inherited from the ARPS simulations setup, DEVINE appears to be an efficient downscaling tool whose minimalist architecture, low input data requirements (NWP wind fields and high-resolution topography), and competitive computing times may be attractive for operational applications. Significance Statement Wind largely influences the spatial distribution of snow in mountains, with direct consequences on hydrology and avalanche hazard. Most operational models predicting wind in complex terrain use a grid spacing on the order of several kilometers, too coarse to represent the real patterns of mountain winds. We introduce a novel method based on deep learning to increase this spatial resolution while maintaining acceptable computational costs. Our method mimics the behavior of a complex model that is able to represent part of the complexity of mountain winds by using topographic information only. We compared our results with observations collected in complex terrain and showed that our model improves the representation of winds, notably at the most elevated and exposed observation stations.

Spatial patterns of snow distribution in the sub-Arctic

Abstract. The spatial distribution of snow plays a vital role in sub-Arctic and Arctic climate, hydrology, and ecology due to its fundamental influence on the water balance, thermal regimes, vegetation, and carbon flux. However, the spatial distribution of snow is not well understood, and therefore, it is not well modeled, which can lead to substantial uncertainties in snow cover representations. To capture key hydro-ecological controls on snow spatial distribution, we carried out intensive field studies over multiple years for two small (2017–2019; ∼ 2.5 km2) sub-Arctic study sites located on the Seward Peninsula of Alaska. Using an intensive suite of field observations (> 22 000 data points), we developed simple models of the spatial distribution of snow water equivalent (SWE) using factors such as topographic characteristics, vegetation characteristics based on greenness (normalized different vegetation index, NDVI), and a simple metric for approximating winds. The most successful model was random forest, using both study sites and all years, which was able to accurately capture the complexity and variability of snow characteristics across the sites. Approximately 86 % of the SWE distribution could be accounted for, on average, by the random forest model at the study sites. Factors that impacted year-to-year snow distribution included NDVI, elevation, and a metric to represent coarse microtopography (topographic position index, TPI), while slope, wind, and fine microtopography factors were less important. The characterization of the SWE spatial distribution patterns will be used to validate and improve snow distribution modeling in the Department of Energy's Earth system model and for improved understanding of hydrology, topography, and vegetation dynamics in the sub-Arctic and Arctic regions of the globe.

Development and parameter estimation of snowmelt models using spatial snow-cover observations from MODIS

Abstract. Given the importance of snow on different land and atmospheric processes, accurate representation of seasonal snow evolution, including distribution and melt volume, is highly imperative to any water resources development trajectories. The limitation of reliable snowmelt estimation in mountainous regions is, however, further exacerbated by data scarcity. This study attempts to develop relatively simple extended degree-day snow models driven by freely available snow-cover images. This approach offers relative simplicity and a plausible alternative to data-intensive models, as well as in situ measurements, and has a wide range of applicability, allowing for immediate verification with point measurements. The methodology employs readily available MODIS composite images to calibrate the snowmelt models on spatial snow distribution in contrast to the traditional snow-water-equivalent-based calibration. The spatial distribution of snow-cover is simulated using different extended degree-day models with parameters calibrated against individual MODIS snow-cover images for cloud-free days or a set of images representing a period within the snow season. The study was carried out in Baden-Württemberg (Germany) and in Switzerland. The simulated snow-cover data show very good agreement with MODIS snow-cover distribution, and the calibrated parameters exhibit relative stability across the time domain. Furthermore, different thresholds that demarcate snow and no-snow pixels for both observed and simulated snow cover were analyzed to evaluate these thresholds' influence on the model performance and identified for the study regions. The melt data from these calibrated snow models were used as standalone inputs to a modified Hydrologiska Byråns Vattenbalansavdelning (HBV) without the snow component in all the study catchments to assess the performance of the melt outputs in comparison to a calibrated standard HBV model. The results show an overall increase in Nash–Sutcliffe efficiency (NSE) performance and a reduction in uncertainty in terms of model performance. This can be attributed to the reduction in the number of parameters available for calibration in the modified HBV and an added reliability of the snow accumulation and melt processes inherent in the MODIS calibrated snow model output. This paper highlights that the calibration using readily available images used in this method allows for a flexible regional calibration of snow-cover distribution in mountainous areas with reasonably accurate precipitation and temperature data and globally available inputs. Likewise, the study concludes that simpler specific alterations to processes contributing to snowmelt can contribute to reliably identify the snow distribution and bring about improvements in hydrological simulations, owing to better representation of the snow processes in snow-dominated regimes.

FlakeOut: A geometric approach to remove wind-blown snow from terrestrial laser scans

Wind-blown snow particles often contaminate Terrestrial Laser Scanning (TLS) data of snow covered terrain. However, common filtering techniques fail to filter wind-blown snow and incorrectly filter data from the true surface due to the spatial distribution of wind-blown snow and the TLS scanning geometry. We present FlakeOut, a filter designed specifically to filter wind-blown snowflakes from TLS data. A key aspect of FlakeOut is a low false positive rate of 2.8x10-4—an order of magnitude lower than standard filtering techniques—which greatly reduces the number of true ground points that are incorrectly removed. This low false positive rate makes FlakeOut appropriate for applications requiring quantitative measurements of the snow surface in light to moderate blowing snow conditions. Additionally, we provide mathematical and software tools to efficiently estimate the false positive rate of filters applied for the purpose of removing erroneous data points that occur very infrequently in a dataset.

DriftScalarDyFoam: An OpenFOAM-Based Multistage Solver for Drifting Snow and Its Distribution Around Buildings

There is a complex coupling relationship between the airflow and snow cover. In a period of hours or even days, the airflow will cause the redistribution of snow, and the redistribution of snow will cause the airflow to change. This study develops a dynamic mesh technology applied in snow drifting simulation through a real-time dynamic mesh update to depict the snow surface evolution process under long-period snow drifting, and a solver application named driftScalarDyFoam based on OpenFOAM is implemented. This solver divides the long-period snow drifting process into several stages, in each of which a snow transport equation is applied to predict the spatial distribution of snow, and finally, the snow surface evolves according to the erosion–deposition model. This method that we have proposed has been validated for several measured cases, including snow distribution on a flat roof and snow distribution around a building.

Wind-tunnel experimental studies of the spatial snow distribution over grass and bush surfaces

Wind-tunnel experimental studies of the spatial snow distribution over grass and bush surfaces

Statistical characteristics of the spatial distribution of wind and snow in the Xinjiang Uygur Autonomous Region

At present, the wind and snow loads are calculated independently when determining the design specifications for building structures. Yet, when snow drifting occurs, the basic wind and snow pressures exist almost simultaneously. Therefore, building specifications based on the independent calculations of the wind and snow loads cannot be used effectively in areas that are severely impacted by snow drifting. Some parts of the Xinjiang Uygur Autonomous Region of China have suffered particularly severe snow drifting natural disasters. In this paper, we offer an analysis of the statistical characteristics of the spatial distribution of wind and snow in this region. In particular, we extract the values of the combined distribution of the wind and snowfall in the Xinjiang Uygur Autonomous Region of China by using parameters such as snowfall intensity, snowfall amount, wind speed, and wind direction as indicators that show the snow drifting disasters. This study found that, after heavy snowfall events, the accumulative wind scale is high and the accumulative snowfall is large in the Altay, the Bortala Mongol Autonomous, and the Tacheng Prefectures of northern Xinjiang, and the Kizilsu Kirghiz Autonomous Prefecture of western Xinjiang. It has an important practical significance for the design of building structures and the treatment of snow drifting disasters.

The sensitivity of runoff generation to spatial snowpack uniformity in an alpine watershed: Green Lakes Valley, Niwot Ridge Long‐Term Ecological Research station

AbstractSeasonal water storage in high‐elevation alpine catchments are critical sources of water for mountainous regions like the western U.S. The spatial distribution of snow in these topographically complex catchments is primarily governed by orography, solar radiation, and wind redistribution. While the effect of solar shading is relatively consistent from year‐to‐year, the redistribution of snow due to wind is more variable – capable of producing snowpacks that have varying degrees of uniformity across these hydrologically‐important catchments. A reasonable hypothesis is that a warmer climate will cause snowfall to become more dense (i.e. wetter and heavier), possibly leading to less wind redistribution and thus produce a more uniformly distributed snowpack across the landscape. In this study, we investigate the role of increasingly uniform spatial snowpack distributions on streamflow generation in the Green Lakes Valley Niwot Ridge Long Term Ecological Research station, within the headwaters of the Boulder Creek watershed in Colorado. A set of idealized hydrologic simulation experiments driven by reconstructed snowpacks spanning 2001–2014 show that more a more uniform spatial snowpack distribution leads to an earlier melt‐out of 31 days on average and tends to produce less total streamflow, with maximum decreases as large as 7.5%. Isolating the role of snowpack heterogeneity from melt‐season precipitation, we find that snowpack uniformity reduces total streamflow by as much as 13.2%. Reductions in streamflow are largely explained by greater exposure to solar radiation in the uniformly distributed case relative to a more heterogeneous snowpack, with this exposure driving shifts towards earlier snowmelt and changes in soil water storage. Overall, we find that the runoff efficiency from shallower snowpacks is more sensitive to the effects of uniformity than deeper snowpacks, which has potential implications for a warming climate where shallower snowpacks and enhanced sensitivities may be present.

Observed snow depth trends in the European Alps: 1971 to 2019

Abstract. The European Alps stretch over a range of climate zones which affect the spatial distribution of snow. Previous analyses of station observations of snow were confined to regional analyses. Here, we present an Alpine-wide analysis of snow depth from six Alpine countries – Austria, France, Germany, Italy, Slovenia, and Switzerland – including altogether more than 2000 stations of which more than 800 were used for the trend assessment. Using a principal component analysis and k-means clustering, we identified five main modes of variability and five regions which match the climatic forcing zones: north and high Alpine, north-east, north-west, south-east, and south and high Alpine. Linear trends of monthly mean snow depth between 1971 and 2019 showed decreases in snow depth for most stations from November to May. The average trend among all stations for seasonal (November to May) mean snow depth was −8.4 % per decade, for seasonal maximum snow depth −5.6 % per decade, and for seasonal snow cover duration −5.6 % per decade. Stronger and more significant trends were observed for periods and elevations where the transition from snow to snow-free occurs, which is consistent with an enhanced albedo feedback. Additionally, regional trends differed substantially at the same elevation, which challenges the notion of generalizing results from one region to another or to the whole Alps. This study presents an analysis of station snow depth series with the most comprehensive spatial coverage in the European Alps to date.

Fractional snow-covered area: scale-independent peak of winter parameterization

Abstract. The spatial distribution of snow in the mountains is significantly influenced through interactions of topography with wind, precipitation, shortwave and longwave radiation, and avalanches that may relocate the accumulated snow. One of the most crucial model parameters for various applications such as weather forecasts, climate predictions and hydrological modeling is the fraction of the ground surface that is covered by snow, also called fractional snow-covered area (fSCA). While previous subgrid parameterizations for the spatial snow depth distribution and fSCA work well, performances were scale-dependent. Here, we were able to confirm a previously established empirical relationship of peak of winter parameterization for the standard deviation of snow depth σHS by evaluating it with 11 spatial snow depth data sets from 7 different geographic regions and snow climates with resolutions ranging from 0.1 to 3 m. An enhanced performance (mean percentage errors, MPE, decreased by 25 %) across all spatial scales ≥ 200 m was achieved by recalibrating and introducing a scale-dependency in the dominant scaling variables. Scale-dependent MPEs vary between −7 % and 3 % for σHS and between 0 % and 1 % for fSCA. We performed a scale- and region-dependent evaluation of the parameterizations to assess the potential performances with independent data sets. This evaluation revealed that for the majority of the regions, the MPEs mostly lie between ±10 % for σHS and between −1 % and 1.5 % for fSCA. This suggests that the new parameterizations perform similarly well in most geographical regions.

Interannual and Seasonal Variability of Snow Depth Scaling Behavior in a Subalpine Catchment

AbstractUnderstanding and characterizing the spatial distribution of snow are critical to represent the energy balance and runoff production in mountain environments. In this study, we investigate the interannual and seasonal variability in snow depth scaling behavior at the Izas experimental catchment of the Spanish Pyrenees (2,000 to 2,300 m above sea level). We conduct variogram analyses of 24 snow depth maps derived from terrestrial light detection and ranging scans, acquired during six consecutive snow seasons (2011–2017) that span a range of hydroclimatic conditions. We complement our analyses with bare ground topography data and wind speed and direction measurements. Our results show temporal consistency in the spatial variability of snow depth, with short‐range fractal behavior and scale break lengths that are similar to the optimal search distance (25 m) previously reported for the topographic position index, a terrain‐based predictor of snow depth. Beyond the 25‐m scale break, there is little to no fractal structure. We report a long‐range scale break of the order of 185–300 m for most dates—aligned with the dominant wind direction—and patterns between anisotropies in scale break lengths of shallow snow cover and directional terrain scaling behavior. The temporal consistency of snow depth scaling patterns suggests that, in addition to guiding the spatial configuration of physically based models, fractal analysis could be used to inform the design of independent variables for statistical models used to predict snow depth and its variability.

Prediction of snow water equivalent using artificial neural network and adaptive neuro-fuzzy inference system with two sampling schemes in semi-arid region of Iran

Direct measurement of snow water equivalent (SWE) in snow-dominated mountainous areas is difficult, thus its prediction is essential for water resources management in such areas. In addition, because of nonlinear trend of snow spatial distribution and the multiple influencing factors concerning the SWE spatial distribution, statistical models are not usually able to present acceptable results. Therefore, applicable methods that are able to predict nonlinear trends are necessary. In this research, to predict SWE, the Sohrevard Watershed located in northwest of Iran was selected as the case study. Database was collected, and the required maps were derived. Snow depth (SD) at 150 points with two sampling patterns including systematic random sampling and Latin hypercube sampling (LHS), and snow density at 18 points were randomly measured, and then SWE was calculated. SWE was predicted using artificial neural network (ANN), adaptive neuro-fuzzy inference system (ANFIS) and regression methods. The results showed that the performance of ANN and ANFIS models with two sampling patterns were observed better than the regression method. Moreover, based on most of the efficiency criteria, the efficiency of ANN, ANFIS and regression methods under LHS pattern were observed higher than the systematic random sampling pattern. However, there were no significant differences between the two methods of ANN and ANFIS in SWE prediction. Data of both two sampling patterns had the highest sensitivity to the elevation. In addition, the LHS and the systematic random sampling patterns had the least sensitivity to the profile curvature and plan curvature, respectively.

Snow cover estimation from MODIS and Sentinel-1 SAR data using machine learning algorithms in the western part of the Tianshan Mountains

Obtaining the spatial distribution of snow cover in mountainous areas using the optical image of remote sensing technology is difficult because of cloud and fog. In this study, the object-based principle component analysis-support vector machine (PCA-SVM) method is proposed for snow cover mapping through the integration of moderateresolution imaging spectroradiometer (MODIS) snow cover products and the Sentinel-1 synthetic aperture radar (SAR) scattering characteristics. First, derived from the Sentinel-1A SAR images, the feature parameters, including VV/VH backscatter, scattering entropy, and scattering alpha, were used to describe the variations of snow and non-snow covers. Second, the optimum feature combinations of snow cover were formed from the feature parameters using the principle component analysis (PCA) algorithm. Finally, using the optimum feature combinations, a snow cover map with a 20 m spatial resolution was extracted by means of an object-based SVM classifier. This method was applied in the study area of the Xinyuan County, which is located in the western part of the Tianshan Mountains in Xinjiang, China. The accuracies in this method were analyzed according to the data observed at different experimental sites. Results showed that the snow cover pixels of the extraction were less than those in the actual situation (FB1=93.86, FB2=59.78). The evaluation of the threat score (TS), probability of detection (POD), and false alarm ratio (FAR) for the snow-covered pixels obtained from the two-stage SAR images were different (TS1=86.84, POD1=90.10, FAR1=4.01; TS2=56.40, POD2=57.62, FAR2=3.62). False and misclassifications of the snow cover and non-snow cover pixels were found. Although the classifications were not highly accurate, the approach showed potential for integrating different sources to retrieve the spatial distribution of snow covers during a stable period.

Sensing Area‐Average Snow Water Equivalent with Cosmic‐Ray Neutrons: The Influence of Fractional Snow Cover

AbstractCosmic‐ray neutron sensing (CRNS) is a promising non‐invasive technique to estimate snow water equivalent (SWE) over large areas. In contrast to preliminary studies focusing on shallow snow conditions (SWE 130 mm), more recently the method was shown experimentally to be sensitive also to deeper snowpacks providing the basis for its use at mountain experimental sites. However, hysteretic neutron response has been observed for complex snow cover including patchy snow‐free areas. In the present study we aimed to understand and support the experimental findings using a comprehensive neutron modeling approach. Several simulations have been set up in order to disentangle the effect on the signal of different land surface characteristics and to reproduce multiple observations during periods of snow melt and accumulation. To represent the actual land surface heterogeneity and the complex snow cover, the model used data from terrestrial laser scanning. The results show that the model was able to accurately reproduce the CRNS signal and particularly the hysteresis effect during accumulation and melting periods. Moreover, the sensor footprint was found to be anisotropic and affected by the spatial distribution of liquid water and snow as well as by the topography of the nearby mountains. Under fully snow‐covered conditions the CRNS is able to accurately estimate SWE without prior knowledge about snow density profiles or other spatial anomalies. These results provide new insights into the characteristics of the detected neutron signal in complex terrain and support the use of CRNS for long‐term snow monitoring in high elevated mountain environments.

A Framework to Delineate Precipitation‐Runoff Regimes: Precipitation Versus Snowpack in the Western United States

AbstractSnowpack accumulation/ablation affects the runoff response to precipitation by modulating the water flux reaching the surface. Previous studies mostly focused on “rain‐on‐snow” events. Here we propose a framework to extend the scope and classify precipitation events accompanied by snow accumulation/ablation (precipitation‐and‐snow, or PAS, events) into five regimes. This framework is applied to a regional climate simulation over the western United States for 1981–2015 to reveal regions where daily changes in snowpack alter the surface hydrologic responses to precipitation. Over the western United States, PAS events account for 50–90% of all the precipitation events. Compared to the broad spatial distribution of snow accumulation‐type PAS events, snowmelting‐type PAS events are limited to coastal high‐elevation areas. Atmospheric rivers, a key driver of heavy precipitation in the region, account for only 2% of the PAS events, but they trigger significant snowmelt, accounting for 20% and 11% of light and heavy snowmelting events, respectively.

Snow-driven uncertainty in CryoSat-2-derived Antarctic sea ice thickness – insights from McMurdo Sound

Abstract. Knowledge of the snow depth distribution on Antarctic sea ice is poor but is critical to obtaining sea ice thickness from satellite altimetry measurements of the freeboard. We examine the usefulness of various snow products to provide snow depth information over Antarctic fast ice in McMurdo Sound with a focus on a novel approach using a high-resolution numerical snow accumulation model (SnowModel). We compare this model to results from ECMWF ERA-Interim precipitation, EOS Aqua AMSR-E passive microwave snow depths and in situ measurements at the end of the sea ice growth season in 2011. The fast ice was segmented into three areas by fastening date and the onset of snow accumulation was calibrated to these dates. SnowModel captures the spatial snow distribution gradient in McMurdo Sound and falls within 2 cm snow water equivalent (s.w.e) of in situ measurements across the entire study area. However, it exhibits deviations of 5 cm s.w.e. from these measurements in the east where the effect of local topographic features has caused an overestimate of snow depth in the model. AMSR-E provides s.w.e. values half that of SnowModel for the majority of the sea ice growth season. The coarser-resolution ERA-Interim produces a very high mean s.w.e. value 20 cm higher than the in situ measurements. These various snow datasets and in situ information are used to infer sea ice thickness in combination with CryoSat-2 (CS-2) freeboard data. CS-2 is capable of capturing the seasonal trend of sea ice freeboard growth but thickness results are highly dependent on what interface the retracked CS-2 height is assumed to represent. Because of this ambiguity we vary the proportion of ice and snow that represents the freeboard – a mathematical alteration of the radar penetration into the snow cover – and assess this uncertainty in McMurdo Sound. The ranges in sea ice thickness uncertainty within these bounds, as means of the entire growth season, are 1.08, 4.94 and 1.03 m for SnowModel, ERA-Interim and AMSR-E respectively. Using an interpolated in situ snow dataset we find the best agreement between CS-2-derived and in situ thickness when this interface is assumed to be 0.07 m below the snow surface.