Jarvis Stewart

Jarvis Stewart#

class abcmodel.land.biosphere.jarvis_stewart.JarvisStewartState(rs, wl, cliq, wCO2, cveg, wltend=<factory>)[source]#

Bases: AbstractBiosphereState

Jarvis-Stewart biosphere state.

rs: Array#

Surface resistance [s m-1].

wl: Array#

Canopy water content [m].

cliq: Array#

Wet fraction of canopy [-].

wCO2: Array#

Kinematic CO2 flux [mol m-2 s-1].

cveg: Array#

Vegetation fraction [-].

wltend: Array#

Canopy water content tendency [m].

class abcmodel.land.biosphere.jarvis_stewart.JarvisStewartModel(rsmin=110.0, lai=2.0, gD=0.0, cveg=0.85, wmax=0.0002, wwilt=0.171, wfc=0.323, w2=0.21)[source]#

Bases: AbstractBiosphereModel[JarvisStewartState]

Jarvis-Stewart biosphere model with empirical surface resistance.

Parameters:

  • rsmin (float) – minimum stomatal resistance [s m-1]. Default is 110.0.

  • lai (float) – leaf area index [m2 m-2]. Default is 2.0.

  • gD (float) – canopy rad extinction coefficient [-]. Default is 0.0.

  • wmax (float) – maximum water storage capacity of the canopy [m]. Default is 0.0002.

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

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

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

init_state(rs=1000000.0, wl=0.0, cliq=0.0, wCO2=0.0)[source]#

Initialize the biosphere state.

Parameters:

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

  • wl (float) – Canopy water content [m]. Default is 0.0.

  • cliq (float) – Wet fraction of canopy [-]. Default is 0.0.

  • wCO2 (float) – Kinematic CO2 flux [mol m-2 s-1]. Default is 0.0.

  • cveg – vegetation fraction [-]. Default is 0.85.

Returns:

The initialized JarvisStewartState.

run(state)[source]#

Compute biosphere surface resistance and canopy wet fraction.

compute_f1(in_srad)[source]#

Compute rad factor f1.

compute_f2(wg)[source]#

Compute soil moisture factor f2.

compute_f3(esat, e)[source]#

Compute VPD factor f3.

compute_f4(theta)[source]#

Compute temperature factor f4.

compute_cliq(wl)[source]#

Compute wet fraction of canopy cliq.

compute_wltend(le_liq)[source]#

Compute canopy water storage tendency.

run_tends(state, surf_state)[source]#

Compute biosphere tendencies that depend on surface fluxes.

integrate(state, dt)[source]#

Integrate canopy water content forward in time.