hydromt_wflow.workflows.soilgrids#

hydromt_wflow.workflows.soilgrids(ds: ~xarray.core.dataset.Dataset, ds_like: ~xarray.core.dataset.Dataset, ptfKsatVer: str = 'brakensiek', soil_fn: str = 'soilgrids', wflow_layers: ~typing.List[int] = [100, 300, 800], logger=<Logger hydromt_wflow.workflows.soilgrids (WARNING)>)[source]#

Return soil parameter maps at model resolution.

Based on soil properties from SoilGrids datasets. Both soilgrids 2017 and 2020 are supported. Soilgrids 2017 provides soil properties at 7 specific depths, while soilgrids_2020 provides soil properties averaged over 6 depth intervals. Ref: Hengl, T., Mendes de Jesus, J., Heuvelink, G. B. M., Ruiperez Gonzalez, M., Kilibarda, M., Blagotic, A., et al.: SoilGrids250m: Global gridded soil information based on machine learning, PLoS ONE, 12, https://doi.org/10.1371/journal.pone.0169748, 2017. Ref: de Sousa, L.M., Poggio, L., Batjes, N.H., Heuvelink, G., Kempen, B., Riberio, E. and Rossiter, D., 2020. SoilGrids 2.0: producing quality-assessed soil information for the globe. SOIL Discussions, pp.1-37. https://doi.org/10.5194/soil-2020-65.

A soil_mapping table can optionally be provided to derive parameters based on soil texture classes. A default table soil_mapping_default is available to derive the infiltration capacity of the soil.

The following soil parameter maps are calculated:

  • theta_s : average saturated soil water content [m3/m3]

  • theta_r : average residual water content [m3/m3]

  • **ksat_vertical ** : vertical saturated hydraulic conductivity at soil

surface [mm/day] - soil_thickness : soil thickness [mm] - f : scaling parameter controlling the decline of ksat_vertical [mm-1] (fitted with curve_fit (scipy.optimize)), bounds are checked - soil_f_ : scaling parameter controlling the decline of ksat_vertical [mm-1]

(fitted with numpy linalg regression), bounds are checked

  • soil_brooks_corey_c_ map: Brooks Corey coefficients [-] based on pore size distribution index for the wflow_sbm soil layers.

  • meta_{soil_fn}_ksat_vertical_[z]cm : ksat vertical [mm/day] at soil depths [z] of SoilGrids data [0.0, 5.0, 15.0, 30.0, 60.0, 100.0, 200.0]

  • meta_soil_texture : USDA Soil texture based on percentage clay, silt,

sand mapping: [1:Clay, 2:Silty Clay, 3:Silty Clay-Loam, 4:Sandy Clay, 5:Sandy Clay-Loam, 6:Clay-Loam, 7:Silt, 8:Silt-Loam, 9:Loam, 10:Sand, 11: Loamy Sand, 12:Sandy Loam]

Parameters:

  • ds (xarray.Dataset) – Dataset containing soil properties.

  • ds_like (xarray.DataArray) – Dataset at model resolution.

  • ptfKsatVer (str) – PTF to use for calculation ksat_vertical .

  • soil_fn (str) – soilgrids version {‘soilgrids’, ‘soilgrids_2020’}

  • wflow_layers (list) – List of soil layer depths [cm] for which soil_brooks_corey_c is calculated.

Returns:

ds_out – Dataset containing gridded soil parameters.

Return type:

xarray.Dataset