RESEARCH ARTICLE The Duality of Reforestation Impacts on Surface and 10.1029/2019JG005543 Air Temperature Key Points: Kimberly A. Novick1 and Gabriel G. Katul2 • Evidence is mounting that reforestation cools the land surface 1O'Neill School of Public and Environmental Affairs, Indiana University, Bloomington, Bloomington, IN, USA, 2Nicholas in many places, but its been challenging to understand how School of the Environment, Duke University, Durham, NC, USA reforestation affects air temperature • A novel approach for detecting the fl in uence of land cover on multiple Abstract Evidence is mounting that temperate‐zone reforestation cools surface temperature (Tsurf), fl metrics of air temperature using ux mitigating deleterious effects of climate warming. While T drives many biophysical processes, air tower observations is presented surf • The analysis shows that temperature (Ta) is an equally important target for climate mitigation and adaptation. Whether reductions in reforestation in the southeastern Tsurf translate to reductions in Ta remains complex, fraught by several nonlinear and intertwined processes. United States cools the near‐surface ‐ – In particular, forest canopy structure strongly affects near surface temperature gradients, complicating air temperature by 1 3 °C during ‐ fl daytime but not nighttime cross site comparison. Here the in uence of reforestation on Ta is assessed by targeting temperature metrics that are less sensitive to local canopy effects. Specifically, we consider the aerodynamic temperature (Taero), estimated using a novel procedure that does not rely on the assumptions of Monin‐Obukhov similarity theory, as well as the extrapolated temperature into the surface layer (T ). The approach is tested with Correspondence to: extrap K. A. Novick, flux tower data from a grass field, pine plantation, and mature hardwood stand co‐located in the Duke Forest [email protected] (North Carolina, USA). During growing season daytime periods, Tsurf is 4–6 °C cooler, and Taero and near‐surface Textrap are 2–3 °C cooler, in the forests relative to the grassland. During the dormant season, Citation: daytime differences are smaller but still substantial. At night, differences in Taero are small, and near‐surface Novick, K. A., & Katul, G. G. (2020). Textrap is warmer over forests than grasslands during the growing season (by 0.5 to 1 °C). Finally, the The duality of reforestation impacts on influence of land cover on T at the interface between the surface and mixed layer is small. Overall, surface and air temperature. Journal of extrap Geophysical Research: Biogeosciences, reforestation appears to provide a meaningful opportunity for adaption to warmer daytime Ta in the 124. https://doi.org/10.1029/ southeastern United States, especially during the growing season. 2019JG005543 Plain Language Summary Reforestation—the process of reestablishing trees where they once Received 26 OCT 2019 dominated—has long been viewed as a strategy to remove CO2 from the atmosphere. Recently, attention Accepted 10 JAN 2020 has focused on understanding if reforestation also offers a direct temperature cooling benefit. By using more Accepted article online 13 MAR 2020 water (a cooling process) and increasing the transfer of heat energy away from the surface, forests may offer a meaningful opportunity for local climate mitigation and adaptation. Evidence is mounting that indeed, in the temperature and tropical zones, the surface of forests is cooler than grasslands and croplands. However, due to confounding effects of forest canopies on wind and temperature profiles near the surface, it has previously been hard to assess if forests also cool the air. Here we present a new approach that accounts for canopy effects, allowing for a more direct assessment of the potential for reforestation to cool near‐surface air temperature. Using a case study from the North Carolina Piedmont, we find that while the air cooling effect of forests is not a large as the surface cooling effect, it is still on the order of 2–3°C during summer daytime periods—times when the need for climate adaptation strategies are particularly pressing.
1. Introduction
Reforestation has long been viewed as an instrument for mitigating the pace of climate change, particularly in the temperate and tropical zone where much of the historical forest cover was lost to harvest within the last 200–300 years (Williams, 1989). This view is largely linked to the carbon sequestration potential of these
forests. Observations of the net ecosystem exchange of CO2 from regional and global networks of flux towers, as well as forest inventory data, reveal that forests in the temperate regions are indeed strong carbon sinks (Jung et al., 2011; Pan et al., 2011) and that even maturing temperate forests are capable of assimilating sub-
stantially more CO2 than expected from conventional ecological theory (Novick et al., 2015; Stoy et al., 2008).
©2020. American Geophysical Union. However, the fate of the future forest carbon sink is less certain (Friedlingstein et al., 2014). As atmospheric All Rights Reserved. CO2 continues to rise, forest carbon uptake potential may saturate due to a number of limitations. Some are
NOVICK AND KATUL 1of15 Journal of Geophysical Research: Biogeosciences 10.1029/2019JG005543
imposed by nutrient and energy availability (Baldocchi & Penuelas, 2019; Oren et al., 2001), while others are
intrinsic to leaf‐level photosynthesis and its saturating behavior with increased CO2. Moreover, expected increases in air temperature, drought, and insect and fire regimes will likely decrease the magnitude, and increase the variability, of forest carbon uptake across much of the world (Frank et al., 2015; McDowell & Allen, 2015; Wear & Coulston, 2015).
While this uncertainty about future forest CO2 uptake continues to motivate research, substantial atten- tion is now focused on the potential for forests to mitigate temperature conditions through alterations to the ecosystem energy balance. Theoretical links between land cover and temperature have long been recognized and incorporated into models (Avissar & Werth, 2005; Foley et al., 2003; Pielke et al., 1998; Raupach, 1991). Modeling work has shown that in the tropics, relatively large evapotranspiration causes forests to be cooler than nonforested ecosystems (Costa, 2005). On the other hand, in boreal cli- mates, relatively low forest albedo likely causes forests to be warmer than nonforested ecosystems (Lee et al., 2011; Swann et al., 2010). In the temperate zone, the overall impact of temperate reforestation on surface temperature was, for a long time, not clear (Bonan, 2008; South et al., 2011). Evaporative cool- ing and emitted longwave radiation from the surface both act to reduce surface temperature, whereas net shortwave radiative load acts to warm the surface. The sensible heat flux, whose efficiency varies with the mean wind and turbulence conditions overlying the surface, plays a dual role and may contri- bute to warming or cooling (Huang et al., 2015). With rapid advancements and proliferation of remote sensing products, observational evidence has emerged to suggest that the combined influence of increased sensible and latent heat in temperate forest ecosystems has an overall surface cooling effect that outweighs albedo‐driving warming by a magnitude of 1–2 °C, annually averaged, across a wide range of temperate ecosystems (Bright et al., 2017; Burakowski et al., 2018; Juang, Katul, et al., 2007; Wickham et al., 2012; Zhang et al., 2020). This work is encouraging, as it suggests a direct and substantial climate mitigation benefit of reforestation in the temperate zone that is mechanistically quite different from the carbon sequestration benefit. However, thus far, observation‐driven studies of land cover effects on microclimate have largely been focused on the response of surface temperature, and not the response of air temperature (e.g. Bright et al., 2017; Juang, Katul, et al., 2007; Luyssaert et al., 2014; Wickham et al., 2012; but see Baldocchi & Ma, 2013). This focus is not surprising for three reasons: (i) Unlike air temperature, surface temperature does not vary with height making it a more logical reference to compare land cover temperature patterns; (ii) likewise, surface tem- perature, representing the integrated radiometric temperature of all canopy elements, is closely coupled to the temperature experienced by foliage in dense canopies, or by microbes near the soil surface of sparse canopies, and is thus biologically relevant; and (iii) operationally, surface temperature is convenient to esti- mate from meteorological towers that report patterns of radiation, albedo, and energy fluxes necessary to attribute variations in surface temperature between ecosystems to specific mechanisms (Juang, Katul, et al., 2007; Lee et al., 2011; Luyssaert et al., 2014). However, when considering direct impacts of reforestation and other land cover changes on climate, air tem-
perature (hereafter Ta), as opposed to surface temperature (Tsurf), is an equally important metric. This rele- vance is certainly true for boundary layer dynamics and rainfall initiation (Juang, Katul, et al., 2007, Juang, Porporato, et al., 2007; Manoli et al., 2016; Siqueira et al., 2009) as well as a plethora of associated “hand‐shakes” between the climate system and the land surface (Baldocchi & Ma, 2013; Luyssaert et al., 2014). Ultimately, climate change is driven by long‐term increases in the temperature of the air (or kinetic temperature) due to increases in greenhouse gases and not by surface temperature.
Linking land cover impacts on Tsurf and subsequent changes in Ta is not straightforward. Since H transfers heat from the surface to the atmosphere, ecosystems with relatively cool surfaces may underlie relatively warm air (Baldocchi & Ma, 2013), and vice versa. These heat transfer mechanisms are further mediated by land cover‐driven changes in the height of the planetary boundary layer (Luyssaert et al., 2014), which is gen- erally greater over forests, creating more “room” for heat energy transferred from the surface to the atmo- sphere. Finally, air temperature gradients near the surface can be steep but also depend on the influence of rough canopy elements on near‐surface turbulence regimes. Thus, uncertainty can arise from a straightfor-
ward comparison of observed Ta above a forested and nonforested canopy if steps are not taken to control for the influence of canopy structure on near‐surface temperature gradients.
NOVICK AND KATUL 2of15 Journal of Geophysical Research: Biogeosciences 10.1029/2019JG005543
The objective of this work is develop and use a “coarse‐grained” model of how land cover change affects air temperature near the surface while ensuring proper matching with surface temperature. The proposed
approach is then applied to the specific question of how reforestation affects Ta in the southeastern United States, across seasons, and over the diurnal cycle. The study leverages eddy covariance flux tower observations from a co‐located grassland, pine forest, and hardwood forest in central North Carolina (USA), which together represent the three primary phases of secondary succession in the region (grassland to pine forest to mature deciduous forest). The proposed approach relies only on data commonly reported by meteorological towers and thus is transferable to studies on land cover change impacts on temperature in other parts of the world.
2. Approach 2.1. Study Sites The study sites were located in the Blackwood Division of the Duke Forest near Durham, North Carolina (35°58′41″N, 79°05′59″W, 163 m above sea level). The old field grassland, hereafter OF (AmeriFlux site code US‐Dk1), was dominated by the C3 grass Festuca arundinacea Schreb., with minor contributions from forbs and other C3 and C4 grass species (Novick et al., 2004). It was harvested at least once a year to prevent refor- estation. The pine forest, hereafter PP (AmeriFlux site code US‐Dk2), was established in 1983 following a clear cut and a burn. Pinus taeda L. (Loblolly pine) seedlings were planted at a 2.0‐m by 2.4‐m spacing with pine density reduced to approximately 1,100 trees/hectare when the site reached maturity. Canopy height increased from 16 m in 2001 to over 20 m in 2008. The hardwood forest, hereafter HW (AmeriFlux site code US‐Dk3), was classified as an uneven‐aged (90–110 year old) oak (Quercus)–hickory (Carya) forest. The stand was dominated by hickories (Carya tomentosa (Poir.) Nutt., C. glabra (P. Mill). Sweet.), yellow poplar (Liriodendron tulipifera L.), sweetgum (Liquidambar styraciflua L.), and oaks (Quercus alba L., Q. michauxii Nutt., Q. phellos L.). The forest was not managed after establishment, and mean canopy height was 25 m. All ecosystems have little topographic variation and lie on Enon silt loam. A clay pan at a depth of ~ 35–50 cm underlies the sites, thereby imposing similar constraints on root‐water access for both ecosystems. Because they are all co‐located to within 1 km, they experience nearly identical macroclimate conditions, character- ized by long‐term mean annual temperature and precipitation of 15.5 °C and 1,146 mm, respectively. More details on the study sites are available elsewhere (Novick et al., 2004; Stoy et al., 2008). Much of the analysis is conducted separately for the dormant season (Julian day of year [DOY] <100 or DOY >300) and growing season (150 < DOY < 270). Some analysis focuses exclusively on daytime (9:00–17:00 hr) or nighttime (2100:2400 and 0:600 hr). The analysis relied on data from 2005 to 2008, as this is the period of record for which the towers supported measurements of outgoing long‐wave radiation necessary to infer
Tsurf. Static canopy heights of 0.5, 19, and 25 m for the OF, PP, and HW sites, respectively, were assumed for the duration during this period. 2.2. Temperature Metrics: Definitions and Observations Approaches
2.2.1. Surface Temperature (Tsurf) This study considers four metrics for temperature at or near the surface (see Figure 1). The first is the radio-
metric surface temperature (Tsurf) or canopy “skin temperature” (Jin & Dickinson, 2010). Conceptually, it represents the aggregated temperature of solid canopy and soil elements projected to a single location in
the vertical dimension. The Tsurf is highly relevant for biophysical processes occurring within the canopy, including respiration and photosynthesis (Luyssaert et al., 2014). It can be inferred directly at the ecosystem scale from observations of outgoing longwave radiation from a solid surface using the Stefan‐Boltzmann law 1/4 —Tsurf =(σεsRL. out) —where σ is the Stefan‐Boltzmann constant and ϵs is the surface emissivity. Because of its dependency on solid surface properties, ϵs may be empirically related to albedo (Juang, Katul, et al., 2007) though no “causal explanations” are to be implied. Albedo itself can be estimated on a daily time step as the ratio of midday incident to outgoing short‐wave radiation measured from the towers, with midday being the time of day when the incident radiation is most orthogonal to the reflecting surface.
Multiple studies have investigated the question of where and why Tsurf varies between forests and grasslands in the study region (Burakowski et al., 2018; Juang, Katul, et al., 2007; Zhang et al., 2020). At the annual time
scale, the Tsurf of forests is substantially cooler than the Tsurf of grasslands, by approximately 1–2 °C in gen- eral, and by approximately 1 °C in the study sites (described later). The cooling effect of forests is enhanced
NOVICK AND KATUL 3of15 Journal of Geophysical Research: Biogeosciences 10.1029/2019JG005543
Figure 1. Conceptual overview of the key variables and processes. Over forests, the tower observations typically occur in the “roughness sublayer,” where wind and temperature dynamics are heavily influenced by canopy elements, and the actual Ta gradients (here shown for daytime conditions, yellow dashed line) diverge substantially from those predicted by Monin‐Obukhov Similarity Theory (MOST; orange profiles). In contrast, in the grassland, Ta is usually measured in the surface layer, where wind speeds are greater and vertical temperature gradients agree well with MOST predictions. The aerodynamic temperature (Taero) and the radiometric surface temperature (Tsurf) can both be conceptually surrogated to a representative source/sink height in the canopy (i.e., the sum of the roughness height for heat and the zero plane dis- placement, blue line). Air temperature extrapolated into the surface layer (Textrap) is expected to be well aligned with predictions from MOST.
in the summer when albedo differences between forests and grasslands are smaller but evapotranspiration (and the associated evaporative cooling) is enhanced over forests (Juang, Katul, et al., 2007; Zhang et al., 2020). Given this substantial body of preexisting work, this study will not investigate in detail the
mechanisms driving differences in Tsurf in the study sites here. Instead, this study focuses on the relation between Tsurf and the multiple “metrics” of air temperature. 2.2.2. Air Temperature (Ta) For our purposes, this is the temperature of the air measured at some height on a meteorological tower in
an aspirated housing unit. Differences in elevation at the Ta measurement height were corrected for using the barometric formula to transform raw Ta into potential temperature. As illustrated in Figure 1, the links between Ta, surface temperature and the temperature of the atmospheric surface layer are compli- cated by the presence of a “roughness sublayer”—the very lowest layer of the atmosphere where the flow and temperature statistics are affected by canopy elements. The height of the roughness sublayer remains a subject of inquiry but is usually approximated as 2–5 times the height of the canopy (Raupach & Thom, 1981), in agreement with higher‐order closure modeling studies as well experimental and simulation stu- dies of flow over complex terrain. The lower limit is generally associated with mean momentum exchange, whereas the upper limit is representative of scalar exchanges (Poggi & Katul, 2007; Siqueira & Katul, 2010).
In short stature grasslands, the observation height for Ta usually exceeds the roughness sublayer height and is instead located in the so‐called “surface layer.” The surface layer represents the lower ~10% of atmospheric boundary layer within which gradients of mean air temperature and mean wind speed are represented by their logarithmic shapes corrected for thermal stratification using Monin‐Obukhov Similarity Theory (here-
after MOST; Monin & Obukhov, 1954). In contrast, except in the case of very tall towers, the Ta from forested flux towers is observed at a height that almost always falls within the roughness sublayer (as defined above).
More importantly, flow immediately above rough surfaces experiences higher friction velocity (u*), higher turbulent diffusivity, and smaller vertical gradients when compared to surfaces characterized by smaller roughness.
NOVICK AND KATUL 4of15 Journal of Geophysical Research: Biogeosciences 10.1029/2019JG005543
2.2.3. Aerodynamic Air Temperature (Taero) This variable represents the air temperature at the “at the apparent source/sink of heat within the canopy” (Chehbouni et al., 2001), which is an idealized plane often specified as the combined height of a zero‐plane
displacement (zd) for momentum or heat and a roughness length for heat (zo,h). Both variables are key MOST parameters. The magnitude of Taero cannot be measured directly, but it is constrained to be within the bounds of Ta measured near the surface and Tsurf (Jin & Dickinson, 2010, and see Figure 1). In bulk aerody- namic representation, it is a key driver of sensible heat flux (H) toward (or away from) the surface, given as H ¼ ðÞ′ − ′ ; ga;h z Taero Ta z (1) cp
where cp is the specific heat capacity for dry air at constant pressure and ga,h(z′) is the aerodynamic conduc- ′ tance at height z = z − zd. The ga,h(z′) can be expressed as
ðÞ′ 2 ðÞ¼′ hi Uzhik ; ga;h z (2) z′ z′ ln − ΨmðÞζ ln − ΨhðÞζ zo;m zo;h
where U is the mean wind speed at z′, k=0.4 is the von Karman constant, ρ is the mean air density, and zo,m is the roughness length for momentum. The Ψm(ζ) and Ψh(ζ) are stability correction functions that depend ′ on the atmospheric stability parameter ζ = z /L, which is a measure of buoyancy production (or destruction) to mechanical production of turbulent kinetic energy (TKE), and L is the Obukhov length. Their form is pre- sented in Appendix A for unstable (i.e., when buoyancy is a source of TKE, ζ < 0) and stably stratified (i.e., when buoyancy is a sink of TKE, ζ > 0) atmosphere.
The zd can be reasonably estimated as 0.6–0.7 times the height of the canopy for dense canopies when treated as the centroid of the drag force within the canopy (Jackson, 1981; Thom, 1971). The zo,m is commonly approximated as 0.1 times the height of the canopy (Campbell & Norman, 1998); alternatively, it can be inferred with reasonable confidence from the diabatic profile equation for wind speed: u* z′ UzðÞ¼′ ln − ΨmðÞζ ; (3) 0:4 zo;m
where u* is friction velocity. Finally, by combining equations (2) and (3), the form of the diabatic profile
equation for temperature can be obtained that depends on the friction velocity, but not zo,m or Ψm(ζ).