Harmonised biomass map from Global Forest Watch

Code

library(terra)
library(sf)

The forest carbon map from Harris et al. (2021) is a global, high-resolution (≈30 m) spatial dataset showing annual forest carbon fluxes—that is, where forests are absorbing (sinks) or releasing (sources) carbon by balancing removals into biomass against emissions from deforestation and disturbance from 2001–2019. Here we only used the initial carbon map estimated for the year 2000.

Code

if (!dir.exists("data/biomass")) dir.create("data/biomass")
paste0(
  "https://data-api.globalforestwatch.org/dataset/",
  "whrc_aboveground_woody_biomass_stock_2000/v1.4/download/",
  "geotiff?grid=10/40000&tile_id=20N_120E&pixel_meaning=Mg_ha-1&",
  "x-api-key=2d60cd88-8348-4c0f-a6d5-bd9adb585a8c"
) |>
  download.file("data/biomass/harris_carbon_2000.tif", "curl")

Code

area <- list.files("data/umrbpl/", full.names = TRUE, pattern = "shp$") |>
  read_sf() |>
  st_transform(crs = "EPSG:4326")
"data/biomass/harris_carbon_2000.tif" |>
  rast() |>
  crop(area) |>
  plot()

Harris, Nancy L., David A. Gibbs, Alessandro Baccini, Richard A. Birdsey, Sytze de Bruin, Mary Farina, Lola Fatoyinbo, et al. 2021. “Global maps of twenty-first century forest carbon fluxes.” Nature Climate Change 11 (3): 234–40. https://doi.org/10.1038/s41558-020-00976-6.