Soil#

Submodules#

Module contents#

class abcmodel.land.soil.StandardSoilModel(a=0.219, b=4.9, p=4.0, cgsat=3.56e-06, wsat=0.472, wfc=0.323, wwilt=0.171, w2=0.21, d1=0.1, c1sat=0.132, c2ref=1.8, rssoilmin=50.0)[source]#

Bases: AbstractSoilModel[StandardSoilState]

Standard soil model with moisture and temperature dynamics.

Parameters:
  • a (float) – Clapp and Hornberger (1978) retention curve parameter. Default is 0.219.

  • b (float) – Clapp and Hornberger (1978) retention curve parameter. Default is 4.90.

  • p (float) – Clapp and Hornberger (1978) retention curve parameter. Default is 4.0.

  • cgsat (float) – saturated soil heat capacity [J m-3 K-1]. Default is 3.56e-6.

  • wsat (float) – saturated soil moisture content [m3 m-3]. Default is 0.472.

  • wfc (float) – soil moisture content at field capacity [m3 m-3]. Default is 0.323.

  • wwilt (float) – soil moisture content at wilting point [m3 m-3]. Default is 0.171.

  • w2 (float) – soil moisture content at the second layer [m3 m-3]. Default is 0.21.

  • d1 (float) – depth of the top soil layer [m]. Default is 0.1.

  • c1sat (float) – saturated soil hydraulic conductivity parameter [-]. Default is 0.132.

  • c2ref (float) – reference soil hydraulic conductivity parameter [-]. Default is 1.8.

  • rssoilmin (float) – minimum soil resistance [s m-1]. Default is 50.0.

init_state(wg=0.21, temp_soil=285.0, temp2=286.0, rssoil=1000000.0)[source]#

Initialize the soil state.

Parameters:
  • wg (float) – Volumetric soil moisture [m3 m-3]. Default is 0.21.

  • temp_soil (float) – Soil temperature [K]. Default is 285.0.

  • temp2 (float) – Deep soil temperature [K]. Default is 286.0.

  • rssoil (float) – Soil resistance [s m-1]. Default is 1.0e6.

Returns:

The initialized StandardSoilState.

run(state)[source]#

Compute soil surface resistance.

compute_soil_resistance(wg)[source]#

Compute the soil resistance rssoil.

compute_temp_soil_tend(gf, temp_soil, temp2)[source]#

Compute the soil temperature tendency.

compute_wgtend(wg, le_soil)[source]#

Compute the soil moisture tendency.

run_tends(state, surf_state)[source]#

Compute soil tendencies that depend on surface fluxes.

integrate(state, dt)[source]#

Integrate soil moisture and temperature forward in time.

class abcmodel.land.soil.StandardSoilState(wg, temp_soil, temp2, rssoil=<factory>, temp_soil_tend=<factory>, wgtend=<factory>)[source]#

Bases: AbstractSoilState

Standard soil state.

wg: Array#

Soil moisture content in the root zone [m3 m-3].

temp_soil: Array#

Soil temperature [K].

temp2: Array#

Deep soil temperature [K].

rssoil: Array#

Soil resistance [m s-1].

temp_soil_tend: Array#

Soil temperature tendency [K s-1].

wgtend: Array#

Soil moisture tendency [m3 m-3 s-1].