scholarly journals Effect of snow microstructure variability on Ku-band radar snow water equivalent retrievals

2019 ◽  
Author(s):  
Nick Rutter ◽  
Melody J. Sandells ◽  
Chris Derksen ◽  
Joshua King ◽  
Peter Toose ◽  
...  

Abstract. Spatial variability in snowpack properties negatively impacts our capacity to make direct measurements of snow water equivalent (SWE) using satellites. A comprehensive data set of snow microstructure (94 profiles at 36 sites) and snow layer thickness (9000 vertical profiles across 9 trenches) collected over two winters at Trail Valley Creek, NWT, Canada, were applied in synthetic radiative transfer experiments. This allowed robust assessment of the impact of first guess information of snow microstructural characteristics on the viability of SWE retrievals. Depth hoar layer thickness varied over the shortest horizontal distances, controlled by subnivean vegetation and topography, while variability of total snowpack thickness approximated that of wind slab layers. Mean horizontal correlation lengths were sub-metre for all layers. Depth hoar was consistently ~ 30 % of total depth, and with increasing total depth the proportion of wind slab increased at the expense of the decreasing surface snow layer. Distinct differences were evident between distributions of layer properties; a single median value represented density and SSA of each layer well. Spatial variability in microstructure of depth hoar layers dominated SWE retrieval errors. A depth hoar SSA estimate of around 7 % under the median value was needed to accurately retrieve SWE. In shallow snowpacks

2019 ◽  
Vol 13 (11) ◽  
pp. 3045-3059 ◽  
Author(s):  
Nick Rutter ◽  
Melody J. Sandells ◽  
Chris Derksen ◽  
Joshua King ◽  
Peter Toose ◽  
...  

Abstract. Spatial variability in snowpack properties negatively impacts our capacity to make direct measurements of snow water equivalent (SWE) using satellites. A comprehensive data set of snow microstructure (94 profiles at 36 sites) and snow layer thickness (9000 vertical profiles across nine trenches) collected over two winters at Trail Valley Creek, NWT, Canada, was applied in synthetic radiative transfer experiments. This allowed for robust assessment of the impact of estimation accuracy of unknown snow microstructural characteristics on the viability of SWE retrievals. Depth hoar layer thickness varied over the shortest horizontal distances, controlled by subnivean vegetation and topography, while variability in total snowpack thickness approximated that of wind slab layers. Mean horizontal correlation lengths of layer thickness were less than a metre for all layers. Depth hoar was consistently ∼30 % of total depth, and with increasing total depth the proportion of wind slab increased at the expense of the decreasing surface snow layer. Distinct differences were evident between distributions of layer properties; a single median value represented density and specific surface area (SSA) of each layer well. Spatial variability in microstructure of depth hoar layers dominated SWE retrieval errors. A depth hoar SSA estimate of around 7 % under the median value was needed to accurately retrieve SWE. In shallow snowpacks <0.6 m, depth hoar SSA estimates of ±5 %–10 % around the optimal retrieval SSA allowed SWE retrievals within a tolerance of ±30 mm. Where snowpacks were deeper than ∼30 cm, accurate values of representative SSA for depth hoar became critical as retrieval errors were exceeded if the median depth hoar SSA was applied.


2020 ◽  
Vol 21 (11) ◽  
pp. 2713-2733 ◽  
Author(s):  
Graham A. Sexstone ◽  
Colin A. Penn ◽  
Glen E. Liston ◽  
Kelly E. Gleason ◽  
C. David Moeser ◽  
...  

AbstractThis study evaluated the spatial variability of trends in simulated snowpack properties across the Rio Grande headwaters of Colorado using the SnowModel snow evolution modeling system. SnowModel simulations were performed using a grid resolution of 100 m and 3-hourly time step over a 34-yr period (1984–2017). Atmospheric forcing was provided by phase 2 of the North American Land Data Assimilation System, and the simulations accounted for temporal changes in forest canopy from bark beetle and wildfire disturbances. Annual summary values of simulated snowpack properties [snow metrics; e.g., peak snow water equivalent (SWE), snowmelt rate and timing, and snow sublimation] were used to compute trends across the domain. Trends in simulated snow metrics varied depending on elevation, aspect, and land cover. Statistically significant trends did not occur evenly within the basin, and some areas were more sensitive than others. In addition, there were distinct trend differences between the different snow metrics. Upward trends in mean winter air temperature were 0.3°C decade−1, and downward trends in winter precipitation were −52 mm decade−1. Middle elevation zones, coincident with the greatest volumetric snow water storage, exhibited the greatest sensitivity to changes in peak SWE and snowmelt rate. Across the Rio Grande headwaters, snowmelt rates decreased by 20% decade−1, peak SWE decreased by 14% decade−1, and total snowmelt quantity decreased by 13% decade−1. These snow trends are in general agreement with widespread snow declines that have been reported for this region. This study further quantifies these snow declines and provides trend information for additional snow variables across a greater spatial coverage at finer spatial resolution.


2014 ◽  
Vol 8 (2) ◽  
pp. 417-437 ◽  
Author(s):  
C. M. Carmagnola ◽  
S. Morin ◽  
M. Lafaysse ◽  
F. Domine ◽  
B. Lesaffre ◽  
...  

Abstract. In the SURFEX/ISBA-Crocus multi-layer snowpack model, the snow microstructure has up to now been characterised by the grain size and by semi-empirical shape variables which cannot be measured easily in the field or linked to other relevant snow properties. In this work we introduce a new formulation of snow metamorphism directly based on equations describing the rate of change of the optical diameter (dopt). This variable is considered here to be equal to the equivalent sphere optical diameter, which is inversely proportional to the specific surface area (SSA). dopt thus represents quantitatively some of the geometric characteristics of a porous medium. Different prognostic rate equations of dopt, including a re-formulation of the original Crocus scheme and the parameterisations from Taillandier et al. (2007) and Flanner and Zender (2006), were evaluated by comparing their predictions to field measurements carried out at Summit Camp (Greenland) in May and June 2011 and at Col de Porte (French Alps) during the 2009/10 and 2011/12 winter seasons. We focused especially on results in terms of SSA. In addition, we tested the impact of the different formulations on the simulated density profile, the total snow height, the snow water equivalent (SWE) and the surface albedo. Results indicate that all formulations perform well, with median values of the RMSD between measured and simulated SSA lower than 10 m2 kg−1. Incorporating the optical diameter as a fully fledged prognostic variable is an important step forward in the quantitative description of the snow microstructure within snowpack models, because it opens the way to data assimilation of various electromagnetic observations.


2021 ◽  
Author(s):  
Leung Tsang ◽  
Michael Durand ◽  
Chris Derksen ◽  
Ana P. Barros ◽  
Do-Hyuk Kang ◽  
...  

Abstract. Seasonal snow cover is the largest single component of the cryosphere in areal extent, covering an average of 46 million square km of Earth's surface (31 % of the land area) each year, and is thus an important expression of and driver of the Earth’s climate. In recent years, Northern Hemisphere spring snow cover has been declining at about the same rate (~ −13 %/decade) as Arctic summer sea ice. More than one-sixth of the world’s population relies on seasonal snowpack and glaciers for a water supply that is likely to decrease this century. Snow is also a critical component of Earth’s cold regions' ecosystems, in which wildlife, vegetation, and snow are strongly interconnected. Snow water equivalent (SWE) describes the quantity of snow stored on the land surface and is of fundamental importance to water, energy, and geochemical cycles. Quality global SWE estimates are lacking. Given the vast seasonal extent combined with the spatially variable nature of snow distribution at regional and local scales, surface observations will not be able to provide sufficient SWE information. Satellite observations presently cannot provide SWE information at the spatial and temporal resolutions required to address science and high socio-economic value applications such as water resource management and streamflow forecasting. In this paper, we review the potential contribution of X- and Ku-Band Synthetic Aperture Radar (SAR) for global monitoring of SWE. We describe radar interactions with snow-covered landscapes, characterization of snowpack properties using radar measurements, and refinement of retrieval algorithms via synergy with other microwave remote sensing approaches. SAR can image the surface during both day and night regardless of cloud cover, allowing high-frequency revisit at high spatial resolution as demonstrated by missions such as Sentinel-1. The physical basis for estimating SWE from X- and Ku-band radar measurements at local scales is volume scattering by millimetre-scale snow grains. Inference of global snow properties from SAR requires an interdisciplinary approach based on field observations of snow microstructure, physical snow modelling, electromagnetic theory, and retrieval strategies over a range of scales. New field measurement capabilities have enabled significant advances in understanding snow microstructure such as grain size, densities, and layering. We describe radar interactions with snow-covered landscapes, the characterization of snowpack properties using radar measurements, and the refinement of retrieval algorithms via synergy with other microwave remote sensing approaches. This review serves to inform the broader snow research, monitoring, and applications communities on progress made in recent decades, and sets the stage for a new era in SWE remote-sensing from SAR measurements.


2013 ◽  
Vol 7 (5) ◽  
pp. 4443-4500 ◽  
Author(s):  
C. M. Carmagnola ◽  
S. Morin ◽  
M. Lafaysse ◽  
F. Domine ◽  
B. Lesaffre ◽  
...  

Abstract. In the SURFEX/ISBA-Crocus multi-layer snowpack model, the snow microstructure was up to now characterized by the grain size and by semi-empirical shape variables which cannot be measured easily in the field or linked to other relevant snow properties. In this work we introduce a new formulation of snow metamorphism directly based on equations describing the rate of change of the optical diameter (dopt). This variable is considered here to be equal to the equivalent sphere optical diameter, which is inversely proportional to the specific surface area (SSA). dopt thus represents quantitatively some of the geometric characteristics of a porous medium. Different prognostic rate equations of dopt, including a re-formulation of the original Crocus scheme and the parametrizations from Taillandier et al. (2007) and Flanner and Zender (2006), were evaluated by comparing their predictions to field measurements carried out at Summit Camp (Greenland) in May and June 2011 and at Col de Porte (French Alps) during the 2009/10 and 2011/12 winter seasons. We focused especially on results in terms of SSA. In addition, we tested the impact of the different formulations on the simulated density profile, the total snow height, the snow water equivalent (SWE) and the surface albedo. Results indicate that all formulations perform well, with median values of the RMSD between measured and simulated SSA lower than 10 m2 kg−1. Incorporating the optical diameter as a fully-fledged prognostic variable is an important step forward in the quantitative description of the snow microstructure within snowpack models, because it opens the way to data assimilation of various electromagnetic observations.


2021 ◽  
Author(s):  
Ilaria Clemenzi ◽  
David Gustafsson ◽  
Jie Zhang ◽  
Björn Norell ◽  
Wolf Marchand ◽  
...  

&lt;p&gt;Snow in the mountains is the result of the interplay between meteorological conditions, e.g., precipitation, wind and solar radiation, and landscape features, e.g., vegetation and topography. For this reason, it is highly variable in time and space. It represents an important water storage for several sectors of the society including tourism, ecology and hydropower. The estimation of the amount of snow stored in winter and available in the form of snowmelt runoff can be strategic for their sustainability. In the hydropower sector, for example, the occurrence of higher snow and snowmelt runoff volumes at the end of the spring and in the early summer compared to the estimated one can substantially impact reservoir regulation with energy and economical losses. An accurate estimation of the snow volumes and their spatial and temporal distribution is thus essential for spring flood runoff prediction. Despite the increasing effort in the development of new acquisition techniques, the availability of extensive and representative snow and density measurements for snow water equivalent estimations is still limited. Hydrological models in combination with data assimilation of ground or remote sensing observations is a way to overcome these limitations. However, the impact of using different types of snow observations on snowmelt runoff predictions is, little understood. In this study we investigated the potential of assimilating in situ and remote sensing snow observations to improve snow water equivalent estimates and snowmelt runoff predictions. We modelled the seasonal snow water equivalent distribution in the Lake &amp;#214;veruman catchment, Northern Sweden, which is used for hydropower production. Simulations were performed using the semi-distributed hydrological model HYPE for the snow seasons 2017-2020. For this purpose, a snowfall distribution model based on wind-shelter factors was included to represent snow spatial distribution within model units. The units consist of 2.5x2.5 km&lt;sup&gt;2&lt;/sup&gt; grid cells, which were further divided into hydrological response units based on elevation, vegetation and aspect. The impact on the estimation of the total catchment mean snow water equivalent and snowmelt runoff volume were evaluated using for data assimilation, gpr-based snow water equivalent data acquired along survey lines in the catchment in the early spring of the four years, snow water equivalent data obtained by a machine learning algorithm and satellite-based fractional snow cover data. Results show that the wind-shelter based snow distribution model was able to represent a similar spatial distribution as the gpr survey lines, when assessed on the catchment level. Deviations in the model performance within and between specific gpr survey lines indicate issues with the spatial distribution of input precipitation, and/or need to include explicit representation of snow drift between model units. The explicit snow distribution model also improved runoff simulations, and the ability of the model to improve forecast through data assimilation.&lt;/p&gt;


2013 ◽  
Vol 17 (7) ◽  
pp. 2781-2796 ◽  
Author(s):  
S. Shukla ◽  
J. Sheffield ◽  
E. F. Wood ◽  
D. P. Lettenmaier

Abstract. Global seasonal hydrologic prediction is crucial to mitigating the impacts of droughts and floods, especially in the developing world. Hydrologic predictability at seasonal lead times (i.e., 1–6 months) comes from knowledge of initial hydrologic conditions (IHCs) and seasonal climate forecast skill (FS). In this study we quantify the contributions of two primary components of IHCs – soil moisture and snow water content – and FS (of precipitation and temperature) to seasonal hydrologic predictability globally on a relative basis throughout the year. We do so by conducting two model-based experiments using the variable infiltration capacity (VIC) macroscale hydrology model, one based on ensemble streamflow prediction (ESP) and another based on Reverse-ESP (Rev-ESP), both for a 47 yr re-forecast period (1961–2007). We compare cumulative runoff (CR), soil moisture (SM) and snow water equivalent (SWE) forecasts from each experiment with a VIC model-based reference data set (generated using observed atmospheric forcings) and estimate the ratio of root mean square error (RMSE) of both experiments for each forecast initialization date and lead time, to determine the relative contribution of IHCs and FS to the seasonal hydrologic predictability. We find that in general, the contributions of IHCs to seasonal hydrologic predictability is highest in the arid and snow-dominated climate (high latitude) regions of the Northern Hemisphere during forecast periods starting on 1 January and 1 October. In mid-latitude regions, such as the Western US, the influence of IHCs is greatest during the forecast period starting on 1 April. In the arid and warm temperate dry winter regions of the Southern Hemisphere, the IHCs dominate during forecast periods starting on 1 April and 1 July. In equatorial humid and monsoonal climate regions, the contribution of FS is generally higher than IHCs through most of the year. Based on our findings, we argue that despite the limited FS (mainly for precipitation) better estimates of the IHCs could lead to improvement in the current level of seasonal hydrologic forecast skill over many regions of the globe at least during some parts of the year.


2016 ◽  
Author(s):  
Jean M. Bergeron ◽  
Mélanie Trudel ◽  
Robert Leconte

Abstract. The potential of data assimilation for hydrologic predictions has been demonstrated in many research studies. Watersheds over which multiple observation types are available can potentially further benefit from data assimilation by having multiple updated states from which hydrologic predictions can be generated. However, the magnitude and time span of the impact of the assimilation of an observation varies according not only to its type, but also to the variables included in the state vector. This study examines the impact of multivariate synthetic data assimilation using the Ensemble Kalman Filter (EnKF) into the spatially distributed hydrologic model CEQUEAU for the mountainous Nechako River located in British-Columbia, Canada. Synthetic data includes daily snow cover area (SCA), daily measurements of snow water equivalent (SWE) at three different locations and daily streamflow data at the watershed outlet. Results show a large variability of the continuous rank probability skill score over a wide range of prediction horizons (days to weeks) depending on the state vector configuration and the type of observations assimilated. Overall, the variables most closely linearly linked to the observations are the ones worth considering adding to the state vector. The performance of the assimilation of basin-wide SCA, which does not have a decent proxy among potential state variables, does not surpass the open loop for any of the simulated variables. However, the assimilation of streamflow offers major improvements steadily throughout the year, but mainly over the short-term (up to 5 days) forecast horizons, while the impact of the assimilation of SWE gains more importance during the snowmelt period over the mid-term (up to 50 days) forecast horizon compared with open loop. The combined assimilation of streamflow and SWE performs better than its individual counterparts, offering improvements over all forecast horizons considered and throughout the whole year, including the critical period of snowmelt. This highlights the potential benefit of using multivariate data assimilation for streamflow predictions in snow-dominated regions.


2020 ◽  
Author(s):  
Bertrand Cluzet ◽  
Matthieu Lafaysse ◽  
Emmanuel Cosme ◽  
Clément Albergel ◽  
Louis-François Meunier ◽  
...  

Abstract. Monitoring the evolution of the snowpack properties in mountainous areas is crucial for avalanche hazard forecasting and water resources management. In-situ and remotely sensed observations provide precious information on the snowpack but usually offer a limited spatio-temporal coverage of bulk or surface variables only. In particular, visible-near infrared (VIS-NIR) reflectance observations can inform on the snowpack surface properties but are limited by terrain shading and clouds. Snowpack modelling enables to estimate any physical variable, virtually anywhere, but is affected by large errors and uncertainties. Data assimilation offers a way to combine both sources of information, and to propagate information from observed areas to non observed areas. Here, we present CrocO, (Crocus-Observations) an ensemble data assimilation system able to ingest any snowpack observation (applied as a first step to the height of snow (HS) and VIS-NIR reflectances) in a spatialised geometry. CrocO uses an ensemble of snowpack simulations to represent modelling uncertainties, and a Particle Filter (PF) to reduce them. The PF is known to collapse when assimilating a too large number of observations. Two variants of the PF were specifically implemented to ensure that observations information is propagated in space while tackling this issue. The global algorithm ingests all available observations with an iterative inflation of observation errors, while the klocal algorithm is a localised approach performing a selection of the observations to assimilate based on background correlation patterns. Experiments are carried out in a twin experiment setup, with synthetic observations of HS and VIS-NIR reflectances available in only a 1/6th of the simulation domain. Results show that compared against runs without assimilation, analyses exhibit an average improvement of snow water equivalent Continuous Rank Probability Score (CRPS) of 60 % when assimilating HS with a 40-member ensemble, and an average 20 % CRPS improvement when assimilating reflectance with a 160-member ensemble. Significant improvements are also obtained outside the observation domain. These promising results open a way for the assimilation of real observations of reflectance, or of any snowpack observations in a spatialised context.


2012 ◽  
Vol 4 (1) ◽  
pp. 13-21 ◽  
Author(s):  
S. Morin ◽  
Y. Lejeune ◽  
B. Lesaffre ◽  
J.-M. Panel ◽  
D. Poncet ◽  
...  

Abstract. A quality-controlled snow and meteorological dataset spanning the period 1 August 1993–31 July 2011 is presented, originating from the experimental station Col de Porte (1325 m altitude, Chartreuse range, France). Emphasis is placed on meteorological data relevant to the observation and modelling of the seasonal snowpack. In-situ driving data, at the hourly resolution, consist of measurements of air temperature, relative humidity, windspeed, incoming short-wave and long-wave radiation, precipitation rate partitioned between snow- and rainfall, with a focus on the snow-dominated season. Meteorological data for the three summer months (generally from 10 June to 20 September), when the continuity of the field record is not warranted, are taken from a local meteorological reanalysis (SAFRAN), in order to provide a continuous and consistent gap-free record. Data relevant to snowpack properties are provided at the daily (snow depth, snow water equivalent, runoff and albedo) and hourly (snow depth, albedo, runoff, surface temperature, soil temperature) time resolution. Internal snowpack information is provided from weekly manual snowpit observations (mostly consisting in penetration resistance, snow type, snow temperature and density profiles) and from a hourly record of temperature and height of vertically free ''settling'' disks. This dataset has been partially used in the past to assist in developing snowpack models and is presented here comprehensively for the purpose of multi-year model performance assessment. The data is placed on the PANGAEA repository (http://dx.doi.org/10.1594/PANGAEA.774249) as well as on the public ftp server ftp://ftp-cnrm.meteo.fr/pub-cencdp/.


Sign in / Sign up

Export Citation Format

Share Document