Simulate soil carbon decomposition rates, calculate nutrient depletion indices, and model land degradation erosion variables.
Soil carbon storage is modeled using multi-compartment kinetics (such as the RothC model). Carbon inputs decompose through pools (active, slow, passive) regulated by local clay fractions, hourly temperatures (ERA5), and moisture indices (SAVI).
Where Y is residual soil carbon pool weight, Y_0 is initial baseline carbon, k is decay constant, and a is local moisture/temperature factor.
Macronutrient states (Nitrogen, Phosphorus, Potassium) are modeled across space using multi-spectral reflectance regression models. Spectral Red-Edge and SWIR bands correlate to canopy chlorophyll and biomass states, providing proxy markers of nitrogen reserves.
Erosion rates are solved using the Universal Soil Loss Equation (USLE), combining terrain slope/aspect with soil erodibility indexes and satellite vegetation indices (like NDVI) to calculate pixel-by-pixel degradation values.
Where A is computed annual soil loss, R is rainfall erosivity, K is soil erodibility factor, LS is slope aspect factor, and C, P are land management variables.