The aim of this work is to investigate the soil water budget across China by means of hydrological modeling under current and future climate conditions and to evaluate the sensitivity to soil parameters. For this purpose, observed precipitation and temperature data (1981–2010) and climate simulations (2021–2050, 2071–2100) at high resolution (about 14 km) on a large part of China are used as weather forcing. The simulated weather forcing has been bias corrected by means of the distribution derived quantile mapping method to eliminate the effects of systematic biases in current climate modeling on water cycle components. As hydrological models, two 1D models are tested: TERRA-ML and HELP. Concerning soil properties, two datasets, provided respectively by Food and Agriculture Organization and U.S. Department of Agriculture, are separately tested. The combination of two hydrological models, two soil parameter datasets and three weather forcing inputs (observations, raw and bias corrected climate simulations) results in five different simulation chains.
The study highlights how the choice of some approaches or soil parameterizations can affect the results both in absolute and in relative terms and how these differences could be highly related to weather forcing in inputs or investigated soil. The analyses point out a decrease in average water content in the shallower part of the soil with different extents according to climate zone, concentration scenario and soil/cover features. Moreover, the projected increase in temperature and then in evapotranspirative demand do not ever result in higher actual evapotranspiration values, due to the concurrent variations in precipitation patterns.
Water balance models ; Model comparison ; Bias correction ; Hydrological impacts ; Soil parameters ; China
In recent years, some relevant studies have shown that the global warming could affect the precipitation patterns and the water cycle at the global and regional scale and at different time scales (Allen and Ingram, 2002 ; Held and Soden, 2006 ; Lenderink and van Meijgaard, 2010 ). However, under the current observational and modeling uncertainties, the Intergovernmental Panel on Climate Change (IPCC) assesses only with “medium confidence” the anthropogenic contribution to changes in precipitation patterns, although it recognizes as “virtual certain” the ongoing processes of global warming and the emerging of new evidences about the variations in global precipitation (IPCC, 2014 ).
This topic and its huge implications for the development of effective adaptation strategies in crucial sectors (agriculture, natural disaster mitigation, water resources, water supply) have meant numerous researches aimed to investigate how climate changes may affect the entire water cycle (Fischer et al., 2007 ). The main studies focus on:
Of course, since the significance, extension and extreme geomorphological and climate heterogeneity, numerous studies have investigated such issues for China. China played indeed a key political and economic role in the mitigation of climate change and it is considered as one of the most vulnerable societies in terms of water resources and agriculture sustainability (IPCC, 2014 ; SC (State Council), 2011 ). Moreover, significant resources have been invested in climate change adaptation strategies and actions (SC, 2014 ). Concerning the topics already introduced, Wu et al. (2015) adopted the Soil and Water Assessment Tool (SWAT) model forced by available historical weather observation to assess the ongoing variations in the water cycle potentially induced by climate change for the Xiangjiang River Basin. On the same river basin, Luo et al. (2013) adopted the same tools for investigating the effect of current and planned scenarios of land use. Zhang et al. (2009) simulated evapotranspiration and soil water balance for different tillage systems on the 2010–2040 time span using as driver the monthly precipitation data provided by statistical downscaling of the Global Climate Model (GCM) HadCM3 projections. Wang et al. (2013) estimated the impacts of climate change on water demand for agriculture in ten basins covering the entire Chinese territory. They coupled the climate inputs provided by the Regional Climate Model (RCM) PRECIS (at 0.44° horizontal resolution; Yang et al., 2010 ) under A1B, A2 and B2 emission scenarios to the China Water Simulation Model (CWSM), assessing water demand and supply and allocating water among its users within the river basins. They showed how the effects could vary in function of severity of emission scenarios but, at the same time, proper adaptation strategies may satisfactorily cope with these issues.
However, for China, relatively much fewer studies have had as specific goal the investigation on potential changes induced by climate change in the soil moisture at local or regional scale. At local scale, e.g. Zhou et al. (2011) (again through the SWAT approach) assessed the effect of climate change on soil moisture dynamics and hydrological variables in an intact-forested watershed in southern China. Through the analysis of observed data, they show how current soil drying processes could be associated to relevant changes in weather forcing. Similarly, Jiang et al. (2007) perform an exhaustive comparison work analyzing the response of six different hydrological models in assessment of variations of soil hydrological balance induced by climate change (using hypothetical future scenarios) in the Dongjiang Basin (southern China).
On the other side, at country level (not only for China), this topic has been hardly investigated until now (Wu et al., 2015 ) also probably due to the limited availability of reliable and continuous soil moisture datasets around the world (Keshta et al., 2012 ). Among available studies, Komuscu et al. (1998) investigated the variations of soil moisture availability in the Southeast Anatolia (Turkey) both under the prescribed variations in weather forcing, and under those provided by the GCM, using a modified version of Thornthwaite water balance model. Naden and Watts (2001) evaluated the variations of monthly soil moisture in five typical areas of ecological interest in the United Kingdom taking into account the actual land covers and soil types. Bormann (2009) performed a sensitivity study for five different regions in Germany, aimed to understand the joined role of soil type and lower boundary conditions in response to variation of weather forcing induced by climate change. Keshta et al. (2012) adopted the Generic System Dynamics Watershed (GSDW) model (Keshta et al., 2009 ) to simulate the potential future soil moisture and evapotranspiration dynamics; the weather input forcing are obtained through the Statistical DownScaling Model (SDSM) proposed by Wilby et al. (2002) of Coupled General Circulation Model (CGCM3) under A2 and B1 emission scenarios. Finally, Bormann (2012) utilized the SIMULAT model (Bormann, 2008 ), a Soil Vegetation Atmosphere Transfer scheme (SVAT) for analysis of virtual soil columns, to assess the soil texture-specific sensitivity of simulated soil moisture to projected climate change. To this aim, he also took into account, in a simplified way, the role played by upper (land cover) and lower (shallow groundwater/free drainage) hydraulic boundary conditions of soil specimens.
Given the state of the art, the availability of climate projections at a very high resolution (0.125°, about 14 km) up to 2100 from RCM COSMO_CLM (Bucchignani et al., 2014 ) under RCP4.5 and RCP8.5 scenarios (Meinshausen et al., 2011 ) for a large part of the Chinese territory, grants an opportunity to attempt partly bridging such gap estimating the main variations assessed in soil water budget for the principal climate zones retrievable in the investigated domain.
The analysis are carried out adopting two hydrological models, two datasets for retrieving hydraulic and thermal soil properties, three weather forcing inputs. In this way, the effects of parameterizations and parameters on the assessment of seasonal and yearly soil water budget considering also the potential effect of soil type and land-use can be evaluated. Moreover, a first estimation about how current biases in weather forcing returned by climate simulations could affect the assessment of the soil water budget and its variation under climate change is also provided. The investigation is performed on five reference “boxes” located in four main climate areas; every box is selected as representative of the most widespread land cover and soil type in that area. Essentially, two time horizons are investigated: near time horizon 2021–2050 and long time horizon 2071–2100 under RCP4.5 and RCP8.5 scenarios. Then in the rest of this paper, firstly material and methods are displayed. Afterwards, attempting to estimate the main variations assessed in the soil water budget for the main climate zones identified in China, under the effect of climate change, and to quantify the reliability of the results under climate and hydrological modeling uncertainties, results are reported and commented in Section 3 . Finally, in conclusions, the main findings are reviewed and the possible future steps are identified.
Given the great extension of China, climate features vary considerably among different areas. According to the Köppen–Geiger classification (Kottek et al., 2006 ), the Chinese territory is characterized by six climates (Fig. 1 a): cold desert climate (BWk) in the Northwest (NW), tundra climate (ET) in the Southwest (SW), humid subtropical climate (Cfa and Cwa) in the Southeast (SE), humid continental climate (Dwb) and subarctic climate (Dwc) in the Northeast (NE). In this study, we refer to the climate areas BWk, ET, Cfa and Cwa. Within each of these areas, a box, characterized by the dominant soil type (Fig. 1 b) and land cover (Fig. 1 c), is selected as representative of the area with the exception of the Cfa climate zone where the almost equal presence of two main different soil types (clay and loamy clay, see following sections for details) suggests analyzing two different boxes. As results, the effects of climate change on the water cycle are investigated on five boxes belonging to four different climate zones.
Koppen–Geiger classification (a), soil type (b) and land cover (c) for the continental part of China.
Concerning climate projections, in general, climate simulations over China return further growing temperatures (3–5 °C) while the average precipitation patterns are expected to vary in different ways. In this regard, further valuable details can be retrieved by the analysis performed by Bucchignani et al. (2014) . The adopted simulation chain in Bucchignani et al. (2014) is composed by the GCM CMCC_CM (Scoccimarro et al., 2011 ) with an horizontal resolution of about 0.75° and the RCM COSMO_CLM (Doms and Schättler, 2002 ; Rockel et al., 2008 ) at a very high resolution (0.125°, about 14 km); the former driving the latter as lateral boundary conditions. Both climate models are forced, as defined in IPCC-AR5, by the Coupled Model Intercomparison Project Phase 5 (CMIP5, Taylor et al., 2012 ) historical experiment up to 2006, and by the RCP4.5 and RCP8.5 emission scenarios from 2006 to 2100 (Meinshausen et al., 2011 ). Improvements associated to their optimized configuration and resolution increase are highlighted in Bucchignani et al. (2014) , where it is also pointed out how the lack of high resolution observational datasets may affect a proper validation, especially for high resolution simulations. According to the Bucchignani et al. (2014) projections, a general warming is expected over the whole area with values consistent with previous studies while an increase in precipitations could occur in the northern part of the domain (except in summer) with the exception of SE area where a precipitation reduction is expected, mainly in winter.
Fig. 2 shows the projected variations, according to Bucchignani et al. (2014) simulations, in monthly cumulative precipitations (as ratio) and temperatures (as difference) for the five boxes on the near (2021–2050) and long (2071–2100) time horizons under the RCP4.5 and RCP8.5 scenarios. In particular, concerning temperatures, a generalized increase is detected: on near term, slight variations are recognized between RCPs with growth for all locations of about 1 °C; on long term, more pronounced differences are returned: for the RCP4.5 scenario, the increase stands between 2 °C and 3 °C while it exceeds 3 °C in RCP8.5 (reaching 6 °C for the SE). Variations in precipitation patterns display more complex trends depending on seasonality and location. Box 1 (in the NW area) is characterized by scarce yearly precipitation values and slight future estimated variations. For box 2 (in the SW area), decreases are detected under RCP8.5 mainly in November and December (corresponding to the dry season in the area) while increases are estimated in the wet season (May–August mainly under RCP4.5). For boxes 3 and 4 (in the SE area), again, large decreases occur in the dry season (November–February) especially under RCP8.5 (up to 60%) while increases are obtained in the wetter time span. Finally, for box 5 (in the NE area), high variations (increases up to 100%) are assessed in the drier season (October–March) while they appear to be limited and fluctuating in the remaining months.
Projected variations in monthly cumulative precipitation and temperature for the five boxes on near (2021–2050) and long (2071–2100) time horizon (relative to 1981–2010).
Beyond the relevance of such data, the “strong non linearity” between climate forcing and soil water balance components and the presence of systematic biases in climate simulations altering quantitatively the expected weather forcing, entails that the climate projections cannot directly provide insights on future variations at the local scale of the soil water balance.
Although the improvements in the evaluation of weather forcing on the area achieved by Bucchignani et al. (2014) are mainly associated to higher horizontal resolutions and better representation of physical processes, biases remain as displayed in Fig. 3 , showing, for the five boxes and for the 1981–2010 control period, the ratio (difference) in monthly cumulative precipitations (temperatures) between observed (Wu and Gao, 2013 ) and simulated data (Bucchignani et al., 2014 ). Results show that precipitation biases are limited in the SE area and relevant (essentially in terms of overestimation in boxes 2 and 5) in the other cases. For temperatures in areas characterized by an overestimation of precipitations, an underestimation (1–2 °C) is detected. For other reference areas, biases stood on about +2 °C.
Ratio in monthly cumulative precipitation, and difference in monthly average temperature for simulated against observed (grey line) and bias corrected simulated against observed (black curve) in 1981–2010.
In this perspective, the adoption of a bias correction approach (Ehret et al., 2012 ; Piani et al., 2010 ; Teutschbein and Seibert, 2012 ; Teutschbein et al., 2011 ; Villani et al., 2015 ) is needed to investigate the soil water balance. To this aim, the non-parametric quantile mapping approach (Gudmundsson et al., 2012 ; Teutschbein and Seibert, 2012 ) is used to correct the precipitation and temperature values provided by numerical climate models. Bias correction approaches are based on a transfer function that imposes the equality in probability between observed and simulated data over a control period:
where FOBS and FCTRL are the cumulative distribution functions of daily observed and simulated variable X . Assuming that the transformation defined by Eq. (1) is valid over the whole simulation time period, the bias corrected value X∗ of the simulated variable X is obtained as
The main differences between the distribution mapping methods are due to h transformation. To this aim, theoretical distribution, parametric relationships quantile–quantile or directly the empirical quantiles are generally adopted. In this work, the last mentioned approach has been preferred ( Gudmundsson et al., 2012 ; Villani et al., 2015 ); it allows greater versatility and requires less a-priori assumptions about statistical distributions reproducing observed and simulated daily weather data. On the other side, higher caution should be paid in their adoption when used beyond the calibration intervals ( Gudmundsson et al., 2012 ; Villani et al., 2015 ).
As it is evident in Fig. 3 , showing a comparison between bias corrected values and the observed data, the non-parametric quantile mapping approach succeeds in almost all cases in nullifying the biases of initial climate simulations. However, as it represents a further step in the simulation chain, additional assumptions should be taken into account (Ehret et al., 2012 ; Maraun, 2013 ). For example, its adoption could partly modify the climate signal returned by RCM simulations, as reported in Fig. 4 where the data for temperature and precipitation (again respectively in terms of difference and ratio) are shown for the two time horizons under RCP4.5 scenario.
Climate signal returned by RCM simulations: ratio in monthly cumulative precipitation and difference in monthly average temperature for simulated and bias corrected simulated for the 2021–2050 and 2071–2100 time horizons under RCP4.5 scenario.
Thus in this study, to take into account the spatial resolution of the RCM, each box includes nine grid cells (square 3 × 3) of the RCM. Fig. 4 illustrates the climate data that will be used to feed the hydrological models described in the following Section, in particular, the climate data that will be used are:
The climate models described in the previous section represent the first element of the modeling chain. The second component is constituted by the hydrological model allowing the evaluation of soil water balance for current and future periods for a control volume in a given time interval. Generally, the main terms composing soil water balance, for a reference time, can be expressed as:
where : water storage variation (mm); P : precipitation (mm); R : runoff (mm); ET : evapotranspiration (mm); DD : deep drainage (mm).
On near time interval, additional terms can be considered. For example, during precipitation events a part can remain on the soil surface due to interception or ponding and after it is removed under evaporation or delayed infiltration. Moreover, on yearly scale, the variation in soil water storage on average tends to be small.
In this work, two hydrological models, TERRA_ML (ML means multi-layer) and HELP (Hydrologic Evaluation of Landfill Performance), are adopted and tested. Both models work on the strong assumption of 1D conditions, hence neglecting actual geomorphology or lateral superficial lost within soil intakes or leakages. In the following, they are briefly described.
TERRA_ML is the module currently adopted in the regional climate model COSMO_CLM (and also for its meteorological version: COSMO-LM) to assess the soil surface conditions in terms of water mass and temperature, constituting the lower boundaries for the atmospheric part (in this way, not only the effects of weather forcing on the soil are considered but also the corresponding soil-atmosphere feedbacks). To this aim, the assessment of the soil water/heat hydrological balance is required.
In the model, the water mass equation is linked to that for soil heat balance through the evaporation term regulating the latent heat flux. Furthermore, also the effect of snow cover assumed as a further upper layer is taken into account. The model works according to the time steps adopted for the atmospheric part.
The code HELP (Schroeder et al., 1994 ) provides a simplified modeling of surface (such as snowmelt, interception by vegetation, surface runoff and surface evaporation) and subsurface processes (evaporation from soil profile, plant transpiration, unsaturated vertical drainage) returning at a daily scale the main components of the soil water balance at surface and within the shallower soil layer.
In order to make the soil water balance estimation obtained through TERRA_ML and that obtained through HELP comparable, the vertical dimension of the domain, the soil type (and then the thermal–hydraulic properties), the land cover and the bottom conditions (free drainage) are equal.
Regarding the soil column dimension, TERRA_ML adopts two different discretization schemes for soil heat and water budget. In the first case, the depth of the domain is 21.87 m, in the second case it is 2.43 m. In particular, soil thermal behavior is required to be assessed to provide surface temperatures constituting a lower boundary for atmospheric part of climate models. Specifically, it (with required soil thermal properties) contributes to assessment of evapotranspiration component (Dickinson, 1984 ). However, in HELP solving the water budget in isothermal conditions, the domain has been set equal to 2.43 m; indeed, HELP estimates evapotranspiration component through literature models (Ritchie, 1972 ; Penman, 1963 ) neglecting soil thermal behavior.
The predominant soil type and the land cover for the five reference boxes (Table 1 ) are retrieved by the datasets implemented in TERRA_ML, respectively Digital Soil Map of the World Food and Agriculture Organization and GLC2000 database (Global Land Cover, 2000; European Commission, Joint Research Centre). Related soil-thermo hydraulic parameters and land cover features (roughness length z0 , depth, min and max plant cover, min and max LAI) are borrowed from the TERRA_ML databases (Doms et al., 2011 ); however, as the HELP model also requires a value for curve number, such value is retrieved based on the soil texture (Table 1 ). Furthermore, as large uncertainties could affect the proper estimation of soil parameters as only function of the soil texture, the variations induced by adopting different databases for soil parameters have been investigated, considering for HELP also the USDA database (included in the software, values are reported in parentheses in Table 1 ) for simulations on current period 1981–2010, driven by raw climate simulations.
|Box||Climate zone||Soil type||Land cover||ωPV||ωFC||ωPWP||ωADP||ksat (cm h−1 )||Ik2 (kg m−2 s−1 )||D0 (cm2 h−1 )||CN|
|Box 1||BWk||Sand||Bare soil||0.364 (0.437)||0.196 (0.062)||0.042 (0.024)||0.012 (−)||17.24 (20.88)||0.0035 (−)||662.40 (−)||− (66)|
|Box 2||ET||Loam||Herbaceous vegetation closed/open||0.455 (0.463)||0.340 (0.232)||0.110 (0.116)||0.035 (−)||1.91 (1.35)||0.001 (−)||128.52 (−)||− (68)|
|Box 3||Cfa||Clay||Cultivated and managed||0.507 (0.475)||0.463 (0.378)||0.257 (0.265)||0.065 (−)||0.006 (0.06)||0.0001 (−)||15.19 (−)||− (85)|
|Box 4||Cfa||Loamy clay||Cultivated and managed||0.475 (0.464)||0.370 (0.310)||0.185 (0.187)||0.060 (−)||0.28 (0.23)||0.0006 (−)||42.48 (−)||− (81)|
|Box 5||Cwa||Loam||Cultivated and managed||0.455 (0.463)||0.340 (0.232)||0.110 (0.116)||0.035 (−)||1.91 (1.35)||0.001 (−)||128.52 (−)||− (53)|
ωPV = saturated water content (volume of voids); ωFC = field capacity; ωPWP = permanent wilting point; ωADP = residual water content (air dryness point); ksat = saturated hydraulic conductivity; Ik2 = minimum infiltration rate; D0 = hydraulic diffusivity parameter; CN = curve number. The parameters ωPV , ωFC , ωPWP , ωADP and CN are dimensionless.
Based on the elements mentioned above, five modeling chains are considered; they are summarized in Fig. 5 where they are identified through capital letters as follows.
Schematic representation of the five modeling chains adopted.
Each modeling chain is tested for the five boxes over 1981–2010. Chains B, C and E are adopted to investigate the climate change effects on the soil water balance for the near and long time horizon under RCP4.5 and RCP8.5.
Fig. 6 shows, in the upper part, the percentages of average yearly surface runoff and deep drainage against cumulative precipitation and, in the lower part, the percentage of average yearly actual evapotranspiration against cumulative potential evapotranspiration. Results are carried out, for the five modeling chains and for each reference box, over 1981–2010. Furthermore, the cumulative values of atmospheric inward (precipitation) and outward (potential evapotranspiration) flows are marked; however, lacking observed data, comparisons can be only in relative terms between the chains.
In the upper part, percentages of average yearly surface runoff and deep drainage against the precipitation (total amount indicated in mm per year); in the lower part percentages of average actual evapotranspiration against potential evapotranspiration (total amount indicated in mm per year). Results refer to 1981–2010.
Firstly, results provided by modeling chain A driven by the observed data are considered. In box 1 the entire precipitation (however low) is assumed to fully enter the soil while only a tiny percentage of atmospheric evaporative demand (<10%) is fulfilled, probably due to low precipitation and coarse grained soils (sand). In box 2, the deep drainage takes values close to 20% while the actual evaporation covers about 75% of the atmospheric demand. The two values could be explained recalling, in the first case, the coarse soil texture that eases the water movement within the soil pores and, in the second one, the snowpack covering for a large part of the year (average temperature about −3 °C) the surface and that is able to satisfy the (albeit small) evaporative demand. In box 3, the clay soils (ksat = 0.006 cm h−1 ) induce runoff percentages higher than 60% (about 1100 mm) and quite similar values of actual evapotranspiration (about 650 mm) maintaining a large part of water in the shallower layer available to return in atmosphere through evapotranspiration. In the same climate areas and for similar weather forcing, but for more pervious (about 50 times) loamy clay soils (box 4), the runoff decreases to 10% with a considerable deep drainage (about 35%, 600 mm) while the actual evapotranspiration exceeds 900 mm (83%). Finally, for reference box 5 where more pervious loam soil is detected, surface runoff is assumed to be zero while about 10% of the precipitation reaches the deeper layers through bottom drainage where actual evapotranspiration is estimated as 70% of potential demand.
Very similar results are returned when the same hydrological tool is forced by the bias corrected weather data (modeling chain E). It reflects the capabilities of quantile mapping approaches in totally removing biases in the weather inputs provided by climate models. On the other hand, as the main components of soil water balance (e.g. actual evapotranspiration or infiltration) are a strong function of antecedent soil water conditions, such results display also how the climate modeling chain (in spite of quantitative biases) is able to satisfactorily reproduce weather patterns (e.g. temperature seasonal cycles, number of wet days, dry interarrival periods).
The comparison between values provided by E and C allows to evaluate the effects of weather inputs (either bias corrected or not) on the soil hydrological balance. Differences are small for box 1. For box 2, the bias in precipitation (+15%) entails a higher deep drainage for relatively pervious loam soil but comparable actual evaporation. For clay soil in box 3, the relative percentage between the surface runoff and overall precipitation is maintained in both analyses (about 64%) despite the relevant differences in average cumulative precipitation (1301 mm vs 1784 mm) while the actual evapotranspiration rate is greatly diminished (−20%), probably due to lower precipitation inputs sustaining the evaporative demand. For loamy clay (box 4), the lower water input results in smaller deep drainage. Finally, for box 5, the components of simulated and modeled soil water balance are similar since average water input and atmospheric demand were comparable. In general terms, biases in precipitation induce higher variations in all components of the soil water balance, compared to those associated to temperature, i.e. the main forcing of potential evapotranspiration. In this perspective, for example, the variations mentioned above in box 3 appear to be meaningful. On the other hand, the potential evapotranspiration estimated by the raw climate data range in between 92% and 97% of that computed through modeling chain A for boxes where the underestimation of the temperature stood at about 1 °C, while the potential discrepancy does not exceed 5% for boxes where the temperature is overestimated on average (up to 2.5 °C in box 3).
The impact of the adopted physical parameterizations in hydrological modules is taken into account considering values from C and B modeling chains. In box 2, for example, TERRA_ML returns a non-negligible surface runoff component (11%) not estimated through the CN approach used in HELP, where the deep drainage rate reaches values higher than 30%; moreover, the former estimates actual evapotranspirative losses covering 75% of the atmospheric demand while, in the second case it totaled 62%. For climate areas in the SE part of the country, while the models return similar results for clay soil (box 3), again they tend to diverge in reproducing the surface runoff components for loamy clay (box 4): 40% versus 10%. The higher infiltrated rate entails greater actual evaporation and deep drainage components. Analogous findings are then retrieved also for box 5, where the surface runoff reaches 24% in TERRA_ML and 1% for HELP model. Of course, the “intermediate” soils (in term of hydraulic properties) appear to be more sensitive to the approach adopted, while coarse (fine) grained soils for which the overwhelming part of precipitation infiltrate (runoff), are less affected.
At last, modeling chains C and D differ for the soil parameter datasets adopted. As it is clear from Fig. 6 and Table 2 , beyond the differences in intrinsic soil parameters accounting for porosity or field content, the higher the differences in soil permeability, the greater the variations in soil hydrological balance. For this reason, higher spreads are detectable for clay soil, due to a difference of one order of magnitude in ksat (box 3); the higher hydraulic conductivity suggested by USDA for clay soil induces a much lower estimated surface runoff (<10%) but, at the same time, higher evaporative losses and deep drainage components.
|Period||Chain||P (mm)||T (°C)||AE (mm)||R (mm)||DD (mm)|
|Box 1 1981–2010||E||46.6||12.6||54.0||0.0||0.0|
|Box 1 RCP4.5 2021–2050||E||1.4||1.2||1.3||0.0||0.0|
|Box 1 RCP4.5 2071–2100||E||1.6||2.1||1.5||0.0||0.0|
|Box 1 RCP8.5 2021–2050||E||1.6||1.4||1.5||0.0||0.0|
|Box 1 RCP8.5 2071–2100||E||1.5||4.1||1.4||0.0||0.0|
|Box 2 1981–2010||E||506.9||−2.9||393.8||1.9||111.6|
|Box 2 RCP4.5 2021–2050||E||1.1||1.1||1.1||1.7||1.1|
|Box 2 RCP4.5 2071–2100||E||1.0||2.7||1.2||0.1||0.6|
|Box 2 RCP8.5 2021–2050||E||1.1||1.3||1.1||1.4||1.1|
|Box 2 RCP8.5 2071–2100||E||1.0||5.4||1.1||0.1||0.7|
|Box 3 1981–2010||E||1784.4||15.8||644.1||1144.3||0.0|
|Box 3 RCP4.5 2021–2050||E||1.0||0.8||1.0||1.0||0.0|
|Box 3 RCP4.5 2071–2100||E||1.1||1.5||1.0||1.1||0.0|
|Box 3 RCP8.5 2021–2050||E||1.0||1.0||1.0||1.0||0.0|
|Box 3 RCP8.5 2071–2100||E||1.0||3.1||0.9||1.1||0.0|
|Box 4 1981–2010||E||1779.3||20.3||964.6||180.6||633.1|
|Box 4 RCP4.5 2021–2050||E||1.0||1.0||1.2||1.2||0.8|
|Box 4 RCP4.5 2071–2100||E||1.0||2.1||1.0||1.0||0.9|
|Box 4 RCP8.5 2021–2050||E||1.0||1.2||1.0||1.2||1.0|
|Box 4 RCP8.5 2071–2100||E||0.9||3.8||1.1||1.0||0.7|
|Box 5 1981–2010||E||745.7||13.6||708||0.7||34.7|
|Box 5 RCP4.5 2021–2050||E||1.0||1.1||1.0||2.5||0.7|
|Box 5 RCP4.5 2071–2100||E||0.9||2.0||1.0||0.1||0.0|
|Box 5 RCP8.5 2021–2050||E||1.0||1.2||1.0||1.7||1.1|
|Box 5 RCP8.5 2071–2100||E||1.2||3.3||1.1||1.2||1.5|
The discussed values are reported for modeling chains B, C and E in Table 2 . In the same Table, future anomalies in main weather forcing (T and P) and principal components of the soil water budget are also given. For box 1 (NW), some issues for the bias correction approach in keeping the climate signal retrieved from climate models are detected not only for precipitations for which the cumulative values are very small, but also for mean temperatures (e.g. −1 °C for RCP8.5 in 2071–2100). On the other hand, much smaller differences between raw and bias corrected data arise if the temperature climate signal is evaluated in terms of relative variation; in this case, the discrepancies do not exceed 3%. However, due to the higher relative availability of soil water, a greater increase (in percentage) in actual evapotranspiration is returned using bias corrected weather inputs. On the other hand, for box 2 (SW), raw and bias corrected values are quite similar and the climate signal is wholly analogous; however, in this case, the discrepancies in hydrological modeling involve significant spreads in anomalies. For example, the actual evaporation is expected to be increasing (+10%) using HELP (driven by raw or bias corrected) while it is fluctuating for TERRA_ML. Such differences are much more evident for surface runoff and deep drainage components: in the first case, the abrupt variations between near and long time horizons using modeling chain E, induced by low runoff contributions, are not detected in the other analysis; for the deep drainage, TERRA_ML returns consistently a growth (up to 60%) between near and long time horizon, while through HELP, a sign variation is also detected (averagely from about 10% to about −20%). In this regard, it is worth noting that the average temperature values in box 2 are currently lower than zero while, under the effect of climate change, on the long time horizon they could move towards positive values, changing dramatically the snow duration and depth values and then with large variations of the whole soil surface water balance. For box 3 and box 4 (SE), in spite of the great differences in cumulative precipitation between raw and bias corrected values, the climate signal is maintained in a proper way, while for the temperature a decrease is evaluated (higher for box 3). Beyond such variations, components of the soil water balance are expected not to be affected in a significant way; a slight decrease in actual evapotranspiration (about 10% for the long time horizon) and increase in surface runoff (up to 10% for box 3 and 20% for box 4) are estimated while the deep drainage for box 4 is assessed as decreasing mainly on the long time horizon. For these boxes, it clearly arises how the compensation and smoothing effects could make the steps from the variations in weather forcing to the corresponding main components of the soil water balance non-trivial. Indeed, considering for example the farthest time horizons, the increases in temperature (in relative terms, 10% under RCP4.5 and 20% under RCP8.5) do not result in variations of equal magnitude in the actual evapotranspiration (which are often unchanged or slightly decreasing). On the other side, for box 4, under a quite similar cumulative precipitation amount, more consistent variations are returned in surface runoff and deep drainage probably due to variations in precipitation patterns (less wet days but characterized by higher intensities). At last, for box 5 (NE), remarks quite similar to those carried out for box 2 can be reiterated: in both cases, probably, the differences in reproducing surface runoff are associated to the features assumed for loam soil (more pervious according HELP modeling and then returning a greater infiltration rate), clearly revealing the crucial role and uncertainties potentially associated to soil parameters.
The overall effects of such variations are summed up in Fig. 7 . In the first two columns, one for each RCP considered, the seasonal cycles of cumulative soil water storage are displayed for the five boxes on three different time spans, control period 1981–2010, near term 2021–2050 and long term 2071–2100. The cumulative soil water storage is computed as the algebraic sum of cumulative values of water budget components from the initial instant t0 , corresponding to the beginning of the hydrological year (1st September) to the first day of the i -th month. Furthermore, in the third column, the relative frequencies of occurrence for daily soil moisture, ranging from the wilting point to saturation are displayed for the five boxes on 1981–2010 and 2071–2100 under both scenarios adopting modeling chain E. Concerning the RCP4.5 scenario, a general decrease of water availability in shallower layers is detected for all boxes except box 3 (more evident on the long time horizon); again, it is worth noting how grain texture can affect the soil response comparing cycles for box 3 and box 4 subjected to similar variations in weather forcing: while in clayey soils small variations arise, in the more pervious loamy clay a large reduction is estimated occurring (mainly on 2071–2100) for most of the year. The greater variations are localized in drier seasons (box 2, 4 and 5). Regarding RCP8.5, a minor decrease is expected in box 2 against a greater increase in temperature but where snow hydrology dynamics could play a relevant role. A slight decrease is estimated for box 3 while for box 4 (loamy clay) the variations between wet and dry seasons appear to be clearer. Finally, for box 5 the increase in precipitation is expected to partly counteract the correspondent higher actual evaporation, resulting in a less pronounced reduction of soil water storage during the entire year. In addition, the trends of the cumulative water storage curves point out clearly how larger variations in weather forcing could not result in larger differences in soil moisture due to the compensation and smoothing effects. In this regard, the land cover and soil act as a filter reducing the amplitude of variations of forcing at surface (precipitation and evapotranspiration demand). Finally, in terms of relative frequencies of occurrence of daily soil moisture, the distributions related to future time spans are characterized by a “backward shift” in the peak value in some cases (mainly for boxes 3 and 4) and generally by a higher data dispersion, consistent with what was shown above.
Seasonal cycles of cumulative soil water storage on two different time spans (first two columns); relative frequencies of occurrence for daily soil moisture ranging from wilting point to saturation on 1981–2010 and 2071–2100 (third column; RCP4.5 and RCP8.5).
The present work is aimed to investigate the main variations assessed in soil water budget for the main climate zones identified in China under the effect of climate change, and to quantify the reliability of the results under climate and hydrological modeling uncertainties. For this purpose, climate simulations at a high resolution (about 14 km) on a large part of China (Bucchignani et al., 2014 ) are used as weather forcing; furthermore, the adoption of two hydrological models, two soil parameter datasets, observed or modeled inputs (provided by raw or bias corrected climate simulations) result in five different simulation chains. The analyses are carried out on five boxes selected in the four main climate areas detected in the investigated domain; however, they should not be considered as strictly related to their location but as examples of potential future dynamics occurring in the different climate zones. Indeed, it is worth remembering that, generally, the soil datasets return parameters applicable at a global scale and estimated from the average features; the local geomorphological peculiarities could induce remarkable differences, in grain size and mainly in soil properties. In this perspective, this work is meant to display how the future impacts induced by variations of weather patterns could vary for different soils and land covers without expecting to represent the actual performance of current and future water balance in the boxes. Concerning the issue of investigating the role played by hydrological modeling and parameters in the evaluation of soil water budget components and their variation under the effect of climate change, the analysis highlights how the choice of some approaches or soil parameterizations can affect the results both in absolute and in relative terms (e.g., as variations between different time spans); moreover, the research shows how the differences could as well be highly related to weather forcing in inputs or investigated soil. Of course, the availability of observed data related to such components could allow for the evaluation of performances and then the selection between the different models; however, at the moment, the available datasets are few and often not adequate. Under these constraints, as currently carried out for climate models, the solution could be the adoption of ensembles multi or single-model with different soil parameter datasets. Moreover, as the processes regulating the soil water budget result as highly non-linear, the adoption of raw or bias corrected inputs can produce assessments different for magnitude and sign. Furthermore, the comparison between the analyses driven by observed and bias corrected values allows to appreciate the great capabilities in removing biases assumed as systematic in climate simulations. On the other hand, the further assumptions related to their use and the possible induced variations in the climate signal with respect to those provided by physically based climate modeling chains should be kept in mind.
Finally, the most evident findings retrieved by the analysis concern the substantial expected decrease in average water content in the shallower part of the soil with different extents according to climate zone, concentration scenario and soil/cover features. Moreover, the assessed increase in temperature and then in evapotranspirative demand do not ever result in higher actual evapotranspiration values, due to the concurrent variations in precipitation patterns. Despite the displayed current uncertainties, the results could represent a further tool in attempting to identify the most vulnerable future areas and what might be the priorities in future adaptation policies. To this aim, the adoption of a single climate simulation could represent a strong constraint while the ensemble simulation could be preferred.
Moreover, the performed analyses are essentially aimed to estimate variations in impacts on the soil, produced by changes in the average values of weather forcing due to climate change that could affect relevant sectors, as agriculture or water supply. On the other hand, special attention should be paid also to the effects of changes in extreme weather events that could entail, in different areas, a substantial increase in occurrence and magnitude of hydro meteorological hazards.
This work was developed within the framework of the GEMINA project, Work Package 7.1.6, “B action” (Italy–China cooperation on climate changes), funded by the Italian Ministry of Education, University, and Research and the Italian Ministry of the Environment, Land, and Sea .