hydromt_wflow.wflow.WflowModel.setup_soilmaps

WflowModel.setup_soilmaps(soil_fn='soilgrids', ptf_ksatver='brakensiek')[source]

This component derives several (layered) soil parameters based on a database with physical soil properties using available point-scale (pedo)transfer functions (PTFs) from literature with upscaling rules to ensure flux matching across scales.

Currently, the only supported soil_fn is “soilgrids” and the default ptf_ksatver (PTF for the vertical hydraulic conductivity) is “brakensiek”.

The following maps are added to staticmaps:

  • thetaS map: average saturated soil water content [m3/m3]

  • thetaR map: average residual water content [m3/m3]

  • KsatVer map: vertical saturated hydraulic conductivity at soil surface [mm/day]

  • SoilThickness map: soil thickness [mm]

  • SoilMinThickness map: minimum soil thickness [mm] (equal to SoilThickness)

  • M map: model parameter [mm] that controls exponential decline of KsatVer with soil depth

    (fitted with curve_fit (scipy.optimize)), bounds of M are checked

  • M_ map: model parameter [mm] that controls exponential decline of KsatVer with soil depth

    (fitted with numpy linalg regression), bounds of M_ are checked

  • M_original map: M without checking bounds

  • M_original_ map: M_ without checking bounds

  • f map: scaling parameter controlling the decline of KsatVer [mm-1]

  • c_0 map: Brooks Corey coefficient [-] based on pore size distribution index at

    depth of 1st soil layer (100 mm) wflow_sbm

  • c_1 map: idem c_0 at depth 2nd soil layer (400 mm) wflow_sbm

  • c_2 map: idem c_0 at depth 3rd soil layer (1200 mm) wflow_sbm

  • c_3 map: idem c_0 at depth 4th soil layer (> 1200 mm) wflow_sbm

  • KsatVer_[z]cm map: KsatVer [mm/day] at soil depths [z] of SoilGrids data [0.0, 5.0, 15.0, 30.0, 60.0, 100.0, 200.0]

Parameters

  • soil_fn ({'soilgrids'}) –

    Name of data source for soil parameter maps, see data/data_sources.yml. Should contain info for the 7 soil depths of soilgrids.

    • Required variables: [‘bd_sl*’, ‘clyppt_sl*’, ‘sltppt_sl*’, ‘oc_sl*’, ‘ph_sl*’, ‘sndppt_sl*’, ‘soilthickness’, ‘tax_usda’]

  • ptf_ksatver ({'brakensiek', 'cosby'}) – Pedotransfer function (PTF) to use for calculation KsatVer (vertical saturated hydraulic conductivity [mm/day]). By default ‘brakensiek’.