Investigating the Effects of Subgrid Cell Dynamic Heterogeneity on the Large-Scale Modeling of Albedo in Boreal Forests*

Stand-clearing disturbances, which remove most of the tree cover but are followed by forest regrowth, affect extensive areas annually, yet each event is usually much smaller than a typical grid cell in Earth system climate models. This study argues that the approach taken to account for the resulting subgrid cell dynamic heterogeneity substantially affects the computation of land–atmosphere exchanges. The authors investigated in a simplified model the effects of three such approaches on the computation of albedo over boreal forests. It was found that the simplest approach—in which any new disturbance-created patch was immediately merged with the rest of the grid cell—underestimated the annual reflected solar radiation by ~3 W m−2 on average (a relative error of 15%) compared with the most accurate approach—in which albedo computations were performed for each individual subgrid patch. This study also investigated an intermediate approach, in which each patch was tracked individually, but albedo was estimated from a much smaller number of subgrid tiles grouping patches having a similar amount of tree cover. Results from this third approach converged quickly toward the most accurate results as the number of tiles increased and were robust to changes in the thresholds used to assign patches to specific tiles. When computing time prevents implementing the most accurate approach in Earth system climate models, the results advocate for using strategies similar to the intermediate approach in order to avoid biasing the net radiative forcing of stand-clearing disturbances toward a warming impact, at least over boreal forests.