Please wait a minute...
Frontiers of Earth Science

ISSN 2095-0195

ISSN 2095-0209(Online)

CN 11-5982/P

Postal Subscription Code 80-963

2018 Impact Factor: 1.205

Front. Earth Sci.    2018, Vol. 12 Issue (4) : 683-692    https://doi.org/10.1007/s11707-018-0714-z
RESEARCH ARTICLE
Distribution of snow depth variability
S.R. FASSNACHT1,2,3,4(), K.S.J. BROWN1, E.J. BLUMBERG5, J.I. LÓPEZ MORENO6, T.P. COVINO1,3, M. KAPPAS4, Y. HUANG7, V. LEONE8, A.H. KASHIPAZHA1
1. ESS-Watershed Science, Colorado State University, Fort Collins, CO 80523-1476, USA
2. Cooperative Institute for Research in the Atmosphere, Fort Collins, CO 80523-1375, USA
3. Natural Resources Ecology Laboratory, Fort Collins, CO 80523-1499, USA
4. Geographisches Institut, Georg-August-Universität Göttingen, 37077 Göttingen, Deutschland
5. Geosciences, Colorado State University, Fort Collins, CO 80523-1482, USA
6. Instituto Pirenaico de Ecología, CSIC, Campus de Aula Dei, 50080 Zaragoza, España
7. State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China
8. Seconda Università degli Studi di Napoli, Viale Abramo Lincoln, 5, 81100 Caserta CE, Italia
 Download: PDF(1061 KB)   HTML
 Export: BibTeX | EndNote | Reference Manager | ProCite | RefWorks
Abstract

Snow depth is the easiest snowpack property to measure in the field and is used to estimate the distribution of snow for quantifying snow storage. Often the mean of three snow depth measurements is used to represent snow depth at a location. This location is used as a proxy for an area, typically a digital elevation model (DEM) or remotely sensed pixel. Here, 11, 17, or 21 snow depth measurements were used to represent the mean snow depth of a 30-m DEM pixel. Using the center snow depth measurement for each sampling set was not adequate to represent the pixel mean, and while the use of three snow depth measurements improved the estimate of mean, there is still large error for some pixels. These measurements were then used to determine the variability of snow depth across a pixel. Estimating variability from few points rather than all in a measurement was not sufficient. The sampling size was increased from one to the total per pixel (11, 17, or 21) to determine how many point samples were necessary to approximate the mean snow depth per pixel within five percent. Binary regression trees were constructed to determine which terrain and canopy variables dictated the spatial distribution of the snow depth, the standard deviation of snow depth, and the sample size to within 5% of the mean per pixel. One location was measured in two years just prior to peak accumulation, and it is shown that there was little to no inter-annual consistency in the mean or standard deviation.

Keywords uncertainty      sampling      binary regression trees      snow telemetry     
Corresponding Author(s): S.R. FASSNACHT   
Just Accepted Date: 14 May 2018   Online First Date: 21 June 2018    Issue Date: 20 November 2018
 Cite this article:   
S.R. FASSNACHT,K.S.J. BROWN,E.J. BLUMBERG, et al. Distribution of snow depth variability[J]. Front. Earth Sci., 2018, 12(4): 683-692.
 URL:  
https://academic.hep.com.cn/fesci/EN/10.1007/s11707-018-0714-z
https://academic.hep.com.cn/fesci/EN/Y2018/V12/I4/683
Fig.1  (a) Location map for the i. Togwotee Pass and ii. Joe Wright SNOTEL stations in Wyoming and Colorado, respectively. (b) Canopy density in percent for the area surveyed: i. Togwotee Pass was surveyed in 2009, and ii. Joe Wright was surveyed in 2009 (red symbols) and iii. 2010 (black symbols). (c) The snow depth sampling design at each sampling plot is i) 21 points for Togwotee Pass (2009), ii. 11 points for Joe Wright 2009, and iii. 17 points for Joe Wright 2010.
Fig.2  Comparison of (a) mean and (b) standard deviation comparison between 3 point statistics (grey diamonds) and plot statistics for three snow depth surveys about SNOTEL stations: i. Togwotee Pass in March 2009, ii. Joe Wright in May 2009, and iii. Joe Wright in May 2010. The center or 1 point snow depth is also compared to the pixel mean (open box). The maximum number of points per station-date are in parentheses.
Fig.3  Mean absolute difference from the pixel mean snow depth for all samples per pixel as a function of the number of points per pixel for the three sampling dates. The number in the (a) parentheses is the number of samples per pixel, while (b) is for only the 11 points in a row (as per Figure 1cii). The 5% difference threshold is shown as the dotted line.
Comparison Togwotee Pass 2009 Joe Wright 2009 Joe Wright 2010
NCSE R2 slope NCSE R2 slope NCSE R2 slope
1 to 11 0.61 0.73 1.02 0.62 0.78 1.14 0.44 0.62 0.96
1 to 17 N/A N/A 0.37 0.61 0.98
1 to 21 0.58 0.72 1.02 N/A N/A
3 to 11 0.90 0.90 0.93 0.93 0.94 0.94 0.89 0.89 0.96
3 to 17 N/A N/A 0.88 0.89 1.0
3 to 21 0.83 0.84 0.91 N/A N/A
5 to 21 0.90 0.91 0.92 N/A N/A
11 to 17 N/A N/A 0.98 0.94 0.98
11 to 21 0.94 0.94 0.98 N/A N/A
Tab.1  Statistics mean summarizing the comparison of fewer sampling points (1, 3, or 5) to all sampling points (11, 17, or 21)
Comparison Togwotee Pass 2009 Joe Wright 2009 Joe Wright 2010
NCSE R2 Slope NCSE R2 Slope NCSE R2 Slope
3 to 11 0.22 0.55 0.96 0.19 0.58 1.05 –0.20 0.39 0.87
3 to 17 N/A N/A –0.30 0.40 0.92
3 to 21 0.03 0.34 0.65 N/A N/A
5 to 21 0.46 0.64 0.95 N/A N/A
11 to 17 N/A N/A 0.76 0.79 0.94
11 to 21 0.35 0.46 0.58 N/A N/A
Tab.2  Statistics standard deviation summarizing the comparison of fewer sampling points (1, 3, or 5) to all sampling points (11, 17, or 21)
Fig.4  Joe Wright SNOTEL inter-annual (May 2009 and 2010) comparison of (a) mean snow depth, and (b) standard deviation for data at the 99 pixels within the same pixel (grey diamonds) and 70 sets of pixels that were within one pixel. A pixel is derived from the 30-meter U.S. Geological Survey (2018) digital elevation model<nationalmap.usgs.gov>. The Nash-Sutcliffe coefficient of efficiency statistic is not included as it is less than zero in all cases.
Fig.5  Togwotee Pass 2009 binary regression trees for depth to 5% of the mean per pixel. Boxes represent decision nodes with the true being the node to the left below. The value below the variable represents the mean value of all pixels remaining. Circles represent terminal nodes with the mean of the number of points to within 5% of snow depth, for the number of points (n) in that category.
Variable (a) Measurements to five 5% of mean per pixel (b) Mean snow depth per pixel (c) Standard deviation per pixel
Togwotee
Pass
2009
Joe Wright 2009 Joe Wright 2010 Togwotee
Pass
2009
Joe Wright 2009 Joe Wright 2010 Togwotee
Pass
2009
Joe Wright 2009 Joe Wright 2010
Canopy density N/A 2,4 4 3 1 3,5 1 2,5 1
Latitude N/A 2,3 4,5 N/A 6 N/A N/A 4,5,5 4,5,5
Longitude 4,4 N/A 6,7 4 N/A 1,3,4,5 4,4 N/A 5
Elevation 1,4,5 1,5,6 1,2,8 3,5,9 2,3,6 2,4,6 3,4 N/A 4
Sine of slope 3,3 N/A 5 5,6,7,9,10 3 N/A 2,3,5 1,3 2,3.4
Max upwind slope 2,3 4 7,8 8 N/A N/A 3,6 6 3,5
Eastness 3 3,7,8 3 1 4,5,6 5,5 2,5 5,7 3
Northness 4 3,4 N/A 4 N/A N/A N/A 3,8 2
Solar radiation 2 N/A N/A 2 N/A 4 4 4 N/A
Relative standard error 0.0321 0.0121 0.0229 0.0354 0.0283 0.0184 0.0708 0.0409 0.0182
Tab.3  Summary of the order of appearance in the binary regression tree for the (a) points to within 5% of the pixel mean from top (1) to last node (e.g., Fig. 2), (b) mean, and (c) standard deviation, and. The same number in the order list identifies that variable occurring more than once in the particular level. The fit statistic is the relative?standard errors using a minimum leaf size of 5% of total number of samples
1 Bales R C, Molotch N P, Painter T H, Dettinger M, Rice R, Dozier J (2006). Mountain hydrology of the western United States. Water Resour Res, 42(8): W08432
https://doi.org/10.1029/2005WR004387
2 Balk B, Elder K (2000). Combining binary decision tree and geostatistical methods to estimate snow distribution in a mountain watershed. Water Resour Res, 36(1): 13–26
https://doi.org/10.1029/1999WR900251
3 Blöschl G (1999). Scaling issues in snow hydrology. Hydrol Processes, 13(14–15): 2149–2175
https://doi.org/10.1002/(SICI)1099-1085(199910)13:14/15<2149::AID-HYP847>3.0.CO;2-8
4 Debarbieux B, Price M F (2012). Mountain regions: a global common good? Mt Res Dev, 32(S1): S7–S11
https://doi.org/10.1659/MRD-JOURNAL-D-11-00034.S1
5 Deems J S, Fassnacht S R, Elder K J (2006). Fractal distribution of snow depth from LiDAR data. J Hydrometeorol, 7(2): 285–297
https://doi.org/10.1175/JHM487.1
6 Deems J S, Fassnacht S R, Elder K J (2008). Interannual consistency in fractal snow depth patterns at two Colorado mountain sites. J Hydrometeorol, 9(5): 977–988
https://doi.org/10.1175/2008JHM901.1
7 Dingman S, Seely-Reynolds D, Reynold D (1988). Application of Kriging to estimating mean annual precipitation in a region of orographic influence. J Am Water Resour Assoc, 24(2): 329–339
https://doi.org/10.1111/j.1752-1688.1988.tb02991.x
8 Dozier J, Frew J (1990). Rapid calculation of terrain parameters for radiation modeling from digital elevation data. IEEE Trans Geosci Remote Sens, 28(5): 963–969
https://doi.org/10.1109/36.58986
9 Elder K, Dozier J, Michaelsen J (1991). Snow accumulation and distribution in an alpine watershed. Water Resour Res, 27(7): 1541–1552
https://doi.org/10.1029/91WR00506
10 Elder K J, Cline D W, Liston G E, Armstrong R (2009). NASA cold land processes experiment (CLPX 2002–2003): field measurements of snowpack properties and soil moisture. J Hydrometeorol, 10(1): 320–329
https://doi.org/10.1175/2008JHM877.1
11 Erickson T, Williams M, Winstral A (2005). Persistence of topographic controls on the spatial distribution of snow in rugged mountain terrain, Colorado, United States. Water Resour Res, 41(4): W04014
https://doi.org/10.1029/2003WR002973
12 Erxleben J, Elder K, Davis R (2002). Comparison of spatial interpolation methods for estimating snow distribution in the Colorado Rocky Mountains. Hydrol Processes, 16(18): 3627–3649
https://doi.org/10.1002/hyp.1239
13 Fassnacht S R, Deems J S (2006). Measurement sampling and scaling for deep montane snow depth data. Hydrol Processes, 20(4): 829–838
https://doi.org/10.1002/hyp.6119
45 Fassnacht S R, Deitemeyer D C, Venable N B H (2014). Capitaliziing on the daily time step of snow telemetry data to model the snowmelt components of the hydrograph for small watersheds. Hydrological Processes, 28(16): 4654–4668
14 Fassnacht S R, Dressler K A, Bales R C (2003). Snow water equivalent interpolation for the Colorado River Basin from snow telemetry (SNOTEL) data. Water Resour Res, 39(8): 1208
https://doi.org/10.1029/2002WR001512
15 Fassnacht S R, Dressler K A, Hultstrand D M, Bales R C, Patterson G G (2012). Temporal inconsistencies in coarse-scale snow water equivalent patterns: Colorado River Basin snow telemetry-topography regressions. Pirineos, 167(0): 165–185
https://doi.org/10.3989/Pirineos.2012.167008
16 Fassnacht S R, Heun C, López Moreno J I, Latron J (2010). Variability of snow density measurements in the Esera valley, Pyrenees mountains, Spain. Cuadernos de Investigación Geográfica, 36(1): 59–72
https://doi.org/10.18172/cig.1227
17 Fassnacht S R, Williams M W, Corrao M V (2009). Changes in the surface roughness of snow from millimetre to metre scales. Ecol Complex, 6(3): 221–229
https://doi.org/10.1016/j.ecocom.2009.05.003
18 Fassnacht S R, Wyss D, Heering S M (2017). Räumliches Denken – ein Beispiel im Schnee / A Spatial Thinking Research-Didactic Example in Snow. Geoöko, 37(3-4)
19 Hultstrand D, Fassnacht S R, Stednick J (2006). Geostatistical methods for estimating snowmelt contribution to the alpine water balance. Proceedings of the Annual Western Snow Conference, (Las Cruces NM), 74: 149–154
20 Kashipazha A H (2012). Practical snow depth sampling around six snow telemetry (SNOTEL) stations in Colorado and Wyoming, United States. Dissertation for M.S. degree. Watershed Science Program, Colorado State University, Fort Collins, Colorado USA
21 Logan L (1973). Basin-wide water equivalent estimation from snowpack depth measurements. In: Role of Snow and Ice in Hydrology. IAHS AIHS Publ. 107: 864–884
22 López–Moreno J I, Fassnacht S R, Begueria S, Latron J (2011). Variability of snow depth at the plot scale: implications for mean depth estimation and sampling strategies. Cryosphere, 5(3): 617–629
https://doi.org/10.5194/tc-5-617-2011
23 López-Moreno J I, Fassnacht S R, Heath J T, Musselman K, Revuelto J, Latron J, Morán-Tejeda E, Jonas T (2013). Small scale spatial variability of snow density and depth over complex alpine terrain: implications for estimating snow water equivalent. Adv Water Resour, 55: 40–52
https://doi.org/10.1016/j.advwatres.2012.08.010
24 López-Moreno J I, Nogués Bravo D (2006). Interpolating local snow depth data: an evaluation of methods. Hydrol Processes, 20(10): 2217–2232
https://doi.org/10.1002/hyp.6199
25 López-Moreno J I, Revuelto J, Fassnacht S R, Azorín-Molina C, Vicente-Serrano S M, Morán-Tejeda E, Sexstone G A (2015). Snowpack variability across various spatio-temporal resolutions. Hydrol Processes, 29(6): 1213–1224
https://doi.org/10.1002/hyp.10245
26 Meiman J R (1968). Snow accumulation related to elevation, aspect, and forest canopy. In Proceedings of Snow Hydrology Workshop Seminar, Canadian National Committee of the International Hydrological Decade, Fredericton NB, pp. 35–47
27 Meromy L, Molotch N, Link T, Fassnacht S R, Rice R (2013). Subgrid variability of snow water equivalent at operational snow stations in the western USA. Hydrol Processes, 27(17): 2383–2400
https://doi.org/10.1002/hyp.9355
28 Molotch N, Brooks P, Burns S, Litvak M, Monson R, McConnell J, Musselman K (2009). Ecohydrological controls on snowmelt partitioning in mixed conifer sub alpine forests. Ecohydrology, 2(2): 129–142
https://doi.org/10.1002/eco.48
29 Molotch N P, Bales R C (2005). Scaling snow observations from the point to the grid-element: implications for observation network design. Water Resour Res, 41: W11421
https://doi.org/https://doi.org/10.1029/2005WR004229
30 Moriasi D N, Arnold J G, Van Liew M W, Bingner R L, Harmel R D, Veith T L (2007). Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans ASABE, 50(3): 885–900
https://doi.org/10.13031/2013.23153
31 Mould B (1985). New Day Rising. Chapter 1 in New Day Rising, Solid State Tuners, Taylor, Texas
32 Nash J E, Sutcliffe J V (1970). River flow forecasting through conceptual models part I- A discussion of principles. J Hydrol (Amst), 10(3): 282–290
https://doi.org/10.1016/0022-1694(70)90255-6
33 Neumann N, Derksen C, Smith C, Goodison B (2006). Characterizing local scale snow cover using point measurements during the winter season. Atmos-ocean, 44(3): 257–269
https://doi.org/10.3137/ao.440304
34 Rice R, Bales R (2010). Embedded sensor network design for snow cover measurements around snow pillow and snow course sites in the Sierra Nevada of California. Water Resour Res, 46(3): W03537
https://doi.org/10.1029/2008WR007318
35 Sexstone G A, Fassnacht S R (2014). What drives basin scale spatial variability of snowpack properties in the Front Range of northern Colorado? Cryosphere, 8(2): 329–344
https://doi.org/10.5194/tc-8-329-2014
36 Sturm M, Wagner A (2010). Using repeated patterns in snow distribution modeling: an Arctic example. Water Resour Res, 46(12): W12549
https://doi.org/10.1029/2010WR009434
37 Trujillo E, Ramírez J A, Elder K J (2007). Topographic, meteorologic, and canopy controls on the scaling characteristics of the spatial distribution of snow depth fields. Water Resour Res, 43(7): W07409
https://doi.org/10.1029/2006WR005317
46 U.S. Department of Agriculture (2011). Snow Survey and Water Supply Forecasting Program. National Engineering Handbook Part 622, Water and Climate Center, Natural Resources Conservation Service
38 U.S. Department of Agriculture (2016). Snow Survey and Water Supply Forecasting Program. National Water and Climate Center, Natural Resources Conservation Service
39 U.S. Geological Survey (2018). The National Map. U.S. Department of Interior
40 Veatch W, Brooks P, Gustafson J, Molotch N (2009). Quantifying the effects of forest canopy cover on net snow accumulation at a continental, mid latitude site. Ecohydrology, 2(2): 115–128
https://doi.org/10.1002/eco.45
41 Wallace C, Gass L (2008). Elevation derivatives for Mojave Desert tortoise habitat models. US Geological Survey Open File Report 2008-1283, 7pp
42 Webb R (2017). Using ground penetrating radar to assess the variability of snow water equivalent and melt in a mixed canopy forest, Northern Colorado. Front Earth Sci, 11(3): 482–495
https://doi.org/10.1007/s11707-017-0645-0
43 Winstral A, Elder K, Davis R (2002). Spatial snow modeling of wind redistributed snow using terrain based parameters. J Hydrometeorol, 3(5): 524–538
https://doi.org/10.1175/1525-7541(2002)003<0524:SSMOWR>2.0.CO;2
44 Winstral A, Marks D (2002). Simulating wind fields and snow redistribution using terrain based parameters to model snow accumulation and melt over a semi arid mountain catchment. Hydrol Processes, 16(18): 3585–3603
https://doi.org/10.1002/hyp.1238
[1] D.M. HULTSTRAND, S.R. FASSNACHT. The sensitivity of snowpack sublimation estimates to instrument and measurement uncertainty perturbed in a Monte Carlo framework[J]. Front. Earth Sci., 2018, 12(4): 728-738.
[2] Igor Appel. Uncertainty in satellite remote sensing of snow fraction for water resources management[J]. Front. Earth Sci., 2018, 12(4): 711-727.
[3] Grey S. NEARING, Hoshin V. GUPTA. Ensembles vs. information theory: supporting science under uncertainty[J]. Front. Earth Sci., 2018, 12(4): 653-660.
[4] Huitao SHEN,Wanjun ZHANG,Jiansheng CAO,Xiang ZHANG,Quanhong XU,Xue YANG,Dengpan XIAO,Yanxia ZHAO. Carbon concentrations of components of trees in 10-year-old Populus davidiana stands within the Desertification Combating Program of Northern China[J]. Front. Earth Sci., 2016, 10(4): 662-668.
[5] Wei LI,Changchun FENG,Chao DAI,Yongping LI,Chunhui LI,Ming LIU. An inexact risk management model for agricultural land-use planning under water shortage[J]. Front. Earth Sci., 2016, 10(3): 419-431.
[6] C. DAI,W. SUN,Q. TAN,Y. LIU,W.T. LU,H.C. GUO. Risk management for sulfur dioxide abatement under multiple uncertainties[J]. Front. Earth Sci., 2016, 10(1): 87-107.
[7] Yue LIU,Qiuming CHENG,Qinglin XIA,Xinqing WANG. The use of evidential belief functions for mineral potential mapping in the Nanling belt, South China[J]. Front. Earth Sci., 2015, 9(2): 342-354.
[8] Ran WANG,Yin LI,Qian TAN. A review of inexact optimization modeling and its application to integrated water resources management[J]. Front. Earth Sci., 2015, 9(1): 51-64.
[9] Fangbing MA,Chunhui LI,Xuan WANG,Zhifeng YANG,Chengchun SUN,Peiyu LIANG. A Bayesian method for comprehensive water quality evaluation of the Danjiangkou Reservoir water source area, for the middle route of the South-to-North Water Diversion Project in China[J]. Front. Earth Sci., 2014, 8(2): 242-250.
[10] Cong DONG, Guohe HUANG, Qian TAN, Yanpeng CAI. Coupled planning of water resources and agricultural land-use based on an inexact-stochastic programming model[J]. Front Earth Sci, 2014, 8(1): 70-80.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed