RothC

satterc.pipeline.models.rothc

evaporation_monthly

evaporation_monthly(actual_evapotranspiration_monthly: DataArray) -> DataArray

Extract evaporation data for RothC model.

Parameters:

actual_evapotranspiration_monthly (DataArray) –

Monthly actual evapotranspiration (mm).

Returns:

DataArray –

Monthly evaporation data.

farmyard_manure_input_monthly

farmyard_manure_input_monthly(
    plant_type: DataArray, dates_monthly: DatetimeIndex
) -> DataArray

For now, return array of zeros for farmyard manure input.

In future, could be determined by pft (non-zero if crop) and month of year.

inert_organic_matter

inert_organic_matter(organic_carbon_stocks: DataArray) -> DataArray

Calculate inert organic matter from organic carbon stocks.

Parameters:

organic_carbon_stocks (DataArray) –

Organic carbon stocks (tC/ha).

Returns:

DataArray –

Inert organic matter (tC/ha).

plant_cover_monthly

plant_cover_monthly(plant_type: DataArray, dates_monthly: DatetimeIndex) -> DataArray

Temporary bridge to the boolean plant cover data required by RothC.

Just returns an array of ones with shape (n_months, n_pixels).

rothc

rothc(
    temperature_celcius_monthly: DataArray,
    precipitation_mm_monthly: DataArray,
    evaporation_monthly: DataArray,
    plant_cover_monthly: DataArray,
    dpm_rpm_ratio_monthly: DataArray,
    soil_carbon_input_monthly: DataArray,
    farmyard_manure_input_monthly: DataArray,
    clay_content: DataArray,
    inert_organic_matter: DataArray,
    soil_depth: DataArray,
    dates_monthly: Index,
    rothc_parameters: tuple[int],
) -> dict[str, DataArray]

Rothamsted Carbon model.

Monthly resolution input data.

Parameters:

temperature_celcius_monthly (DataArray) –

Monthly mean temperature in degrees Celsius.
precipitation_mm_monthly (DataArray) –

Monthly precipitation in mm.
evaporation_monthly (DataArray) –

Monthly evaporation in mm.
plant_cover_monthly (DataArray) –

Monthly plant cover as boolean (True = covered).
dpm_rpm_ratio_monthly (DataArray) –

Ratio of decomposable to resistant plant material.
soil_carbon_input_monthly (DataArray) –

Carbon input in tC/ha/month.
farmyard_manure_input_monthly (DataArray) –

Farmyard manure input in tC/ha/month.
clay_content (DataArray) –

Clay content percentage.
soil_depth (DataArray) –

Soil depth in cm.
inert_organic_matter (DataArray) –

Inert organic matter in tC/ha.
rothc_parameters (tuple[int]) –

Tuple of parameters.

Returns:

dict –

Dictionary containing monthly model outputs: - decomposable_plant_material_monthly: DPM pool (tC/ha) - resistant_plant_material_monthly: RPM pool (tC/ha) - microbial_biomass_monthly: Microbial biomass pool (tC/ha) - humified_organic_matter_monthly: HUM pool (tC/ha) - soil_organic_carbon_monthly: Total SOC (tC/ha)

Notes

All outputs have units tC/ha (tonnes of Carbon per hectare). All outputs are at monthly resolution.

rothc_parameters

rothc_parameters(n_years_spinup: int = 1) -> tuple[int]

Static parameters for the Rothamsted Carbon model.

Parameters:

n_years_spinup (int, default: 1 ) –

Number of years to use for model spin-up.

Returns:

Tuple containing these parameters. –

soil_carbon_input_monthly

soil_carbon_input_monthly(litter_pool_monthly: DataArray) -> DataArray

Temporary bridge to map litter input to soil carbon input.

Parameters:

litter_pool_monthly (DataArray) –

Monthly litter input (tC/ha).

Returns:

DataArray –

Monthly soil carbon input (tC/ha).