=== versions === python 3.12.0 | rasterio 1.5.2 | numpy 2.5.4 | pyproj 3.8.0 | pystac-client 0.9.0 | planetary-computer 1.0.0 === Global Forest Change layers === GFC-2025-v1.13: 4000 x 4000 pixels of 0.00025 deg over lon -64.0..-63.0, lat -11.0..-10.0 lossyear values: 0..25 (0 = no loss, 1..25 = 2001..2025) datamask: land 99.99%, water 0.01%, no data 0.00% === pixel area === at 10.5 S one pixel is 27.37 m wide x 27.65 m tall = 756.8 m2 (0.0757 ha) row areas range 755.6..757.9 m2 across the AOI pixel area at 0.0 deg: 769.3 m2 (0.85 x a 30 m square) pixel area at 10.5 deg: 756.8 m2 (0.84 x a 30 m square) pixel area at 30.0 deg: 668.5 m2 (0.74 x a 30 m square) pixel area at 50.0 deg: 498.4 m2 (0.55 x a 30 m square) pixel area at 60.0 deg: 388.5 m2 (0.43 x a 30 m square) pixel area at 70.0 deg: 266.3 m2 (0.30 x a 30 m square) AOI land area 1,210,727 ha; forest (>= 30% canopy) in 2000 1,007,753 ha (83.2% of land) pixel count x 0.09 ha would give 1,198,664 ha, +18.9% too much === forest loss by year (canopy >= 30% in 2000) === year loss_ha pct_of_2000_forest moving_avg_3yr 2001 27,360.5 2.7 NaN 2002 22,565.0 2.2 23,398.2 2003 20,269.1 2.0 22,426.4 2004 24,445.3 2.4 28,665.0 2005 41,280.8 4.1 28,187.9 2006 18,837.7 1.9 23,180.7 2007 9,423.7 0.9 11,388.7 2008 5,904.6 0.6 7,142.8 2009 6,100.2 0.6 7,063.9 2010 9,187.1 0.9 6,951.3 2011 5,566.7 0.6 8,134.4 2012 9,649.4 1.0 6,445.2 2013 4,119.6 0.4 8,148.4 2014 10,676.3 1.1 8,713.9 2015 11,345.8 1.1 12,084.4 2016 14,231.0 1.4 14,036.1 2017 16,531.4 1.6 12,873.1 2018 7,856.9 0.8 10,926.6 2019 8,391.5 0.8 7,936.8 2020 7,562.0 0.8 9,573.5 2021 12,767.1 1.3 11,023.6 2022 12,741.8 1.3 10,063.4 2023 4,681.3 0.5 9,967.1 2024 12,478.2 1.2 8,931.0 2025 9,633.4 1.0 NaN total 2001-2025: 333,606 ha = 33.1% of the 2000 forest; peak 2005 (41,281 ha) 2001-2010 mean 18,537 ha/yr; 2011-2020 mean 9,593 ha/yr; 2021-2025 mean 10,460 ha/yr canopy threshold 10%: forest 2000 1,023,246 ha, loss 2001-2025 335,496 ha canopy threshold 30%: forest 2000 1,007,753 ha, loss 2001-2025 333,606 ha canopy threshold 50%: forest 2000 987,944 ha, loss 2001-2025 329,503 ha canopy threshold 75%: forest 2000 942,816 ha, loss 2001-2025 313,231 ha loss pixels outside the 30% forest mask: 7,487 ha === zoom area === zoom: lon -63.520..-63.400, lat -10.720..-10.600 (480 x 480 GFC pixels) forest 2000 17,172 ha; loss 2001-2017 877 ha; loss 2018-2024 4,807 ha; loss 2025 884 ha === before/after imagery === Landsat 7 before: LE07_L2SP_232067_20000612_02_T1 cloud 0.0% Sentinel-2 before: S2B_MSIL2A_20170720T142749_R053_T20LMP_20210210T120110 cloud 0.0% baseline 02.12 Sentinel-2 after: S2C_MSIL2A_20250703T142811_R053_T20LMP_20250703T194403 cloud 0.0% baseline 05.11 === Sentinel-2 dry-season NBR, 2017-2025 === 2017: 3 scenes (20 Jul, 15 Jun, 25 Jul), valid 100.0% 2018: 3 scenes (25 Jun, 20 Jul, 10 Jul), valid 100.0% 2019: 3 scenes (15 Jun, 10 Jul, 05 Jun), valid 100.0% 2020: 3 scenes (08 Aug, 18 Aug, 29 Jul), valid 100.0% 2021: 3 scenes (19 Jul, 09 Jul, 04 Jul), valid 100.0% 2022: 3 scenes (09 Jul, 19 Jul, 19 Jul), valid 100.0% 2023: 3 scenes (04 Jul, 24 Jul, 23 Aug), valid 100.0% 2024: 3 scenes (18 Jul, 13 Jul, 23 Jul), valid 100.0% 2025: 3 scenes (03 Jul, 15 Jul, 27 Aug), valid 100.0% median NBR by year (rows: GFC class, columns: 2017..2025) stable forest 0.66 0.64 0.68 0.66 0.66 0.67 0.67 0.66 0.65 lost in 2020 0.62 0.62 0.66 0.50 0.50 0.46 0.35 0.22 0.31 lost in 2021 0.64 0.63 0.67 0.64 0.52 0.53 0.49 0.39 0.34 lost in 2022 0.65 0.64 0.67 0.65 0.64 0.54 0.54 0.42 0.37 lost in 2023 0.65 0.63 0.67 0.64 0.64 0.63 0.52 0.48 0.39 lost in 2024 0.62 0.61 0.65 0.61 0.60 0.59 0.60 0.56 0.41 lost in 2025 0.65 0.63 0.67 0.64 0.64 0.65 0.65 0.64 0.45 loss pixels 2018-2025 with a full NBR record: 75,238 largest annual NBR drop in the GFC year: 35.6%; one year later: 15.9%; one year earlier: 2.9%; further off: 45.6% largest drop >= 0.2 for 55.6% of GFC loss pixels vs 3.1% of stable-forest pixels dNBR 2017->2025 median: stable forest 0.005, lost 2018-2025 0.225, lost 2001-2017 0.109 === MODIS burned area vs GFC loss === year modis_months gfc_loss_ha loss_in_modis_burn_ha modis_burned_ha share_burned 2019 12 8,392 354 5,714 0.04 2020 12 7,562 610 4,648 0.08 2021 12 12,767 3,166 9,742 0.25 2022 12 12,742 1,307 5,565 0.10 2023 11 4,681 46 280 NaN 2024 10 12,478 1,996 19,309 NaN 2025 8 9,633 0.00 0.00 NaN complete years 2019-2022: 13.1% of GFC loss sat inside a same-year MODIS burn scar wrote figures 01-06, output/loss_by_year.csv and output/loss_vs_modis_burned.csv