Next Article in Journal
A Comparison of Three Sediment Acoustic Models Using Bayesian Inversion and Model Selection Techniques
Next Article in Special Issue
Satellite-Based Evapotranspiration in Hydrological Model Calibration
Previous Article in Journal
Editorial for the Special Issue “Microwave Indices from Active and Passive Sensors for Remote Sensing Applications”
Previous Article in Special Issue
Long-Term Dynamic of Poyang Lake Surface Water: A Mapping Work Based on the Google Earth Engine Cloud Platform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Climate Variabilities and Human Activities on Surface Water Extents in Reservoirs of Yongding River Basin, China, from 1985 to 2016 Based on Landsat Observations and Time Series Analysis

1
College of Resource Environment and Tourism, Capital Normal University, Beijing 100048, China
2
China Institute of Water Resources and Hydropower Research, Beijing 100038, China
3
Laboratory Cultivation Base of Environment Process and Digital Simulation, Beijing 100048, China
4
Atmospheric Sciences and Global Change Division, Pacific Northwest National Laboratory, Richland, WA 99352, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(5), 560; https://doi.org/10.3390/rs11050560
Submission received: 28 January 2019 / Revised: 21 February 2019 / Accepted: 1 March 2019 / Published: 7 March 2019
(This article belongs to the Special Issue Remote Sensing and Modeling of Land Surface Water)

Abstract

:
Yongding River is the largest river flowing through Beijing, the capital city of China. In recent years, Yongding River Basin (YDRB) has witnessed increasing human impacts on water resources, posing serious challenges in hydrological and ecological health. In this study, remote sensing techniques and statistical time series approaches for hydrological studies were combined to characterize the dynamics and driving factors of reservoir water extents in YDRB during 1985–2016. First, 107 Landsat 4, 5, 7 and 8 images were used to extract surface water extents in YDRB during 1985–2016 using a combination of water indices and Otsu threshold algorithm. Significant positive correlation was found between water extents and the annual inflow for the two biggest reservoirs, the downstream Guanting and upstream Cetian reservoirs, proving their representativeness of surface water availability in this basin. Then, statistical time series approaches including trend-free pre-whitening Mann-Kendall trend test, Pettit change-point test and double mass curve method, which are frequently used in hydrological studies, were adopted to quantify the trend of reservoir water extents dynamics and the relative contributions of climate variability and human activities. Results showed that the water extents in both reservoirs exhibited significant downward trend with change point occurring in 2001 and 2005 for Guanting and Cetian, respectively. About 74%~75% of the shrinkage during the post-change period can be attributed to human activities, among which GDP, population, electricity power production, raw coal production, steel and crude iron production, value of agriculture output, and urban area were the major human drivers. Hydrological connectivity between the upstream Cetian and downstream Guanting reservoirs declined during the post-change period. Since 2012, water extents in both reservoirs recovered as a result of various governmental water management policies including the South-to-North Water Diversion Project. The methodology presented in this study can be used for analyzing the dynamics and driving mechanism of surface water resources, especially for un-gauged or poorly-gauged watersheds.

Graphical Abstract

1. Introduction

Water-related issues, including floods, droughts and water shortages have been identified as major global challenges in the 21st century [1]. In particular, the conflict between freshwater supply and demand has been widely recognized to be a source of strategic rivalry nationally, regionally and locally. While climate change might alter the hydrological cycle and the availability of water resources, human activities and socioeconomic development increasingly put pressure on freshwater supply, and may lead to water scarcity [2].
Quantifying the spatiotemporal dynamics of surface water resources and understanding their driving mechanism over the past decades will provide decision-makers with information for feasible restoration and management strategies and to further to evaluate their effects. Existing research in the hydrology community has generally utilized in-situ streamflow data obtained from hydrological gauge stations to analyze the temporal dynamics of surface water resources [3,4,5]. Compared to gauge station observations, remote sensing techniques provide direct and synoptic observations on spatiotemporal dynamics of surface water extents at regional and even global scales that can be used as indicators of water resources availability. In recent years, Landsat-series satellite datasets have been increasingly used to monitor surface water dynamics, given the advantage of their public accessibility, high resolution, and long-term availability (i.e., featuring a three-decade long data archive). Using these datasets, global open-water products have been developed, such as the Global 3 arc-second Water Body Map (G3WBM) at a five-year interval during 1990–2010 [6] at ~90 m resolution at the equator (http://hydro.iis.u-tokyo.ac.jp/~yamadai/G3WBM/), the 30 m-resolution Global Surface Water Explorer (GSWE) dataset from 1984 to 2015 [7] (https://global-surface-water.appspot.com/), and the 90m-resolution Global Inundation Extent from Multi-Satellite (GIEMS-D3) by combining Landsat-series and Synthetic Aperture Radar (SAR) datasets during the period of 1993–2007 [8,9]. Although these products provide good descriptions of surface water dynamics with regard to global coverage, they may not satisfy all the needs for regional-scale studies. Specifically, G3WBM and GIEMS-D3 have relatively short temporal coverages; GSWE do not provide annual time series of surface water extents, and only maximum water extent and water occurrence frequency were available. For time series analysis, annual or even seasonal water extents need to be extracted. Several remote-sensing-based indices are available for extracting water surfaces, including the Normalized Difference Water Index (NDWI) [10], the Modified Normalized Difference Water Index [11], the water index developed by Danaher and Collett (2006) [12], the Automated Water Extraction Index (AWEI) [13] and the water index developed by Fisher et al. (2016) [14]. All these indices allow water pixels to be identified by specifying a simple threshold, which can be adjusted for specific images. Compared to supervised classification methods, the water index method is easy to be implemented and facilitates fully automated water extraction.
In this study, we chose Yonding River Basin (YDRB) as our study region. Yongding River, the largest tributary of the Haihe River Basin (HRB), is commonly known as the mother river of the city of Beijing, China. The river basin is situated in a typical temperate semi-arid continental climate with long-term annual mean precipitation around 400 mm. The river flows through Shanxi Province, Inner Mongolia Autonomous Region, Hebei Province, and down to Beijing and Tianjin. Historically, wood harvest in the headwater region of the basin led to massive and frequent flood disasters prior to the 20th century [15]. For flood mitigation, dams and various water conservancy facilities were built in the YDRB, such as Guanting Reservoir built in 1954 and Cetian Reservoir built in 1960 [16]. In addition to flood control, these reservoirs have also played important roles in electricity power generation, and agricultural, domestic, and industrial supply. Over the past three decades, YDRB has become a region with intensive human activity. In the upper reach of YDRB, Shanxi Province produces about one-third of coal in China; Hebei Province is an important industrial base for steel and iron production [17]. Beijing, in the lower reach of YDRB, has been seeing massive developments at an unprecedented pace with rapid growth of the population and urban expansion. Increasing water abstraction has put more pressure on resources, causing serious water scarcity [1,3]. Water shortage has not only led to degradation of water quality and ecological security, hindering sustainable development, but it has also been reported that the lower reaches (the Beijing-Tianjin section) of Yongding River have been perennially dry [4]. In recent years, the state and local governments have implemented several inter-basin water transfer projects aiming to alleviate the water shortage in north China, including the South-to-North Water Diversion Project transporting water from the Yangtze River basin to north China and the Yellow River Diversion Project conveying water from the Yellow River to Shanxi Province. Intra-basin water transfer projects have also been implemented in YDRB, including water conveyance from upstream reservoirs to downstream reservoirs. Nevertheless, the degradation of aquatic ecosystems has not been fully mitigated. In 2017, the Development and Reform Commission of China announced the “The General Plan of Comprehensive Treatment and Ecological Restoration for the Yongding River Basin”, a cross-provincial program designed to restore riverine ecosystems in five to ten years.
The objectives of this study therefore are to: (1) investigate the annual dynamics and statistical trends of surface water extents in YDRB over the past 32 years (1985–2016) using Landsat-based datasets. Due to the narrow river channel width in YDRB from Landsat-series imagery, we hence mainly focus on water extents in reservoirs, as they represent the majority of the surface water storage in YDRB; (2) quantify the contribution of climate variability and human activities to variations in surface water extents in key reservoirs at the annual time scale; and (3) assess natural, socioeconomic factors and governmental policies that drive the temporal trends of reservoir water extents.

2. Study Area

The Yongding River (112°~117°45′ E, 39°~41°20′ N) is a major tributary of the Haihe River in North China (Figure 1). It is 747 km long, flowing through Shanxi Province, Inner Mongolia Autonomous Region, Hebei Province, Beijing and Tianjin, and drains into the Bohai Bay. Yanghe and Sanggan rivers are the two major tributaries of the Yongding River, originating from the Inner Mongolia Autonomous Region and Shanxi Province, respectively. Both rivers flow through the Guanting reservoir. The YDRB is located in a semi-humid and semi-arid climate transition zone with a drainage area of 47,016 km2, annual average temperatures ranging between 5.1–7.1 °C and a mean annual precipitation of ~400 mm. The climate of YDRB features a strong seasonality with 70–80% of its rainfall occurring in the summer (i.e., June to August).
As aforementioned, rapid industrial development, population growth and urban expansion have led to intensive water abstraction during the past three decades. Water management policies have been implemented in recent years to alleviate the water shortage. In this study, we focus our analyses on the Guanting reservoir and its upstream drainage area. Currently, 14 large and medium reservoirs exist in the study domain, all of which were constructed between 1950 and 1980 [16]. These reservoirs have been mostly used for flood control and water supply, while some of them are designed for electricity generation [18]. Among these reservoirs, the Guanting reservoir in Beijing has the largest area, followed by the Cetian reservoir located along the Sangan river in Shanxi Province. Sedimentation in both reservoirs is marginal, according to regional water authority bureaus. Two major prefecture-level cities upstream to the Guanting reservoir are Datong in Shanxi Province and Zhangjiakou in Hebei Province. Major land cover types in the study area are forest (13%), grassland (26%), cropland (56%) and urban (5%).

3. Materials and Methods

3.1. Remote Sensing Imagery and Additional Datasets

In this study, a total of 107 images from Landsat 4 Multispectral Scanner (MSS) at 60-m resolution (1 scene), Landsat 4-5 Thematic Mapper (TM) (87 scenes), Landsat Enhanced Thematic Mapper (ETM+) (8 scenes) and Landsat 8 Operational Land Imager (OLI) images (11 scenes) at 30-m resolution acquired between 1985 and 2016 were used to monitor the dynamics of reservoir water extent in the YDRB (Figure A1). In order to minimize the impact of seasonal variability, only images acquired from late-August to October were selected because the water extents reach the maximum in a year as a result of summer precipitation. All images are cloudless over the reservoirs. Landsat 4 MSS level-1 terrain-corrected images (1985) and Landsat TM/ETM+/OLI surface reflectance products (1986–2015) were downloaded from the United States Geological Survey’s remote sensing image database (http://earthexplorer.usgs.gov/).
Climate, hydrological and anthropogenic variables were collected and used to quantify the drivers of surface water extent variation. Annual precipitation observed at 50 precipitation stations (Figure 1) was acquired from the Hydrological Year Books of the HRB maintained by Haihe River Water Conservancy Commission and the National Meteorological Information Center of China. Annual streamflow observed at the Gudingqiao hydrological station located upstream to the Cetian Reservoir, Shixiali and Xiangshuipu stations upstream to the Guanting Reservoir were also obtained from the Hydrological Year Books. Anthropogenic variables (Table 1) were collected from the regional statistical yearbooks, including annual population, gross domestic production (GDP), electricity production, raw coal production, steel and crude iron production, and total value of agriculture output from Datong and Zhangjiakou. Urban areas were calculated from the Climate Change Initiative Land Cover (CCI-LC) dataset that provide land cover maps during long time periods (1992–2015) (http://maps.elie.ucl.ac.be/CCI/viewer/download.php). In order to justify the use of this dataset, the CCI-LC urban cover was evaluated by comparing with urban areas defined by 30-m GlobeLand30 2000 and 2010 datasets, which was reported to have high accuracies [19,20]. Comparison results showed that the difference between the two datasets are within 12%, indicating the feasibility of CCI-LC data.

3.2. Methods

Figure 2 summarizes the workflow of this study. First, annual reservoir water extents were extracted based on Landsat-series imagery using the combination of AWEI/NDWI and Otsu threshold algorithms. In order to characterize the dynamics of reservoir water extents, trend and change-point analyses were then conducted on the time series water extents; variations in the hydrological connectivity between upstream and downstream reservoirs were then analyzed. Third, the contributions of climate variabilities and human activities on the variations of reservoir water extent were quantified using Double Mass Curve (DMC) method. Finally, natural and anthropogenic driving factors of reservoir water extent dynamics were analyzed using correlation analysis. The influence of water management policies on surface water extents were also examined. The software used to implement the procedures is listed in Figure 2.

3.2.1. Image Preprocessing and Extraction of Reservoir Water Extent

For Landsat TM, ETM+ and Landsat OLI, surface reflectance data products were generated using the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) and Landsat 8 Surface Reflectance Code (LASRC), respectively [21], which are both atmospheric correction procedures officially developed by NASA Goddard Space Flight Center and the University of Maryland. These products can therefore be directly used for water body extraction. Because there were no surface reflectance products from Landsat MSS, we implemented the Fast Line-of-sight Atmospheric Analysis of Hypercubes (FLAASH) algorithm embedded in the Environment for Visualizing Images (ENVI) software to correct for atmospheric effect [22]. The data gaps in the Landsat ETM+ images caused by the Scan Line Corrector failure issue were also corrected using Landsat Gapfill algorithm embedded in ENVI software. For each image, the quality assessment band was examined to remove clouds/cloud shadows.
Previous inter-comparison studies suggest that AWEI performs well in both urban area and mountainous areas, and the optimal threshold value is fairly stable across images [13,14]. It has been successfully applied in many recent studies to extract water bodies from Landsat TM/ETM+ images [23,24,25]. Therefore, in this study, we adopted the AWEI approach to extract reservoir water extent from Landsat TM/ETM+ images within the YDRB. Because it can be challenging to separate water pixels from dark surfaces due to their similar reflectance, the AWEI approach provides two indices, namely AWEIsh and AWEInsh. AWEIsh is suitable for scenes with shadows and/or other dark surfaces, typically found in mountainous regions, while AWEInsh is suitable for scenes without shadows. As none of the reservoirs in our study area are located in mountainous regions, only AWEInsh is calculated from the Landsat TM/ETM+/OLI images in this study (Equation (1)) [13].
AWEI n s h = 4 × ( ρ G r e e n ρ S W I R 1 ) ( 0.25 × ρ N I R + 2.75 × ρ S W I R 2 )
where ρ G r e e n , ρ N I R , ρ S W I R 1 , and ρ S W I R 2 are the surface reflectance of green, near infrared (NIR), shortwave infrared 1 (SWIR1) and shortwave infrared 2 (SWIR2) bands, respectively.
For the Landsat MSS images, NDWI was used instead because they only consist of four spectral bands, namely green, red, and two NIR bands (700–800 um, 800–1100 um). The NDWI is defined as [10]:
NDWI = ( ρ G r e e n ρ N I R 2 ) / ( ρ G r e e n + ρ N I R 2 )
After the water index map was created for each image, a threshold value is applied to convert it into a binary map, with “1” denoting water and “0” denoting non-water. In this study, we adopted the Otsu threshold algorithm [26], which is a nonparametric and unsupervised approach for automatic threshold selection. This method determines optimal threshold as the one that maximizes the inter-class variance based on the AWEI or NDWI histogram [27,28]. In our study, the optimal thresholds for the AWEI and NDWI maps ranged from 0.09 to 0.13. Finally, reservoir water extent was extracted over the subset of the binary maps around each reservoir.
The water extent results were then evaluated using high spatial resolution images from Google Earth and Sentinel-2 MSI images with 10-m resolution. Google Earth images were carefully selected so that the acquisition dates were within 10 days from the Landsat images. In the study area, two scenes of 2-m QuickBird images covering Guanting Reservoir (2005) and Yiliuhe Reservoir (2015) were selected as reference image. Sentinel-2 images covering Guanting and Cetian reservoirs (2016) with acquisition dates close to Landsat scenes were also used as reference images. Two thousand sampling points were randomly generated within the 5-km buffer area around each reservoir, visually interpreted as water or non-water, and used as reference datasets.

3.2.2. Trend and Change-point Analyses of Reservoir Water Extents

Once the annual reservoir water extents were extracted, the trend-free pre-whitening (TFPW) Mann-Kendall (M-K) test and Pettit mutation test were adopted to determine the trend and change point of water extents from 1985 to 2016. The Mann-Kendall (M-K) trend test is a nonparametric test method that has been extensively used for time-series hydrological datasets such as runoff [29]. The samples for the M-K test do not need to obey a prescribed distribution and is not sensitive to outliers. In order to eliminate autocorrelation and improve accuracy of trend detection, the TFPW strategy [30] was applied prior to the M-K test. Meanwhile, the Theil-Sen method was used to determine the magnitude of trends in reservoir water surface area and precipitation [31,32].
Change point detection is an important tool for determining if and when a change in a time-series data set has occurred, and has proven to be useful in analyzing the effects of natural factors and human activities on hydrological time series. The Pettitt test [33] is a rank-based nonparametric statistical test method that poses no assumption on the distribution of the data. It adapts the Mann-Whitney test that allows the identification of the time at which the change occurs. Given a time series represented by ( x 1 ,   x 2 ,   ,   x t ,   x t + 1 ,   ,   x T ) , which represents the time-series reservoir water extents in this study, the statistical index U(t,T) is defined as follows [33]:
U t , T = i = 1 t j = t + 1 T sgn ( x j x i ) , 1 t < T
where
sgn ( θ ) = { 1 0 1 θ > 0 θ = 0 θ < 0 ,
When the time series follows a continuous distribution, the test statistics U(t,T) can also be calculated by the following recursive iteration:
U t , T = U t 1 , T + V t , T , 2 t T ,
V t , T = j = 1 T sgn ( x j x t ) , 1 t T ,
The statistics of the most probable change point τ and the significance probability of Kτ are defined as follows:
K τ = | U τ , T | = max | U t , T | ,
p = 2 exp ( 6 K τ 2 T 2 + T 3 )
When p is less than a certain significance level α, we consider that xτ is a significant change point at level α.

3.2.3. Quantifying Relative Contributions of Climate Variability and Human Activities on Reservoir Water Extents

Climate variability and human interventions are commonly considered two major contributors to variations of water resources in a river basin. In this study, the double mass curve (DMC) method, which was frequently used in hydrological studies, was adopted to quantify their relative contributions to reservoir water extents dynamics. The DMC is a simple visual method composed of cumulative values of two datasets plotted against one another during a given time period [34]. It has been used to evaluate the contributions of potential drivers to changes in streamflow or sediment load [35,36]. Generally, the plot of cumulative streamflow against cumulative precipitation is a straight line if the effect of human activities on streamflow is constant or negligible. Significant changes in DMC slopes indicates departure from the original relationship, suggesting changes in contributions of human activities.
In this study, we plotted the cumulative reservoir water extent (A) against the cumulative upstream precipitation (P). In order to evaluate the change point detected using methods described in Section 3.2, the nonparametric Mann-Whitney U test was used to compare DMC slopes before and after each change point on the curve statistically [37]. Once the significance of the change point is confirmed, regression equations were established based on the data points in the pre-change (reference) and post-change (disturbance) periods. The baseline regression equation in the reference period was then employed to predict the cumulative reservoir water extents for the disturbance period. The difference between the observed and predicted cumulative water extents in the disturbance period was considered as the deviation caused by intensified human activities (∆Ah). Thus, the deviation caused by climate variability (∆Ac) can be calculated as [34]:
∆Ac = ∆A − ∆Ah
where ∆A, ∆Ac, and ∆Ah are the deviations of annual water extent between disturbance and reference periods, annual water extent deviation caused by climate change, and annual water extent deviation caused by human activities, respectively. The relative contributions of climate variability and human activities to water extents can be determined as:
R c = | Δ A c | | Δ A c | +   | Δ A h |   ×   100 %
R h =   | Δ A h | | Δ A c | +   | Δ A h |   ×   100 %
where R c and R h are the relative contributions of climate variability and human activities to water extents, respectively.

3.2.4. Correlation Analysis of Driving Factors and Reservoir Water Extents

For both pre-change and post-change periods, Pearson correlation coefficients were calculated to examine the effect of annual precipitation and socio-economic factors on reservoir water extent dynamics. Here precipitation was treated as a natural driving factor. Human factors include GDP, population, electricity production, raw coal, steel and crude iron productions, total value of agriculture output and urban area, as they are all associated with water usage. For the Guanting reservoir, the values of human factors were calculated as the sum of Datong and Zhangjiakou, as they are the two major upstream cities. For the Cetian Reservoir, the values for Datong was used. Steel and crude iron production was not included for the Cetian reservoir because of its low production (average annual production < 40 × 104 ton). High correlation coefficients indicate stronger possible impacts on variations in reservoir water extent. The differences between pre- and post- change coefficients suggest that the main driving factors of the reservoir area might have changed. The correlation coefficient between upstream and downstream reservoir water extents can be used to assess the hydrological connectivity in the river basin. Higher correlation coefficients indicate stronger connections, and can thus be interpreted as stronger influence of the upstream on downstream.

4. Results

4.1. Spatiotemporal Dynamics of Surface Water Extents in Reservoirs

The water extents of the 14 reservoirs in YDRB were extracted based on the 107 Landsat-series images from 1985 to 2016 using AWEI or NDWI approaches. Using Google Earth images as references, the overall accuracies vary from 98.8% to 99.7%, and the kappa coefficients vary from 0.95 to 0.98 (Table 2).
Figure 3 shows the frequency of surface water occurrence in the study area, calculated as the number of years that a pixel has been detected as open water, divided by the total number of years (i.e., 32 years in this study). The statistics of annual surface water extents in the 14 reservoirs are shown in Figure 4. The Guanting and Cetian reservoirs are the two biggest reservoirs in the YDRB. The average water extents in the Guanting and Cetian reservoirs are 68.6 km2 and 16.8 km2, respectively, which account for over 75% of the total water extent in the basin. Both reservoirs show great variations in water extents, ranging from 31.9 km2 to 121.3 km2 in Guanting and 6.5 km2 to 33.8 km2 in Cetian. Other reservoirs are very small, with an average area of less than 5 km2, and several of them less than 1 km2. Nonetheless, almost all reservoirs showed similar temporal dynamics during the whole study period, with increase in their water extent at the beginning and decrease afterwards. Most reservoirs reached their maximum areas in 1995 or 1996, including Guanting, Cetian and Youyi.
Figure 5 showed inter-annual variability in water extents in both Guanting and Cetian reservoirs (see Figure A2 and Figure A3 for the maps of annual water extent). The water extent dynamics of Guanting can be divided into four stages: slow rising period, severe atrophy period, stable period and recovery period. In the slow rising period of 1985–1995, its area increased from 66.7 km2 to 121.3 km2, with an expansion of 81.9%. In the severe atrophy period of 1996–2007, its area decreased sharply from 113.2 km2 to 31.9 km2, shrinking by 71.8%. From 2008–2012, the area was relatively stable, with a mean area of 36.2 km2. During 2013–2016, the reservoir area rebounded slightly, increasing from 41.6 km2 to 77.1 km2, with an average area of 62.9 km2. The area of Cetian is much smaller than the Guanting reservoir, but varied significantly over the study period figuring three states (Figure 4): a slow rising period (1985–1995), a shrinking period (1996–2011) and a recovery period (2012–2016). During the slow rising period of 1985–1995, the water extent increased from 11.8 km2 to 33.8 km2, with an average annual increase rate of 2.2 km2/year. In the shrinking period from 1996 to 2011, the reservoir area decreased from 33.3 km2 to 9.3 km2, reducing by 72.1%. During the 2012–2016 period, its area recovered from 10.2 km2 to 15.6 km2. Given that the total surface area of all other reservoirs is too small, Guanting and Cetian were selected as reservoirs for analyzing variability of surface water dynamics in this basin.

4.2. Trend and Change-points of Dynamic Water Extents

The TFPW-MK trend test showed that the water extents in Guanting and Cetian had a significant downward trend from 1985 to 2016 (Table 3), and the annual average declination rate were 1.56 km2/year for Guanting and 0.116 km2/year for Cetian, respectively. The Pettitt test showed that the significant change-point of time-series water extent occurred in 2001 (p < 0. 01) for Guanting and 2005 (p < 0.05) for Cetian (Table 3). For the Guanting Reservoir, the average water extent was 88.5 km2 in the pre-change period of 1985 to 2001 and 46.1 km2 in the post-change period of 2002 to 2016, decreasing by 47.9%. The average water extent of Cetian was 18.8 km2 in the pre-change period of 1985 to 2005 and reduced to 12.9 km2 during the post-change period of 2006 to 2016, decreasing by 31.4%. Interestingly, from Figure 5 we can visually detect an apparent upward trend in Guanting since 2012. Nonetheless, the Pettit test on the water extents during the post-change period (2002–2016) indicated 2012 as the most probable change-point with p-value of 0.18 due to the small sample size.
As illustrated in Figure 6a,b, both Guanting and Cetian reservoirs demonstrated similar temporal trends in their annual inflow observed from their upstream hydrological stations. Note that the Guanting inflow was calculated as the sum of streamflow observed at the Shixiali and Xiangshuipu hydrological stations (Figure 1). The TFPW-MK trend test showed that the streamflow also demonstrated significant downward trend with change points occurring in 1999 for the inflow of Guanting and 2005 for the inflow of Cetian. For the Cetian Reservoir, the correlation analysis showed a significantly positive relationship between annual streamflow and reservoir water extent. The Pearson correlation coefficients are 0.81 (p < 0.01) for the whole study period from 1985 to 2016, 0.77 (p < 0.01) for the pre-change period (1985–2005), and 0.72 (p < 0.05) for the post-change period (2006–2016). For the Guanting reservoir, the positive correlations are statistically significant for the whole study period (r = 0.83, p < 0.01) and pre-change period (r = 0.78, p < 0.01) (1985–2001), while the correlation coefficient for the post-change period (2002–2016) is 0.49 and p-value is 0.07. Correlations between water extents in Guanting and Cetian extents are 0.55 (p < 0.01) for the 32 years, and 0.72 (p < 0.01) during the Guanting pre-change period (1985–2001) and 0.25 (p > 0.05) during the post-change period (2002–2016).

4.3. Relative Contributions of Climate Variability and Human Activity to Reservoir Water Extents

As shown in Figure 7, the linear regression lines between the cumulative water extents and the cumulative annual precipitation have different slopes before and after the change point. By performing Mann-Whitney U tests [36], the slopes are confirmed to be statistically different, and the change points detected by the Pettit test are reasonable (p < 0.01 for both reservoirs).
Water extents in Guanting and Cetian reservoirs increased 5.8 ± 16.4 km2 and 0.8 ± 3.1 km2, respectively, in response to climate variability. The relative contributions to the changes in water extents being 25.9 ± 22.4% and 24.9 ± 17.3%. Human intervention led to decreases in annual water extents in Guanting and Cetian reservoirs of 48.1 ± 20.0 km2 and 6.7 ± 3.0 km2, respectively, contributing 74.1 ± 22.4% and 75.1 ± 17.3% to the changes in reservoir water extents (Table 4). Evidently, contributions of human activities were much higher than those from climate variability in both reservoirs after the change point. For Guanting, the disturbance period can be further divided into two sub-periods, i.e., 2002–2012 and 2013–2016, based on the distinct expansion rates of water extents (Figure 5). The negative contribution of human activities during 2002–2012 was higher than that during 2013–2016. Our results suggest that human activities and climate variability played opposite roles in variations of water extents in both reservoirs, while the negative effect of human activities seems to decrease after 2012.

4.4. Correlations between Driving Factors and Reservoir Water Extents

The graphical trends of human factors and water extents (Figure 8 and Figure 9) showed that GDP, population, urban area, electricity energy production and value of agriculture output in the upstream regions showed a steady positive trend since 1985, with a rapid increase in the slope after 2000 for both Guanting and Cetian reservoirs. The raw coal production peaked in 1995, decreased from 1995 to 2000, followed by a sharp increase since 2001, matching the expansion and shrinkage of the water extents in both reservoirs (Figure 8e and Figure 9e). For the Guanting reservoir, the upstream regions exhibited sharp increase in steel and crude iron production since 2000, which becomes stable after 2012.
Table 5 lists the Pearson’s correlation coefficients between water extents and the drivers for both reservoirs. Precipitation had significantly positive relationships (p < 0.01) with water extents during the pre-change periods for both reservoirs, while after the break points the positive correlations become insignificant. It is interesting that none of the factors had significant correlations with water extents during the post-change periods. As Figure 6 indicates evident water extent expansions after 2012, we further analyzed the correlations between the break points and 2012. It can be seen that GDP, population, electricity energy production, raw coal, steel and crude iron productions, value of agriculture output and urban area are all negatively correlated to water extents in Guanting reservoir during this period. The scatterplots of the water extents against the driving factors (Figure A4 and Figure A5) showed negative slope of the linear regression line. In Cetian, the negative correlation coefficients are statistically significant (p < 0.05), and the regression line has R2 greater than 0.55, suggesting negative effects of human activities on water extents.

5. Discussions

5.1. Characteristics of Reservoir Water Extents Dynamics

5.1.1. Trend of Reservoir Water Extents and Their Representativeness of Water Resources in YDRB

The Guanting and Cetian reservoirs are the two biggest surface water bodies in the YDRB. From 1985 to 2016, the water extents in both reservoirs experience similar temporal patterns. They both expanded from 1985 to 1995, reached the maximum extent in 1995, and then shrank until 2011. From 2012, the water extents of both reservoirs started to rebound. The combination of the TFPW M-K and Pettit tests showed that both reservoirs experienced downward trend over the 32 years, and the abrupt change points occurred in 2000s. Our results for Guanting reservoir are consistent with a recent study, which reported that the water level in its hydro-fluctuation zone had dropped 8.19 m from 1996 to 2007 [38].
For both reservoirs, the time-series of water extents showed consistent temporal patterns with the streamflow observed from upstream hydrological gauging stations. The positive correlation between streamflow and water extents were statistically significant throughout the study period (p < 0.01). Trend analyses on annual streamflow also showed that abrupt changes occurred closely with the reservoir water extents. This consistency suggests that surface water extents can be used to represent variations in water resources in YDRB.
Previous studies have suggested that surface water extents can be used as a proxy for streamflow in hydrological modelling. Revilla-Romero et al. (2016) [39] assimilated satellite-derived surface water extents into a rainfall-runoff model to update simulated streamflow, and demonstrated that such an integration improved streamflow simulations. Similarly, Revilla-Romero et al. (2015) [40] employed remotely-sensed surface water extents from the Global Flood Detection System (GFDS) to calibrate a hydrological model and demonstrated the potential of the satellite-derived water extent for improving hydrologic model calibration in ungauged or poorly gauged watersheds. Huang et al. (2014) [41] reported high correlation coefficients between satellite-derived wetland water extents and the stream flow, and suggested that inundation maps can be used for improving model studies of drought and stream runoff. The findings of these studies, as well as our results, suggest that the satellite-derived water extents can be considered as feasible indicators for time-series water resource dynamics in YDRB.

5.1.2. Hydrological Connectivity Discovered from Reservoir Water Extents Dynamics

Correlation analyses between upstream and downstream streamflow observed at gauge stations have been frequently used to examine hydrological connectivity in a watershed [4,42]. The correlation coefficient of reservoir water extents between Cetian and Guanting reservoirs was statistically significant over the past three decades (p < 0.01), indicating strong connections between upstream and downstream reservoirs in YDRB. That is, variability in upstream reservoirs has a strong influence on downstream reservoirs.
However, we found that the correlation in the pre-change period (r = 0.715, p < 0.01) was stronger than that in the post-change period (r = 0.254, p > 0.05). Although showing similar temporal trends, the downstream Guanting reservoir experienced more rapid shrinkage than the upstream Cetian reservoir. In addition, the abrupt change in Guanting water extent occurred four years earlier than that of Cetian. This indicates that the faster drying of Guanting can be attributed to more intensified upstream human activities.
Jiang et al. (2014) [4] studied the hydrological connections in the Beijing section of Yongding River, which is located in the lower reach starting from Guanting. They reported that the hydrological connections between Guanting and the downstream Sanjiadian had declined since 2000 due to the prolonged zero flows below Sanjiadian. In our study, we found that the hydrological connections above the Guanting Reservoir also declined during the post-change period, which is likely to contribute to the zero flows in the lower reach of YDRB.

5.2. Relative Contributions of Climate Variability and Human Activities to Reservoir Water Extents

Our study showed that the relative contributions of human activities to water extents in both reservoirs were around 74%~75% compared to those in the pre-change period for both reservoirs. Previous studies have investigated the contributions of climate variability and human activities to runoff variations in HRB based on the streamflow observations at gauge stations. Liu et al. (2017) [5] and Xu et al. (2014) [43] reported that 66.2%~73.1% of runoff reduction in HRB was attributed to human activities based on the Budyko hypothesis model. Wang et al. (2013) [44] evaluated four catchments in HRB, and reported that human activities contributed to 63.1% of runoff reduction. Xia et al. (2014) [45] studied the runoff variations in YDRB from 1960 to 2010, and reported that the contributions from human activities to runoff declination accounted for 71.0%~92.7%. Although our study focused on the reservoir water extents instead of runoff, we have shown that the reservoir water extents correlated significantly with runoff, and the relative contributions of climate variability and human activities to water extents were similar. Most previous studies have focused on periods before 2010. Our analysis updates the knowledge on changes in water resources and their driving factors in YDRB, highlighting the prolonged intensification of human activities in the recent decades.

5.3. Drivers for Reservoir Water Extent Dynamics

Precipitation had significant positive influences on the reservoir water extents variations during the pre-change period. During the post-change period, the positive influence become weaker, confirming that precipitation had a minor influence on the declination of reservoir water extents. From the change point to the year 2012, almost all socioeconomic factors were negatively associated with water extents. Since 2000, YDRB has witnessed increasingly rapid socioeconomic development. For the upstream regions of Guanting Reservoir, the growth rate of GDP was 2.4 billion Yuan RMB per year from 1985 to 2001, which increased to 14.3 billion CNY per year from 2002 to 2016. During the post-change period, the growth rate of GDP in the upstream region of Cetian reservoir increased more than twice than that of the pre-change period. Population and urban area also increased more rapidly after the change point. The fast economic development, population growth and urban expansion led to increase of domestic and industrial water demands in the basin. For both reservoirs, human-induced water extents changes were also negatively correlated with crop production. Increasing crop production increases water usage for agriculture irrigation and hence decreases reservoir water storage. In YDRB, the primary industrial water consumption comes from coal power generation, metallurgy, and food production [46]. Steel and crude iron are the mainstay industry in Zhangjiakou city in Hebei Province, which is located in the upper reaches of the Yonding River. Datong city in Shanxi Province, which is close to the Cetian reservoir, is home to one of the largest coal deposits in China. Coal mining has been the mainstay economic industry in Shanxi Province, and is also a major surface water consumer [47]. Therefore, the sharp increase in productivity of mainstay industries in the upstream regions, i.e., raw coal for Cetian and steel and crude iron for Guanting led to water extents declination in the reservoirs.

5.4. Effects of Policies on Reservoir Water Extents

Table 5 showed that human activities correlated with reservoir water extents negatively during the period from the change-point year to 2012, while their effects were not significant in the entire post-change period. The water extents in both reservoirs have recovered since 2012, despite continuous rise of socioeconomic indicators. While this can be partially attributed to increasing precipitation from 2012, the implementation of water-related policies also contributed to the recovery of water extents.
In 2011, the Shanxi Wanjiazhai Yellow River Diversion Project in the northwest region of Shanxi Province was completed [48]. It was one of the strategic measures by the ShanXi Province to transfer water from the Yellow River and alleviate water shortages in the three industrial cities of Datong, Taiyuan, and Pingsuo. The project brought 160 million cubic meters of water to Datong each year. The supply of the Yellow River water reduced water deficit in Datong, therefore allowing more water to drain downstream, leading to the rebound of water extent in the Cetain reservoir. Meanwhile, the middle route of the South-to-North Water Diversion Project, which originates from Hanjiang river, began its operation in December 2014. By the end of 2016, the volume of water coming to Beijing from the South Water has reached 2.3 billion cubic meters, of which 1.7 billion cubic meters was supplied to tap water plants, and the remainder directed to reservoirs, emergency water sources, rivers and lakes (http://www.bjwater.gov.cn). All these policies help alleviate the pressure of water demand from the reservoirs; hence, the surface areas of Guanting and Cetian were able to recover.
These policies affected not only the inter-annual but also intra-annual variability of reservoir water extents. After 1996, the storage of the Guanting reservoir continued to decrease. In order to rationally allocate limited water resources in the YDRB, in 2001, the Government of Beijing and the Ministry of Water Resources of the People’s Republic of China authorized the "Capital Water Resources Sustainable Utilization Plan in the Early 21st Century" [49]. To achieve this goal, from 2003, the Cetian reservoir began to transfer water to the Guanting reservoir (Table 6). Except for 2009, 2011 and 2012, water was delivered during early-October to mid-November each year. The total amount of water delivered reached 393 million m3 since 2003. The water delivery caused distinct variations in reservoir water extents during the two months. Figure 10 demonstrates the water extent before and after water delivery from 2002 to 2008 observed from Landsat-series imagery following the same procedure in Section 3.2.1. Starting from 2003, the area of Guanting increased obviously after water delivery, and the water extent expansion was proportional to the amount of water inflows. The area of Cetian decreased after water delivery, and the water extent shrinkage was also proportional to the water output. As a reference, the water extents in both reservoirs did not change during September to November in 2002, when centralized water delivery was not initiated.

6. Conclusions

In this study, the dynamics of surface water extents in reservoirs of YDRB was characterized based on time-series of Landsat 4 MSS, Landsat 4-5 TM, Landsat 7 ETM+ and Landsat 8 OLI images during 1985-2016 using a combination of AWEI, NDWI and Otsu threshold algorithms. All reservoirs demonstrated considerable variations in water extents over the past 32 years. The downstream Guanting and upstream Cetian reservoirs were identified as the biggest open water bodies, with a combined water area of over 75% of the total reservoir areas in the basin.
We proved that the reservoir water extents in YDRB are representative of variations in surface water availability, as the annual inflow to each reservoir and the reservoir water extents showed consistent temporal trend and had significant positive correlation. The time series approaches, namely TFPW M-K test and Pettit test, which were previously used for streamflow analysis, were then applied to examine the temporal trend and change-point of water extents dynamics. The DMCs, which were previously used to assess the trend of river discharges, were further adopted to analyze the relative contribution of human activities and climate change on water extent variability. Results showed that both reservoirs experienced statistically significant downward trend in water extent, and the most probable change points occurred in 2001 and 2005 for the Guanting reservoir and Cetian reservoir, respectively. The relative contribution of human activities to shrinkage of water extents was 74%~75% for both reservoirs during the post-change period.
Precipitation, socio-economic factors and policies are all drivers for variations in reservoir water extents, while they played different roles in different periods. Precipitation is the major driving factor during the pre-change periods. From the change-point year (2006) to 2012, human drivers including GDP, population, electricity power production, raw coal production, value of agriculture output, and urban area played more important roles and negatively influenced water extents in Cetian. For the downstream Guanting reservoir, steel and crude iron production was also an important factor.
From the remotely-sensed water extents, we identified a reduced hydrological connectivity between the upstream Cetian reservoir and downstream Guanting reservoir during the post-change period. The effects of inter-basin and intra-basin water transfer projects implemented by regional governments, including the South-to-North Water Diversion Project, were also revealed by the recovery of water extents since 2012 and the intra-annual change of the water extents.
Existing studies in the remote-sensing community mainly focused on methodological development for water extent extraction, while studies in hydrology communities usually relied on streamflow measurements at hydrological stations to analyze the dynamics and driving factors of water resources. The framework presented in this study linked the remote-sensing techniques and hydrological time series analysis approaches to assess the dynamics and driving factors of water resource variations. Compared to streamflow observations in hydrological stations, remote-sensing technique features high efficiency and low cost. This is potentially useful in un-gauged or poorly-gauged watersheds. The hydrological time series approaches help quantify the contributions of natural and anthropogenic factors to variations of water availability in the basin, which could provide guidance for better decision-making for water resources management in the basin.

Author Contributions

M.W. performed the experiment, analyzed the data, and wrote the manuscript draft. Y.K. provided crucial guidance and support during the research, and revised the manuscript. L.D. provided crucial guidance and some research data during the research. M.H. provided crucial guidance and revised the manuscript. J.Z. provided crucial research data. Y.Z., X.L. and H.G. made important suggestions for data processing.

Funding

This research was supported by the National Key R&D Program of China [2017YFC0505903], the Beijing Natural Science Foundation under Grant [5172002], by the National Natural Science Foundation of China under Grant [41401493], and by Capacity Building for Sci-Tech Innovation—Fundamental Scientific Research Funds.

Acknowledgments

This research was supported by the National Key R&D Program of China [2017YFC0505903], the Beijing Natural Science Foundation under Grant [5172002], by the National Natural Science Foundation of China under Grant [41401493], and by Capacity Building for Sci-Tech Innovation—Fundamental Scientific Research Funds. M.H. was supported by the U.S. Department of Energy, Office of Science, as part of research in Multi-Sector Dynamics, Earth and Environmental System Modeling Program. We also thank the reviewers for their valuable suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Acquisition dates of Landsat-series images.
Figure A1. Acquisition dates of Landsat-series images.
Remotesensing 11 00560 g0a1
Figure A2. Water extent in Guanting reservoir superimposed on Landsat-series images.
Figure A2. Water extent in Guanting reservoir superimposed on Landsat-series images.
Remotesensing 11 00560 g0a2
Figure A3. Water extent in Cetian reservoir superimposed on Landsat-series images.
Figure A3. Water extent in Cetian reservoir superimposed on Landsat-series images.
Remotesensing 11 00560 g0a3
Figure A4. Scatterplots of Guanting reservoir water extents against (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) steel production, (g) crude iron production, (h) value of agriculture output and (i) urban area during 1985–2001 and 2002–2012.
Figure A4. Scatterplots of Guanting reservoir water extents against (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) steel production, (g) crude iron production, (h) value of agriculture output and (i) urban area during 1985–2001 and 2002–2012.
Remotesensing 11 00560 g0a4aRemotesensing 11 00560 g0a4b
Figure A5. Scatterplots of Cetian reservoir water extents against (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) value of agriculture output and (g) urban area during 1985–2005 and 2006–2012.
Figure A5. Scatterplots of Cetian reservoir water extents against (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) value of agriculture output and (g) urban area during 1985–2005 and 2006–2012.
Remotesensing 11 00560 g0a5aRemotesensing 11 00560 g0a5b

References

  1. Palmer, M.A.; Reidy Liermann, C.A.; Nilsson, C.; Flörke, M.; Alcamo, J.; Lake, P.S.; Bond, N. Climate change and the world’s river basins: Anticipating management options. Front. Ecol. Environ. 2008, 6, 81–89. [Google Scholar] [CrossRef]
  2. Veldkamp, T.I.E.; Wada, Y.; Aerts, J.C.J.H.; Doll, P.; Gosling, S.N.; Liu, J.; Masaki, Y.; Oki, T.; Ostberg, S.; Pokhrel, Y.; et al. Water scarcity hotspots travel downstream due to human interventions in the 20th and 21st century. Nat. Commun. 2017, 8, 15697. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, L.; Wang, Z.J.; Koike, T.; Yin, H.; Yang, D.W.; He, S. The assessment of surface water resources for the semi-arid Yongding River Basin from 1956 to 2000 and the impact of land use change. Hydrol. Process. 2010, 24, 1123–1132. [Google Scholar] [CrossRef]
  4. Jiang, B.; Wong, C.P.; Lu, F.; Ouyang, Z.Y.; Wang, Y.J. Drivers of drying on the Yongding River in Beijing. J. Hydrol. 2014, 519, 69–79. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, J.Y.; Zhang, Q.; Singh, V.P.; Shi, P. Contribution of multiple climatic variables and human activities to streamflow changes across China. J. Hydrol. 2017, 545, 145–162. [Google Scholar] [CrossRef]
  6. Yamazaki, D.; Trigg, M.A.; Ikeshima, D. Development of a global similar to 90 m water body map using multi-temporal Landsat images. Remote Sens. Environ. 2015, 171, 337–351. [Google Scholar] [CrossRef]
  7. Pekel, J.F.; Cottam, A.; Gorelick, N.; Belward, A.S. High-resolution mapping of global surface water and its long-term changes. Nature 2016, 540, 418–422. [Google Scholar] [CrossRef] [PubMed]
  8. Aires, F.; Papa, F.; Prigent, C.; Crétaux, J.-F.; Berge-Nguyen, M. Characterization and Space–Time Downscaling of the Inundation Extent over the Inner Niger Delta Using GIEMS and MODIS Data. J. Hydrometeorol. 2014, 15, 171–192. [Google Scholar] [CrossRef] [Green Version]
  9. Prigent, C.; Papa, F.; Aires, F.; Rossow, W.B.; Matthews, E. Global inundation dynamics inferred from multiple satellite observations, 1993–2000. J. Geophys. Res. Atmos. 2007, 112. [Google Scholar] [CrossRef] [Green Version]
  10. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  11. Xu, H.Q. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. Int. J. Remote Sens. 2006, 27, 3025–3033. [Google Scholar] [CrossRef]
  12. Danaher, T.; Collett, L. Development, optimisation and multi-temporal application of a simple Landsat based water index. In Proceedings of the 13th Australasian Remote Sensing and Photogrammetry Conference, Canberra, Australia, 29 September–3 October 2006. [Google Scholar]
  13. Feyisa, G.L.; Meilby, H.; Fensholt, R.; Proud, S.R. Automated Water Extraction Index: A new technique for surface water mapping using Landsat imagery. Remote Sens. Environ. 2014, 140, 23–35. [Google Scholar] [CrossRef]
  14. Fisher, A.; Flood, N.; Danaher, T. Comparing Landsat water index methods for automated water classification in eastern Australia. Remote Sens. Environ. 2016, 175, 167–182. [Google Scholar] [CrossRef]
  15. Li, H.Y. Impact analysis of historical forest change on Yongding River. China Water Resour. 2005, 18, 23. (In Chinese) [Google Scholar]
  16. Haihe River Water Conservancy Commission, Ministry of Water Resources of the People′s Republic of China. Maps of Locations of Hydraulic Structures in Haihe River Basin. 2011. (In Chinese) [Google Scholar]
  17. Wang, D.Y. Retrospect and Prospect of Hebei Iron and Steel Industry Development During the 40 Years of Reform and Opening up. Econ. Manag. 2018, 32, 1–7. (In Chinese) [Google Scholar]
  18. Fan, S. The Position and Function of Cetian Reservoir in Comprehensive Utilization of Water Resources in Datong City. Shanxi Water Resour. 2009, 25, 11–12. (In Chinese) [Google Scholar]
  19. Chen, J.; Chen, J.; Liao, A.; Cao, X.; Chen, L.; Chen, X.; He, C.; Han, G.; Peng, S.; Lu, M.J.; et al. Global land cover mapping at 30 m resolution: A POK-based operational approach. ISPRS J. Photogramm. Remote Sens. 2015, 103, 7–27. [Google Scholar] [CrossRef]
  20. Yang, Y.; Xiao, P.; Feng, X.; Li, H.; Sensing, R. Accuracy assessment of seven global land cover datasets over China. ISPRS J. Photogramm. Remote Sens. 2017, 125, 156–173. [Google Scholar] [CrossRef]
  21. Masek, J.G.; Vermote, E.F.; Saleous, N.E.; Wolfe, R.; Hall, F.G.; Huemmrich, K.F.; Gao, F.; Kutler, J.; Lim, T.-K.; Letters, R.S. A Landsat surface reflectance dataset for North America, 1990–2000. IEEE Geosci. Remote Sens. Lett. 2006, 3, 68–72. [Google Scholar] [CrossRef]
  22. Ke, Y.H.; Im, J.; Lee, J.; Gong, H.L.; Ryu, Y. Characteristics of Landsat 8 OLI-derived NDVI by comparison with multiple satellite sensors and in-situ observations. Remote Sens. Environ. 2015, 164, 298–313. [Google Scholar] [CrossRef]
  23. Huang, C.; Chen, Y.; Zhang, S.Q.; Wu, J.P. Detecting, Extracting, and Monitoring Surface Water from Space Using Optical Sensors: A Review. Rev. Geophys. 2018, 56, 333–360. [Google Scholar] [CrossRef]
  24. Tulbure, M.G.; Broich, M.; Stehman, S.V.; Kommareddy, A. Surface water extent dynamics from three decades of seasonally continuous Landsat time series at subcontinental scale in a semi-arid region. Remote Sens. Environ. 2016, 178, 142–157. [Google Scholar] [CrossRef]
  25. Xie, H.; Luo, X.; Xu, X.; Pan, H.Y.; Tong, X.H. Evaluation of Landsat 8 OLI imagery for unsupervised inland water extraction. Int. J. Remote Sens. 2016, 37, 1826–1844. [Google Scholar] [CrossRef]
  26. Otsu, N.; Cybernetics, A. Threshold Selection Method from Gray-Level Histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef]
  27. Guo, Q.D.; Pu, R.L.; Li, J.L.; Cheng, J. A weighted normalized difference water index for water extraction using Landsat imagery. Int. J. Remote Sens. 2017, 38, 5430–5445. [Google Scholar] [CrossRef]
  28. Xie, H.; Luo, X.; Xu, X.; Pan, H.Y.; Tong, X.H. Automated Subpixel Surface Water Mapping from Heterogeneous Urban Environments Using Landsat 8 OLI Imagery. Remote Sens. 2016, 8, 584. [Google Scholar] [CrossRef]
  29. Choi, W.; Nauth, K.; Choi, J.; Becker, S.J. Urbanization and rainfall–runoff relationships in the Milwaukee River basin. Prof. Geogr. 2016, 68, 14–25. [Google Scholar] [CrossRef]
  30. Von Storch, H. Misuses of Statistical Analysis in Climate Research. In Analysis of Climate Variability; Springer: New York, NY, USA, 1995; pp. 11–26. [Google Scholar]
  31. Theil, H. A rank-invariant method of linear and polynomial regression analysis. In Henri Theil’s Contributions to Economics and Econometrics; Springer: Dordrecht, The Netherlands, 1992; pp. 345–381. [Google Scholar]
  32. Sen, P.K. Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  33. Pettitt, A.N. A non-parametric approach to the change-point problem. Appl. Stat. 1979, 28, 126–135. [Google Scholar] [CrossRef]
  34. Searcy, J.K.; Hardison, C.H. Double-Mass Curves, Manual of Hydrology: Part 1. General Surface-Water Techniques. 1960. Available online: https://pubs.usgs.gov/wsp/1541b/report.pdf (accessed on 31 January 2019).
  35. Costigan, K.H.; Daniels, M.D.; Dodds, W.K. Fundamental spatial and temporal disconnections in the hydrology of an intermittent prairie headwater network. J. Hydrol. 2015, 522, 305–316. [Google Scholar] [CrossRef]
  36. Wang, X.; Hao, G.; Yang, Z.; Liang, P.; Cai, Y.; Li, C.; Sun, L.; Zhu, J. Variation analysis of streamflow and ecological flow for the twin rivers of the Miyun Reservoir Basin in northern China from 1963 to 2011. Sci. Total Environ. 2015, 536, 739–749. [Google Scholar] [CrossRef] [PubMed]
  37. Gao, P.; Li, P.; Zhao, B.; Xu, R.; Zhao, G.; Sun, W.; Mu, X. Use of double mass curves in hydrologic benefit evaluations. Hydrol. Process. 2017, 31, 4639–4646. [Google Scholar] [CrossRef]
  38. Gong, Z.N.; Chen, Y.Z.; Zhao, W.H. Analysis of fluctuation characteristics and subzone in hydro-fluctuation belt in the Guanting Reservoir during 1979–2013. Acta Ecol. Sin. 2017, 37, 5294–5304. [Google Scholar]
  39. Revilla-Romero, B.; Wanders, N.; Burek, P.; Salamon, P.; de Roo, A. Integrating remotely sensed surface water extent into continental scale hydrology. J Hydrol. 2016, 543, 659–670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Revilla-Romero, B.; Beck, H.E.; Burek, P.; Salamon, P.; de Roo, A.; Thielen, J. Filling the gaps: Calibrating a rainfall-runoff model using satellite-derived surface water extent. Remote Sens. Environ. 2015, 171, 118–131. [Google Scholar] [CrossRef] [Green Version]
  41. Huang, C.; Peng, Y.; Lang, M.; Yeo, I.-Y.; McCarty, G. Wetland inundation mapping and change monitoring using Landsat and airborne LiDAR data. Remote Sens. Environ. 2014, 141, 231–242. [Google Scholar] [CrossRef]
  42. Zhang, Z.X.; Chen, X.; Xu, C.Y.; Yuan, L.F.; Yong, B.; Yan, S.F. Evaluating the non-stationary relationship between precipitation and streamflow in nine major basins of China during the past 50 years. J. Hydrol. 2011, 409, 81–93. [Google Scholar] [CrossRef]
  43. Xu, X.Y.; Yang, D.W.; Yang, H.B.; Lei, H.M. Attribution analysis based on the Budyko hypothesis for detecting the dominant cause of runoff decline in Haihe basin. J. Hydrol. 2014, 510, 530–540. [Google Scholar] [CrossRef]
  44. Wang, W.G.; Shao, Q.X.; Yang, T.; Peng, S.Z.; Xing, W.Q.; Sun, F.C.; Luo, Y.F. Quantitative assessment of the impact of climate variability and human activities on runoff changes: A case study in four catchments of the Haihe River basin, China. Hydrol. Process. 2013, 27, 1158–1174. [Google Scholar] [CrossRef]
  45. Xia, J.; Zeng, S.D.; Du, H.; Zhan, C.S. Quantifying the effects of climate change and human activities on runoff in the water source area of Beijing, China. Hydrol. Sci. J. 2014, 59, 1794–1807. [Google Scholar] [CrossRef] [Green Version]
  46. Yang, Z.F.; Zhang, W.G. Upstream-downstream cooperation approach in Guanting Reservoir watershed. J. Environ. Sci. China 2005, 17, 676–680. (In Chinese) [Google Scholar] [PubMed]
  47. Cheng, D.Z.; Chen, M.; Shi, S.P.; Ma, Z.Z. Impact of upstream human activities on the rainfall–runoff relationship in Yongding River. Des. Water Resour. Hydroelectr. Eng. 2001, 2, 19–21. (In Chinese) [Google Scholar]
  48. Yu, M.; Wang, C.R.; Liu, Y.; Olsson, G.; Wang, C.Y. Sustainability of mega water diversion projects: Experience and lessons from China. Sci. Total Environ. 2018, 619, 721–731. [Google Scholar] [CrossRef] [PubMed]
  49. Zhang, J.; Cao, X.S.; Meng, X.Z. Sustainable urban sewerage system and its application in China. Resour. Conserv. Recycl. 2007, 51, 284–293. [Google Scholar] [CrossRef]
Figure 1. Location of the Yongding River Basin (YDRB) study area with elevation from Shuttle Radar Topography (SRTM) Mission Digital Elevation Model (DEM).
Figure 1. Location of the Yongding River Basin (YDRB) study area with elevation from Shuttle Radar Topography (SRTM) Mission Digital Elevation Model (DEM).
Remotesensing 11 00560 g001
Figure 2. Workflow of this study.
Figure 2. Workflow of this study.
Remotesensing 11 00560 g002
Figure 3. Landsat-derived annual surface water extent dynamics during 1985–2016 in YDRB.
Figure 3. Landsat-derived annual surface water extent dynamics during 1985–2016 in YDRB.
Remotesensing 11 00560 g003
Figure 4. Statistics of reservoir water extents in YDRB.
Figure 4. Statistics of reservoir water extents in YDRB.
Remotesensing 11 00560 g004
Figure 5. Water extents in the Guanting and Cetian reservoirs.
Figure 5. Water extents in the Guanting and Cetian reservoirs.
Remotesensing 11 00560 g005
Figure 6. Water extent and inflow for (a) Guanting and (b) Cetian reservoirs.
Figure 6. Water extent and inflow for (a) Guanting and (b) Cetian reservoirs.
Remotesensing 11 00560 g006
Figure 7. Double mass curve of cumulative annual water extents against cumulative precipitation for (a) Guanting reservoir and (b) Cetian reservoir.
Figure 7. Double mass curve of cumulative annual water extents against cumulative precipitation for (a) Guanting reservoir and (b) Cetian reservoir.
Remotesensing 11 00560 g007
Figure 8. Temporal trends of (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) steel production, (g) crude iron production, (h) value of agriculture output and (i) urban area, compared to reservoir water extents in the Guanting reservoir.
Figure 8. Temporal trends of (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) steel production, (g) crude iron production, (h) value of agriculture output and (i) urban area, compared to reservoir water extents in the Guanting reservoir.
Remotesensing 11 00560 g008aRemotesensing 11 00560 g008b
Figure 9. Temporal trends of (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) value of agriculture output and (g) urban area, compared to reservoir water extents in the Cetian reservoir.
Figure 9. Temporal trends of (a) precipitation, (b) GDP, (c) population, (d) electricity energy production, (e) raw coal production, (f) value of agriculture output and (g) urban area, compared to reservoir water extents in the Cetian reservoir.
Remotesensing 11 00560 g009
Figure 10. Reservoir water extents extracted from Landsat-series imagery acquired before and after water delivery. The dates of the image acquisition are listed above the bar.
Figure 10. Reservoir water extents extracted from Landsat-series imagery acquired before and after water delivery. The dates of the image acquisition are listed above the bar.
Remotesensing 11 00560 g010
Table 1. Anthropogenic variables used in this study.
Table 1. Anthropogenic variables used in this study.
Anthropogenic VariablesUnitYearSource
GDPChinse Yuan (CNY)1985, 1990–2016Regional statistical yearbooks
PopulationPerson1985, 1990–2016
Electricity energy productionKilowatt hour1985, 1990, 1995–2016
Raw coal productionTon1985, 1990, 1995–2016
Steel productionTon1985, 1990, 1995–2016
Crude iron productionTon1985, 1990, 1995–2016
Value of agriculture outputCNY1985–2016
Urbankm21992–2015CCI-LC data
Table 2. Accuracies of the extracted reservoir water extents.
Table 2. Accuracies of the extracted reservoir water extents.
ReservoirYearReference ImagePA (%)UA (%)OA (%)Kappa Coefficient
Guanting2005QuickBird (2 m)93.497.298.80.95
2016Sentinel-2 (10 m)94.696.898.90.95
Yiliuhe2015QuickBird (2 m)98.097.399.70.97
Cetian2016Sentinel-2 (10 m)98.198.799.80.98
PA: producer’s accuracy; UA: user’s accuracy; OA: overall accuracy.
Table 3. Trend and change-point of reservoir water extents and streamflow from 1985 to 2016.
Table 3. Trend and change-point of reservoir water extents and streamflow from 1985 to 2016.
Reservoir TFPW-MK Trend TestPettit Change-Point Test
Slope (km2/year)ZYearP
GuangtingWater extent (km2)−1.56−4.86 **20010.00 **
Streamflow (m3/s)−0.28−4.21 **19990.00 **
CetianWater extent (km2)−0.116−2.24 *20050.03 *
Streamflow (m3/s)−0.067−2.45 *20050.01 **
* Significant at the 0.05 level; ** Significant at the 0.01 level.
Table 4. Temporal variation of the relative contributions of climate variability and human interventions to annual water extents in Guanting and Cetian reservoirs.
Table 4. Temporal variation of the relative contributions of climate variability and human interventions to annual water extents in Guanting and Cetian reservoirs.
ReservoirPeriodsΔA (km2)ΔAc (km2)ΔAh (km2)ΔA/A (%)Rc (%)Rh (%)
Guanting2002–2016−42.45.8 ± 16.4−48.1 ± 20.0−47.825.9 ± 22.474.1 ± 22.4
2002–2012−48.50.9 ± 15.3−49.3 ± 20.5−54.825.3 ± 25.374.7 ± 25.3
2013–2016−25.619.3 ± 12.0−44.9 ± 21.0−28.927.6 ± 13.672.4 ± 13.6
Cetian2006–2016−5.90.8 ± 3.1−6.7 ± 3.0−31.624.9 ± 17.375.1 ± 17.3
2006–2012−6.8−0.6 ± 2.3−6.2 ± 2.9−30.025.5 ± 20.074.5 ± 20.0
2013–2016−4.33.2 ± 3.0−7.5 ± 3.4−22.623.8 ± 14.376.2 ± 14.3
Note: ΔA is the difference in water extent between the reference and disturbance periods. ΔAc and ΔAh are the water extent changes attributed to climate variability and human activities, respectively. Rc and Rh are the relative contributions of climate variability and human activity to reservoir water extents.
Table 5. Pearson correlation coefficients between water extents and precipitation & human activity factors in Guanting reservoir and Cetian reservoir.
Table 5. Pearson correlation coefficients between water extents and precipitation & human activity factors in Guanting reservoir and Cetian reservoir.
Guanting ReservoirCetian Reservoir
1985–2001
(Pre-Change)
2002–2016
(Post-Change)
2002–20121985–2005
(Pre-Change)
2006–2016
(Post-Change)
2006–2012
Precipitation0.680 **0.4130.2290.610 **0.308−0.194
GDP0.2460.336−0.5640.2640.162−0.884 **
Population0.2010.300−0.5190.2020.136−0.755 *
Electricity energy production0.3730.246−0.5160.3730.004−0.786 *
Raw coal production−0.4920.258−0.3040.492 *0.126−0.771 *
Steel production0.020−0.005−0.833 **------
Crude iron production0.3340.082−0.722 **------
Value of agriculture output0.538 *0.43−0.5020.5360.269−0.740
Urban area0.2170.350−0.651 *0.2170.211−0.851 *
* Correlation is significant at the 0.05 level; ** Correlation is significant at the 0.01 level.
Table 6. The outflow and dates of water delivery from Cetian reservoir to Guanting reservoir.
Table 6. The outflow and dates of water delivery from Cetian reservoir to Guanting reservoir.
YearCetian Outflow
(×104 m3)
DateYearCetian Outflow
(×104 m3)
Date
2003501009/26–10/072010150010/15–11/16
2004714010/12–11/152011----
2005670310/17–11/062012----
2006157910/13–11/222013300010/08–11/11
2007260110/10–10/222014300011/07–11/22
2008180010/10–10/202015201809/29–11/23
2009----2016220010/17–10/30

Share and Cite

MDPI and ACS Style

Wang, M.; Du, L.; Ke, Y.; Huang, M.; Zhang, J.; Zhao, Y.; Li, X.; Gong, H. Impact of Climate Variabilities and Human Activities on Surface Water Extents in Reservoirs of Yongding River Basin, China, from 1985 to 2016 Based on Landsat Observations and Time Series Analysis. Remote Sens. 2019, 11, 560. https://doi.org/10.3390/rs11050560

AMA Style

Wang M, Du L, Ke Y, Huang M, Zhang J, Zhao Y, Li X, Gong H. Impact of Climate Variabilities and Human Activities on Surface Water Extents in Reservoirs of Yongding River Basin, China, from 1985 to 2016 Based on Landsat Observations and Time Series Analysis. Remote Sensing. 2019; 11(5):560. https://doi.org/10.3390/rs11050560

Chicago/Turabian Style

Wang, Mingli, Longjiang Du, Yinghai Ke, Maoyi Huang, Jing Zhang, Yong Zhao, Xiaojuan Li, and Huili Gong. 2019. "Impact of Climate Variabilities and Human Activities on Surface Water Extents in Reservoirs of Yongding River Basin, China, from 1985 to 2016 Based on Landsat Observations and Time Series Analysis" Remote Sensing 11, no. 5: 560. https://doi.org/10.3390/rs11050560

APA Style

Wang, M., Du, L., Ke, Y., Huang, M., Zhang, J., Zhao, Y., Li, X., & Gong, H. (2019). Impact of Climate Variabilities and Human Activities on Surface Water Extents in Reservoirs of Yongding River Basin, China, from 1985 to 2016 Based on Landsat Observations and Time Series Analysis. Remote Sensing, 11(5), 560. https://doi.org/10.3390/rs11050560

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop