Long-Term Projection of Water Cycle Changes over China Using RegCM

: The global water cycle is becoming more intense in a warming climate, leading to extreme rainstorms and ﬂoods. In addition, the delicate balance of precipitation, evapotranspiration, and runoff affects the variations in soil moisture, which is of vital importance to agriculture. A systematic examination of climate change impacts on these variables may help provide scientiﬁc foundations for the design of relevant adaptation and mitigation measures. In this study, long-term variations in the water cycle over China are explored using the Regional Climate Model system (RegCM) developed by the International Centre for Theoretical Physics. Model performance is validated through comparing the simulation results with remote sensing data and gridded observations. The results show that RegCM can reasonably capture the spatial and seasonal variations in three dominant variables for the water cycle (i.e., precipitation, evapotranspiration, and runoff). Long-term projections of these three variables are developed by driving RegCM with boundary conditions of the Geophysical Fluid Dynamics Laboratory Earth System Model under the Representative Concentration Pathways (RCPs). The results show that increased annual average precipitation and evapotranspiration can be found in most parts of the domain, while a smaller part of the domain is projected with increased runoff. Statistically signiﬁcant increasing trends (at a signiﬁcant level of 0.05) can be detected for annual precipitation and evapotranspiration, which are 0.02 and 0.01 mm/day per decade, respectively, under RCP4.5 and are both 0.03 mm/day per decade under RCP8.5. There is no signiﬁcant trend in future annual runoff anomalies. The variations in the three variables mainly occur in the wet season, in which precipitation and evapotranspiration increase and runoff decreases. The projected changes in precipitation minus evapotranspiration are larger than those in runoff, implying a possible decrease in soil moisture.


Introduction
The global water cycle is becoming more intense in a warming climate; increases in precipitation, evapotranspiration, and runoff can be widely observed over the world [1,2]. The resulting extreme rainstorms and heavy runoff can themselves lead to losses of life and damage to infrastructure (such as urban drainage systems), not to mention that they may also cause flood events that are of more devastating consequences [3]. On the other hand, the delicate balance of the three variables also deserves attention, as, according to the surface water budget equation dS/dt = P − E − R (where S denotes the subsurface storage of water substances, P is precipitation, E is evapotranspiration, and R is runoff) [2,4], they are closely associated with variations in soil moisture. In cases where meteorological drought occurs (i.e., long-term rainfall deficit), the interactions among the water cycle components could affect how the meteorological drought is propagated to hydrological drought (i.e.,

Methodology and Data
Dynamical downscaling of water cycle components over China is developed using RegCM, which is a regional climate model developed by the International Center for Theoretical Physics [26]. The Community Land Model version 4.5 (CLM4.5) coupling is enabled in RegCM simulations to provide an improved description of land surface processes (e.g., carbon cycle, vegetation dynamics, and river routing) [27,28]. The detailed representation of water vapor fluxes for both non-vegetated and vegetated surfaces in CLM is expected to help with the simulation of evapotranspiration. In addition, CLM is embedded with SIMTOP (simple TOPMODEL-based runoff model) [29], which can take into account the influence of topological information in runoff generation. The simulated total runoff is then routed to active ocean or marginal seas through a river transport model [30]. More details about RegCM parameterization scheme configuration can be found in Lu et al. [22].
Two rounds of RegCM hindcast simulations are conducted for model validation purposes; one of them is driven by the ERA-Interim reanalysis data developed by the European Centre for Medium-Range Weather Forecasts [31], which provide the realistic historical climate over China; the other one is driven by the historical climate scenario of the Earth System Model developed by the Geophysical Fluid Dynamics Laboratory (GFDL) [32,33], which is used to provide the baseline for projections. The baseline period is 1986 to 2005.
The RegCM performance is validated through comparisons of the annual and seasonal averages of model results with those of the observations, remote sensing data, and reconstructed data. The months included in each season are as follows: December (of the previous year), January, and February for winter; March, April, and May for spring; June, July, and August for summer; September, October, and November for autumn. Spatial correlation is employed as a quantitative metric to reflect the similarity between the annual and seasonal averaged observation/remote sensing data/reconstructed data and the simulation results.
In this study, version 4 of the high-resolution gridded observation dataset from the Climate Research Unit (CRU) [34] is used for verifying the model-generated temperature and precipitation. This dataset is generated through the interpolation of extensive networks of gauge station observations into a 0.5 • regular grid [34], and is widely applied for the calibration/validation of global and regional climate models [35].
For evapotranspiration, two sets of remote sensing data are employed, specifically, version 6 of the Resolution Imaging Spectroradiometer (MODIS) terrestrial evapotranspiration product (MOD16A2GF v006) [36], and the latest version of the Global Land Evapotranspiration Amsterdam Model (GLEAM, v3.5b) [37,38]. The MOD16A2GF dataset is created from remotely sensed data products of MODIS (e.g., land cover and albedo) based on the Penman-Monteith equation [39]. On the other hand, GLEAM assembles various satellite-based observations (e.g., radiation from GERES, precipitation from TMPA and MSWEP, air temperature from AIRS, soil moisture from SMOS and ERA CCI SM) and derives global evaporation variables with the Priestley and Taylor model [37]. Both datasets were demonstrated to be able to reasonably represent the actual evapotranspiration over China [40,41]. The evapotranspiration from MOD16A2GF and the actual evaporation from GLEAM are used for validating the RegCM-generated evapotranspiration. It is worth noting that the start years of these two remote sensing datasets are 2000 and 2003, respectively; therefore, RegCM results averaged over the same periods (i.e., 2000 to 2005 and 2003 to 2005) are used for comparison. To assist the validation for evapotranspiration over the entire baseline period, station-based observed tank evaporation from the National Meteorological Information Center of China (NIMC) is also used (data available at data.cma.cn; accessed on 9 April 2019). The locations of the stations are shown in Figure S1 in the Supplementary Materials.
For runoff, the validation was undertaken through the comparison of the model outputs with the Global Runoff Reconstruction (GRUN), which is constructed through machine learning techniques based on runoff and meteorological observations [42]. This dataset is widely used in weather and climate research and is shown to have a reasonable performance over China [43][44][45].
For future climate, GFDL projections of two representative concentration pathways (RCPs) are employed, which are RCP4.5 and RCP8.5, respectively, for intermediate and heavy emissions. Simulations are conducted for the entire twenty-first century. Three time-slices are considered in result analysis: 2020 to 2039 (or 2030s), 2040 to 2069 (or 2050s), and 2070 to 2099 (or 2080s); time averages and trends are calculated with respect to these periods. In addition, since the water cycle components are closely related to atmospheric moisture contents, and the saturation vapor pressure is related to the temperature following the Clausius-Clapeyron equation [46], different warming periods are also considered in this study. The warming periods are defined as twenty-year periods in which the domain average temperature increases by 1, 1.5, 2, 3, and 4 • C compared with the baseline average. Table 1 lists the respective periods for each warming level under the two emission scenarios.
The numerical values of future trends are obtained based on Sen's slope estimator, and their statistical significance is examined by Mann-Kendall tests [47][48][49][50].

Model Validation
Validation results for near-surface temperature, precipitation, evapotranspiration, and runoff are shown in Figures 1-4, respectively. The columns of each figure are for different datasets (i.e., CRU, MODIS, GLEAM, GRUN, RegCM driven by ERA-Interim, and RegCM driven by GFDL); the rows are for annual and seasonal averages. It can be observed that RegCM can reasonably reproduce the spatial distribution and seasonal variations of temperature over China. As shown by the CRU observations (Figure 1a), above-zero baseline-average temperature can be found in most parts of the domain, except for the Tibetan Plateau and a small part of northeastern China. This spatial feature is realistically generated in the two sets of RegCM results, although with underestimations of various degrees (Figure 1b,c). Such underestimations, as discussed by Lu et al. [22], are partly caused by the model setup and partly due to the driving GCM data. Temperature over China demonstrates clear seasonality, i.e., hot summer, cold winter, and mild spring and autumn. RegCM is able to capture the seasonal variations, although underestimations can still be observed. The spatial correlation between the observed and modeled temperature can be found in Table S1 in the Supplementary Materials. The high correlation (ranges between 0.94 and 0.98) indicates RegCM's good performance in temperature. In addition, the results of RegCM driven by GFDL have higher correlations with the observations than those of the raw GCM data (please refer to Table S2 in the Supplementary Materials), which suggests that RegCM is able to correct some biases in GFDL's temperature results. The performance of RegCM with respect to precipitation is less satisfactory than fo temperature. Although it is able to generate the observed wet-in-the-southeast and dry in-the-northwest precipitation pattern, underestimations can be found over the southeas ern part of the domain and overestimations over the northwestern part (Figure 2a-c). Thi over-and underestimation pair was shown to be related to the simulation bias in vapo pressure by Lu et al. [22], who further argued that the bias in vapor pressure could b associated with the bias in temperature. There is also an apparent dry bias near the S chuan Basin and a wet bias near the southeastern edge of the Tibetan Plateau, which coul be related to the configuration of the cumulus convective scheme [22]. The spatial corre lation (Table S1) is lower for precipitation than for temperature, which is consistent wit The performance of RegCM with respect to precipitation is less satisfactory than for temperature. Although it is able to generate the observed wet-in-the-southeast and dry-inthe-northwest precipitation pattern, underestimations can be found over the southeastern part of the domain and overestimations over the northwestern part (Figure 2a-c). This overand underestimation pair was shown to be related to the simulation bias in vapor pressure by Lu et al. [22], who further argued that the bias in vapor pressure could be associated with the bias in temperature. There is also an apparent dry bias near the Sichuan Basin and a wet bias near the southeastern edge of the Tibetan Plateau, which could be related to the configuration of the cumulus convective scheme [22]. The spatial correlation (Table S1) is lower for precipitation than for temperature, which is consistent with the above results. The seasonal precipitation over China shows clear monsoon features (more precipitation in summer and less in winter), which is shared by the two sets of RegCM results.
Remote Sens. 2021, 13, 3832 6 the above results. The seasonal precipitation over China shows clear monsoon featu (more precipitation in summer and less in winter), which is shared by the two set RegCM results. The spatial pattern of actual evapotranspiration from MODIS and GLEAM (Fig  3a,d) shows more regional details than the observed precipitation pattern from C which is in part due to the higher resolution of remote sensing data than that of the g reasonable performance for RegCM in evapotranspiration. The spatial correlation between RegCM results and the observed tank evaporation from NIMC is 0.50 and 0.42, respectively, for the two rounds of hindcast simulations. The lower correlation between RegCM and NIMC could be related to the difference between tank evaporation and actual evapotranspiration. In terms of the seasonal variations, RegCM demonstrates overestimations in spring and winter, and underestimations in summer and autumn. The RegCM-generated spatial patterns for runoff show considerably larger biases than for other variables. Apparent overestimations can be spotted near the southeastern The spatial pattern of actual evapotranspiration from MODIS and GLEAM (Figure 3a,d) shows more regional details than the observed precipitation pattern from CRU, which is in part due to the higher resolution of remote sensing data than that of the gridded observation, and in part due to evapotranspiration's closer relationship with the geophysical characteristics of the domain than precipitation. The two sets of remote sensing data show similar annual patterns; subtle differences exist sporadically over the domain, which could be explained by the different skills of the two datasets over different land-use types [51]. RegCM shows better skills in simulating evapotranspiration than precipitation (Figure 3b,c,e,f), although minor overestimations can be noticed. As shown in Table S1, for annual average evapotranspiration, the spatial correlations between RegCM results and GLEAM are higher than 0.8, and those for MODIS are higher than 0.6, indicating a reasonable performance for RegCM in evapotranspiration. The spatial correlation between RegCM results and the observed tank evaporation from NIMC is 0.50 and 0.42, respectively, for the two rounds of hindcast simulations. The lower correlation between RegCM and NIMC could be related to the difference between tank evaporation and actual evapotranspiration. In terms of the seasonal variations, RegCM demonstrates overestimations in spring and winter, and underestimations in summer and autumn.
entire Tibetan Plateau, where overestimations in evapotranspiration can be identified; the latter is likely to be the cause of the former. The spatial correlations between the annual average runoff of the two sets of RegCM results and GRUN are 0.67 and 0.66, respectively (Table S1). From the seasonal perspective, the GRUN reconstructed runoff is high in summer and low in winter. This monsoon feature is well captured by RegCM. As shown by the spatial correlation, for runoff, RegCM shows better skills in summer/autumn (ranges between 0.65 and 0.69) than in winter/spring (between 0.37 and 0.53).  The RegCM-generated spatial patterns for runoff show considerably larger biases than for other variables. Apparent overestimations can be spotted near the southeastern corner of the Tibetan Plateau, which is likely to be related to the wet bias in precipitation that occurs at the same location. Slight underestimations in runoff can be found over the entire Tibetan Plateau, where overestimations in evapotranspiration can be identified; the latter is likely to be the cause of the former. The spatial correlations between the annual average runoff of the two sets of RegCM results and GRUN are 0.67 and 0.66, respectively (Table S1). From the seasonal perspective, the GRUN reconstructed runoff is high in summer and Remote Sens. 2021, 13, 3832 9 of 21 low in winter. This monsoon feature is well captured by RegCM. As shown by the spatial correlation, for runoff, RegCM shows better skills in summer/autumn (ranges between 0.65 and 0.69) than in winter/spring (between 0.37 and 0.53).
The RegCM-simulated domain average annual cycles for precipitation, evapotranspiration, and runoff are also examined (shown in Figure 5). The two sets of RegCM results present similar features (please refer to Figure S2 in the Supplementary Materials for a direct comparison of the two sets of simulation results). All three variables show a peak in their annual cycles during the monsoon period (May to September), which is consistent with previous results. On domain average, a considerable part of the precipitation is balanced by evapotranspiration, and a smaller portion is attributed to the runoff. The difference between the simulated precipitation and evapotranspiration is also plotted in Figure 5. As indicated by the surface water budget equation, the amount of precipitation that is not balanced by the other two variables contributes to the moisture storage in soil. It can be observed that the difference between precipitation and evapotranspiration is larger than runoff in early spring, and this relationship reverses in autumn, which indicates water storage in spring and dissipation in autumn.
Remote Sens. 2021, 13, 3832 9 of 22 The RegCM-simulated domain average annual cycles for precipitation, evapotranspiration, and runoff are also examined (shown in Figure 5). The two sets of RegCM results present similar features (please refer to Figure S2 in the Supplementary Materials for a direct comparison of the two sets of simulation results). All three variables show a peak in their annual cycles during the monsoon period (May to September), which is consistent with previous results. On domain average, a considerable part of the precipitation is balanced by evapotranspiration, and a smaller portion is attributed to the runoff. The difference between the simulated precipitation and evapotranspiration is also plotted in Figure  5. As indicated by the surface water budget equation, the amount of precipitation that is not balanced by the other two variables contributes to the moisture storage in soil. It can be observed that the difference between precipitation and evapotranspiration is larger than runoff in early spring, and this relationship reverses in autumn, which indicates water storage in spring and dissipation in autumn.

Precipitation
Having reasonable skills in reproducing the historical climate over China, RegCM is subsequently used to project future changes in water cycle components. The projected changes in precipitation over China in three future periods and under two emission scenarios are shown in Figure 6. The left three columns are for RCP4.5 and the right three for RCP8.5; for each scenario, the three columns, respectively, indicate 2030s, 2050s, and 2080s. The rows in the figure are the annual and seasonal averages. On annual average, increases in precipitation can be found in most parts of the domain. Precipitation changes of the two RCPs show certain similarities. For example, in the 2080s, precipitation increases of larger than 0.3 mm/day are expected in parts of the Tibetan Plateau, Yellow River Basin, Haihe River Basins, Yangtze Plain, and southern coastal hilly regions under both scenarios. In general, the area experiencing increased precipitation is larger under RCP8.5, especially in the Tibetan Plateau, where more pronounced increases (over 0.9 mm/day in the southeastern corner) can be observed as well.   Table 2). Some seasonal trends are also statistically significant; summer precipitation shows increasing trends of 0.02 and 0.05 mm/day per decade under the two scenarios, and autumn precipitation increases at a rate of 0.02 mm/day per decade under RCP8.5. Precipitation changes demonstrate apparent seasonal variations. In winter, precipitation decreases of over 0.3 mm/day can be found in parts of the Yunnan-Guizhou Plateau in the 2050s and 2080s under RCP4.5 and in the 2030s and 2050s under RCP8.5. In the 2080s under RCP8.5, more severe changes over larger areas can be noticed; parts of the Pearl River Basin are projected with precipitation decreases of over 0.3 mm/day and parts of the Yunnan-Guizhou Plateau of over 0.6 mm/day. Summer precipitation exhibits similar pat-terns of changes in the 2050s under the two scenarios, where decreases of over 0.3 mm/day are to be found in the middle and lower reaches of the Yangtze River Basin, and increases of over 0.9 mm/day are expected near the Hengduan Mountains located in the southeastern corner of the Tibetan Plateau. In the 2080s, the spatial distributions of summer precipitation change under the two scenarios are quite different. Under RCP4.5, precipitation decreases of over 0.3 mm/day are likely to occur in the very north of northeastern China, parts of the Yangtze Plain, and parts of the Pearl River Basin, while increases of over 0.9 mm/day are projected in the Hengduan Mountains and the southeastern coastal hilly regions. In comparison, under RCP8.5, decreases in summer precipitation mainly occur in the middle and lower reaches of the Yangtze River Basin (over 0.3 mm/day), while increases of over 0.9 mm/day are to be found in the southern parts of the Tibetan Plateau, Hengduan Mountains, and parts of the Haihe River Basin. The changes in summer precipitation may be related to variations in its major moisture transport branches, which are the transportations by the Indian monsoon, Southeast Asian monsoon, and midlatitude westerlies, as shown by Simmonds [52]. The increases in summer precipitation in southeastern China and decreases in central south and southeastern China under both scenarios could indicate an enhanced Indian summer monsoon and a subsided Southeast Asian summer monsoon. In spring and autumn, precipitation increases are to be seen in the Yangtze Plain (can reach over 1.8 mm/day) and the Yellow River Basin (over 0.9 mm/day), respectively.
The annual series of domain average precipitation under both RCPs are shown in Figure 7. For both scenarios, precipitation demonstrates an evident increasing trend (although a decreasing trend can be observed between 2050 and 2060 under RCP4.5). The Mann-Kendall test confirms the statistical significance of the trends at an α level of 0.05 (when the period of 2010 to 2100 is considered as a whole). The magnitudes of trends, given by Sen's slope estimator, are 0.02 and 0.03 mm/day per decade under RCP4.5 and RCP8.5, respectively (as shown in Table 2). Some seasonal trends are also statistically significant; summer precipitation shows increasing trends of 0.02 and 0.05 mm/day per decade under the two scenarios, and autumn precipitation increases at a rate of 0.02 mm/day per decade under RCP8.5.     Precipitation changes (with respect to the baseline period) in the three future periods under both scenarios are listed in Table 3. Under RCP4.5, the annual average precipitation is projected to increase by 0.06, 0.08, and 0.16 mm/day, respectively, in the three future periods (statistically significant at an α level of 0.05). Under RCP8.5, precipitation change in the 2030s is not statistically significant, and the increases in the 2050s and 2080s are 0.12 and 0.2 mm/day, respectively. Spring and summer precipitation is expected to undergo larger increases than that in winter and autumn. These numbers are consistent with the changes in the annual cycles of precipitation as shown in Figure 8, in which large precipitation increases can be found from April to September in the 2050s and 2080s under both scenarios.  Precipitation changes with respect to different warming levels under the two scenarios are shown in Table 4. The magnitudes of changes under the same warming level are similar under different scenarios. For example, with a domain average warming of 2 • C, the annual average precipitation is likely to increase by 0.12 and 0.13 mm/day under RCP4.5 and RCP8.5, respectively. This phenomenon is reasonable because, according to the Clausius-Clapeyron equation, the increase in the water holding capacity of the atmosphere is the same given the same temperature increase. The change in precipitation is not necessarily the same since the actual amount of water vapor available can be different. At a warming level of 2 • C, spring and summer precipitation is also projected to increase by 0.15 and 0.17 mm/day under RCP4.5, and by 0.21 and 0.17 mm/day under RCP8.5. When domain average warming reaches 4 • C under RCP8.5, the annual, spring, and summer precipitation are projected to increase by 0.19, 0.24, and 0.27 mm/day, respectively. Table 4. Projected changes (mm/day) in precipitation, evapotranspiration, and runoff at different warming levels (P, E, and R denote precipitation, evapotranspiration, and runoff, respectively).

Evapotranspiration
The projected changes in evapotranspiration over China are shown in Figure 9. Similar to the spatial pattern of precipitation changes, the annual average evapotranspiration is likely to increase over most parts of the domain. The area experiencing increased evapotranspiration is larger than that for precipitation, but the magnitude of the increase is smaller when compared with precipitation. Annual average evapotranspiration changes are within ±0.3 mm/day for most of the time under both scenarios, except that increases of over 0.3 mm/day are projected in the Hengduan Mountains, Yunnan-Guizhou Plateau, and southeastern coastal hilly regions in the 2080s under RCP8.5. Intra-annual variations can also be observed for evapotranspiration changes. The most pronounced changes are to occur in summer, in which evapotranspiration increases of over 0.3 mm/day can be found over the entire domain except parts in northern and northeastern China in the 2080s under RCP8.5, and increases of over 0.6 mm/day are to be seen in the Hengduan Mountains. In spring, evaporation increases of over 0.3 mm/day are projected in areas between the Yangtze River and Pearl River Basins in the 2080s under RCP4.5, and in most of the southern parts of the domain in the same period under RCP8.5. In autumn, evapotranspiration increases of over 0.3 mm/day can be found in the Hengduan Mountains. Evapotranspiration changes in winter are within ±0.3 mm/day for all future periods under both scenarios. The annual average evapotranspiration time series is shown in Figure 7, which appears to be less fluctuating than those of the other two variables. Evident increasing trends can be observed under both scenarios, the magnitudes of which are 0.01 and 0.03 mm/day per decade (statistically significant at an  level of 0.05), respectively ( Table 2). All seasonal trends are exclusively statistically significant, which are 0.01, 0.02, 0.01, and 0.01 mm/day per decade, respectively, for winter, spring, summer, and autumn under RCP4.5 and 0.01, 0.03, 0.04, and 0.02 mm/day per decade under RCP8.5. Annual and seasonal evapotranspiration changes in the three future periods under both scenarios are all statistically significant, except for winter in the 2030s under RCP4.5 (Table 3). Under RCP4.5, annual average evapotranspiration is projected to increase by 0.06, 0.1, and 0.12 mm/day, respectively, for the 2030s, 2050s, and 2080s, while, under RCP8.5, the increases are 0.08, 0.14, and 0.22 mm/day. Such increases in evapotranspiration can also be observed in Figure 8, in which evapotranspiration changes are always above zero and larger increases can be found in the monsoon months. At a domain average warming of 2 °C, annual average evapotranspiration is projected to increase by 0.12 and 0.14 mm/day under the two scenarios (Table 4), which are close to the amount of precipitation increases. For the 4 °C warming period under RCP8.5, annual evapotranspiration is likely to increase by 0.21 mm/day, and seasonal evapotranspiration by 0.08, 0.23, 0.36, and 0.18 mm/day for winter, spring, summer, and autumn, respectively. The annual average evapotranspiration time series is shown in Figure 7, which appears to be less fluctuating than those of the other two variables. Evident increasing trends can be observed under both scenarios, the magnitudes of which are 0.01 and 0.03 mm/day per decade (statistically significant at an α level of 0.05), respectively ( Table 2). All seasonal trends are exclusively statistically significant, which are 0.01, 0.02, 0.01, and 0.01 mm/day per decade, respectively, for winter, spring, summer, and autumn under RCP4.5 and 0.01, 0.03, 0.04, and 0.02 mm/day per decade under RCP8.5. Annual and seasonal evapotranspiration changes in the three future periods under both scenarios are all statistically significant, except for winter in the 2030s under RCP4.5 (Table 3). Under RCP4.5, annual average evapotranspiration is projected to increase by 0.06, 0.1, and 0.12 mm/day, respectively, for the 2030s, 2050s, and 2080s, while, under RCP8.5, the increases are 0.08, 0.14, and 0.22 mm/day. Such increases in evapotranspiration can also be observed in Figure 8, in which evapotranspiration changes are always above zero and larger increases can be found in the monsoon months. At a domain average warming of 2 • C, annual average evapotranspiration is projected to increase by 0.12 and 0.14 mm/day under the two scenarios (Table 4), which are close to the amount of precipitation increases. For the 4 • C warming period under RCP8.5, annual evapotranspiration is likely to increase by 0.21 mm/day, and seasonal evapotranspiration by 0.08, 0.23, 0.36, and 0.18 mm/day for winter, spring, summer, and autumn, respectively. Figure 10 shows the projected changes in future runoff. In general, the annual and seasonal variations in runoff share a certain resemblance with those in precipitation, and only the area experiencing increased runoff is considerably smaller. Runoff changes in most parts of the domain are within ±0.3 mm/day in the three future periods under both scenarios. Under RCP4.5, runoff increases of over 0.3 mm/day can be found in the Yangtze Plain in the 2080s. Under RCP8.5, over 0.3 mm/day decreases in runoff are likely to occur in parts of the Yunnan-Guizhou Plateau in the 2050s and also in parts of the Yangtze River Basin in the 2080s. Seasonal changes in future runoff also demonstrate distinct characteristics. Among all seasons, runoff changes in summer are most severe. For example, under RCP4.5, runoff increases of over 1.2 mm/day are projected in the southeastern coastal hilly regions in the 2080s, which can be related to the increased precipitation in this area. Under RCP8.5, a considerable part of the Yangtze River Basin is projected to experience a runoff reduction of over 0.6 mm/day in the 2080s, which is likely to be caused by the simultaneous decrease in precipitation and increase in evapotranspiration. In winter, changes in runoff are within ±0.3 mm/day in most parts of the domain, with decreases in the central and southwestern parts and increases elsewhere. Spring runoff is projected to increase in the Yangtze Plain by at least over 0.3 mm/day in the future. In autumn, parts of the Yellow River and Haihe River Basins are to receive runoff increases of over 0.3 mm/day. The observed changes in runoff are closely related to the corresponding changes in precipitation.

Runoff
As shown in Figure 7, the interannual variation in runoff is highly correlated with those in precipitation and precipitation minus evapotranspiration (hereafter as P − E). In addition, the changes in P − E are slightly larger than those in runoff, which are consistent with the observations of Zhang et al. [2]. They offered two possible explanations: problems with model spinup or water balance closure, and changes in terrestrial water storage as a result of global warming [2]. The latter case indicates a reduction in future soil moisture, which may bring a negative influence on agriculture. The time series of annual average runoff do not exhibit apparent trends, which is consistent with the lack of a significant trend as shown in Table 2. Previous studies also noticed that the annual runoff series does not demonstrate a trend in the historical period; for example, no statistically significant trend can be detected in the Huaihe River Basin, according to Yu et al. [53]. Significant weak trends of 0.01 mm/day per decade are projected in autumn under RCP4.5 and in winter and autumn under RCP8.5. In comparison with precipitation and evapotranspiration, there is no statistically significant change in annual average runoff except for the 2080s under RCP4.5 (0.05 mm/day). Statistically significant runoff reductions are projected in summer, which are −0.08 and −0.09 mm/day in the 2030s and 2050s under RCP4.5 and −0.12, −0.13, and −0.16 mm/day in the three future periods under RCP8.5. Such a runoff reduction in summer is also evident in Figure 8. These observations are partially consistent with those of Zhang et al. [2], who noticed that changes in precipitation, evapotranspiration, runoff, and P − E mainly occur in the wet season with simultaneous increases in all variables. In this study, the most pronounced changes also occur in the wet season. The projected reduction in runoff and P − E could be related to the difference in the domain selection between this and Zhang et al.'s studies. There is no significant change in annual average runoff under all warming levels and emission scenarios, except that an increase of 0.04 mm/day is projected at a warming level of 3 • C under RCP8.5. At a warming level of 2 • C, runoff changes of 0.02 and 0.08 mm/day are likely to occur in winter and autumn under RCP4.5, and 0.03, −0.11, and 0.08 mm/day in winter, summer, and autumn under RCP8.5. With a domain average warming of 4 • C, runoff in winter, summer, and autumn is to change by 0.03, −0.15, and 0.01 mm/day under RCP8.5.
regions in the 2080s, which can be related to the increased precipitation in this area. Under RCP8.5, a considerable part of the Yangtze River Basin is projected to experience a runoff reduction of over 0.6 mm/day in the 2080s, which is likely to be caused by the simultaneous decrease in precipitation and increase in evapotranspiration. In winter, changes in runoff are within ±0.3 mm/day in most parts of the domain, with decreases in the central and southwestern parts and increases elsewhere. Spring runoff is projected to increase in the Yangtze Plain by at least over 0.3 mm/day in the future. In autumn, parts of the Yellow River and Haihe River Basins are to receive runoff increases of over 0.3 mm/day. The observed changes in runoff are closely related to the corresponding changes in precipitation.

Discussions and Conclusions
In this study, the long-term variations in the water cycle over China are studied using RegCM. The performance of RegCM in terms of temperature, precipitation, evapotranspiration, and runoff is validated through comparisons of the model-generated results with gridded observations, remote sensing data, and reconstructed data. The results show that RegCM can reasonably capture the spatial and seasonal variations in these variables, although certain biases exist, such as a cold bias in the entire domain, a dry and wet bias pair in the southeastern and northwestern parts of the domain, some over-and underestimations of evapotranspiration, respectively, in winter/spring and summer/autumn, and some over-and underestimations of runoff near the Tibetan Plateau.
Long-term projections of precipitation, evapotranspiration, and runoff under two emission scenarios are then developed. The results show that increased annual average precipitation and evapotranspiration can be found in most parts of the domain, while a smaller part of the domain is projected with increased runoff.
For precipitation, the regions most affected by global warming are the Yangtze Plain, Yellow River and Haihe River Basins, and southeastern parts of the Tibetan Plateau, where over 0.3 mm/day increases are expected in the 2080s under both scenarios. The projected increase in precipitation in the Yellow River and Haihe River Basins can also be observed in a CMIP6 GCM ensemble according to Tian et al. [54]. It is worth noting that although the projected precipitation increase in the Tarim Basin is within 0.3 mm/day, the percentage change could be large considering its low annual average total precipitation, which is why several studies (e.g., [55,56]) identified it as an area vulnerable to climate change. Precipitation increase has been shown to be among the driving factors for the increased flood frequency in the Tarim River Basin since the 1980s [57]; the increase in future precipitation as projected by this and the previous studies may indicate increased flood risks in this area, which suggests the need for relevant flood prevention measures.
For evapotranspiration, areas experiencing evident increases are the Hengduan Mountains, Yunnan-Guizhou Plateau, and southeastern coastal hilly regions; the magnitude of change is 0.3 mm/day in the 2080s under RCP8.5. The apparent evapotranspiration increase in southeastern coastal hilly regions is also noted by Su et al. [58]. In terms of seasonal variations, summer and spring are likely to see larger increases in evapotranspiration; this feature is consistent with Ma et al.'s observation [59]. The developed projections in evapotranspiration can be used to evaluate climate change impacts on drought conditions through the calculation of evapotranspiration deficit; this index, compared with those that are based on precipitation and soil moisture, could more effectively reflect moisture deficiency in ecosystems [60].
For runoff, the regions most affected are the Yangtze Delta, Yunnan-Guizhou Plateau, and parts of the Yangtze River Basin, with increases of 0.3 mm/day for the first region in the 2080s under RCP4.5 and decreases of over 0.3 mm/day for the latter two regions in the 2080s under RCP8.5. The projected decrease in runoff in the middle reaches of the Yangtze River Basin is consistent with the results from Xing et al.'s study [61]. In addition, the runoff reduction in the Yunnan-Guizhou Plateau and the upper reaches of the Yangtze River Basin is also noticed by Zhai et al. in their ensemble projection of runoff [62]. Extreme high and low runoff are often related to flood and drought hazards [63], and the runoff projection developed in this study could help identify regions vulnerable to increased flood and drought risks, and thus support flood mitigation and water resource management [64].
In summary, future precipitation and evapotranspiration are likely to increase over China in the wet season, while runoff decreases. The projected changes in precipitation minus evapotranspiration are larger than those in runoff, implying a possible decrease in soil moisture. It is important that future variations in the water cycle components be considered when designing flood and drought mitigation measures.
Extensions of this study can be conducted with respect to the current limitations. For example, more sophisticated bias correction techniques can be applied to the model results. In addition, more GCMs can be used to drive RegCM so that an ensemble can be constructed for more robust projections.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/rs13193832/s1, Figure S1: Location of NIMC weather stations. Note: Red dots indicate selected stations for evapotranspiration data. The NIMC tank evaporation data contains a large number of missing data; the stations are selected if the percentage of missing data is less than 50% for each season. Figure S2: Simulated annual cycles for precipitation, evapotranspiration, runoff, and the difference between precipitation and evapotranspiration in the baseline period. The R2's for the two sets of RegCM results are 0.88, 0.92, 0.83, and 0.82, respectively for precipitation, evapotranspiration, runoff, and the difference between precipitation and evapotranspiration. Table S1: Spatial correlation between RegCM-generated results and observations, remote sensing data, and reconstructed data. Table S2: Spatial correlation between raw GFDL data and observations.