Nationwide native forest structure maps for Argentina based on forest inventory data, SAR Sentinel-1 and vegetation metrics from Sentinel-2 imagery

SILVEIRA, EDUARDA M.O.; RADELOFF, VOLKER C.; MARTINUZZI, SEBASTIÁN; MARTINEZ PASTUR, GUILLERMO J.; BONO, JULIETA; POLITI, NATALIA; LIZARRAGA, LEONIDAS; RIVERA, LUIS O.; CIUFFOLI, LUCIA; ROSAS, YAMINA M.; OLAH, ASHLEY M.; GAVIER-PIZARRO, GREGORIO I.; PIDGEON, ANNA M.

REMOTE SENSING OF ENVIRONMENT

ELSEVIER SCIENCE INC

Año: 2023 vol. 285

0034-4257

Detailed maps of forest structure attributes are crucial for sustainable forest management, conservation, and forest ecosystem science at the landscape level. Mapping the structure of broad heterogeneous forests is challenging, but the integration of extensive field inventory plots with wall-to-wall metrics derived from synthetic aperture radar (SAR) and optical remote sensing offers a potential solution. Our goal was to map forest structure attributes (diameter at breast height, basal area, mean height, dominant height, wood volume and canopy cover) at 30-m resolution across the diverse 463,000 km2 of native forests of Argentina based on SAR Sentinel-1, vegetation metrics from Sentinel-2 and geographic coordinates. We modelled the forest structure attributes based on the latest national forest inventory, generated uncertainty maps, quantified the contribution of the predictors, and compared our height predictions with those from GEDI (Global Ecosystem Dynamics Investigation) and GFCH (Global Forest Canopy Height). We analyzed 3788 forest inventory plots (1000 m2 each) from Argentina´s Second Native Forest Inventory (2015–2020) to develop predictive random forest regression models. From Sentinel-1, we included both VV (vertical transmitted and received) and VH (vertical transmitted and horizontal received) polarizations and calculated 1st and 2nd order textures within 3 × 3 pixels to match the size of the inventory plots. For Sentinel-2, we derived EVI (enhanced vegetation index), calculated DHIs (dynamic habitat indices (annual cumulative, minimum and variation) and the EVI median, then generated 1st and 2nd order textures within 3 × 3 pixels of these variables. Our models including metrics from Sentinel-1 and 2, plus latitude and longitude predicted forest structure attributes well with root mean square errors (RMSE) ranging from 23.8% to 70.3%. Mean and dominant height models had notably good performance presenting relatively low RMSE (24.5% and 23.8%, respectively). Metrics from VH polarization and longitude were overall the most important predictors, but optimal predictors differed among the different forest structure attributes. Height predictions (r = 0.89 and 0.85) outperformed those from GEDI (r = 0.81) and the GFCH (r = 0.66), suggesting that SAR Sentinel-1, DHIs from Sentinel-2 plus geographic coordinates provide great opportunities to map multiple forest structure attributes for large areas. Based on our models, we generated spatially-explicit maps of multiple forest structure attributes as well as uncertainty maps at 30-m spatial resolution for all Argentina´s native forest areas in support of forest management and conservation planning across the country.