Spatial variability in degree-day factorsin Yarlung Zangbo River Basin in China*

LIU Jinping, ZHANG Wanchang

(1Key Laboratory of Digital Earth Science,Institute of Remote Sensing and Digital Earth,Chinese Academy ofSciences,Beijing100094,China; 2University of Chinese Academy of Sciences,Beijing100049,China) (Received 22 June 2017; Revised 25 October 2017)

Abstract Degree-day model, based on a linear relationship between glacier/snow ablation and positive accumulated temperature, is widely used for estimating glacier and snow melt for its simplicity and outstanding performance. However, degree-day factor (DDF), as the most important parameter in this model, is difficult to obtain in an extensive region, especially in Yarlung Zangbo River Basin, which is the highest basin in the world and characterized by rugged and complex terrain. With the development of distributed hydrological model, increasing demand for spatial dataset of DDFs has become a challenging issue in the relevant fields. The objective of this study is to develop an algorithm to spatially estimate the positive DDFs based on available observed DDFs. Topographical factors are fully taken into account to correlate with DDFs, and spatial DDFs in moderate resolution are generated by stepwise and neighboring mean filtering improved in this study. Finally, the produced DDFs are indirectly validated by snow lines in different orientations, and the correlation coefficient of up to 0.82 proves that the DDFs generated in the present study are reasonable in applications. The spatial distribution of DDFs in Yarlung Zangbo River Basin shows that DDFs in downstream were lower than those in upstream and midstream, which is probably because there is less precipitation and more solar radiation in upstream and midstream than in downstream, due to monsoon from India Ocean. To better estimate glacial meltwater runoff and manage water resources in Yarlung Zangbo River Basin, further studies on spatial variability in DDFs are yet needed.

Keywords glacier; degree-day factor; Yarlung Zangbo River Basin; topographical factors; spatial variability

In the context of global warming, a lot of the world’s small mountain glaciers are losing mass at an accelerated rate, which causes a rise in sea level[1]. Much effort has been made to develop various methods to estimate the amount of glacier ablation throughout the world[2]. In these methods, degree-day model was considered to be the most common approach for melt modelling due to its simplicity and good performance[3]. The model depends on an assumed relationship between ablation and air temperature, which is expressed, in general, in the form of positive mean daily temperature sum. Degree-day factor (hereafter, DDF), being the factor of proportionality, represents the glacier ablation for a unit positive degree-day[4]. It simplifies the complex physical processes that are more properly described by energy balance models, which means that the factor itself has a certain physical mechanism. Accordingly, there is various in DDF caused by the energy partitioning with different seasons and spatial locations[5]. Several studies have proved that the seasonal variation in DDFs of glaciers can be ignored for its no significant[3,6]. However, a large difference of DDFs in the spatial distribution was revealed[7], which greatly affect the accuracy of glaciers melt model[3].

In general, the DDF can be derived from measurement of ablation stick[8], as well as calculation by energy balance model[7]. In the Yarlung Zangbo River Basin (hereafter, YZRB), characterized by high altitude and extreme topographical conditions[9], direct and continuous measurements of DDFs are not available due to the inaccessibility of glaciers and the logistical difficulties associated with maintaining a DDF record. Energy balance model is characterized by complex physical processes, which results difficulties in application due to the restricted data availability in this basin. In order to improve the accuracy of runoff yield in the glacierized basin for distributed hydrologic model, a proper algorithm need to be developed to generate a distributed DDF map of glaciers scale in YZRB on moderate resolution (1 km).

With different topography (i.e. elevation, slope, aspect, latitude and longitude) and heat fluxes, DDFs will vary considerably in spatial distribution[10]. In general, the DDFs are considered to be positive correlation with elevation and direct solar radiation input, but negative correlation with sensible heat flux[3]. Pre-existing study suggested that high shares of sensible heat flux, resulting from high air temperature and wind speed, usually associated with low DDFs[11]. It is well known that air temperature has a strong relationship with topography. Therefore, a multiple regression equation with statistical optimum choosing between observed DDFs of different glaciers around the YZRB and topographic factors is established to produce distributed DDFs in YZRB, and the method of improved neighborhood mean filtering is used to revise some abnormal pixels in DDFs map in this study.

1 Study area

Originating from the Chemayungdung Glacier in the northern slope with an average elevation of 4 600 m. a.s.l. over the Himalayas, the River in YZRB is the highest river on the earth stretching across the south part of the Tibetan Plateau from the west to the east with total length of about 2 057 km, a drainage area of about 257 680 km2[9,12]. Geographically, the YZRB is situated between 82°1′-97°6′E and 27°49′-31°17′ N (Fig. 1), encompassing a narrow strip with the east-west length exceeding the north-south width, stretching northward to Gangdisê Mountains and the Nyainqêntanglha Mountains and southward to the Himalayas, where about 41 countries (cities) located. The YZRB is a part of Indian Ocean water system. Because of its special geographic location in the Yarlung Zangbo Suture Zone, contributing to the moisture movement from the Indian Ocean to the inner regions of the Qinghai-Tibet Plateau, the YZRB has a significant effect on the rainfall amount of Qinghai-Tibet Plateau[13]. The annual precipitation has an obvious east-west declining gradient from above 2 000 mm to less than 300 mm for the decaying effect of the monsoon[14]. In this basin, there are 10 816 glaciers that covers about 14 493 km2in area and 1 293 km3in ice volume, accounting for 29% of the total number of glaciers (36 918) on the Qinghai-Tibet Plateau, as well as 23% of the glacial total in China[15].

Fig.1 Fundamental geographic information of Yarlung Zangbo River Basin and the observed DDFs of glaciers around this basin

2 Data utilization

Due to sparse data availability in the west of China, some short-term observed value of DDFs in different glaciers were collected from literatures and measurement obtained by ourselves in the Yanlongba glacier, the DDFs of ten glaciers around the YZRB were calculated, which was summarized in Table 1 and their geographic locations were shown in Fig.1.

Topography information is important to estimate spatial distribution of ice-DDFs in this paper. Shuttle Radar Topography Mission (SRTM) was hosted by U.S. National Aeronautics and Space Administration (NASA), U.S. National Geospatial-Intelligence Agency (NGA), German Aerospace Center (DLR) and Italian Space Agency (ASI). The version 4 of SRTM DEM (digital elevation model) used in this research is seamless. The DEM data, covering the scope of 60°N-56°S, includes three levels in resolutions, at 1 000, 90 and 30 m, respectively. The spatial resolution of 1 000 m, called moderate resolution image from the aspect of remote sensing technology, is adopted in this paper, which involves two points: one is that the specified resolution of DDFs is enough for simulating runoff in YZRB, the other is that it can properly reduce heterogeneity of underlying surface in DDFs of glaciers. With the 1 000 m resolution of DEM data, ArcGIS software is used to calculate elevation, latitude, longitude, aspect and slope, constituting the topography data in the YZRB.

In addition, the datasets of first Glacier Inventory of China (CGI) derived from Cold and Arid Regions Science Data Center (WestDC)[10]was used as the extent of glaciers in YZRB.

Table 1 The average DDFs for different glaciers

3 Methodology

3.1 Definition and calculation of glacier DDF

DDF is defined as the ratio between the amount of glacier ablation and the positive accumulated temperature in the same period[20]. Therefore, the amount of glacier ablation in a specific period can be expressed as

wherekgstandsfortheglacierDDF(mm·d-1·℃-1); Agdenotestheamountofglacierablation(mmw.e.),whichisobtainedbythedifferenceoftwomeasurementsusingsightingrodmultiplybytheglacierdensity(～900kg·m-3)[21];PDDgrepresentsthepositiveaccumulatedtemperatureinacertainperiod(℃).ThePDDgcanbedescribedas

whereTtdenotesthedailymeantemperature(℃); Ht,alogicalvariable,is0whentheTtisnegativevalueandis1whenTtispositivevalue.PDDsstandsforthepositiveaccumulatedtemperatureusedforsnow-melt.ItisimportanttokeepinmindthatPDDsisequal0whenthereisnosnowcoveronthespecificglacier-coverpixel,andthepositiveaccumulatedtemperatureisfirstusedforsnow-meltandthenfortheglacier-melt[16].However,themeasurementsdatainthefieldarethetotalmassbalanceandaccumulation,therefore,thetotalablationofglacierandsnowisobtained,whichleadstothefactthatthePDDsandPDDgcannotbeseparatedinthetotalpositiveaccumulatedtemperature.Therefore,actuallythecalculatedkgbasedonmeasurementsdataistheintegratedDDFofsnowandglacier[16].

3.2 Estimation of the distributed glacier DDFs

The stepwise regression is widely used in climate and environment prediction for its low computational load and effective variable selection[22-23], which was introduced into estimation of spatial ice-DDFs in this study. In an actual application, the dependent variable DDF was selected as a function of parts or all of independent variables that included elevation, latitude, longitude, slope and aspect. The concrete procedure of stepwise regression was entering one variable at a time into the independent variable set (hereafter, x-set) according to the influence degree of independent variable on dependent variable, and the previous variable in x-set would be removed if it became not significant for dependent variable resulting from a new variable entering. This process would end when no available variables satisfying the tolerance criterion remains or when the reduction in the residual sum of squares failed to exceed a specified threshold[24].

According to previous study, the minimum DDF of glaciers in western China was 2.6 mm·d-1·℃-1[10], and maximum value 16.9 mm·d-1·℃-1referring to Table 1. Therefore, the value of DDFs derived by stepwise regression without this range was considered to be abnormal pixels (hereafter, it was defined as outliers). Aiming to correct outliers in DDFs and keep normal value unchanged, the method of improved neighboring mean filtering[25]was applied in this paper. The specific process of the improved neighboring mean filtering was shown in Fig. 2. In this proposed algorithm, only the outliers were processed as

whereODDFstandsforcorrectedresultofthecurrentoutlierbycalculatingthemeanofnormalneighboringDDFs, NDDFthenormalneighboringDDFandnthenumberofnormalpixels.ThepixelwillbeskippedifthenumberofnormalDDFsisequalto0andrecalculatedinnextiteration.Itshouldbenotedthatthe5×5filteringwindowisselectedallowingfortheareaofsingleglacier.

Fig.2 Flowchart of the neighboring mean filtering improved in this work

4 Results

With the method of stepwise regression, variables of longitude, slope and aspect in topographical factors were excluded because they were not significant in regression model. The multiple linear regression model between DDFs and significant topographical factors was shown in Equation (2) (r2=0.8) and it was significant at the level of 0.01. Many studies have already shown that higher elevation, in general, results in higher DDFs[2-3,10]. In addition, Jones[26]argued that there was large difference in melting rate at various altitude and latitude, which indicated that the DDFs have strong relationship with elevation and latitude. These related studies all proved the reasonability of Equation (4).

DDF=0.009×ELE-0.934×LAT-8.1.

(4)

Fig.3 Spatially distributed DDFs of glaciers in YZRB

ThefinalspatialdistributionofDDFsinYZRBwasshowninFig. 3.Actually,inordertogettheresult,therewerethreestepsincluded:thefirststepwastocalculatethepreliminaryDDFsinYZRBusingstepwiseregressionmethod.ThesecondstepwasthattheDDFsinglacialextentwasextractedwiththefirstGlacierInventoryofChina(CGI).ThelaststepwastorevisetheincorrectedDDFsbythealgorithmofimprovedneighboringmeanfiltering.DuetolackofobservedglacialDDFs,directvalidationforthisresultseemedtobedifficult.Fortunately,previousstudyshowedthathighersnowlineusuallymeanslargernegativeglacialmassbalance[15],aswellaslargerglacialDDF.Inordertovalidatethereliabilityofthisresult,themeanDDFandsnowlineindifferentorientationsofglaciers(Table2)werecalculatedbyArcGISsoftwareandderivedfromYaoetal[15],respectively.ThecorrelationcoefficientbetweenDDFsandsnowlinesindifferentorientationswasupto0.82 (p<0.05),whichsuggestedthattheDDFsproducedinthispaperwasreasonablefortheapplicationinhydrology.

FromFig. 3,wecouldfindthatDDFsindownstreamofYZRBwaslowerthanthatinmidstreamandupstream.Oneofpossibleexplanationswasthattherewerelessprecipitationandmoresolarradiationinupstreamandmidstreamcomparedtodownstream.AccordingtotheresultofDDFs,theaveragevalueofglacialDDFsoverYZRBwas10.97mm·d-1· ℃-1,whichcouldbeusedasareferenceforsomeapplicationsofcryosphereandhydrology.

Table 2 Mean DDFs and snow lines

5 Discussion

The DDFs data produced in this study is proven to be reasonable in the different glacier orientations by indirect validation method. However, in addition to the influence of terrain factors presented in this study (i.e. elevation, longitude, latitude, aspect, and slope), the DDF will also be affected by the conditions such as mountain sheltering, wind speed, and solar radiation[16,21], therefore, no doubt errors do exist in the distributed glacier DDFs data, which may also result from few observed samples, simple regression model, spatial scale, and the neglected seasonal fluctuation. Still, it is very significant to study on this in the ungauged basin under the development of the distributed hydrological modeling, but it requires further studies on the accurate estimation of the distributed glacier DDF in future.

A follow-up study will make an attempt at retrieval of the distributed glacier DDFs based on remotely sensed data. The monthly glacier mass change can be computed by GRACE (gravity recovery and climate experiment) data, which is described detailedly in these literatures[27-28]. Then the monthly distributed glacier DDFs can be obtained with the combination of the glacier mass change and positive accumulated temperature. In addition, the distributed snow DDFs also can be computed by snow water equivalent (SWE) change, which was derived from the special sensor microwave/imager (SSM/I) and scanning multi-channel microwave radiometer (SMMR)[29], and the positive accumulated temperature. This is a preliminary idea of using remotely sensing (RS) data to produce the distributed glacier/snow DDFs, therefore, further studies need to be done in future. Also, the RS-based glacier DDFs can provide an opportunity to make a comparison with the distributed glacier DDFs presented in this study.

6 Conclusion

Because of insensitivity to time dimension, DDFs were only considered to change in spatial distribution. In this study, observed glacial DDFs were collected from literatures and measurement around YZRB. Topographical factors, including elevation, longitude, latitude, aspect and slope derived from SRTM DEM, were fully taken into account in spatial variability of DDFs. With the method of stepwise regression and the datasets of CGI, distributed DDFs were produced in the glacial extent of YZRB by regression model, as well as correction of outliers through improved neighboring mean filtering proposed in this study. The result was indirectly validated by snow lines in different orientations due to lack of available glacial DDFs for validation, which indicated that the reliability of produced DDFs was comparable because the correlation coefficient between snow lines and mean DDFs was up to 0.82.

In the YZRB, DDFs in downstream was lower than that in upstream and midstream, one of possible explanations could be that there were drier and more solar radiation in upstream and midstream compared to downstream influenced by monsoon from India Ocean. According to the result of DDFs, the mean value of DDFs in glacial extent of YZRB was 10.97 mm·d-1·℃-1, which was important to cryosphere and hydrology, especially in the region of sparse data availability.

The spatial dataset in DDFs was produced in this paper, and it was very significant in the distributed hydrological model. However, owning to little studies focusing on the variability of DDFs in the western China, especially in high-mountain regions, future studies on validation of this result are still needed.