Heat stress is one of the most severe climate threats to human society in a future warmer world. The situation is further exacerbated
in urban areas by urban heat islands (UHIs). Because the majority of
world's population is projected to live in cities, there is a pressing need
to find effective solutions for the heat stress problem. We use a climate
model to investigate the effectiveness of various urban heat mitigation
strategies: cool roofs, street vegetation, green roofs, and reflective
pavement. Our results show that by adopting highly reflective roofs, almost
all the cities in the United States and southern Canada are transformed into
“white oases” – cold islands caused by cool roofs at midday, with an
average oasis effect of
Heat stress associated with climate change is projected to cause a substantial increase in human mortality (Patz et al., 2005; Huang et al., 2011) and a large reduction in workplace productivity (Dunne et al., 2013; Zander et al., 2015). These risks are further amplified for urban residents because of the urban heat island (UHI) effect, a phenomenon in which surface temperatures are higher in urban areas than in surrounding rural areas (Grimmond, 2007). Because more than 50 % of the world's population currently lives in cities, and that number is projected to increase to 70 % by year 2050 (Heilig, 2012), there is a pressing need to find effective solutions for urban heat stress. In recent years, urban climate agendas have broadened beyond carbon management, which brings marginal heat relief to urban residents (Revi et al., 2014), to include urban climate adaptation. It is now recognized that in addition to the traditional emphasis on preparedness to cope with heat stress (Stone et al., 2012), urban adaptation should include active modifications to the urban landscape to reduce urban temperatures (Rizwan et al., 2008; Castleton et al., 2010; Santamouris, 2013; Bowler et al., 2010; Jacobson, 2009).
Field experiments have demonstrated the cooling benefits of street vegetation (Shashua-Bar and Hoffman, 2000; Mackey et al., 2012), bright pavement (Rosenfeld et al., 1995; Mackey et al., 2012), green roofs (Jim and Peng, 2012), solar roofs (Dominguez et al., 2011), and cool roofs (Ismail et al., 2011) at the scale of individual buildings. Specifically, street vegetation and green roofs provide localized cooling through enhanced evaporation, and for street vegetation by increasing surface roughness and therefore convection efficiency. Bright pavement and cool roofs, both characterized by a higher albedo than traditional materials, cool the urban environment by reflecting a higher portion of solar radiation away from the surface. Solar roofs produce localized cooling through their conversion of solar energy to electricity.
Although it is financially challenging to conduct field experiments on these methods at the city scale, modeling studies have evaluated the potential of citywide implementation of these methods on temperature reduction. A majority of the modeling studies deploy a mesoscale weather forecast model, such as the Weather Research Forecasting Model (WRF; Skamarock et al., 2008), and model simulations are performed for a business-as-usual case and for one or more mitigation strategies, including tree planting (Stone et al., 2013, 2014), green roof (Li et al., 2014; Stone et al., 2014), cool roof (Georgescu et al., 2014), and rooftop solar panels (Salamanca et al., 2016).
The high computational demand of mesoscale modeling limits long-term and large-scale simulations, with most of these studies restricting their analysis to high-impact heatwaves in a single city (Rosenzweig et al., 2009; Li et al., 2014). Recently, studies using a mesoscale modeling approach have begun to scale up the geographic scope of analysis, comparing and evaluating urban heat mitigation methods across multiple (up to 3) cities in contrasting regional climates (Stone et al., 2014; Vargo et al., 2016) and even over the continental United States (Georgescu et al., 2014).
A second modeling strategy, using subgrid output from a global climate model (GCM), overcomes the computational constraints of mesoscale models, allowing for the continental or global-scale analysis of UHI and heat mitigation strategies over climatological timescales (Oleson et al., 2010c; Hu et al., 2016; Jacobson and Ten Hoeve, 2012). Although climate models cannot resolve local-scale processes, such as advection between adjacent urban and rural land in the same grid cell, they can capture the urban effects on the large-scale dynamics of the atmosphere. Additionally, GCMs can simulate the interactions of greenhouse gas (GHG)-induced warming and the biophysical drivers of urban heat islands and heat mitigation methods. As urban heat mitigation is becoming an integral part of the climate adaptation agenda, GCMs can provide useful knowledge to city planners about the effectiveness of various mitigation strategies across a range of climate conditions and under different climate scenarios.
In this study, we use an urban climate model to assess the effectiveness of urban climate mitigation methods in future warmer climates. Our specific objectives are to (1) investigate the urban heat island intensity under current and future climate scenarios for three climate regions in the US and southern Canada, (2) quantify the effectiveness of several mitigation strategies (street vegetation, green roof, cool roof and bright pavement) on offsetting urban warming, and (3) estimate aggregated temperature reduction potential of these strategies using a UHI mitigation wedge approach. Our research complements the studies by Georgescu et al. (2014), Rosenzweig et al. (2009), and Stone et al. (2014) in that we are also interested in quantifying the effectiveness of multiple UHI mitigation strategies, both individually and collectively, across a diverse range of climate conditions. Here, we expand the analysis to several climate scenarios and a large number of cities (57 cities) to understand the interactions between the UHI biophysical drivers, GHG-induced warming, and influence of local background climate. We propose a new method to assess the mitigation strategies, which is based on a theoretical understanding of the surface energy balance and is unconstrained by computational demand.
We used a global climate model to simulate the UHI and to quantify the cooling potential of urban heat mitigation strategies: cool roofs, street vegetation, green roofs, and reflective pavement. As will be discussed in more detail in the following subsections, the simplified representation of urban areas in the model does not directly allow for online assessment of each of the mitigation strategies. Therefore, in addition to the online simulations to assess the mitigation potential of cool roofs, we used an offline attribution method to quantify the mitigation potential of green roofs and reflective pavement as well as a two-member interpolation method to quantify the cooling potential of street vegetation (Table 1) for three climate scenarios. The offline calculations are based on the diagnostic variables of the surface energy balance produced by the model.
The simulations run and the methods used to assess the temperature mitigation strategies. The online simulation method calculated the mitigation potential of cool roofs as the difference between the WHT and CTR simulations (WHT – CTR). The offline attribution and two-endmember interpolation methods used diagnostic data from the CTR simulations.
We used the Community Earth System Model (CESM; Hurrell et al., 2013) to simulate the UHI for cities in the United States and southern Canada. The land surface processes are represented by the Community Land Model (CLM, version 4.0), which is the land component of CESM (Oleson et al., 2010b). In CLM 4.0, the land surface heterogeneity is represented as a nested hierarchy of subgrid levels. The model grid cell consists of up to five land units: vegetated land, glacier, wetland, urban, and lake, all of which are driven by the same atmospheric forcing. In the CESM architecture, fluxes and state variables are computed at each land unit level and then area-weighted to grid cell means before being passed to the atmosphere model. We use data from the urban and vegetated land units in this analysis to represent urban and rural land cover types, respectively. The urban land unit is represented in a conceptual canyon structure that consists of a roof, sunlit wall, shaded wall, and pervious (bare soil) and impervious (road, sidewalk, and parking lot) ground (Oleson et al., 2010a). The thermal (such as heat capacity and thermal conductivity), radiative (such as albedo and emissivity), and morphological (such as height-to-width ratio, roof areal fraction, average building height, and pervious ground fraction) characteristics of these canyon components are provided by Jackson et al. (2010) for each grid cell. The vegetated land unit is comprised of up to 15 natural plant functional types and bare soil.
The CLM was run under three climate scenarios: current climate,
Representative Concentration Pathway (RCP) 4.5, and RCP8.5. These three
simulations comprise the control (CTR) group of simulations, where all the
urban parameters were kept as the default values prescribed by the model. For
the current climate, the model was run for 33 years after a 60-year spin-up,
driven by a reconstructed climatology from 1972 to 2004 (Qian et al., 2006).
For the two future scenarios, the model was forced by atmospheric outputs
from fully coupled runs of the CESM and run for 96 years from 2005 to 2100
after a 600 year spin-up. Dynamic land surface input data corresponding to
each future scenario are used in the simulations. It should be noted here
that although the vegetated land unit changes dynamically in the climate
scenarios, the urban land is kept fixed in the default configuration of the
current version of CESM. This model setup is a shortcut to the fully coupled
mode and can be considered as a simplified retrieval of the surface climate
variables from the fully coupled runs. The results should be nearly identical
to those obtained from the coupled simulations because the impacts of
large-scale feedbacks represented in a coupled CESM, such as large-scale
dynamics and ocean–air feedbacks, are preserved by the atmospheric variables.
All the simulations were run at a horizontal resolution of
This model configuration cannot simulate the dynamic impact of changes to the urban land on the atmosphere. However, because the urban land unit in CLM comprises only a small areal fraction of each grid cell, changes in urban temperature would lead to negligible changes in grid cell mean temperature. Therefore, the dynamic feedbacks between changes in urban land and the atmosphere are negligible. Studies that fully couple CLM to the atmosphere conclude that modification of the urban albedo has negligible effects on the regional and global climate (Hu et al., 2016; Zhang et al., 2016).
We used 33 years (1972–2004) of data in the current climate and
the last 30 years (2071–2100) of data in the two future scenarios to
compute climatological mean temperatures and the UHI intensity. We analyzed
the UHI intensity for 57 selected cities, representing three Köppen–Geiger
climate zones: temperate climate (24 cities, eastern and southern US),
continental climate (20 cities, northern US and southern Canada), and dry
climate (13 cities, arid and semiarid western US). These cities comprise
less than 20 % of the grid cell area. Model outputs at 13:00 and 01:00
local time were used to represent daytime and nighttime conditions,
respectively. These two times were selected for two reasons. First,
they are close to the times of the daily maximum and minimum temperature, thus
giving a better representation of the diurnal range of the UHI intensity
(
In this analysis, we use the surface temperature instead of the screen-height air temperature in our assessment of UHI and the temperature mitigation potential of cooling strategies. In other words, the UHI intensity is defined as the surface UHI, which is the difference in radiative surface temperature between urban and nonurban subgrid land units in the grid cell where the city is located. This is different from the air UHI which is defined using screen-height air temperature. These two UHI definitions differ in a variety of aspects. First, the air UHI can be affected by the local air and landscape conditions and thus is susceptible to the inhomogeneous urban landscape. The surface UHI, instead, is a spatial average. Although the surface temperature does not necessarily match the human experience of warmth while walking across the rural–urban boundary, it is more stable as a metric of the city-scale microclimate and can be directly validated against satellite observations over large areas, especially when compared across different cities over diverse local climate conditions. Second, the air UHI is more tightly related to heat stress and heat exposure assessment (Oleson et al., 2015) than surface UHI. However, the surface UHI has a firm theoretical basis that underpins the wedge method for assessment of multiple mitigation strategies. It has been shown from the surface energy balance principle that the overall surface UHI can be estimated by linear supposition of different biophysical contributions (Zhao et al., 2014). The linear supposition property, that the overall temperature reduction due to simultaneous implementation of several mitigation methods can be approximated by the sum of cooling benefits from implementation of these methods one at a time, is the underlying basis for the UHI mitigation wedge approach used in this study. In addition, this linear supposition seems to be a robust property for the air UHI as well. We will present evidence from published UHI mitigation studies using the air UHI (Taha et al., 1997; Georgescu et al., 2014; Stone et al., 2014; Salamanca et al., 2016).
Defining UHI by surface temperature allows us to estimate the contribution
of different biophysical mechanisms to the overall surface UHI from surface
energy balance principles (Zhao et al., 2014). In addition to
To investigate the mitigation potential of cool roofs, an additional set of
simulations was conducted for each of the three climate scenarios (including
current climate, RCP 4.5, and RCP 8.5). In this group of simulations (white roofs, here abbreviated as WHT),
the roof albedo was increased to 0.88, the value of the US EPA Energy Star
Solareflect coating material after 3 years of wear and aging (Georgescu
et al., 2014). This value is slightly lower than the initial albedo of 0.89,
indicating the material's high capability of maintaining high solar
reflectance
(
Attribution of summer mean UHI intensity during 2071–2100
under the RCP 4.5 scenario for the control (CTR) run.
The morphological properties were kept unchanged in the CTR and WHT simulations. Of the 57 selected cities, the average prescribed roof areal fraction in CLM 4.0 is 48.8 %, the average default roof albedo is 0.29, and the average citywide albedo is 0.18. After the deployment of white roofs, the average citywide albedo is increased to 0.47.
To help the reader visualize the changes brought by the cool roofs to the
urban landscape, we created two animations using a 3-D point cloud and image
mosaic dataset acquired by an unmanned aerial vehicle over an urban
neighborhood in Switzerland (
It is not possible to directly evaluate the cooling benefit of street
vegetation with CLM because vegetation is not explicitly represented in the
urban land unit in CLM. Here, we used a simple two-endmember interpolation
method to calculate the surface temperature change
This method yields the net cooling of all the biophysical effects associated with street vegetation, including changes in albedo, convection efficiency, evaporation, and storage. It should be noted here that this method may slightly overestimate the vegetation cooling benefit for two reasons. First, the UHIs from the CTR simulations also include anthropogenic heat contribution. According to the offline attribution results, the anthropogenic heat contribution is about 20 % (Fig. 1). Second, street vegetation is inherently not identical to its rural counterparts due to the radiative trapping by urban canyons via the “canyon effect” (Wang, 2014; Ryu et al., 2016). Therefore, the same amount of vegetation inside urban canyons may not induce as much cooling as it would in the adjacent rural land.
We calculated the cooling benefit of green roofs using the diagnostic
surface energy balance data as
The Bowen ratio of green roofs was the average Bowen ratio value for
grassland calculated at the subgrid level by the CLM for the grass plant
functional type with its own soil column (Schultz et al., 2016). This
implies that the water-conserving native grass is used for green roofs. The
citywide Bowen ratio change was the area-weighted average of the green roof
Bowen ratio and the Bowen ratio of other surface components in the urban
land unit. The Bowen ratio change (
The green roofs were assigned the average warm-season grassland albedo value
of 0.20 (Bonan, 2008). The citywide albedo was the area-weighted average of
the green roof albedo and the albedo of other surface components in the
urban land unit. The albedo change (
The change in surface roughness induced by green roofs is minimal and thus was omitted in this analysis.
The same offline method was used to calculate the effect of reflective
pavement as
An example of the attribution result is given in Fig. 1 for the period
2071–2100 under RCP4.5. The offline
This linear supposition property of different biophysical contributions is
the theoretical basis underpinning the wedge approach that we use to
quantify the overall effectiveness of multiple UHI mitigation strategies.
Different mitigation strategies affect the urban surface temperature through
these biophysical contributions. Combining different UHI abatement actions
is actually adding up the perturbed biophysical contributions in the process
level. Linear supposition of various biophysical contributions states that
different mitigation strategies add up linearly as well. Therefore, good
agreement between online
The validity of our offline method is further supported by the excellent
agreement between the online and offline
Scatter plot between the daytime UHI intensity computed online with the CLM and that calculated offline using the surface energy balance diagnostic data. Filled and open symbols denote results from control (CTR) and white roof (WHT) simulations, respectively. Red, green, and blue colors denote dry (dry), continental (cont) and temperate (temp) climate regions, respectively. Circle, triangle, and square symbols denote current climate (cur), RCP4.5 (RCP45), and RCP8.5 (RCP85) scenarios, respectively.
We find strong urban climate change signals under both RCP scenarios. Near
the end of this century, the average summer daytime
In the current climate, the spatial variations of daytime
Urban heat island intensity from the control (CTR;
The mean anthropogenic heat flux for all the selected
cities and the cities in the three climate zones in the current climate, RCP
4.5, and RCP 8.5 scenarios from both the control (CTR) and the white roof (WHT)
simulations. (Units: W m
Adoption of the bright, reflective Solareflect roof material transforms
almost all the cities into cold islands in both the current and the future
climates, as indicated by the negative daytime
Maps of summer mean urban heat island intensity during
2071–2100 under the RCP4.5 scenario.
Our results indicate a stronger cooling effect of reflective roofs than those estimated by previous research. Studies using the WRF model showed a cooling effect of 0.2–2.0 K in screen-height air temperature with reflective roof materials of different albedos (Georgescu et al., 2014; Li et al., 2014; Stone et al., 2014). A case study of the California South Coast Air Basin area using a mesoscale model coupled with an urban airshed model showed that a citywide albedo increase by 0.30 can reduce the air temperature by up to 4.5 K in the mid-afternoon (Taha et al., 1997). However, Rosenzweig et al. (2009) reported very modest cooling effects (0.3–0.6 K) from cool roofs with albedo of 0.5 in New York City using a regional climate model. Studies using global climate models also report local cooling effects, ranging from a global average of 0.02 K (Jacobson and Ten Hoeve, 2012) to 0.4 K (Oleson et al., 2010c) decrease in air temperature with different prescribed cool roof albedos. Apparently, these estimates vary with the choice of cool roof albedo as well as the roof space available for cool roof installation (Akbari et al., 2012; Rosenzweig et al., 2009).
There are three reasons why we observe a stronger cooling effect compared to
previous studies. First, the albedo of the cool roof in our study (0.88) is
at the higher end of the range (0.3–0.9) used in previous modeling
studies. Second, we define
Our results are based on the albedo value of roof coating material after 3 years of use. The long-term weathering and aging problem of the reflective material has not been taken into account. Although the albedo of the Solareflect coating material decreases only slightly from 0.89 to 0.88 after 3 years, its long-term performance of retaining solar reflectance is unknown. In addition, none of the current roofing materials can last 96 years (2005–2100) while preserving high albedos. Re-coating or reroofing seems necessary for long-term use. Therefore, our results should be interpreted as an upper bound of the cooling potential of the cool roof strategy.
Street vegetation is less effective than cool roofs in terms of temperature
reduction. The daytime surface temperature reduction is 1.3
As mentioned previously, our method for calculating the cooling potential of street vegetation may slightly overestimate the cooling benefit because of the anthropogenic heat contribution to the UHI and the radiative trapping effects by urban canyons. Our offline attribution results show that anthropogenic heat contribution is about 20 % (Fig. 1). Ryu et al. (2016) showed that the trees in the urban canyon decreased the sensible heat flux and increased the latent heat flux by roughly a similar amount. Another potential source of uncertainty to this method is the street vegetation composition. Cities that have different vegetation composition from their surrounding rural landscapes may generate slightly different cooling than predicted here. However, integrating all these considerations into our results would not change the point that street vegetation reduces urban surface temperature less effectively in comparison to cool roofs.
The cooling effects of street vegetation estimated in this study are generally consistent with previous studies. Rosenzweig et al. (2009) found that planting trees in open spaces and along streets reduces the mid-afternoon air temperature by 0.2–0.6 K in neighborhoods of New York City. Luley and Bond (2002) reported that in their maximum scenario, in which all urban grass is replaced by trees, an up to 1.0 K reduction in air temperature is expected in Manhattan, New York City, on a summer afternoon. By forcing a mesoscale weather forecast model with atmospheric conditions generated by a global climate model, Stone et al. (2014) found that, under the business-as-usual climate change scenario, increasing green vegetation in urban public areas reduces the urban daily average air temperature by 0.1–0.3 K. Consistent with our results, they also reported a weaker vegetation cooling effect over a city in the dry climate compared to the cities in the temperate climate. In extreme scenarios, in which the urban center of Atlanta, Georgia, is fully substituted by forest, Stone et al. (2013) found a reduction in air temperature of 0.5–1.5 K. Generally, the cooling benefit varies with different cities because of the area available for tree planting. Our estimates fall at the higher end of the range reported previously, partly because the temperature difference is measured by surface radiative temperature.
Green roofs generate modest cooling. The average daytime surface temperature
reduction is 1.6
These results are comparable with the green roof cooling effects reported by previous studies. Building-scale observational studies have shown a relatively large range of temperature reduction. A case study evaluating the “green” and “white” policy in Chicago, Illinois, found that green roofs can reduce the building-top temperature by 0.3–2.6 K according to the LANDSAT satellite observations (Mackey et al., 2012). Ismail et al. (2011) used an experimental approach on a test building in Malaysia and showed a 4.6 K reduction in the surface temperature around noontime between a green roof and a black bare roof. Mesoscale modeling studies have also shown cooling estimates ranging from 0.2 to 0.8 K if measured by air temperature (Rosenzweig et al., 2009; Stone et al., 2014) and 3 to 4 K if measured by surface temperature (Li et al., 2014). All these results, together with our study, confirm that a major factor in determining the mitigation potential of green roof is water availability.
Our results favor cool roofs as the preferred method for urban heat
mitigation. Although empirical data are not yet available to validate the
white oasis effect, some lessons can be drawn from studies of green oases.
Local cooling has been observed for green oases as small as several hectares
in area (Kai et al., 1997). A typical green oasis effect ranges from
An explanation to this phenomenon can be achieved using the big-leaf model
(Sellers et al., 1996) solution of surface temperature shown below:
On the global (Pacala and Socolow, 2004) and urban (Creutzig et al., 2015) scales, the highest carbon mitigation potential is achieved by using the strategy of mitigation wedges or combining the incremental benefits of a number of abatement actions. Because nonlinear interactions among the various biophysical UHI contributors are small (Zhao et al., 2014; Fig. 1), a similar strategy can be used to estimate the aggregated potential of multiple urban temperature reduction methods. Figure 5 shows the additive benefit of three UHI mitigation wedges: cool roofs (or green roofs), street vegetation, and reflective pavement. Because cool roofs and green roofs are mutually exclusive under the assumption of a 100 % penetration rate, we did not combine the cooling benefit of cool roofs and green roofs in a single scenario. The urban daytime surface temperature is on average 2.4 K greater than the rural background temperature in the current climate. Near the end of the century, the urban temperature will be 7.3 K greater under the RCP4.5 scenario than the current background temperature due to the GHG warming (4.2 K) and the UHIs (3.1 K). The combined wedges with cool roofs (Fig. 5a) provide greater cooling benefit than the ones with green roofs (Fig. 5b). The total cooling effect of the three mitigation wedges with cool roofs is 8.0 K, essentially eliminating all the UHI effect plus the GHG warming.
A UHI strategy consisting of three mitigation wedges
under the RCP4.5 scenario.
The example given in Fig. 5 is the most optimistic scenario. Because roof
space is a “precious” resource that must accommodate other competing
needs, a 100 % conversion to white roofs is probably impractical. However,
we can linearly scale the endmember results (Table S1) to
estimate the cooling benefits of other wedge combinations. For example, a
wedge strategy consisting of 50 % cool roof, 50 % green roof, 100 %
street vegetation, and 100 % reflective pavement will bring a total
temperature reduction of 5.7 K (cool roof 3.3 K
Conversion of roof space to solar photovoltaics (PVs) can also cool the urban
air (Hu et al., 2016; Salamanca et al., 2016). A typical solar PV has an
albedo of 0.10 and a futuristic electricity conversion efficiency of 25 %
(Hu et al., 2016), giving a heat removal rate equivalent to having a surface
albedo of 0.33, which is 0.04 more than the average default roof albedo
value. Scaling the WHT simulation result proportionally by the albedo
increase, we estimate that the solar PV will lower the surface temperature
by an average amount of 1.2 K in the daytime. This estimate is made using
the offline energy balance diagnostics and Eq. (6), where
The above discussion focuses on mitigation of the daytime temperature.
Unfortunately, none of the strategies are effective in eliminating the
nighttime UHI. Conversion to cool roofs reduces building storage of solar
radiation in the daytime and subsequent heat release at night, contributing
to nighttime cooling by an average amount of 0.7
Comparison between combined cooling benefit (combined
simulation results) and sum of the components (sum of individual simulation
results). Both are measured by screen-height air temperature. Black symbols: cool roof
The above wedge strategy is applied to the surface UHI. To find out if the
linear supposition can also be extended to air UHI, we also conducted a
meta-analysis based on the published studies that presented modeled air UHI
results of both combined mitigation strategies and individual strategies. We
calculated the sum of the cooling benefit from each simulation of individual
strategy (as is shown in the
The UHI intensity is projected to increase in future warmer climates, partly due to higher anthropogenic heat release from AC energy use to cope with the GHG-induced warming. Our modeling analyses favor cool roofs as the preferred method for urban heat mitigation in comparison to green roofs, street vegetation, and reflective pavement. By adopting highly reflective roofs citywide, almost all the selected cities in the USA and in southern Canada are transformed into white oases. Cool roofs also bring large daytime AC energy savings in future climate scenarios. A UHI mitigation wedge strategy consisting of 50 % cool roof, 50 % green roof, street vegetation, and reflective pavement has the potential to reduce the urban daytime surface temperature by 5.7 K in the summer from the unmitigated urban scenario. Unfortunately, none of the UHI mitigation methods are effective in eliminating the nighttime UHI.
Cities are in fact engineered landscapes. The above UHI wedge strategies amount to a re-engineering of these landscapes (Animation Movies S1 and S2). A key distinction between urban engineering and geoengineering is the scale. Unlike planetary-scale geoengineering, urban engineering impacts a much smaller areal extent (about 2.7 % of the terrestrial land; Schneider et al., 2009). Planetary albedo modification can temporarily cool the global climate but at potentially large environmental costs (Mcnutt et al., 2015). In contrast, re-engineering of the urban land should have minimal negative consequences. Reorientation of the discussion from the global-scale albedo intervention to the small-scale temperature modification can motivate local actions because the payback is immediate and direct.
Data used in this study are CESM simulation results conducted by Lei Zhao on the US National Center for Atmospheric Research's cluster “Yellowstone” and can be provided upon request to the corresponding author, Lei Zhao (lei.zhao@princeton.edu).
The authors declare that they have no conflict of interest.
This research was supported by the Ministry of Education of China (grant PCSIRT), a Princeton STEP Postdoctoral Fellowship (to Lei Zhao), and a Yale University Graduate Fellowship (to Natalie Schultz). We thank Chang Cao for making the animation movies. We acknowledge high-performance computing support from Yellowstone (ark:/85065/d7wd3xhc) provided by NCAR's Computational and Information Systems Laboratory, sponsored by the US National Science Foundation. Edited by: Yun Qian Reviewed by: three anonymous referees