Unraveling the Role of Human Activities and Climate Variability in Water Level Changes in the Taihu Plain Using Artificial Neural Network

Water level, as a key indicator for the floodplain area, has been largely affected by the interplay of climate variability and human activities during the past few decades. Due to a nonlinear dependence of water level changes on these factors, a nonlinear model is needed to more realistically estimate their relative contribution. In this study, the attribution analysis of long-term water level changes was performed by incorporating multilayer perceptron (MLP) artificial neural network. We took the Taihu Plain in China as a case study where water level series (1954–2014) were divided into baseline (1954–1987) and evaluation (1988–2014) periods based on abrupt change detection. The results indicate that climate variables are the dominant driver for annual and seasonal water level changes during the evaluation period, with the best performance of the MLP model having precipitation, evaporation, and tide level as inputs. In the evaluation period, the contribution of human activities to water level changes in the 2000s is higher than that in the 1990s, which indicates that human activities, including the rapid urbanization, are playing an important role in recent years. The influence of human activities, especially engineering operations, on water level changes in the 2000s is more evident during the dry season (March-April-May (MAM) and December-January-February (DJF)).


Introduction
It is widely recognized that climate variability and human activities have exerted considerable impacts on hydrological processes globally during the past decades [1][2][3][4].Climate variability can influence these processes through altering the spatial and temporal patterns of precipitation and evaporation [5][6][7].The influence of human activities, such as landscape change, hydraulic construction, and river channelization, make hydrological processes more complex [8][9][10].In order to effectively manage regional water resources, it is important to investigate the influence of both climate variability and human activities on hydrological processes.
The water level is an important indicator of surface water studies and flood management across the floodplain areas [11,12].Recently, the variations in water level have attracted increasing attention because of a worldwide increased risk of flood and drought/water scarcity [13,14].Some researchers have investigated the influence of either climate variability or anthropogenic activities on water level variations [15,16].However, both climatic and anthropogenic factors should be considered for the attribution analyses of water level variations.As human activities and the interventions in catchments have increased in recent decades, especially in fast-developing countries, such as China, the need to recognize and address the related impact becomes even more urgent [17][18][19][20].Despite this urgency, the number of studies quantifying the relative contribution of climate variability and human activities to water level changes in China is still limited.Zhang et al. [21] reported that sand excavation was one of the most intensive human activities after the 1980s in the Pearl River Delta, which affected extreme water levels by altering the morphology of the river channels and estuary.In the Hang-Jia-Hu plain, rapid urbanization since the 1980s has changed the temporal and spatial variation of water levels [22].In the Mekong River basin, the temporal and spatial trends in water infrastructure development after 1991 have caused alterations in water level patterns [23].
The Taihu Plain is one of the typical Chinese floodplains located in the center of the Yangtze River Delta with complex fluvial systems and heavy industrialization [24].Due to the flat topography and gentle slope, the plain has a complex flow direction which makes runoff monitoring quite difficult.It, moreover, experiences continuous rainfall with high flood risks in summer [25].Thus, water level becomes a prime indicating factor for hydraulic construction as well as flood control.The Taihu Plain is also one of the most prosperous areas in China and has experienced a rapid urban development over the past three decades, which has triggered the occurrence of different kinds of human activities in this region [26].With further urbanization and a projected increase for future summer precipitation in this region, the flood risks are expected to increase by approximately 4-15 times by 2050 under different climatic and socio-economic scenarios [27].The largest flood disaster in the region in 1999 resulted in a direct economic loss of $16 billion.Therefore, it is necessary to conduct a quantitative attribution analysis of each driving force on regional hydrological processes to reduce the increasing flood risks and implement water management projects in the floodplain area.
The previous studies have mostly performed linear attribution analyses for hydrological changes [21,23].However, due to an interaction of different human activities and climate variability influencing water level, water level varies nonlinearly.Modeling of such complex process requires a robust model as otherwise, results can be biased.The artificial neural network (ANN) is one of the mathematical models with high abilities to find the nonlinear relationship between input and output parameters [28].Until now, although the ANN method has been widely used in many hydrological and meteorological applications [29][30][31][32], there is a lack of study using ANN for the attribution analysis of hydrological changes.To address this knowledge gap, this study applied multilayer perceptron (MLP) artificial neural network for the attribution analysis of long-term water level changes during 1954-2014 in the Taihu Plain.The relative contributions of climate variability and human activities to water level changes are quantified.

Study Area
The Taihu Plain covers a total land area of 7929 km 2 , with an elevation of 2-4 m above sea level (Figure 1).The region is controlled by the East Asia Monsoon with four distinct seasons: spring (March-April-May: MAM), summer (June-July-August: JJA), autumn (September-October-November: SON), and winter (December-January-February: DJF).The average annual rainfall is 1180 mm, and the average annual temperature is 16 • C. The precipitation is mainly concentrated in summer, which accounts for almost 60% of the annual precipitation [33,34].
Rapid urbanization and economic growth have taken place in this region since the 1980s [35].The population of the Taihu Plain increased from 9.08 million in 1962 to 14.65 million in 2010.Industrial production in 2010 increased by 1000 times than that in 1960 and reached 4500 billion RMB [36].The proportion of farmland areas decreased sharply, while industrial and urban areas increased by 4.7% annually during the same period, which indicates that the urbanization in this region is a process of farmland occupation.River system has also been severely disturbed during the rapid urbanization, with the decrease in drainage density and water surface ratio by 12.6% and 20.2% from the 1980s to the 2010s, respectively [24].The hydrological processes of the Taihu Plain are also affected by tide level changes to some extent [37,38].Meanwhile, more and more hydraulic structures, including sluices and pumps, have been constructed for water resources management and flood control [39].Overall, the dense population, rapid industrialization, and heavy urbanization have put considerable pressure on the hydrology of the region.affected by tide level changes to some extent [37,38].Meanwhile, more and more hydraulic structures, including sluices and pumps, have been constructed for water resources management and flood control [39].Overall, the dense population, rapid industrialization, and heavy urbanization have put considerable pressure on the hydrology of the region.

Data Preparation
Monthly water level data from January 1954 to December 2014 were collected from eight water level stations within the Taihu Plain (Figure 1).Monthly precipitation data of the same period were obtained from ten rainfall stations.Monthly tidal level data at Jiangyin station located in the estuary of the Yangtze River were also collected.What's more, monthly meteorological data of Dongshan station obtained from the China Administration of Meteorology, including air temperature, sunshine hours, relative humidity, and wind speed, were used to estimate potential evaporation using the Food and Agriculture Organization (FAO) Penman-Monteith method.
The geographic information of the study stations and some statistical metrics, including longterm mean (LTM) and coefficient of variation (CV), are listed in Table 1.Considering similar characteristics of water level and rainfall at each station, regional mean time series for the Taihu Plain was calculated by the Thiessen polygon in ArcGIS 9.3.In addition to climate data, human activity related data, including land use and hydraulic structure within the Taihu Plain, were collected.The land use data with a resolution of 1 km for five periods during 1990-2010 were provided by Data Center for Resources and Environmental Sciences (RESDC, http://www.resdc.cn).The hydraulic structures, including sluices and pumps, were acquired from the Taihu Basin Authority.

Data Preparation
Monthly water level data from January 1954 to December 2014 were collected from eight water level stations within the Taihu Plain (Figure 1).Monthly precipitation data of the same period were obtained from ten rainfall stations.Monthly tidal level data at Jiangyin station located in the estuary of the Yangtze River were also collected.What's more, monthly meteorological data of Dongshan station obtained from the China Administration of Meteorology, including air temperature, sunshine hours, relative humidity, and wind speed, were used to estimate potential evaporation using the Food and Agriculture Organization (FAO) Penman-Monteith method.
The geographic information of the study stations and some statistical metrics, including long-term mean (LTM) and coefficient of variation (CV), are listed in Table 1.Considering similar characteristics of water level and rainfall at each station, regional mean time series for the Taihu Plain was calculated by the Thiessen polygon in ArcGIS 9.3.In addition to climate data, human activity related data, including land use and hydraulic structure within the Taihu Plain, were collected.The land use data with a resolution of 1 km for five periods during 1990-2010 were provided by Data Center for Resources and Environmental Sciences (RESDC, http://www.resdc.cn).The hydraulic structures, including sluices and pumps, were acquired from the Taihu Basin Authority.

Statistical Analysis
To determine the change point in the water level time series, the Pettitt [40], Sequential Cluster [41], and Mann-Kendall [42,43] tests were used.They are commonly used methods to identify the occurrence of a change point in hydrological series [44][45][46].Then, the most probable change point is selected based on the overall results of the three methods.In this way, the detected change points by different methods can be validated.Once the change point is detected, the time series is divided into two sub-series and the series, which are expected to have different statistics (e.g., mean and standard deviation).
The Mann-Whitney test was used to test the null hypothesis that the two samples come from the same population; whether observations in one sample tend to be larger or smaller than observations in the other sample.Given the time series ).The Mann-Whitney test statistic is given as where r(x i ) is the rank of the observations, n1 and n2 are the sizes of samples.The null hypothesis H 0 of no steep change between two periods is accepted if , where α is the significance level.Z 1−α/2 for the selected α = 0.01 and α = 0.05 is equal to 2.58 and 1.96, respectively.

ANN-Based Attribution Analysis
In order to separate the impacts of climate variability and human activities on the changes in water level, the sub-series before and after the detected change points were referred to as baseline and evaluation periods, respectively.An assumption was made that there is no significant human activity during the baseline period  [47].The total change (∆H) in mean water level between two independent periods can be estimated as: where H obs 1 and H obs 2 denote the observed mean water level during the baseline and evaluation periods, respectively.The total changes in water level caused by the combined impacts of climate variability and human activities can be expressed as: (5) where ∆H c is the water level change due to the climate variability between two periods, ∆H h is the water level change caused by the changes in human activities, H sim 2 is natural water level during the evaluation period, and η c and η h are the contributions of climate and human activity to the water level change, respectively.To partition water level impacts induced by climate variability and human activities, ∆H c and ∆H h need to be quantified.In this study, a monthly scale non-linear model between mean monthly water and climate variables (i.e., monthly rainfall, evaporation, and tide level) was built based on ANN.The H sim 2 and ∆H c can be estimated using the validated ANN model.Then, the percentage contribution of climate-induced (η c ) and human-activities-induced (η h ) to water level changes can be calculated.The used ANN model is briefly explained below.In this study, the multilayer perceptron (MLP) model was applied, which is probably the most popular ANN model in hydrological studies [48][49][50].The MLP is a feed-forward network that consists of three layers of neurons: the input layer, the hidden layers, and the output layer.The hidden layer provides the nonlinearity between the input and output sets.Most complex problems can be solved by increasing the number of hidden layers and/or the neurons in the hidden layer.The MLP neural network can be expressed as follows.
where Y t is the output vector, x i is the input vector, w i and w j are the weights vector between the neurons of the input and hidden layers and between the hidden and output layers, respectively, f 1 and f 2 are the transfer functions for the hidden layer and output layer, respectively.In order to lie all the data in the central zone of the sigmoid function, all inputs were normalized between 0.1 and 0.9 before training the MLP model, as recommended by [51].The normalization prevents the problem of output signal saturation that sometimes happens in ANN applications [52].The model was trained using the Levenberg-Marquardt algorithm.
To find the optimal hidden neurons, a trial and error process was used, starting with two hidden neurons and gradually increased up to 10 with a step size of 1 in each trial.For each set of hidden neurons, the network was trained in a batch mode to minimize the mean square error at the output layer.In order to avoid overfitting, early stop technique and cross-validation were performed by keeping track of the efficiency of the fitted model [53].Once no significant improvement in the efficiency was noticed, the training was stopped.The best MLP model was selected based on Nash-Sutcliffe efficiency (NSE) and the coefficient of determination (R 2 ).

Variability of Annual and Seasonal Water Level
The variation of regional annual water level and rainfall from 1954 to 2014 was analyzed by linear regression.It can be seen from Figure 2a that an upward trend was detected for annual water level series.Minimum and maximum annual water level occurred in 1978 and 2010, respectively.Further visual inspection of Figure 2a shows the 1970s as a period with a lower water level than other decades, which is related to the shortage of rainfall in this decade.There is an increase in water level during the 1990s, which is related to the abundant rainfall in this decade.On the seasonal scale (Figure 2b-e), there are increasing linear trends in water level for each season.The evolution of water level in JJA similar to that of the annual water level.
The change point detection results show a significant change point at the 95% confidence level in annual water level at around 1988, 1990, and 1988 using the Pettitt, Mann-Kendall, and Sequential Cluster tests, respectively (Table 2).As for seasonal water level, significant abrupt changes are identified in the late 1980s.The change point identified for annual and seasonal water level in the late 1980s is more or less consistent with the beginning of the extensive rapid urbanization in the Taihu Plain.H 0 represents null hypothesis.R represents that H 0 is rejected at the 95% confidence level.MAM: March-April-May, JJA: June-July-August, SON: September-October-November, DJF: December-January-February.According to the results, it can be concluded that the Taihu Plain experienced two different hydrological periods with a change point at around the late 1980s.For comparison, the entire study period was split into two sub-periods (baseline and evaluation periods) by the change point year.The Mann-Whitney test was applied for comparing the water level in the two sub-periods (Table 3).The results show that the null hypothesis H 0 of the Mann-Whitney test is rejected at the 95% confidence level, implying significant changes in the annual and seasonal water level series between the two sub-periods.The largest shift changes in water level are observed for winter and spring, respectively.To better understand the characteristics of the water level changes, the means and standard deviations of annual and seasonal water level in the two sub-periods were compared.There is an increasing ratio of 9.9% for the mean annual water level in the evaluation period relative to the baseline period.For the seasonal scale, the increasing ratios of water level between the two sub-periods range from 7.5% to 13.3%.In contrast to water level, annual and seasonal rainfall series do not show a significant increase (Figure 2).Thus, we can infer that the significant changes in water level in the Taihu Plain might be influenced by a combination of human activity and climate variability.

MLP Model Setup
To quantify the relative contribution of climate variability and human activity to water level changes, a monthly scale MLP model was established based on the data in the baseline period.In order to obtain the optimal model, the MLP models with three input combinations were established to estimate monthly water levels (Table 4).Input variables include mean monthly rainfall, evaporation, and tide level.The 70% and 30% of the data in the baseline period were chosen as the training and validation periods of the MLP models, respectively.The optimal number of neurons in the hidden layer is six for all the MLP models.It is evident that the best performed MLP model is the one with all the three variables as input (i.e., MLP3).The MLP3 simulated water level values in both training and validation periods have generally a good consistency with observations.The R 2 , NSE, and relative error (RE) values are 0.85, 0.84, and 0.026 for the MLP3 for the training period and 0.81, 0.80, and 0.029 for the validation period, respectively (Figure 3).Hence, the validated MLP3 model was selected for the subsequent contribution analysis.

Contribution Analysis to Water Level Changes
Based on the validated MLP3 model, the climate-induced water level series in the evaluation period was simulated.The simulated and observed water level in the evaluation period is exhibited in Figure 4. Comparing the two-time series, there is a high consistency for the peak values, while the low water level is underestimated by the MLP3 model.The difference in low water level during the 2000s is significantly larger than that during the 1990s.The contribution of climate variability and human activities to water level changes is presented in Table 5. Regarding the evaluation period, the contributions of climate variability and human activities to annual water level changes are 69% and 31%, respectively.For each season, the contribution of climate variability is also higher than that of human activities.The contribution of climate variability is beyond 80% for summer and autumn.Thus, it indicates that climate variables play dominant roles in water level variations during the whole evaluation period.To further investigate the potential influence of climate variability and human activities on water level changes at a decadal scale, the evaluation period was divided into two shorter periods: sub-period I and II referring to the periods before and after 2000, respectively.The results of the attribution analysis become slightly different when we focus on two shorter periods.As shown in Table 5, the contribution of climate variability during the sub-period I is higher than that during the sub-period II for annual and seasonal water level.In other words, the contribution of human activities in the sub-period II becomes larger than that in the sub-period I.As for annual water level, the contribution of human activities (49%) is almost equal to that of climate variability (51%) during the sub-period II.As for each season in the sub-period II, the contribution of human activities during MAM (73%) and DJF (53%) is even beyond that of climate variability.Similar results also can be found in Figure 4 that the observed water level in MAM and DJF is larger than the simulated water level, especially after 2000.Overall, although climate variability plays dominant roles in water level variations, human activities also have an increasing impact in recent years, especially for MAM and DJF.
To further investigate the potential influence of climate variability and human activities on water level changes at a decadal scale, the evaluation period was divided into two shorter periods: subperiod I and II referring to the periods before and after 2000, respectively.The results of the attribution analysis become slightly different when we focus on two shorter periods.As shown in Table 5, the contribution of climate variability during the sub-period I is higher than that during the sub-period II for annual and seasonal water level.In other words, the contribution of human activities in the sub-period II becomes larger than that in the sub-period I.As for annual water level, the contribution of human activities (49%) is almost equal to that of climate variability (51%) during the sub-period II.As for each season in the sub-period II, the contribution of human activities during MAM (73%) and DJF (53%) is even beyond that of climate variability.Similar results also can be found in Figure 4 that the observed water level in MAM and DJF is larger than the simulated water level, especially after 2000.Overall, although climate variability plays dominant roles in water level variations, human activities also have an increasing impact in recent years, especially for MAM and DJF.

Impacts of Human Activities on Water Level Changes
In this study, three popular methods were employed to identify change points in annual and seasonal water level across the Taihu Plain.It is clear that abrupt changes of water level occurred in the late 1980s, which are consistent with previous studies using the Mann-Kendall test by Yin et al. [54] and Xu et al. [22].Although previous studies showed that the increasing anthropogenic perturbations have significant effects on water level alteration since the 1980s in the Taihu Plain, few studies have so far investigated their relative quantitative contribution.In our study, a validated MLP-based model was used to investigate the individual impact of climatic and anthropogenic factors on water level changes during the evaluation period .The evaluation period was also divided into two shorter periods (1990s and 2010s) in order to investigate the contribution of recent human activities to water level variation.
Due to urban development, intensive human activities have taken place within the Taihu Plain since the late 1980s, including hydraulic construction, land use change, and river reduction [9,35,55].For the sake of managing water resources and flood, a large number of hydraulic structures (pump, sluices, and dikes) have been constructed along the rivers after devastating floods in 1991 and 1999 [56][57][58].The increasing trend of hydraulic structure construction since the 1980s can be seen in Figure

Impacts of Human Activities on Water Level Changes
In this study, three popular methods were employed to identify change points in annual and seasonal water level across the Taihu Plain.It is clear that abrupt changes of water level occurred in the late 1980s, which are consistent with previous studies using the Mann-Kendall test by Yin et al. [54] and Xu et al. [22].Although previous studies showed that the increasing anthropogenic perturbations have significant effects on water level alteration since the 1980s in the Taihu Plain, few studies have so far investigated their relative quantitative contribution.In our study, a validated MLP-based model was used to investigate the individual impact of climatic and anthropogenic factors on water level changes during the evaluation period (1988-2014).The evaluation period was also divided into two shorter periods (1990s and 2010s) in order to investigate the contribution of recent human activities to water level variation.
Due to urban development, intensive human activities have taken place within the Taihu Plain since the late 1980s, including hydraulic construction, land use change, and river reduction [9,35,55].For the sake of managing water resources and flood, a large number of hydraulic structures (pump, sluices, and dikes) have been constructed along the rivers after devastating floods in 1991 and 1999 [56][57][58].The increasing trend of hydraulic structure construction since the 1980s can be seen in Figure 5.The total number of hydraulic structures has reached 2468 in the 2000s.The capacity of water abstraction for water management across the entire region has been enhanced significantly after 1989, which is almost consistent with the identified change point in the water level series.The alteration of the hydrological regime in the Taihu Plain due to these hydraulic structures has been shown by several studies [25,39].The water abstraction from the Yangtze River to the Taihu Plain through pumping stations has been mainly concentrated during the dry season (from September to April) since 2000 [59,60].Besides, it is worthy of noting that many polders are distributed across the Taihu Plain, and most of them are set up with high-power pumping stations and sluices since the devastating flood in 1991 [39].A large number of abstracted water from the Yangtze River is first stored in the river network of the Taihu Plain and then pumped into polders according to their individual demand in the dry season.Although a part of water resources is abstracted into polders for agricultural irrigation and water quality improvement, a high water level is kept in the river network.Thus, we can conclude that engineering operation plays an important role in water level changes during the dry season.The same results are confirmed by our contribution analysis showing a larger contribution of human activities during MAM and DJF.
Water 2019, 11, x FOR PEER REVIEW 11 of 16 abstraction for water management across the entire region has been enhanced significantly after 1989, which is almost consistent with the identified change point in the water level series.The alteration of the hydrological regime in the Taihu Plain due to these hydraulic structures has been shown by several studies [25,39].The water abstraction from the Yangtze River to the Taihu Plain through pumping stations has been mainly concentrated during the dry season (from September to April) since 2000 [59,60].Besides, it is worthy of noting that many polders are distributed across the Taihu Plain, and most of them are set up with high-power pumping stations and sluices since the devastating flood in 1991 [39].A large number of abstracted water from the Yangtze River is first stored in the river network of the Taihu Plain and then pumped into polders according to their individual demand in the dry season.Although a part of water resources is abstracted into polders for agricultural irrigation and water quality improvement, a high water level is kept in the river network.Thus, we can conclude that engineering operation plays an important role in water level changes during the dry season.The same results are confirmed by our contribution analysis showing a larger contribution of human activities during MAM and DJF.In addition to the hydraulic structure construction, the land use of the Taihu Plain has also altered dramatically due to the rapid urbanization.The general trend of land use in this region was characterized mainly by a reduced cultivated land and increased urbanized land, causing a significant increase in the effective impervious area and thereby a larger runoff coefficient [35].As shown in Figure 6, the proportion of the urbanized land increased by 1751 km 2 with an annual growth rate of 56.48 km 2 per year, while that of paddy land decreased from 74% in 1990 to 42% in 2010.It should be noted that the increase rate of urbanized land from 2000 to 2010 is much higher than that from 1990 to 2000, which might have accelerated the rainfall-runoff process in the 2000s compared with the 1990s [61,62].In addition, the reduction of river and lake also has an influence on hydrological regimes to some extent.Urbanization has accelerated the disappearance of lakes and river networks, which reduces river storage capacity [63][64][65].The river density and water surface ratio decreased by 17% and 19% over the past 50 years, respectively [24].The reduction of rivers and lakes could lead to excess runoff, exceeding the present channel capacity, and increase the possibility of flood disasters.In addition to the hydraulic structure construction, the land use of the Taihu Plain has also altered dramatically due to the rapid urbanization.The general trend of land use in this region was characterized mainly by a reduced cultivated land and increased urbanized land, causing a significant increase in the effective impervious area and thereby a larger runoff coefficient [35].As shown in Figure 6, the proportion of the urbanized land increased by 1751 km 2 with an annual growth rate of 56.48 km 2 per year, while that of paddy land decreased from 74% in 1990 to 42% in 2010.It should be noted that the increase rate of urbanized land from 2000 to 2010 is much higher than that from 1990 to 2000, which might have accelerated the rainfall-runoff process in the 2000s compared with the 1990s [61,62].In addition, the reduction of river and lake also has an influence on hydrological regimes to some extent.Urbanization has accelerated the disappearance of lakes and river networks, which reduces river storage capacity [63][64][65].The river density and water surface ratio decreased by 17% and 19% over the past 50 years, respectively [24].The reduction of rivers and lakes could lead to excess runoff, exceeding the present channel capacity, and increase the possibility of flood disasters.
Overall, the impact of these human activities on the water level changes is coexistent and complicated.However, it can be inferred from the results that the construction of hydraulic structures is the most important driving factor in recent years, which can cause dramatic changes in water level.Especially, land use change and river reduction are gradual processes whose impacts on water level can accumulate gradually.Therefore, scientific considerations of engineering operations, land development, and river protection are very necessary for future basin management practices.Overall, the impact of these human activities on the water level changes is coexistent and complicated.However, it can be inferred from the results that the construction of hydraulic structures is the most important driving factor in recent years, which can cause dramatic changes in water level.Especially, land use change and river reduction are gradual processes whose impacts on water level can accumulate gradually.Therefore, scientific considerations of engineering operations, land development, and river protection are very necessary for future basin management practices.

Limitations and Recommendations for Future Research
In this study, the sub-series before the identified change point was used as the baseline period to analyze the relative impacts of climate variability and human activity in the evaluation period, which is because the intensity of human activity is relatively small in the baseline period.Some uncertainties exist in the attribution analysis when separating the impacts of climate variability and human activity on water level changes.Firstly, due to limited long-term observation data in the Taihu Plain, the performance of attribution analysis depends on the water level of the baseline period, with no of human activity on the model's calibration.In reality, some human disturbances, such as agricultural irrigation, land use changes, and river regulations, still existed during the baseline period.Hence, the contribution of human activities can be slightly overestimated, although the calibrated ANN model well reflects the water level evolution during the baseline period.Secondly, the employed attribution analysis method assumes that climate and human activity are independent of each other, while climate and human activity interact with each other especially in a catchment scale [66].In addition, the ANN-based attribution analysis highly depends on the quality of rainfall, tide, and evaporation measurements and estimates.The tide (evaporation) data from Jiangyin (Dongshan) station in the study area might not be representative of the entire region.The possible influences of the assumptions and data error on the attribution results are yet to be investigated in future research.

Conclusions
In this study, the relative impacts of climate variability and human activities on water level changes across the Taihu Plain were investigated based on statistical analysis and modeling approach.Several popular methods were used to detect change points in annual and seasonal water level time series to determine the baseline period, and then a monthly scale MLP model was established and used for the attribution analysis.
Significant upward trends were detected for annual and seasonal water level.The abrupt change of water level mainly occurred in the late 1980s, which corresponds to the beginning of intensive human activities in this region.For the whole evaluation period, climate variability is the dominant driver for water level changes, accounting for a 69% contribution to annual water level changes

Limitations and Recommendations for Future Research
In this study, the sub-series before the identified change point was used as the baseline period to analyze the relative impacts of climate variability and human activity in the evaluation period, which is because the intensity of human activity is relatively small in the baseline period.Some uncertainties exist in the attribution analysis when separating the impacts of climate variability and human activity on water level changes.Firstly, due to limited long-term observation data in the Taihu Plain, the performance of attribution analysis depends on the water level of the baseline period, with no effect of human activity on the model's calibration.In reality, some human disturbances, such as agricultural irrigation, land use changes, and river regulations, still existed during the baseline period.Hence, the contribution of human activities can be slightly overestimated, although the calibrated ANN model well reflects the water level evolution during the baseline period.Secondly, the employed attribution analysis method assumes that climate and human activity are independent of each other, while climate and human activity interact with each other especially in a catchment scale [66].In addition, the ANN-based attribution analysis highly depends on the quality of rainfall, tide, and evaporation measurements and estimates.The tide (evaporation) data from Jiangyin (Dongshan) station in the study area might not be representative of the entire region.The possible influences of the assumptions and data error on the attribution results are yet to be investigated in future research.

Conclusions
In this study, the relative impacts of climate variability and human activities on water level changes across the Taihu Plain were investigated based on statistical analysis and modeling approach.Several popular methods were used to detect change points in annual and seasonal water level time series to determine the baseline period, and then a monthly scale MLP model was established and used for the attribution analysis.
Significant upward trends were detected for annual and seasonal water level.The abrupt change of water level mainly occurred in the late 1980s, which corresponds to the beginning of intensive human activities in this region.For the whole evaluation period, climate variability is the dominant driver for water level changes, accounting for a 69% contribution to annual water level changes versus a 31% contribution of human activities.However, further attribution analysis for two sub-periods of the evaluation period reveals that the contribution of human activities to water level changes in the 2000s is higher than that in the 1990s.Especially for the dry season (MAM and DJF), human activities have become the dominant factors during the 2000s.
The present study focused on evaluating individual impacts of climatic and anthropogenic factors on long-term water level variations.Compared with previous studies, this study provided some new insights on the impact of human activities on annual and seasonal water level variations and proposed the use of nonlinear ANN models for attribution analyses.Although different kinds of human activities are interconnected across the Taihu Plain, yet abrupt changes in water level are basically synchronized with the time of hydraulic construction, with more noticeable effects of engineering operations during the dry season (MAM and DJF).Moreover, the individual impact of each driving force could be more informative for designing countermeasures compared to the combined impacts.Such information is helpful for water resources managers for risk mitigation planning for flood and drought disasters.Additionally, considering the impacts of multiple human activities on the water level, further studies are needed to analyze the contribution of individual anthropogenic activities on water level changes based on coupling hydrological and hydraulic models.Such studies will be beneficial to identify the most important driver for water level changes in order to achieve sustainable watershed management under climate change and extensive human activities.

Figure 1 .
Figure 1.Location of the study area and the considered gauge stations along with the digital elevation model.

Figure 1 .
Figure 1.Location of the study area and the considered gauge stations along with the digital elevation model.

Figure 2 .
Figure 2. Evolution of regional annual and seasonal water level and rainfall in the Taihu Plain from 1954 to 2014.Trendlines were fitted separately to the time series before and after the identified change points (baseline and evaluation periods).

Figure 2 .
Figure 2. Evolution of regional annual and seasonal water level and rainfall in the Taihu Plain from 1954 to 2014.Trendlines were fitted separately to the time series before and after the identified change points (baseline and evaluation periods).

Figure 3 .
Figure 3.Comparison of the observed and the MLP3 simulated monthly water level for training and validation phases for the baseline period (1954-1987).MLP: multilayer perceptron.

Figure 3 .
Figure 3.Comparison of the observed and the MLP3 simulated monthly water level for training and validation phases for the baseline period (1954-1987).MLP: multilayer perceptron.

Figure 4 .
Figure 4. Comparison between the observed and the MLP3 simulated monthly water level during the evaluation period (1988-2014).MLP: multilayer perceptron.

Figure 4 .
Figure 4. Comparison between the observed and the MLP3 simulated monthly water level during the evaluation period (1988-2014).MLP: multilayer perceptron.

Figure 5 .
Figure 5. (a) Changes in hydraulic construction and (b) water abstraction capacity in the Taihu Plain.

Figure 5 .
Figure 5. (a) Changes in hydraulic construction and (b) water abstraction capacity in the Taihu Plain.

Figure 6 .
Figure 6.Land use changes across the Taihu Plain from 1990 to 2010.

Figure 6 .
Figure 6.Land use changes across the Taihu Plain from 1990 to 2010.

Table 1 .
Selected stations and their basic properties of annual data.

Station Name Measuring Variable Longitude ( • E) Latitude ( • N) Mean Value 1 CV Period of Record
1Rainfall and evaporation are in mm/year, and water and tide levels in m.CV: coefficient of variation.

Table 2 .
Results of change point detection using the Pettit, Mann-Kendall, and Sequential Cluster tests.

Table 3 .
Shift change test using the Mann-Whitney test, and statistics summary before and after the identified change point for annual and seasonal water level.indicates significant changes at the 95% level.Z c represent Mann-Whitney test statistic.CV and SD represent the coefficient of variation and standard deviation, respectively. *

Table 3 .
Shift change test using the Mann-Whitney test, and statistics summary before and after the identified change point for annual and seasonal water level.indicates significant changes at the 95% level.Zc represent Mann-Whitney test statistic.CV and SD represent the coefficient of variation and standard deviation, respectively. *

Table 4 .
Best MLP architecture and performance criteria for proposed models.

Table 5 .
Contributions of climate variability and human activities to water level changes in the Taihu Plain.
R: monthly rainfall; E: evaporation; T: tide; W: water level.η c and η h are the contributions of climate and human activity to the water level change, respectively.Baseline and evaluation periods cover 1954-1987 and 1988-2014, respectively.Evaluation sub-period I and II cover the 1990s (1988-1999) and the 2000s (2000-2014), respectively.MAM: March-April-May; JJA: June-July-August; SON: September-October-November; DJF: December-January-February.