Soil

This module computes soil storage and related fluxes for all grid cells based on section 4.4 of Müller Schmied et al 2021 [1].

soil.soil_water_balance(soil_water_content, pet_to_soil, current_landarea_frac, landareafrac_ratio, max_temp_elev, canopy_evap, effective_precipitation, precipitation, immediate_runoff, land_storage_change_sum, sublimation, daily_storage_transfer, snow_freeze_temp, gamma, max_daily_pet, humid_arid, soil_texture, drainage_direction, max_groundwater_recharge, groundwater_recharge_factor, critcal_gw_precipitation, max_soil_water_content, areal_corr_factor, minstorage_volume, x, y)

Compute daily soil balance.

Parameters:

soil_water_contentfloat

Soil water content, Units: [mm]

pet_to_soilfloat

Remaining energy for addtional evaporation from soil, Units: [mm]

current_landarea_fracfloat

Land area fraction of current time step, Units: [-]

landareafrac_ratiofloat

Ratio of land area fraction of previous to current time step,

Units: [-]

max_temp_elevfloat

Maximum temperature from the 1st(lowest) elevation from snow algorithm. , Units: [K]

canopy_evapfloat

Canopy evaporation, Units: [mm/day]

effective_precipitationfloat

Effective precipitation based on Müller Schmied et al 2021, Units: [mm/day]

precipitationfloat

Daily precipitation, Units: [mm/day]

immediate_runofffloat

Immediate runoff, Units: [mm/day]

land_storage_change_sumfloat

Sum of change in vertical balance storages, Units: [mm]

sublimationfloat

Sublimation, Units: [mm/day]

daily_storage_transferfloat

Storage to be transfered to runoff when land area fraction of current time step is zero, Units: [mm]

snow_freeze_tempfloat

Snow freeze temperature , Units: [K]

gammafloat

Runoff coefficient , Units: [-]

max_daily_petfloat

Maximum daily potential evapotranspiration, Units: [mm/day]

humid_aridfloat

Humid-arid calssification based on Müller Schmied et al. 2021

soil_texturefloat

Soil texture classification based on Müller Schmied et al. 2021

drainage_directionfloat

Drainage direction based on Müller Schmied et al. 2021

max_groundwater_rechargefloat

Maximum groundwater recharge from soil, Units: [mm/day]

groundwater_recharge_factorfloat

Groundwater recharge factor, Units: [-]

critcal_gw_precipitationfloat

critical precipitation for groundwater recharge, Units: [mm/day]

max_soil_water_contentfloat

Maximum soil water content , Units: [mm]

areal_corr_factorfloat

Areal correction factor-CFA for correcting runoff, Units: [-]

minstorage_volumefloat

Volume at which storage is set to zero, Units: [km3]

x, y : Latititude and longitude indexes of grid cells.

Returns:

soil_water_contentfloat

Updated soil water content , Units: [mm]

groundwater_recharge_from_soil_mmfloat

Groundwater recharge from soil, Units: [mm/day]

actual_soil_evapfloat

Actual evapotranspiration from the soil, Units: [mm/day]

soil_saturationfloat

Soil saturation, Units: [-]

surface_runofffloat

Surface runoff from land, Units: [mm/day]

daily_storage_transferfloat

Updated storage to be transfered to runoff when land area fraction of current time step is zero, Units: [mm]

total_daily_runofffloat

Total daily runoff from land (RL) (runoff + immediate runoff + soil water overflow), Units: [mm/day]

daily_runofffloat

daily runoff (R3), Bergström (1995), Units: [mm/day]

soil_water_overflowfloat

Soil water overflow (R2), Units: [mm/day]

immediate_runofffloat

runoff from urban areas or immediate runoff (R1), Units: [mm/day]

Note

The computation order for the soil storage module is as follows: First, immediate runoff (R1) is calculated. After, effective precipitation is reduced by the immediate runoff. Then, runoff from soil water overflow (R2) is computed. After, daily runoff (R3) is calculated followed by actual evapotranspiration. Soil storage is updated. Afterwards, ground water recharge is calculated based on daily runoff. Then, total daily runoff from land (RL) is computed as daily runoff (R3) + immediate runoff (R1) + soil water overflow (R2).

Note: Total daily runoff from land is corrected with an areal correction factor (CFA) (if gamma is insufficient to fit simulated discharge). To conserve water balance, actual evapotranspiration is also corrected with CFA. After evapotranspiration correction, soil storage, total daily runoff from land and groundwater recharge are adjusted as well. Finally, Surface runoff \(({R}_{s})\) is calculated as total daily runoff from land minus groundwater recharge.

Water balance

Soil storage \(S_s\) \([mm]\) is calculated as:

\[\frac{dS_s}{d_t} = {P}_{eff} − {R}_{l}− {E}_{s}\]

where \({P}_{eff}\) is effective precipitation \([mm/d]\), \({R}_{l}\) is total runoff from land \([mm/d]\) and \({E}_{s}\) is the actual evapotranspiration from soil \([mm/d]\).

Inflows

Effective precipitation \({P}_{eff}\) is calculated as:

\[{P}_{eff} = {P}_{t} − {P}_{sn} + {M}\]

where \({P}_{t}\) is throughfall \([mm/d]\), \({P}_{sn}\) is snowfall \([mm/d]\) and \({M}\) is snow melt \([mm/d]\).

Actual evapotranspiration \({E}_{s}\) from soil \([mm/d]\) is calculated as:

\[{E}_{s} = min\biggl(({E}_{pot} - E_c) , ({E}_{pot,max} - E_c) \times \frac{S_s}{S_s,max} \biggr)\]

where \({E}_{pot}\) is potential evapotranspiration \([mm/d]\), \({E}_{c}\) is canopy evaporation \([mm/d]\) and \({S}_{s,max}\) is the maximum soil water content \([mm]\) derived as a product of total available water capacity in the upper meter of the soil [2] and land-cover-specific rooting depth (Table C2 [1]). The maximum potential evapotranspiration \({E}_{pot,max}\) is set to \(15 {mm}/{d}\) globally.

Total daily runoff from land (RL) is calculated as:

\[RL = R1 + R3 + R2\]

where immediate runoff from urban areas (R1) is computed as:

\[{immediate \: runoff} = 0.5 \times {P}_{eff} \times fraction \: of \: build \: up \: area\]

and where soil water overflow (R2) is calculated as:

\[\begin{split}{P}_{eff} = \begin{cases} {P}_{eff} +{S}_{s,p} - {S_s,max}, & \text{if } ({P}_{eff} + {S}_{s,p})>{S_s}_{,max} \\ 0, & \text{otherwise} \end{cases}\end{split}\]

where \({P}_{eff}\) is effective precipitation, \({S}_{s,p}\) and \({S}_{s,max}\) is soil storage of the previous day and maximum soil storage respectively.

Daily runoff from soil (\({R3}\)) \([mm/day]\) is calculated following Bergström (1995) [3] as:

\[R3 = {P}_{eff} \biggl(\frac{S_s}{S_s,max}\biggr)^\Gamma\]

where Gamma is the runoff coefficient \([–]\). This parameter, which varies between :math`0.1` and \(5.0\), is used for calibration. Together with soil saturation, it determines the fraction of \({P}_{eff}\) that becomes \({R3}\).

Note

If the sum of \({P}_{eff}\) and \({S}_{s}\) of the previous day exceed \({S}_{s,max}\), the overflow \({R2}\) is added to daily runoff \({R3}\) and immediate runoff \({R1}\) to total daily runoff from land \({R}_{L}\).

\({R3}\) is partitioned into fast surface and subsurface runoff \({R}_{s}\) and diffuse groundwater recharge \({R}_{g}\) according to the heuristic scheme.

\[{R}_{g} = min\biggl({R}_{gmax} , ({f}_{g} \times {R3} \biggr)\]

where \({R}_{gmax}\) is soil-texture-specific maximum groundwater recharge with values of 7, 4.5 and 2.5 \([mm/d]\) for sandy, loamy and clayey soils, respectively, and \({f}_{g}\) is the groundwater recharge factor ranging between 0 and 1. \({f}_{g}\) is determined based on relief, soil texture, aquifer type, and the existence of permafrost or glaciers [4].

Note

If a grid cell is defined as (semi)arid and has coarse (sandy) soil, groundwater recharge will only occur if precipitation exceeds a critical value of 12.5 \([mm/d]\), otherwise the water remains in the soil. The fraction of \({R}_{3}\) that does not recharge the groundwater becomes \({R}_{s}\), which recharges surface water bodies and the river compartment.

If the gamma parameter is not enough to match the observed discharge, the total daily runoff from land is adjusted by multiplying it with a real correction factor (CFA). To conserve mass balance actual total evaporation from land (\({E}_{s}\)) is corrected such that when a real correction factor is increased or reduced to increase runoff, actual total evaporation will also be reduced or increased respectively.

\[\begin{split}RL &= P - E_s - \frac{dS}{dt}, (dt=1), eqn.1, \\ RL \times CFA &= P - {E}_{s,corr} - dS,eqn2,\end{split}\]

substituting \(RL\) from equation 2 into equation 1:

\[\begin{split}P - {E}_{s,corr} - ds &= CFA (P - E_s - ds), \\ {E}_{s,corr} &= dS(CFA-1) - P(CFA-1) + E_s(CFA)\end{split}\]

where \({P}\) is precipitation \([mm/day]\).

Note

Surface runoff \(({R}_{s})\) is finally calculated as total daily runoff from land minus groundwater recharge.

References

[1] (1,2)

Müller Schmied, H., Cáceres, D., Eisner, S., Flörke, M., Herbert, C., Niemann, C., Peiris, T. A., Popat, E., Portmann, F. T., Reinecke, R., Schumacher, M., Shadkam, S., Telteu, C.E., Trautmann, T., & Döll, P. (2021). The global water resources and use model WaterGAP v2.2d: model description and evaluation. Geoscientific Model Development, 14(2), 1037–1079. https://doi.org/10.5194/gmd-14-1037-2021

[2]

Batjes, N. H.: ISRIC-WISE derived soil properties on a 5 by 5 arc-minutes global grid (ver. 1.2), Tech. Rep. 2012/01, ISRIC–World Soil Information, Wageningen, 2012.

[3]

Bergström, S.: The HBV model, in: Computer models of watershed hydrology, edited by: Singh, V., Water Resources Publications, Lone Tree, USA, 443–476, 1995

[4]

Döll, P. and Fiedler, K.: Global-scale modeling of groundwater recharge, Hydrol. Earth Syst. Sci., 12, 863–885, https://doi.org/10.5194/hess-12-863-2008, 2008.