为了深入分析气温对长江上游年径流的影响和解释青藏高原冰川融水再冻结现象的物理机制,采用对位置、尺度、形状的广义可加模型(简称GAMLSS)建立控制因素降水、气温、ATD与年径流量之间的关系。在GAMLSS框架下,气温影响因子可以用两种形式表示,一种是直接采用气温,另一种是采取ATD指数(累积气温亏损值)。通过比较不同解释变量组合下的GAMLSS模型,进而研究气温对长江上游巴塘站1960~2012年的年径流影响。结果表明:基于ATD的回归模型,在年径流序列服从对数正态分布假设的条件下拟合效果最优。与气温值相比,ATD指数能更有效地解释长江上游径流变化的特征和冰川产流的物理机制。研究成果对长江上游年径流预报、高原气候下的产流特征分析具有理论意义。
To analyze the impacts of air temperature on annual runoff in the headstream of Yangtze River and explain the physical mechanism behind the glacier melt from Qinghai-Tibetan Plateau,the Generalized Additive Model of Location,Scale and Shape(GAMLSS)is employed to detect the non-stationarity of the annual runoff series(1960-2012)and to quantify the relationships of both the mean and variance of the annual runoff variable to physical factors such as precipitation and air temperature.In the GAMLSS model,the impact of air temperature can be represented in two different ways,one is to directly use the annual average temperature as a predictor and the other way is to use the index of accumulated temperature deficit(ATD)as a predictor,so two kinds of regression models are established in this paper under four different distribution assumptions for annual runoff variable.By comparing the efficiencies of the two kinds of GAMLSS models with different combinations of predictors using a variety of criteria such as the generalized Akaike information criterion,both normal QQ and worm plots of the residuals,Filliben correlation coefficient and kernel density estimation,we investigated the impacts of air temperature on annual runoff in the upper Yangtze River.We found that annual runoff series from the Batang station generally had an unstable trend,with a slight decline during the period 1960-1988 followed by a rise for over twenty years later,while annual precipitation and air temperature values showed continuous,steady and slow increasing trends.The fact that there existed a decline in annual runoff for nearly 30 years can be partly explained by the phenomenon called refreezing of meltwater of glaciers.The fitted model with ATD as a predictor under the lognormal distribution assumption for annual runoff performed better than the direct use of air temperature as the predictor in estimating annual runoff values.The ATD series,as an indirect reflection of heat energy deficit,showed roughly the same change trend as annua