Attributions of Evapotranspiration and Gross Primary Production Changes in Semi-Arid Region: A Case Study in the Water Source Area of the Xiong’an New Area in North China

Investigating the attributions of evapotranspiration (ET) and gross primary production (GPP) changes is of great importance for regional, sustainable water resources and ecological management in semi-arid regions. Based on the simulation conducted during 2000–2019 by improving water-carbon coupling Distributed Time Variant Gain Model, the trends of ET and GPP were estimated and the driving factors were identified via 10 experimental scenarios in the water source area of the Xiong’an New Area in North China. The results show significant increases both in ET and GPP by 2.4 mm/a and 6.0 gC/m2/a in the region, respectively. At the annual scale, increasing precipitation dominates the ET uptrend. Air temperature, humidity and the interactive effects also contribute to the ET uptrend, and the contributions are 12.8%, 2.0% and 2.3%, respectively, while elevated atmospheric CO2 concentration (eCO2) and solar dimming lead to ET changes of about −7.2% and −12.4%, respectively. For the GPP changes, the increase in GPP is mainly caused by eCO2, increasing precipitation and rising temperature with the contributions of 56.7%, 34.8% and 27.8%, respectively. Solar dimming, humidity and windspeed contribute −6.8%, −4.8% and −3.5% of the GPP changes. Compared to climate change, land use and cover change has smaller effects on both ET and GPP for the few changes in land coverage. At the seasonal scale, ET and GPP increase to a greater extent during the growing season in spring and summer than in autumn and winter. Precipitation, temperature and eCO2 are generally the main causes for ET and GPP changes. Meanwhile, the decreasing humidity and rising temperature are dominant factors for ET and GPP increases, respectively, in winter. Furthermore, solar dimming has strong effects on ET reduction in autumn. The contribution of the interactive effects is much higher on a seasonal scale than annual scale, contributing to considerable decreases in ET and GPP in spring, increases in ET in autumn and winter, and an increase in GPP in winter. This study highlights the importance of considering water-carbon coupling on the attributions of ET and GPP changes and the differentiation of the effects by the abovementioned influential factors at annual and seasonal scales.


Introduction
Terrestrial evapotranspiration (ET) and gross primary production (GPP) are the key variables in global water and carbon cycles, respectively [1]. ET has a strong effect on catchment runoff, soil water content and atmospheric moisture feedback [2][3][4]. GPP, the total amount of carbon dioxide fixed by vegetation through photosynthesis, is of mainly caused by climate change [33]. However, to our knowledge, there are still only a few studies conducted on the attributions of both ET and GPP changes considering multiple factors in this region.
In this study, the improved coupling water-carbon DTVGM ecohydrological model was used to simulate the ET and GPP in the six catchments of the water source area of the Xiong'an New Area and then the contributions of multiple influence factors including climate variables, eCO 2 and underlying land cover change and their interactive effects on ET and GPP were quantified at annual and monthly timescales based on ten experimental scenarios. The specific objectives of this paper include: (1) evaluating the model applicability in the study region; (2) estimating the changes of ET and GPP in the catchments; and (3) investigating the dominant factors for the ET and GPP changes both at annual and monthly timescales.

Study Area
The water-source region of the Xiong'an New Area is located between 38 • 3 ~40 • 4 N and 113 • 39 ~116 • 12 E in the semi-arid region in North China, and covers an area of 3.12 × 10 4 km 2 with the mountainous areas in the northwest and plain area in the southeast ( Figure 1). Because the mountainous areas contribute to more than 75% of the water-yield generation for the Xiong'an New Area, the major six catchments of the water-source region were selected as the study area, including the Baigouhe Catchment (BGC), Jumahe Catchment (JMC), Baishahe Catchment (BSC), Tanghe Catchment (THC), Dashahe Catchment (DSC) and Cihe Catchment (CHC), ordered from the northeast to southwest. Details of the six catchments are shown in Table 1. The region has a semi-arid continental monsoon climate, with a mean annual temperature of 8.4 • C and it increases gradually from southwest to northeast across the mountains. The annual precipitation is about 576 mm and is mainly concentrated in the summer season. The major land-cover types in the region are forest (42.8%), grassland (39.9%), cropland (12.0%) and the rest of the land types, including water area and urban land, are relatively small, only accounting for about 5.3%. The dominant land-use type in the northeastern part is forest (83%) in the BGC and the area of forest decreases gradually to 22% in the THC and then increases again to 47% in the southwestern part (CHC), while the grassland and cropland increases firstly and then decreases from northeast to southwest. The forest is a mix of temperate coniferous and broad-leaved, mainly composing of oak, hornbeam, Chinese pie, larch and poplar. The grassland consists of white grass, yellow back grass and artemisia, etc. The crops include grains, beans, tubers and oil crops, etc.

Dataset
In order to establish the ecohydrological DTVGM model, several datasets including the topographic and underlying surface data, climate forcing data and remotely-sensed based datasets were used. The topographic and underlying surface data were collected, including the SRTM 90 m Digital Elevation Database, the homogenized soil data at 1000 m resolution from FAO and land-use coverage data at 1000 m resolution from Resource and Environment Science and Data Center. The meteorological data includes the daily atmospheric pressure, precipitation, maximum and minimum temperature, windspeed, relative humidity and sunshine duration in seven meteorological stations obtained from the China Meteorological Data Sharing Service System and daily precipitation data in 28 rainfall gauges during 2000-2019. The daily precipitation data were interpolated to the catchments using both the rainfall and weather stations, while other climate variables were interpolated using the weather stations, both via the Inverse Distance Weighting method (The evaluation of this method is available in the supplementary file). The XCO2 products at Level 3 were used to obtain the carbon dioxide in the study region (https://cds.climate. copernicus.eu/, accessed on 8 October 2021). The remotely-sensed based ET and GPP data were used to calibrate and validate the DTVGM model from MOD16A and MOD17A (500 m resolution) at monthly and yearly timestep (https://lpdaac.usgs.gov/products/, accessed on 13 October 2021). Moreover, the annual ET data obtained from Global Land Evaporation Amsterdam Model (GLEAM v3) on a 0.25 • × 0.25 • grid are also used for annual ET evaluation (https://www.gleam.eu/, accessed on 15 November 2021).

The Improved Water-Carbon Coupling DTVGM Ecohydrological Model
The Distributed Time Variant Gain Model (DTVGM) has been developed and successfully applied to many basins in China [34][35][36]. So far, it has been extended to simulating land-surface processes including water, energy and carbon fluxes by coupling with the CASACNP model [37]. The physically based coupled model could capture the mechanism of the interactions between vegetation dynamics and the water cycle quite well and could be used to investigate the responses of ecohydrological processes to changing climate and vegetation, and consequently the attributions of ET and GPP.
In this study, the DTVGM is further improved based on the coupled model by solving the ecohydrological processes iteratively according to the following equations to ensure water balance. At each daily time step, the water balance equation can be expressed as: The above equation could be solved using the Newton iteration method as follows: where P is the daily effective rainfall arriving at the ground surface, AW is the daily mean soil water content, ET is the daily actual evapotranspiration, RS, RI and RG are the daily surface runoff, interflow runoff and groundwater runoff, respectively. The subscripts i and i+1 represent the time steps. The superscript j is the iteration number. The evapotranspiration was calculated with three components including canopy-intercepted water evaporation in the wet leaves, transpiration in the dry leaves and soil evaporation based on the two-source method [37]. For gross primary production, the daily canopy photosynthesis model [38] is selected. The photosynthesis rate are calculated as follows: where A c and A j are the Rubiso-limited and light-limited gross photosynthesis, R d is the daytime leaf dark respiration. Details can be found in [38]. The major calculation procedure to iteratively solve of the ecohydrological processes are shown in Figure 2. Firstly, the initial mean soil water content was set to calculate the soil resistance. The stomatal conductance, photosynthesis and the canopy resistance were calculated, and were affected by soil moisture. Then, the ET with its three components and runoff could be estimated. After that, ET, effective rainfall and runoff were input to Equation (1) to solve the mean soil water content by comparing it with the initial mean soil water content until it converges using the Newton iteration. With these calculation steps, the mean soil moisture was iteratively calculated, and then the ET, GPP and other fluxes were estimated.  Figure 2. The major calculation procedure of the simulation.

Quantification of the Effects on ET and GPP
For a given catchment, the changes of certain ecohydrological indicator Y (ET or GPP) can be affected by meteorological variables (M), elevated atmospheric carbon dioxide (eCO2, C), land cover change (L) and the interactive effects (I) between climate, eCO2 and

Quantification of the Effects on ET and GPP
For a given catchment, the changes of certain ecohydrological indicator Y (ET or GPP) can be affected by meteorological variables (M), elevated atmospheric carbon dioxide (eCO 2 , C), land cover change (L) and the interactive effects (I) between climate, eCO 2 and land cover.
According to [39], the changes in Y (∆Y) at a catchment scale under changing environment can be expressed using Taylor approximation as follows: where ∆M, ∆C, ∆L are the changes of climate variables, eCO 2 and land cover, ∆Y M , ∆Y C , ∆Y L , ∆Y I represent the changes in Y caused by climate change, eCO 2 , land cover changes and the interactive effects between the above three factors.
In order to quantify the contributions of climate change, eCO 2 and land-use and cover change on GPP and ET in the six catchments from 2000 to 2019, ten experimental scenarios were designed, including the baseline (S0), climate variables detrend including atmospheric pressure-detrend (S1), temperature-detrend (S2), humidity-detrend (S3), precipitationdetrend (S4), windspeed-detrend (S5), solar radiation-detrend (S6), carbon dioxide-detrend (S7), land-use and cover changes (S8) and all factor changes scenario (S9). The details of the scenarios are shown in Table 2. The detrend of climate variables is based on the removal of the long-term linear trend during 2000 to 2019 as follows: where X Y i,d and X Y i,d denote the daily original and detrended climate variable at the d day in the year Y i ; X Y i is the average values of the climate variable; a represents the linear trend of the values of the climate variable; Y i refers to the year from 2000 to 2019 and Y 0 is the base year. Table 2. Details of the ten experimental scenarios settings.

S0 Baseline
All the driving factors including historical meteorological variables, CO 2 and land coverage are used to drive the model. S1 Atmospheric pressure-detrend The atmospheric pressure is detrended using historical records. S2 Temperature-detrend The temperature is detrended using historical records. S3 Humidity-detrend The relative humidity is detrended using historical records. S4 Precipitation-detrend The precipitation is detrended using historical records. S5 Windspeed-detrend The windspeed is detrended using historical records. S6 Solar radiation-detrend The sunshine duration is detrended using historical records. S7 Carbon dioxide-detrend The carbon dioxide is detrended using historical records. S8 Land use and cover changes The land use is set at 2000 using historical records.

S9
All-detrend and no land use change All the meteorological variables and CO 2 are detrended and the land use is set at 2000. Each of the ten experimental scenarios were used to simulate the ecohydrological processes during 2000-2019 within the control variate method, i.e., only one influence factor is changed and others are fixed for S1 to S8. Furthermore, S9 denotes all the influence-factor changes. By comparing the simulated ET and GPP of the detrend experiment with the baseline, the contributions of each of the factors on ET and GPP could be quantified in the six catchments of the water-source region of the Xiong'an New Area.

Model Performance
Prior to using the model for quantifying the contributions of different driving factors, the DTVGM model was evaluated in simulating the ET and GPP in the six catchments. The DTVGM simulated ET and GPP were evaluated using the remote sensing product. The model was then calibrated using the data during 2000-2006 and validated for 2007-2010 on a monthly scale. Figure 3 shows the monthly ET by the DTVGM model and the MODIS data from 2000 to 2010 in the six catchments. Generally, the evaluation shows that the seasonal patterns of the simulated ET are consistent with the MODIS data. Comparing to the MODIS data, the ET was found to be underestimated in the winter season. The relative error (RE) is below 7.0%, and the Pearson correlation coefficient (r) is above 0.73 except for the THC, and the root mean squares error (RMSE) is below 0.67 mm in the six catchments as shown in the statistical results in Table 3. Figure 4 displays the estimated monthly GPP by DTVGM and MODIS in the catchments. The results show that the simulated GPP captures the MODIS GPP patterns and magnitudes reasonably well. The RE values are almost all below 11.0% except for the BSC, the r values even reach up to 0.90 except for the BGC about 0.86, and the RMSE is below 1.24 g C/m 2 (Table 3). Figure 5 shows the comparison of the annual simulated and remote-sensing based ET and GPP in the region. The least-squares linear fits to the simulated ET and GPP is near to the 1:1 line, while slightly underestimated, and the R 2 value of linear goodness of fit test is about 0.51 and 0.71 for ET and GPP, respectively. Furthermore, the simulated annual ET was also compared with the GLEAM annual ET data. The results show that the RE is −4.3%, 4.1% −5.7%, −17.2%, −12.7% and −20.9% in the GBC, JMC, BSC, THC, DSC and CHC, respectively and the average RE is −8.7%. Figure 5c shows that the simulated annual ET by DTVGM is mostly between the values of MODIS and GLEAM. Overall, the results are acceptable for ET and GPP simulations in the region. compared with the GLEAM annual ET data. The results show that the RE is −4.3 −5.7%, −17.2%, −12.7% and −20.9% in the GBC, JMC, BSC, THC, DSC and CHC, tively and the average RE is −8.7%. Figure 5c shows that the simulated annua DTVGM is mostly between the values of MODIS and GLEAM. Overall, the res acceptable for ET and GPP simulations in the region.

Changes in Climate Variables and Land Coverage
(1) Changes in the climate variables Climate variables are supposed to be the most important driving factors that affect hydrological processes, the vegetation dynamics and their interactions. Therefore, the changing trends of the climate variables, including atmospheric pressure, temperature, humidity, precipitation, windspeed, sunshine duration and carbon dioxide (CO 2 ) were detected using the linear regression method and Mann-Kendall [40,41] test, which are both commonly used in the hydrometeorological trend detection. The monthly results in Table 4 show that the atmospheric pressure decreased slightly in the BGC and JMC but did not exceed the significance level, while it shows no obvious changing trends in the BSC, THC and DSC. The temperature increased significantly in the THC and DSC, while it decreased significantly in the BSC. The humidity decreased in five catchments except for the BSC and the decreasing trends are significant in the THC and DSC. The precipitation increased significantly in the five catchments except for the THC. The windspeed shows significant reduction in the BGC, DSC and CHC. The sunshine duration also displays a decreased trend, particularly in the BSC and CHC. The CO 2 presents a significant increasing trend in the region within the slope of 2.24 ppmv/a. The monthly Mann-Kendall changing trends of the climate variables in the six catchments are shown in Figure 6. The color bar indicates the changing trends in the variables. The atmospheric pressure displayed a significant decrease in summer and autumn, particularly in August and September. The humidity shows a dropping trend, particularly in January. The monthly precipitation increased mostly in February, April, July, August and October, particularly in the wet season, while it decreased in the dry season, especially in January, November and December. For the sunshine duration, it decreased almost throughout the year, especially for July. The windspeed decreased in most of the months, particularly during March, April, October and November. Unlike the changes of windspeed, the monthly temperature increased gradually in most of the months, and fell mostly in February and July.  (2) Land-use and cover change Land-use and cover change (LUCC) is recognized as important factors that influe (2) Land-use and cover change Land-use and cover change (LUCC) is recognized as important factors that influence evapotranspiration and vegetation carbon assimilation. Therefore, the land use changes were also investigated in the six catchments as shown in Figure 7. In the six catchments, the dominant land coverage is forest, grassland and cropland, covering almost nearly 95% of the whole basin. In the BGC, the major land cover is forest (83%) and grassland (9%). The forest decreases to 64% in the JMC and the rest are mainly grassland (22%), cropland (10%). While in the BSC, the main land cover is grassland (49%), the forest decreases to 36%, cropland is about 10%, respectively. In the THC, the grassland increases to 51% and the cropland increases to 21%, while the forest decreases to 22%. In the DSC, the chief land cover is still grassland (57%), while the cropland decreases to 7% and the forest increase to 32%. In the CHC, the main land covers are forest (48%), grassland (38%) and the cropland constitutes about 8%. Figure 7 displays the land-coverage changes from 2000 to 2020 in the six catchments. Generally, the study area is a mountainous region with relatively few human activities, and the urban area increased by only 0.8% due to urbanization. A growth of 0.2% in the forest was found in the BGC. While in the JMC, the cropland and forest both decreased by 0.2% and the urban area increased by 0.4%. In the THC, the forest increased by 0.2%, the cropland and grassland area decreased by 0.5% and 0.7%, respectively. In the DSC, only the grassland was reduced by 2.3% and the cropland and forest both increased by 0.2%. In the CHC, the forest and grassland decreased by 1.5% and 0.6% as the result of urban expansion. Table 5 shows the changing trends of land coverage using the linear regression method and Mann-Kendall test. The result shows that most of the land-use and coverage changed slightly, except for urban land, during the past 20 years. However, it should be mentioned that the urban land covers small areas.   Figure 7. Land use and cover change in the six catchments(a-f).

Changes in ET and GPP
The changes in ET and GPP were examined from 2000 to 2019 based on the DTVGM simulation using historical observations in the six catchments ( Table 6). The ET and GPP in the region show a significant increase (α = 0.01) of 2.40 mm/a and 6.00 gC/m 2 /a, respectively. Figure 8 shows simulated annual ET in the six catchments during the past years. The ET increased at the rate of 6.92 mm/a, 3.49 mm/a, 0.44 mm/a, 1.59 mm/a, 1.54 mm/a and 1.20 mm/a for GBC, JMC, BSC, THC, DSC and CHC, respectively. The Mann-Kendall test reveals that the ET increased significantly in the BGC, JMC, DSC (α = 0.01) and THC (α = 0.05), while it is not obvious for the BSC and CHC. Figure 9 displays the modeled annual GPP in the six catchments. The GPP shows significant increasing trends in all of the six catchments at a significance level of 0.01. Additionally, the changing trends are 14.30 gC/m 2 /a, 7.73 gC/m 2 /a, 6.41 gC/m 2 /a, 4.25 gC/m 2 /a, 4.38 gC/m 2 /a, 5.28 gC/m 2 /a for GBC, JMC, BSC, THC, DSC and CHC, respectively.
Furthermore, the variations in ET within the year were investigated as shown in Figure 10a. In the BGC, the ET increases in all the months, especially in the spring and summer season from April to September. However, the ET shows decreasing trends in some months in the autumn and winter seasons and increased mostly in spring and summer seasons in the other five catchments. Particularly, the ET decreased from autumn to early spring in the following year in the BSC. Figure 10b displays the changing pattern in GPP within the year. The GPP increased in nearly all the months except for February. The major increase in monthly GPP mainly occurred in the spring and summer seasons, from April to September. Figure 10c summarizes the mean changes in ET and GPP in the water source region, and both ET and GPP increased mainly in spring and summer, while the uptrends were relatively smaller in autumn and winter. spring in the following year in the BSC. Figure 10b displays the changing pattern in G within the year. The GPP increased in nearly all the months except for February. The ma increase in monthly GPP mainly occurred in the spring and summer seasons, from Ap to September. Figure 10c summarizes the mean changes in ET and GPP in the water sou region, and both ET and GPP increased mainly in spring and summer, while the uptren were relatively smaller in autumn and winter.

(1) Attributing changes in ET
To identify the dominant factors that affect the changes in ET and GPP, ten experimental scenarios including the baseline and the detrended scenarios were conducted. The contributions of the climate variables, land-use and cover change (LUCC), elevated atmospheric carbon dioxide (eCO 2 ) and their interactive effects on ET are shown in Table 7. The results show that the contribution rate of precipitation to the growth ET in all the six catchments is the largest, with an average contribution of 104.8%. The contributions of the remaining influence factors are temperature (12.8%), solar dimming (−12.4%) and carbon dioxide (−7.2%), and then the interactive effects (2.3%), humidity (2.0%), atmospheric pressure (−1.3%), wind speed (−0.8%), and LUCC (−0.2%). In the BGC, the results show that precipitation, temperature, windspeed, atmospheric pressure, humidity and LUCC contributed positive effects, and the effects of sunshine duration, eCO 2 and interactions were negative. However, in the JMC, only precipitation, temperature and interactive effects show positive effects and the other climate variables all contributed negative effects. It should be noticed that the interactive effects could reach up to 21.2% in the JMC, which contributes more than temperature and the other factors. In the BSC, the precipitation could affect the ET increase by up to 164.1%, and sunshine duration contributed to a considerable ET decrease, by 42.5%. In the THC and DSC, all the influence factors have similar contributions. The precipitation, temperature, humidity and windspeed had positive effects on the ET increase, while the sunshine duration, eCO 2 , LUCC, atmospheric pressure and interactions led to ET reductions. In the CHC, the precipitation, temperature and windspeed continue to contribute to the ET increase, while the interactions and atmospheric pressure also increase the ET, whereas humidity, sunshine duration, eCO 2 and LUCC contribute negative effects. Figures 11 and 12a display the contributions of the factors on monthly and seasonal ET changes. Generally, precipitation is the dominant factor for ET changes. The positive effects of precipitation occur throughout the year except for the negative effects that occur in January and December in winter. The temperature has a negative effect on ET changes in summer, while it leads to an ET increase in other three seasons, particularly from November to January in winter. The effect of sunshine duration on ET changes is negative in summer and autumn, particularly from September to November, which is the opposite in winter and spring. The eCO 2 contributes to ET increase in autumn and ET decrease in other seasons, especially in winter. The humidity leads to an ET decrease in autumn and ET increase in other seasons, particularly in January and December in winter and it is the dominant factor contributing to ET change in winter. The windspeed has a positive effect on ET increase in most of the months except for June to August in summer. The effects of LUCC and atmospheric pressure on ET are almost rare throughout the year. Additionally, the interactive effects lead to ET reduction, mostly in spring, and its increase in autumn and winter, particularly in September. Overall, precipitation is the dominant factor that increases ET in spring, summer and autumn and decrease ET in winter. In spring, the interactive effects contribute to most of the ET reduction. In autumn, the humidity and sunshine duration are the main factors that reduce ET. However, humidity and temperature are the major reasons for ET increase in winter.  Figure 11. Attributions of monthly ET changes in the six catchments.
(2) Attributing changes in GPP Based on the ten experimental scenarios, the effects of the climate variables, LUCC, eCO2 and their interactions on GPP are also quantified in Table 8. The results show that eCO2 contribute the largest GPP increases in most of region within average contribution about 56.7%. The precipitation and temperature also have significant influences on the increase in GPP with the effects of 34.8% and 27.8%, while the sunlight duration, humidity, wind speed, atmospheric pressure, LUCC and the interactive effects have negative effects on GPP growth with the contributions of −6.8%, −4.8%, −3.5%, −2.4%, −1.0% and −0.8%, respectively. We could also find that the effects of the influence factors on GPP changes are various from each catchment in the region. In the BGC and JMC, the precipitation contributes the largest GPP increase, and eCO2 and temperature also contribute significant increase in GPP. Meanwhile, the interactive effects (39.5%) should also be focused on in the JMC. The effects of atmospheric pressure, windspeed and LUCC are different in the two catchments. In the BSC, the eCO2 have a large effect increasing the GPP by 108.3%, the humidity and precipitation also have positive effects, while temperature, sunshine duration and the interactions have negative effects. In the THC and DSC, the dominant factor is eCO2 with the contributions of 66.2% and 73.8%, and the followed factors are temperature and precipitation. The interactive effects could reduce the GPP by 39.7% in the DSC. The effects of the other impact factors are generally small. In the CHC, the most important factors are eCO2 and temperature increasing the GPP by 66.8% and 45.7%, while precipitation has a negative effect of −6.6% which is different from the results in other catchments. Meanwhile, the sunshine duration make a large reduction on GPP by 21.3%. The LUCC have a larger negative effect about −4.9% in the CHC comparing to other catchments. Figures 12b and 13 presents the contributions of the factors on monthly and seasonal GPP changes. The eCO2 contribute a large GPP increase throughout the year except for small decrease in February. Additionally, the eCO2 is the dominant factor contributing to GPP increase in summer and autumn, specifically from June to October. Furthermore, Figure 11. Attributions of monthly ET changes in the six catchments.
(2) Attributing changes in GPP Based on the ten experimental scenarios, the effects of the climate variables, LUCC, eCO 2 and their interactions on GPP are also quantified in Table 8. The results show that eCO 2 contribute the largest GPP increases in most of region within average contribution about 56.7%. The precipitation and temperature also have significant influences on the increase in GPP with the effects of 34.8% and 27.8%, while the sunlight duration, humidity, wind speed, atmospheric pressure, LUCC and the interactive effects have negative effects on GPP growth with the contributions of −6.8%, −4.8%, −3.5%, −2.4%, −1.0% and −0.8%, respectively. We could also find that the effects of the influence factors on GPP changes are various from each catchment in the region. In the BGC and JMC, the precipitation contributes the largest GPP increase, and eCO 2 and temperature also contribute significant increase in GPP. Meanwhile, the interactive effects (39.5%) should also be focused on in the JMC. The effects of atmospheric pressure, windspeed and LUCC are different in the two catchments. In the BSC, the eCO 2 have a large effect increasing the GPP by 108.3%, the humidity and precipitation also have positive effects, while temperature, sunshine duration and the interactions have negative effects. In the THC and DSC, the dominant factor is eCO 2 with the contributions of 66.2% and 73.8%, and the followed factors are temperature and precipitation. The interactive effects could reduce the GPP by 39.7% in the DSC. The effects of the other impact factors are generally small. In the CHC, the most important factors are eCO 2 and temperature increasing the GPP by 66.8% and 45.7%, while precipitation has a negative effect of −6.6% which is different from the results in other catchments. Meanwhile, the sunshine duration make a large reduction on GPP by 21.3%. The LUCC have a larger negative effect about −4.9% in the CHC comparing to other catchments. Figures 12b and 13 presents the contributions of the factors on monthly and seasonal GPP changes. The eCO 2 contribute a large GPP increase throughout the year except for small decrease in February. Additionally, the eCO 2 is the dominant factor contributing to GPP increase in summer and autumn, specifically from June to October. Furthermore, precipitation has positive effect on GPP change in all months, and it is the dominant factor leading to an increased GPP in spring. The temperature increases the GPP during a year, particularly from January to March, while it has a slight negative effect in summer. The temperature is the most important factor that contributes to the GPP increase in winter and it also has a great contribution in spring. The windspeed increased the GPP in spring and autumn and leads to a slight decrease in the GPP in summer and winter. The sunshine duration has a negative effect in summer and positive effects in the other seasons. Generally, the effects of atmospheric pressure, humidity and LUCC are relatively small compared to the above influence factors. The effects of the interactions of the factors on GPP should be of concern, as they can lead to a considerable decrease in GPP in spring, from March to May, and an increase in winter, in January and February. Overall, in spring, the precipitation is the dominant factor that increases GPP, while the temperature and eCO 2 also make great contributions, while the interactions have considerable effects on GPP reduction. The eCO 2 has a dominant effect on GPP increases in summer and autumn, followed by precipitation. Conversely, in winter, the temperature is the major factor that increases GPP, and the interactions of the factors should also be considered.

Discussion
In this study, the improved water-carbon coupled DTVGM ecohydrologica was used to simulate the ET and GPP in the six catchments of the water-source re the Xiong'an New Area. Overall, the model was able to capture the ET and GPP va quite well in the simulation periods, with relatively error of mostly below 10.0% high correlation coefficient of up to around 0.9. Compared to the MODIS da DTVGM simulation results also show some underestimations in the dry season while GPP is underestimated in the wet season and overestimated in the dry which is relevant to the vegetation, runoff and soil parameterization. The simulate also consistent with the GLEAM ET Product [42] with a relative error of about which further proved the reliability of the DTVGM results. Some researches show large seasonal amplitudes MODIS GPP data is overestimated in the wet season a values are underestimated. The reason may be due to the artificially high values of FPAR and the MODIS light-use efficiency approach, which does not account for th ration of canopy photosynthesis under clear sky conditions [43][44][45]. Based on the D simulation, ET decreases from northeast to southwest and ranges from arou mm/year to 500 mm/year, with an average value of 446.67 mm/year. GPP rang around 600 gCm 2 /year to 800 gCm 2 /year from northeast to southwest, with a mea of 702.57 gCm 2 /year. These results also agreed with previous studies that included and GPP simulations in this region [28,46]. Overall, we suppose that the DTVGM

Discussion
In this study, the improved water-carbon coupled DTVGM ecohydrological model was used to simulate the ET and GPP in the six catchments of the water-source region of the Xiong'an New Area. Overall, the model was able to capture the ET and GPP variations quite well in the simulation periods, with relatively error of mostly below 10.0% and a high correlation coefficient of up to around 0.9. Compared to the MODIS data, the DTVGM simulation results also show some underestimations in the dry season for ET, while GPP is underestimated in the wet season and overestimated in the dry season, which is relevant to the vegetation, runoff and soil parameterization. The simulated ET is also consistent with the GLEAM ET Product [42] with a relative error of about −8.7%, which further proved the reliability of the DTVGM results. Some researches show that the large seasonal amplitudes MODIS GPP data is overestimated in the wet season and low values are underestimated. The reason may be due to the artificially high values of MODIS FPAR and the MODIS light-use efficiency approach, which does not account for the saturation of canopy photosynthesis under clear sky conditions [43][44][45]. Based on the DTVGM simulation, ET decreases from northeast to southwest and ranges from around 400 mm/year to 500 mm/year, with an average value of 446.67 mm/year. GPP ranges from around 600 gCm 2 /year to 800 gCm 2 /year from northeast to southwest, with a mean value of 702.57 gCm 2 /year. These results also agreed with previous studies that included ET [18] and GPP simulations in this region [28,46]. Overall, we suppose that the DTVGM ecohydrological model could provide reasonable simulations for ET and GPP in the experimental scenario analysis to identify the dominant impact factors.
The attributions of ET and GPP changes were evaluated based on ten experimental scenarios. At the annual scale, ET increased by about 2.4 mm/a in the region. The increasing precipitation dominated the ET uptrend absolutely, probably on account of the water limit in the semi-arid region [47]. Additionally, similar results were also found in other semi-arid regions [48]. Meanwhile, rising air temperature and decreasing humidity also contribute to the ET uptrend by increasing the vapor-pressure deficit and then accelerating ET [49]. The elevated atmospheric CO 2 concentration (eCO 2 ) and solar dimming are the main causes for the decrease in ET, which is similar to the findings in [19]. It has been widely recognized that eCO 2 tends to induce the stomatal closure, resulting in a decrease in transpiration [50]. Additionally, solar dimming reduces the solar radiation on the land surface and limits the available energy allocation to the latent heat. The decreasing windspeed also contribute to a small decrease in ET. The GPP increases to about 6.0 gC/m 2 /a in the region. The eCO 2 , increasing precipitation and rising temperature are responsible for the GPP increases, while solar dimming, decreasing humidity and windspeed contribute to decreases in GPP. These results are also supported in previous studies [21,51]. The eCO 2 and rising temperature may directly lead to increasing carboxylation rates and GPP. Increasing precipitation promotes the available water resources for plant growth and also contributes to a GPP increase with a sigmoid relationship [52]. Meanwhile, decreasing humidity leads to vapor-pressure deficit increase, which would result in stomatal conductance declines and a slowing of photosynthesis [53]. Besides the negative effects of solar dimming on ET, it would also reduce the photosynthetically active radiation for carbon fix [54]. The results also show that the effects of LUCC on both ET and GPP are slightly negative, at about −0.2% and −1.0%, respectively. This may be related to the slight land-use change, less urbanization and human activities in the mountainous regions.
On the seasonal scale, the increases in ET and GPP mainly occur in spring and summer and the increases are relatively small particularly in winter. Similar results are also found in the northeastern part of China [55]. The major reason for the increase in ET and GPP in spring and summer is obvious rising temperature and increasing precipitation in the two seasons, which may be due to the availability of more accumulated temperature and water during the growing season. In autumn, significant solar dimming should lead to lower ET and GPP uptrends, while the humidity decreases significantly and lower increases in precipitation and temperature cause little ET and GPP increase in winter. Besides, the interactive effects of the climate variables, eCO 2 and LUCC are significant in spring and winter, and have positive effects on the ET (16.2%) and GPP (27.8%) uptrend in spring and negative effects on the ET (−24.4%) and GPP (−38.0%) uptrend in winter. Meanwhile, the interactive effects are also significant in autumn in ET increase (42.9%). This finding is consistent with the conclusion that interactive effects are important in water-limited catchments, whereby the ecosystem structure and function are mainly constrained by water [39,56]. Additionally, this may be related to the strategy mechanism of the ecosystem in adapting to changing circumstances.
The ET and GPP increase in the water-source region of the Xiong'an New Area, and increase more during the growing season from April to September, which can cause adverse impacts more easily on agriculture and the ecosystem [57,58]. As we know, water resources were relatively scarce in semi-arid area, the increased GPP was helpful for regional-vegetation recovery and more water consumption may have led to runoff reduction under global climate change in this region, subsequently resulting in a water shortage for social-economic usage and river environmental requirements. Besides, the increased ET mainly caused by increased transpiration will lead to more soil water-consumption of the soil layer and intensify the drought, resulting in high ecological sensitivity and serious vulnerability in water-limited zones [6]. The attributions of ET and GPP changes reveal that the multitude of factors dominating water and carbon fluxes are different in the region over different time scales. Precipitation is responsible for both ET and GPP changes, which proved the determination of the vegetation by water in the semi-arid region. Rising temperature and eCO 2 also contribute greatly to both ET and GPP changes.
Increasing precipitation, rising temperature and eCO 2 lead to greater GPP and ET for vegetation. Meanwhile, more precipitation may lead to greater water availability and fewer droughts, but a warmer temperature could enhance evaporation from the soil, producing more droughts [59]. The complex mechanisms of ET and GPP changes may bring more uncertainties of the drought in future and pose challenges to drought management.
The attributions of ET and GPP changes are quantified in this study by considering climate variables, eCO 2 and LUCC, at both annual and seasonal timescales using the water-carbon coupling DTVGM ecohydrological model. We consider the water-carbon interactions by better capturing the mechanisms of ET and GPP changes, and it could especially reflect the interactive effects of hydrological processes on GPP by the climate forces and the vegetation dynamics on ET by the elevated atmospheric carbon dioxide, which cannot be well understood by the ecological model or rainfall-runoff relationship based on hydrological model. Besides, as global climate change may increase drought frequency and intensity in the future [60], more effective models to evaluate ET and GPP are needed considering anthropogenic influences such as irrigation [61]. The improved DTVGM ecohydrological model in this study could characterize the ET and GPP variations reasonably, providing a useful tool to study drought prediction and also the drought effects on vegetation carbon-uptake. However, it should also be further improved by considering human activities and more studies should be focused on the model parameterization and deeper insights into the spatiotemporal patterns. Furthermore, future work should also continue to assess the driving factors.

Conclusions
In this study, the attributions of the evapotranspiration (ET) and gross primary production (GPP) changes including climate variables, elevated atmospheric carbon dioxide (eCO 2 ), land-use and coverage change (LUCC) were quantified in the six catchments in the water source region of the Xiong'an New Area using the improved water-carbon coupled DTVGM ecohydrological model. The ecohydrological model could capture the ET and GPP variations reasonably in the study area. The modelling experimental scenarios show that annual ET and GPP increased significantly, particularly in the spring and summer seasons. At the annual scale, increasing precipitation dominated the ET uptrend absolutely and air temperature also contributed to the ET uptrend. The eCO 2 concentration and solar dimming were found to be the main causes for the ET decrease. For the GPP changes, eCO 2 , increasing precipitation and rising temperature dominant the GPP increases. Compared to climate change, the LUCC has small effects on both ET and GPP in terms of the few changes in land coverage. At the seasonal scale, the modelling results show different effects in the four seasons. Overall, precipitation, temperature and eCO 2 are the main causes for ET and GPP changes. However, the decreasing humidity was found to be the dominant factor for ET increase in winter, while solar dimming imparts strong effects on ET decreases in autumn. It should also be noted that the contribution of the interactive effects cannot be neglected at the seasonal scale, particularly in terms of the considerable decreases in ET and GPP in spring, increases in ET in autumn and winter and increase in GPP in winter. This study highlights the importance of considering water-carbon coupling on the attribution of ET and GPP and the differentiation of the effects by the abovementioned influence factors at the annual and seasonal scales.