Protected planet (protected areas), forests and intact forest landscapes at 100 m, 250 m to 1 km resolution

Protected planet (protected areas; version Oct 2024) and intact forest landscapes (2000, 2013, 2016 and 2020) rasterized to 100 m, 250 m and 1 km resolutions. The aggregated map contains all pixels that are either protected or intacts. To use these resources please refer to original data producers:

• Defourny, P., Lamarche, C., Bontemps, S., De Maet, T., Van Bogaert, E., Moreau, I., Brockmann, C., Boettcher, M., Kirches, G., Wevers, J., Santoro, M., Ramoino, F., & Arino, O. (2017). Land Cover Climate Change Initiative - Product User Guide v2. Issue 2.0. http://maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdf
• Olsson, E., Albrecht, R., & Golden Kroner, R.E. (2021). PADDDtracker Data Release Version 2.1: Technical Notes. Conservation International, Arlington, VA. DOI: 10.5281/zenodo.4749615.
• Potapov, P., Hansen, M. C., Laestadius L., Turubanova S., Yaroshenko A., Thies C., Smith W., Zhuravleva I., Komarova A., Minnemeyer S., Esipova E. The last frontiers of wilderness: Tracking loss of intact forest landscapes from 2000 to 2013. Science Advances, 2017; 3:e1600821
• UNEP-WCMC and IUCN (2024), Protected Planet: The World Database on Protected Areas (WDPA) [Online], October 2024, Cambridge, UK: UNEP-WCMC and IUCN. Available at: www.protectedplanet.net.

The time-series of forest areas (forest.areas_esa.cci_p) are based on the ESA CCI Land Cover time-series (2000–2022) 300-m resolution data; also available at 1-km resolution based on "average" resampling. Two maps (forest.cover.sum_esa.cci_p_250m and forest.cover.diff_esa.cci_p_250m) show long term cumulative forest cover and difference in forest cover for 2022 vs 2000.

The protected planet areas and intact forest landscapes were rasterized using:

## https://www.protectedplanet.net/en/thematic-areas/wdpa?tab=WDPA
for(j in 0:2){
  system(paste0('gdal_rasterize -ot Byte -a_nodata 0 -burn 100 -where "IUCN_CAT LIKE \'I%\'" /data/CCI_LandCover/WDPA_Oct2024_Public_shp_', j, '/WDPA_Oct2024_Public_shp-polygons.shp WDPA_Oct2024_Public_shp_', j, '_1km.tif -tr 0.008333333 0.008333333 -te -180 -65.00208 180 87.37 -co COMPRESS=DEFLATE -a_srs EPSG:4326'))
}
s = sds(rast("WDPA_Oct2024_Public_shp_ALL_0_1km.tif"), rast("WDPA_Oct2024_Public_shp_ALL_1_1km.tif"), rast("WDPA_Oct2024_Public_shp_ALL_2_1km.tif"))
dg.x = app(s, fun=max, na.rm=TRUE, cores = 32)
dg.x0 = terra::ifel(is.na(dg.x), 0, dg.x,
           filename="protected.areas_wdpa.all_p_1km_s_2023_2024_go_epsg4326_v20241025.tif",
           wopt=list(gdal=c("COMPRESS=DEFLATE"), datatype='INT2S'), overwrite=TRUE)
## https://intactforests.org/data.ifl.html
for(j in c(2000,2013,2016,2020)){
  system(paste0('gdal_rasterize -ot Byte -a_nodata 0 -burn 100 -l \"ifl_', j, '\" /mnt/lacus/raw/protectedplanet/ifl_', j, '.shp intact.forest_gfw_p_1km_s_', j, '0101_', j, '1231_go_epsg4326_v20241025.tif -tr 0.008333333 0.008333333 -te -180 -65.00208 180 87.37 -co COMPRESS=DEFLATE -a_srs EPSG:4326'))
}
## Combination IFL & WPDA
b = sds(rast("protected.areas_wdpa.all_p_1km_s_2023_2024_go_epsg4326_v20241025.tif"), rast("intact.forest_gfw_p_1km_s_20200101_20201231_go_epsg4326_v20241025.tif"))
bg.x = app(b, fun=max, na.rm=TRUE, cores = 32)
bg.x0 = terra::ifel(is.na(bg.x), 0, bg.x,
                    filename="protected.intact.areas_wdpa.ifl_p_1km_s_2020_2024_go_epsg4326_v20241025.tif",
                    wopt=list(gdal=c("COMPRESS=DEFLATE"), datatype='INT2S'), overwrite=TRUE)