A Modiﬁed ABCD Model with Temperature-Dependent Parameters for Cold Regions: Application to Reconstruct the Changing Runo ﬀ in the Headwater Catchment of the Golmud River, China

: The runo ﬀ changes due to global warming in hydrological basins in the Qinghai–Tibetan Plateau (QTP) have received worldwide attention. The headwater catchment of the Golmud River, located in the northern QTP, is the main source of water resources for the Golmud city in an arid region but has been poorly known for the hydroclimatological behaviors. In this study, a widely-used hydrological model, the ABCD model (Thomas, H.A., Washington, DC, USA), is modiﬁed by incorporating temperature-dependent hydrological processes and groundwater evapotranspiration in cold regions with a few additional parameters. The new model is used to reconstruct the monthly runo ﬀ in the past decades for the headwater catchment of the Golmud River and performs better than other comparable models. As indicated, the annual snowmelt runo ﬀ increased with the increasing air temperature and became more concentrated in April than in May. The frozen soil degradation could increase the hydraulic conductivity of soils and lead to a rise in cold season runo ﬀ . The groundwater level in the Golmud city was positively correlated to the annual runo ﬀ in the headwater catchment of the Golmud River, indicating that an increase of the groundwater level could be triggered by the rising trend in the streamﬂow of the Golmud River. This study suggests a useful hydrological model for the groundwater management in the Golmud city. data agreed well with the observed precipitation. This agreement illustrates


Introduction
Catchments in mountains generally role as headwater region for adjacent lowlands, providing variable water resources with the changing runoff that have great impacts on humanity and ecosystems in the middle and lower reaches area [1]. Most of the large mountains in the world are cold with a local climate that is sensitive to the rapid global warming in the last decades. Climate change caused by global warming may significantly alter the hydrological processes in headwater catchments, leading to a remarkable shift in the water balance, available water resources and extremes events (floods and droughts) [2][3][4]. Therefore, it is a great challenge to assess or predict the change in the runoff of headwater catchments for better managements of water use and ecosystem protections.
The Qinghai-Tibetan Plateau (QTP) roles as the water tower of Asia, which is the source area of many large Asian rivers, such as the Yellow River, Yangtze River, and Brahmaputra River, providing water resources for more than 1.4 billion of people [5]. The QTP has experienced significant warming The leakage of water from the Golmud River is the major recharge of groundwater in the area of the Golmud city. This recharge highly depends on the runoff of the headwater catchment of the Golmud River. As a consequence, the groundwater level in the Golmud city is subject to the streamflow in the river. Due to changes in streamflow, the depth to subsurface water table in Golmud City has experienced large oscillations in the past 40 years. In 1981In , 1983, and 2010, significant rise of the groundwater level caused serious disasters in the city area such as house collapses and grassland immersion [51,52]. Previous studies on this kind of groundwater disasters mainly focused on the inter-annual change in groundwater level [53,54], but it has been highlighted that groundwater in this area is closely related to surface water [55]. However, one of the difficulties in the management of the Golmud groundwater is the lack of understanding of long-term runoff changes in the headwater catchment of the Golmud River, where frozen soils exist and streamflow is sensitive to climate change in the QTP.

Data and Treatment
Several datasets were used in this study to prepare the input of the hydrological models, including the distributed precipitation (P), air temperature (T), and potential evaporation (E0). As shown in Figure 1, there was no national meteorological station in the headwater catchment of the Golmud River, although some meteorological stations exist around the study area. To estimate the basin-averaged monthly precipitation, the recently released global precipitation product TERRACLIMATE data [56] was used, which has a spatial resolution of 4 km. Among the nearby national meteorological stations (Figure 1), Golmud (with the lowest altitude) and Qumarleb (with the highest altitude), were selected to examine the accuracy of the TERRACLIMATE precipitation data. The observed mean monthly precipitation during 1965-2015 at the two meteorological stations are shown in Figure 2 and are compared with the TERRACLIMATE data at pixel elements where the stations located in. It indicates that the TERRACLIMATE data agreed well with the observed The leakage of water from the Golmud River is the major recharge of groundwater in the area of the Golmud city. This recharge highly depends on the runoff of the headwater catchment of the Golmud River. As a consequence, the groundwater level in the Golmud city is subject to the streamflow in the river. Due to changes in streamflow, the depth to subsurface water table in Golmud City has experienced large oscillations in the past 40 years. In 1981In , 1983, and 2010, significant rise of the groundwater level caused serious disasters in the city area such as house collapses and grassland immersion [50,51]. Previous studies on this kind of groundwater disasters mainly focused on the inter-annual change in groundwater level [52,53], but it has been highlighted that groundwater in this area is closely related to surface water [54]. However, one of the difficulties in the management of the Golmud groundwater is the lack of understanding of long-term runoff changes in the headwater catchment of the Golmud River, where frozen soils exist and streamflow is sensitive to climate change in the QTP.

Data and Treatment
Several datasets were used in this study to prepare the input of the hydrological models, including the distributed precipitation (P), air temperature (T), and potential evaporation (E 0 ). As shown in Figure 1, there was no national meteorological station in the headwater catchment of the Golmud River, although some meteorological stations exist around the study area. To estimate the basin-averaged monthly precipitation, the recently released global precipitation product TERRACLIMATE data [55] was used, which has a spatial resolution of 4 km. Among the nearby national meteorological stations (Figure 1), Golmud (with the lowest altitude) and Qumarleb (with the highest altitude), were selected to examine the accuracy of the TERRACLIMATE precipitation data. The observed mean monthly precipitation during 1965-2015 at the two meteorological stations are shown in Figure 2 and are compared with the TERRACLIMATE data at pixel elements where the stations located in. It indicates that the TERRACLIMATE data agreed well with the observed precipitation. This agreement illustrates Water 2020, 12, 1812 5 of 23 that the TERRACLIMATE precipitation data are reliable in the study area and can be used to estimate the average precipitation from this distributed data.
Water 2020, 12, x FOR PEER REVIEW 5 of 23 precipitation. This agreement illustrates that the TERRACLIMATE precipitation data are reliable in the study area and can be used to estimate the average precipitation from this distributed data. The distributed potential evaporation and air temperature were interpolated with distance weighting method from six national meteorological stations around the study area ( Figure 1). The correlation between the air temperature and the elevation [51] was incorporated in generating the distributed air temperature. The potential evaporation was calculated as the pan evaporation (Epan) observed from a standard pan (601 mm in diameter) at each meteorological station multiplying an empirical coefficient [57]. Daily T, P and Epan at the meteorological stations were obtained from the dataset released by the China Meteorological Data Service Center (http://data.cma.cn). These data were then used to estimate the monthly and annual T, P, and E0. A brief summary of the meteorological data was provided in Table 1. The monthly and annual T and E0 data at these meteorological stations were interpolated to produce 4-km resolution raster data in the study area. Applying the raster data of P extracted from TERRACLIMATE and the raster data of T and E0 interpolated from the meteorological stations, the basin-average monthly and annual T, P, and E0 The distributed potential evaporation and air temperature were interpolated with distance weighting method from six national meteorological stations around the study area ( Figure 1). The correlation between the air temperature and the elevation [50] was incorporated in generating the distributed air temperature. The potential evaporation was calculated as the pan evaporation (E pan ) observed from a standard pan (601 mm in diameter) at each meteorological station multiplying an empirical coefficient [56]. Daily T, P and E pan at the meteorological stations were obtained from the dataset released by the China Meteorological Data Service Center (http://data.cma.cn). These data were then used to estimate the monthly and annual T, P, and E 0 . A brief summary of the meteorological data was provided in Table 1. The monthly and annual T and E 0 data at these meteorological stations were interpolated to produce 4-km resolution raster data in the study area. Applying the raster data of P extracted from TERRACLIMATE and the raster data of T and E 0 interpolated from the meteorological stations, the basin-average monthly and annual T, P, and E 0 were obtained for the study catchment. Figure 3 shows the annual data between 1965 and 2015 at the basin-scale. As exhibited, there were increasing trends of P and T in the period 1965-2015 with Water 2020, 12, 1812 6 of 23 different fluctuations. The trend of E 0 was negative in this period and the fluctuation seemed to be also in contrast to that of T. This decreasing trend in E 0 is consistent with the decreasing trend in E pan indicated by meteorological observations in the past 50 years around China [57,58]. Existing analyses found that this negative trend in E pan was linked with the decreasing wind speed and solar radiation as well as the increasing relative humidity [57,58], even though the increasing air temperature was expected to be an positive factor for the evaporation process. Indeed, such a "paradox" did not only exist in China but also was found to be a global phenomenon, which has been interpreted as a result of significant decreases in sunlight due to increased cloud coverage and aerosol concentration [59]. Linear slope in trends of T, P, and E 0 at different meteorological stations are listed in Table 1. It is indicated that the slope of T was relatively uniform for different places but the slope of P and E 0 show significant spatial variance. In particular, the trend of E 0 was positive in some places, not always negative. These differences imply the complexity of the climate change in the QTP.
Water 2020, 12, x FOR PEER REVIEW 6 of 23 were obtained for the study catchment. Figure 3 shows the annual data between 1965 and 2015 at the basin-scale. As exhibited, there were increasing trends of P and T in the period 1965-2015 with different fluctuations. The trend of E0 was negative in this period and the fluctuation seemed to be also in contrast to that of T. This decreasing trend in E0 is consistent with the decreasing trend in Epan indicated by meteorological observations in the past 50 years around China [58,59]. Existing analyses found that this negative trend in Epan was linked with the decreasing wind speed and solar radiation as well as the increasing relative humidity [58,59], even though the increasing air temperature was expected to be an positive factor for the evaporation process. Indeed, such a "paradox" did not only exist in China but also was found to be a global phenomenon, which has been interpreted as a result of significant decreases in sunlight due to increased cloud coverage and aerosol concentration [60]. Linear slope in trends of T, P, and E0 at different meteorological stations are listed in Table 1. It is indicated that the slope of T was relatively uniform for different places but the slope of P and E0 show significant spatial variance. In particular, the trend of E0 was positive in some places, not always negative. These differences imply the complexity of the climate change in the QTP.  before 1965 at the outlet of the headwater catchment of the Golmud River but destroyed by a heavy flood in 1990. The second hydrological station, W-IV, was established in 1990 at a lower reach to replace the W-III station, where the streamflow changed since 1990 as shown in Figure 4c. Several reservoirs have been constructed in the river between W-III and W-IV since 1977. A significant leakage of water from the reservoirs and riverbeds exists between W-III and W-IV because the groundwater level over there is lower than the surface water level. Thus, the streamflow at W-IV does not represent the actual runoff in the study catchment. However, the annual streamflow during 1990-2015 at W-III has been estimated from the W-IV data using a restoration method in a previous study [54]. The combined annual streamflow data at W-III are shown in Figure 4d. As indicated, the annual runoff shows a significant fluctuation, partly following the variations of the local precipitation. The relationship between the annual streamflow and groundwater level is discussed in Section 4.5.
Water 2020, 12, x FOR PEER REVIEW 7 of 23 were monitored in observation wells G1-G4 (Figure 4a), with available data from 1990 to 2012 as shown in Figure 4b. The regime of groundwater levels was significantly influenced by the streamflow in the Golmud River observed at hydrological stations. The first hydrological station, W-III, was established before 1965 at the outlet of the headwater catchment of the Golmud River but destroyed by a heavy flood in 1990. The second hydrological station, W-IV, was established in 1990 at a lower reach to replace the W-III station, where the streamflow changed since 1990 as shown in Figure 4c. Several reservoirs have been constructed in the river between W-III and W-IV since 1977. A significant leakage of water from the reservoirs and riverbeds exists between W-III and W-IV because the groundwater level over there is lower than the surface water level. Thus, the streamflow at W-IV does not represent the actual runoff in the study catchment. However, the annual streamflow during 1990-2015 at W-III has been estimated from the W-IV data using a restoration method in a previous study [55]. The combined annual streamflow data at W-III are shown in Figure 4d. As indicated, the annual runoff shows a significant fluctuation, partly following the variations of the local precipitation. The relationship between the annual streamflow and groundwater level is discussed in Section 4.5. denotes the restored streamflow at W-III obtained by considering the leakage of river water between W-III and W-IV [55]. The annual precipitation at the Golmud city is also shown as histograms in (c) and (d).

ABCD Model
The ABCD model divides the whole catchment water storage into soil water and groundwater [42], as that shown in Figure 5a. The rainfall (P) is partitioned into the evapotranspiration (E), the surface runoff (D), and the change in soil water storage. A delayed infiltration from the soil water recharges the groundwater. The total runoff is a combination of the surface runoff and the base flow that released from groundwater. . Streamflow observed at W-IV (c) and W-III (d) in different periods. Dashed line in (c) denotes the restored streamflow at W-III obtained by considering the leakage of river water between W-III and W-IV [54]. The annual precipitation at the Golmud city is also shown as histograms in (c) and (d).

ABCD Model
The ABCD model divides the whole catchment water storage into soil water and groundwater [41], as that shown in Figure 5a. The rainfall (P) is partitioned into the evapotranspiration (E), the surface runoff (D), and the change in soil water storage. A delayed infiltration from the soil water recharges the groundwater. The total runoff is a combination of the surface runoff and the base flow that released from groundwater. For a monthly simulation, the change in soil water storage in the ith month is calculated as follows: where −1 and are, respectively, the effective soil water storages at the end of (i − 1)th and ith month, , , and are the precipitation, actual evapotranspiration, groundwater recharge and surface runoff in the ith month, respectively. The change in groundwater storage is calculated as follows: where and −1 are, respectively, the effective groundwater storages at the at the end of (i − 1)th and ith month, is the base flow in the ith month. The ABCD model estimates , , and as follows: where and are two dimensionless parameters for the monthly simulation, whereas and are two state variables defined, respectively, as follows: For a monthly simulation, the change in soil water storage in the ith month is calculated as follows: where S i−1 and S i are, respectively, the effective soil water storages at the end of (i − 1)th and ith month, P i , E i , R i and D i . are the precipitation, actual evapotranspiration, groundwater recharge and surface runoff in the ith month, respectively. The change in groundwater storage is calculated as follows: where G i and G i−1 are, respectively, the effective groundwater storages at the at the end of (i − 1)th and ith month, F i is the base flow in the ith month. The ABCD model estimates R i , G i , and F i as follows: where c and d are two dimensionless parameters for the monthly simulation, whereas W i and Y i are two state variables defined, respectively, as follows: Water 2020, 12, 1812 9 of 23 A nonlinear relationship between Y i and W i was proposed [42]: where a (dimensionless) and b (mm) are two parameters. The values of S i and G i are finally estimated as: where E 0i is the potential evaporation in the ith month. In the ABCD model, the summation of the surface runoff, D i , and the base flow, F i, yields the total runoff, Q i , which can be expressed as follows:

Modified ABCD Model for Cold Regions (ABCD-CR)
The first essential condition that should be considered in cold regions is the existing of frozen soils, which highly depends on the air temperature. In general, frozen soils are low-permeability materials that can decrease the effective hydraulic conductivity of shallow soils and aquifers for the subsurface flow. Therefore, development of shallow frozen soils when the air temperature dropped below the freezing point could reduce both of groundwater recharge and discharge [60]. This mechanism was not considered in the ABCD model but the parameters c and d were related to groundwater recharge and discharge. We suggest a modification of c and d from constants to temperature-dependent parameters as follows: where α is an additional parameter (1/ • C) representing the influence of negative air temperature, T i is the basin-average monthly mean air temperature, c 0 and d 0 are normal values of c and d when the air temperature is positive. Another condition that should be considered in cold regions is the snow cover (possibly overlying glaciers). Snowmelt water (possibly including the release of glaciers) could be an additional supply to the soil water in non-frozen seasons, particularly in April and May. This was also ignored in the ABCD model. To modify the model, we divide the total precipitation, P i , into the rainfall, P i , and the snowfall, P * i , as follows: Here we did not follow the snow model supposed by Martinez and Gupta [47] who introduced two temperature thresholds and required daily minimum temperature data for input. The practicability of Equation (13) is higher because it does not include additional parameters and only require the average monthly mean air temperature. The snowfall increases the solid water storage in the snow cover and/or glaciers, S * i , as: When T i is higher than 0 • C, it is considered that a basin-scale snowmelt water depth (possibly mixed with melting glaciers), D * i , will comes into being as follows [61]: where β is also an additional parameter (1/ • C). Equation (15) is different from the snowmelt formula proposed by Martinez and Gupta [47], who assumed that the released runoff from the snowmelt is independent on temperature when the air temperature is higher than the "rain" threshold. Their assumption is not reasonable because the snowmelt rate is highly dependent on the air temperature.
In the frozen season with T i < 0, we assume that D * i = 0. A positive D* will cause a decrease in S*: Note that D * i should be restricted to S * i−1 when the result of Equation (15) is higher than S * i−1 . In this study we include the smelt water in the equivalent precipitation P i (The role of P i is the same of that in the ABCD model.) on the land surface, as follows: In addition to the previous modification, we introduce the groundwater evapotranspiration in the model. Wang et al. [45] highlighted that the groundwater evapotranspiration in zones with a small depth to water table could significantly influence the catchment water balance. In cold regions, shallow groundwater with relatively high temperature maintains unfrozen soils and can directly contribute to the land surface evapotranspiration. The groundwater evapotranspiration occurs in a small part of the catchment and roles as an additional loss of water from the normal soil water evapotranspiration, which depends on the depth to water table with a complicate process. To simplify the simulation, in this study we assume that the groundwater evapotranspiration at the basin scale, E * i , is proportional to both of the groundwater storage, G i , and the potential evaporation, E 0i , in excess of the soil water evapotranspiration, E i , as follows: where G max is the potential maximum groundwater storage (L). The modifications on the frozen soils, snow cover and groundwater evapotranspiration are incorporated into the original ABCD model to develop a new conceptual hydrological model, named ABCD-CR, for catchments in cold regions. The outline of ABCD-CR model is shown in Figure 5b. Equations (3)-(9) in the ABCD model are still applied but includes the temperature-dependent parameters defined in Equation (12). The change in groundwater storage is estimated by: Substituting Equation (17) into Equation (18), we have: The total runoff is the summation of the surface runoff, D i , the base flow, F i , and the snowmelt runoff, D * i , as follows: The ABCD-CR model has 7 parameters: a, b, c 0 , d 0 , G max , α and β.

Calibration and Performance Assessment of Models
In this study, not only the ABCD and ABCD-CR models are used. To check the efficiency of considering different processes (frozen soils, snow cover, and groundwater evapotranspiration) in this cold region, we add two models, ABCD-I and ABCD-II, in the comparison. The ABCD-I model is close to the ABCD model, and the only difference is that the ABCD-I model includes the groundwater evapotranspiration. In the ABCD-II model, the temperature-dependent behaviors are incorporated as the same in the ABCD-CR model, but the groundwater evapotranspiration is not considered. It is an equivalent model of ABCD-CR when we specify G max = ∞ in Equation (17). After preheating of the models on runoff simulation of the study catchment, we presumed the ranges of parameter values that should be satisfied in the model calibration (Table 2). Table 2. Presumed parameter ranges for model calibration and optimal results. The models were calibrated with the input monthly meteorological data in the period 1982-1985, to produce monthly runoff results matching observation data at the hydrological station W-III in the same period. Modeling performance was assessed using the Nash-Sutcliffe efficiency (NSE) [62]. Better performance is exhibited when NSE is closer to 1. For monthly simulation, the performance would be quite well when NSE > 0.7 and satisfactory when NSE > 0.5. For a test model, the NSE depends on the parameter values and was taken as the objective function in the calibration. The error of the modeling results was also evaluated with the Mean Absolute Relative Error (MARE).
The data in the period 1986-1989 were used to verify the calibrated models. This validation has to be terminated in 1990 when the hydrological station W-III was destroyed. The restored streamflow data [54] from the hydrological station W-IV was not used for the model verification due to uncertainties in estimating the leakage loss between W-III and W-IV. However, the calibrated models can be applied to produce new restored streamflow for the period since 1990 and the result is presented in Section 4.

Sensitivity Analysis on Runoff to Temperature
In the study area, the change in air temperature is a critical force that triggering changes in climatic and hydrological conditions. To check the unique impact of the change in air temperature on the change in runoff, we perform a special sensitivity analysis using the ABCD-CR model. A hypothesized scenario is used as the baseline for the period 1975-2015 when the air temperature showed a significant increasing trend. This scenario is characterized by a repeat of rising and falling temperature in every year with the same patterns from 1965 to 1974. However, the monthly precipitation and potential evaporation in this scenario are not changed from the actual data. Such a scenario is organized for the ABCD-CR model to produce a baseline of the runoff in the study area and denoted as baseline simulation. The comparison between this baseline simulation and the normal simulation with the actual air temperature is discussed in Section 5.1.

Validation and Calibration Results of the ABCD-CR Model
We calibrated and validated the ABCD-CR model using the observed runoff data, for periods of 1982-1985 and 1986-1989, respectively. The parameters were optimized by maximizing the Nash-Sutcliffe efficiency (NSE) in the calibration period. The time from January 1963 to December 1981 was used as the preheating period to minimize the influence of uncertainties in initial conditions on the model simulation. The length of this period is enough to eliminate the effect of initial conditions on the comparison of different models. The identified parameters of the ABCD-CR model are shown in Table 2. Figure 6 shows the simulated and observed monthly runoff curves at the W-III station for a comparison. The modeling results basically agree with the observation. As shown in Table 3, the NSE of the ABCD-CR model is greater than 0.7 in both of the calibration and validation periods, indicating that this model performs well. The MARE is less than 20%, also indicating a high accuracy. To check the robustness of the model calibration, an experiment with a longer validation period was also operated. In the experiment, the period 1963-1973 was chosen for model preheating and the period 1974-1977 with mean annual T of 0.99 • C was chosen as the calibration period. A 12 years period from 1978 to 1989 with mean annual T of 1.49 • C, which was warmer than that in the calibration period, was treated as the validation period. Result of the ABCD-CR model in this experiment is shown in Figure S1 in the Supplementary Material, which indicates that the simulated monthly runoff also agrees well with observations at the W-III station. The NSE was 0.70 and 0.69 for the calibration and validation periods, respectively, similar to that of the ABCD-CR model ( Table 3) using periods 1982-1985 and 1986-1989 for calibration and validation, respectively. These results illustrate that the calibrated ABCD-CR model is sufficient for long-term runoff simulation in the headwater catchment of Golmud River.
Water 2020, 12, x FOR PEER REVIEW 12 of 23 of the ABCD-CR model is greater than 0.7 in both of the calibration and validation periods, indicating that this model performs well. The MARE is less than 20%, also indicating a high accuracy. To check the robustness of the model calibration, an experiment with a longer validation period was also operated. In the experiment, the period 1963-1973 was chosen for model preheating and the period 1974-1977 with mean annual T of 0.99 °C was chosen as the calibration period. A 12 years period from 1978 to 1989 with mean annual T of 1.49 °C, which was warmer than that in the calibration period, was treated as the validation period. Result of the ABCD-CR model in this experiment is shown in Figure S1 in the Supplementary Material, which indicates that the simulated monthly runoff also agrees well with observations at the W-III station. The NSE was 0.70 and 0.69 for the calibration and validation periods, respectively, similar to that of the ABCD-CR model ( Table 3) using periods 1982-1985 and 1986-1989 for calibration and validation, respectively. These results illustrate that the calibrated ABCD-CR model is sufficient for long-term runoff simulation in the headwater catchment of Golmud River.

Comparison Among Different Models
We firstly compared the ABCD Model with the ABCD-CR model using the same methods in preheating, calibration, and validation. The calibrated parameters and simulated runoff for the ABCD model are shown in Table 2 and Figure 7a, respectively. The NSE of the ABCD model is less than 0.5 (Table 3) in the calibration period, indicating an unsatisfactory result. Although the NSE is slightly higher than 0.5 in the calibration period, the MARE is quite large (higher than 30%). The ABCD model is not successful in capture the runoff patterns, especially in the summer (Figure 7a). It significantly underestimated the peaks of runoff in summer and overestimated the low flow during the cold season.

Comparison among Different Models
We firstly compared the ABCD Model with the ABCD-CR model using the same methods in preheating, calibration, and validation. The calibrated parameters and simulated runoff for the ABCD model are shown in Table 2 and Figure 7a, respectively. The NSE of the ABCD model is less than 0.5 (Table 3) in the calibration period, indicating an unsatisfactory result. Although the NSE is slightly higher than 0.5 in the calibration period, the MARE is quite large (higher than 30%). The ABCD model is not successful in capture the runoff patterns, especially in the summer (Figure 7a). It significantly underestimated the peaks of runoff in summer and overestimated the low flow during the cold season. The simulated monthly runoff from the ABCD-I model is shown in Figure 7b. The ABCD-I model yields higher NSE values than that of the ABCD model. However, compared with the ABCD-CR model, the ABCD-I model still has lower NSE values and lager MARE ( Table 3). The errors are mainly caused by the poor performance in the cold season (October to March). As shown in Figure  7b, the observed runoff decreased from above 4 mm in October to about 2 mm in January when the temperature decreased, while the simulated runoff could not drop so much due to using constant c and d. As shown in Table 2, the b value of the ABCD-I model (less than 50 mm) is quite smaller than the b value in the ABCD model (close to 100 mm) because the limit of Gmax at 5200 mm. It allows the ABCD-I model to capture the high runoff in the summer. Figure 7c shows the simulated monthly runoff by the ABCD-II model. The NSE of the ABCD-II model is higher than the NSE values of the ABCD and ABCD-I models but still smaller than 0.7 (Table  3). For the runoff in the cold season, the ABCD-II model performs better than the ABCD and ABCD-I models because it incorporates the temperature-dependent parameters. However, it remains a large b value (100 mm) that is similar to the ABCD model and requires a highest d value among that of these models (Table 2), causing a significant change in groundwater storage. We found that the annual groundwater storage obtained by the ABCD and ABCD-II model increased rapidly and leaded to unreasonable rises in runoff during 1990-2015 (Section 4.3). This unreasonable result is caused by ignoring the groundwater evapotranspiration. In comparison, the groundwater storage simulated from the ABCD-CR model was relatively stable with a fluctuation that agrees with the variations in the annual precipitation.
Consequently, the different performances of the ABCD-I and ABCD-II models indicate that both of the groundwater evapotranspiration and temperature-dependent runoff generation are necessary in consideration for the study area. The ABCD-CR model significantly improves the runoff simulation from the original ABCD model, with a relatively low number of additional parameters. The simulated monthly runoff from the ABCD-I model is shown in Figure 7b. The ABCD-I model yields higher NSE values than that of the ABCD model. However, compared with the ABCD-CR model, the ABCD-I model still has lower NSE values and lager MARE ( Table 3). The errors are mainly caused by the poor performance in the cold season (October to March). As shown in Figure 7b, the observed runoff decreased from above 4 mm in October to about 2 mm in January when the temperature decreased, while the simulated runoff could not drop so much due to using constant c and d. As shown in Table 2, the b value of the ABCD-I model (less than 50 mm) is quite smaller than the b value in the ABCD model (close to 100 mm) because the limit of Gmax at 5200 mm. It allows the ABCD-I model to capture the high runoff in the summer. Figure 7c shows the simulated monthly runoff by the ABCD-II model. The NSE of the ABCD-II model is higher than the NSE values of the ABCD and ABCD-I models but still smaller than 0.7 (Table 3). For the runoff in the cold season, the ABCD-II model performs better than the ABCD and ABCD-I models because it incorporates the temperature-dependent parameters. However, it remains a large b value (100 mm) that is similar to the ABCD model and requires a highest d value among that of these models (Table 2), causing a significant change in groundwater storage. We found that the annual groundwater storage obtained by the ABCD and ABCD-II model increased rapidly and leaded to unreasonable rises in runoff during 1990-2015 (Section 4.3). This unreasonable result is caused by ignoring the groundwater evapotranspiration. In comparison, the groundwater storage simulated from the ABCD-CR model was relatively stable with a fluctuation that agrees with the variations in the annual precipitation.
Consequently, the different performances of the ABCD-I and ABCD-II models indicate that both of the groundwater evapotranspiration and temperature-dependent runoff generation are necessary in consideration for the study area. The ABCD-CR model significantly improves the runoff simulation from the original ABCD model, with a relatively low number of additional parameters.

Restored Runoff at W-III for the Period 1990-2015
As previously mentioned, the hydrological station W-III had been out of service since it was destroyed in 1990. As a consequence, the river discharge from the headwater catchment of the Golmud Rver has been unknown since 1990. This causes a difficulty in assessing the change in water resources for the Golmud city. However, the unknown streamflow can be restored by the hydrological models developed in this study.
From the comparison among different models, we were sure of that the ABCD-CR model can provide the most reliable result of the restored river discharge at the hydrological station W-III. We used the calibrated models introduced in Section 4.2 to simulate the monthly runoff at the hydrological station W-III for the period 1990-2015 with the same meteorological forcing data. Results are summarized in Figure 8. In the period 1990-2015, the annual P showed an increasing trend before 2010 and then fell rapidly, whereas the annual E 0 showed a slight decreasing trend (Figure 8a). The annual runoff reproduced by the ABCD and ABCD-II model show evident increasing trends and are much larger than that reproduced by ABCD-CR (Figure 8b) because of ignoring the groundwater evaporation, which seem unreasonable. It also demonstrates that the improvement of considering groundwater evaporation is necessary. In comparison, the simulated changes in the annual runoff at W-III from ABCD-I and ABCD-CR are closer to that reconstructed by Wang and Wang [54] with an estimation of the river water leakage between W-III and W-IV. It can be seen that the result in Wang and Wang [54] is generally smaller than that from ABCD-I and ABCD-CR, especially in the period 2000-2004 and after 2013. The reconstructed runoff from the ABCD-CR model follows well with the fluctuation in annual P, which seems more reasonable. It is clear that the climate change was the major force of the increase in the annual runoff in this study catchment with an extremely low population density (<3 people/km 2 , 2018 census data) and a little change in the land use.

Restored Runoff at W-III for the Period 1990-2015
As previously mentioned, the hydrological station W-III had been out of service since it was destroyed in 1990. As a consequence, the river discharge from the headwater catchment of the Golmud Rver has been unknown since 1990. This causes a difficulty in assessing the change in water resources for the Golmud city. However, the unknown streamflow can be restored by the hydrological models developed in this study.
From the comparison among different models, we were sure of that the ABCD-CR model can provide the most reliable result of the restored river discharge at the hydrological station W-III. We used the calibrated models introduced in Section 4.2 to simulate the monthly runoff at the hydrological station W-III for the period 1990-2015 with the same meteorological forcing data. Results are summarized in Figure 8. In the period 1990-2015, the annual P showed an increasing trend before 2010 and then fell rapidly, whereas the annual E0 showed a slight decreasing trend (Figure 8a). The annual runoff reproduced by the ABCD and ABCD-II model show evident increasing trends and are much larger than that reproduced by ABCD-CR (Figure 8b) because of ignoring the groundwater evaporation, which seem unreasonable. It also demonstrates that the improvement of considering groundwater evaporation is necessary. In comparison, the simulated changes in the annual runoff at W-III from ABCD-I and ABCD-CR are closer to that reconstructed by Wang and Wang [55] with an estimation of the river water leakage between W-III and W-IV. It can be seen that the result in Wang and Wang [55] is generally smaller than that from ABCD-I and ABCD-CR, especially in the period 2000-2004 and after 2013. The reconstructed runoff from the ABCD-CR model follows well with the fluctuation in annual P, which seems more reasonable. It is clear that the climate change was the major force of the increase in the annual runoff in this study catchment with an extremely low population density (<3 people/km 2 , 2018 census data) and a little change in the land use. The mean annual water balance data in the study catchment were summarized for the periods [1990][1991][1992][1993][1994][1995][1996][1997][1998][1999][2000][2001][2002][2003][2004][2005][2006][2007][2008][2009][2010][2011][2012][2013][2014][2015], and compared to that in the period from 1982 to 1989 (Table 4). From the period 1982-1989 to period 1990-2000, the mean annual T increased with 0.55 °C, whereas the mean annual P decreased with 29.7 mm and the mean annual E0 changed a little. This led to a reduction in recharge from P to G and a small change in the mean annual Q (+2.4 mm) as estimated from the ABCD-CR model, with slight releases in S * and D *. However, the ABCD model resulted in a higher increase in The mean annual water balance data in the study catchment were summarized for the periods [1990][1991][1992][1993][1994][1995][1996][1997][1998][1999][2000][2001][2002][2003][2004][2005][2006][2007][2008][2009][2010][2011][2012][2013][2014][2015], and compared to that in the period from 1982 to 1989 (Table 4). From the period 1982-1989 to period 1990-2000, the mean annual T increased with 0.55 • C, whereas the mean annual P decreased with 29.7 mm and the mean annual E 0 changed a little. This led to a reduction in recharge from P to G and a small change in the mean annual Q (+2.4 mm) as estimated from the ABCD-CR model, with slight releases in S* and D*. However, the ABCD model resulted in a higher increase in the mean annual Q (+5.18 mm). In the period 2001-2015, both the mean annual T and P increased significantly in comparison with that in the period 1982-1989, whereas the mean annual E0 decreased with 8.3 mm. It led to a substantial increase in the mean annual Q (+16.4 mm from the period 1982-1989) with increases in E and E* while decreases in S* and D*, as estimated from the ABCD-CR model. In comparison, the ABCD model overestimated the increase in the mean annual Q. This is mainly caused by the overestimation of the groundwater storage with a much higher ∆G and the underestimation of the mean annual E. Wang and Wang [54] only restored the annual streamflow at the hydrological station W-III in the period from 1990 to 2015. In this study, we restored the monthly runoff with respect to W-III using the ABCD-I model and ABCD-CR model, and their differences are compared in Figure 9. As that illustrated in Figure 9a, the average monthly T from April to September (warm season) in every year is higher than zero, followed by a rise and fall in the average monthly P. In the other months (cold season), the average monthly T is below zero and P is generally smaller than 20 mm. Both the ABCD-I model and ABCD-CR model could capture the response of runoff to the precipitation in the warm season, resulting in similar monthly runoff data. For the cold season, the runoff reconstructed from the ABCD-CR model soundly follows the change in the air temperature, but that reconstructed from the ABCD-I model shows too small variation. This result implies that the temperature-effect should be considered in a hydrological model for reconstructing the streamflow at W-III in the cold season.
Water 2020, 12, x FOR PEER REVIEW 15 of 23 the mean annual Q (+5.18 mm). In the period 2001-2015, both the mean annual T and P increased significantly in comparison with that in the period 1982-1989, whereas the mean annual E0 decreased with 8.3 mm. It led to a substantial increase in the mean annual Q (+16.4 mm from the period 1982-1989) with increases in E and E * while decreases in S * and D *, as estimated from the ABCD-CR model. In comparison, the ABCD model overestimated the increase in the mean annual Q. This is mainly caused by the overestimation of the groundwater storage with a much higher ΔG and the underestimation of the mean annual E. Wang and Wang [55] only restored the annual streamflow at the hydrological station W-III in the period from 1990 to 2015. In this study, we restored the monthly runoff with respect to W-III using the ABCD-I model and ABCD-CR model, and their differences are compared in Figure 9. As that illustrated in Figure 9a, the average monthly T from April to September (warm season) in every year is higher than zero, followed by a rise and fall in the average monthly P. In the other months (cold season), the average monthly T is below zero and P is generally smaller than 20 mm. Both the ABCD-I model and ABCD-CR model could capture the response of runoff to the precipitation in the warm season, resulting in similar monthly runoff data. For the cold season, the runoff reconstructed from the ABCD-CR model soundly follows the change in the air temperature, but that reconstructed from the ABCD-I model shows too small variation. This result implies that the temperature-effect should be considered in a hydrological model for reconstructing the streamflow at W-III in the cold season.

Effect of Temperature Changes on Runoff
The increasing air temperature may directly lead to changes in the snowmelt and frozen soils, which are worth checking for the study on changing streamflow in a cold region. In the study catchment, the annual T increased rapidly at a rate of 0.5 • C/dec during the period of 1975-2015, as shown by the normal curve in Figure 10a. After 1994, the annual T became significantly higher than the baseline extracted from the data during 1965-1974. The snowfall also displayed a slight increase from 1975 to 2015 at a rate of 0.07 mm/yr (Figure 10b). In particular, we use the method of sensitivity analysis (Section 3.4) to investigate the change in snowmelt runoff in April and May. When the increasing trend in the air temperature is replaced by the baseline, the snowmelt simulated from the ABCD-CR model will slightly increase in both April (Figure 10c) and May (Figure 10d). This rise of snowmelt runoff is a unique effect of the increasing snowfall. In comparison, the normal simulation (with the true data of T) results of the ABCD-CR model yield an increasing snowmelt in April with a higher rate (Figure 10e) than the baseline result and display a decreasing snowmelt in May (Figure 10f). Such opposite trends for the April and May snowmelt runoff were a unique effect of the increasing T. The snowmelt becomes more concentrated in April as triggered by the increasing air temperature and it leads to a reduction in the May snowmelt because the remained snow cover storage in May becomes smaller and smaller. Note that the total snowmelt runoff still exhibits an increasing trend as a result of the increasing snowfall. Thus, the period of spring flood in the Golmud River will be shorten by the air temperature rising while the peak of the flood will increase.
Water 2020, 12, x FOR PEER REVIEW 16 of 23

Effect of Temperature Changes on Runoff
The increasing air temperature may directly lead to changes in the snowmelt and frozen soils, which are worth checking for the study on changing streamflow in a cold region. In the study catchment, the annual T increased rapidly at a rate of 0.5 °C/dec during the period of 1975-2015, as shown by the normal curve in Figure 10a. After 1994, the annual T became significantly higher than the baseline extracted from the data during 1965-1974. The snowfall also displayed a slight increase from 1975 to 2015 at a rate of 0.07 mm/yr (Figure 10b). In particular, we use the method of sensitivity analysis (Section 3.4) to investigate the change in snowmelt runoff in April and May. When the increasing trend in the air temperature is replaced by the baseline, the snowmelt simulated from the ABCD-CR model will slightly increase in both April (Figure 10c) and May (Figure 10d). This rise of snowmelt runoff is a unique effect of the increasing snowfall. In comparison, the normal simulation (with the true data of T) results of the ABCD-CR model yield an increasing snowmelt in April with a higher rate (Figure 10e) than the baseline result and display a decreasing snowmelt in May ( Figure  10f). Such opposite trends for the April and May snowmelt runoff were a unique effect of the increasing T. The snowmelt becomes more concentrated in April as triggered by the increasing air temperature and it leads to a reduction in the May snowmelt because the remained snow cover storage in May becomes smaller and smaller. Note that the total snowmelt runoff still exhibits an increasing trend as a result of the increasing snowfall. Thus, the period of spring flood in the Golmud River will be shorten by the air temperature rising while the peak of the flood will increase. Landscape changes in the study area mainly resulted from frozen soil changes. The hydraulic conductivity is affected by temperature, which is controlled by soil freezing and thawing. In this study, this mechanism is described by linking air temperature to parameter c and d in the ABCD-CR model. As air temperature rises, c and d increase, which describes the increase in hydraulic conductivity due to changes in frozen soil. The trends in temperature are almost the same at the rate of +0.4 °C/dec in the warm season and cold season, as shown in Figures 11a and b. Figure 11c shows Landscape changes in the study area mainly resulted from frozen soil changes. The hydraulic conductivity is affected by temperature, which is controlled by soil freezing and thawing. In this study, this mechanism is described by linking air temperature to parameter c and d in the ABCD-CR model. As air temperature rises, c and d increase, which describes the increase in hydraulic conductivity due to changes in frozen soil. The trends in temperature are almost the same at the rate of +0.4 • C/dec in the warm season and cold season, as shown in Figure 11a,b. Figure 11c shows that the runoff in the baseline simulation almost equals to the normal simulation in the warm season. However, for the cold season the runoff in the normal simulation was larger than that in the baseline simulation, as shown in Figure 11d. This result implies that the frozen soil degradation was a positive change for the cold season runoff, leading to the increase trend in the annual cold season runoff during the period from 1975 to 2015. Differently, the impact of the frozen soil on the warm season runoff was quite limited. Thus, the change in the warm season runoff should be mainly caused by the changes in the precipitation and evapotranspiration.
Water 2020, 12, x FOR PEER REVIEW 17 of 23 that the runoff in the baseline simulation almost equals to the normal simulation in the warm season. However, for the cold season the runoff in the normal simulation was larger than that in the baseline simulation, as shown in Figure 11d. This result implies that the frozen soil degradation was a positive change for the cold season runoff, leading to the increase trend in the annual cold season runoff during the period from 1975 to 2015. Differently, the impact of the frozen soil on the warm season runoff was quite limited. Thus, the change in the warm season runoff should be mainly caused by the changes in the precipitation and evapotranspiration.

Streamflow Versus Groundwater Level
It has been considered that the change in streamflow of the Golmud River triggered the fluctuation of groundwater in the Golmud city [55] that shown in Figure 4b. To check the relationship between the stream flow and groundwater, we plotted the annual runoff in the headwater catchment of the Golmud river (estimated from the ABCD-CR model) and the observed groundwater level in G1-G4 in Figure 12.

Streamflow Versus Groundwater Level
It has been considered that the change in streamflow of the Golmud River triggered the fluctuation of groundwater in the Golmud city [54] that shown in Figure 4b. To check the relationship between the stream flow and groundwater, we plotted the annual runoff in the headwater catchment of the Golmud river (estimated from the ABCD-CR model) and the observed groundwater level in G1-G4 in Figure 12.
that the runoff in the baseline simulation almost equals to the normal simulation in the warm season. However, for the cold season the runoff in the normal simulation was larger than that in the baseline simulation, as shown in Figure 11d. This result implies that the frozen soil degradation was a positive change for the cold season runoff, leading to the increase trend in the annual cold season runoff during the period from 1975 to 2015. Differently, the impact of the frozen soil on the warm season runoff was quite limited. Thus, the change in the warm season runoff should be mainly caused by the changes in the precipitation and evapotranspiration.

Streamflow Versus Groundwater Level
It has been considered that the change in streamflow of the Golmud River triggered the fluctuation of groundwater in the Golmud city [55] that shown in Figure 4b. To check the relationship between the stream flow and groundwater, we plotted the annual runoff in the headwater catchment of the Golmud river (estimated from the ABCD-CR model) and the observed groundwater level in G1-G4 in Figure 12.  As indicated, a strong positive correlation between the streamflow and the groundwater level exists for G1-G3 (Figure 12a-c). It implies that the groundwater level in the area of the Golmud city is mainly controlled by the runoff in the study catchment, even other factors such as the local precipitation could also influence the groundwater regime. In the Golmud city, the annual P is less than the 4% of the annual E 0 (Table 1), resulting in a tiny groundwater recharge from the local precipitation in comparison to the leakage recharge from the Golmud river. More streamflow in the river brings more leakage into groundwater and causes higher groundwater level so that the groundwater level is positively related to the runoff in the headwater catchment. The groundwater utilization and management in the city could also influence the groundwater regime. However, the observation wells G1-G3 are far enough from the center of the city to avoid the significant impact from human activities. The correlation between the streamflow and the groundwater level in G4 (Figure 12d) is also positive but weak (R 2 < 0.4). G4 is located in the eastern bank of the Golmud River and close to the urban zone (Figure 4a), where the disturbance from the groundwater utilization and management is significant. Nevertheless, the positive correlation between the groundwater level and streamflow demonstrates that the runoff rise in the headwater catchment of the Golmud river was the major cause for the rising groundwater level in the area of the Golmud city over the past 20 years.
The ABCD-CR model developed in this study can role as a useful tool for predicting the streamflow changes in the upper Golmud River basin with the changing climate in the QTP. And then, the results of the ABCD-CR model can be used to assess the trend in the groundwater level near the Golmud city. This is a worthy work in the future study.

Comparison with Previous Similar Studies
The ABCD model is a widely used hydrological model due to its simplicity and efficiency. In recent years, the ABCD model was mainly applied for analyzing changes in runoff at the annual scale [44][45][46][47][48], although it was proposed at monthly scale [42]. On the other hand, the ABCD model was rarely examined in cold regions. This study shows that the difference between the annual runoff simulated by the ABCD model and the ABCD-CR model was limited (Figure 8b). However, the ABCD model has unreasonable results at monthly scale, and ABCD-CR model greatly improves the performance in runoff simulation by considering groundwater evaporation and the temperature effect on model parameter. This finding is partly consistent with that of Wang and Zhou [45] who found that groundwater evaporation must be considered in the simulation of the ABCD model when the groundwater level was high at monthly scale. Moreover, this study found that the temperature effect on the runoff should also be considered in cold regions in addition to the groundwater evaporation.
The basic finding with the ABCD-CR model in this study is consistent with some other studies in cold regions using distributed hydrological models [29][30][31]. For example, Qi et al. [30] incorporated a three-layer snowmelt model and frozen soil parameterization to the Water and Energy Budget Based Distributed Hydrological Model (WEB-DHM). They found that the modified WEB-DHM improved the performance for simulating river discharge in the source region of the Yangtze River, especially in spring and in the cold season. Gao et al. [31] simulated cryospheric processes including snow and soil freezing/thawing with the Geomorphology-Based Eco-Hydrological Model (GBEHM) for the upper Heihe River basin, northeastern QTP. They also found that an increasing air temperature will cause an increase in runoff in spring and winter, which were attributed to the increased snowmelt and increased hydraulic conductivity triggered by frozen soil degradation. These models required a large amount of observed data for model input and parameter estimation. Therefore, they are not easy to be applied in the regions with scare observations such as the headwater catchment of the Golmud River. This study considered the key processes in cold regions and incorporated them in the ABCD-CR model with a simplified way. As indicated, the temperature-dependent processes can be sufficiently captured with a few additional parameters to the original ABCD model, which is an advantage for the application of hydrological models in cold regions with limited observations.
In recent years, researchers proposed remote sensing approaches to retrieve river discharges in the QTP [25,26]. Wang et al. [63] proposed a method to retrieve the river discharge using the area of water surface estimated from remote sensing images. However, this method was mainly applied for large rivers, and was difficult to be used for the narrow river segment (width < 400 m), such as the Golmud River, due to the limitation in remote sensing images. For the narrow rivers, they argued that a rating curve need to be developed based on in-situ discharges and river width (or precipitation) estimated from satellite data [63]. Similar with hydrological models, the rating curve needs to be fitted and may have large uncertainties in parameter. On the other hand, the remote sensing approaches could only retrieve runoff data in the period after the launching of the satellite. Compared with the remote sensing approaches, the ABCD-CR model provides a more flexible and efficiency method to reconstruct the runoff data.
This study found that the runoff from the headwater catchment of the Golmud River dominated the changes in groundwater level in the middle and lower Golmud River basin. This finding is consistent with Kou [64] and Wang and Wang [54]. Kou [64] studied the interactions between surface water and groundwater in the middle and lower Golmud River basin. He found that the lateral recharge in deep aquifers could be ignored and the leakage from the Golmud River was the main source of the recharge of groundwater. He also found that groundwater utilization rate was less than 10%, indicating water extraction from human activities had limited effect on the changes in groundwater level. Wang and Wang [54] also argued that the upstream surface water was the main source of groundwater recharge.
Changes in runoff in the headwater catchment of the Golmud River were rarely studied in the previous studies. Kong et al. [65] analyzed the isotopic characteristics of the river water in the Qaidam Basin including the whole Golmud River basin and analyzed the changes in streamflow in recent 20 year. They found an increasing trend in precipitation and snowmelt, and considered climatic factor was the main factor increasing river discharge in the rivers. These findings agreed with the results of model simulation in this study. Wen et al. [66] used three trend analysis methods to investigate variations in monthly and annual runoff in the upper reaches of the Bayin River basin located near the headwater catchment of the Golmud River. They found that annual runoff showed a significant increase from 1959 to 2013 and the river experienced a period of large positive anomaly since the beginning of 21st century. This trend is similar with the changes in runoff retrieved by the ABCD-CR model in the study area.

Uncertainties and Limitations
Results in this study could be influenced by uncertainties in the observation. Because no national meteorological stations were located in the study catchment, it is difficult to fully reproduce the actual spatial distributions of the precipitation, causing uncertainties in the model input. More meteorological observations in the study area are expected to improve hydrological studies. In this study we assumed that the change in the air temperature has the same hydrologic effect as that from the change in ground temperature, which is not reasonable if the difference between the air temperature and ground temperature was significant. It could be fixed if the ABCD-CR model uses the monthly data of the ground temperature, while, this kind of observations was absent in the study area.
Additional limitations in the ABCD-CR model are linked with the simplifications. For example, the effect of frozen soil changes on runoff was simply represented by the temperature-dependent parameters, c and d. Complex soil freezing and thawing processes were not fully incorporated. It reduces the computational cost but may ignore more complicated impacts of frozen soil changes on hydrological processes. This limitation should be examined in further investigations. Equations or formulas describing the heat transfer in soil may be incorporated to the model to improve the model performance in a future study. On the other hand, the runoff simulated by the ABCD-CR model could provide boundary conditions for a groundwater model, leading to a coupling of surface and subsurface processes which could be applied to predict groundwater changes in future studies.

Conclusions
In this study, the ABCD model was modified with temperature-dependent parameters and groundwater evapotranspiration to simulate monthly runoff for catchments in a cold region. The new hydrological model, namely ABCD-CR, was used in analyzing the headwater catchment of the Golmud River located in the QTP, China, with significantly limited observations. This catchment has missed hydrological gauge data since 1990 because the damage of the hydrological station in a heavy flood. We reconstructed the monthly streamflow patterns for the catchment and analyzed the temperature effect. The following conclusions can be drawn: (1) Incorporating the temperature and groundwater-evapotranspiration effects could significantly enhance the performance of the hydrological model for catchments in cold regions. Among several comparable hydrological models, the ABCD-CR model provides the best simulation on the monthly runoff in the headwater catchment of the Golmud River observed before 1990. For the period after 1990, the ABCD model overestimated the mean annual runoff and underestimated the mean annual evapotranspiration. (2) With increasing in the air temperature during 1975-2015, the annual snowmelt runoff in the study catchment showed an increasing trend and a change in seasonal distribution. The concentrated release period of the snowmelt runoff gradually shifted from May to April. The spring flood in the Golmud River will shorten the time while raise the peak. (3) The ABCD-CR model can capture the increasing effect of the hydraulic conductivity of soils with the frozen soil degradation caused by the global warming. As revealed, the frozen soil degradation led to an increase in the cold season runoff whereas did not affect the warm season runoff. (4) The annual runoff in the headwater catchment of the Golmud River was positively correlated with the groundwater level in the area of the Golmud city. The increasing streamflow in the Golmud River was the major cause for the rising groundwater level in the Golmud city over the past 20 years.
The model developed in this study may provide a useful tool for analyzing the changes in the catchment scale runoff in cold regions triggered by the climate change.