hydromt_sfincs.workflows.curve_number_with_recovery#

hydromt_sfincs.workflows.curve_number_with_recovery(da_lulc: DataArray, da_hsg: DataArray, da_ksat: DataArray, df_map: DataFrame, *, effective: float, ksat_max: float | None = 360.0, da_mask: DataArray | None = None) → Dataset[source]#

Estimate SCS curve-number-with-recovery parameters.

Curve numbers are derived from land use and HSG, converted to maximum soil retention, then combined with Ksat to calculate effective retention and recovery rate. Inputs must already be-read xarray objects and a DataFrame.

Parameters:

  • da_lulc (xr.DataArray) – Land-use or land-cover class raster.

  • da_hsg (xr.DataArray) – Hydrologic soil group class raster.

  • da_ksat (xr.DataArray) – Saturated hydraulic conductivity raster [mm/hr].

  • df_map (pd.DataFrame) – Curve-number table indexed by land-use class, with HSG classes as columns.

  • effective (float) – Fraction of smax used for effective retention, e.g. 0.5 for 50%.

  • ksat_max (float or None, optional) – Upper cap for Ksat in mm/hr; None disables the cap. By default 360.0 (100 um/s), high enough to only remove outliers.

Returns:

smax and seff in metres and ks in mm/hr. Cells without valid Ksat stay NaN in ks.

Return type:

xr.Dataset

On this page