Estimation of future precipitation change in the Yangtze River basin by using statistical downscaling method

In this study, the applicability of the statistical downscaling model (SDSM) in downscaling precipitation in the Yangtze River basin, China was investigated. The investigation includes the calibration of the SDSM model by using large-scale atmospheric variables encompassing NCEP/NCAR reanalysis data, the validation of the model using independent period of the NCEP/NCAR reanalysis data and the general circulation model (GCM) outputs of scenarios A2 and B2 of the HadCM3 model, and the prediction of the future regional precipitation scenarios. Selected as climate variables for downscaling were measured daily precipitation data (1961–2000) from 136 weather stations in the Yangtze River basin. The results showed that: (1) there existed good relationship between the observed and simulated precipitation during the calibration period of 1961–1990 as well as the validation period of 1991–2000. And the results of simulated monthly and seasonal precipitation were better than that of daily. The average R ² values between the simulated and observed monthly and seasonal precipitation for the validation period were 0.78 and 0.91 respectively for the whole basin, which showed that the SDSM had a good applicability on simulating precipitation in the Yangtze River basin. (2) Under both scenarios A2 and B2, during the prediction period of 2010–2099, the change of annual mean precipitation in the Yangtze River basin would present a trend of deficit precipitation in 2020s; insignificant changes in the 2050s; and a surplus of precipitation in the 2080s as compared to the mean values of the base period. The annual mean precipitation would increase by about 15.29% under scenario A2 and increase by about 5.33% under scenario B2 in the 2080s. The winter and autumn might be the more distinct seasons with more predicted changes of precipitation than in other seasons. And (3) there would be distinctive spatial distribution differences for the change of annual mean precipitation in the river basin, but the most of Yangtze River basin would be dominated by the increasing trend.     

Comparison of direct statistical and indirect statistical-deterministic frameworks in downscaling river low-flow indices

ABSTRACT

This work explores the ability of two methodologies in downscaling hydrological indices characterizing the low flow regime of three salmon rivers in Eastern Canada: Moisie, Romaine and Ouelle. The selected indices describe four aspects of the low flow regime of these rivers: amplitude, frequency, variability and timing. The first methodology (direct downscaling) ascertains a direct link between large-scale atmospheric variables (the predictors) and low flow indices (the predictands). The second (indirect downscaling) involves downscaling precipitation and air temperature (local climate variables) that are introduced into a hydrological model to simulate flows. Synthetic flow time series are subsequently used to calculate the low flow indices. The statistical models used for downscaling low flow hydrological indices and local climate variables are: Sparse Bayesian Learning and Multiple Linear Regression. The results showed that direct downscaling using Sparse Bayesian Learning surpassed the other approaches with respect to goodness of fit and generalization ability.

Editor D. Koutsoyiannis; Associate editor K. Hamed     

Assessment of change in design flood frequency under climate change using a multivariate downscaling model and a precipitation-runoff model

Precipitation and runoff are key elements in the hydrologic cycle because of their important roles in water supply, flood prevention, river restoration, and ecosystem management. Global climate change, widely accepted to be happening, is anticipated to have enormous consequences on future hydrologic patterns. Studies on the potential changes in global, regional, and local hydrologic patterns under global climate change scenarios have been an intense area of research in recent years. The present study contributes to this research topic through evaluation of design flood under climate change. The study utilizes a weather state-based, stochastic multivariate model as a conditional probability model for simulating the precipitation field. An important premise of this study is that large-scale climatic patterns serve as a major driver of persistent year-to-year changes in precipitation probabilities. Since uncertainty estimation in the study of climate change is needed to examine the reliability of the outcomes, this study also applies a Bayesian Markov chain Monte Carlo scheme to the widely used SAC-SMA (Sacramento soil moisture accounting) precipitation-runoff model. A case study is also performed with the Soyang Dam watershed in South Korea as the study basin. Finally, a comprehensive discussion on design flood under climate change is made.     

A comparison of three methods for downscaling daily precipitation in the Punjab region

Many downscaling techniques have been developed in the past few years for projection of station‐scale hydrological variables from large‐scale atmospheric variables simulated by general circulation models (GCMs) to assess the hydrological impacts of climate change. This article compares the performances of three downscaling methods, viz. conditional random field (CRF), K‐nearest neighbour (KNN) and support vector machine (SVM) methods in downscaling precipitation in the Punjab region of India, belonging to the monsoon regime. The CRF model is a recently developed method for downscaling hydrological variables in a probabilistic framework, while the SVM model is a popular machine learning tool useful in terms of its ability to generalize and capture nonlinear relationships between predictors and predictand. The KNN model is an analogue‐type method that queries days similar to a given feature vector from the training data and classifies future days by random sampling from a weighted set of K closest training examples. The models are applied for downscaling monsoon (June to September) daily precipitation at six locations in Punjab. Model performances with respect to reproduction of various statistics such as dry and wet spell length distributions, daily rainfall distribution, and intersite correlations are examined. It is found that the CRF and KNN models perform slightly better than the SVM model in reproducing most daily rainfall statistics. These models are then used to project future precipitation at the six locations. Output from the Canadian global climate model (CGCM3) GCM for three scenarios, viz. A1B, A2, and B1 is used for projection of future precipitation. The projections show a change in probability density functions of daily rainfall amount and changes in the wet and dry spell distributions of daily precipitation. Copyright © 2011 John Wiley & Sons, Ltd.     

Projection of future rainfall for the North China Plain using two statistical downscaling models and its hydrological implications

This study projected the future rainfall (2046–2065 and 2081–2100) for the North China Plain (NCP) using two stochastic statistical downscaling models, the non-homogeneous hidden Markov model and the generalized linear model for daily climate time series, conditioned by the large-scale atmospheric predictors from six general circulation models for three emission scenarios (A1B, A2 and B1). The results indicated that the annual total rainfall, the extreme daily rainfall and the maximum length of consecutive wet/dry days would decline, while the number of annual rainfall days would slightly increase (correspondingly rainfall intensity would decrease) in the NCP, in comparison with the base period (1961–2010). Moreover, the summer monsoon rainfall, which accounted for 50–75 % of the total annual rainfalls in NCP, was projected to decrease in the latter half of twenty-first century. The spatial patterns of change showed generally north–south gradients with relatively larger magnitude decrease in the northern NCP and less decrease (or even slightly increase) in the southern NCP. This could result in decline of the annual runoff by ?5.5 % (A1B), ?3.3 % (A2) and ?4.1 % (B1) for 2046–2065 and ?5.3 % (A1B), ?4.6 % (A2) and ?1.9 % (B1) decrease for 2081–2100. These rainfall changes, combined with the warming temperature, could lead to drier catchment soil profiles and further reduce runoff potential, would hence provide valuable references for the water availability and related climate change adaption in the NCP.     

Comparison of transfer functions in statistical downscaling models for daily temperature and precipitation over Canada

This study compares three linear models and one non-linear model, specifically multiple linear regression (MLR) with ordinary least squares (OLS) estimates, robust regression, ridge regression, and artificial neural networks (ANNs), to identify an appropriate transfer function in statistical downscaling (SD) models for the daily maximum and minimum temperatures (T_max and T_min) and daily precipitation occurrence and amounts (P_occ and P_amount). This comparison was made over twenty-five observation sites located in five different Canadian provinces (British Columbia, Saskatchewan, Manitoba, Ontario, and Québec). Reanalysis data were employed as atmospheric predictor variables of SD models. Predictors of linear transfer functions and ANN were selected by linear correlations coefficient and mutual information, respectively. For each downscaled case, annual and monthly models were developed and analysed. The monthly MLR, annual ANN, annual ANN, and annual MLR yielded the best performance for T_max, T_min, P_occ and P_amont according to the modified Akaike information criterion (AICu). A monthly MLR is recommended for the transfer functions of the four predictands because it can provide a better performance for the T_max and as good performance as the annual MLR for the T_min, P_occ, and P_amount. Furthermore, a monthly MLR can provide a slightly better performance than an annual MLR for extreme events. An annual MLR approach is also equivalently recommended for the transfer functions of the four predictands because it showed as good a performance as monthly MLR in spite of its mathematical simplicity. Robust and ridge regressions are not recommended because the data used in this study are not greatly affected by outlier data and multicollinearity problems. An annual ANN is recommended only for the T_min, based on the best performance among the models in terms of both the RMSE and AICu.     

Prediction of daily rainfall state in a river basin using statistical downscaling from GCM output

Conventional statistical downscaling techniques for prediction of multi-site rainfall in a river basin fail to capture the correlation between multiple sites and thus are inadequate to model the variability of rainfall. The present study addresses this problem through representation of the pattern of multi-site rainfall using rainfall state in a river basin. A model based on K-means clustering technique coupled with a supervised data classification technique, namely Classification And Regression Tree (CART), is used for generation of rainfall states from large-scale atmospheric variables in a river basin. The K-means clustering is used to derive the daily rainfall state from the historical daily multi-site rainfall data. The optimum number of clusters in the observed rainfall data is obtained after application of various cluster validity measures to the clustered data. The CART model is then trained to establish relationship between the daily rainfall state of the river basin and the standardized, dimensionally-reduced National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis climatic data set. The relationship thus developed is applied to the General Circulation Model (GCM)-simulated, standardized, bias free large-scale climate variables for prediction of rainfall states in future. Comparisons of the number of days falling under different rainfall states for the observed period and the future give the change expected in the river basin due to global warming. The methodology is tested for the Mahanadi river basin in India.     

Recent trends of extreme precipitation indices in the Iberian Peninsula using observations and WRF model results

Spatial and temporal distributions of the trends of extreme precipitation indices were analysed between 1986 and 2005, over the Iberian Peninsula (IP). The knowledge of the patterns of extreme precipitation is important for impacts assessment, development of adaptation and mitigation strategies. As such, there is a growing need for a more detailed knowledge of precipitation climate change.This analysis was performed for Portuguese and Spanish observational datasets and results performed by the Weather Research and Forecast (WRF) model forced by the ERA-Interim reanalysis. Extreme precipitation indices recommended by the Expert Team for Climate Change Detection Monitoring and Indices were computed, by year and season. Then, annual and seasonal trends of the indices were estimated by Theil-Sen method and their significance was tested by the Mann-Kendal test. Additionally, a second simulation forced by the Max Planck Institute Earth System Model (MPI-ESM), was considered. This second modelling configuration was created in order to assess its performance when simulating extremes of precipitation.The annual trends estimated for the 1986–2005, from the observational datasets and from the ERA-driven simulation reveal: 1) negative statistically significant trends of the CWD index in the Galicia and in the centre of the IP; 2) positive statistically significant trends of the CDD index over the south of the IP and negative statistically significant trends in Galicia, north and centre of Portugal; 3) positive statistically significant trends of the R75p index in some regions of the north of the IP; 4) positive statistically significant trends in the R95pTOT index in the Central Mountains Chain, Leon Mountains and in the north of Portugal.Seasonally, negative statistically significant trends of the CWD index were found in Galicia, in winter and in the south of the IP, in summer. Positive statistically significant trends of the CWD index were identified in the Leon Mountains, in spring, and in Galicia, in autumn. For the CDD index, negative statistically significant trends were seen in Valencia, in the spring, and, in Galicia and Portugal (north and centre), in summer. Positive statistically significant trends of the CDD index were found: in the east of the IP, in the winter; in the Cantabrian Mountain, in the spring; and, in the south of the IP, in summer. Regarding to the R75p index, negative statistically significant trends were found in Galicia, in winter and positive statistically significant trends in the north of Portugal, in spring and in the Central Mountains Chain and north of Portugal, in autumn. For the R95pTOT index, negative statistically significant trends were found over the Sierra Cuenca and Sierra Cazorla, in winter and positive statistically significant trends were found over the Sierra Cebollera, in winter and in Castile-la Mancha region, in spring.The results of the annual and seasonal trends of the extreme precipitation indices performed for observational datasets and the simulation forced by ERA-Interim, are similar. The results obtained for the simulation forced by MPI-ESM are not satisfactory, and can be a source of criticism for the use of simulation forced by MPI-ESM in this type of climate change studies. Even for the relatively short period used, the WRF model, when properly forced is a useful tool due to the similar results of Portuguese and Spanish observational datasets and the simulation forced by ERA-Interim.     

Prediction of variability of precipitation in the Yangtze River Basin under the climate change conditions based on automated statistical downscaling

Many impact studies require climate change information at a finer resolution than that provided by general circulation models (GCMs). Therefore the outputs from GCMs have to be downscaled to obtain the finer resolution climate change scenarios. In this study, an automated statistical downscaling (ASD) regression-based approach is proposed for predicting the daily precipitation of 138 main meteorological stations in the Yangtze River basin for 2010–2099 by statistical downscaling of the outputs of general circulation model (HadCM3) under A2 and B2 scenarios. After that, the spatial–temporal changes of the amount and the extremes of predicted precipitation in the Yangtze River basin are investigated by Mann–Kendall trend test and spatial interpolation. The results showed that: (1) the amount and the change pattern of precipitation could be reasonably simulated by ASD; (2) the predicted annual precipitation will decrease in all sub-catchments during 2020s, while increase in all sub-catchments of the Yangtze River Basin during 2050s and during 2080s, respectively, under A2 scenario. However, they have mix-trend in each sub-catchment of Yangtze River basin during 2020s, but increase in all sub-catchments during 2050s and 2080s, except for Hanjiang River region during 2080s, as far as B2 scenario is concerned; and (3) the significant increasing trend of the precipitation intensity and maximum precipitation are mainly occurred in the northwest upper part and the middle part of the Yangtze River basin for the whole year and summer under both climate change scenarios and the middle of 2040–2060 can be regarded as the starting point for pattern change of precipitation maxima.     

Uncertainty analysis of statistical downscaling methods using Canadian Global Climate Model predictors

Three downscaling models, namely the Statistical Down‐Scaling Model (SDSM), the Long Ashton Research Station Weather Generator (LARS‐WG) model and an artificial neural network (ANN) model, have been compared in terms of various uncertainty attributes exhibited in their downscaled results of daily precipitation, daily maximum and minimum temperature. The uncertainty attributes are described by the model errors and the 95% confidence intervals in the estimates of means and variances of downscaled data. The significance of those errors has been examined by suitable statistical tests at the 95% confidence level. The 95% confidence intervals in the estimates of means and variances of downscaled data have been estimated using the bootstrapping method and compared with the observed data. The study has been carried out using 40 years of observed and downscaled daily precipitation data and daily maximum and minimum temperature data, starting from 1961 to 2000. In all the downscaling experiments, the simulated predictors of the Canadian Global Climate Model (CGCM1) have been used. The uncertainty assessment results indicate that, in daily precipitation downscaling, the LARS‐WG model errors are significant at the 95% confidence level only in a very few months, the SDSM errors are significant in some months, and the ANN model errors are significant in almost all months of the year. In downscaling daily maximum and minimum temperature, the performance of all three models is similar in terms of model errors evaluation at the 95% confidence level. But, according to the evaluation of variability and uncertainty in the estimates of means and variances of downscaled precipitation and temperature, the performances of the LARS‐WG model and the SDSM are almost similar, whereas the ANN model performance is found to be poor in that consideration. Further assessment of those models, in terms of skewness and average dry‐spell length comparison between observed and downscaled daily precipitation, indicates that the downscaled daily precipitation skewness and average dry‐spell lengths of the LARS‐WG model and the SDSM are closer to the observed data, whereas the ANN model downscaled precipitation underestimated those statistics in all months. Copyright © 2006 John Wiley & Sons, Ltd.     

Statistical downscaling of extreme daily precipitation,evaporation, and temperature and construction of future scenarios

Generally, the statistical downscaling approaches work less perfectly in reproducing precipitation than temperatures, particularly for the extreme precipitation. This article aimed to testify the capability in downscaling the extreme temperature, evaporation, and precipitation in South China using the statistical downscaling method. Meanwhile, the linkages between the underlying driving forces and the incompetent skills in downscaling precipitation extremes over South China need to be extensively addressed. Toward this end, a statistical downscaling model (SDSM) was built up to construct future scenarios of extreme daily temperature, pan evaporation, and precipitation. The model was thereafter applied to project climate extremes in the Dongjiang River basin in the 21st century from the HadCM3 (Hadley Centre Coupled Model version 3) model under A2 and B2 emission scenarios. The results showed that: (1) The SDSM generally performed fairly well in reproducing the extreme temperature. For the extreme precipitation, the performance of the model was less satisfactory than temperature and evaporation. (2) Both A2 and B2 scenarios projected increases in temperature extremes in all seasons; however, the projections of change in precipitation and evaporation extremes were not consistent with temperature extremes. (3) Skills of SDSM to reproduce the extreme precipitation were very limited. This was partly due to the high randomicity and nonlinearity dominated in extreme precipitation process over the Dongjiang River basin. In pre‐flood seasons (April to June), the mixing of the dry and cold air originated from northern China and the moist warm air releases excessive rainstorms to this basin, while in post‐flood seasons (July to October), the intensive rainstorms are triggered by the tropical system dominated in South China. These unique characteristics collectively account for the incompetent skills of SDSM in reproducing precipitation extremes in South China. Copyright © 2011 John Wiley & Sons, Ltd.     

Intercomparison of multiple statistical downscaling methods: multi-criteria model selection for South Korea

A number of statistical downscaling methodologies have been introduced to bridge the gap in scale between outputs of climate models and climate information needed to assess potential impacts at local and regional scales. Four statistical downscaling methods [bias-correction/spatial disaggregation (BCSD), bias-correction/constructed analogue (BCCA), multivariate adaptive constructed analogs (MACA), and bias-correction/climate imprint (BCCI)] are applied to downscale the latest climate forecast system reanalysis (CFSR) data to stations for precipitation, maximum temperature, and minimum temperature over South Korea. All methods are calibrated with observational station data for 19 years from 1973 to 1991 and validated for the more recent 19-year period from 1992 to 2010. We construct a comprehensive suite of performance metrics to inter-compare methods, which is comprised of five criteria related to time-series, distribution, multi-day persistence, extremes, and spatial structure. Based on the performance metrics, we employ technique for order of preference by similarity to ideal solution (TOPSIS) and apply 10,000 different weighting combinations to the criteria of performance metrics to identify a robust statistical downscaling method and important criteria. The results show that MACA and BCSD have comparable skill in the time-series related criterion and BCSD outperforms other methods in distribution and extremes related criteria. In addition, MACA and BCCA, which incorporate spatial patterns, show higher skill in the multi-day persistence criterion for temperature, while BCSD shows the highest skill for precipitation. For the spatial structure related criterion, BCCA and MACA outperformed BCSD and BCCI. From the TOPSIS analysis, we found that MACA is the most robust method for all variables in South Korea, and BCCA and BCSD are the second for temperature and precipitation, respectively. We also found that the contribution of the multi-day persistence and spatial structure related criteria are crucial to ranking the skill of statistical downscaling methods.     

Statistical downscaling of extremes of precipitation and temperature and construction of their future scenarios in an elevated and cold zone

Reliable projections of extremes at finer spatial scales are important in assessing the potential impacts of climate change on societal and natural systems, particularly for elevated and cold regions in the Tibetan Plateau. This paper presents future projections of extremes of daily precipitation and temperature, under different future scenarios in the headwater catchment of Yellow River basin over the 21st century, using the statistical downscaling model (SDSM). The results indicate that: (1) although the mean temperature was simulated perfectly, followed by monthly pan evaporation, the skill scores in simulating extreme indices of precipitation are inadequate; (2) The inter-annual variabilities for most extreme indices were underestimated, although the model could reproduce the extreme temperatures well. In fact, the simulation of extreme indices for precipitation and evaporation were not satisfactory in many cases. (3) In future period from 2011 to 2100, increases in the temperature and evaporation indices are projected under a range of climate scenarios, although decreasing mean and maximum precipitation are found in summer during 2020s. The findings of this work will contribute toward a better understanding of future climate changes for this unique region.     

Retrospective forecasting test of a statistical physics model for earthquakes in Sichuan-Yunnan region

Pattern informatics (PI) model is one of the recently developed predictive models of earthquake phys- ics based on the statistical mechanics of complex systems. In this paper, retrospective forecast test of the PI model was conducted for the earthquakes in Sichuan-Yunnan region since 1988, exploring the possibility to apply this model to the estimation of time-dependent seismic hazard in continental China. Regional earthquake catalogue down to ML3.0 from 1970 to 2007 was used. The ‘target magnitude’ for the forecast test was MS5.5. Fifteen-year long ‘sliding time window’ was used in the PI calculation, with ‘anomaly training time window’ being 5 years and ‘forecast time window’ being 5 years, respectively. Receiver operating characteristic (ROC) test was conducted for the evaluation of the forecast result, showing that the PI forecast outperforms not only random guess but also the simple number counting approach based on the clustering hypothesis of earthquakes (the RI forecast). If the ‘forecast time window’ was shortened to 3 years and 1 year, respectively, the forecast capability of the PI model de- creased significantly, albeit outperformed random forecast. For the one year ‘forecast time window’, the PI result was almost comparable to the RI result, indicating that clustering properties play a more important role at this time scale.     

Urbanization and changes in precipitation, a statistical approach

Monthly rainfall data from urbanized and rural areas in the western part of The Netherlands were used to investigate the urban effects upon the rainfall regime. A trend test in the sequence of differences between urban and rural rainfall amounts is described. For the urban areas of Amsterdam and Rotterdam some evidence was found for an increase in precipitation.     

Spatial downscaling of TRMM precipitation data based on the orographical effect and meteorological conditions in a mountainous area

The lack of high resolution precipitation data has posed great challenges to the study and management of extreme rainfall events. Satellite-based rainfall products with large areal coverage provide a potential alternative source of data where in situ measurements are not available. However, the mismatch in scale between these products and model requirements has limited their application and demonstrates that satellite data must be downscaled before being used. This study developed a statistical spatial downscaling scheme based on the relationships between precipitation and related environmental factors such as local topography and pre-storm meteorological conditions. The method was applied to disaggregate the Tropical Rainfall Measuring Mission (TRMM) 3B42 products, which have a resolution of 0.25° × 0.25°, to 1 × 1 km gridded rainfall fields. The TRMM datasets in accord with six rainstorm events in the Xiao River basin were used to validate the effectiveness of this approach. The downscaled precipitation data were compared with ground observations and exhibited good agreement with r² values ranging from 0.612 to 0.838. In addition, the proposed approach provided better results than the conventional spline and kriging interpolation methods, indicating its promise in the management of extreme rainfall events. The uncertainties in the final results and the implications for further study were discussed, and the needs for additional rigorous investigations of the rainfall physical process prior to institutionalizing the use of satellite data were highlighted.     

Comparison of SDSM and LARS-WG for simulation and downscaling of extreme precipitation events in a watershed

Future climate projections of Global Climate Models (GCMs) under different emission scenarios are usually used for developing climate change mitigation and adaptation strategies. However, the existing GCMs have only limited ability to simulate the complex and local climate features, such as precipitation. Furthermore, the outputs provided by GCMs are too coarse to be useful in hydrologic impact assessment models, as these models require information at much finer scales. Therefore, downscaling of GCM outputs is usually employed to provide fine-resolution information required for impact models. Among the downscaling techniques based on statistical principles, multiple regression and weather generator are considered to be more popular, as they are computationally less demanding than the other downscaling techniques. In the present study, the performances of a multiple regression model (called SDSM) and a weather generator (called LARS-WG) are evaluated in terms of their ability to simulate the frequency of extreme precipitation events of current climate and downscaling of future extreme events. Areal average daily precipitation data of the Clutha watershed located in South Island, New Zealand, are used as baseline data in the analysis. Precipitation frequency analysis is performed by fitting the Generalized Extreme Value (GEV) distribution to the observed, the SDSM simulated/downscaled, and the LARS-WG simulated/downscaled annual maximum (AM) series. The computations are performed for five return periods: 10-, 20-, 40-, 50- and 100-year. The present results illustrate that both models have similar and good ability to simulate the extreme precipitation events and, thus, can be adopted with confidence for climate change impact studies of this nature.     

基于随机森林模型的GRACE数据3种空间降尺度对比

陆地水储量是赋存在陆地上各种形式水的综合体现,研究其时空变化对认识区域水循环过程和水资源调控等具有重要意义。然而现有陆地水储量变化数据实际分辨率较低,限制了其在中小流域或地区中的应用。针对这一问题,本文基于GRACE重力卫星和其后续卫星GRACE-FO反演的陆地水储量变化数据,首先采用随机森林模型,分别基于格点、区域(流域)和区域(全国)3种空间降尺度思路将GRACE数据降尺度至0.25°×0.25°,后结合GLDAS模型数据,基于水量平衡原理计算得到地下水储量变化数据,最后基于降尺度模型模拟效果和实测地下水位数据评估3种降尺度思路在全国的适用性。结果表明:随机森林模型能够较好地模拟驱动数据(降水、气温、植被条件指数和土壤水储量)与GRACE数据的统计关系,验证期格点降尺度思路的平均相关系数总体在0.6左右,区域降尺度思路的平均纳什效率系数、相关系数和均方根误差分别>0.5、>0.75和<6.6 cm,3种空间降尺度思路的模拟精度均满足基本要求;2003—2021年间,GRACE数据、格点降尺度、区域降尺度(流域)和区域降尺度(全国)得到的我国陆地水储量亏缺量分别约为...