Next Article in Journal
Soil Hydrology for a Sustainable Land Management: Theory and Practice
Next Article in Special Issue
Changes in the Water Temperature of Rivers Impacted by the Urban Heat Island: Case Study of Suceava City
Previous Article in Journal
Intertwining Observations and Predictions in Vadose Zone Hydrology: A Review of Selected Studies
Previous Article in Special Issue
An Improved Model for the Evaluation of Groundwater Recharge Based on the Concept of Conservative Use Potential: A Study in the River Pandeiros Watershed, Minas Gerais, Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tracking Lake and Reservoir Changes in the Nenjiang Watershed, Northeast China: Patterns, Trends, and Drivers

1
Key Laboratory of Wetland Ecology and Environment, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun 130102, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
National Earth System Science Data Center, Beijing 100101, China
4
School of Management Engineering, Qingdao University of Technology, Qingdao 266033, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(4), 1108; https://doi.org/10.3390/w12041108
Submission received: 17 March 2020 / Revised: 30 March 2020 / Accepted: 31 March 2020 / Published: 13 April 2020

Abstract

:
In terms of evident climate change and human activities, investigating changes in lakes and reservoirs is critical for sustainable protection of water resources and ecosystem management over the Nenjiang watershed (NJW), an eco-sensitive semi-arid region and the third-largest inland waterbody cluster in China. In this study, we established a multi-temporal dataset documenting lake and reservoir (area ≥ 1 km2) changes in this region using an object-oriented image classification method and Landsat series images from 1980 to 2015. Using the structural equation model (SEM), we analyzed the diverse impacts of climatic and anthropogenic variables on lake changes. Results indicated that lakes experienced significant changes with fluctuations over the past 35 years including obvious declines in the total area (by 42%) and number (by 51%) from 1980 to 2010 and a slight increase in the total lake area and number from 2010 to 2015. More than 235 lakes in the size class of 1–10 km2 decreased to small lakes (area < 1 km2), while 59 lakes covering 243.75 km2 disappeared. Total reservoir area and number had continuous increases during the investigated 35 years, with an areal expansion of 54.9% from 919 km2 to 1422 km2, and a number increase by 65.3% from 78 to 129. The SEM revealed that the lake area in the NJW had a significant correlation with the mean annual precipitation (MAP), suggesting that the MAP decline clarified most of the lake shrinkage in the NJW. Furthermore, agricultural consumption of water had potential impacts on lake changes, suggested by the significant relationship between cropland area and lake area.

1. Introduction

Inland waterbodies (lakes and reservoirs) are closely related to human life due to their diverse ecosystem services [1,2]. Lakes and reservoirs in the drylands, although covering only a small proportion of the landscape, play irreplaceable roles in fragile environments and for local residents [3]. Meanwhile, they are sensitive to climate change and human disturbances [4]. Monitoring changes in lakes and reservoirs and investigating their driving forces are thus of great significance to sustainable water resource management and regional economic development.
Due to diverse driving forces, lakes around the world have experienced changes in both area and number during the past decades. Thaw and “breaching” of permafrost have caused a widespread decline in the Arctic lake number and area [5]. Because of warmer air temperatures, which allow higher evaporation rates, as the world’s largest surface freshwater system, the total area of the Great Lakes of North America experienced a continuous decrease from the 1980s to 2005 [6]. Under the background of global warming, glaciers in China’s Tibetan Plateau (TP) showed a general retreat, while the glacial lakes expanded notably [7]. However, a drier climate and degraded permafrost have led to lake shrinkage in some basins of the TP [8]. Lakes over Mongolia experienced areal decline mostly related to a drier climate, while lakes in Inner Mongolia had notable shrinkage which was attributed mainly to human consumption of water resources, particularly coal mining [9]. Moreover, the lake area in the Jianghan Plain of China decreased dramatically from the 1950s to 1998, which was mainly caused by agricultural cultivation.
During the past 60 years, the number of reservoirs has increased notably, and there are about 50,000 large reservoirs nowadays around the world [10]. For example, China constructed nearly 45,000 reservoirs in the Yangtze River Basin to meet the large demand for water resources. By 2013, 98,000 reservoirs had been built in China [11]. Previous studies [12,13] indicated that most of large river systems and lakes in China were affected by reservoirs. Due to the construction of large reservoirs in upstream areas, water was impounded upstream, which seriously affected the water supply of lakes in the middle and lower reaches [14]. Therefore, it is necessary to track and analyze the changes in lakes and reservoirs at different scales and regions to form a scientific management response [15,16].
Optical images from different satellite sensors were widely used to monitor spatiotemporal variations of inland waterbodies at multiple geographic scales [17,18,19,20,21,22]. The Moderate Resolution Imaging Spectroradiometer (MODIS) data have been widely used to assess the water extent at daily to 16-day timescales despite a spatial resolution at 250 or 500 m [23]. Yet, areal changes of small lakes or reservoirs with irregular shapes were not accurately delineated due to the coarse resolution of source data. Other optical sensors, such as Quickbird (DigitalGlobe, Longmont, Colorado, USA) and IKONOS (DigitalGlobe, Longmont, Colorado, USA), provide finer images comparable to aerial photography for the extraction of lake or reservoir boundaries [24]. However, those data are limited in application at broad scale due to the high costs, narrow swath size, and so on [25,26]. Landsat series images have a fine spatial resolution (30 or 80 m) compared with previously mentioned data and have provided the longest temporal and spatial records for surface observations since their first launch in 1972 [27]. Consequently, Landsat imagery has been widely-used remote sensing data in examining changes in lakes or reservoirs.
As one of the third-largest waterbody clusters, an eco-environmentally fragile area, and an important base for grain production, the Nenjiang watershed (NJW) plays an important role in ecological conservation and national grain security in China [28]. However, the changes in lakes and reservoirs across the NJW during recent decades have rarely been examined and their correlations with climate change and human disturbances have rarely been quantitatively investigated. Doing so is critical to understanding the regional water cycle and sustainable water resource management. Therefore, this study aims to (1) employ Landsat 8 images to investigate the current status of lakes and reservoirs in the NJW, (2) evaluate changes in area and number of lakes and reservoirs using consistent Landsat images (i.e., 1980, 1990, 2000, 2010, and 2015), and (3) quantify the roles of climatic factors and artificial variables in driving lake changes.

2. Materials and Methods

2.1. Study Area

The NJW, covering an area of 297 × 103 km2, is located in the core region of Northeast China (Figure 1), with latitudes from 44°1′48″ N to 51°42′1″ N and longitudes from 119°12′1″ E to 127°54′2″ E [28]. Obvious terrain variances can be found in this area with elevations ranging from 120 to 1740 m above sea level. The highest elevation is in the northwest at the Greater and Lesser Khingan Mountains, and the lowest elevation is in the southeast at the Songnen Plain. This is why the lakes and reservoirs are mainly observed in the southeast. The NJW is dominated by a temperate semi-arid continental climate with an mean annual precipitation of around 470 mm and a mean annual air temperature of 4 °C [29]. Most precipitation (82%) occurs in months from June to September which is a critical growth period for vegetation in this region. As an important grain base in China, the NJW plays an critical role in promoting regional economic development and ensuring national food security [30,31]. Moreover, this region provides important habitats for threatened species and migratory waterfowls.

2.2. Data Source and Processing

2.2.1. Satellite Data

In this study, multi-temporal Landsat images from multispectral scanner (MSS), thematic mapper (TM), enhanced thematic mapper plus (ETM+), and operational land imager (OLI) sensors were acquired from the United States Geological Survey (USGS) to investigate the changes of lakes and reservoirs from 1980 to 2015. Specifically, we used 28, 27, 27, 29, and 26 scenes of images to extract the lakes and reservoirs for five dates of 1980, 1990, 2000, 2010, and 2015, respectively, (Figure 2). A total of 137 images acquired for months from June to September were used for the classification because the vegetation in this period is easier to be identified and the lakes and reservoirs have the largest amount. All of the images have a little cloud cover (less than 5%) and lakes/reservoirs on these images are clearly visible. In addition, in order to analyze the impact of floods on the area and the number of lakes and reservoirs in 1998 and 2013, we used 30 scenes of TM and 30 scenes of OLI images to extract the lakes and reservoirs in 1998 and 2013, respectively. Prior to image classification, data preprocessing including geometric, topographic, and radiometric corrections was performed for all the images using the ENVI 5.3 (Exelis Visual Information Solutions, Inc. Boulder, USA) software package. To ensure the data consistency, all images were re-projected to the 1984 WGS UTM zone 51N projection.

2.2.2. Meteorological Data

In order to analyze the relationships of lake/reservoir area with climatic factors, the daily climatic data, including extreme air temperature and precipitation, mean wind speed, and sunshine hours during 1980–2015, were collected from the meteorological records of the China Meteorological Data Service Center. Figure 1 shows the locations of the meteorological stations in the NJW. Spatial patterns of the mean annual air temperature (MAAT) and mean annual precipitation (MAP) were interpolated from those meteorological records using an Anusplin software considering the elevation differences [32]. We used the daily meteorological data and the Food and Agriculture Organization of the United Nations (FAO) Penman–Monteith model [8,33] to calculate the actual evapotranspiration (ET). ET is calculated as follows:
E T 0 = 0.408 Δ ( R n G ) + γ 900 T + 273 U 2 ( e s e a ) Δ + γ ( 1 + 0.34 U 2 )
E T = 9.78 + 0.0072 × E T 0 × P P T + 0.051 × P P T × L A I
where E T 0 is the potential evapotranspiration, R n is the net radiation, G is the soil heat flux, e s is the saturation vapor pressure, e a is the actual vapor pressure, e s e a is the saturation vapor pressure deficit of the air, T is the air temperature at 2 m height, U 2 is the wind speed at 2 m height, Δ is the slope of the saturation vapor pressure temperature relationship, γ is the psychrometric constant, P P T is precipitation, and the calculation of L A I is referenced from the method of Lu et al. [34]. These parameters are directly calculated or derived from the average daily maximum and minimum temperature, the daily average temperature, the daily actual vapor pressure, the daily wind speed data, the actual duration of sunshine hours, the relative humidity data, and other empirical metrics.
In order to analyze the potential impact of climate change in the next 35 years (2015–2050) on regional lake area and number changes, we downloaded the Representative Concentration Pathways 5 (RCP 5) datasets from the China Agrosys Platform (http://stdown.agrivy.com). The future climate change data were generated based on the Fifth Generation Coupled Global Climate Model (CGCM 5) from the Canadian Centre for Climate Modeling and Analysis. The datasets constructed using original meteorological observations for each station included daily mean temperatures, daily precipitation, and daily potential evaporation. We converted daily mean temperature, daily precipitation, and daily potential evapotranspiration into annual mean temperature, annual mean precipitation, and annual potential evapotranspiration for our analysis.

2.2.3. Other Data

The areal data of the cropland area and gross domestic product (GDP) index collected from local statistical yearbook (http://www.stats.gov.cn).

2.3. Data Analysis

2.3.1. Extracting Lakes and Reservoirs

We focused on lakes and reservoirs greater than 1 km2 in terms of the 30 and 80 m resolutions of the satellite images [35]. An object-based image analysis (OBIA) was used rather than the traditional pixel-based classification method, because OBIA classifies objects instead of individual pixels [36,37,38]. In the process of OBIA classification, the spectrum, spatial information, texture, and geometric features characterized by remote sensing images were fully utilized. The lake and reservoir extents of the NJW in 1980, 1990, 1998, 2000, 2010, 2013, and 2015 were extracted in eCognition Developer 8.6 [14,39]. The normalized difference water index (NDWI) is the most popular used index for automated inland waterbodies delineation [35,40,41]. In addition, we tested different versions of NDWI (mNDWI, NDRW, and NDWI) [21,42] and found that the selected NDWI is more effective in our study area. Therefore, NDWI was applied to extract lakes and reservoirs, which was defined as:
N D W I = ( G r e e n N I R ) / ( G r e e n + N I R )
where G r e e n and N I R represent the reflectance of the green and Near Infrared (NIR) bands, respectively [25]. The data processing steps for the lake and reservoir inventory are shown in Figure 3. The extraction of lakes and reservoirs consisted of three major steps: image multi-resolution segmentation, NDWI threshold testing, and classification rule designing. By defining the length/width index and the rectangular fit index built into the software, rivers and artificial ponds were removed. Then, we used visual interpretation to extract the reservoirs.
To validate the accuracy of extracted lakes and reservoirs, we calculated error matrices based on 1492 validation samples selected from Google Earth and Landsat sensors. The detailed number of validation samples for the five dates is given in Table 1. The overall accuracies of the lake and reservoir classification for the five dates were respectively evaluated and the standard errors bars (uncertainties) were estimated [43]. The overall accuracies of lakes and reservoirs were larger than 90%, and the kappa coefficients were larger than 0.87.

2.3.2. Temporal Analysis of Lakes and Reservoirs

To fully understand the changes of lakes and reservoirs of different sizes, the lakes and reservoirs were categorized into four classes: 1–10 km2, 10–50 km2, 50–100 km2, and >100 km2. In order to examine the areal change rate of the lake or reservoir in different periods, an indicator of the lake or reservoir area dynamic degree, as shown Equation (4), was used to analyze their changes [44,45].
K = U b U a U a × 1 T × 100 %
where K is the dynamic indicator for the lake or reservoir area; U a and U b are the area of the lake and reservoir at the start date and end date, respectively; and T is the time scale under consideration.

2.3.3. Assessing the Roles of Climatic Factors and Anthropogenic Causes in Lake Changes

The roles of climatic factors and artificial variables in changes of lakes during the study period were quantified using structural equation modeling (SEM) by the AMOS 22 software (IBM, Armonk, NY, USA) [39]. This paper analyzed the changes in MAAT, MAP, and ET to examine the influences of climate change on lake changes. The statistical data of cropland area and GDP were selected to investigate the impacts of artificial variables on lake changes [46,47]. Table 2 illustrates the optimum values for these indicators necessary for the SEM.

3. Results

3.1. Spatial Pattern of Lakes and Reservoirs in 2015

Figure 4 shows the lake and reservoir distribution in 2015 across the NJW. Lakes and reservoirs were identified dominantly in the southeastern part of the NJW. A total of 233 lakes (area ≥ 1 km2) covering an area of 2110 ± 53 km2 in 2015 were extracted from satellite images. Most of the lakes (89.4%) had an area of smaller than 10 km2. There were three lakes with area greater than 100 km2, and their accumulated area accounted for approximately 33.8% of the total lake area in the NJW. The Chagan Lake was the largest lake with an estimated area of 292 km2.
There were 129 reservoirs (area ≥ 1 km2) in the NJW with a total area of 1422 ± 56 km2 in 2015. Of these, 108 reservoirs had an area smaller than 10 km2. There were two reservoirs with an area larger than 100 km2, with the accumulated area accounting for approximately 41.5% of the total reservoir area in the NJW. The Nierji Reservoir was the largest reservoir with an estimated area of 360 km2.

3.2. Temporal Changes of Lakes and Reservoirs from 1980 to 2015

Figure 5 suggests that dramatic changes in the total area and number of lakes and reservoirs over the NJW occurred from 1980 to 2015. During the observed 35 years, the total lake area decreased by 38.7% from 3440 ± 72 km2 to 2110 ± 53 km2, whilst the total number of lakes decreased by 233 from the original 468 in 1980. Specifically, lake changes had significant fluctuations over the past 35 years, including obvious declines in total area (42%) and number (51%) from 1980 to 2010 and slight increases in the total lake area and number from 2010 to 2015. Reservoirs in the NJW experienced continuous expansion during 1980–2015. The total number of reservoirs increased from 78 to 129 with the total area expansion being 55%, from 919 ± 53 km2 to 1422 ± 56 km2.

3.2.1. Lake Changes

For a further understanding of the temporal changes of lakes, the spatial heterogeneity of lake area changes across the NJW was investigated (Figure 6). Changes in the lake area presented clear variations. During the first period, 1980–1990, most of the expanded lakes were mainly distributed in the southeastern part of the NJW, while the lakes that shrunk were mainly distributed in the central part (Figure 6a). A rapid shrinkage of the total lake area and a decline in lake number occurred during 1990–2000. The larger lake shrinkage occurred in the southern NJW (Figure 6b), while lakes in the Chagan Lake, Yangsha Lake, and Tumuji Lake zones exhibited expanding trends during this period. During 2000–2010, lakes in the NJW with large areal loss were distributed mainly in the Dalong Lake zones, while expanded lakes were identified mainly at the intersection point between the Nenjiang River and Taoerhe River (Figure 6c). During 2010–2015, it is noteworthy that the characteristic of lake expansion was relatively obvious in the NJW compared to the other periods. Lakes with shrinkage were mainly distributed in the southeastern part, whereas lake expansions mostly occurred in the eastern part (Figure 6d).
The detailed lake changes in each class are shown in Table 3. From 1980 to 1990, the largest areal decline of lakes was observed for the 50–100 km2 size class. Interestingly, the Dalong Lake, which was in the size class of 50–100 km2, separated into five lakes in the size class of 10–50 km2 during this period. Specifically, 10 lakes with a size of 1–10 km2 covering a total area of 34 km2 in the central region disappeared in this period. During 1990–2000, lake shrinkage occurred with the total lake area declining by 18%, and 53 lakes (area ≤ 50 km2) vanished during this decade. During 2000–2010, both the total number and area of lakes in the size class of 50–100 km2 increased, while both the total number and area of lakes in other size classes decreased. From 2010 to 2015, both the total area and number of lakes in the size classes 1–10, 10–50, and 50–100 km2, increased.

3.2.2. Reservoir Changes

Figure 7 shows the spatiotemporal patterns of reservoir changes in the NJW during different periods. During 1980–1990, the reservoirs that shrunk were distributed mainly in the central and northeast parts, whereas the expanded reservoirs were mainly distributed in the southern and eastern parts of the NJW (Figure 7a). During 1990–2000, the expanded reservoirs were mainly distributed in the northeast and southern parts, whereas the reservoirs that shrunk were mainly distributed in the central and southern parts of the NJW (Figure 7b). During 2000–2010, significant reservoir expansions were identified across the NJW, while the larger reservoirs that shrunk were mainly distributed in the southern part of the NJW (Figure 7c). After 2010, the reservoirs that shrunk were mainly distributed in the southern and eastern parts. Larger expanded reservoirs were mainly identified in the northern part and the Taoer River basin (Figure 7d).
During the period of 1980–1990, the total area of reservoirs in size classes of 10–50 km2 and 50–100 km2 increased by 121%, with an areal increase from 149 ± 10 to 328 ± 16 km2, and by 133% from 52 ± 14 to 121 ± 17 km2, respectively. The total areas of reservoirs in the size classes of 1–10 km2 and >100 km2 decreased from 205 ± 9 to 200 ± 6 km2 (−2%) and from 513 ± 20 to 151 ± 14 km2 (−71%), respectively (Table 4). During 1990–2000, the total areas of reservoirs in the size classes of 1–10 km2 and 50–100 km2 increased by 23%, with an areal increase from 200 ± 6 to 247 ± 5 km2, and by 123% from 121 ± 17 to 270 ± 15 km2, respectively. The total area of reservoirs in the size classes of 10–50 km2 decreased from 328 ± 16 to 263 ± 8 km2 (−20%). During the period of 2000–2010, the total area of reservoirs increased from 780 ± 28 to 1086 ± 58 km2. The total areas of reservoirs in the size classes of 1–10 km2 and 10–50 km2 increased by 30% and 24%, respectively. The total area of reservoirs in the size class of 50–100 km2 decreased by 50%. During the period of 2010–2015, the total reservoir area increased by 31% with an increase in area from 1086 ± 58 to 1422 ± 56 km2. Detailed reservoir changes for these different classes are shown in Table 4.

3.3. Roles of Climatic Factors and Artificial Variables in Driving Lake Changes

As shown in Figure 8, the model passed the reliability test, convergent validity test, and discriminant validity test. The dominant climatic factor that influenced lake changes was the MAP (β = 0.66, p < 0.001), followed by the ET (β = 0.39, p < 0.05) and MAAT (β = 0.15, p < 0.1). Agricultural consumption of water had a significant effect on lake changes, suggested by the significant relationship between cropland area and lake area (β = 0.17, p < 0.1).

4. Discussion

Lake and reservoir mapping is affected by the resolution of the used satellite data [26]. We focused only on lakes and reservoirs with area greater than 1 km2 to reduce as much as possible the errors induced by a coarse resolution of 30 and 80 m. Although some images out of the optimal season were used, these data mainly covered the Greater and Lesser Khingan Mountains where few lakes and reservoirs were identified. This did not yield large uncertainties in our analysis. This study integrated OBIA and visual interpretation instead of automatic classification to extract lakes and reservoirs, which ensured data accuracy and effective analysis.
Both climate change and human activities contributed to the changes of lakes in the NJW. On the one hand, the climate over the study area is a semi-arid continental climate. Water supply and output for these lakes dominated with precipitation and evapotranspiration, respectively. SEM analysis revealed that the lake shrinkage in the NJW had a significant correlation with the MAP (β = 0.66, p < 0.001) (Figure 8), followed by the ET (β = 0.39, p < 0.05), and the MAAT (β = 0.15, p < 0.1). This suggests that precipitation had a significant statistical relationship with lake area and potentially had the largest impacts on the lake area. During the investigated 35 years, a warmer climate was identified for the NJW with a significant increase of MAAT (p < 0.05) (Figure 9a). A reduced water supply from precipitation and output by increased ET characterized most of the lake shrinkages in the NJW. Specifically, lakes in the NJW showed areal changes with a decrease from 1980 to 2010 and then an increase from 2010 to 2015 (Figure 5), which is consistent with the changed trend of MAP during the 35 years (Figure 9a). The MAP could be regarded as the main climatic factor to explain the lake shrinkage in the NJW from 1980 to 2015 in terms of its decline.
Future projections of climate change using the Intergovernmental Panel on Climate Change (IPCC) models (RCP 5) indicate that the warming and drying trend will continue in the NJW (Figure 9b). Both the MAAT and PET show increasing trends with a rate of 0.03 ℃ yr−1 and 1.64 mm yr−1 (p < 0.05), while the MAP exhibits a significant decline with a rate of −1.63 mm yr−1 (p < 0.01). If this trend continues, the total area of lakes in the NJW may continue to decrease, and some small lakes will disappear.
This study found that flood events can markedly affect the lake area and number, especially for the small lakes (1–10 km2). According to the hydrological records, the most serious floods during the recent century occurred in 1998 and 2013 [51]. Compared with the total lake area in 2010 (normal flow year), the total lake area in the NJW respectively increased by 2946 km2 in 1998 and 895 km2 in 2013 during the flooding events. It is clear that extreme precipitation events (floods) have accelerated the expansion of the total lake area (Figure 10). In arid and semi-arid regions, a multidimensional view on the prevention and exploitation of floods is required. Lakes and reservoirs have huge water storage capacity, and thus they could serve as hydrological buffers to prevent floods and provide irrigation water for sustainable agricultural development [10].
On the other hand, the study area is an important grain production base in China. Agricultural demand for water had potential impacts on lakes. Besides climate change, human activities have imposed marked impacts on lakes and reservoirs [52]. While our study revealed that climatic factors drive striking lake and reservoir changes, human-induced changes should be responded to as quickly as possible. In this study, a clear correlation (β = 0.17, p < 0.1) was observed between lake area and crop area (Figure 8). Lake shrinkage caused by artificial variables is mainly attributed to the agricultural water consumption in the NJW [50,53]. Due to the construction of large reservoirs upstream water was therefore impounded upstream, which seriously affected the water supply of lakes in the middle and lower reaches. In particular, the largest reservoir, Nierji Reservoir, constructed over the upstream area of the Nenjiang River exerted a marked influence on downstream lakes. Due to reclamation from lakes and agricultural development during the past 35 years, some lakes in the size class of 1–10 km2 decreased to small lakes with areas smaller than 1 km2. For example, a large area of shallow waters was reclaimed for planting rice (Figure 11). Since 1990, the area of paddy fields in the NJW has increased significantly by 90.26 km2 from 3.46 to 93.73 km2. The total cropland area increased from 459.42 to 948.18 km2. With the population increase and agricultural land expansion in this region, more and more open water resources have been applied to agricultural irrigation. It is evident that the human impacts on lake and reservoir changes could not be ignored.
Lake shrinkage and disappearance lead to ecological and environmental degradation, such as aggravating the degree of sandstorms and desertification and reducing the number of wild animals [4,54]. To optimize water distribution programs through scientific water conservancy projects is very important. Therefore, appropriate measures need to be implemented by managers to further reduce the decline of lake areas [55], such as enhancing the drainage capacity of reservoirs in drought years but the storage capacity of reservoirs in flood years. In addition, we should control the expansion of paddy fields to relieve water stress and guarantee regional sustainable development in the NJW [56].

5. Conclusions

In this paper, we established a multi-temporal dataset of lakes and reservoirs in the NJW using long time-series Landsat images from 1980 to 2015 to document their changes on a decadal scale and quantified the contribution degree of MAAT, MAP, ET, cropland, and GDP to change in lake area. A notable decline in the total lake area by 1330 km2 in the period of 1980 to 2015 was identified, while the lake number decreased contemporaneously. In contrast to lake shrinkage, the total area and number of reservoirs in the NJW experienced continuous increases. We identified 51 newborn reservoirs with total area of 504 km2 in 2015 compared to 1980. SEM analysis revealed that decrease of the MAP is the dominant factor driving the changes of lakes, followed by ET and MAAT, especially for those of small size. Furthermore, the human impacts on lake and reservoir changes could not be ignored. Timely and appropriate policies and measures are required to reduce lake shrinkage and respond to the degraded environment. The results and analysis in this study are expected to provide guidance for sustainable management of water resource in the NJW.

Author Contributions

B.D. and D.M. designed the analytical framework of this study. B.D. performed the data analysis and drafted the manuscript. H.L. and H.X. provided methodological advice. D.M. and Z.W. made major revisions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was jointly funded by the National Key Research and Development Program of China (2019YFA0607100, 2016YFC0500408, 2016YFC0500201, and 2016YFA0602301), the National Natural Science Foundation of China (41771383), and the funding from Youth Innovation Promotion Association, Chinese Academy of Sciences (2017277, 2012178) and National Earth System Science Data Center of China (http://northeast.geodata.cn).

Acknowledgments

We appreciate the facility that make Landsat 8 OLI images accessible through the USGS. We thank the three anonymous reviewers for the constructive comments and suggestions, which help improve the quality of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bai, J.; Chen, X.; Li, J.; Yang, L.; Fang, H. Changes in the area of inland lakes in arid regions of central Asia during the past 30 years. Environ. Monit. Assess. 2011, 178, 247–256. [Google Scholar] [CrossRef] [PubMed]
  2. Li, B.; Yang, G.; Wan, R.; Zhang, L.; Zhang, Y.; Dai, X. Using fuzzy theory and variable weights for water quality evaluation in Poyang Lake, China. Chin. Geogr. Sci. 2017, 27, 39–51. [Google Scholar] [CrossRef] [Green Version]
  3. Duan, Z.; Bastiaanssen, W.G.M. Estimating water volume variations in lakes and reservoirs from four operational satellite altimetry databases and satellite imagery data. Remote Sens. Environ. 2013, 134, 403–416. [Google Scholar] [CrossRef]
  4. Williamson, C.E.; Saros, J.E.; Vincent, W.F.; Smol, J.P. Lakes and reservoirs as sentinels, integrators, and regulators of climate change. Limnol. Oceanogr. 2009, 54, 2273–2282. [Google Scholar] [CrossRef]
  5. Smith, L.C.; Sheng, Y.; MacDonald, G.; Hinzman, L. Disappearing arctic lakes. Science 2005, 308, 1429. [Google Scholar] [CrossRef] [Green Version]
  6. McBean, E.; Motiee, H. Assessment of impact of climate change on water resources: A long term analysis of the Great Lakes of North America. Hydrol. Earth Syst. Sci. 2008, 12, 239–255. [Google Scholar] [CrossRef] [Green Version]
  7. Zhiguo, L. Glacier and lake changes across the Tibetan Plateau during the past 50 years of climate change. J. Resour. Ecol. 2014, 5, 123–132. [Google Scholar] [CrossRef]
  8. Mao, D.; Wang, Z.; Yang, H.; Li, H.; Thompson, R.J.; Li, L.; Song, K.; Chen, B.; Gao, H.; Wu, J. Impacts of Climate Change on Tibetan Lakes: Patterns and Processes. Remote Sens. 2018, 10, 358. [Google Scholar] [CrossRef] [Green Version]
  9. Tao, S.; Fang, J.; Zhao, X.; Zhao, S.; Shen, H.; Hu, H.; Tang, Z.; Wang, Z.; Guo, Q. Rapid loss of lakes on the Mongolian Plateau. Proc. Natl. Acad. Sci. USA 2015, 112, 2281–2286. [Google Scholar] [CrossRef] [Green Version]
  10. Lehner, B.; Liermann, C.R.; Revenga, C.; Vörösmarty, C.; Fekete, B.; Crouzet, P.; Döll, P.; Endejan, M.; Frenken, K.; Magome, J.; et al. High-resolution mapping of the world’s reservoirs and dams for sustainable river-flow management. Front. Ecol. Environ. 2011, 9, 494–502. [Google Scholar] [CrossRef] [Green Version]
  11. Ministry of Water Resources. Bulletin of First National Census for Water; China Water and Power Publisher: Beijing, China, 2013. [Google Scholar]
  12. Yang, X.; Lu, X. Drastic change in China’s lakes and reservoirs over the past decades. Sci. Rep. 2014, 4, 6041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mao, D.; Wang, Z.; Wu, J.; Wu, B.; Zeng, Y.; Song, K.; Yi, K.; Luo, L. China’s wetlands loss to urban expansion. Land Degrad. Dev. 2018, 29, 2644–2657. [Google Scholar] [CrossRef]
  14. Chang, B.; Li, R.; Zhu, C.; Liu, K. Quantitative Impacts of Climate Change and Human Activities on Water-Surface Area Variations from the 1990s to 2013 in Honghu Lake, China. Water 2015, 7, 2881–2899. [Google Scholar] [CrossRef] [Green Version]
  15. Zheng, Z.; Li, Y.; Guo, Y.; Xu, Y.; Liu, G.; Du, C. Landsat-Based Long-Term Monitoring of Total Suspended Matter Concentration Pattern Change in the Wet Season for Dongting Lake, China. Remote Sens. 2015, 7, 13975–13999. [Google Scholar] [CrossRef] [Green Version]
  16. Wu, G.; Liu, Y. Mapping Dynamics of Inundation Patterns of Two Largest River-Connected Lakes in China: A Comparative Study. Remote Sens. 2016, 8, 560. [Google Scholar] [CrossRef] [Green Version]
  17. 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. [Google Scholar] [CrossRef] [Green Version]
  18. 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]
  19. 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]
  20. Beaton, A.; Whaley, R.; Corston, K.; Kenny, F. Identifying historic river ice breakup timing using MODIS and Google Earth Engine in support of operational flood monitoring in Northern Ontario. Remote Sens. Environ. 2019, 224, 352–364. [Google Scholar] [CrossRef]
  21. Carroll, M.; Wooten, M.; DiMiceli, C.; Sohlberg, R.; Kelly, M. Quantifying Surface Water Dynamics at 30 Meter Spatial Resolution in the North American High Northern Latitudes 1991–2011. Remote Sens. 2016, 8, 622. [Google Scholar] [CrossRef] [Green Version]
  22. Yamazaki, D.; Trigg, M.A.; Ikeshima, D. Development of a global ~90m water body map using multi-temporal Landsat images. Remote Sens. Environ. 2015, 171, 337–351. [Google Scholar] [CrossRef]
  23. Li, Q.; Lu, L.; Wang, C.; Li, Y.; Sui, Y.; Guo, H. MODIS-Derived Spatiotemporal Changes of Major Lake Surface Areas in Arid Xinjiang, China, 2000–2014. Water 2015, 7, 5731–5751. [Google Scholar] [CrossRef]
  24. Karlsson, J.M.; Jaramillo, F.; Destouni, G. Hydro-climatic and lake change patterns in Arctic permafrost and non-permafrost areas. J. Hydrol. 2015, 529, 134–145. [Google Scholar] [CrossRef] [Green Version]
  25. Wendleder, A.; Friedl, P.; Mayer, C. Impacts of Climate and Supraglacial Lakes on the Surface Velocity of Baltoro Glacier from 1992 to 2017. Remote Sens. 2018, 10, 1681. [Google Scholar] [CrossRef] [Green Version]
  26. Khadka, N.; Zhang, G.; Thakuri, S. Glacial Lakes in the Nepal Himalaya: Inventory and Decadal Dynamics (1977–2017). Remote Sens. 2018, 10, 1913. [Google Scholar] [CrossRef] [Green Version]
  27. Mishra, N.; Haque, M.O.; Leigh, L.; Aaron, D.; Helder, D.; Markham, B. Radiometric Cross Calibration of Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper Plus (ETM+). Remote Sens. 2014, 6, 12619–12638. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, Z.; Wang, Z.; Zhang, B.; Lu, C.; Ren, C. Impact of land use/land cover changes on ecosystem services in the Nenjiang River Basin, Northeast China. Ecol. Process. 2015, 4, 11. [Google Scholar] [CrossRef] [Green Version]
  29. Feng, X.; Zhang, G.; Yin, X. Hydrological responses to climate change in Nenjiang river basin, northeastern China. Water Resour. Manag. 2011, 25, 677–689. [Google Scholar] [CrossRef]
  30. Qiu, B.; Lu, D.; Tang, Z.; Chen, C.; Zou, F. Automatic and adaptive paddy rice mapping using Landsat images: Case study in Songnen Plain in Northeast China. Sci. Total Environ. 2017, 598, 581–592. [Google Scholar] [CrossRef]
  31. Zhang, Y.; Zang, S.; Sun, L.; Yan, B.; Yang, T.; Yan, W.; Meadows, M.E.; Wang, C.; Qi, J. Characterizing the changing environment of cropland in the Songnen Plain, Northeast China, from 1990 to 2015. J. Geogr. Sci. 2019, 29, 658–674. [Google Scholar] [CrossRef] [Green Version]
  32. Liu, Z.-H.; McVicar, T.R.; LingTao, L.; Van Niel, T.G.; Yang, Q.-K.; Li, R.; Mu, X.-M. Interpolation for time series of meteorological variables using ANUSPLIN. J. Northwest A F Univ. 2008, 36, 227–234. [Google Scholar]
  33. Yoo, S.-H.; Choi, J.-Y.; Jang, M.-W. Estimation of design water requirement using FAO Penman–Monteith and optimal probability distribution function in South Korea. Agric. Water Manag. 2008, 95, 845–853. [Google Scholar] [CrossRef]
  34. Lü, Y.; Fu, B.; Feng, X.; Zeng, Y.; Liu, Y.; Chang, R.; Sun, G.; Wu, B. A policy-driven large scale ecological restoration: Quantifying ecosystem services changes in the Loess Plateau of China. PLoS ONE 2012, 7, e31782. [Google Scholar] [CrossRef] [PubMed]
  35. Korzeniowska, K.; Korup, O. Object-Based Detection of Lakes Prone to Seasonal Ice Cover on the Tibetan Plateau. Remote Sens. 2017, 9, 339–361. [Google Scholar] [CrossRef] [Green Version]
  36. Walter, V. Object-based classification of remote sensing data for change detection. ISPRS J. Photogramm. 2004, 58, 225–238. [Google Scholar] [CrossRef]
  37. Blaschke, T.; Lang, S.; Hay, G. Object-Based Image Analysis: Spatial Concepts for Knowledge-Driven Remote Sensing Applications; Springer: Berlin, Germany, 2008. [Google Scholar]
  38. Jia, M.; Wang, Z.; Wang, C.; Mao, D.; Zhang, Y. A new vegetation index to detect periodically submerged Mangrove forest using single-tide sentinel-2 imagery. Remote Sens. 2019, 11, 2043. [Google Scholar] [CrossRef] [Green Version]
  39. Li, H.; Mao, D.; Li, X.; Wang, Z.; Wang, C. Monitoring 40-Year Lake Area Changes of the Qaidam Basin, Tibetan Plateau, Using Landsat Time Series. Remote Sens. 2019, 11, 343. [Google Scholar] [CrossRef] [Green Version]
  40. Nie, Y.; Sheng, Y.; Liu, Q.; Liu, L.; Liu, S.; Zhang, Y.; Song, C. A regional-scale assessment of Himalayan glacial lake changes using satellite observations from 1990 to 2015. Remote Sens. Environ. 2017, 189, 1–13. [Google Scholar] [CrossRef] [Green Version]
  41. Li, H.; Gao, Y.; Li, Y.; Yan, S.; Xu, Y. Dynamic of Dalinor Lakes in the Inner Mongolian Plateau and Its Driving Factors during 1976–2015. Water 2017, 9, 749. [Google Scholar] [CrossRef] [Green Version]
  42. Mueller, N.; Lewis, A.; Roberts, D.; Ring, S.; Melrose, R.; Sixsmith, J.; Lymburner, L.; McIntyre, A.; Tan, P.; Curnow, S.; et al. Water observations from space: Mapping surface water from 25 years of Landsat imagery across Australia. Remote Sens. Environ. 2016, 174, 341–352. [Google Scholar] [CrossRef] [Green Version]
  43. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  44. Yuan, F.; Sawaya, K.E.; Loeffelholz, B.C.; Bauer, M.E. Land cover classification and change analysis of the Twin Cities (Minnesota) Metropolitan Area by multitemporal Landsat remote sensing. Remote Sens. Environ. 2005, 98, 317–328. [Google Scholar] [CrossRef]
  45. Zhou, W.; Troy, A.; Grove, M. Object-based land cover classification and change analysis in the Baltimore metropolitan area using multitemporal high resolution remote sensing data. Sensors 2008, 8, 1613–1636. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Enders, C.K.; Mansolf, M. Assessing the fit of structural equation models with multiply imputed data. Psychol. Methods 2018, 23, 76. [Google Scholar] [CrossRef]
  47. Meydan, C.; Sesen, H. Structural Equation Modeling AMOS Applications; Detay Yayıncılık: Ankara, Turkey, 2011. [Google Scholar]
  48. Schreiber, J.B. Core reporting practices in structural equation modeling. Res. Soc. Adm. Pharm. 2008, 4, 83–97. [Google Scholar] [CrossRef]
  49. Melucci, M.; Paggiaro, A. Evaluation of information retrieval systems using structural equation modeling. Comput. Sci. Rev. 2019, 31, 1–18. [Google Scholar] [CrossRef]
  50. Gefen, D.; Straub, D.; Boudreau, M.-C. Structural equation modeling and regression: Guidelines for research practice. Commun. Assoc. Inf. Systems 2000, 4, 7. [Google Scholar] [CrossRef] [Green Version]
  51. Li, S.; Gang, A. A diagnostic study of northeast cold vortex heavy rain over the Songhuajiang-Nenjiang River Basin in the Summer of 1998. Chin. J. Atmos. Sci. 2001, 342–353. [Google Scholar]
  52. Zhang, F.; Tiyip, T.; Johnson, V.C.; Kung, H.-t.; Ding, J.-l.; Sun, Q.; Zhou, M.; Kelimu, A.; Nurmuhammat, I.; Chan, N.W. The influence of natural and human factors in the shrinking of the Ebinur Lake, Xinjiang, China, during the 1972–2013 period. Environ. Monit. Assess. 2014, 187, 4128. [Google Scholar] [CrossRef]
  53. Zhang, G.; Yao, T.; Shum, C.; Yi, S.; Yang, K.; Xie, H.; Feng, W.; Bolch, T.; Wang, L.; Behrangi, A. Lake volume and groundwater storage variations in Tibetan Plateau’s endorheic basin. Geophys. Res. Lett. 2017, 44, 5550–5560. [Google Scholar] [CrossRef]
  54. Zhong, L.; Liu, L.; Liu, Y. Natural Disaster Risk Assessment of Grain Production in Dongting Lake Area, China. Agric. Sci. Procedia 2010, 1, 24–32. [Google Scholar] [CrossRef] [Green Version]
  55. Mao, D.; He, X.; Wang, Z.; Tian, Y.; Xiang, H.; Yu, H.; Man, W.; Jia, M.; Ren, C.; Zheng, H. Diverse policies leading to contrasting impacts on land cover and ecosystem services in Northeast China. J. Clean. Prod. 2019, 240, 117961. [Google Scholar] [CrossRef]
  56. Mao, D.; Luo, L.; Wang, Z.; Wilson, M.C.; Zeng, Y.; Wu, B.; Wu, J. Conversions between natural wetlands and farmland in China: A multiscale geospatial analysis. Sci. Total Environ. 2018, 634, 550–560. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Geographic location and general situation of the Nenjiang watershed.
Figure 1. Geographic location and general situation of the Nenjiang watershed.
Water 12 01108 g001
Figure 2. Landsat images selection in optimum months for different dates.
Figure 2. Landsat images selection in optimum months for different dates.
Water 12 01108 g002
Figure 3. Flowchart for the lake and reservoir mapping process.
Figure 3. Flowchart for the lake and reservoir mapping process.
Water 12 01108 g003
Figure 4. Spatial distribution of lakes and reservoirs in the Nenjiang watershed (NJW) in 2015.
Figure 4. Spatial distribution of lakes and reservoirs in the Nenjiang watershed (NJW) in 2015.
Water 12 01108 g004
Figure 5. Decadal variations in the total area and number of lakes and reservoirs (area ≥ 1 km2) in the NJW from 1980 to 2015. The error bars represent the 95% confidence level.
Figure 5. Decadal variations in the total area and number of lakes and reservoirs (area ≥ 1 km2) in the NJW from 1980 to 2015. The error bars represent the 95% confidence level.
Water 12 01108 g005
Figure 6. Spatial variations of lake changes in the Nenjiang watershed during different periods. Red and blue dots present lake shrinkage and expansion, respectively, while the dot extent denotes change proportion in the corresponding period. (ad) represent the periods 1980–1990, 1990–2000, 2000–2010, and 2010–2015, respectively.
Figure 6. Spatial variations of lake changes in the Nenjiang watershed during different periods. Red and blue dots present lake shrinkage and expansion, respectively, while the dot extent denotes change proportion in the corresponding period. (ad) represent the periods 1980–1990, 1990–2000, 2000–2010, and 2010–2015, respectively.
Water 12 01108 g006
Figure 7. Spatial distribution of reservoir changes in the Nenjiang watershed during different periods. Purple and blue dots present reservoir shrinkage and expansion, respectively, while the dot extent denotes change proportion in the corresponding period. (ad) represent the periods 1980–1990, 1990–2000, 2000–2010, and 2010–2015, respectively.
Figure 7. Spatial distribution of reservoir changes in the Nenjiang watershed during different periods. Purple and blue dots present reservoir shrinkage and expansion, respectively, while the dot extent denotes change proportion in the corresponding period. (ad) represent the periods 1980–1990, 1990–2000, 2000–2010, and 2010–2015, respectively.
Water 12 01108 g007
Figure 8. Structural model and path coefficient. Arrows show the effect path and direction. Numbers adjacent to arrows are path coefficients, (β). The path coefficients (β) characterize the extent of the effect of the examined variable on the lake changes, while the larger value indicates a strong positive or negative effect. P represents a significant level. * represents p-value < 0.1, ** represents p-value < 0.05, and *** represents p-value < 0.01. MAP denotes the mean annual precipitation, MAAT represents the mean annual air temperature, ET represents the evapotranspiration, and GDP represents the gross domestic product. The “e” represents the error (residual term) of path analysis of observed variables in the structural model. The numbers on the arrows are the values of the standardized regression weights of the model.
Figure 8. Structural model and path coefficient. Arrows show the effect path and direction. Numbers adjacent to arrows are path coefficients, (β). The path coefficients (β) characterize the extent of the effect of the examined variable on the lake changes, while the larger value indicates a strong positive or negative effect. P represents a significant level. * represents p-value < 0.1, ** represents p-value < 0.05, and *** represents p-value < 0.01. MAP denotes the mean annual precipitation, MAAT represents the mean annual air temperature, ET represents the evapotranspiration, and GDP represents the gross domestic product. The “e” represents the error (residual term) of path analysis of observed variables in the structural model. The numbers on the arrows are the values of the standardized regression weights of the model.
Water 12 01108 g008
Figure 9. Past changes and future projections of climatic factors in the NJW. (a) Changes of mean annual air temperature (MAAT), mean annual precipitation (MAP), and evapotranspiration (ET) from 1980 to 2015 and (b) changes of future MAAT, MAP, and PET from 2015 to 2050.
Figure 9. Past changes and future projections of climatic factors in the NJW. (a) Changes of mean annual air temperature (MAAT), mean annual precipitation (MAP), and evapotranspiration (ET) from 1980 to 2015 and (b) changes of future MAAT, MAP, and PET from 2015 to 2050.
Water 12 01108 g009
Figure 10. Increased lake area due to floods in 1998 and 2013.
Figure 10. Increased lake area due to floods in 1998 and 2013.
Water 12 01108 g010
Figure 11. Image examples for lakes and reservoirs affected by the expanded paddy fields.
Figure 11. Image examples for lakes and reservoirs affected by the expanded paddy fields.
Water 12 01108 g011
Table 1. Collected validation samples for the five dates.
Table 1. Collected validation samples for the five dates.
Date19801990200020102015
SourcesGoogle Earth Image/MSSTMTM/ETM+TMOLI
Lakes126158170244256
Reservoirs3847506265
Non-water bodies3540796557
Total Samples199245299371378
MSS—multispectral scanner, TM—thematic mapper, ETM+—enhanced thematic mapper plus, OLI—operational land imager.
Table 2. Measures used to test the goodness of model.
Table 2. Measures used to test the goodness of model.
MeasureOptimum ValuesReference
RMSEA (root mean square error of approximation)Less than 0.08Li et al. (2019) [39]
λ2/df (chi-square/degree of freedom)Less than 3James (2007) [48]
GFI (goodness of fit index)0.90 and aboveMelucci et al. (2019) [49]
CFI (comparative fit index)0.90 and aboveDavid et al. (2000) [50]
Table 3. Lake number, area, and changes in different classes between 1980 and 2015.
Table 3. Lake number, area, and changes in different classes between 1980 and 2015.
Year or PeriodArea Classes (km2)Total
1–1010–5050–100>100
Number of lakes19804194333468
19904144013458
20003433904386
20101982723230
20151992923233
Change in number (%)1980–1990−1−7−670−2
1990–2000−17−3033−15
2000–2010−42−31100−25−41
2010–201507001
1980–2015−53−33−330−50
Lake area (km2)19801369 ± 11916 ± 17197 ± 20958 ± 243440 ± 72
19901145 ± 10884 ± 1581 ± 17900 ± 173010 ± 59
2000934 ± 8729 ± 190.0804 ± 132467 ± 40
2010548 ± 7535 ± 16111 ± 17767 ± 131961 ± 53
2015595 ± 7656 ± 14145 ± 15714 ± 172110 ± 53
Change in area (%)1980–1990−16 **−4 **−59 **−6.0 **−13 **
1990–2000−18 **−18 **−100 **−10.7 **−18 **
2000–2010−41 **−27 **0−4.5 **−20 **
2010–20159 **23 **30 **−7.0 **8 **
1980–2015−57 **−28 **−26 **−25.5 **−39 **
Note: * denotes change at a significant level of 0.05, ** denotes significant at 0.01.
Table 4. Reservoir number, area, and changes in different classes between 1980 and 2015.
Table 4. Reservoir number, area, and changes in different classes between 1980 and 2015.
Year or PeriodArea Classes (km2)Total
1–1010–5050–100>100
Number of reservoirs19806961278
199068152186
200073194096
2010841322101
20151081722129
Change in number (%)1980–1990−1150100−5010
1990–2000727100013
2000–201015−32−501004
2010–201529310028
1980–201557183100065
Reservoir area (km2)1980205 ± 9149 ± 1052 ± 14513 ± 20919 ± 53
1990200 ± 6328 ± 16121 ± 17151 ± 14800 ± 53
2000247 ± 5263 ± 8270 ± 150780 ± 28
2010320 ± 6324 ± 16134 ± 17308 ± 191086 ± 58
2015259 ± 4397 ± 16175 ± 18591 ± 181422 ± 56
Change in area (%)1980–1990−2121 **133 **−71 **−13 **
1990–200023 **−20 **123 **−100 **−3
2000–201030 **24 **−50 **039 **
2010–2015−19 **23 **31 **92 **31 **
1980–201527 **167 **237 **15 **55 **
Note: * denotes change at a significant level of 0.05, ** denotes significant at 0.01.

Share and Cite

MDPI and ACS Style

Du, B.; Wang, Z.; Mao, D.; Li, H.; Xiang, H. Tracking Lake and Reservoir Changes in the Nenjiang Watershed, Northeast China: Patterns, Trends, and Drivers. Water 2020, 12, 1108. https://doi.org/10.3390/w12041108

AMA Style

Du B, Wang Z, Mao D, Li H, Xiang H. Tracking Lake and Reservoir Changes in the Nenjiang Watershed, Northeast China: Patterns, Trends, and Drivers. Water. 2020; 12(4):1108. https://doi.org/10.3390/w12041108

Chicago/Turabian Style

Du, Baojia, Zongming Wang, Dehua Mao, Huiying Li, and Hengxing Xiang. 2020. "Tracking Lake and Reservoir Changes in the Nenjiang Watershed, Northeast China: Patterns, Trends, and Drivers" Water 12, no. 4: 1108. https://doi.org/10.3390/w12041108

APA Style

Du, B., Wang, Z., Mao, D., Li, H., & Xiang, H. (2020). Tracking Lake and Reservoir Changes in the Nenjiang Watershed, Northeast China: Patterns, Trends, and Drivers. Water, 12(4), 1108. https://doi.org/10.3390/w12041108

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