Quantifying the Individual Contributions of Climate Change, Dam Construction, and Land Use/Land Cover Change to Hydrological Drought in a Marshy River

Hydrological drought for marshy rivers is poorly characterized and understood. Our inability to quantify hydrological drought in marshy river environments stems from the lack of understanding how wetland loss in a river basin could potentially change watershed structure, attenuation, storage, and flow characteristics. In this study, hydrological drought in a marshy river in far Northeast China at a higher latitude was assessed with a streamflow drought index (SDI). A deterministic, lumped, and conceptual Rainfall–Runoff model, the NAM (Nedbor Afstromnings Model), was used to quantify the individual contributions of climate change, land use/land cover (LULC) change, and river engineering to hydrological drought. We found that in the last five decades, the frequency of hydrological droughts has been 55% without considering LULC change and reservoir construction in this wetland-abundant area. The frequency of hydrological drought increased by 8% due to land use change and by 19% when considering both the impacts of LULC change and a reservoir construction (the Longtouqiao Reservoir). In addition to the more frequent occurrence of hydrological droughts, human activities have also increased drought intensity. These findings suggest that LULC and precipitation changes play a key role in hydrological drought, and that the effect can be significantly modified by a river dam construction.


Introduction
Hydrological drought is one of the four major types of drought occurrences, alongside meteorological drought, agricultural drought, and socio-economic drought [1,2]. Hydrological drought is a natural hazard defined as a significant deficiency of water availability in a particular period and over a particular area [3,4], caused by multifaceted interactions between climate and catchment processes [5]. Thus, hydrological droughts can have direct impacts on a number of stakeholders, such as irrigation, power generation, and recreational purposes within the affected river basin [6,7]. Understanding hydrological drought and its relationship with climate change, land use change, and river engineering can be helpful for sustainable surface and groundwater water resource management [8].
Several recent studies [9][10][11][12] have found that hydrological droughts have sharply increased in both frequency and intensity in many parts of the world due to climate change and/or due to increased water demand caused by population growth and rapid expansion of the agricultural, energy, and industrial sectors. In the most recent decade, increasing attention has been devoted to identifying site-specific hydrological droughts in different basin environments. For instance, Lopes The headwater area of the NR, upstream of the Baoqing hydrological station (Figure 1), was chosen for this study. The headwater area has a total drainage area of approximately 3767 km 2 , providing crucial water resources for the downstream wetlands. In 2002, a 750 m long, 25.7 m tall dam across the river was completed, creating a reservoir known as the Longtouqiao Reservoir, with a maximum storage capacity of 615 million m 3 . The reservoir receives discharge from the uppermost 1730 km 2 of the NRB (Figure 1), causing an immediate streamflow change downstream [24,25].

Change Point Detection and Division of Sub-Periods
In order to separate the individual contributions of climate change and human activities to streamflow changes that occurred from 1961 to 2014, we made a time separation of the 54-year study period by testing a changing point of annual streamflow in the river basin using the Mann-Kendall method [26][27][28]. The significance of a change was verified at the α = 0.05 level, with rejection of the null hypothesis of no trend if |Z| > 1.96 [29,30]. Additionally, the year of 2002 was also considered separately in the time separation because of the Longtouqiao Reservoir construction. Based on the change point detection and the reservoir building, we divided the 54-year study period into three sub-periods: (1) 1961-1966, when the area had little in the way human activities; (2) from 1967 to 2001, during which time-intensive LULC change occurred; and (3) after 2002, when the dam construction was completed.

Impacts of Climate Change and Human Activity on Streamflow
In this study, we used a slope change ratio of cumulative quantity (SCRCQ) to quantify the contribution of climate change and human activities to the streamflow change [31,32]. Climate change was represented by precipitation and potential evapotranspiration (PET) changes, while human  The headwater area of the NR, upstream of the Baoqing hydrological station (Figure 1), was chosen for this study. The headwater area has a total drainage area of approximately 3767 km 2 , providing crucial water resources for the downstream wetlands. In 2002, a 750 m long, 25.7 m tall dam across the river was completed, creating a reservoir known as the Longtouqiao Reservoir, with a maximum storage capacity of 615 million m 3 . The reservoir receives discharge from the uppermost 1730 km 2 of the NRB (Figure 1), causing an immediate streamflow change downstream [24,25].

Change Point Detection and Division of Sub-Periods
In order to separate the individual contributions of climate change and human activities to streamflow changes that occurred from 1961 to 2014, we made a time separation of the 54-year study period by testing a changing point of annual streamflow in the river basin using the Mann-Kendall method [26][27][28]. The significance of a change was verified at the α = 0.05 level, with rejection of the null hypothesis of no trend if |Z| > 1.96 [29,30]. Additionally, the year of 2002 was also considered separately in the time separation because of the Longtouqiao Reservoir construction. Based on the change point detection and the reservoir building, we divided the 54-year study period into three sub-periods: (1) 1961-1966, when the area had little in the way human activities; (2) from 1967 to 2001, during which time-intensive LULC change occurred; and (3) after 2002, when the dam construction was completed.

Impacts of Climate Change and Human Activity on Streamflow
In this study, we used a slope change ratio of cumulative quantity (SCRCQ) to quantify the contribution of climate change and human activities to the streamflow change [31,32]. Climate change Sustainability 2020, 12, 3777 4 of 16 was represented by precipitation and potential evapotranspiration (PET) changes, while human activities were represented by land use/land cover change and the Longtouqiao Reservoir construction. The slope of the linear relationship between the year and cumulative streamflow before and after the changing year was assumed to be S Sb and S Sa (10 8 m 3 /year), respectively. The slope of the linear relationship between the year and cumulative precipitation before and after the changing year was S Pb and S pa (mm), respectively. The streamflow variation ratio and precipitation variation ratio were expressed as (S Sa − S Sb )/|S Sb | and (S pa − S pb )/|S Pb |, respectively. The contribution of precipitation (C P , unit: %) to streamflow variation after the changing year was quantified as follows: Similarly, the slope of the linear relationship between the year and cumulative PET before and after the changing year was assumed to be S Eb and S Ea (mm), respectively. The variation ratio of PET can be expressed as ((S Ea − S Eb )/|S Eb |). Therefore, the contribution of PET (C E , unit: %) to streamflow variation after the changing year was determined by: Based on the water balance, the contribution of human activities (C H , unit: %) to the streamflow variation was expressed as: Based on Equations (1) and (2), the individual contributions of LULC (C LULC , unit: %) and the Longtouqiao Reservoir (C Res , unit: %) to streamflow variation were estimated by:

Hydrological Model
(1) Model Introduction In this study, we used a deterministic, lumped, and conceptual Rainfall-Runoff model, the NAM (Nedbor Afstromnings Model), to simulate streamflow for the headwater area of the NR. The model is part of the Mike Basin software package developed by the Danish Hydraulic Institute (DHI) in Denmark, and has been widely used for different geographical regions in the world [33][34][35][36][37][38]. The model uses time series of precipitation and evapotranspiration as input, and estimates the rainfall-runoff process using the linkage rule between the four different storages, which are connected together; each is the representative of different physical specifications ( Figure 2). These four storages are snow storage, surface water storage, root zone storage, and groundwater storage. However, in this study, snow storage was not used due to a lack of corresponding monitoring data. The NAM was prepared with nine parameters (Table 2) representing the rest of the storage types. More details on the NAM can be found in the "MIKE Hydro" Reference Guide [39].  (2) Data Input The basic input data for the NAM model included the catchment boundary conditions and meteorological data, such as precipitation and potential evapotranspiration. In this study, we collected (1) DEM (Digital Elevation Model) data for this region from the Chinese Geospatial Data Cloud (http://www.gscloud.cn) to set up spatial boundary parameters, (2) daily precipitation, as well as daily air temperature, humidity, and radiation data from the Baoqing weather station for estimation of potential evapotranspiration, (3) discharge records from the Baoqing hydrological station for model calibration and validation, and (4) LULC data, which was collected by remote sensing interpretation, as well as from previous studies [20,[23][24][25]36,37].
(3) Calibration and Validation of the Model We calibrated and validated the NAM model for two periods: 1967-2014 and 2002-2014, that is, prior to and after the Longtouqiao Reservoir construction, in order to quantify the river engineering effect on streamflow. The model performance for the calibration and validation was assessed using the Nash-Sutcliffe efficiency (ENS) [39], as follows in Table 3. Table 3. Illustration of Nash-Sutcliffe efficiency (ENS).   (2) Data Input The basic input data for the NAM model included the catchment boundary conditions and meteorological data, such as precipitation and potential evapotranspiration. In this study, we collected (1) DEM (Digital Elevation Model) data for this region from the Chinese Geospatial Data Cloud (http://www.gscloud.cn) to set up spatial boundary parameters, (2) daily precipitation, as well as daily air temperature, humidity, and radiation data from the Baoqing weather station for estimation of potential evapotranspiration, (3) discharge records from the Baoqing hydrological station for model calibration and validation, and (4) LULC data, which was collected by remote sensing interpretation, as well as from previous studies [20,[23][24][25]36,37].

Nash-Sutcliffe Efficiency
(3) Calibration and Validation of the Model We calibrated and validated the NAM model for two periods: 1967-2014 and 2002-2014, that is, prior to and after the Longtouqiao Reservoir construction, in order to quantify the river engineering effect on streamflow. The model performance for the calibration and validation was assessed using the Nash-Sutcliffe efficiency (E NS ) [39], as follows in Table 3.   (Table 4). The natural streamflow without the impacts of human activities during 1967-2014 (Q S1 (2002-2014) ) was simulated with the calibrated parameters (Table 5) (Table 5). Note: " # " means the calibration, "*"means the validation.    The validation period (1965)(1966) in the headwater area of the Naoli River in Northeast China. Note: " # " means the calibration, "*"means the validation. The NAM was first calibrated with the meteorological and streamflow data from 1961-1964 (Figure 3a), and the observed streamflow in 1965-1966 ( Figure 3b) was used for validating the model. The results during the periods of calibration and validation were assessed using ENS (Table 4). The natural streamflow without the impacts of human activities during 1967-2014 (QS1 (2002-2014)) was simulated with the calibrated parameters (Table 5) (Table 5). Note: " # " means the calibration, "*"means the validation.    Comparing the simulated streamflow using NAM with the observed streamflow, it was found that the cumulative streamflow, impacted by LULC change and the Longtouqiao Reservoir, has obviously declined ( Figure 5). Comparing the simulated streamflow using NAM with the observed streamflow, it was found that the cumulative streamflow, impacted by LULC change and the Longtouqiao Reservoir, has obviously declined ( Figure 5).

Analysis of Hydrological Drought
A number of drought indices have been presented over the past decades to distinguish and quantify drought events [40]. These include the Surface Water Supply Index [41], the Standardized Precipitation Index (SPI) [42], the Standardized Precipitation Evapotranspiration Index (SPEI) [43], the Reconnaissance Drought Index (RDI) [44], the Rainfall Anomaly Index (RAI) [45], the Standardized Palmer Drought Index (SPDI) [46], and the Streamflow Drought Index (SDI) [47]. Information on hydrological drought indicators is important for the design and operation of many sectors that serve society and the economy, including public water supply, irrigation, energy, navigation, and industry [1]. As such, there is an urgent need to advance hydrological drought research and operational applications. In this study, the hydrological drought was assessed by using the Streamflow Drought Index (SDI) developed by Nalbantis and Tsakiris [47], because this index (SDI) can fully represent the conditions of hydrological drought in the river.
We calculated SDI series (SDI-12) with a running series of monthly observed and simulated streamflows in a 12-month interval. Application of this approach was mainly due to the fact that a mild drought or a dependent drought can be removed or combined into larger independent drought events at such a time scale [48]. Moreover, there have been studies showing that the ability of the many drought indices, using a 12-month time interval, to recognize an anomaly in hydrological variables is much stronger than on a short time scale [9,49]. By comparing the observed SDI-12 with the simulated SDI-12, the impacts of climate change and human activities on the hydrological drought were identified. Based on the SDI, hydrological droughts can be divided into five levels ( Table 6).
In this study, the frequency of hydrological droughts was also used for quantifying the impacts of climate change and human activities on hydrological droughts. This can be calculated using the following formulas:

Analysis of Hydrological Drought
A number of drought indices have been presented over the past decades to distinguish and quantify drought events [40]. These include the Surface Water Supply Index [41], the Standardized Precipitation Index (SPI) [42], the Standardized Precipitation Evapotranspiration Index (SPEI) [43], the Reconnaissance Drought Index (RDI) [44], the Rainfall Anomaly Index (RAI) [45], the Standardized Palmer Drought Index (SPDI) [46], and the Streamflow Drought Index (SDI) [47]. Information on hydrological drought indicators is important for the design and operation of many sectors that serve society and the economy, including public water supply, irrigation, energy, navigation, and industry [1]. As such, there is an urgent need to advance hydrological drought research and operational applications. In this study, the hydrological drought was assessed by using the Streamflow Drought Index (SDI) developed by Nalbantis and Tsakiris [47], because this index (SDI) can fully represent the conditions of hydrological drought in the river.
We calculated SDI series (SDI-12) with a running series of monthly observed and simulated streamflows in a 12-month interval. Application of this approach was mainly due to the fact that a mild drought or a dependent drought can be removed or combined into larger independent drought events at such a time scale [48]. Moreover, there have been studies showing that the ability of the many drought indices, using a 12-month time interval, to recognize an anomaly in hydrological variables is much stronger than on a short time scale [9,49]. By comparing the observed SDI-12 with the simulated SDI-12, the impacts of climate change and human activities on the hydrological drought were identified. Based on the SDI, hydrological droughts can be divided into five levels ( Table 6). In this study, the frequency of hydrological droughts was also used for quantifying the impacts of climate change and human activities on hydrological droughts. This can be calculated using the following formulas: where P is the frequency of hydrological droughts, n is the number of samples of hydrological droughts, and N is the total number of statistical samples. Due to human activity being relatively simple in this region, only LULC and reservoir construction were included. Based on the results of the NAM model, the individual contributions of LULC change and the Longtouqiao Reservoir to the changed volume of streamflow caused by human activities can be calculated as follows: where  where P is the frequency of hydrological droughts, n is the number of samples of hydrological droughts, and N is the total number of statistical samples. Due to human activity being relatively simple in this region, only LULC and reservoir construction were included. Based on the results of the NAM model, the individual contributions of LULC change and the Longtouqiao Reservoir to the changed volume of streamflow caused by human activities can be calculated as follows:      Table 7). The contributions of LULC change and the Longtouqiao Reservoir construction to the changed volume of the streamflow caused by human activities were simulated with Equations (1) and (2), respectively ( Table 6). The results showed that the streamflow, once affected by human activities, decreased by 5.  Table 9). According to Table 9, the slopes of the linear relationships of the accumulative streamflows decreased Similarly, the changes in observed cumulative precipitation (SP) and PET (SE) are shown in Figure 9. It was illustrated that the correlation coefficients are all higher than 0.99 during the three periods mentioned above. The slope change rates of the accumulative precipitation (RSP) and PET (RSE) are also listed in Table 8. Compared with the baseline condition, the slopes of the linear relationship of the accumulative precipitation decreased by 91. 46 Table 7). The contributions of LULC change and the Longtouqiao Reservoir construction to the changed volume of the streamflow caused by human activities were simulated with Equations (1) and (2), respectively ( Table 6). The results showed that the streamflow, once affected by human activities, decreased by 5.9 m 3 s −1 when compared with the natural streamflow (Q S 1 (2002−2014) ), and the contributions of LULC change (P LULC ) and the Longtouqiao Reservoir (P Res ) to this change were 38.3% and 61.7%, respectively. The fitted linear relationships between the year and observed cumulative streamflow during 1961-1966, 1967-2001, and 2002-2014 are shown in Figure 8. As shown in Figure 9, all of the relationships have high correlation coefficients above 0.98. The corresponding slopes (S R ) are listed in Table 8.  Table 9). According to Table 9, the slopes of the linear relationships of the accumulative streamflows decreased by 3.34 ×    Note: "-"means the change ratios are 0 compared with the baseline period.

The Characteristics of Streamflow Change
The contribution of climate change and human activities to the change of streamflow was simulated with Equations (3)-(7) ( Table 9). The quantitative assessment results showed that both climate change and human activities had a negative effect on streamflow in the headwater area of the NR. For the 1967-2001 period, the contributions of precipitation, potential evapotranspiration (PET), and LULC change were 35%, 12.7%, and 52.3%, respectively. During the period of 2002-2014, the contributions of precipitation, PET, LULC change, and the Longtouqiao Reservoir were 32.5%, 12.0%, 21.1% and 34.4%, respectively.    Note: "-"means the change ratios are 0 compared with the baseline period.
The contribution of climate change and human activities to the change of streamflow was simulated with Equations (3)-(7) ( Table 9). The quantitative assessment results showed that both climate change and human activities had a negative effect on streamflow in the headwater area of the NR. For the 1967-2001 period, the contributions of precipitation, potential evapotranspiration (PET), and LULC change were 35%, 12.7%, and 52.3%, respectively. During the period of 2002-2014, the contributions of precipitation, PET, LULC change, and the Longtouqiao Reservoir were 32.5%, 12.0%, 21.1% and 34.4%, respectively.  Note: "-"means the change ratios are 0 compared with the baseline period. Table 9. Contributions of climate and human activities to streamflow changes in the headwater area of the Naoli River in Northeast China. Similarly, the changes in observed cumulative precipitation (S P ) and PET (S E ) are shown in Figure 9. It was illustrated that the correlation coefficients are all higher than 0.99 during the three periods mentioned above. The slope change rates of the accumulative precipitation (R SP ) and PET (R SE ) are also listed in Table 8. Compared with the baseline condition, the slopes of the linear relationship of the accumulative precipitation decreased by 91. 46  The contribution of climate change and human activities to the change of streamflow was simulated with Equations (3)-(7) ( Table 9). The quantitative assessment results showed that both climate change and human activities had a negative effect on streamflow in the headwater area of the NR. For the 1967-2001 period, the contributions of precipitation, potential evapotranspiration (PET), and LULC change were 35%, 12.7%, and 52.3%, respectively. During the period of 2002-2014, the contributions of precipitation, PET, LULC change, and the Longtouqiao Reservoir were 32.5%, 12.0%, 21.1% and 34.4%, respectively.

The Characteristics of Hydrological Drought Changes
(1) Long-Term Hydrological Drought Hydrological droughts in the 1961-2014 period were analyzed using SDI for twelve months (SDI-12) (Figure 10). It was found that annual hydrological droughts did not occur during 1961-1966. However, the frequency of hydrological droughts has increased substantially from 1967 to 2014. Successive hydrological droughts have occurred in the following six periods: Two-year droughts in [1967][1968][1993][1994]

The Characteristics of Hydrological Drought Changes
(1) Long-Term Hydrological Drought Hydrological droughts in the 1961-2014 period were analyzed using SDI for twelve months (SDI-12) (Figure 10). It was found that annual hydrological droughts did not occur during 1961-1966. However, the frequency of hydrological droughts has increased substantially from 1967 to 2014. Successive hydrological droughts have occurred in the following six periods: Two-year droughts in [1967][1968][1993][1994]  (2) Contribution of Climate Change and Human Activities to Hydrological Drought The SDI-12 was calculated based on the simulated monthly streamflow without LULC and the reservoir during 1967-2014 ( Figure 11a) and 2002-2014 (Figure 11b). The frequency of hydrological droughts was calculated with Equation (8), and an obvious increase under the impacts of precipitation and PET over the past 50 years was noted with the effect of human activities. From 1967 to 2001, hydrological drought frequency increased from 53% to 58% under the impact of LULC changes. During the period of 2002-2014, the frequency of hydrological droughts without the impacts of LULC change and the Longtouqiao Reservoir was 55%. It increased to 63% with the impact of LULC change, and eventually increased to 74% when considering both the impacts of LULC change and the Longtouqiao Reservoir (Table 10). In addition to a more frequent occurrence of hydrological droughts, human activities also increased drought intensity. For instance, with the influence of LULC change, the hydrological droughts changed from mild drought to extreme drought in 1978, from mild drought to severe drought in 1979, and from moderate drought to extreme drought in 1993. Meanwhile, the hydrological droughts changed from mild drought to severe drought in 2002 after the Longtouqiao Reservoir started storing the upstream discharge. (2) Contribution of Climate Change and Human Activities to Hydrological Drought The SDI-12 was calculated based on the simulated monthly streamflow without LULC and the reservoir during 1967-2014 ( Figure 11a) and 2002-2014 (Figure 11b). The frequency of hydrological droughts was calculated with Equation (8), and an obvious increase under the impacts of precipitation and PET over the past 50 years was noted with the effect of human activities. From 1967 to 2001, hydrological drought frequency increased from 53% to 58% under the impact of LULC changes. During the period of 2002-2014, the frequency of hydrological droughts without the impacts of LULC change and the Longtouqiao Reservoir was 55%. It increased to 63% with the impact of LULC change, and eventually increased to 74% when considering both the impacts of LULC change and the Longtouqiao Reservoir (Table 10). In addition to a more frequent occurrence of hydrological droughts, human activities also increased drought intensity. For instance, with the influence of LULC change, the hydrological droughts changed from mild drought to extreme drought in 1978, from mild drought to severe drought in 1979, and from moderate drought to extreme drought in 1993. Meanwhile, the hydrological droughts changed from mild drought to severe drought in 2002 after the Longtouqiao Reservoir started storing the upstream discharge.

Streamflow Change of a Marshy River
Several studies have reported a decrease in streamflow in the headwater area of the NR in the past 54 years, and have attributed it to a precipitation decline in the region [22,24,50]. Overall, our findings showed a similar trend in the headwater area of the NR. Climate change plays an important role in streamflow change [51]. Compared with the period of 1961-1966, we found that there was a significant decreasing trend of precipitation in 1967-2001, and moreover, the trend was more obvious in 2002-2014. In addition, PET increased during the above two periods. The combination of these two factors has led to the negative effect of climate change on streamflow.
In addition, the exploitation and utilization of water resources are direct influencing factors for streamflow flowing into wetlands decreasing in the headwater area of the NR [51]. In order to develop the economic potential in the headwater area of the NR, a large amount of land was

Streamflow Change of a Marshy River
Several studies have reported a decrease in streamflow in the headwater area of the NR in the past 54 years, and have attributed it to a precipitation decline in the region [22,24,50]. Overall, our findings showed a similar trend in the headwater area of the NR. Climate change plays an important role in streamflow change [51]. Compared with the period of 1961-1966, we found that there was a significant decreasing trend of precipitation in 1967-2001, and moreover, the trend was more obvious in 2002-2014. In addition, PET increased during the above two periods. The combination of these two factors has led to the negative effect of climate change on streamflow.
In addition, the exploitation and utilization of water resources are direct influencing factors for streamflow flowing into wetlands decreasing in the headwater area of the NR [51]. In order to develop the economic potential in the headwater area of the NR, a large amount of land was cultivated into farmland (Table 1) [24,25]. Farmland increased from 3.5% (1961) to 34% (2001), while forest decreased from 75.6% (1961) to 54.6% (2001) across the entire region. This may have destroyed the natural hydrological environment of the underlying surface, and the rapid development of irrigated agriculture aggravated the exploitation of water resources in the headwater area of the NR. Hydraulic engineering has also had an important role in the decrease of streamflow. Specifically, the trend of streamflow reduction became more obvious after the Longtouqiao Reservoir was built. According to our study, the contribution of reservoirs to streamflow changes may be significantly higher than that of LULC changes. The primary function of this reservoir is to control floods and supply irrigation water. The average amount of water used for irrigation is about 2.2 × 10 8 m 3 yr −1 [25]. The irrigation area is still increasing, and the water supply will therefore likely continue to increase, which may lead to a further reduction in streamflow downstream due to the Longtouqiao Reservoir in the future.

Hydrological Drought of a Marshy River
Hydrological drought of the marshy Naoli River has not been thoroughly investigated in the past. Our findings from this study showed that the frequency of hydrological droughts has been increasing over the past 54 years under the impacts of climate change and human activities. This indicates that the Naoli River Basin has been on a trend towards drier conditions. In order to find the reasons for the increase of hydrologic drought, first of all, we eliminated the influence of human activities using a hydrological model (NAM). The results show that climate change was the main reason for frequent hydrological droughts in the headwater area of the NR. Second, the frequencies of hydrological droughts under the influence of human activities and without human activities were compared, and the results showed that the frequency was exacerbated, but not significantly. However, according to Figure 11, the greatest effect of human activities was that the intensity of hydrological droughts in this marshy river was exacerbated. The frequency of severe drought increased throughout 1967-2014, and extreme drought had occurred in 1978 and 1993. This may be due to increased water demand for irrigation and reservoir storage in dry years, which would exacerbate hydrological droughts downstream of the reservoir. This situation is rather unfavorable for the protection of riparian wetland.
There has been a lot of evidence demonstrating the importance of mitigating hydrological droughts for wetland protection [52][53][54][55]. The Naoli River Basin has a large area of wetlands, but most have been converted into farmland in the past half century [20], with the area of wetlands in the headwater area of the NR decreasing from 9.3% in 1961 to 0.8% in 2014. Therefore, the study of hydrological droughts will be beneficial to the protection of the remaining wetlands in the Naoli River Basin.
In spite of the uncertainties and limitations of methods, our study sheds light on the conditions of hydrological droughts in the headwater area of this marshy river, which is critical for the protection of ecosystems of river and wetlands. In addition, we plan to undertake more work on these uncertainties in future studies to improve the quantified results.

Conclusions
This study is the first to assess the individual contributions of precipitation change, potential evapotranspiration change, land use/land cover change, and a river dam construction to hydrological droughts in a marshy river in far Northeast China at higher latitudes. Based on our results, the following conclusions can be drawn: (1) Climate change, land use change (primarily wetland loss), and river dam construction have all contributed to streamflow decline in this marshy river, with the greatest contribution coming from the dam construction, followed by LULC change and meteorological factors. (2) While the streamflow of the marshy Naoli River continued declining, hydrological drought in the river has occurred more frequently and more intensely since 1967. The frequency of hydrological drought increased most obviously during the period of 2002-2014.
(3) Climate change was the main factor for the increase of drought events, while LULC change and the dam construction contributed to the severity of the droughts. These findings suggest that future water resource management in this and other agriculture-intensive river basins of high-latitude needs to develop effective strategies in reducing river drought occurrence and intensity.