Evaluating methods for reconstructing large gaps in historic snow depth time series
Open Access
- 24 November 2021
- journal article
- research article
- Published by Copernicus GmbH in Geoscientific Instrumentation, Methods and Data Systems
- Vol. 10 (2), 297-312
- https://doi.org/10.5194/gi-10-297-2021
Abstract
Historic measurements are often temporally incomplete and may contain longer periods of missing data, whereas climatological analyses require continuous measurement records. This is also valid for historic manual snow depth (HS) measurement time series, for which even whole winters can be missing in a station record, and suitable methods have to be found to reconstruct the missing data. Daily in situ HS data from 126 nivo-meteorological stations in Switzerland in an altitudinal range of 230 to 2536 m above sea level are used to compare six different methods for reconstructing long gaps in manual HS time series by performing a “leave-one-winter-out” cross-validation in 21 winters at 33 evaluation stations. Synthetic gaps of one winter length are filled with bias-corrected data from the best-correlated neighboring station (BSC), inverse distance-weighted (IDW) spatial interpolation, a weighted normal ratio (WNR) method, elastic net (ENET) regression, random forest (RF) regression and a temperature index snow model (SM). Methods that use neighboring station data are tested in two station networks with different density. The ENET, RF, SM and WNR methods are able to reconstruct missing data with a coefficient of determination (r2) above 0.8 regardless of the two station networks used. The median root mean square error (RMSE) in the filled winters is below 5 cm for all methods. The two annual climate indicators, average snow depth in a winter (HSavg) and maximum snow depth in a winter (HSmax), can be reproduced by ENET, RF, SM and WNR well, with r2 above 0.85 in both station networks. For the inter-station approaches, scores for the number of snow days with HS>1 cm (dHS1) are clearly weaker and, except for BCS, positively biased with RMSE of 18–33 d. SM reveals the best performance with r2 of 0.93 and RMSE of 15 d for dHS1. Snow depth seems to be a relatively good-natured parameter when it comes to gap filling of HS data with neighboring stations in a climatological use case. However, when station networks get sparse and if the focus is set on dHS1, temperature index snow models can serve as a suitable alternative to classic inter-station gap filling approaches.Keywords
This publication has 45 references indexed in Scilit:
- 21st century Wisconsin snow projections based on an operational snow model driven by statistically downscaled climate dataInternational Journal of Climatology, 2011
- Long-term changes in annual maximum snow depth and snowfall in Switzerland based on extreme value statisticsClimatic Change, 2011
- Reconstructing missing daily precipitation data using regression trees and artificial neural networks for SWAT streamflow simulationJournal of Hydrology, 2010
- Regime shift of snow days in SwitzerlandGeophysical Research Letters, 2008
- Regularization and Variable Selection Via the Elastic NetJournal of the Royal Statistical Society Series B: Statistical Methodology, 2005
- Long-term variability in reconstructed and observed snow cover over the last 100 winter seasons in Cracow and Zakopane (southern Poland)Climate Research, 2002
- Spatial interpolation of climatic Normals: test of a new method in the Canadian boreal forestAgricultural and Forest Meteorology, 1998
- A Three-Way Model for Interpolating for Monthly Precipitation ValuesMonthly Weather Review, 1992
- Estimating Missing Daily Maximum and Minimum TemperaturesJournal of Climate and Applied Meteorology, 1983
- Ridge Regression: Biased Estimation for Nonorthogonal ProblemsTechnometrics, 1970