scholarly journals Multilayer observation and estimation of the snowpack cold content in a humid boreal coniferous forest of eastern Canada

2021 ◽  
Vol 15 (12) ◽  
pp. 5371-5386
Author(s):  
Achut Parajuli ◽  
Daniel F. Nadeau ◽  
François Anctil ◽  
Marco Alves

Abstract. Cold content (CC) is an internal energy state within a snowpack and is defined by the energy deficit required to attain isothermal snowmelt temperature (0 ∘C). Cold content for a given snowpack thus plays a critical role because it affects both the timing and the rate of snowmelt. Measuring cold content is a labour-intensive task as it requires extracting in situ snow temperature and density. Hence, few studies have focused on characterizing this snowpack variable. This study describes the multilayer cold content of a snowpack and its variability across four sites with contrasting canopy structures within a coniferous boreal forest in southern Québec, Canada, throughout winter 2017–2018. The analysis was divided into two steps. In the first step, the observed CC data from weekly snowpits for 60 % of the snow cover period were examined. During the second step, a reconstructed time series of modelled CC was produced and analyzed to highlight the high-resolution temporal variability of CC for the full snow cover period. To accomplish this, the Canadian Land Surface Scheme (CLASS; featuring a single-layer snow model) was first implemented to obtain simulations of the average snow density at each of the four sites. Next, an empirical procedure was used to produce realistic density profiles, which, when combined with in situ continuous snow temperature measurements from an automatic profiling station, provides a time series of CC estimates at half-hour intervals for the entire winter. At the four sites, snow persisted on the ground for 218 d, with melt events occurring on 42 of those days. Based on snowpit observations, the largest mean CC (−2.62 MJ m−2) was observed at the site with the thickest snow cover. The maximum difference in mean CC between the four study sites was −0.47 MJ m−2, representing a site-to-site variability of 20 %. Before analyzing the reconstructed CC time series, a comparison with snowpit data confirmed that CLASS yielded reasonable bulk estimates of snow water equivalent (SWE) (R2=0.64 and percent bias (Pbias) =-17.1 %), snow density (R2=0.71 and Pbias =1.6 %), and cold content (R2=0.93 and Pbias =-3.3 %). A snow density profile derived by utilizing an empirical formulation also provided reasonable estimates of layered cold content (R2=0.42 and Pbias =5.17 %). Thanks to these encouraging results, the reconstructed and continuous CC series could be analyzed at the four sites, revealing the impact of rain-on-snow and cold air pooling episodes on the variation of CC. The continuous multilayer cold content time series also provided us with information about the effect of stand structure, local topography, and meteorological conditions on cold content variability. Additionally, a weak relationship between canopy structure and CC was identified.

2021 ◽  
Author(s):  
Achut Parajuli ◽  
Daniel F. Nadeau ◽  
François Anctil ◽  
Marco Alves

Abstract. Cold content (CC) is an internal energy state within a snowpack and is defined by the energy deficit required to attain isothermal snowmelt temperature (0 °C). For any snowpack, fulfilling the cold content deficit is a pre-requisite before the onset of the snowmelt. Cold content for a given snowpack thus plays a critical role because it affects both the timing and the rate of snowmelt. Estimating the cold content is a labour-intensive task as it requires extracting in-situ snow temperature and density. Hence, few studies have focused on characterizing this snowpack variable. This study describes the multilayer cold content of a snowpack and its variability across four sites with contrasting canopy structures within a coniferous boreal forest in southern Québec, Canada, throughout winter 2017–18. The analysis was divided into two steps. In the first step, the observed CC data from weekly snowpits for 60 % of the snow cover period were examined. During the second step, a reconstructed time series of CC was produced and analyzed to highlight the high-resolution temporal variability of CC for the full snow cover period. To accomplish this, the Canadian Land Surface Scheme (CLASS; featuring a single-layer snow model) was first implemented to obtain simulations of the average snow density at each of the four sites. Next, an empirical procedure was used to produce realistic density profiles, which, when combined with in situ continuous snow temperature measurements from an automatic profiling station, provides a time series of CC estimates at half-hour intervals for the entire winter. At the four sites, snow persisted on the ground for 218 days, with melt events occurring on 42 of those days. Based on snowpit observations, the largest mean CC (−2.62 MJ m−2) was observed at the site with the thickest snow cover. The maximum difference in mean CC between the four study sites was −0.47 MJ m−2, representing a site-to-site variability of 20 %. Before analyzing the reconstructed CC time series, a comparison with snowpit data confirmed that CLASS yielded reasonable estimates of the snow water equivalent (SWE) (R2 = 0.64 and percent bias (Pbias) = −17.1 %), bulk snow density (R2 = 0.71 and Pbias = 1.6 %), and bulk cold content (R2 = 0.90 and Pbias = −2.0 %). A snow density profile derived by utilizing an empirical formulation also provided reasonable estimates of cold content (R2 = 0.42 and Pbias = 5.17 %). Thanks to these encouraging results, the reconstructed and continuous CC series could be analyzed at the four sites, revealing the impact of rain-on-snow and cold air pooling episodes on the variation of CC. The continuous multilayer cold content time series also provided us with information about the effect of stand structure, local topography, and meteorological conditions on cold content variability. Additionally, a weak relationship between canopy structure and CC was identified.


Biologia ◽  
2014 ◽  
Vol 69 (11) ◽  
Author(s):  
Martin Bartík ◽  
Roman Sitko ◽  
Marek Oreňák ◽  
Juraj Slovik ◽  
Jaroslav Škvarenina

AbstractIn the presented paper we deal with the impact of the mature spruce stand on the accumulation and melting of snow cover at Červenec research area located in the Western Tatras at an elevation of 1420 m a.s.l. The work analyses the data obtained from the monitoring of snow cover during the period 2009–2014 (6 seasons). Since the season 2012/2013 the measurements have been also performed in a dead part of the stand and in a meadow. The results proved significant impact of the spruce stand on hydro-physical characteristics of snow cover — snow water equivalent, snow density, as well as on their change due to the dieback of the stand. The data measured at individual locations (open space in the forest, open meadow area, living and dead forest) were tested with the paired t-test for the significance of average differences. Average snow water equivalent in the living forest, dead forest and meadow was 42%, 47% and 83% of the reference value measured at the open space in the forest, respectively. The process of snow accumulation and melting was fastest at the open space, followed by the dead forest. In the living forest, the processes were the slowest.


2017 ◽  
Vol 11 (4) ◽  
pp. 1647-1664 ◽  
Author(s):  
Emmy E. Stigter ◽  
Niko Wanders ◽  
Tuomo M. Saloranta ◽  
Joseph M. Shea ◽  
Marc F. P. Bierkens ◽  
...  

Abstract. Snow is an important component of water storage in the Himalayas. Previous snowmelt studies in the Himalayas have predominantly relied on remotely sensed snow cover. However, snow cover data provide no direct information on the actual amount of water stored in a snowpack, i.e., the snow water equivalent (SWE). Therefore, in this study remotely sensed snow cover was combined with in situ observations and a modified version of the seNorge snow model to estimate (climate sensitivity of) SWE and snowmelt runoff in the Langtang catchment in Nepal. Snow cover data from Landsat 8 and the MOD10A2 snow cover product were validated with in situ snow cover observations provided by surface temperature and snow depth measurements resulting in classification accuracies of 85.7 and 83.1 % respectively. Optimal model parameter values were obtained through data assimilation of MOD10A2 snow maps and snow depth measurements using an ensemble Kalman filter (EnKF). Independent validations of simulated snow depth and snow cover with observations show improvement after data assimilation compared to simulations without data assimilation. The approach of modeling snow depth in a Kalman filter framework allows for data-constrained estimation of snow depth rather than snow cover alone, and this has great potential for future studies in complex terrain, especially in the Himalayas. Climate sensitivity tests with the optimized snow model revealed that snowmelt runoff increases in winter and the early melt season (December to May) and decreases during the late melt season (June to September) as a result of the earlier onset of snowmelt due to increasing temperature. At high elevation a decrease in SWE due to higher air temperature is (partly) compensated by an increase in precipitation, which emphasizes the need for accurate predictions on the changes in the spatial distribution of precipitation along with changes in temperature.


2020 ◽  
Vol 24 (8) ◽  
pp. 4061-4090 ◽  
Author(s):  
Silvia Terzago ◽  
Valentina Andreoli ◽  
Gabriele Arduini ◽  
Gianpaolo Balsamo ◽  
Lorenzo Campo ◽  
...  

Abstract. Snow models are usually evaluated at sites providing high-quality meteorological data, so that the uncertainty in the meteorological input data can be neglected when assessing model performances. However, high-quality input data are rarely available in mountain areas and, in practical applications, the meteorological forcing used to drive snow models is typically derived from spatial interpolation of the available in situ data or from reanalyses, whose accuracy can be considerably lower. In order to fully characterize the performances of a snow model, the model sensitivity to errors in the input data should be quantified. In this study we test the ability of six snow models to reproduce snow water equivalent, snow density and snow depth when they are forced by meteorological input data with gradually lower accuracy. The SNOWPACK, GEOTOP, HTESSEL, UTOPIA, SMASH and S3M snow models are forced, first, with high-quality measurements performed at the experimental site of Torgnon, located at 2160 m a.s.l. in the Italian Alps (control run). Then, the models are forced by data at gradually lower temporal and/or spatial resolution, obtained by (i) sampling the original Torgnon 30 min time series at 3, 6, and 12 h, (ii) spatially interpolating neighbouring in situ station measurements and (iii) extracting information from GLDAS, ERA5 and ERA-Interim reanalyses at the grid point closest to the Torgnon site. Since the selected models are characterized by different degrees of complexity, from highly sophisticated multi-layer snow models to simple, empirical, single-layer snow schemes, we also discuss the results of these experiments in relation to the model complexity. The results show that, when forced by accurate 30 min resolution weather station data, the single-layer, intermediate-complexity snow models HTESSEL and UTOPIA provide similar skills to the more sophisticated multi-layer model SNOWPACK, and these three models show better agreement with observations and more robust performances over different seasons compared to the lower-complexity models SMASH and S3M. All models forced by 3-hourly data provide similar skills to the control run, while the use of 6- and 12-hourly temporal resolution forcings may lead to a reduction in model performances if the incoming shortwave radiation is not properly represented. The SMASH model generally shows low sensitivity to the temporal degradation of the input data. Spatially interpolated data from neighbouring stations and reanalyses are found to be adequate forcings, provided that temperature and precipitation variables are not affected by large biases over the considered period. However, a simple bias-adjustment technique applied to ERA-Interim temperatures allowed all models to achieve similar performances to the control run. Regardless of their complexity, all models show weaknesses in the representation of the snow density.


2014 ◽  
Vol 15 (2) ◽  
pp. 551-562 ◽  
Author(s):  
Jiarui Dong ◽  
Mike Ek ◽  
Dorothy Hall ◽  
Christa Peters-Lidard ◽  
Brian Cosgrove ◽  
...  

Abstract Understanding and quantifying satellite-based, remotely sensed snow cover uncertainty are critical for its successful utilization. The Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover errors have been previously recognized to be associated with factors such as cloud contamination, snowpack grain sizes, vegetation cover, and topography; however, the quantitative relationship between the retrieval errors and these factors remains elusive. Joint analysis of the MODIS fractional snow cover (FSC) from Collection 6 (C6) and in situ air temperature and snow water equivalent measurements provides a unique look at the error structure of the MODIS C6 FSC products. Analysis of the MODIS FSC dataset over the period from 2000 to 2005 was undertaken over the continental United States (CONUS) with an extensive observational network. When compared to MODIS Collection 5 (C5) snow cover area, the MODIS C6 FSC product demonstrates a substantial improvement in detecting the presence of snow cover in Nevada [30% increase in probability of detection (POD)], especially in the early and late snow seasons; some improvement over California (10% POD increase); and a relatively small improvement over Colorado (2% POD increase). However, significant spatial and temporal variations in accuracy still exist, and a proxy is required to adequately predict the expected errors in MODIS C6 FSC retrievals. A relationship is demonstrated between the MODIS FSC retrieval errors and temperature over the CONUS domain, captured by a cumulative double exponential distribution function. This relationship is shown to hold for both in situ and modeled daily mean air temperature. Both of them are useful indices in filtering out the misclassification of MODIS snow cover pixels and in quantifying the errors in the MODIS C6 product for various hydrological applications.


2014 ◽  
Vol 11 (11) ◽  
pp. 12531-12571 ◽  
Author(s):  
S. Gascoin ◽  
O. Hagolle ◽  
M. Huc ◽  
L. Jarlan ◽  
J.-F. Dejoux ◽  
...  

Abstract. The seasonal snow in the Pyrenees is critical for hydropower production, crop irrigation and tourism in France, Spain and Andorra. Complementary to in situ observations, satellite remote sensing is useful to monitor the effect of climate on the snow dynamics. The MODIS daily snow products (Terra/MOD10A1 and Aqua/MYD10A1) are widely used to generate snow cover climatologies, yet it is preferable to assess their accuracies prior to their use. Here, we use both in situ snow observations and remote sensing data to evaluate the MODIS snow products in the Pyrenees. First, we compare the MODIS products to in situ snow depth (SD) and snow water equivalent (SWE) measurements. We estimate the values of the SWE and SD best detection thresholds to 40 mm water equivalent (we) and 105 mm respectively, for both MOD10A1 and MYD10A1. Kappa coefficients are within 0.74 and 0.92 depending on the product and the variable. Then, a set of Landsat images is used to validate MOD10A1 and MYD10A1 for 157 dates between 2002 and 2010. The resulting accuracies are 97% (κ = 0.85) for MOD10A1 and 96% (κ = 0.81) for MYD10A1, which indicates a good agreement between both datasets. The effect of vegetation on the results is analyzed by filtering the forested areas using a land cover map. As expected, the accuracies decreases over the forests but the agreement remains acceptable (MOD10A1: 96%, κ = 0.77; MYD10A1: 95%, κ = 0.67). We conclude that MODIS snow products have a sufficient accuracy for hydroclimate studies at the scale of the Pyrenees range. Using a gapfilling algorithm we generate a consistent snow cover climatology, which allows us to compute the mean monthly snow cover duration per elevation band. We finally analyze the snow patterns for the atypical winter 2011–2012. Snow cover duration anomalies reveal a deficient snowpack on the Spanish side of the Pyrenees, which seems to have caused a drop in the national hydropower production.


2021 ◽  
Author(s):  
Colleen Mortimer ◽  
Lawrence Mudryk ◽  
Chris Derksen ◽  
Kari Luojus ◽  
Pinja Venalainen ◽  
...  

<p>The European Space Agency Snow CCI+ project provides global homogenized long time series of daily snow extent and snow water equivalent (SWE). The Snow CCI SWE product is built on the Finish Meteorological Institute's GlobSnow algorithm, which combines passive microwave data with in situ snow depth information to estimate SWE. The CCI SWE product improves upon previous versions of GlobSnow through targeted changes to the spatial resolution, ancillary data, and snow density parameterization.</p><p>Previous GlobSnow SWE products used a constant snow density of 0.24 kg m<sup>-3</sup> to convert snow depth to SWE. The CCI SWE product applies spatially and temporally varying density fields, derived by krigging in situ snow density information from historical snow transects to correct biases in estimated SWE. Grid spacing was improved from 25 km to 12.5 km by applying an enhanced spatial resolution microwave brightness temperature dataset. We assess step-wise how each of these targeted changes acts to improve or worsen the product by evaluating with snow transect measurements and comparing hemispheric snow mass and trend differences.</p><p>Together, when compared to GlobSnow v3, these changes improved RMSE by ~5 cm and correlation by ~0.1 against a suite of snow transect measurements from Canada, Finland, and Russia. Although the hemispheric snow mass anomalies of CCI SWE and GlobSnow v3 are similar, there are sizeable differences in the climatological SWE, most notably a one month delay in the timing of peak SWE and lower SWE during the accumulation season. These shifts were expected because the variable snow density is lower than the former fixed value of 0.24 kg m<sup>-3</sup> early in the snow season, but then increases over the course of the snow season. We also examine intermediate products to determine the relative improvements attributable solely to the increased spatial resolution versus changes due to the snow density parameterizations. Such systematic evaluations are critical to directing future product development.</p>


2021 ◽  
pp. 117-127
Author(s):  
M. V. GEORGIEVSKY ◽  
◽  
N. I. GOROSHKOVA ◽  
V. A. KHOMYAKOVA ◽  
A. V. STRIZHENOK

The article presents an analysis of the impact of climate change on the main characteristics of ice phenomena, snow cover and the water regime in the Small Northern Dvina River basin occurring in recent decades. Recently, a significant climate warming has been observed in the basin. As a result, winters are getting warmer and shorter. There is also an increase in winter precipitation and the number of thaws. Climate warming directly affects the duration of snow cover, which decreases both due to the later formation and to the earlier destruction of snow. There is also a slight downward trend in the annual values of the maximum snow water equivalent, which may be the result of an increase in the number of thaws in winter, when a part of the snow cover melts contributing to the winter river runoff. The analysis of the main characteristics of the ice cover on the rivers of the studied basin shows that their changes are similarly to changes in the snow cover: there is a reduction in the freeze-up period due to its later formation and earlier complete destruction. The maximum ice thickness on the rivers of the basin also tends to decrease. There is an increase in winter and a decrease in spring runoff. Predictive estimates of changes in the observed trends in the future are presented in the fi nal part of the article based on the CMIP5 project data.


2020 ◽  
Author(s):  
Gabriele Schwaizer ◽  
Lars Keuris ◽  
Thomas Nagler ◽  
Chris Derksen ◽  
Kari Luojus ◽  
...  

<p>Seasonal snow is an important component of the global climate system. It is highly variable in space and time and sensitive to short term synoptic scale processes and long term climate-induced changes of temperature and precipitation. Current snow products derived from various satellite data applying different algorithms show significant discrepancies in extent and snow mass, a potential source for biases in climate monitoring and modelling. The recently launched ESA CCI+ Programme addresses seasonal snow as one of 9 Essential Climate Variables to be derived from satellite data.</p><p>In the snow_cci project, scheduled for 2018 to 2021 in its first phase, reliable fully validated processing lines are developed and implemented. These tools are used to generate homogeneous multi-sensor time series for the main parameters of global snow cover focusing on snow extent and snow water equivalent. Using GCOS guidelines, the requirements for these parameters are assessed and consolidated using the outcome of workshops and questionnaires addressing users dealing with different climate applications. Snow extent product generation applies algorithms accounting for fractional snow extent and cloud screening in order to generate consistent daily products for snow on the surface (viewable snow) and snow on the surface corrected for forest masking (snow on ground) with global coverage. Input data are medium resolution optical satellite images (AVHRR-2/3, AATSR, MODIS, VIIRS, SLSTR/OLCI) from 1981 to present. An iterative development cycle is applied including homogenisation of the snow extent products from different sensors by minimizing the bias. Independent validation of the snow products is performed for different seasons and climate zones around the globe from 1985 onwards, using as reference high resolution snow maps from Landsat and Sentinel- 2as well as in-situ snow data following standardized validation protocols.</p><p>Global time series of daily snow water equivalent (SWE) products are generated from passive microwave data from SMMR, SSM/I, and AMSR from 1978 onwards, combined with in-situ snow depth measurements. Long-term stability and quality of the product is assessed using independent snow survey data and by intercomparison with the snow information from global land process models.</p><p>The usability of the snow_cci products is ensured through the Climate Research Group, which performs case studies related to long term trends of seasonal snow, performs evaluations of CMIP-6 and other snow-focused climate model experiments, and applies the data for simulation of Arctic hydrological regimes.</p><p>In this presentation, we summarize the requirements and product specifications for the snow extent and SWE products, with a focus on climate applications. We present an overview of the algorithms and systems for generation of the time series. The 40 years (from 1980 onwards) time series of daily fractional snow extent products from AVHRR with 5 km pixel spacing, and the 20-year time series from MODIS (1 km pixel spacing) as well as the coarse resolution (25 km pixel spacing) of daily SWE products from 1978 onwards will be presented along with first results of the multi-sensor consistency checks and validation activities.</p>


2016 ◽  
Vol 17 (5) ◽  
pp. 1467-1488 ◽  
Author(s):  
Reinel Sospedra-Alfonso ◽  
Lawrence Mudryk ◽  
William Merryfield ◽  
Chris Derksen

Abstract The ability of the Canadian Seasonal to Interannual Prediction System (CanSIPS) to provide realistic forecast initial conditions for snow cover is assessed using in situ measurements and gridded snow analyses. Forecast initial conditions for snow in CanCM3 and CanCM4 employed by CanSIPS are determined by the response of the Canadian Land Surface Scheme (CLASS) used in both models to forcing from model atmospheric fields constrained by assimilation of 6-hourly reanalysis data. These snow initial conditions are found to be representative of the daily climatology of snow water equivalent (SWE) as well as interannual variations in maximum SWE and the timing of snow onset and snowmelt observed at eight in situ measurement sites located across Canada. The level of this agreement is similar to that of three independent gridded snow analyses (MERRA, the European Space Agency’s GlobSnow, and an offline forced version of CLASS). Total Northern Hemisphere snow mass generated by the CanSIPS initialization procedure is larger for both models (especially CanCM3) than in MERRA, mostly because of higher SWE in regions of common snow cover. Globally, the interannual variability of initial SWE is found to correlate highly with that of MERRA in locations with appreciable snow. These initial values are compared to SWE in freely running CanCM3 and CanCM4 simulations produced without data assimilation of atmospheric fields. Differences in climatological SWE relative to MERRA are similar in the freely running and assimilating CanCM3 and CanCM4 simulations, suggesting that inherent model biases are a major contributor to biases in CanSIPS snow initial conditions.


Sign in / Sign up

Export Citation Format

Share Document