Journal cover Journal topic
Earth System Dynamics An interactive open-access journal of the European Geosciences Union
Journal topic
Earth Syst. Dynam., 10, 473–484, 2019
https://doi.org/10.5194/esd-10-473-2019
Earth Syst. Dynam., 10, 473–484, 2019
https://doi.org/10.5194/esd-10-473-2019

Research article 19 Jul 2019

Research article | 19 Jul 2019

# Different response of surface temperature and air temperature to deforestation in climate models

Different response of surface temperature and air temperature to deforestation in climate models
Johannes Winckler1,2, Christian H. Reick1, Sebastiaan Luyssaert3, Alessandro Cescatti4, Paul C. Stoy5, Quentin Lejeune6,a, Thomas Raddatz1, Andreas Chlond1, Marvin Heidkamp1,2, and Julia Pongratz1,7 Johannes Winckler et al.
• 1Max Planck Institute for Meteorology, Hamburg, Germany
• 2International Max Planck Research School on Earth System Modeling, Hamburg, Germany
• 3Faculty of Science, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands
• 4Institute for Environment and Sustainability, Joint Research Centre, European Commission, Ispra, Italy
• 5Department of Land Resources and Environmental Sciences, Montana State University, Bozeman, MT, USA
• 6Institute for Atmospheric and Climate Science, ETH-Zürich, Zurich, Switzerland
• 7Ludwig-Maximilians-Universität München, Munich, Germany
• acurrently at: Climate Analytics, Berlin, Germany

Correspondence: Julia Pongratz (julia.pongratz@geographie.uni-muenchen.de)

Abstract

When quantifying temperature changes induced by deforestation (e.g., cooling in high latitudes, warming in low latitudes), satellite data, in situ observations, and climate models differ concerning the height at which the temperature is typically measured/simulated. In this study the effects of deforestation on surface temperature, near-surface air temperature, and lower atmospheric temperature are compared by analyzing the biogeophysical temperature effects of large-scale deforestation in the Max Planck Institute Earth System Model (MPI-ESM) separately for local effects (which are only apparent at the location of deforestation) and nonlocal effects (which are also apparent elsewhere). While the nonlocal effects (cooling in most regions) influence the temperature of the surface and lowest atmospheric layer equally, the local effects (warming in the tropics but a cooling in the higher latitudes) mainly affect the temperature of the surface. In agreement with observation-based studies, the local effects on surface and near-surface air temperature respond differently in the MPI-ESM, both concerning the magnitude of local temperature changes and the latitude at which the local deforestation effects turn from a cooling to a warming (at 45–55 N for surface temperature and around 35 N for near-surface air temperature). Subsequently, our single-model results are compared to model data from multiple climate models from the Climate Model Intercomparison Project (CMIP5). This inter-model comparison shows that in the northern midlatitudes, both concerning the summer warming and winter cooling, near-surface air temperature is affected by the local effects only about half as strongly as surface temperature. This study shows that the choice of temperature variable has a considerable effect on the observed and simulated temperature change. Studies about the biogeophysical effects of deforestation must carefully choose which temperature to consider.

1 Introduction

Afforestation has been proposed as a tool to mitigate climate change globally (UNFCCC2011), mainly because forests can store large amounts of carbon . In addition, changes in forest cover can cause a warming or cooling via an alteration of the exchange of energy and water between the Earth's surface and the atmosphere, i.e., the so-called biogeophysical effects (Bonan2008). Earth System models have been employed to assess how these biogeophysical effects affect the temperature of the surface and the temperature of the near-surface air (usually air temperature 2 m above zero-plane displacement height) . The different temperature variables that are considered in studies about deforestation effects are relevant for different questions and applications. Satellite-based studies on changes in radiometric surface temperature provide important information about the biophysical mechanisms of surface energy partitioning and thereby surface–atmosphere interactions . Compared to changes in surface temperature, changes in air temperature may be considered more relevant for human living conditions because of their importance, e.g., for the perceived temperature (e.g., Staiger et al.2011). Within- and below-canopy air temperature (which is not included in this study) is the most relevant variable for many organisms that live within forests (e.g., De Frenne et al.2013, 2019). The coupling between ground temperature and air temperature is strongly influenced by the type of vegetation that covers the surface , but it remains unclear whether surface temperature and near-surface air temperature respond differently to deforestation in climate models. This is the focus of the present study.

An answer to this question could help to reconcile apparent inconsistencies in observation-based studies on the effects of deforestation on surface temperature and air temperature. Studies based on satellite observations investigated changes in radiometric surface temperature resulting from deforestation . These studies reported that deforestation results in a local cooling in the boreal regions (north of approximately 45–55N) and a warming in lower latitudes, highlighting that the contributions from changes in surface albedo and other surface properties vary with latitude . Studies based on observations of air temperature from weather stations and FLUXNET towers also reported a deforestation-induced boreal local cooling and a warming for lower latitudes, but they indicated that the transition between cooling and warming is located further south (at approximately 35 N). It remains unclear whether part of this apparent inconsistency can be attributed to the different heights above the surface at which temperature changes are considered. In contrast to observations, climate models allow us to assess the biogeophysical effects on surface temperature, on near-surface air temperature, and on the temperature of the atmosphere within a single consistent framework, rendering climate models suitable tools to investigate this question.

Both air and surface temperature are affected by local and nonlocal biogeophysical effects of deforestation. We define local effects as effects that are only apparent in deforested locations and nonlocal effects as effects that are also apparent in non-deforested locations (see Sect. 2.1 and Winckler et al.2017). Local effects can, for example, be caused by a redistribution of heat between the surface and the atmosphere (e.g., Vanden Broucke et al.2015), while the nonlocal effects can be caused by advection or by changes in global circulation . In this study, local and nonlocal effects are analyzed separately for three reasons. First, the difference between local and nonlocal effects matters for decision makers: the local effects may be relevant for policies that aim at adapting to a warming climate locally because they link the climate effects to the areas where policies are implemented . The nonlocal effects are also relevant for international policies that aim at mitigating global climate change because the nonlocal effects may dominate the global mean biogeophysical temperature response to deforestation . Second, the observation-based data sets only record the local effects when comparing temperature between nearby locations with and without forest, or between locations with and without deforestation. The nearby locations share the same background climate, and thus the nonlocal effects cancel out when temperature differences between the locations are considered . For a consistent comparison to observation-based data sets, the local effects need to be separated from the nonlocal effects when analyzing climate model results. The third reason to consider local and nonlocal temperature changes separately is that different mechanisms trigger local and nonlocal temperature changes . If surface and air temperature respond differently to deforestation, it is unclear whether this difference arises from mechanisms that trigger the local temperature changes (predominantly via changes in the turbulent heat fluxes; ), from mechanisms that trigger the nonlocal temperature changes (predominantly via the incoming radiation that reaches the surface; ), or from both. A separate analysis of local and nonlocal temperature changes facilitates an investigation of the mechanisms that may cause a different response of surface and air temperature to deforestation.

Here, we investigate how deforestation in the MPI-ESM (Max Planck Institute Earth System Model) affects surface and air temperature differently and analyze this separately for the local and nonlocal effects. A previous study contrasted the response of surface temperature only with the response of near-surface air temperature and found mainly differences between surface and air temperature for the local effects (Appendix C in Winckler et al.2017). We go beyond this previous study by additionally analyzing the effects on temperature in the lowest atmospheric layer and by using simulations with an interactive ocean which enables us to better capture the nonlocal temperature effects of deforestation from the surface to the lower atmosphere . To further analyze the mechanisms that are responsible for differences in these three temperature variables, we investigate the local effects separately for the response in mean daily minimum and maximum temperature. To test the robustness of our results for this particular climate model, we compare the simulation results of the MPI-ESM to existing simulation results from multiple climate models from the Climate Model Intercomparison Project CMIP5. In this inter-model comparison, we contrast the response of the local effects on near-surface air temperature and surface temperature.

2 Methods

## 2.1 Simulations of large-scale deforestation in the MPI-ESM

Using the fully coupled climate model MPI-ESM , the temperature response to deforestation at the surface, at 2 m and the lowest layer of the atmosphere are obtained from simulations of large-scale deforestation. Simulations are performed at T63 atmospheric resolution (about 1.9) for 550 years, and the last 200 years (which are free of substantial trends in the investigated variables, not shown) are used for the analysis. Two simulations are performed: a first simulation (“forest world”) with forest plant functional types on all areas that can potentially be covered with vegetation (i.e., forests do not exist in deserts, etc., Fig. S1 in the Supplement). These vegetated areas are taken from a previous study , which derived a map of potential vegetation from remote sensing . The non-forest plant functional types were replaced by the forest types occurring in that grid cell (preserving the relative fraction of the different forest types). In a second simulation forests in the forest world are completely replaced by grasslands in three of four grid boxes in a regular spatial pattern (Fig. S1, equivalent to simulation “3∕4” in a previous study; ). In both simulations, atmospheric CO2 concentrations are prescribed at preindustrial level in order to obtain only the biogeophysical effects of deforestation. The total, i.e., local plus nonlocal, biogeophysical deforestation effects are then computed as the differences (e.g., in temperature) between these two simulations.

Following the approach in , the total effects can be separated into the local and nonlocal effects of deforestation as follows:

$\begin{array}{}\text{(1)}& \mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{total}}=\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{local}}^{\prime }+\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{nonlocal}},\end{array}$

where Δ Ttotal are the temperature changes that are simulated at a deforested grid box and $\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{local}}^{\prime }=\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{local}}+\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{local}×\mathrm{nonlocal}}$ includes both the local effects and possible interactions between local and nonlocal effects. However, such interactions were found to be small across a wide range of deforestation scenarios , and in the following we refer $\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{local}}^{\prime }$ as “local effects”. The nonlocal effects are determined from non-deforested grid boxes, where only the nonlocal effects are present. The nonlocal effects are spatially interpolated to the deforested grid boxes using bilinear interpolation. The local effects at deforested grid boxes can thus be obtained by subtracting the nonlocal effects from the simulated total effects:

$\begin{array}{}\text{(2)}& \mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{local}}^{\prime }=\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{total}}-\mathrm{\Delta }\phantom{\rule{0.125em}{0ex}}{T}_{\mathrm{nonlocal}}.\end{array}$

The local effects that are obtained this way are similar to the local effects that were obtained in previous studies by comparing temperatures in sub-grid tiles . In contrast to their methods, the method that is applied in this study includes local feedbacks between the surface and the atmosphere leading, for example, to local changes in clouds or precipitation. In addition, the method used in this study allows us to assess the local effects on the temperature of the lowest layer of the atmosphere (which in most models is not calculated separately for sub-grid tiles) and to assess nonlocal effects on temperature. The choice of deforesting three of four grid boxes is to some extent arbitrary; varying the spatial extent of deforestation influences the magnitude of the nonlocal effects on surface temperature , but the local effects on surface temperature within a grid box are largely insensitive to deforestation elsewhere . A detailed description and discussion of the separation approach can be found in .

## 2.2 Temperature of the surface, the lowest atmospheric layer, and near-surface air in the MPI-ESM

This study investigates the response of three types of temperature to deforestation in the MPI-ESM: surface temperature (Tsurf), the temperature of the lowest atmospheric layer (Tatm), and near-surface air temperature (T2 m, called “tas” in CMIP5). Although these temperature variables are part of the standard output of climate model simulations, different models may calculate them differently. In the following, details are provided on the calculations of these variables in the MPI-ESM.

The surface temperature Tsurf in the MPI-ESM is determined by solving the surface energy balance equation in a bulk canopy layer. For simplicity this layer has a heat capacity that is independent of the vegetation type. This bulk canopy layer exchanges heat with deeper soil layers via the ground heat flux.

It is not possible to assign one geometrical height to the surface layer in the MPI-ESM because there is an internal inconsistency between the two different aspects that are involved in the process of solving the surface energy balance equation: the calculation of the surface radiative budget (absorption of solar radiation and emission of terrestrial radiation to the atmosphere) and the calculation of the turbulent heat fluxes (latent and sensible heat). From the perspective of the radiative budget, the surface is where this radiative budget is calculated (i.e., where the energy balance is solved). In the presence of vegetation this is somewhere in the canopy, but geometrically its exact height cannot be specified. From the perspective of turbulent fluxes, the geometrical height d+z0 above the surface is where the wind speed would become zero in the wind profile based on Monin–Obukhov theory . Here, z0 denotes the aerodynamic roughness length and d is the zero-plane displacement height. This d takes into account the displacement effect exerted by vegetation . Geometrically the height d+z0 may differ from the height where the radiative budget is calculated.

What does this inconsistency imply for the comparison between Tsurf in the MPI-ESM and satellite-based products? For comparison with satellite observations only the radiative perspective is relevant – because satellites estimate temperature based on the emissions of terrestrial radiation. That the Monin–Obukhov theory provides a different definition of surface height must be considered as a special approximation to solve the energy balance but has no consequences for comparison with satellite observations of Tsurf.

The “atmospheric temperature” Tatm is defined here as the temperature of the lowest of the 47 atmospheric layers in the MPI-ESM . The thickness of this layer is around 60 m (at 15 C), and the temperature is volume-averaged in this layer. This temperature is used for the calculation of the turbulent heat fluxes and Tsurf.

The near-surface air temperature T2 m is estimated in the MPI-ESM as temperature of the air at 2 m above d+z0. Because it is unclear (and irrelevant for the calculations) where within the canopy this d+z0 is, a comparison of T2 m between the MPI-ESM and observations is challenging, especially in forests (see Sect. 4). The MPI-ESM does not have a representation of within-canopy air temperature or separate temperatures of the surface and the vegetation canopy.

In the MPI-ESM, following , T2 m is obtained via a procedure based on Monin–Obukhov similarity theory that uses the values at the surface and the lowest atmospheric layer. This procedure employs dry static energy instead of temperature because dry static energy is a conserved quantity in an adiabatic process.

$\begin{array}{}\text{(3)}& {s}_{{\mathrm{z}}_{\mathrm{aero}}}={c}_{p}\phantom{\rule{0.125em}{0ex}}{T}_{{\mathrm{z}}_{\mathrm{aero}}}+g\phantom{\rule{0.125em}{0ex}}{z}_{\mathrm{aero}},\end{array}$

where cp is the heat capacity of moist air, and g is the gravitational acceleration of the Earth, and ${s}_{{z}_{\mathrm{aero}}}$ and ${T}_{{z}_{\mathrm{aero}}}$ are the dry static energy and temperature at the aerodynamic height ${z}_{\mathrm{aero}}=z-\left(d+{z}_{\mathrm{0}}\right)$ where z is the height above the surface.

At 2 m above d+z0, the dry static energy is then obtained as follows:

$\begin{array}{}\text{(4)}& {s}_{\mathrm{2}\phantom{\rule{0.125em}{0ex}}\mathrm{m}}={s}_{\mathrm{surf}}+\left({s}_{\mathrm{atm}}-{s}_{\mathrm{surf}}\right)\phantom{\rule{0.125em}{0ex}}\mathit{\gamma }\left(\frac{\mathrm{2}\phantom{\rule{0.125em}{0ex}}\mathrm{m}}{{z}_{\mathrm{0}}},{R}_{i}\right),\end{array}$

where ssurf and satm denote the dry static energy at the surface and the lowest atmospheric layer, and γ is a nonlinear function based on Monin–Obukhov similarity theory with values ranging between 0 and 1 that depends on the roughness length z0 and on the bulk Richardson number Ri, which is closely related to the temperature gradient TsurfTatm. Different functions for γ are used for near-surface neutral (Ri≈0), stable (Ri<0), and unstable conditions (Ri>0) (Geleyn1988; ECMWF Research Department1991, Sect. 3.1.3). Note that both ssurf and satm, but also Ri and z0 are affected by deforestation. After this procedure, T2 m is obtained as Taero from Eq. (3) by setting zaero=2 m.

## 2.3 Isolation of local effects across CMIP5 models

The response of Tsurf and T2 m to deforestation in the Northern Hemisphere midlatitudes is compared across a wide range of climate models from CMIP5: CanESM2, CCSM4, CESM1-CAMS, GFDL-CM3, HadGEM2-ES, IPSL-CM5A-LR, MPI-ESM-LR, and NorESM1-M. Given that the these models did not simulate the 3/4 deforestation (see Sect. 2.1), we invoke the difference between “historical” and “piControl” simulations to isolate the local temperature response from the CMIP5 ensemble . The historical simulations are subject to all forcings including changes in greenhouse gases and land use, while the piControl simulations are subject to constant boundary conditions and no forcings. To isolate the local effects of deforestation, we use a method that was already applied and validated on these simulations . This method assumes that temperature in neighboring grid boxes can be affected differently by the local effects of deforestation, depending on the forest cover change in each grid box, whereas other climate forcings (like greenhouse gases, but also the nonlocal effects) influence neighboring grid boxes in a similar way. Linear regressions are fitted between temporal changes in temperature (the so-called predictand) and forest cover change (the so-called predictor) within a spatially moving window encompassing 5×5 model grid boxes. In the center of this moving window, the local effects are then defined as the temperature change for a hypothetical conversion of 100 % forest into 100 % open land (given by the slope of the regression) and is by construction largely independent of the changes due to the nonlocal greenhouse gas forcing and nonlocal deforestation effects (given by the y intercept of the regression).

We consider here the difference between the last 30 years (1971–2000) of historical simulations for which data in all models are available and 30 years of the preindustrial control simulations (piControl), from which the temporal changes in both temperature variables and forest fraction since 1860 are computed. The historical simulations consist of several ensemble members for each model, where each ensemble member experiences the same forcings but starts from different initial conditions. The moving-window method is applied to several combinations of ensemble members from the historical simulations and time slices from the piControl simulations for each model, and the number of analyzed ensemble members is shown in Table S1 in the Supplement.

For this inter-model comparison, we focus on the local effects for two reasons. (1) The nonlocal deforestation effects cannot be isolated from the analyzed set of simulations because the nonlocal deforestation effects cannot be distinguished from nonlocal greenhouse gas forcing in these simulations. (2) The local effects exhibit a better signal  noise ratio compared to the nonlocal effects (e.g., Lejeune et al.2017). This is important because the climate variability can be large compared to the nonlocal effects for the short time spans (30 years) that are analyzed here . Furthermore, climate variability is especially large in the midlatitudes that are analyzed in the inter-model comparison.

3 Results

## 3.1 Different temperature response of surface and air temperature in the MPI-ESM

In the MPI-ESM, deforestation in three of four grid boxes triggers substantial nonlocal cooling in most regions (Fig. 1). This happens because deforestation locally reduces the input of latent and sensible heat from the surface to the atmosphere . Thus, the atmosphere becomes cooler and drier (not shown), and this leads to a reduction in Tsurf in many regions, mainly because of reduced longwave incoming radiation . The spatial pattern of these nonlocal effects is very similar for Tsurf, T2 m, and Tatm.

Figure 1The deforestation-induced annual mean response of temperature in the lowest atmospheric layer (Tatm), near-surface air (T2 m), and at the surface (Tsurf). Deforestation was applied to three of four grid boxes (Fig. S1). Stippling indicates where results are not statistically significant at a 5 % level for a Student's t test accounting for lag-1 autocorrelation . Zonal averages are shown in Fig. S2.

In contrast to the nonlocal effects, the local effects differ strongly between Tsurf, T2 m, and Tatm. Deforestation strongly influences the local surface energy balance: the imposed changes in surface properties in the model (surface albedo, evapotranspirative efficiency and surface roughness) cause a surface warming for the local effects in most regions, except for the high northern latitudes where the local effects cause a surface cooling (Fig. 1). The changes in surface properties influence not only the local Tsurf but also the flux of sensible heat from the surface into the lower boundary layer (not shown). Intuitively one would expect that the change in sensible heat flux alters Tatm; e.g., an increased input of sensible heat into the atmosphere could raise the temperature of the atmospheric air above a deforested location. However, in our model results Tatm is largely unaffected by the local effects of deforestation (Fig. 1). We interpret this lack of local effects in Tatm as follows: it takes long enough for the lowest atmospheric layer to warm up (due to the deforestation-induced increase in sensible heat flux) for the heated air to be transported to higher atmospheric layers and the neighboring grid boxes. Due to the advection, the change in Tatm is hence not only seen in a deforested location but also in nearby grid boxes that are not deforested. Thus, this warming or cooling is accounted for in the nonlocal effects. In the nearby grid boxes, the change in Tatm and/or atmospheric moisture can then also influence Tsurf via changes in longwave incoming radiation , which could explain why the nonlocal effects are similar for Tatm and Tsurf. While in Tatm, advection can lead to a direct exchange of heat between neighboring grid cells, the same is not possible for Tsurf; there is no direct horizontal exchange of heat between the surface of neighboring grid cells, and this difference (advection for Tatm but not Tsurf) may explain why local effects can be seen in Tsurf but not in Tatm.

Because T2 m is an interpolation between Tsurf and Tatm, we expected that the local response of T2 m would also lie in between the response of Tsurf and Tatm. In a lot of regions this is the case, but in other regions, most notably those that show a cooling, the local effects on T2 m seem to cool more than Tsurf, and in some regions even the sign differs between ΔTsurf and ΔT2 m (e.g., parts of the US and regions in the southern extratropics; Fig. S3). The different response of Tsurf and T2 m in relation to the observation-based findings is discussed in Sect. 4.

To better understand the apparent discrepancy between ΔTsurf and ΔT2 m, we separately analyze the local temperature response for boreal winter (DJF) and summer (JJA) and the response of mean daily minimum temperature (Tmin, which approximately corresponds to nighttime conditions) and maximum temperature (Tmax, which approximately corresponds to daytime conditions). For Tsurf, the response to deforestation locally differs strongly between DJF, JJA, Tmin, and Tmax values (Fig. 2). For Northern Hemisphere DJF and Tmin, deforestation leads to a local Tsurf cooling, while for JJA and Tmax, deforestation leads to a local Tsurf warming. This is qualitatively in good agreement with observation-based studies that show a local cooling in the boreal regions in DJF and in agreement with the local increase in the diurnal amplitude due to deforestation . Similarly as in the case of long-term mean temperature (Fig. 1), Tatm locally shows little response to deforestation, neither for DJF, JJA, Tmin, nor Tmax (Fig. 2). For both Tsurf and T2 m, the annual mean response then depends on the balance between the daytime and nighttime response and on the balance between the responses in different seasons.

Figure 2Seasonal and diurnal temperature response to the local effects of deforestation, separately for boreal winter (DJF) and summer (JJA) and daily maximum (Tmax) and minimum temperature (Tmin). Stippling indicates where results are not statistically significant at a 5 % level for a Student's t test accounting for lag-1 autocorrelation .

For near-surface air temperature, the Tmax response is substantially weaker, and in many areas of opposite sign, than for the surface, similar to the lowest atmospheric layer (land mean absolute changes for Tmax: 1.62 K for Tsurf, 0.19 K for T2 m, and 0.10 K for Tatm). On the contrary, most regions exhibit a strong Tmin cooling of T2 m, similar to the Tmin response of Tsurf (land mean absolute changes for Tmin: 0.67 K for Tsurf, 0.48 K for T2 m, and 0.10 K for Tatm).

Tmax responds differently for Tsurf and T2 m not only for annual mean Tmax, but in some tropical or subtropical regions also for Tmax during DJF and JJA (Fig. S9). We interpret this as follows: during the daytime when Tmax occurs, the Tsurf is higher than the temperature of the lowest atmospheric layer (Fig. S4) because the surface is heated by incoming radiation. In accordance with previous studies (e.g., Li et al.2015), deforestation further increases Tsurf (Fig. 2) and thus also the difference between Tsurf and Tatm (illustrated in Fig. S5b). Intuitively, one would expect that an increase in this difference would result also in an increase in T2 m. But by the increased surface temperature, the well-mixed zone in the boundary layer not only extends to larger heights, but also extends further down. Accordingly, the cooler air from above mixes further down, and if this affects heights below 2 m, T2 m is lowered. In the model's calculation of T2 m, a deforestation-induced increase in Tsurf increases the difference between Tatm and Tsurf and thus leads to an even more negative Richardson number Ri. This Richardson number enters the similarity function for the calculation of s2 m (γ in Eq. 4; see underlying report; ) such that the vertical profile of s2 m becomes more nonlinear and the calculated T2 m approaches Tatm. As a result, daily maximum T2 m in the model may decrease although Tsurf increases.

During nighttime when Tmin occurs, the surface loses energy via outgoing longwave radiation, and thus the surface is often cooler than the overlying atmosphere (Fig. S4). In accordance with previous studies (e.g., Schultz et al.2017), deforestation further decreases Tsurf and thus Tsurf and Tatm further diverge. For T2 m, one expects that less vertical mixing than during the daytime results in T2 m tending towards Tsurf. In the model's calculation of T2 m, this behavior is a result of the γ in Eq. (4) being less nonlinear in stable compared to unstable conditions (for the concrete form of the γ function, see ) with the consequence that T2 m at deforested locations stays closer to Tsurf compared to daytime conditions (Fig. S5a). As a result, in most regions the calculated T2 m follows the deforestation-induced nighttime surface cooling but not necessarily the deforestation-induced daytime surface warming.

The different T2 m responses for DJF and JJA in the model may be caused by the same mechanism that may be responsible for the different behaviors for Tmin and Tmax. In JJA, when the surface is often warmer than the overlying atmosphere, deforestation leads to a further local Tsurf warming (Fig. 2). Similar to the case of Tmax, T2 m barely responds or even responds with a cooling in some regions. In DJF, when the surface is often cooler than the overlying atmosphere, deforestation leads to a local Tsurf cooling. Similar to the case of Tmin, T2 m also shows a substantial cooling in most Northern Hemisphere regions. The reason why T2 m responds similarly to Tsurf in Northern Hemisphere DJF but not JJA may be similar to the reason why T2 m responds similar to Tsurf for Tmin but not for Tmax, as described in the previous two paragraphs. Our findings are in qualitative agreement with the findings of (their Fig. 9), who found a strong deforestation-induced daytime Tsurf warming and a moderate T2 m cooling in many regions using the CLM (Community Land Model). Whether or not Tsurf and T2 m respond similarly to deforestation may strongly depend on how T2 m is calculated in the respective climate model. For the summer and winter response in the northern midlatitudes, the responses of Tsurf and T2 m are compared across climate models from CMIP5 in the following.

## 3.2 Different temperature response of surface temperature and air temperature across climate models

While the above results refer only to the MPI-ESM climate model, quantitatively the local effects on Tsurf and T2 m also differ in other climate models (Fig. 3). We analyze the local effects of historical deforestation and the average over midlatitude areas (40–60 N) that have experienced intense deforestation (≥15 %) since 1860. We choose the northern midlatitudes for two reasons. (1) This is where most historical deforestation happened, and regions with intense deforestation are required in the moving-window approach for isolating the local effects across models (see Sect. 2.3), and (2) the midlatitudes are a suitable test case because there the local effects on Tsurf have a different sign in the winter (DJF) and summer season (JJA) (Fig. 2).

In the areas considered, in most models (with the exception of CanESM2 and GFDL-CM3) Tsurf and T2 m respond similarly for changes in annual means (Fig. 3a). This is also true for the MPI-ESM – Tsurf and T2 m respond similarly when averaged over the respective regions (midlatitude areas (40–60 N) that have experienced intense deforestation (≥15 %) since 1860). A difference in response between Tsurf and T2 m becomes apparent also for the midlatitudes when analyzing seasons (DJF and JJA separately) instead of annual means. Almost all of the tested models show substantial differences between Tsurf and T2 m at a seasonal scale (Fig. 3b–c).

Figure 3Local effects on near-surface air temperature (T2 m) and surface temperature (Tsurf) for CMIP5 models. Values are averaged over midlatitude areas (40–60 N) that have experienced intense deforestation (≥15 %) since 1860. Positive values indicate a deforestation-induced warming. Each transparent marker denotes one combination of ensemble members from the historical and piControl experiments, respectively. The solid markers denote the mean values. The corresponding maps are shown in Figs. S6–S8. The local effects are isolated as in the study by .

In JJA (Fig. 3c), all but one model show a surface warming locally, with the Tsurf responding more strongly than T2 m by a factor of around 2 (Table S1). Only the HadGEM2-ES climate model is an outlier: there, Tsurf responds to deforestation with a local cooling (−0.13 K), which is not in agreement with observation-based studies , and in the HadGEM2-ES T2 m cools even more strongly (−0.27 K) than Tsurf. This inter-model comparison confirms that the local deforestation responses of Tsurf and T2 m quantitatively differ strongly in JJA, but in contrast with the MPI-ESM results shown above (Fig. 2), all models in Fig. 3 show the same sign of the JJA responses for Tsurf and T2 m.

In DJF (Fig. 3b), all but one model show a surface cooling locally, again with Tsurf responding stronger than T2 m in most models. An exception is the CanESM2 model, which locally responds to deforestation with a strong Tsurf warming and T2 m cooling. In some of the other climate models (CCSM4, CESM1-CAM5, NorESM1-M, all sharing the same land surface model CLM4), the Tsurf cooling is approximately twice the cooling of T2 m, analogous to the JJA response. In other models (GFDL-CM3, HadGEM2-ES, IPSL-CM5A-LR, MPI-ESM-LR), the two variables respond more similarly in DJF compared to JJA.

Overall, the inter-model comparison suggests that the quantitatively different response of Tsurf and T2 m is not specific to the MPI-ESM model. In agreement with previous studies , the inter-model spread in the temperature response in Fig. 3 is large (e.g., in JJA inter-model (excluding the HadGEM2) standard deviation 0.20 K, inter-model mean 0.38 K). However, the investigated models agree better concerning the ratio between the T2 m and Tsurf response (JJA inter-model (excluding the HadGEM2) standard deviation 0.11, inter-model mean 0.50). Both for DJF and JJA (Table S1), and for most of the investigated models, the ratio of changes in T2 m and Tsurf of 0.5 (range 0.35–0.65, excluding HadGEM2) is largely independent of the magnitude and sign of the Tsurf response. The temperature response of T2 m is between zero and the response of Tsurf with a ratio that depends on the exact way in which T2 m is calculated in the respective models (but most likely all models use a Monin–Obukhov approach). Further studies may investigate to what extent the calculation of T2 m differs across models.

4 Discussion and conclusions

This study shows that in climate models, surface temperature (Tsurf) and near-surface air temperature (T2 m) respond differently to deforestation. In the MPI-ESM, the nonlocal cooling in most regions (present also in locations that were not deforested) of Tsurf and T2 m is similar, while the local warming in the tropics and cooling in the higher latitudes (present only in locations that were deforested) differs. In the northern mid- and high latitudes, the annual mean local cooling of T2 m can be stronger than the local cooling of Tsurf, but in most regions Tsurf responds more strongly than T2 m. Across most models, the local effects of deforestation on Tsurf and T2 m in the midlatitudes differ by a factor of 2, for both local warming during summer and local cooling during winter.

This study illustrates that the conclusions concerning the effects of deforestation can depend on the considered temperature measure. The differences in magnitude and pattern between ΔTsurf and ΔT2 m (e.g., Lee et al.2011; Zhang et al.2014) obtained from observation-based studies largely agree with our findings (more cooling for T2 m than for Tsurf in the midlatitudes for the local effects in the MPI-ESM (see Figs. 1 and S2)), and thus our results make it seem plausible that the consideration of different temperature measures can explain some of the discrepancies between the satellite-based and in situ-based studies. A consistent comparison between satellite-based and in situ-based studies can be challenging because they may report different variables. Because of the heterogeneous emissivity of the land surface , satellite-based data sets usually report changes in radiometric surface temperature, which represents a combination of temperature at the top of the vegetation and the soil (through gaps in the canopy). Satellite-based direct estimates of air temperature (based on the intensity of upwelling microwave radiation from atmospheric oxygen) are available only for broad vertical layers of the atmosphere and at a coarse spatial scale . Instead of direct observations, air temperature was derived from surface temperature by empirical methods or process-oriented models (i.e., by solving the surface energy budget) . More direct observational investigations on the effects of deforestation on air temperature were based on recordings from weather stations and FLUXNET towers, which measure temperature at different heights. For instance, weather stations, e.g., in forest clearings, recorded temperatures at a height of between 1.2 and 2.0 m above ground level (WMO2008), while FLUXNET sites recorded temperatures typically 2–15 m above forest canopies . The different measurement height may lead to systematic differences because of the steep vertical temperature profile that develops near the surface under stable atmospheric conditions (e.g., at night) especially over open land . In contrast to satellite-based products, which are available at a high spatial resolution, the spatial distribution of FLUXNET towers and weather stations is biased toward developed countries and there is a relatively poor geographical coverage of rural areas in developing countries where deforestation has occurred recently . To perform a meaningful comparison, near-surface air temperature would have to be available at the same height above canopy top (preferably multiple heights) for the various land cover types and with a good geographical coverage.

The comparison of deforestation effects in observations and climate models is even more challenging. First, the respective variables in the models are only a proxy of the variables that were recorded in observation-based data sets (see Sect. 2.2 for the MPI-ESM). Second, model-based studies usually analyzed the combination of local and nonlocal effects (especially relevant for simulations of large-scale deforestation where the nonlocal effects can be large), while observation-based studies only analyzed local effects, for which Tsurf and T2 m respond differently (Fig. 1). Any nonlocal effects are excluded from the observations because possible nonlocal effects are present both in forest locations and nearby open land, and thus the nonlocal effects cancel out when looking at the difference between forests and open land, which is acknowledged by these studies . Note that Earth system models consider further climate effects when simulating deforestation-induced releases of land carbon into the atmosphere (e.g., Pongratz et al.2010; Le Quéré et al.2017). Because CO2 is a well-mixed greenhouse gas, the resulting warming can be expected to act essentially nonlocally and likely influences surface and air temperature similarly.

The different response of surface temperature (Tsurf) and near-surface air temperature (T2 m) is relevant for climate policies. Strategies that aim at adapting locally to warming air temperature may focus on perceived temperature and thus T2 m, but this study shows that the local effects on T2 m may substantially differ from those on Tsurf. Our results for the MPI-ESM suggest that the difference between T2 m and Tsurf is particularly strong for mean daily maximum temperature (see Fig. 2). Further studies may investigate whether this is also true for other climate models and observation-based data sets. Consequently, strategies in the agricultural sector that aim at adapting locally to warming soil and canopy temperatures may focus on the local effects on surface temperature because this variable is relevant for the organisms that live there. On the other hand, for international policies that aim at mitigating global warming, what matters is not only temperature at the location of deforestation but also in nearby and in remote regions. Thus, international policies may additionally consider the nonlocal effects. For the nonlocal effects, the responses of Tsurf and T2 m are rather similar (Fig. 1), and a distinction between the two temperature measures is therefore less relevant. To sum up, this study emphasizes that the local biogeophysical effects of deforestation influence Tsurf and T2 m differently, and thus, a careful choice based on the respective application has to be made regarding whether a study should focus on changes in surface temperature or near-surface air temperature.

Data availability
Data availability.

Data, scripts used in the analysis, and other supporting information that may be useful in reproducing the authors' work can be found here: http://hdl.handle.net/21.11116/0000-0002-CCA1-2 .

Supplement
Supplement.

Author contributions
Author contributions.

JW, CHR, and JP designed the research; JW performed the simulations with the MPI-ESM, analyzed the data, and drafted the paper. All authors contributed to the interpretation of the data and revision of the paper.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

Our simulations were performed at the German Climate Computing Center (DKRZ). We want to thank all groups who provided data for the inter-model comparison and Vivek Arora for discussions of the results for the CanESM2.

Financial support
Financial support.

This research was supported by the Deutsche Forschungsgemeinschaft (grant no. PO1751). Paul C. Stoy acknowledges support by the US National Science Foundation (grant nos. DEB-1552976 and EF-1702029).

The article processing charges for this open-access
publication were covered by the Max Planck Society.

Review statement
Review statement.

This paper was edited by Somnath Baidya Roy and reviewed by four anonymous referees.

References

Alkama, R. and Cescatti, A.: Biophysical climate impacts of recent changes in global forest cover, Science, 351, 600–604, https://doi.org/10.1126/science.aac8083, 2016. a, b, c, d, e, f, g, h

Bala, G., Caldeira, K., Wickett, M., Phillips, T. J., Lobell, D. B., Delire, C., and Mirin, A.: Combined climate and carbon-cycle effects of large-scale deforestation, P. Natl. Acad. Sci. USA, 104, 6550–6555, https://doi.org/10.1073/pnas.0608998104, 2007. a

Baldocchi, D.: How will land use affect air temperature in the surface boundary layer? Lessons learned from a comparative study on the energy balance of an oak savanna and annual grassland in California, USA, Tellus B, 65, 19994, https://doi.org/10.3402/tellusb.v65i0.19994, 2013. a

Bathiany, S., Claussen, M., Brovkin, V., Raddatz, T., and Gayler, V.: Combined biogeophysical and biogeochemical effects of large-scale forest cover changes in the MPI earth system model, Biogeosciences, 7, 1383–1399, https://doi.org/10.5194/bg-7-1383-2010, 2010. a

Boisier, J. P., de Noblet-Ducoudré, N., Pitman, A. J., Cruz, F. T., Delire, C., van den Hurk, B. J. J. M., van der Molen, M. K., Müller, C., and Voldoire, A.: Attributing the impacts of land-cover changes in temperate regions on surface temperature and heat fluxes to specific causes: Results from the first LUCID set of simulations, J. Geophys. Res.-Atmos., 117, 1–16, https://doi.org/10.1029/2011JD017106, 2012. a, b

Bonan, G. B.: Forests and climate change: forcings, feedbacks, and the climate benefits of forests, Science, 320, 1444–1449, https://doi.org/10.1126/science.1155121, 2008. a

Bright, R. M., Davin, E. L., O'Halloran, T. L., Pongratz, J., Zhao, K., and Cescatti, A.: Local temperature response to land cover and management change driven by non-radiative processes, Nat. Clim. Change, 7, 296–302, https://doi.org/10.1038/nclimate3250, 2017. a, b, c, d

Campbell, G. S. and Norman, J. M.: An Introduction to Environmental Biophysics, Springer-Verlag, New York, https://doi.org/10.1007/978-1-4612-1626-1, 1998. a

Claussen, M., Brovkin, V., and Ganopolski, A.: Biogeophysical versus biogeochemical feedbacks of large-scale land cover change, Geophys. Res. Lett., 28, 1011–1014, 2001. a

Davin, E. L. and de Noblet-Ducoudré, N.: Climatic impact of global-scale deforestation: radiative versus nonradiative processes, J. Climate, 23, 97–112, https://doi.org/10.1175/2009JCLI3102.1, 2010. a, b, c, d

De Frenne, P., Rodriguez-Sanchez, F., Coomes, D. A., Baeten, L., Verstraeten, G., Vellend, M., Bernhardt-Romermann, M., Brown, C. D., Brunet, J., Cornelis, J., Decocq, G. M., Dierschke, H., Eriksson, O., Gilliam, F. S., Hedl, R., Heinken, T., Hermy, M., Hommel, P., Jenkins, M. A., Kelly, D. L., Kirby, K. J., Mitchell, F. J. G., Naaf, T., Newman, M., Peterken, G., Petrik, P., Schultz, J., Sonnier, G., Van Calster, H., Waller, D. M., Walther, G.-R., White, P. S., Woods, K. D., Wulf, M., Graae, B. J., and Verheyen, K.: Microclimate moderates plant responses to macroclimate warming, P. Natl. Acad. Sci. USA, 110, 18561–18565, https://doi.org/10.1073/pnas.1311190110, 2013. a

De Frenne, P., Zellweger, F., Rodríguez-Sánchez, F., Scheffers, B. R., Hylander, K., Luoto, M., Vellend, M., Verheyen, K., and Lenoir, J.: Global buffering of temperatures under forest canopies, Nature Ecology & Evolution, 3, 744–749, https://doi.org/10.1038/s41559-019-0842-1, 2019. a

de Noblet-Ducoudré, N., Boisier, J.-P., Pitman, A., Bonan, G. B., Brovkin, V., Cruz, F., Delire, C., Gayler, V., van den Hurk, B. J. J. M., Lawrence, P. J., van der Molen, M. K., Müller, C., Reick, C. H., Strengers, B. J., and Voldoire, A.: Determining robust impacts of land-use-induced land cover changes on surface climate over North America and Eurasia: Results from the first set of LUCID experiments, J. Climate, 25, 3261–3281, https://doi.org/10.1175/JCLI-D-11-00338.1, 2012. a, b

Deser, C., Knutti, R., Solomon, S., and Phillips, A. S.: Communication of the role of natural variability in future North American climate, Nat. Clim. Change, 2, 775–779, https://doi.org/10.1038/nclimate1562, 2012. a

Devaraju, N., Bala, G., and Modak, A.: Effects of large-scale deforestation on precipitation in the monsoon regions: Remote versus local effects, P. Natl. Acad. Sci. USA, 112, 3257–3262, https://doi.org/10.1073/pnas.1423439112, 2015. a, b

Duveiller, G., Hooker, J., and Cescatti, A.: The mark of vegetation change on Earth's surface energy balance, Nat. Commun., 9, 1–12, https://doi.org/10.1038/s41467-017-02810-8, 2018. a, b, c, d, e, f, g, h

ECMWF Research Department: Research Manual 3, ECMWF forecast model, Physical parametrization, 3rd Edition, ECMWF, Reading, England, 1991. a, b, c

Findell, K. L., Knutson, T. R., and Milly, P. C. D.: Weak simulated extratropical responses to complete tropical deforestation, J. Climate, 19, 2835–2850, https://doi.org/10.1175/JCLI3737.1, 2006. a

Geleyn, J.-F.: Interpolation of wind, temperature and humidity values from model levels to the height of measurement, Tellus 40A, 4, 347–351, 1988. a, b

Gibbard, S., Caldeira, K., Bala, G., Phillips, T. J., and Wickett, M.: Climate effects of global land cover change, Geophys. Res. Lett., 32, 1–4, https://doi.org/10.1029/2005GL024550, 2005. a

Giorgetta, M. A., Roeckner, E., Mauritsen, T., Bader, J., Crueger, T., Esch, M., Rast, S., Kornblueh, L., Schmidt, H., Kinne, S., Hohenegger, C., Möbis, B., Krismer, T., Wieners, K.-H., and Stevens, B.: The atmospheric general circulation model ECHAM6 – Model description, Tech. Rep. 135, Max-Planck-Institute for Meteorology, Hamburg, 2013. a

Hansen, M. C., Potapov, P. V., Moore, R., Hancher, M., Turubanova, S. A., and Tyukavina, A.: High-resolution global maps of 21st-century forest cover change, Science, 342, 850–853, https://doi.org/10.1126/science.1244693, 2013. a

Hou, P., Chen, Y., Qiao, W., Cao, G., Jiang, W., and Li, J.: Near-surface air temperature retrieval from satellite images and influence by wetlands in urban region, Theor. Appl. Climatol., 111, 109–118, https://doi.org/10.1007/s00704-012-0629-7, 2013. a

Jin, M. and Dickinson, R. E.: Land surface skin temperature climatology: Benefitting from the strengths of satellite observations, Environ. Res. Lett., 5, 044004, https://doi.org/10.1088/1748-9326/5/4/044004, 2010. a

Jones, A. D., Collins, W. D., and Torn, M. S.: On the additivity of radiative forcing between land use change and greenhouse gases, Geophys. Res. Lett., 40, 4036–4041, https://doi.org/10.1002/grl.50754, 2013. a

Lague, M. and Swann, A.: Progressive Midlatitude Afforestation: Impacts on Clouds, Global Energy Transport, and Precipitation, J. Climate, 29, 5561–5573, https://doi.org/10.1175/JCLI-D-15-0748.1, 2016. a

Leclerc, M. and Foken, T.: Footprints in Micrometeorology and Ecology, Springer Berlin Heidelberg, https://doi.org/10.1007/978-3-642-54545-0, 2014. a, b

Lee, X., Goulden, M. L., Hollinger, D. Y., Barr, A., Black, T. A., Bohrer, G., Bracho, R., Drake, B., Goldstein, A., Gu, L., Katul, G., Kolb, T., Law, B. E., Margolis, H., Meyers, T., Monson, R., Munger, W., Oren, R., Paw U, K. T., Richardson, A. D., Schmid, H. P., Staebler, R., Wofsy, S., and Zhao, L.: Observed increase in local cooling effect of deforestation at higher latitudes, Nature, 479, 384–387, https://doi.org/10.1038/nature10588, 2011. a, b, c, d

Lejeune, Q., Seneviratne, S. I., and Davin, E. L.: Historical land-cover change impacts on climate: Comparative assessment of LUCID and CMIP5 multimodel experiments, J. Climate, 30, 1439–1459, https://doi.org/10.1175/JCLI-D-16-0213.1, 2017. a, b

Lejeune, Q., Davin, E. L., Gudmundsson, L., Winckler, J., and Seneviratne, S. I.: Historical deforestation increased the risk of heat extremes in northern mid-latitudes, Nat. Clim. Change, 8, 386–390, https://doi.org/10.1038/s41558-018-0131-z, 2018. a, b

Le Quéré, C., Andrew, R. M., Canadell, J. G., Sitch, S., Korsbakken, J. I., Peters, G. P., Manning, A. C., Boden, T. A., Tans, P. P., Houghton, R. A., Keeling, R. F., Alin, S., Andrews, O. D., Anthoni, P., Barbero, L., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Currie, K., Delire, C., Doney, S. C., Friedlingstein, P., Gkritzalis, T., Harris, I., Hauck, J., Haverd, V., Hoppema, M., Klein Goldewijk, K., Jain, A. K., Kato, E., Körtzinger, A., Landschützer, P., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Melton, J. R., Metzl, N., Millero, F., Monteiro, P. M. S., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., O'Brien, K., Olsen, A., Omar, A. M., Ono, T., Pierrot, D., Poulter, B., Rödenbeck, C., Salisbury, J., Schuster, U., Schwinger, J., Séférian, R., Skjelvan, I., Stocker, B. D., Sutton, A. J., Takahashi, T., Tian, H., Tilbrook, B., van der Laan-Luijkx, I. T., van der Werf, G. R., Viovy, N., Walker, A. P., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2016, Earth Syst. Sci. Data, 8, 605–649, https://doi.org/10.5194/essd-8-605-2016, 2016. a, b

Li, Y., Zhao, M., Motesharrei, S., Mu, Q., Kalnay, E., and Li, S.: Local cooling and warming effects of forests based on satellite observations, Nat. Commun., 6, 1–8, https://doi.org/10.1038/ncomms7603, 2015. a, b, c, d, e, f, g

Li, Y., De Noblet-Ducoudré, N., Davin, E. L., Motesharrei, S., Zeng, N., Li, S., and Kalnay, E.: The role of spatial scale and background climate in the latitudinal temperature response to deforestation, Earth Syst. Dynam., 7, 167–181, https://doi.org/10.5194/esd-7-167-2016, 2016. a

Luyssaert, S., Schulze, E.-D., Börner, A., Knohl, A., Hessenmöller, D., Law, B. E., Ciais, P., and Grace, J.: Old-growth forests as global carbon sinks, Nature, 455, 213–215, https://doi.org/10.1038/nature07276, 2008. a

Luyssaert, S., Marie, G., Valade, A., Chen, Y.-y., Djomo, S. N., Ryder, J., Otto, J., Naudts, K., Lansø, A. S., Ghattas, J., and Mcgrath, M. J.: Trade-offs in using European forests to meet climate objectives, Nature, 562, 259–262, https://doi.org/10.1038/s41586-018-0577-1, 2018. a

Malyshev, S., Shevliakova, E., Stouffer, R. J., and Pacala, S. W.: Contrasting local versus regional effects of land-use-change-induced heterogeneity on historical climate: Analysis with the GFDL Earth system model, J. Climate, 28, 5448–5469, https://doi.org/10.1175/JCLI-D-14-00586.1, 2015. a

Meier, R., Davin, E. L., Lejeune, Q., Hauser, M., Li, Y., Martens, B., Schultz, N. M., Sterling, S., and Thiery, W.: Evaluating and improving the Community Land Model's sensitivity to land cover, Biogeosciences, 15, 4731–4757, https://doi.org/10.5194/bg-15-4731-2018, 2018. a, b

Melo-Aguilar, C., González-Rouco, J. F., García-Bustamante, E., Navarro-Montesinos, J., and Steinert, N.: Influence of radiative forcing factors on ground–air temperature coupling during the last millennium: implications for borehole climatology, Clim. Past, 14, 1583–1606, https://doi.org/10.5194/cp-14-1583-2018, 2018. a

Pitman, A. J., de Noblet-Ducoudré, N., Cruz, F. T., Davin, E. L., Bonan, G. B., Brovkin, V., Claussen, M., Delire, C., Ganzeveld, L., Gayler, V., van den Hurk, B. J. J. M., Lawrence, P. J., van der Molen, M. K., Müller, C., Reick, C. H., Seneviratne, S. I., Strengers, B. J., and Voldoire, A.: Uncertainties in climate responses to past land cover change: First results from the LUCID intercomparison study, Geophys. Res. Lett., 36, 1–6, https://doi.org/10.1029/2009GL039076, 2009. a, b

Pongratz, J., Reick, C. H., Raddatz, T., and Claussen, M.: A reconstruction of global agricultural areas and land cover for the last millennium, Global Biogeochem. Cy., 22, 1–16, https://doi.org/10.1029/2007GB003153, 2008. a

Pongratz, J., Reick, C. H., Raddatz, T., and Claussen, M.: Biogeophysical versus biogeochemical climate response to historical anthropogenic land cover change, Geophys. Res. Lett., 37, 1–5, https://doi.org/10.1029/2010GL043010, 2010. a, b

Ramankutty, N. and Foley, J. A.: Estimating historical changes in global land cover: Croplands from 1700 to 1992, Global Biogeochem. Cy., 13, 997–1027, 1999. a

Schultz, N. M., Lawrence, P. J., and Lee, X.: Global satellite data highlights the diurnal asymmetry of the surface temperature response to deforestation, J. Geophys. Res.-Biogeo., 122, 903–917, https://doi.org/10.1002/2016JG003653, 2017. a, b, c, d

Staiger, H., Laschewski, G., and Graetz, A.: The perceived temperature – a versatile index for the assessment of the human thermal environment. Part A: Scientific basics, Int. J. Biometeorol., 56, 165–176, https://doi.org/10.1007/s00484-011-0409-6, 2011. a

Stevens, B., Giorgetta, M., Esch, M., Mauritsen, T., Crueger, T., Rast, S., Salzmann, M., Schmidt, H., Bader, J., Block, K., Brokopf, R., Fast, I., Kinne, S., Kornblueh, L., Lohmann, U., Pincus, R., Reichler, T., and Roeckner, E.: Atmospheric component of the MPI-M Earth System Model: ECHAM6, J. Adv. Model. Earth Sy., 5, 146–172, https://doi.org/10.1002/jame.20015, 2013.  a

Swann, A. L. S., Fung, I. Y., and Chiang, J. C. H.: Mid-latitude afforestation shifts general circulation and tropical precipitation, P. Natl. Acad. Sci. USA, 109, 712–716, https://doi.org/10.1073/pnas.1116706108, 2012. a

Taylor, K. E., Stouffer, R. J., and Meehl, G. A.: An Overview of CMIP5 and the Experiment Design, B. Am. Meteorol. Soc., 93, 485–498, https://doi.org/10.1175/BAMS-D-11-00094.1, 2012. a

UNFCCC: The Cancun Agreements: Land use, land-use change and forestry in Report of the Conference of the Parties serving as the meeting of the Parties to the Kyoto Protocol on its sixth session, United Nations Framework Convention on Climate Change (UNFCCC), Cancun, 1–32, 2011. a

Vanden Broucke, S., Luyssaert, S., Davin, E. L., Janssens, I., and Lipzig, N.: New insights in the capability of climate models to simulate the impact of LUC based on temperature decomposition of paired site observations, J. Geophys. Res.-Atmos., 120, 1–20, https://doi.org/10.1002/2015JD023095, 2015. a

Von Engeln, A. and Bühler, S.: Temperature profile determination from microwave oxygen emissions in limb sounding geometry, J. Geophys. Res.-Atmos., 107, 1–15, https://doi.org/10.1029/2001JD001029, 2002. a

Winckler, J., Reick, C. H., and Pongratz, J.: Robust identification of local biogeophysical effects of land-cover change in a global climate model, J. Climate, 30, 1159–1176, https://doi.org/10.1175/JCLI-D-16-0067.1, 2017. a, b, c, d, e, f, g

Winckler, J., Reick, C. H., Lejeune, Q., and Pongratz, J.: Nonlocal effects dominate the global mean surface temperature response to the biogeophysical effects of deforestation, Geophys. Res. Lett., 46, 745–755, https://doi.org/10.1029/2018gl080211, 2019a. a, b, c, d, e, f, g, h, i, j

Winckler, J., Reick, C. H., Raddatz, T., Pongratz, J., Kauhs, C.: Supplementary material and data for the manuscript `Different response of surface temperature and air temperature to deforestation in climate models', MPG Publication Repository, http://hdl.handle.net/21.11116/0000-0002-CCA1-2 (last access: July 2019), 2019b. a

WMO: Guide to Meteorological Instruments and Methods of observation, vol. I & II, World Meteorological Organization (WMO), Geneva, Switzerland, ISBN 978-92-63-100085, 2008. a

Zhang, M., Lee, X., Yu, G., Han, S., Wang, H., Yan, J., Zhang, Y., Li, Y., Ohta, T., Hirano, T., Kim, J., Yoshifuji, N., and Wang, W.: Response of surface air temperature to small-scale land clearing across latitudes, Environ. Res. Lett., 9, 034002, https://doi.org/10.1088/1748-9326/9/3/034002, 2014. a, b, c

Zwiers, F. W. and von Storch, H.: Taking serial correlation into account in tests of the mean, J. Climate, 8, 336–351, https://doi.org/10.1175/1520-0442(1995)008<0336:TSCIAI>2.0.CO;2, 1995. a, b