Introduction
The hard work of removing roof snow has been regarded as a part of life in the heavy snow regions of Japan, and people incur labour and expense for snow removal. During the development of various apparatus for removing roof snow, it became necessary to express snow weight and its daily variance as accurately as possible in order to quantify roof snow to be removed during the winter season.
At present, we can obtain long-term meteorological data on air temperature, precipitation, wind speed and direction, sunlight hours, and snow depth with AMeDAS, which automatically observes normal meteorological conditions at many points over a broad area. We have assumed that estimates of the weight of the snow cover can be obtained using only AMeDAS.
We use a computing method which depends on the difference between the daily quantities of precipitation and melting, in spite of the difficulties in estimating daily snowmelt.
In order to estimate glacier ablation and snowmelt runoff to rivers or streams, methods based on either the energy balance of the snow cover or degree days have been proposed (Reference Kubota, Fukami, Ohmae, Kaneda and YamadaKubota and others, 1978; Reference Takahashi, Sato and NaruseTakahashi and others, 1981; Reference Kojima, Motoyama and YamadaKojima and others, 1983; Reference Ishikawa, Kojima and MotoyamaIshikawa and others, 1985). Heat-balance methods may be accurate, but detailed data over a broad area are difficult to obtain. Functional relationships may estimate long-term snowmelt but can produce unsatisfactory short-term results. An estimate of snowmelt using only degree days is generally improved by increasing the integration period. The suitable integration period depends upon the desired confidence in the estimate of the daily snowmelt variable.
Methods
To estimate the weight of the snow cover using only AMeDAS data, the correlation between long-term measurements for the weight of the snow cover and meteorological factors taken from AMeDAS data observed at the same place must be obtained by statistical methods. The statistical correlation was obtained by the trial and error method using values of the weight of the ground snow cover and meteorological data from 11 winter seasons (1977–87) observed at the Ministry of Agriculture, Forestry and Fishery’s Forestry and Forest Product Research Institute, Tokamachi Experiment Station, in Niigata Prefecture (Tokamachi Experimental Site and Disaster Prevention Forest Laboratory, 1986, 1987).
The total weight of snow cover was obtained by summing the difference between the daily quantities of snowmelt and precipitation. This was computed from the beginning of snowfall each year for observations over a period from 0900 h of the ith day to 0900 h of the (i + 1)th day. We determined the daily quantities of precipitation and melting by the following methods.
First, since precipitation in the air temperature range of 0°–2°C may assume the forms of snow, sleet, hail or rain, it was necessary to evaluate the quantity of precipitation according to its form (Reference TamuraTamura, 1990; Reference HasemiHasemi, 1991). Using the data of Tokamachi Experiment Station during the 1978/79–1984/85 winter seasons, we assumed that the precipitation may be regarded as snow below 0°C and as rain above 2°C. Furthermore, the ratio of snowfall was assumed to decrease linearly with increasing air temperature in the range of 0°–2°C. Consequently, the precipitation quantity of snow at these temperatures was computed for the ith day from
where Tai is daily mean air temperature and Wfi ′ is total daily precipitation quantity at 0°–2°C.
Next, the daily quantities of snowmelt were assumed to be the difference between the total weight of snow cover measured at 0900–1000 h on two sequential days without snowfall or rainfall. Therefore, the relationship between the daily weight reduction of the snow cover and other meteorological factors may be considered statistically. However, it was assumed that the daily precipitation quantity due to rainfall drained away from the snow cover within a day, and that snowmelt did not occur on the day of snowfall. Further, the previous and following days on which snow depth was zero were not included in the data for obtaining daily snowmelt because of the uncertainty of disappearance time, and rain days with precipitation below 5 mm d−1 were adopted to compensate for insufficient data. The timing of the actual runoff from the snow cover is usually delayed by the depth and the stratigraphy of snow cover. Though the daily quantity of snowmelt is not strictly the same as daily quantity of runoff, these values were treated the same for the sake of convenience.
Snowmelt is linked to the energy balance of the snow cover as determined by the various energy fluxes. However, in the case of long-term snowmelt, the quantity of snowmelt may be approximately expressed as a function of daily mean air temperature.
On the basis of the above results, attempts were made to find a functional formula for snowmelt quantity which has significant statistical correlation by using the trial and error method, and incorporating daily mean air temperature, daily maximum air temperature, daily minimum air temperature, daily mean wind speed and daily sunlight hours, etc. As a result, the relationship between the three-day moving accumulation of snowmelt and the three-day moving degree days (moving accumulated positive daily mean air temperature for three days) was obtained with a correlation coefficient of 0.797.
Figure 1 shows the relationship between the three-day moving accumulation of snowmelt and the three-day moving degree days.
The formula obtained approximately from these data is
where wmj is the daily quantity of melt on the ith day. Accordingly, the daily quantity of melt is computed according to
and the weight of the snow cover is expressed for 0 ≤ Tai ≤ 2°C as
where Wsn is the weight of the snow cover on the nth day. When Tai ≤ 0°C
Results and Conclusion
By comparing the computed values using Equations (4) and (5) to the measured weights of the snow cover, good approximations could be obtained throughout the winter season.
Examples of the comparison of measured values with estimated results of the weight of the snow cover for a light snow season and a heavy snow season are shown in Figure 2. Further, Figure 3 shows the relationship between the measured values and the estimated results for daily data through the 11 winter seasons.
The coefficient of correlation among the groups of each value for 1323 days during the 11 winter seasons was 0.986 with a standard deviation of ±590 N m−2. Thus, if the data on daily mean air temperatures and daily precipitation quantities can be obtained for an area, the daily weight of ground snow cover can be estimated with confidence.
The weight of snow on a roof may be estimated from data for the weight of snow cover on the ground. Snow amounts were measured on roofs and on the ground of nearby residences, comparing Aomori City in Aomori Prefecture and Tokamachi City in Niigata Prefecture for several winter seasons (December 1978–March 1991). As a result, it was found that roof snow was 60–80% of that on the ground even during snowfall periods, possibly due to partial snowmelt on roofs because of heat loss through the roof from the residence (Reference Shinojima and HaradaShinojima and Harada, 1992). Although these results are preliminary, we conclude that there is a good possibility of estimating the daily total weight of the snow cover where data exist for daily precipitation and daily mean air temperature.
Acknowledgement
We wish to express appreciation to the Ministry of Agriculture, Forestry and Fishery’s Forestry and Forest Product Research Institute, Tokamachi Experiment Station for permission to use their data.