Quantitative Assessment of Impact of Climate Change and Human Activities on Streamflow Changes Using an Improved Three-Parameter Monthly Water Balance Model

Understanding the impact of climate change and human activities on the hydrological cycle of any watershed can provide a scientific basis for regional water resource planning, flood management, and disaster mitigation. An improved three-parameter hydrological model (CM) based on monthly water balance using an exponential equation to depict the distribution of groundwater storage capacity was developed and evaluated. The model uses Asian Precipitation Highly Resolved Observational Data Integration Towards Evaluation (APHRODITE) rainfall data as input, with the Zhejiang Province as the case application, and the effects of climate change and human activities on streamflow changes were assessed by separating environmental variables in this study. The results indicate that APHRODITE data has excellent monthly accuracy, with a mean correlation coefficient (CC) of more than 0.96 and an average absolute percentage bias (Pbais) of less than 5%. The three models are relatively close in their ability to simulate high flows, but the CM simulated low flow is better than the other two models. Positive and negative Pbais phenomena occur in the CM model in each catchment, and absolute levels are regulated by 5%. Furthermore, the CM model’s average Nash efficiency coefficient (NSE) is greater than 0.9, indicating that it can correctly fulfill the water balance. The results are more consistent throughout multiple catchments in each watershed using Budyko-based and hydrological model technique to evaluate the influence of climate change and human activities on streamflow. Climate change dominated streamflow variations in 18 of the 21 catchments in Zhejiang Province, whereas human activities dominated the rest. The findings of the study will be used to influence the management, development, and usage of water resources in


Introduction
In recent decades, global climate change and human activities have significantly impacted the quantity, quality, and spatial distribution of water resources. Climate change is expected to disrupt the hydrological cycle, resulting in a regional and temporal redistribution of water resources, with significant implications for rainfall, evaporation, runoff, and soil water content [1][2][3][4][5][6]. Water supply, irrigation, navigation, urban drainage, hydropower generation, flood scheduling, soil and water conservation, and water pollution control will all be impacted by changes in the water cycle [7][8][9][10]. The hydrological cycle mechanisms are influenced by both climate change and human activity. Climate change influences watershed runoff in terms of input conditions, such as regional fluctuations in precipitation and temperature changes, while human activities can modify and affect the subsurface conditions directly or indirectly [3]. Long-term water cycle and water resource distribution have become important issues in water resource planning [11], and the monthly scale water balance model is frequently employed in long-term water resource distribution forecasting [12][13][14][15]. The reconstruction of catchment hydrological structures, forecasting of river runoff variations, the consequences of climate change in the catchment, and the evaluation of seasonal and regional water supply and irrigation demand are all examples of current research on monthly water balance models [16][17][18][19].
The hydrological model is a tool for characterizing and simulating the complex spatial and temporal variability in the conversion of rainfall to runoff from a watershed [20,21]. Hydrological models can be categorized into physical and conceptual models [22]. The physical model enables the study of the changing patterns of hydrological processes using physical scenarios, while the conceptual model is a hydrological model based on physical concepts in hydrological phenomena, as well as empirical formulas [11,23]. According to the degree of description of spatial characteristics, models can be classified as lumped, semi-distributed, and distributed hydrological models. Lumped hydrologic models describe the production of runoff using rainfall by taking the entire watershed into account, and the majority of these models lack the capability to mechanically account for the impact of unevenly distributed rainfall and subsurface conditions on the creation of rainfall runoff [24]. The distributed hydrological model divides the basin into several units based on the different information of climate and subsurface characteristics at each part of the basin, and it has the ability to mechanistically take into account the impact of the uneven spatial distribution of rainfall and subsurface conditions on runoff generation [25,26].
Hydrological models generally use mathematical and physical equations to simplify hydrological processes and model parameters representing watershed characteristics [27]. The number of parameters of hydrological models often decreases as the time scale increases. The daily Xinanjiang model, for example, includes 12 parameters, but the monthly ABCD model has only four [28]. The monthly water balance model provides a fundamental generalization of the catchment's hydrological processes, with relatively minimal model input requirements and quick estimates with a few model parameters [29][30][31]. Thornthwaite [32] developed the first monthly water balance model in the 1940s, which was later improved by Thornthwaite and Mather [33]. Many water-balance models based on various assumptions and model architectures were sequentially suggested by researchers to accommodate hydrological simulations with varied aims, such as drought analysis, water quantity estimation, and ecological studies [34][35][36]. For example, considering Thomas's [37] well-known ABCD monthly water balance model, Makhlouf and Michel [38] proposed the GR2M two-parameter monthly water balance model for evaluating and managing water resources, and Xiong and Guo [39] studied a two-parameter monthly water balance model and utilized it to simulate streamflow in Southern China.
Currently, the water balance method based on Budyko's hypothesis, and the hydrological simulation method are the main methods for studying the effects of climate change and human activities on runoff changes in watersheds [40,41]. Hydrological models with varying time scales and levels of complexity have been examined in recent years for hydrological simulations of climate change and human activities in watersheds [42][43][44]. The hydrological model simulation approach can separate the study period into two parts: the base period and the impact phase. The impact period is determined by the time of the abrupt change in streamflow, assuming that the streamflow in the base period is only influenced by climate change and that the streamflow in the impact period is influenced by both climate change and human activities. The hydrological model established in the base period is also used to predict the runoff in the impact period and then compared to the measured streamflow in the impact period in order to quantitatively separate the influence of climate change and human activities [45]. Zhang et al. [46] developed a variable time scale dynamic water balance model (DWBM) based on the Budyko assumption and conducted example research on 265 watersheds in Australia, demonstrating that the model's runoff simulation was sufficient. Wang and Tang [15] computed a single-parameter Budyko water balance model at the multi-year average scale based on the soil conservation service (SCS) model's isopycnic assumption and investigated the model parameters' key affecting aspects. Based on Wang and Tang's work, Zhao et al. [47] derived the maximum entropy generation principle for the same model structure at the multi-year average scale, monthly scale, and field-scale water balance models and elucidated the correlations between the water balance models at three different time scales.
Precipitation data are a critical input for hydrological models, and the hydrological community has struggled to overcome the challenge of hydrological simulation and forecasting in locations where the data is scarce. Precipitation data derived through satellite inversion and model operations offer a solution to the challenge of hydrological modeling in poorly informed areas [48]. The frequently used two-parameter monthly water balance hydrological model has excellent simulation capabilities and is very practical [39]. However, it does not include a baseflow simulation component. Because the storage capacity of subsurface runoff is not uniformly distributed, it is more appropriate to describe baseflow using an exponential equation to depict the distribution of groundwater storage capacity as a parabola. At the same time, reliable quantitative analysis of the impacts of climate change and human activities on streamflow variations is dependent on the quality of the hydrological model simulation. As a result, the primary goal of this research is to develop a simple and accurate monthly water balance model that can simulate the entire hydrological process, including streamflow, surface runoff, and baseflow. The objectives of this study are: (1) to evaluate and validate the APHRODITE satellite precipitation products for use in hydrological modeling, (2) propose and develop a three-parameter hydrological model based on monthly water balance, and (3) use the water balance model to assess the impact of climate change and human activities on streamflow changes. The Wang Model (WM), a two-parameter monthly water balance model, is a conceptually simple monthly water balance model that can be used to estimate monthly flows using monthly rainfall and potential evapotranspiration data [49]. The WM model has been successfully applied to semi-arid and humid watersheds in China with only three parameters [50]. In comparison to other hydrological models such as SIMHYD, TANK, and SMAR, it not only has a more straightforward structure, fewer parameters, and greater flexibility, but it also performs well. Different components of the model are estimated using Equations (1)- (3).
The variable E(t) is the actual monthly evaporation value of the catchment; PET(t) is the monthly potential evapotranspiration value; P(t) is the monthly precipitation; t is the time; Q S (t) is the monthly surface runoff volume; S(t) is the monthly soil water content; Q b (t) is the monthly baseflow volume; S max is the maximum monthly soil water content of the watershed; k s is the runoff coefficient; k g is the baseflow coefficient.

Improved Three-Parameter Monthly Water Balance Model
This study develops an improved three-parameter monthly water balance model called the Chen Model (CM). The model is constructed as shown in Figure 1. According Remote Sens. 2022, 14, 4411 4 of 22 to the innovation of Xiong and Guo [39], the actual evaporation from the watershed is calculated using the Ol'dekop equation [51]. It can be considered as a conversion factor from evaporation dish observations to actual evaporation values, which is a hyperbolic tanh() function of the ratio of precipitation to evaporation dish observations. The tanh(x) function for a given compression, x reflects the greater buffering and regulating capacity of watershed soils compared to air for hydrological phenomena or processes. Based on the monthly potential evapotranspiration values and monthly precipitation, the actual monthly evaporation values of the catchments were obtained using Equation (4) as follows: where k is the first parameter of the model to respond to the evaporation and precipitation changes.
content; is the monthly baseflow volume; is the maximum monthly soil water content of the watershed; is the runoff coefficient; is the baseflow coefficient. This study develops an improved three-parameter monthly water balance model called the Chen Model (CM). The model is constructed as shown in Figure 1. According to the innovation of Xiong and Guo [39], the actual evaporation from the watershed is calculated using the Ol'dekop equation [51]. It can be considered as a conversion factor from evaporation dish observations to actual evaporation values, which is a hyperbolic ℎ function of the ratio of precipitation to evaporation dish observations. The ℎ function for a given compression, reflects the greater buffering and regulating capacity of watershed soils compared to air for hydrological phenomena or processes. Based on the monthly potential evapotranspiration values and monthly precipitation, the actual monthly evaporation values of the catchments were obtained using Equation (4) as follows: where is the first parameter of the model to respond to the evaporation and precipitation changes. If the entire catchment is treated as a reservoir, assuming that the monthly surface runoff is a hyperbolic ℎ function of the soil water content. The monthly surface runoff is obtained from the net soil water content of the month (i.e., the soil water content of the previous period plus the precipitation of the month minus the actual evaporation of the month) using Equation (5) as follows: where is the second parameter of the model, which usually represents the maximum monthly soil water content of the entire watershed.
An exponential equation was utilized in the baseflow calculation (Equation (6)), but in the WM model, a linear equation was used. Because baseflow storage capacity is seldom equally distributed, a linear equation cannot accurately replicate the actual condition of baseflow. To solve the issue of continually changing groundwater storage capacity, the CM model mimics the unequal distribution of tension water by describing groundwater storage capacity distribution as a parabola, which can more precisely imitate the baseflow If the entire catchment is treated as a reservoir, assuming that the monthly surface runoff is a hyperbolic tanh() function of the soil water content. The monthly surface runoff is obtained from the net soil water content of the month (i.e., the soil water content of the previous period plus the precipitation of the month minus the actual evaporation of the month) using Equation (5) as follows: where SC is the second parameter of the model, which usually represents the maximum monthly soil water content of the entire watershed. An exponential equation was utilized in the baseflow calculation (Equation (6)), but in the WM model, a linear equation was used. Because baseflow storage capacity is seldom equally distributed, a linear equation cannot accurately replicate the actual condition of baseflow. To solve the issue of continually changing groundwater storage capacity, the CM model mimics the unequal distribution of tension water by describing groundwater storage capacity distribution as a parabola, which can more precisely imitate the baseflow than the WM model. Compared to linear techniques, using exponential equations to model baseflow is more consistent with groundwater outflow.
The third parameter d of the model represents the baseflow coefficient. The output value of the model is the total monthly runoff flow Q(t) (Equation (7)), and the entire catchment needs to meet the monthly water balance (Equation (8)) as follows: According to the physical significance of the three parameters of this model, the range of values of the three parameters can be given, and the range applies to areas with different geographical features, with k from 0.5 to 1.5, SC from 100 to 2000, and d from 0.01 to 1.00.

Four-Parameter Monthly Water Balance Model
The four-parameter monthly water balance model (ABCD) is a conceptual hydrological model proposed by Thomas [37], which has been applied in several wet or semi-humid basins in China, with the same input data as the WM model and an additional groundwater content in the output data [15]. The water balance in a soil aquifer is represented by Equation (9).
The two-state variables (effective water volume W(t) and possible evapotranspiration Y(t) in the ABCD model are calculated using Equations (10)-(12) as follows: W(t) is related to Y(t) as a nonlinear function as follows: where a is the probability of forming runoff before the soil is fully saturated and takes a value in the range of (0, 1]; b is the upper limit of unsaturated aquifer storage capacity and takes a value in the range of (0, 1000]. The relationship between soil moisture content and potential evapotranspiration is given by Equation (13).
The water balance equations of the model for the groundwater layer are as follows: where G(t) is the amount of underground storage in the current month; GD(t) is the underground runoff; G(t − 1) is the amount of underground storage in the previous month; c is the proportion of the amount of groundwater recharged by the soil aquifer, taking the value range of (0, 1]; d is the rate of groundwater formation outflow, taking the value range of (0, 1].

Performance Indicators
The accuracy of satellite rainfall data was evaluated by percentage bias (Pbais), root mean square error (RMSE), and correlation coefficient (CC). In addition, the effect of the hydrological model runoff simulation was evaluated using the Nash efficiency coefficient (NSE) (Nash and Sutcliff, 1970 [52]) and percentage bias (Pbais). The performance indicators are given by Equations (17)- (20).
where P oi is the observed rainfall; P si satellite rainfall; Q t obs is the observed runoff; Q t sim is the simulated runoff.

The Budyko-Based Method
The elasticity factor approach based on the Budyko framework combines the Budyko model with a watershed water balance model, where human activities and climate change are the main factors affecting runoff changes [53]. Climate change is mainly characterized as changes in precipitation and potential evapotranspiration, and changes in coupled hydrothermal parameters and other effects are described as anthropogenic impacts for quantitative calculations [54]. In a closed catchment, its multi-year water balance is expressed as: where P, E, and Q represent catchment precipitation, evaporation, and runoff, respectively; ∆S is the catchment water storage variable, which is negligible on a multi-year average scale. According to the assumption of Budyko [55], the annual average actual evapotranspiration can be estimated from precipitation and potential evapotranspiration, i.e., E/P = F(Φ), where Φ = PET /P, which describes the dryness of a region. The change in runoff due to climate change can be calculated by the method proposed by Milly and Dunne [56].
where ∆P and ∆PET are the changes in precipitation and potential evapotranspiration, respectively; ∂Q/∂P and ∂Q⁄∂PET are the sensitivity coefficients of runoff to precipitation and potential evapotranspiration, respectively. The sensitivity coefficients of runoff to precipitation and potential evapotranspiration can be calculated from Equation (21) and the Budyko series of hydrothermal equations. When based on the Foucauldian formula [57], the sensitivity coefficient can be expressed as: where w is a model parameter related to vegetation type, soil hydraulic properties, and topography terrain [58]. It can be calculated based on the annual average P, PET, and Q of the basin.

Method of Separating Environment Variables
The hydrological time series (especially the streamflow time series) are generally divided into "natural phase" and "human-influenced phase" according to the change nodes, i.e., "base period" and "impact period", depending on the degree of influence of human activities on the natural environment in the study area. Human activities and climate change have been the fundamental causes of the difference between the observed streamflow in the impact and base periods [46]. Human activities and climate change have not affected the observed streamflow series in the base period. The hydrological time series in the impact period is simulated using the enhanced three-parameter monthly water balance model with the streamflow data as the base period, while the observed streamflow series in the impact period derives from the combined influence of both. The difference between the simulated streamflow in the impact period and the measured streamflow in the base period is used to calculate the climate change component. In contrast, the difference between the simulated streamflow in the impact period and the measured streamflow in the same period is used to calculate the human activities component. The specific calculation formulas are as follows: where ∆R is the total streamflow change; R 1 is the measured streamflow during the impact period; R 0 is the measured streamflow during the base period; R is the simulated streamflow during the impact period; ∆R climate is the impact value of climate change on streamflow; ∆R human is the impact value of human activities on streamflow; η climate and η human are the contribution rates of climate change and human activities on streamflow change, respectively.

Study Area
Zhejiang Province lies in the southeastern section of China's mainland, bordering the East China Sea, a humid and rainy zone with an average annual precipitation of more than 1600 mm and an average annual precipitation day between 140 and 180 days. Zhejiang's geography trends from southwest to northeast, with high points in the southwest and low points in the northeast. The southwest is mostly mountain ranges with elevations of more than 1000 meters; the central part is hilly and catchment with an altitude of about 50-500 meters, and the northeast is plain with an altitude less than 10 meters, with the mountainous hills accounting for 70.4 percent of the province's land area. Many rivers are formed by the abundance of water and the steep terrain. From north to south, Zhejiang Province has major water systems, including the Tiaoxi (TX) River Basin, Canal Basin, Qiantangjiang (QTJ) River Basin, Yongjiang (YJ) River Basin, Jiaojiang (JJ) River Basin, Oujiang (OJ) River Basin, Feiyunjiang (FYJ) River Basin, Aojiang (AJ) River Basin, Zhoushan Island (ZSI) River Basin (Figure 2), as well as a slew of tiny rivers that run into the sea and neighboring provinces on their own. A large plain formed by rivers, lakes, and shallow sea deposits exists in northern Zhejiang and coastal areas, and regional rivers and lakes are joined to build a dense water network, establishing a unique river network water system known as "Jiangnan Water Town".
Canal Basin, Qiantangjiang (QTJ) River Basin, Yongjiang (YJ) River Basin, Jiaojiang (JJ) River Basin, Oujiang (OJ) River Basin, Feiyunjiang (FYJ) River Basin, Aojiang (AJ) River Basin, Zhoushan Island (ZSI) River Basin (Figure 2), as well as a slew of tiny rivers that run into the sea and neighboring provinces on their own. A large plain formed by rivers, lakes, and shallow sea deposits exists in northern Zhejiang and coastal areas, and regional rivers and lakes are joined to build a dense water network, establishing a unique river network water system known as "Jiangnan Water Town."

Data Sources
Streamflow data from 21 hydrological stations in Zhejiang Province were utilized as measured runoff data in this study, and the essential information for each hydrological station can be seen in Table 1. The China Meteorological Science Data Sharing Service (http://data.cma.cn, accessed on 12 May 2020) provided daily precipitation data, evaporation dish measurements, and the average temperature at meteorological stations in the watershed. The Thiessen polygon technique was used to calculate the watershed's monthly surface rainfall.
The APHRODITE was created by the Japan Meteorological Agency (JMA) Research Foundation and utilized to collect satellite rainfall data for this study. APHRODITE gathered rainfall station measurements from around Asia and created gridded daily precipitation data with high precision [59]. The spatial resolution of the 1951-2007 precipitation reanalysis data for Central Asia is 0.25° × 0.25°. (V1101). In 2018, the same technique was used to develop an extended update (V1101EX_R1) covering 2008-2015 for diagnosing climate change and associated hydrological investigations (V1101).

Data Sources
Streamflow data from 21 hydrological stations in Zhejiang Province were utilized as measured runoff data in this study, and the essential information for each hydrological station can be seen in Table 1. The China Meteorological Science Data Sharing Service (http://data.cma.cn, accessed on 12 May 2020) provided daily precipitation data, evaporation dish measurements, and the average temperature at meteorological stations in the watershed. The Thiessen polygon technique was used to calculate the watershed's monthly surface rainfall.
The APHRODITE was created by the Japan Meteorological Agency (JMA) Research Foundation and utilized to collect satellite rainfall data for this study. APHRODITE gathered rainfall station measurements from around Asia and created gridded daily precipitation data with high precision [59]. The spatial resolution of the 1951-2007 precipitation reanalysis data for Central Asia is 0.25 • × 0.25 • . (V1101). In 2018, the same technique was used to develop an extended update (V1101EX_R1) covering 2008-2015 for diagnosing climate change and associated hydrological investigations (V1101).

Assessment of APHRODITE Rainfall Product
Daily rainfall measurements from eight meteorological stations in the QTJ River basin from 1981 to 1995, as well as four meteorological stations in the OJ River basin from 1971 to 1998, were used to test the accuracy of APHRODITE rainfall products at the station scale in this study. The raster data of APHRODITE co-located with the observed rainfall stations were used for comparative analysis, and the results are the average of the comparison results of each station in each catchment. The scatter plots of the surface rainfall on the monthly time scale compared with the measured rainfall data at the rainfall stations are shown in Figure 3. The Pbais between APHRODITE satellite and measured rainfall at the monthly scale between the four catchments is 3.97%, −5.99%, −5.09%, and −3.08%, respectively, while the CC values are 0.97, 0.96, 0.97, and 0.98, and the RMSE values are 31.93, 24.62, 32.95, and 20.97. As a result, the APHRODITE satellite rainfall products offer excellent monthly accuracy.
Taylor diagrams at monthly scales for the catchments are given in Figure 4a to visualize the accuracy of satellite rainfall products at different scales and catchments. The length of its radius reflects the amount of change in the overall standard deviation, the angular relationship with the vertical coordinate represents the correlation coefficient, and the length from the observation point represents the error in Taylor diagrams. As a result, using the Taylor diagram helps to indicate how different precipitation products differ from recorded precipitation data at ground rainfall stations. The closer the simulated value's geographic distribution is to the observed value, the shorter the simulated value's distance from the observation point. Over diverse watersheds, the Taylor diagram displays the degree of superiority and inferiority of the APHRODITE satellite product estimations. At the watershed scale, as shown in Figure 4a, the WR station is closer to the observation point in the Taylor diagram, and the accuracy of the three sites (QZ, JH, and LX) within the QTJ River watershed is similar. As a result, the accuracy of APHRODITE satellite rainfall products on the OJ River basin is slightly greater than that of the QTJ River basin. The greater catchment area has a better hydrological model simulation performance, according to Table 1. The accuracy of APHRODITE satellite rainfall products is much greater in winter than in other seasons and somewhat less applicable in summer than in other seasons, as illustrated in Figure 4b. At the same time, the watershed scale provides a broader range of applications at the WR station. Remote Sens. 2022, 14, x FOR PEER REVIEW 10 of 24 Taylor diagrams at monthly scales for the catchments are given in Figure 4a to visualize the accuracy of satellite rainfall products at different scales and catchments. The length of its radius reflects the amount of change in the overall standard deviation, the angular relationship with the vertical coordinate represents the correlation coefficient, and the length from the observation point represents the error in Taylor diagrams. As a result, using the Taylor diagram helps to indicate how different precipitation products differ from recorded precipitation data at ground rainfall stations. The closer the simulated value's geographic distribution is to the observed value, the shorter the simulated value's distance from the observation point. Over diverse watersheds, the Taylor diagram displays the degree of superiority and inferiority of the APHRODITE satellite product estimations. At the watershed scale, as shown in Figure 4a, the WR station is closer to the observation point in the Taylor diagram, and the accuracy of the three sites (QZ, JH, and LX) within the QTJ River watershed is similar. As a result, the accuracy of APHRODITE satellite rainfall products on the OJ River basin is slightly greater than that of the QTJ River basin. The greater catchment area has a better hydrological model simulation performance, according to Table 1. The accuracy of APHRODITE satellite rainfall products is much greater in winter than in other seasons and somewhat less applicable in summer than in other seasons, as illustrated in Figure 4b. At the same time, the watershed scale provides a broader range of applications at the WR station.

Simulation Analysis of Monthly Water Balance Model
This study's examination of the monthly water balance model's simulation results is separated into two parts: (1) actual rainfall data from two watersheds (four catchments) were used to compare the simulation results of the three-parameter monthly water balance (CM) model with those of the WM and ABCD models; (2) APHRODITE precipitation data were used as the input to the hydrological model, and the CM, WM, and ABCD models for eight watersheds (21 catchments in total) in the Zhejiang Province.

Simulation Analysis of Monthly Water Balance Model
This study's examination of the monthly water balance model's simulation results is separated into two parts: (1) actual rainfall data from two watersheds (four catchments) were used to compare the simulation results of the three-parameter monthly water balance (CM) model with those of the WM and ABCD models; (2) APHRODITE precipitation data were used as the input to the hydrological model, and the CM, WM, and ABCD models for eight watersheds (21 catchments in total) in the Zhejiang Province.

Simulation Results Using the Observed Rainfall
The WM, ABCD, and CM models were used to maximize the values of each model parameter in four catchments of the QTJ and OJ River basins. The NSE value based on simulated and measured streamflow was selected as the objective function for obtaining the optimal model parameters. A genetic algorithm-based solution approach is used to search for the optimal parameter set. Table 2 shows the parameter values for three hydrological models, and Table 3 shows the monthly streamflow simulation accuracy comparison for each model over the calibration and validation period. Tables 2 and 3 show that the established water balance models in each catchment have great simulation accuracy, with the NSE values of the CM model over the validation period above 0.9 and substantially exceeding the simulation results of the WM and ABCD models. The Pbais of the WM model simulations are all positive, showing that the model is overestimated, but the percentage deviations of the CM model in each catchment are both positive and negative. The absolute numbers are regulated by 5%, allowing for proper water balance. Figure 5 depicts the simulated hydrological processes of each model in each catchment, as well as the simulated measured scatter plots. The simulation results of the CM model are closer to the actual hydrological process than the WM and ABCD models, as seen in the above figure, and the simulated and observed scatter spots are closer to the 1:1 line. The simulation ability of the three models for high flows is relatively close, but the CM model performs better in the simulation of low flows, and the simulation of the baseflow through the exponential equation is more consistent with the hydrological processes in the semi-humid and semi-arid regions of southern China. S max = maximum monthly soil water content of the watershed, k s = runoff coefficient, k g = baseflow coefficient; a = probability of forming runoff before the soil is fully saturated, b = upper limit of unsaturated aquifer storage capacity, c = proportion of the amount of groundwater recharged by the soil aquifer, d = rate of groundwater formation outflow; k = evaporation and precipitation changes, SC = maximum monthly soil water content of the entire watershed, d = baseflow coefficient.

Simulation Results Using the APHRODITE Rainfall
The WM, ABCD, and CM models were applied to the QTJ River Basin (QZ, LX, ZG, JH, ZX, XF, SZ, and QSD stations), YJ River Basin (NXK station), OJ River Basin (LQ, WR, JJK, SB, BY, and SZ stations), FYJ River Basin (XK station), AJ River Basin (DT station), JJ River Basin (BZA and SD stations), TX River Basin (QDC station) and ZSI River Basin (CCL station), and the model NSE and Pbais for the calibration and validation period are shown in Table 4. The results show that the simulation accuracy of the established water balance models for each catchment is high, the NSE for each catchment during the calibration and validation period is greater than 0.7, the absolute values of Pbais are within 15%, and each model can satisfy the water balance using APHRODITE rainfall as the model input. The average NSE for the CM model validation period is about 0.85, whereas the average NSE values for the WM and ABCD models are 0.79 and 0.82, respectively. On the other hand, the Pbais value for the CM model in the validation period, on the other hand, is closer to 0 and less than 5%, suggesting that the model matches the monthly water balance within each watershed, showing that the CM model has suitable application in Zhejiang Province's basins.
The streamflow above the 75th percentile is defined as peak flow, while that below the 25th percentile is considered as low flow. The peak flow and low flow simulation performance indicators for each model are shown in Figure 6. The three models are relatively close in their ability to simulate high flows, and this can be attributed to their reasonably accurate calculations for both actual evaporation and soil water content. All models show a good low-flow simulation capability due to the higher accuracy of precipitation input used in these models. It is known from the previous section that APHRODITE precipitation product is more applicable in winter than in other seasons and that low flows of hydrological processes tend to occur in winter. However, the CM simulated low flow is better than the other two models in terms of NSE and Pbais values. The reason for this is that both WM and ABCD models use linear equations to simulate the receding process, while CM uses exponential expressions that are more consistent with the actual baseflow process. To illustrate this phenomenon more closely, Figure 7 depicts the simulation results of hydrological processes in the LQ catchment of the OJ River basin in order to better observe the simulation of hydrological processes by each model. All models can accurately mimic the monthly runoff process for the LQ catchment, and the simulated values of the peak streamflow have high consistency with observations. At the same time, the CM model's simulation of the hydrological process is more accurate than the WM and ABCD models, particularly in terms of peak and minimum values, and the peak-shaped process is closer to the experimental process. According to the results in Tables 3 and 4, the rain gauge product provides better streamflow estimates than the APHRODITE product for the three given models. However, the difference in simulation performance between the two is not significant, and both provide a better representation of the basic hydrologic processes in the catchment.
This part splits the 21 catchments into three ranges by comparing the NSE and Pbais values of the model within various catchment regions (see Figure 8) to explore the application of the three-parameter monthly water balance model at the geographical scale: 1 catchment area < 1000 km 2 ; 2 1000 km 2 < catchment area < 5000 km 2 ; 3 catchment area > 5000 km 2 . When the catchment area exceeds 5000 km 2 , the model's simulation capability becomes critical. Using the NSE as an example, the CM model's average NSE value is greater than 0.9, and the absolute value of the mean Pbais is less than 1% when the catchment area is greater than 5000 km 2 . The reason for this phenomenon is that the bias of APHRODITE rainfall product is smaller in watersheds with large catchment areas, while small catchment areas are more challenging for the simulation capability of the ensemble hydrological model. It should be noted that the accuracy of each model simulation result is slightly poor owing to the effect of rainfall data or human activities in different watersheds.
The following study will investigate the impact of climate change and human activities on the change in watershed runoff in-depth to solve this problem. To illustrate this phenomenon more closely, Figure 7 depicts the simulation results of hydrological processes in the LQ catchment of the OJ River basin in order to better observe the simulation of hydrological processes by each model. All models can accurately mimic the monthly runoff process for the LQ catchment, and the simulated values of the peak streamflow have high consistency with observations. At the same time, the CM model's simulation of the hydrological process is more accurate than the WM and ABCD models, particularly in terms of peak and minimum values, and the peak-shaped process is closer to the experimental process. According to the results in Tables 3 and 4, the rain gauge product provides better streamflow estimates than the APHRODITE product for the three given models. However, the difference in simulation performance between the two is not significant, and both provide a better representation of the basic hydrologic processes in the catchment.
This part splits the 21 catchments into three ranges by comparing the NSE and Pbais values of the model within various catchment regions (see Figure 8) to explore the application of the three-parameter monthly water balance model at the geographical scale: ① catchment area < 1000 km 2 ; ② 1000 km 2 < catchment area < 5000 km 2 ; ③ catchment area > 5000 km 2 . When the catchment area exceeds 5000 km 2 , the model's simulation capability becomes critical. Using the NSE as an example, the CM model's average NSE value is greater than 0.9, and the absolute value of the mean Pbais is less than 1% when the catchment area is greater than 5000 km 2 . The reason for this phenomenon is that the bias of APHRODITE rainfall product is smaller in watersheds with large catchment areas, while small catchment areas are more challenging for the simulation capability of the ensemble hydrological model. It should be noted that the accuracy of each model simulation result is slightly poor owing to the effect of rainfall data or human activities in different watersheds. The following study will investigate the impact of climate change and human activities on the change in watershed runoff in-depth to solve this problem.

Impacts of Climate Change and Human Activities on Streamflow Changes
The hydrological cycle in the basin has changed due to the dual effect of climate change and human activity, causing the rainfall-runoff process in the basin to deviate from the historical hydrological process. In this changing environment, the first step in the hydrological study is to determine if the basin's rainfall-runoff mechanism fits the stationarity assumption. The Mann-Kendall test (Mann, 1945 [60]) was used for each watershed in Zhejiang Province to assess the abrupt change points in the runoff series: under the same coordinate system, two statistical curves, UF and UB, as well as the upper and lower critical lines. If UF and UB overlap within the upper and lower crucial lines, the year when the junction occurs is the mutation year. The significance level α = 0.05 was chosen, and the corresponding critical value U 0.05 = ±1.96. The Mann-Kendall test curves of annual runoff at different time scales for 21 catchments in 8 River basins are shown in Figure 9. The hydrological cycle in the basin has changed due to the dual effect of climate change and human activity, causing the rainfall-runoff process in the basin to deviate from the historical hydrological process. In this changing environment, the first step in the hydrological study is to determine if the basin's rainfall-runoff mechanism fits the stationarity assumption. The Mann-Kendall test (Mann, 1945 [60]) was used for each watershed in Zhejiang Province to assess the abrupt change points in the runoff series: under the same coordinate system, two statistical curves, UF and UB, as well as the upper and lower critical lines. If UF and UB overlap within the upper and lower crucial lines, the year when the junction occurs is the mutation year. The significance level α = 0.05 was chosen, and the corresponding critical value U .
1.96. The Mann-Kendall test curves of annual runoff at different time scales for 21 catchments in 8 River basins are shown in Figure 9. Based on the findings of the streamflow abrupt change study, a Budyko-based and hydrological model technique was used to evaluate the influence of climate change and human activities on streamflow. The research period was divided into two halves of different environmental factors: the base period and the impact period. The CM model was Based on the findings of the streamflow abrupt change study, a Budyko-based and hydrological model technique was used to evaluate the influence of climate change and human activities on streamflow. The research period was divided into two halves of different environmental factors: the base period and the impact period. The CM model was used to investigate the causes of streamflow changes in each watershed in Zhejiang Province. Table 5 shows the results of the two techniques' calculations, and the two methods' estimations of streamflow variations resulting from climate change and human activities are primarily comparable. The contribution of climate change and human activities to streamflow changes using the hydrological model method in each watershed is shown in Figure 10. It can be observed that climatic change dominates streamflow variations in 18 of the 21 catchments, while human activities contributed to those variations in the other three. Furthermore, the results are more consistent throughout multiple catchments in each watershed. In the OJ River basin, all stations show a decrease in streamflow due to human activities, but an increase in streamflow caused by climate change. Climate change contributes 51.39 percent to 87.39 percent in the QTJ River basin, with only declining streamflow at SZ station being a direct result of it. Except for the NXK watershed, climate change dominance spans from 60.97 percent to 80.74 percent throughout the OJ River basin. The contribution rates of human activities impacting yearly streamflow fluctuations in the JJ River basin (BZA and SD catchments) were 60.51 percent and −91.99 percent, respectively. Climate change and human activities impact the geographical distribution of yearly runoff from watersheds in Zhejiang Province (see Figure 11). Climate change has a greater impact on the province's watersheds than human activity. However, human activities have a greater effect on the eastern coastal areas than climate change, which might be supported by land-use changes (rural to urban) or other changes in those watersheds as shown to be impacted by human activities.

Discussions
It is worth pointing out that rainfall data are crucial for creating accurate models of hydrological processes as it is the main input to the models. As all the research regions

Discussions
It is worth pointing out that rainfall data are crucial for creating accurate models of hydrological processes as it is the main input to the models. As all the research regions

Discussions
It is worth pointing out that rainfall data are crucial for creating accurate models of hydrological processes as it is the main input to the models. As all the research regions are mesoscale watersheds with a maximum streamflow propagation time of less than 1 day from the uppermost section of the watershed control area to the streamflow station, a monthly temporal scale was employed to analyze the spatial distribution of satellite rainfall products. Additionally, studies [61,62] have indicated that watersheds with quicker production and confluence times are less affected by the spatial distribution of rainfall. The method for comparing satellite precipitation products with observation precipitation is another crucial issue that has to be considered. Thiessen polygon, a spatial interpolation technique, is not a reasonable way to compare product correctness [63]. As a result, the scale of the meteorological station was chosen for comparison. The APHRODITE dataset is the only high-resolution, long-term gridded information on land rainfall covering the entire Asian region, and the rainfall estimates are highly accurate at the monthly temporal scale [64,65]. Rainfall data from APHRODITE was selected as the model input because the hydrological model utilized in this work operates on a monthly scale. The monthly water balance model is a conceptual hydrological model based on the principle of monthly water balance [66,67], which has the advantages of simple structure and fewer parameters. In this study, a three-parameter monthly water balance model is modified, and the baseflow calculation was carried out using an exponential equation. The three models are fairly comparable in terms of their ability to simulate high flows, but the CM model is more effective at simulating low flows, and simulating baseflow using an exponential equation is more in line with the hydrological processes in the semi-humid and semi-arid regions of southern China. Further research is needed on the applicability of the model in other regions, such as arid regions. On the other hand, the baseflow simulation capability of this hydrological model cannot be evaluated due to the limitation of the accurate baseflow, but the improved model is improved in terms of the simulation capability of low flows.
When using this hydrological model to simulate watershed responses under changing environments, more climatic factors and human activity factors (e.g., land use type and vegetation cover conditions) are not considered in relation to runoff changes due to the limitations of the model characteristics. Different types of human activities may have different effects on runoff changes, so further research can be conducted in the future to separate the effects of different types of human activities on runoff changes. Furthermore, many studies [68][69][70] on attribution analysis are founded on the premise that human activity and climate change are independent of one another. However, in actual watersheds, human activity and climate change interact with one another, so future research is expected to continue to investigate how to conduct attribution analysis in these situations.

Conclusions
In this study, the accuracy of APHRODITE rainfall data is examined and studied, which uses eight river systems in Zhejiang Province as research objects. The research area's runoff process is simulated and tested using an enhanced three-parameter monthly water balance model. Finally, by isolating environmental factors, the impacts of climate change and human activities on runoff changes were assessed. The key findings are stated below: (1) The APHRODITE satellite rainfall products offer great accuracy on a monthly scale, with a CC mean of more than 0.96 and an average absolute Pbais of less than 5%. In addition, the accuracy of APHRODITE rainfall data is higher in the winter than in the other seasons and slightly less so in the summer.
(2) Three monthly water balance models offer good simulation accuracy and hydrological process inversion. Furthermore, the average NSE utilizing the CM model is greater than 0.9, indicating that the water balance may be accurately met. The CM model mimics the hydrological process better than the WM and ABCD models, particularly in the modeling of low flow processes, which is closer to the experimental process.
(3) In their estimations of streamflow variations as impacted by climate change and human activities, the Budyko-based and hydrological model methodologies are typically compatible. Climate change is the primary driver of streamflow variations in 18 of Zhejiang Province's 21 watersheds, while human activities are the primary driver in the other three. Climate change has a more significant impact on most river systems than human activities, although human activities have a greater impact on eastern coastal areas than climate change. Data Availability Statement: The daily precipitation data, evaporation dish observations, and average temperature at meteorological stations in the basin are obtained from the China Meteorological Science Data Sharing Service (http://data.cma.cn, accessed on 10 January 2020). The APHRODITE dataset is available from https://www.chikyu.ac.jp/precip/ (accessed on 12 May 2020).