Abstract It has been recognized that the snowmelt models developed in the past do not fully meet current prediction requirements. Part of the reason is that they do not account for the spatial variation in the dynamics of the spatially heterogeneous snowmelt process. Most of the current physics-based distributed snowmelt models utilize point-location-scale conservation equations which do not represent the spatially varying snowmelt dynamics over a grid area that surrounds a computational node. In this study, to account for the spatial heterogeneity of the snowmelt dynamics, areally averaged mass and energy conservation equations for the snowmelt process are developed. As a first step, energy and mass conservation equations that govern the snowmelt dynamics at a point location are averaged over the snowpack depth, resulting in depth averaged equations (DAE). In this averaging, it is assumed that the snowpack has two layers. Then, the point location DAE are averaged over the snowcover area. To develop the areally averaged equations of the snowmelt physics, we make the fundamental assumption that snowmelt process is spatially ergodic. The snow temperature and the snow density are considered as the stochastic variables. The areally averaged snowmelt equations are obtained in terms of their corresponding ensemble averages. Only the first two moments are considered. A numerical solution scheme (Runge-Kutta) is then applied to solve the resulting system of ordinary differential equations. This equation system is solved for the areal mean and areal variance of snow temperature and of snow density, for the areal mean of snowmelt, and for the areal covariance of snow temperature and snow density. The developed model is tested using Scott Valley (Siskiyou County, California) snowmelt and meteorological data. The performance of the model in simulating the observed areally averaged snowmelt is satisfactory.
[1]
Stephen F. Ackley,et al.
Snow and ice
,
1975
.
[2]
E. Anderson,et al.
A point energy and mass balance model of a snow cover
,
1975
.
[3]
W. Swinbank.
A comparison between predictions of dimensional analysis for the constant‐flux layer and observations in unstable conditions
,
1968
.
[4]
W. Tollmien,et al.
Meteorologische Anwendung der Strömungslehre
,
1961
.
[5]
T. Dunne,et al.
The generation of runoff from subarctic snowpacks
,
1976
.
[6]
Colleagues,et al.
Physical Studies on Deposited Snow. I. ; Thermal Properties.
,
1955
.
[7]
E. F. Bradley,et al.
Flux-Profile Relationships in the Atmospheric Surface Layer
,
1971
.
[8]
D. H. Male.
Snowcover ablation and runoff
,
1981
.
[9]
K. Kojima.
Densification of Seasonal Snow Cover
,
1967
.
[10]
J. L. Smith,et al.
Advanced concepts and techniques in the study of snow and ice resources. An interdisciplinary symposium held at Monterey, California, December 2--6, 1973
,
1974
.
[11]
Engineering Properties of Snow
,
1977
.
[12]
J. Kondo,et al.
A Prediction Model for Snowmelt, Snow Surface Temperature and Freezing Depth Using a Heat Balance Method
,
1990
.
[13]
E. Morris,et al.
Modelling the Flow of Mass and Energy within a Snowpack for Hydrological Forecasting
,
1983,
Annals of Glaciology.
[14]
Roger G. Barry,et al.
Application of Computed Global Radiation for Areas of High Relief
,
1972
.
[15]
Y. Yen.
Effective thermal conductivity of ventilated snow
,
1962
.