TIMBER v0.1: a conceptual framework for emulating temperature responses to tree cover change

Land cover changes have been proposed to play a significant role, alongside emission reductions, in achieving the temperature goals agreed upon under the Paris Agreement. Such changes carry both global implications, pertaining to the biogeochemical effects of land cover change and thus the global carbon budget, and regional or local implications, pertaining to the biogeophysical effects arising within the immediate area of land cover change. Biogeophysical effects of land cover change are of high relevance to national policy and decision makers, and accounting for them is essential for effective deployment of land cover practices that optimise between global and regional impacts. To this end, Earth system model (ESM) outputs that isolate the biogeophysical responses of climate to land cover changes are key in informing impact assessments and supporting scenario development exercises. However, generating multiple such ESM outputs in a manner that allows comprehensive exploration of all plausible land cover scenarios is computationally untenable. This study proposes a framework to explore in an agile manner the local biogeophysical responses of climate under customised tree cover change scenarios by means of a computationally inexpensive emulator, the Tree cover change clIMate Biophysical responses EmulatoR (TIMBER) v0.1. The emulator is novel in that it solely represents the biogeophysical responses of climate to tree cover changes, and it can be used as either a standalone device or as a supplement to existing climate model emulators that represent the climate responses from greenhouse gas (GHG) or global mean temperature (GMT) forcings. We start off by modelling local minimum, mean, and maximum surface temperature responses to tree cover changes by means of a month- and Earth system model (ESM)-specific generalised additive model (GAM) trained over the whole globe; 2 m air temperature responses are then diagnosed from the modelled minimum and maximum surface temperature responses using observationally derived relationships. Such a two-step procedure accounts for the different physical representations of surface temperature responses to tree cover changes under different ESMs whilst respecting a definition of 2 m air temperature that is more consistent across ESMs and with observational datasets. In exploring new tree cover change scenarios, we employ a parametric bootstrap sampling method to generate multiple possible temperature responses, such that the