Climate and Hydrological Change Characteristics and ... - MDPI

7 downloads 0 Views 3MB Size Report
Mar 1, 2018 - State Key Laboratory of Hydraulics and Mountain River Engineering, .... regional research, including the Yangtze River, the Yellow River, the ...
water Article

Climate and Hydrological Change Characteristics and Applicability of GLDAS Data in the Yarlung Zangbo River Basin, China Hong Zhang, Ling Lei Zhang *, Jia Li *, Rui Dong An and Yun Deng State Key Laboratory of Hydraulics and Mountain River Engineering, College of Water Resource & Hydropower, Sichuan University, Chengdu 610065, China; [email protected] (H.Z.); [email protected] (R.D.A.); [email protected] (Y.D.) * Correspondence: [email protected] (L.L.Z.); [email protected] (J.L.); Tel.: +86-028-8540-7075 (L.L.Z.); +86-028-8540-1158 (J.L.) Received: 20 December 2017; Accepted: 26 February 2018; Published: 1 March 2018

Abstract: The hydrological cycle is particularly sensitive to and is greatly affected by global climate change. In addition, runoff change has a strong influence on the hydrological cycle and migration of biogenic substances. The Yarlung Zangbo River basin in China is a typical basin for which climate and hydrological data are lacking. Land surface models can provide data for studying land surface substance and energy circulation, which are meaningful to face climate change. The midstream region of the Yarlung Zangbo River basin, which is strongly affected by climate change, was selected as the study area. First, the observed mean temperature, precipitation and runoff characteristics were analysed. Second, after combining the Global Land Data Assimilation System (GLDAS) and the water balance equation, we simulated climate and hydrological processes for the same time period. Finally, the correlation and error between GLDAS and observed data were analysed to verify applicability of the GLDAS data, and the impacts of climate factors on runoff were discussed. The results revealed that under the background of global warming, precipitation, temperature, and runoff changed significantly and showed strong consistency during the research period. Mean monthly precipitation, temperature and runoff exhibited clear cyclical fluctuations of approximately 12 months, and they all tended to increase. GLDAS is not a good system to describe the land surface conditions of the Yarlung Zangbo River basin all the time. However, within a certain time period, GLDAS data have a good applicability in the basin. Thereinto, the GLDAS mean monthly precipitation was moderately correlated with observed precipitation, with a correlation coefficient of 0.75. GLDAS mean monthly temperature was highly correlated with observed data, with a correlation coefficient of 0.94. Based on the Brunke ranking method, it indicates that GLDAS-Noah-based runoff data were closer to observed runoff data than the three other GLDAS models. Correlation coefficients between precipitation and runoff for the three time scales were higher than those between temperature and runoff. This means that rainfall was the main factor affecting natural runoff change, as opposed to temperature, and it can control the evolution of the river to some degree. This paper indicates the impacts of climate change on runoff and the application of GLDAS for data-limited basins. The results provide a deeper understanding of the Yarlung Zangbo River basin characteristics and can provide a scientific basis for the management of water resources and policy implementation for this basin. Keywords: the Yarlung Zangbo River; climate; runoff; GLDAS; applicability

1. Introduction Global warming is becoming an indisputable truth. The atmospheric carbon dioxide concentration has increased from 280 ppm to 407 ppm since the industrial revolution and will continue to increase

Water 2018, 10, 254; doi:10.3390/w10030254

www.mdpi.com/journal/water

Water 2018, 10, 254

2 of 16

at a rate of 1.8 ppm/a [1]. After five intensive discussions, the Intergovernmental Panel on Climate Change (IPCC) reported the following climatic and hydrological problems [2–5]: Mean temperature of the sea and land surface has exhibited a linear upward trend over the past 100 years, which has led to sea level rise and the snowfall and glaciers decreased. Extreme weather (i.e., heavy rainfall, heat waves, flooding and drought) due to global warming threatens human life, and the occurrence probability may increase. In China, the rates of warming and sea level rise are faster than the mean global rate. The glacier area has decreased by 10.1%, and frozen areas on earth have decreased by approximately 18.6% since the late 1970s. China’s air temperature will continue to increase in the future. Under the background of global climate change, the global water cycle will change inevitably. The increase of greenhouse gas concentrations will have direct first-order effects on the global hydrological cycle, with impacts on water availability and demand [6]. The critical importance of water is undeniable-it is the essence of life. Insufficient rainfall results in a decrease in water supply and quality, crop failure, and famine. Too much water results in flooding, destruction, poor water quality, and famine [7]. Looking forward, water resources will become key and will play a restrictive role in the global environment and economic development on a daily basis. Research on the effects of climate change on water resources started in the late 1970s, and included the following: relationships between climate, climate change and water supply [8], the extent of climate prediction and the impact of human activity on climate [9], impact of climate change and fluctuation on water resources and the sensitivity of water resources to climate change [10,11], and large-scale climate and hydrological simulation technology [12]. The main results of these research are as follows: (1) Higher rainfall causes runoff to increase in high latitude regions, and runoff decreases in low latitude regions due to decreased rainfall and increased evaporation; (2) Mean annual runoff is more sensitive to precipitation than to temperature; (3) In the future, the frequency and intensity of drought may increase in areas with higher temperature and lower precipitation; (4) In basins with seasonal snowfall and melting snow, seasonal runoff is more sensitive to temperature changes than to other climate changes. Since the 1980s, China has also carried out a series of research to determine the impacts of gradual climate change on hydrology and water resources. However, Chinese scholars have focused mainly on regional research, including the Yangtze River, the Yellow River, the Hu River, Ning Xia Province, and other regions [13–16]. This research showed that under the background of global warming, mean temperature, and precipitation in major regions of China have changed significantly in recent years, with significant impacts on water resources. The results can provide a scientific basis for the management of climate and water resources and for the formulation of relevant policies. The Tibetan Plateau is the highest plateau in the world. Due to its unique geographic environment, the energy and water cycles of this area contribute substantially to global energy and water cycles. In addition, the exchange of energy and water between the land surface and atmosphere play an important role in regional climate change and water resource management [17]. Therefore, the Tibetan Plateau is also known as “a driver and amplifier of global climate change”. The Yarlung Zangbo River basin, which is located in the southern Tibetan Plateau, is the main source of fresh water in the Tibet autonomous region, and is also a strategic water resource reserve area in China. Under the background of global warming, the climate in this area has changed significantly. The temperature is increasing significantly, especially in winter. Precipitation fluctuates considerably, and the regional difference in climate change is significant. As global warming intensifies, the water cycle is strengthened, especially in high latitude and alpine regions, and water resources are strongly affected by climate change. Hence, it is necessary to carry out researches on water cycle in the basin under climate change. The weather stations in western China are sparsely distributed, especially on the Tibetan Plateau, and the study of climate change in this region depends only on observed data, which are insufficient. Therefore, we need more accurate and high-resolution grid data to cover the whole area, and the Global Land Data Assimilation System (GLDAS) can provide data support for this requirement.

Water 2018, 10, 254

3 of 16

Midstream of the Yarlung Zangbo River basin is a core region of economic and social development in Tibet. The population in this area is very high, and water resource use is relatively high. However, this area is strongly affected by climate change. Hence, carrying out a study on runoff response under climate change in this area has important practical significance. In this paper, we selected this area as the study area. First, characteristics of observed temperature, precipitation and runoff change of Nuxia hydrologic station from 1979 to 2016 were analysed. Second, four GLDAS models were used to extract temperature, precipitation, soil moisture, and snow equivalent data for the same period. The water balance equation and GLDAS data were used to simulate runoff change. Finally, the GLDAS data in the study area were verified by comparison with observed data to analyse the application, and the impact of climate factors on runoff was discussed. This research is helpful to better understand the impact of climate change on water resources in the midstream region of the Yarlung Zangbo River basin and provides significant guidance for water resource management, related policy formulation, drought disaster, agricultural irrigation, and environmental problems that are caused by the reduction of hydrological flow. This research is also useful to study the response mechanisms of hydrology and ecological processes to climate change at the basin scale. 2. Objective, Study Area, Data and Methods 2.1. Objective Under the background of global climate change, this research targets to the core rule of the applicability of the assimilation system data in data-limited area. With the midstream region of the Yarlung Zangbo River basin taken as study object, this research systematically focus on the analyse of climate and hydrology change characteristics, and the applicability of four models of GLDAS data in this area. 2.2. Study Area The Yarlung Zangbo River, which originates in Jie Ma Zong Yang glacier, runs through southern Tibet and leaves China through Pasighat. It is an important international river that has the highest mean altitude. The main stream in China is approximately 2.06 × 103 km, ranging from 82◦ 010 E to 97◦ 060 E and from 27◦ 400 N to 31◦ 170 N. The catchment area is approximately 2.58 × 105 km2 , and mean annual runoff is approximately 1.40 × 1011 m3 , ranking third in China [18]. This basin is a typical subtropical climate zone that is located at the intersection of the southern and central Asian tropics and is greatly affected by the south Asian summer monsoon. The study area of this paper is the midstream region of this river, which runs from Lizi to Pai town. Its total length is 1300 km, and the control area is 1.64 × 105 km2 . The major tributaries are located in this area, and these huge tributaries not only provide abundant water resources, but also maintain a broad plain, which is the most important and richest agricultural region in Tibet. Therefore, understanding the spatial and temporal dynamics and determining the trend of the water cycle under climate change are meaningful for agricultural production, as well as a precondition to reasonably allot and develop water resources, leading to the sustainable use of regional water resources [19]. The geographic location, drainage distribution, and related hydrological stations of the study area can be observed in Figure 1.

important and richest agricultural region in Tibet. Therefore, understanding the spatial and temporal dynamics and determining the trend of the water cycle under climate change are meaningful for agricultural production, as well as a precondition to reasonably allot and develop water resources, leading to the sustainable use of regional water resources [19]. The geographic location, drainage Water 2018, 10, 254 4 of 16 distribution, and related hydrological stations of the study area can be observed in Figure 1.

Water 2018, 10, x FOR PEER REVIEW

4 of 15

(a)

(b) FigureFigure 1. (a) Geographic location of of the study (b) The Thewater water system map ofYarlung the Yarlung Zangbo 1. (a) Geographic location the studyarea; area; (b) system map of the Zangbo River basin. The green is the midstreamregion region of and thethe blueblue line line is theisupstream and and River basin. The green lineline is the midstream of this thisriver, river, and the upstream downstream region of this river. downstream region of this river. . 2.3. Data 2.3. Data 2.3.1. Observed Meteorological and Hydrological Data

2.3.1. Observed Meteorological and Hydrological Data

Observed meteorological (i.e., precipitation and temperature) data are published by the China meteorological data service(i.e., centre (http://data.cma.cn). The monthly mean Observed meteorological precipitation and temperature) data areprecipitation published (P) by and the China temperature (T) from January 1979 to December 2016 of 20 meteorological stations in and near meteorological data service centre (http://data.cma.cn). The monthly mean precipitationthe (P) and midstream of this basin were downloaded2016 for this The 20 meteorological temperature (T) region from January 1979 to December of paper. 20 meteorological stationsstations in andare near the distributed unevenly because of manpower and resource limitations, as well as the natural environment. midstream region of this basin were downloaded for this paper. The 20 meteorological stations are Hence, to obtain mean monthly values, related measures should be taken. distributedUsing unevenly because of manpower and resource limitations, as well as the natural ArcGIS software, the Thiessen polygon method is used to calculate mean precipitation of a environment. Hence, mean in monthly related basin [20], which to canobtain be expressed terms of values, the formula (1). measures should be taken. UsingThe ArcGIS software, Thiessen polygon method is used to calculate mean precipitation of Thiessen polygonthe formula:

a basin [20], which can be expressed in terms of the formula (1). p1 f 1 + p2 f 2 + . . . + p n f n The Thiessen polygon formula: P = F

(1)

p f + p f 2 +meteorological … + p n f n station; f 1 , f 2 , . . . , fn represent where p1 , p2 , . . . , pn represent the precipitation P = 1 1 of2each (1) the area of each controlled meteorological station; and F F is the area of the whole basin. To calculate mean evapotranspiration of a basin, this paper adopts a traditional but high-precision pn represent precipitation of each meteorological station; f1, f2, …, fn represent the where method, p1, p2, …, called Takahashithe formula [20]. area of each controlled meteorological station; and F is the area of the whole basin. To calculate mean evapotranspiration of a basin, this paper adopts a traditional but highprecision method, called Takahashi formula [20]. The Takahashi formula: ET

3100 · P

= 3100

+ 1 . 8 · P 2 ·exp(

− 34 . 4 ·T ) 235 + T

(2)

Water 2018, 10, 254

5 of 16

The Takahashi formula: ET =

3100· P −34.4· T ) 3100 + 1.8· P2 · exp( 235 +T

(2)

where, P is monthly mean precipitation; T is monthly mean temperature. Observed runoff data for the Yarlung Zangbo River basin are extremely limited. There are four hydrological stations in the main stream of the midstream region, i.e., Lazi, Nugesha, Yongcun, and Nuxia stations. Figure 1b shows the location of the 4 gauging stations. We could only obtain partial observed runoff data from the 4 stations. In this paper, we used observed monthly runoff data from Nuxia gauging station from January 1979 to December 2000 to validate the results that were calculated from the GLDAS data. 2.3.2. Global Land Data Assimilation System (GLDAS) GLDAS is an offline global land simulation system that was jointly developed by the Goddard Space Flight Center (GSFC) of National Aeronautics and Space Administration (NASA) and the National Centers for Environmental Prediction (NCEP) of National Oceanic and Atmospheric Administration (NOAA) (https://hydro1.gesdisc.eosdis.nasa.gov/data/). The system contains 3 land models (i.e., Mosaic-MOS, Noah and a community land surface model—CLM), and 1 hydrology model (i.e., the variable infiltration capacity model—VIC). The Noah model is a land surface process model which regards the soil surface and vegetation as a whole. The CLM contains three important sub-models: land surface model of the National Center for Atmospheric Research (NCAR), Biosphere—Atmosphere Transfer Scheme, and Land Surface model (LSM) of the institute of atmospheric physics of Chinese academy of sciences. The VIC model is a large scale semi-distributed hydrological model based on thermal balance and physical dynamic mechanism. In addition, the system has two versions (i.e., GLDAS-1 and GLDAS-2). GLDAS-1 provides all four models of global land surface data from 1979 to present, with 1◦ × 1◦ of spatial resolutions and 3 h and one month of temporal resolutions. GLDAS-1 also provides data of Noah model from 2000 to present with 0.25◦ × 0.25◦ of spatial resolutions. To GLDAS-2, it only has Noah model to provide data from 1979 to 2010 with 1◦ × 1◦ of spatial resolutions. All four GLDAS models coalesce data from both of the land surface and satellites, and provide 28 optimized real-time land surface state variables (e.g., soil moisture, surface temperature, rainfall rate, snow water equivalent, and pressure) [21–23]. These high quality global land surface data are widely used in weather and climate research, water cycle research and applications of water resource. In this paper, all four GLDAS models from GLDAS-1 were used to provide monthly rainfall rate, air temperature, soil moisture, and snow water equivalent parameters for the 1.0◦ spatial resolution [17]. Precipitation was calculated by rainfall rate. By combining GLDAS data with the water balance equation, the runoff of this region was calculated. The purpose of this research is to analyse the characteristics of climate and hydrology change and to select the most suitable alternative climate and hydrology information for the midstream region of the Yarlung Zangbo River basin. 2.4. Methods 2.4.1. Comparison and Applicability Analysis of Meteorological Data When GLDAS data were compared with observed meteorological data, we first extracted monthly data according to the longitude and latitude of the research region. We selected GLDAS-MOS to assess the precipitation and temperature data because these input data are the same for the four GLDAS models, and Wang wen [24] and Wang Xiaoju [25] also found that precipitation and temperature data from the 4 GLDAS models were strongly consistent. The change trend, correlation coefficient (R), and root mean square error (RMSE) were used to determine the degree of similarity between the two datasets.

Water 2018, 10, 254

6 of 16

2.4.2. Comparison and Applicability Analysis of Runoff Data As the soil moisture input data for the four GLDAS models are different, we used all GLDAS data to compare runoff with observed values. The total runoff of the midstream of the Yarlung Zangbo River basin was calculated using the water balance equation by combining terrestrial water storage variations derived from the four GLDAS models (terrestrial water storage is represented by the sum of soil moisture and the snow water equivalent [20]) with precipitation and evapotranspiration. The regulation of distribution, the trend of variation, and continuity and mutability were analysed using the least square method. Subsequently, four statistical parameters (i.e., the correlation coefficient (R), mean bias (BIAS), standard deviation of the differences (σd ), and the ratio of the standard deviations (σr /σobs )) were calculated to compare the correlation and error between the GLDAS and observed runoff data. In addition, the Brunke ranking method [26] was applied to evaluate the data quality and applicability of the four datasets for the Yarlung Zangbo River basin. The water balance equation is as follows: Water 2018, 10, x FOR PEER REVIEW

6 of 15

P(t) − ET (t) − R(t) = ∆S(t) (3) where P is the areal mean rate of precipitation, ET is the areal mean rate of evapotranspiration, R is the total basin runoff, S represents terrestrial water storage and t is the time. where P is the areal mean rate of precipitation, ET is the arealobserved mean ratebyofGLDAS, evapotranspiration, R is the Figure 2 is a Schematic Diagram for the Water Balance Equation. In this paper, the variations are total basin runoff, S represents terrestrial water storage observed by GLDAS, and t is the time. the monthly values for the Yarlung Zangbo River basin. So, ΔS(t) is the water-storage variation Figure 2 is a Schematic Diagram for the Water Balance Equation. In this paper, the variations are between times t1 and 2 (in which the subscripts 1 and 2 refer to the beginning and the end of the the monthly values for tthe Yarlung Zangbo River basin. So, ∆S(t) is the water-storage variation between month). times t and t (in which the subscripts 1 and 2 refer to the beginning and the end of the month). 1

2

Δ S = S ( t 2 − t1 ) ∆S = S(t2 − t1 )

(4) (4) Because GLDAS can only provide 3 h and one month’s mean land surface state variables. We can only 3 h and month’s mean land surface state We usedBecause a brieflyGLDAS and efficient wayprovide to calculate the one ΔS following the method. The ΔS of avariables. month could used a briefly be denoted as:and efficient way to calculate the ∆S following the method. The ∆S of a month could be denoted as: 1 Δ ∆S S i i ≈≈ 1 ((SSi+i1+ 1−−Si−S1i)− 1 ) (5) (5) 22

Figure 2. A Schematic Diagram for the Water Balance Equation.

The equations of the four statistical parameters are as follows (6)–(9) [27]: The equations of the four statistical parameters are as follows (6)–(9) [27]: R

1 N

N

N ( M − M )( O − O = 1 1 (−MMi −) M)(1Oi − O( O)  N( Mi∑  =1 N sN R= s i

i

N

N

2

i

i=1

N

σ

r



obs

i

− O )

i=1

N 1 1 (NMi − N ∑ BIAS = i=1 Di

σ

)

i=1

d



M)

2

where

1 N

i =1

= (

1 N −1

 1  N −1  =  1

N



2

(6) (6)

N

∑ (Oi − O) iD =i1= M i − O i

2

(7)

(8)

D i2 ) 1 / 2

i =1

1/2

N



i=1 N

(M

i

 − M )2   1/2 

(9)

Water 2018, 10, 254

7 of 16

BI AS =

1 N Di N i∑ =1

where Di = Mi − Oi

1 N 2 σd = ( Di ) N − 1 i∑ =1 

1 N −1

σr /σobs = 

1/2

N

(8)

∑ ( Mi − M )

2

1/2

i =1

1 N −1

N

(7)

∑ (Oi − O)

2

1/2

(9)

i =1

where Mi is runoff calculated from the four models of GLDAS data by using the water balance equation; M is the mean runoff during the study period, based on the four models of GLDAS; Oi is the observed runoff from Nuxia gauging station; O is the mean runoff from Nuxia gauging station during the study period; and, N is the number of months during the study period. We used the Brunke ranking method [28] to rank the four datasets from 1 to 4. A ranking of 1 was assigned to the dataset with the lowest BIAS (or lowest σd , or highest R). A ranking of 4 was assigned to the dataset with the largest BIAS (or largest σd , or lowest R). The σr /σobs value closest to 1 was ranked 1, and the value furthest from 1 was ranked 4. The arithmetic mean of each dataset was calculated from the rankings of the four statistical parameters. We then obtained the total ranking scores of these datasets. 3. Results and Discussion 3.1. Characteristics of Observed Climate and Runoff Changes 3.1.1. Characteristics of Observed Precipitation Change during the Past 40 Years Precipitation changes in the midstream region of the Yarlung Zangbo River basin were analysed for three different time scales. Figure 3a–c show that mean monthly precipitation during the research period exhibited clear cyclical fluctuations of approximately 12 months. However, the peak values differed considerably between years. The mean monthly precipitation was 39.96 mm and fluctuated between 1.61 mm and 125.56 mm. The maximum monthly precipitation occurred in August 1998, whereas no precipitation was recorded in December 1996 and January 2016. In terms of secular mean monthly precipitation, the highest precipitation occurred in July (125.56 mm), and the lowest precipitation occurred in December (1.61 mm). In terms of interannual variation, the mean annual precipitation was 497.50 mm and fluctuated between 318.11 mm and 621.84 mm. The maximum monthly precipitation occurred in 2000, and the minimum monthly precipitation occurred in 1983. The trend of annual precipitation showed a slight increase of approximately 2.3 mm per year during the study period. Based on a 5-year moving average, the lowest mean annual precipitation was 386.53 mm, which occurred from 1981 to 1985, and the highest mean annual precipitation was 580.74 mm, which occurred from 1998 to 2002.

precipitation occurred in December (1.61 mm). In terms of interannual variation, the mean annual precipitation was 497.50 mm and fluctuated between 318.11 mm and 621.84 mm. The maximum monthly precipitation occurred in 2000, and the minimum monthly precipitation occurred in 1983. The trend of annual precipitation showed a slight increase of approximately 2.3 mm per year during the study period. Based on a 5-year moving average, the lowest mean annual precipitation was 386.53 Watermm, 2018,which 10, 254 occurred from 1981 to 1985, and the highest mean annual precipitation was 580.74 mm, 8 of 16 which occurred from 1998 to 2002.

(a)

(b)

(c)

Figure 3. Characteristics of precipitation change in the midstream region of the Yarlung Zangbo River

Figure 3. Characteristics of precipitation change in the midstream region of the Yarlung Zangbo River basin. (a) Time series of mean monthly precipitation; (b) The secular mean monthly precipitation; and, basin. (a) Time series of mean monthly precipitation; (b) The secular mean monthly precipitation; and, (c) The mean annual precipitation (c) The mean annual precipitation

3.1.2. Characteristics of Observed Temperature Change during the Past 40 Years

3.1.2. Characteristics of Observed Temperature Change during the Past 40 Years

Figure 4a shows that the mean monthly temperature during the research period also presented 4a shows that the similar mean monthly temperature duringmean the research period aFigure clear cyclical fluctuation, to precipitation. The monthly temperature wasalso 5.16presented °C and a ◦ C and clearfluctuated cyclical fluctuation, precipitation. The monthly temperature 5.162009, between −6.82similar °C andto 15.57 °C. The maximum monthly mean temperature occurredwas in July ◦ ◦ and thebetween minimum−monthly temperature occurred in January 1983. Intemperature terms of the occurred secular mean fluctuated 6.82 C and 15.57 C. The maximum monthly in July the highest temperatureoccurred occurredininJanuary July (13.57 °C) In and the lowest January 2009,monthly and thetemperature, minimum monthly temperature 1983. terms of the in secular mean ◦ (−3.88 °C). In terms of interannual variation, the mean annual temperature of the study area was monthly temperature, the highest temperature occurred in July (13.57 C) and the lowest in January 5.61◦°C fluctuated between 4.45variation, and 6.85 mm. The maximum monthly temperature occurred (−3.88 C).and In terms of interannual the mean annual temperature of the study areainwas 2009, and the minimum monthly temperature occurred in 1983. The annual temperature trend 5.61 ◦ C and fluctuated between 4.45 and 6.85 mm. The maximum monthly temperature occurred in showed a slight increase of approximately 0.035 °C per year. Based on a five-year moving average, it 2009, and the minimum monthly temperature occurred in 1983. The annual temperature trend showed can be concluded that the lowest mean annual temperature was 4.67 °C, which occurred from 1979 a slight increase of approximately 0.035 ◦ C per year. Based on a five-year moving average, it can be to 1983, and the highest mean annual temperature was 6.49 °C, which occurred from 2006 to 2010.

concluded that the lowest mean annual temperature was 4.67 ◦ C, which occurred from 1979 to 1983, and Water the highest annual temperature was 6.49 ◦ C, which occurred from 2006 to 2010. 8 of 15 2018, 10, xmean FOR PEER REVIEW

(a)

(b)

(c)

Figure 4. Characteristics of temperature changes in the midstream region of the Yarlung Zangbo River

Figure 4. Characteristics of temperature changes in the midstream region of the Yarlung Zangbo River basin. (a) Time series of monthly temperature; (b) The secular mean monthly temperature; and, (c) basin. (a) Time series of monthly temperature; (b) The secular mean monthly temperature; and, (c) The The mean annual temperature. mean annual temperature.

3.1.3. Characteristics of Observed Runoff Change during the Past 20 Years We used the observed monthly runoff data of Nuxia gauging station from January 1979 to December 2000 to validate the results based on the GLDAS. Figure 5a–c show observed runoff for the station for three time scales. It can be concluded that mean monthly runoff presented a clear cyclical fluctuation. Mean monthly runoff during the research period was 25.01 mm and fluctuated between

(b)

(c)

Figure 4. Characteristics of temperature changes in the midstream region of the Yarlung Zangbo River basin. Water 2018, 10,(a) 254Time series of monthly temperature; (b) The secular mean monthly temperature; and, (c)9 of 16 The mean annual temperature.

3.1.3. 3.1.3. Characteristics Characteristics of of Observed ObservedRunoff RunoffChange Changeduring duringthe thePast Past20 20Years Years We We used used the the observed observed monthly monthly runoff runoff data data of of Nuxia Nuxia gauging gauging station station from from January January 1979 1979 to to December December2000 2000to tovalidate validatethe theresults resultsbased basedon onthe the GLDAS. GLDAS. Figure Figure5a–c 5a–cshow show observed observedrunoff runofffor forthe the station station for forthree threetime timescales. scales. ItIt can can be be concluded concluded that that mean mean monthly monthly runoff runoff presented presentedaaclear clearcyclical cyclical fluctuation. fluctuation. Mean Mean monthly monthly runoff runoff during during the the research researchperiod periodwas was25.01 25.01mm mmand andfluctuated fluctuatedbetween between 3.07 3.07 and and 142.62 142.62 mm. mm. The The maximum maximum monthly monthly runoff runoff occurred occurred in in August August 1998, 1998, and and the the minimum minimum monthly 1997 (3.07 mm). In In terms of the secular mean monthly runoff, the monthlyrunoff runoffoccurred occurredininFebruary February 1997 (3.07 mm). terms of the secular mean monthly runoff, highest runoff occurred in August (70.19(70.19 mm),mm), and the runoffrunoff occurred in February (5.00 mm). the highest runoff occurred in August andlowest the lowest occurred in February (5.00 In terms interannual variation, the mean annual runoffrunoff of Nuxia gauging station was 300.08 mm mm). In of terms of interannual variation, the mean annual of Nuxia gauging station was 300.08 and between 201.41201.41 mm and mm. The annual runoffrunoff occurred in 2000, mmfluctuated and fluctuated between mm474.52 and 474.52 mm.maximum The maximum annual occurred in and minimum runoffrunoff occurred in 1983. The trend of annual runoff showed a slight increase of 2000,the and the minimum occurred in 1983. The trend of annual runoff showed a slight increase 0.035 mmmm perper year. This result is is consistent with other of 0.035 year. This result consistent with otherresearch research[29]. [29].Based Basedon onaafive-year five-yearmoving moving average, average,the the lowest lowest runoff runoffwas was245.46 245.46 mm, mm, which which occurred occurredfrom from1982 1982to to1986, 1986, and and the the highest highest runoff runoff was 375.10 mm, which occurred from 1996 to 2000. was 375.10 mm, which occurred from 1996 to 2000.

(a)

(b)

(c)

Figure 5.5. Characteristics Characteristics of of observed observed runoff runoff change change for for Nuxia Nuxia gauging gauging station. station. (a) (a) Time Time series series of of Figure monthlyrunoff; runoff;(b) (b)The Thesecular secularmean meanmonthly monthlyrunoff; runoff;and, and,(c) (c)The Themean meanannual annualrunoff. runoff. monthly

3.2. Applicability Analysis of GLDAS Data 3.2.1. Applicability Analysis of GLDAS Precipitation Data Because all four GLDAS data models have the same precipitation and temperature input data, we used GLDAS-MOS model data for comparison with observed data. Figure 6 indicates that GLDAS mean monthly precipitation during the research period presented the same clear cyclical fluctuation as that for observed precipitation. However, the amplitude of observed precipitation was significantly higher than the GLDAS precipitation. A root mean square error of 34.35 mm was recorded between the GLDAS and observed mean monthly precipitation from 1979 to 2016. Correlation analysis of the two data sets (Figure 7) indicated that GLDAS mean monthly precipitation from 1979 to 2016 was moderately correlated with observed precipitation, and the correlation coefficient was 0.75. Strong correlations were observed between GLDAS and observed precipitation from January to June and from September to December, with respective correlation coefficients of 0.98 and 0.96. However, there were big differences between GLDAS and observed precipitation in July and August; the mean monthly observed precipitation was higher than the GLDAS precipitation, and the correlation coefficient was very low. The above analysis shows that GLDAS precipitation data have a good applicability in the

from 1979 to 2016 was moderately correlated with observed precipitation, and the correlation from 1979 to 2016 was moderately correlated with observed precipitation, and the correlation coefficient was 0.75. Strong correlations were observed between GLDAS and observed precipitation coefficient was 0.75. Strong correlations were observed between GLDAS and observed precipitation from January to June and from September to December, with respective correlation coefficients of from January to June and from September to December, with respective correlation coefficients of 0.98 and 0.96. However, there were big differences between GLDAS and observed precipitation in 0.98 and 0.96. However, there were big differences between GLDAS and observed precipitation in July and August; the mean monthly observed precipitation was higher than the GLDAS precipitation, Water 10, 254 10 of 16 July 2018, and August; the mean monthly observed precipitation was higher than the GLDAS precipitation, and the correlation coefficient was very low. The above analysis shows that GLDAS precipitation and the correlation coefficient was very low. The above analysis shows that GLDAS precipitation data have a good applicability in the Yarlung Zangbo River basin within a certain time period. data have a good applicability in the Yarlung Zangbo River basin within a certain time period. Yarlung Zangbo River basin a certain time period. the quality of the the follow-up July and August However, the quality of thewithin July and August data needs However, to be improved. So, in study, However, the quality of the July and August data needs to be improved. So, in the follow-up study, data needs to be measures improved.toSo, in the follow-up study, we should take measures improve calculation we should take improve calculation method and mothod to extracttodata to improve data we should take measures to improve calculation method and mothod to extract data to improve data method and precision mothod toofextract data to improve data quality and precision of GLDAS. quality and GLDAS. quality and precision of GLDAS.

Figure 6. 6. GLDAS and and observed precipitation precipitation changes from from 1979 to to 2016. Figure Figure 6. GLDAS GLDAS and observed observed precipitation changes changes from1979 1979 to2016. 2016.

Figure 7. Correlation coefficient between Global Land Data Assimilation System (GLDAS) and Figure 7.7. Correlation Correlation coefficient coefficient between between Global Global Land Land Data Data Assimilation Assimilation System System (GLDAS) (GLDAS) and and Figure observed precipitation data. . observed precipitation data. . observed precipitation data. .

3.2.2. Applicability Analysis of GLDAS Temperature Data 3.2.2. Applicability Applicability Analysis Analysis of of GLDAS GLDAS Temperature TemperatureData Data 3.2.2. Figure 8 shows the GLDAS and observed temperature change from 1979 to 2016. It can be Figure 88 shows shows the the GLDAS GLDAS and and observed observed temperature temperature change change from from 1979 1979 to to2016. 2016. It It can can be be Figure concluded that the GLDAS mean monthly temperature during the research period presented the concluded that the GLDAS mean monthly temperature during the research period presented the concluded that thefluctuation GLDAS mean monthly temperature during the research periodofpresented the same clear cyclical as the observed temperature. However, the amplitude the observed same clear clear cyclical cyclical fluctuation fluctuationas as the the observed observed temperature. temperature. However, However,the theamplitude amplitudeof ofthe the observed observed same temperature was significantly lower than that of the GLDAS temperature. The root mean square error temperaturewas wassignificantly significantlylower lowerthan thanthat thatof ofthe the GLDAS GLDAS temperature. temperature. The The root mean mean square square error error temperature between GLDAS and observed mean monthly temperature from 1979 toroot 2016 was 11.57 °C. between GLDAS and observed mean monthly temperature from 1979 to 2016 was 11.57 °C. ◦ between GLDAS and observed mean monthly temperature from 1979that to 2016 11.57 mean C. Correlation Correlation analysis of the two data sets (Figure 9) indicated the was GLDAS monthly Correlation analysis of the two data sets (Figure 9) indicated that the GLDAS mean monthly analysis of the two1979 datato sets (Figure 9) indicated that with the GLDAS mean monthly temperature from temperature from 2016 was highly correlated observed temperature, with a correlation temperature from 1979 to 2016 was highly correlated with observed temperature, with a correlation 1979 to 2016 was highly correlated withGLDAS observed with a correlation coefficient of coefficient of 0.94. This result shows that datatemperature, can reproduce observed temperature changes coefficient of 0.94. This that result showsdata that can GLDAS data can reproduce observed changes temperature changes 0.94. This result shows GLDAS reproduce observed temperature in the study period. However, there was a big difference between the GLDAS and the observed temperature from November to February of the next year. The observed mean monthly temperature was higher than the GLDAS temperature, and the correlation coefficient was very low. The above analysis shows that GLDAS temperature data have a good applicability to the midstream region of the Yarlung Zangbo River basin within a certain time period. However, the quality of data from November to February of the next year needs to be improved. We also should take measures to improve data quality and precision of GLDAS in the follow-up study.

temperature temperature from from November November to to February February of of the the next next year. year. The The observed observed mean mean monthly monthly temperature temperature was higher than the GLDAS temperature, and the correlation coefficient was very low. was higher than the GLDAS temperature, and the correlation coefficient was very low. The The above above analysis shows that GLDAS temperature data have a good applicability to the midstream analysis shows that GLDAS temperature data have a good applicability to the midstream region region of of the Yarlung Zangbo River basin within a certain time period. However, the quality of data from the Yarlung Zangbo River basin within a certain time period. However, the quality of data from November to February of the next year needs to be improved. We also should take measures to November to Water 2018, 10, to 254 February of the next year needs to be improved. We also should take measures 11 of 16 improve improve data data quality quality and and precision precision of of GLDAS GLDAS in in the the follow-up follow-up study. study.

Figure 8. observed temperature changes from 1979 to 2016. Figure Figure 8. 8. GLDAS GLDAS and and observed observed temperature temperature changes changes from from1979 1979to to2016. 2016.

Figure 9. Correlation coefficient between GLDAS and observed temperature data. Figure 9. 9. Correlation Correlation coefficient coefficient between between GLDAS GLDAS and and observed observed temperature temperaturedata. data. Figure

3.2.3. Applicability Analysis of GLDAS Runoff Data 3.2.3. Applicability Applicability Analysis Analysis of of GLDAS GLDAS Runoff Runoff Data Data 3.2.3. Figure 10 presents the time series of mean monthly runoff based on the GLDAS models and Figure 10 10 presents presents the the time time series series of of mean mean monthly monthly runoff runoff based based on on the the 444 GLDAS GLDAS models models and and Figure the observed runoff for Nuxia gauging station. It must point that the runoff values calculated in this the observed runoff for Nuxia gauging station. It must point that the runoff values calculated in this the observed runoff for Nuxia gauging station. It must point that the runoff values calculated in this paper are relative values relative to secular mean monthly runoff. Hence, there exists negative runoff paper are arerelative relativevalues valuesrelative relativeto tosecular secularmean mean monthly monthly runoff. runoff. Hence, Hence, there there exists exists negative negative runoff runoff paper in the text and Figure 10. It can observed that runoff from the 44 GLDAS models presented similar in the text and Figure 10. It can observed that runoff from the GLDAS models presented similar in the text and Figure 10. It can observed that runoff from the 4 GLDAS models presented similar fluctuations and significant seasonal variability when compared with observed runoff. However, fluctuations and and significant significant seasonal seasonal variability variability when when compared compared with with observed observed runoff. runoff. However, However, fluctuations there were also some differences. The peak runoff values from the four GLDAS models were in good there were also some differences. The peak runoff values from the four GLDAS models were in good good there were also some differences. The peak runoff values from the four GLDAS models were in agreement, and the runoff amplitude of each GLDAS model was higher than the observed runoff. To agreement, and the runoff amplitude of each GLDAS model was higher than the observed runoff. To agreement, and the runoff amplitude of each GLDAS model was higher than the observed runoff. quantitatively describe runoff variations, we adopted the least square method to fit the annual quantitatively describe runoff variations, weweadopted the annual annual To quantitatively describe runoff variations, adoptedthe theleast leastsquare squaremethod method to to fit fit the amplitude, cycle, phase and linear trend of the total runoff time series for each dataset. The result amplitude, cycle, phase and linear trend of the total runoff time series for each dataset. The result amplitude, cycle, phase and linear trend of the total runoff time series for each dataset. The result (Table 1) shows that all 44 datasets had similar change cycles of approximately 12 months compared (Table 1) shows that all datasets had similar change cycles of approximately 12 months compared (Table 1) shows that all 4 datasets had similar change cycles of approximately 12 months compared with the observed runoff for Nuxia gauging station. The GLDAS-MOS based result had the highest with the the observed observed runoff runofffor forNuxia Nuxiagauging gaugingstation. station. The The GLDAS-MOS GLDAS-MOS based based result result had had the the highest highest with annual amplitude of 91.91 mm water equivalent. In addition, the GLDAS model-based runoff values annual amplitude of 91.91 mm water equivalent. In addition, the GLDAS model-based runoff values annual amplitude of 91.91 mm water equivalent. In addition, the GLDAS model-based runoff values were considerably different from the observed runoff in terms of the phase. Observed runoff occurred were considerably different from the observed runoff in terms of the phase. Observed runoff occurred were considerably different from the observed runoff in terms of the phase. Observed runoff occurred later than GLDAS-based runoff. The reason for this off that result observed runoff presented is later than thanGLDAS-based GLDAS-basedrunoff. runoff.The Thereason reasonfor forthis thisoff offisis isthat thatresult resultofof ofobserved observedrunoff runoff presented later presented is is a aa value after aa series of complex process, namely rainfall, evapotranspiration, surface and subsurface value after series complexprocess, process,namely namelyrainfall, rainfall,evapotranspiration, evapotranspiration,surface surfaceand andsubsurface subsurface value after a series ofofcomplex runoff, infiltration, al. So, will take long time for all forms of water to converge to the exit section runoff,infiltration, infiltration,etet etal. al.So, So,itit it will will take take aaa long long time time for for all all forms forms of of water water to to converge converge to to the the exit exit section section runoff, of the river basin. However, result of GLDAS presented is relativly simpler. The terrestrial water of the river basin. However, result of GLDAS presented is relativly simpler. The terrestrial water of the river basin. However, result of GLDAS presented is relativly simpler. The terrestrial water storage change (ΔS(t)) was represented only by the sum of soil moisture change and the snow water storage change (ΔS(t)) was represented only by the sum of soil moisture change and the snow water storage change (∆S(t)) was represented only by the sum of soil moisture change and the snow water equivalent change. Hence, there exists one or two month’s off when compared with the observed equivalent change. change. Hence, exists one one or or two two month’s month’s off off when when compared compared with with the the observed observed equivalent Hence, there there exists runoff. With respect to the linear trend, runoff from the Yarlung Zangbo River basin increased runoff. With respect to the linear trend, runoff from the Yarlung Zangbo River basin increased runoff. With respect to the linear trend, runoff from the Yarlung Zangbo River basin increased slightly slightly during the period for all slightlythe during the research research period fordatasets. all five five datasets. datasets. during research period for all five

Table 2 presents four statistical parameters between the 4 datasets and observed runoff, as well as average ranking scores from these datasets. The results show that GLDAS-Noah-based products best represented time series of observed runoff, except for BIAS. The total Brunke ranking scores of the four datasets were 2.75, 1.75, 2.50, and 3.00, respectively. Hence, the result shows that GLDASNoah-based runoff data were closer to observed runoff compared with the three other GLDAS Water 2018, 10, 254 12 of 16 models.

Figure 10. 10. Time Time series series of of monthly monthly runoff runoff calculated calculated from from GLDAS GLDAS and and observed observedrunoff runoffdata. data. Figure Table1.1.Annual Annualamplitude, amplitude,cycle, cycle,phase, phase,and and linear linear trend trend of of runoff runoff based based on on the the least least square square method. method. Table

Dataset Amplitude (mm) Cycle (month) Phase (°) Amplitude (mm) Cycle (month) Phase (◦ ) GLDAS-MOS 91.91 12.00 146 GLDAS-MOS 91.91 12.0012.01 146108 GLDAS-Noah 67.11 GLDAS-Noah 67.11 12.01 108 GLDAS-VIC 77.92 GLDAS-VIC 77.92 11.9911.99 151151 GLDAS-CLM 83.19 GLDAS-CLM 83.19 12.0012.00 107107 Observed runoff 28.81 12.0012.00 220220 Observed runoff 28.81 Dataset

Linear Trend (mm/a) Linear Trend (mm/a) 0.44 0.44 0.33 0.33 0.26 0.26 0.11 0.11 0.44 0.44

Table 2.four The average ranking scores of runoff data the four models. Table 2 presents statistical parameters between the 4from datasets andGLDAS observed runoff, as well as average ranking scores from these datasets. The results show that GLDAS-Noah-based products best Name GLDAS-MOS GLDAS-Noah GLDAS-VIC GLDAS-CLM represented timeRseries of observed0.70/3 runoff, except for0.74/1 BIAS. The total Brunke 0.69/4 ranking scores 0.70/2of the four datasets were BIAS 2.75, 1.75, 2.50, and 18.82/1 3.00, respectively. 19.53/4 Hence, the result 18.87/2 shows that GLDAS-Noah-based 18.96/3 runoff data were runoff compared with the three other GLDAS models. σdcloser to observed 26.13/3 21.41/1 23.53/2 30.92/4

σr/σobs 2.78/4 2.05/1 2.38/2 2.66/3 Table 2. The average ranking scores of runoff data from the four GLDAS models. Total ranking scores 2.75 1.75 2.50 3.00 Number a/number b: a represents the values of 4 statistical parameters; b represents the total ranking Name GLDAS-MOS GLDAS-Noah GLDAS-VIC GLDAS-CLM scores. R 0.70/3 0.74/1 0.69/4 0.70/2 BIAS 18.82/1 19.53/4 18.87/2 18.96/3 3.3. Impact of Climate Change on Runoff σd 26.13/3 21.41/1 23.53/2 30.92/4 σ /σ 2.78/4 2.05/1 2.38/2 2.66/3 r obs Climate change affects hydrological cycles locally and globally [30]. From precipitation and Total ranking scores 2.75 1.75 2.50 3.00

temperature and runoff curves under three time scales, from 1979 to 2000, over midstream of the Number a/number b: a represents the values of 4 statistical parameters; b represents the total ranking scores. Yarlung Zangbo river basin. It can be concluded that changes in runoff were very similar to precipitation and temperature changes: the high flow period was consistent with high precipitation 3.3. of Climate Change on Runoff andImpact temperature, and the low flow period was consistent with low precipitation and temperature. July Climate and August were the peak periods for mean monthly rainfall, [30]. temperature and runoff, and and change affects hydrological cycles locally and globally From precipitation December and January exhibited the lowest mean monthly rainfall, temperature and runoff each temperature and runoff curves under three time scales, from 1979 to 2000, over midstream of year. At the Zangbo annual time of runoff was consistent lowest mean the Yarlung river scale, basin.the It lowest can be amount concluded that changes in runoffwith werethe very similar to annual precipitation and temperature. For instance, in the midstream region of the Yarlung Zangbo precipitation and temperature changes: the high flow period was consistent with high precipitation River basin, 1983 was yearperiod duringwas the research period, a mean annual runoff of 201.41 and temperature, andthe thedriest low flow consistent with with low precipitation and temperature. mm,and mean annualwere precipitation 318.11for mm, andmonthly a mean annual 4.45runoff, °C. Onand the July August the peak of periods mean rainfall,temperature temperatureofand contrary, the year with the highest amount of runoff had the highest mean annual precipitation and December and January exhibited the lowest mean monthly rainfall, temperature and runoff each year. temperature. 1998 wasamount the wettest year in this region during thelowest research period; the At the annual For timeinstance, scale, the lowest of runoff was consistent with the mean annual precipitation and temperature. For instance, in the midstream region of the Yarlung Zangbo River basin, 1983 was the driest year during the research period, with a mean annual runoff of 201.41 mm, mean annual precipitation of 318.11 mm, and a mean annual temperature of 4.45 ◦ C. On the contrary, the year with the highest amount of runoff had the highest mean annual precipitation and temperature. For instance, 1998 was the wettest year in this region during the research period; the mean annual runoff

Water 2018, 10, 254 Water 2018, 10, x FOR PEER REVIEW

13 of 16 12 of 15

was 459.58 mm, the mean annual mm, precipitation 625.17precipitation mm, and thewas mean annual temperature was mean annual runoff was 459.58 the meanwas annual 625.17 mm, and the mean ◦ C. 6.34 annual temperature was 6.34 °C. A relationship was determined between precipitation, temperature, and runoff. Figure 11 shows the correlation between monthly, annual, and secular mean monthly climate factors and runoff, and Table 33 shows shows the the corresponding corresponding correlation correlation coefficients coefficients at at different different time time scales. scales. The correlation precipitation and runoff at monthly, annual and secular mean monthly scales coefficients between precipitation were 0.752, 0.897 and 0.925, respectively, and all three passed the significance significance test test at at the the 0.001 0.001 level. level. The correlation coefficients between temperature and runoff at the three time scales were were 0.704, 0.704, 0.499, 0.499, and 0.806, and all three also passed the significance test at the 0.001 level. These results indicate that the correlation coefficients between precipitation and runoff at the three time scales were higher than the correlation coefficients between temperature and runoff. Secular mean monthly precipitation had mean monthly runoff. ThisThis means that that precipitation is theismost the strongest strongestrelationship relationshipwith withsecular secular mean monthly runoff. means precipitation the important factor that affects natural runoff changes in the midstream region of the Yarlung Zangbo most important factor that affects natural runoff changes in the midstream region of the Yarlung River basin, and rainfall control evolution of the river to some Zangbo River basin, andcan rainfall canthe control the evolution of the river degree. to some degree.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 11. Correlation Correlationanalysis analysisbetween between climate factors runoff 3 different time scales. (a,b) Figure 11. climate factors andand runoff for 3for different time scales. (a,b) Mean Mean monthly values; (c,d) Mean annual values; and, (e,f) Secular mean monthly values. monthly values; (c,d) Mean annual values; and, (e,f) Secular mean monthly values. Table and runoff runoff data data at at different different scales. scales. Table 3. 3. Correlation Correlation analysis analysis between between climate climate and Data Data Mean monthly values Mean values Meanmonthly annual values Mean annual values Secular mean monthly values Secular mean monthly values

Correlation Coefficient between Correlation Coefficient between Precipitation and Runoff Precipitation and Runoff 0.752 0.752 0.897 0.897 0.925 0.925

Correlation Coefficient between Correlation Coefficient between Temperature and Runoff Temperature and Runoff 0.704 0.704 0.499 0.499 0.806 0.806

4. Conclusions 4. Conclusions This paper used observed data to analyse the characteristics of temperature, precipitation and This paperin used data to analyse the characteristics of temperature, precipitation and runoff changes the observed midstream region of the Yarlung Zangbo River basin. In order to enrich data runoff changes the midstream region the Yarlung Zangbo River basin.alternative In order toclimate enrich data material for thisindata-limited basin, theofGLDAS was selected to provide and hydrological data, and expect to provide reference for the fellow-up studies. The main results of this research are as follows:

Water 2018, 10, 254

14 of 16

material for this data-limited basin, the GLDAS was selected to provide alternative climate and hydrological data, and expect to provide reference for the fellow-up studies. The main results of this research are as follows: (1)

(2)

(3)

(4)

Under the background of global warming, precipitation, temperature, and runoff changed significantly, and showed certain regularity during the research period. Based on the three time scales, we observed that mean monthly precipitation, temperature and runoff presented clear cyclical fluctuations of approximately 12 months. The highest secular mean monthly values occurred in July and August, and the lowest values occurred in December and January. The annual precipitation, temperature, and runoff trends showed slight increases during the study period. The maximum mean annual values occurred in 1998 and 1999, and the minimum mean annual values occurred in 1983. These analyses indicate that changes in precipitation, temperature, and runoff during the research period exhibited strong consistency. A comparison of observed precipitation and temperature with GLDAS climate data indicated that the GLDAS data have good applicability in the midstream region of the Yarlung Zangbo River basin. The GLDAS mean monthly precipitation and temperature presented the same clear cyclical fluctuation as the observed climate values. The GLDAS mean monthly precipitation was moderately correlated with observed precipitation, with a correlation coefficient of 0.75. The GLDAS mean monthly temperature was highly correlated with observed temperature, with a correlation coefficient of 0.94. This shows that GLDAS temperature data can be used to reproduce observed temperature in the study period. However, the quality of precipitation data from July to August and of temperature data from November to February of the next year needs to be improved. A comparison of observed and GLDAS-based runoff data indicated that runoff from the four GLDAS models had similar fluctuations and significant seasonal variability compared with observed runoff during the study period. All four datasets had similar change cycles of approximately 12 months compared with observed runoff for Nuxia gauging station, and runoff from the Yarlung Zangbo River basin increased slightly during the study period. However, the amplitudes of each GLDAS model were higher than those of observed runoff, and observed runoff occurred later than GLDAS-based runoff. The Brunke ranking method showed that the total ranking scores of the four GLDAS models were 2.75, 1.75, 2.50, and 3.00, respectively. Hence, it can be concluded that GLDAS-Noah-based runoff data are closer to observed runoff data for Nuxia gauging station compared with the three other GLDAS models. Runoff from the Yarlung Zangbo River basin is strongly affected by climate change. Correlation coefficients between precipitation and runoff for the 3 time scales were higher than those between temperature and runoff. Secular mean monthly precipitation had the strongest relationship with runoff. This means that a change in precipitation is the most important factor that affects natural runoff changes in the midstream region of the Yarlung Zangbo River basin. Rainfall can control the evolution of the river to some degree.

Based on observed and GLDAS climate and hydrology data, the relationship between runoff and climate change characteristics, the applicability of GLDAS data and the impact of climate factors on runoff data were quantitatively analysed. In conclusion, under the background of global warming, the climate and hydrological elements in the study period showed significant regularity changes. Rainfall is the main factor that impacts runoff change in this area. To the whole study period, GLDAS is not a good system to describe the land surface conditions of the Yarlung Zangbo River basin. However, within a certain time period, GLDAS data have a good applicability in the Yarlung Zangbo River basin. To some extent, it could be a good supplement of climate and hydrology data to data limited areas. We need to take measures to improve data quality and the precision of GLDAS in the follow-up study. However, some problems exist. The Yarlung Zangbo River is an important international river that has the highest mean altitude in China. The number of climate and hydrological gauging stations

Water 2018, 10, 254

15 of 16

is limited because of the complex terrain and environment. Certain results and the applicability of GLDAS in this period are limited because the time series are short and the data quality is poor. Hence, the selection of an alternative dataset specifically aimed at the Yarlung Zangbo River basin has certain limitations. This problem results in a degree of uncertainty in this study. In addition, with increasing human interference in the river, how to reduce the impacts of the underlying surface and other human activities will become the focus of future work when discussing the impact of climate change on runoff. Acknowledgments: This work was funded by the National Natural Science Foundation of China (No. 91547211) and the National Key R&D Program of China (No. 2016YFC0502200). The authors also would like to thank China meteorological data service center providing meteorological data and the GES DISC providing GLDAS data. The work of these scientists is the basis for our investigations. Author Contributions: Hong Zhang, Ling Lei Zhang and Jia Li conceived and designed the study; Rui Dong An provided data and analysis tools; Hong Zhang, Ling Lei Zhang and Jia Li analyzed the data; Hong Zhang wrote the paper; and all of the authors polished the writing. Conflicts of Interest: The authors declare no conflict of interest.

References 1. 2. 3. 4. 5. 6.

7. 8. 9. 10. 11. 12. 13. 14. 15. 16.

17.

Guo, J. Research on the Impact of Climate Change on Watershed Water Cycle and Water Resource; Wuhan University: Wuhan, China, 2010. IPCC (Intergovernmental Panel on Climate Change). Synthesis Report; Cambridge University Press: Cambridge, UK, 2001. IPCC (Intergovernmental Panel on Climate Change). Climate Change: The IPCC Impacts Assessment; AGPS Press: Canberra, Australia, 2001. IPCC (Intergovernmental Panel on Climate Change). Impacts, Adaptations and Mitigation of Climate Change: Scientific Technical Analyses; Cambridge University Press: Cambridge, UK, 1995. IPCC (Intergovernmental Panel on Climate Change). Synthesis Report; Cambridge University Press: Cambridge, UK, 2007. Animesh, K.G.; Gain, C.; Fabrice, G. Renaud Climate Change Adaptation and Vulnerability Assessment of Water Resources Systems in Developing Countries: A Generalized Framework and a Feasibility Study in Bangladesh. Water 2012, 4, 345–366. [CrossRef] Lisa, J.L.; Joel, D.G.; Vernon, L.S. Climate Change and Classic Maya Water Management. Water 2011, 3, 479–494. [CrossRef] National Research Council. Climate, Climate Change and Water Supply; National Academy Press: Washington, DC, USA, 1977. Dong, L.H.; Xiong, L.H.; Yu, K.X.; Li, S. Research advances in effects of climate change and human activities on hydrology. Adv. Water Sci. 2012, 23, 278–285. Zhao, X.S.; Zhou, H.; Wang, W.C. Research Progresses of Impacts of the Climate Change on the Regional Hydrological Cycle System. J. North China Inst. Water Conserv. Hydroelectr. Power 2012, 33, 46–49. Zhang, L.P.; Chen, X.F.; Zhao, Z.P.; Hu, Z. Progress in Study of Climate Change Impacts on Hydrology and Water Resources. Prog. Geogr. 2008, 27, 60–67. Wang, G.Q.; Zhang, J.Y.; Liu, J.F.; Li, Y.; He, R.M. Review of impacts of climate change on hydrology and water resources. Adv. Water Sci. 2008, 2, 47–51. Bai, L.Y.; Rong, Y.S. Impacts of climate on water resources in source regions of Yangtze River and Yellow River. Water Resour. Prot. 2012, 28, 46–50. Tang, F.F.; Xun, Z.X.; Zuo, D.P. Response of Runoff to Climate Change in the Upper Yellow River Basin. Resour. Sci. 2012, 34, 1079–1088. Wang, G.Q.; Li, M.; Jin, J.L.; Liu, C. Variation Trend of Runoff in Fujiang River Catchment and Its Responses to Climate Change. J. China Hydrol. 2012, 32, 22–28. Jiao, K.Q.; Ye, B.S.; Han, T.D.; Jin, Z.; Yang, H.A. Response of Runoff to Climate Change in the Glacier NO.1 at the Headwater of Urumpi River, Tianshan Mountains during 1980–2006. J. Glaciol. Geocryol. 2011, 33, 606–611. Li, X.; Gao, Y.H.; Wang, W.Z.; Lan, Y.; Xu, J.; Li, K. Climate changes and applicability of GLDAS in the headwater of the Yellow River Basin. Adv. Earth Sci. 2014, 29, 531–540.

Water 2018, 10, 254

18. 19. 20. 21.

22. 23.

24. 25. 26. 27. 28. 29. 30.

16 of 16

Nie, N.; Zhang, W.C.; Deng, C. Spatial and Temporal Climate Variations from 1978 to 2009 and their Trend Projection over the Yarlung Zangbo River Basin. J. Glaciol. Geocryol. 2012, 34, 64–71. Li, Z.L.; Xu, Z.X.; Gong, T.L. Characteristics of runoff variationin the Yarlung Zangbo River basin. Geogr. Res. 2008, 27, 353–361. Xu, P.K. Study on Glacier&Groundwater Mass Changes of Yarlung Zangbo River Basin in the Recent Years; Nanjing University: Nanjing, China, 2013. Fang, H.L.; Beaudoing, H.K.; Rodell, M.; Teng, W.; Vollmer, B.E. Global Land Data Assimilation System (GLDAS) products, services and application from NASA Hydrology Data and Information Services Center (HDISC). In Proceedings of the ASPRS 2009 Annual Conference, Baltimore, MD, USA, 1–9 March 2009. Wang, Y.J.; Wang, S.D.; Song, W.L.; Yang, S. Application of GLDAS data to the potential evapotranspiration monitoring in Weihe River Basin. J. Arid Land Resour. Environ. 2013, 27, 54–58. Rodell, M.; Houser, P.R.; Jambo, R.U.; Gottschalck, J.; Mitchell, K.; Meng, C.; Arsenault, K.; Cosgrove, B.; Radakovich, J.; Bosilovich, M.; et al. The global land data assimilation system. Bull. Am. Meteorol. Soc. 2004, 85, 381–394. [CrossRef] Wang, W.; Wang, X.J.; Wang, P. Assessing the applicability of GLDAS monthly precipitation data in China. Adv. Water Sci. 2014, 25, 769–778. Wang, X.J.; Wang, W.; Wang, P. Assessing Applicability of GLDAS Monthly Temperature Data in China. Water Res. Power 2014, 32, 10–14. Brunke, M.A.; Fairall, C.W.; Zeng, X.B.; Eymard, L.; Curry, J.A. Which bulk aerodynamic algorithms are least problematic incomputing ocean surface turbulent fluxes. J. Clim. 2003, 16, 619–635. [CrossRef] Wei, F.Y. Modern Climatic Statistical Diagnosis and Prediction Technology; Climate Press: London, UK, 2007; pp. 1–298. Liu, C.; Ye, Y.; Jie, J.; Zhou, X.; Li, J.L.; Ge, J. Applicability of Soil Temperature and Moisture in Several Datasets over Qinghai-Xizang Plateau. Meteorology 2015, 34, 653–665. Jiang, C.; Zhang, L.B.; Li, D.Q.; Li, F. Water Discharge and Sediment Load Changes in China: Change Patterns, Causes, and Implications. Water 2015, 7, 5849–5875. [CrossRef] Zhu, T.J.; Ringler, C. Climate Change Impacts on Water Availability and Use in the Limpopo River Basin. Water 2012, 4, 63–84. [CrossRef] © 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).