Estimating Forest Biomass from Satellite Data: Carbon Stock Mapping from Space
Quick Answer: Forest biomass — the dry weight of living trees — is the largest terrestrial carbon pool after soils. Satellites estimate biomass through three complementary approaches: (1) SAR backscatter intensity correlates with biomass up to saturation limits (C-band: ~50 t/ha, L-band: ~150 t/ha, P-band: ~300+ t/ha), (2) LiDAR-derived canopy height strongly predicts biomass (taller trees = more biomass), (3) Optical indices and texture capture canopy structure. Global biomass maps at 100m resolution now exist (ESA CCI Biomass, NASA GEDI+). Accuracy is typically ±20-40% at pixel level, improving to ±10-15% at landscape scale. The ESA BIOMASS mission (P-band SAR, launching 2025) will be the first satellite specifically designed for forest biomass estimation.
Forests contain roughly 400 billion tonnes of carbon in their living biomass — trunks, branches, roots, and leaves. When forests are cleared or degraded, this carbon enters the atmosphere as CO₂. When forests grow, they absorb CO₂ and store it as biomass. Quantifying this carbon pool — knowing how much biomass is where — is fundamental to climate science, carbon markets, and forest management.
The problem is that you can't weigh a forest. Traditional biomass measurement requires destructive sampling: felling trees, drying them, and weighing the pieces. This has been done for thousands of individual trees worldwide, but extrapolating point measurements to entire continents requires satellites.
How Biomass Relates to What Satellites See
SAR Backscatter and Biomass
Radar signals interact with forest structure in ways that correlate with biomass:
Volume scattering: Radar waves penetrate the forest canopy and scatter off branches, trunks, and leaves. Denser forests with more woody material produce more backscatter — up to a saturation point where adding more biomass doesn't increase backscatter further.
Saturation levels by wavelength:
- C-band (Sentinel-1, 5.6 cm): Saturates at ~50 t/ha. Useful for young forests and plantations but insufficient for mature tropical forests (200-400+ t/ha).
- L-band (ALOS-2, 23.6 cm): Saturates at ~100-150 t/ha. Penetrates deeper into the canopy; sensitive to larger branches and trunks.
- P-band (BIOMASS mission, 69 cm): Saturates at ~300+ t/ha. Penetrates to the forest floor; interacts with trunks and large branches. The only SAR frequency sensitive to the full biomass range of tropical forests.
HV polarization is most sensitive to biomass because volume scattering from randomly oriented branches depolarizes the radar signal, producing strong cross-polarized (HV) returns.
LiDAR Canopy Height
LiDAR measures canopy height directly by timing laser pulses. Forest height is the single strongest predictor of biomass — tall trees have thick trunks and large crowns containing most of the biomass.
GEDI (Global Ecosystem Dynamics Investigation): NASA's spaceborne LiDAR on the ISS measures forest canopy height at 25m footprint diameter, sampling billions of points across the tropics and temperate zones (±51.6° latitude).
ICESat-2: Photon-counting LiDAR providing canopy height profiles along orbital tracks globally (including polar regions that GEDI doesn't cover).
The allometric relationship between height and biomass varies by forest type:
- Tropical rainforest: 50m trees may contain 300-500 t/ha
- Boreal forest: 20m trees may contain 80-150 t/ha
- Plantation: 30m trees in uniform stands may contain 150-250 t/ha
Optical Indicators
Optical satellite data contributes to biomass estimation through:
Canopy cover: NDVI and related indices estimate the fraction of ground covered by tree canopy. Dense canopy cover indicates more trees (but not necessarily taller or thicker trees).
Texture: Image texture measures (variance, entropy, correlation of pixel values within a window) capture canopy structure — rough texture indicates large crowns and gaps; smooth texture indicates uniform canopy.
Phenology: Seasonal NDVI patterns distinguish forest types with different biomass densities.
Estimation Methods
Empirical Models
Statistical relationships between satellite features and field-measured biomass:
- Establish field plots with measured biomass (from allometric equations applied to tree diameter and height measurements)
- Extract satellite features (SAR backscatter, canopy height, optical indices) at plot locations
- Fit regression models (linear, power-law, Random Forest, gradient boosting)
- Apply the model wall-to-wall across the satellite coverage
Physical Models
Radar backscatter models (like the Water Cloud Model) simulate how radar interacts with forest structure. By inverting these models, biomass can be estimated from observed backscatter.
Multi-Sensor Fusion
The most accurate biomass estimates combine multiple data sources:
- LiDAR height (strongest individual predictor, but spatially sparse from GEDI)
- SAR backscatter (continuous spatial coverage, but saturates)
- Optical data (canopy structure, forest type classification)
The approach: use GEDI LiDAR as training data to build wall-to-wall biomass maps from SAR and optical data that have complete spatial coverage.
Global Biomass Maps
ESA CCI Biomass
- Resolution: 100m
- Coverage: Global
- Epochs: 2010, 2017, 2018, 2019, 2020
- Method: Integration of SAR (ALOS PALSAR, Sentinel-1, Envisat ASAR), LiDAR (GEDI, ICESat-2), and optical data
- Product: Above-ground biomass (AGB) in tonnes of dry matter per hectare
NASA GEDI L4B
- Resolution: 1km
- Coverage: ±51.6° latitude (GEDI coverage)
- Method: GEDI-derived height + allometric models
- Strength: Direct height measurement foundation
GlobBiomass
- Resolution: 100m
- Coverage: Global
- Epoch: 2010
- Method: SAR-based (ALOS PALSAR + Envisat)
The ESA BIOMASS Mission
Launching in 2025, BIOMASS will be the first satellite carrying a P-band SAR:
- Wavelength: 69 cm — penetrates entire forest canopy to interact with trunks
- Polarimetry: Full quad-pol (HH, HV, VH, VV) providing complete polarimetric information
- Tomography: Multi-baseline InSAR enabling 3D imaging of forest vertical structure
- Coverage: Global forest areas
BIOMASS addresses the fundamental C-band and L-band saturation limitation. For tropical forests where most forest carbon resides (and where C-band saturates at ~50 t/ha and L-band at ~150 t/ha), P-band sensitivity extends to 300+ t/ha — covering the full biomass range.
Applications
National Greenhouse Gas Inventories
Countries report forest carbon stocks and changes under the UN Framework Convention on Climate Change (UNFCCC). Satellite biomass maps provide:
- Spatially explicit carbon stock estimates
- Change detection for emissions/removals
- Independent verification of reported values
REDD+ (Reducing Emissions from Deforestation and Degradation)
Carbon credit programs that pay countries or communities for avoiding deforestation need:
- Baseline biomass maps (how much carbon is at risk)
- Monitoring of biomass change (did the carbon stay or leave?)
- Measurement, Reporting, and Verification (MRV) systems
Satellite biomass maps are the primary data source for REDD+ MRV at national and project scales.
Forest Management
Commercial forestry uses biomass estimates for:
- Timber volume estimation (biomass relates to merchantable wood volume)
- Growth rate monitoring
- Harvest planning
- Certification compliance (sustainable forest management standards)
Accuracy
| Scale | Typical RMSE |
|---|---|
| Pixel (100m) | ±30-50% of mean biomass |
| Plot (1 ha) | ±20-40% |
| Landscape (10,000 ha) | ±10-20% |
| National | ±10-15% |
Errors decrease with aggregation because positive and negative pixel-level errors partially cancel. For carbon accounting at national scale (where UNFCCC reporting occurs), the accuracy is generally acceptable. For individual project-level REDD+ verification, pixel-level errors remain a concern.
The quest for accurate forest biomass estimation drives some of the most sophisticated satellite remote sensing science. The combination of radar backscatter, LiDAR height, and optical structure information — calibrated against extensive field measurement networks — is producing biomass maps of unprecedented accuracy and coverage. With the upcoming BIOMASS mission adding P-band capability, the remaining saturation limitation for tropical forests will be largely addressed, completing the global carbon inventory that climate science requires.
