Random Forest Development and Modeling of Gross Primary Productivity in the Hudson Bay Lowlands

Jason Beaver,Elyn R. Humphreys,Douglas King

DOI: https://doi.org/10.1080/07038992.2024.2355937

IF: 2.6

2024-06-02

Canadian Journal of Remote Sensing

Abstract:Peatlands are a critical component of the global carbon cycle. Within Canada, the Hudson Bay Lowlands (HBL) has accumulated an estimated 33 Gt of carbon as peat because of a small but persistent difference between gross primary productivity (GPP) and ecosystem respiration over millennia. However, the impacts of disturbance and climate change on GPP are difficult to monitor across the HBL due to its vast size and remote location. This study evaluates the potential for random forest regression models to estimate GPP at five HBL eddy covariance flux monitoring sites using only optical data from MODIS (500 m, 8 day) or harmonized Landsat/Sentinel (HLS; 30 m, 16 day or more frequent). The results show that spatial resolution has less impact on modeled daily GPP compared to temporal resolution across model configurations. Using MODIS data, individual sites' daily GPP could be simulated with minimal bias, R 2 up to 0.89 and mean absolute error as low as 0.37 g C m −2 day −1 . For annual GPP, MODIS (R 2 = 0.84; mean absolute error 40.5 g C m −2 year −1 ) also outperformed the HLS models (R 2 = 0.46; mean absolute error 86.4 g C m −2 year −1 ). Les tourbières sont une composante essentielle du cycle du carbone. Au Canada, les basses terres de la baie d'Hudson (HBL) ont accumulé environ 33 Gt de carbone sous forme de tourbe en raison d'une différence légère, mais persistante entre la productivité primaire brute (GPP) et la respiration des écosystèmes au fil des millénaires. Cependant, les impacts des perturbations et du changement climatique sur la GPP sont difficiles à quantifier sur l'ensemble des HBL en raison de la vaste étendue du territoire et de son éloignement. Cette étude évalue le potentiel des modèles de régression de forêt aléatoire pour estimer la GPP sur cinq sites dans les HBL équipés de tour à covariance des turbulences, en utilisant uniquement les données optiques de MODIS (500 m, 8 jours) ou les données Landsat/Sentinel-2 harmonisées (HLS; 30 m, 16 jours ou plus fréquents). Les résultats montrent que dans toutes les configurations de modèle, la résolution spatiale a moins d'impact sur la GPP journalière modélisé que la résolution temporelle. À l'aide des données MODIS, la GPP journalière de chaque site a pu être simulé avec un faible biais, un R 2 jusqu'à 0.89 et une erreur absolue moyenne aussi faible que 0.37 g C m −2 jour −1 . Pour la GPP annuelle, MODIS (R 2 = 0.84; erreur absolue moyenne 40.5 g C m −2 an −1 ) a également surperformé les modèles issus des données HLS (R 2 = 0.46; erreur absolue moyenne 86.4 g C m −2 an −1 ).

remote sensing