Declining but spatially redistributed global fire nitrogen emissions regulated by vegetation-climate interactive mechanisms

The code combines field-measured emission factors (EFs), combustion factors (CFs), satellite-derived burned area (BA), and aboveground biomass (AGB / fuel load, FL) to produce gridded emission estimates with full uncertainty characterization.</p><p dir="ltr">The analysis employs two complementary uncertainty methods:</p><p dir="ltr">1. **First-order error propagation** — analytical calculation of per-pixel uncertainty and parameter contribution breakdown</p><p dir="ltr">2. **Monte Carlo simulation with Latin Hypercube Sampling (LHS)** — numerical uncertainty estimation with region-specific EF sampling, block-LHS spatial correlation structure, and freeze-one contribution decomposition</p><p dir="ltr">The code supports seven biomes (boreal forest, temperate forest, tropical forest, grassland, shrubland, savanna, cropland), seven N species (NH3, NO, NO2, HNO2, HNO3, NH4+, NO3-), and regional aggregation across six continents.</p>