Urban warming and future air-conditioning use in an Asian megacity: importance of positive feedback

The impact of feedback between urban warming and air-conditioning (AC) use on temperatures in future urban climates is explored in this study. Pseudo-global warming projections are dynamically downscaled to 1 km using a regional climate model (RCM) coupled to urban canopy and building energy models for current and six future global warming (ΔTGW) climates based on IPCC RCP8.5. Anthropogenic heat emissions from AC use is projected to increase almost linearly with ΔTGW, causing additional urban warming. This feedback on urban warming reaches 20% of ΔTGW in residential areas. This further uncertainty in future projections is comparable in size to that associated with: a selection of emission scenarios, RCMs, and urban planning scenarios. Thus this feedback should not be neglected in future urban climate projections, especially in hot cities with large AC use. The impact of the feedback during the July 2018 Japanese heat waves is calculated to be 0.11 °C.


INTRODUCTION
With the global proportion of people living in cities expected to exceed two-thirds by 2050, 1 the urban climate affects many aspects of human life, including energy demand, health and the economy.Therefore, urban climate projections that include extreme weather events, such as heat waves from climate change, are of interest, given the potentially wide range of social and scientific impacts on human activities.These projections are also important when producing strategies to mitigate urban heat and for adapting to urban climate change.For example, given that US urban temperatures are predicted to increase by 1-2 °C due to climate change, energy demand is predicted to increase by 5-25%.Therefore, effective adaptation strategies that include measures to cool the urban climate and built environment that do not result in further emissions of heat or greenhouse gases (GHGs) will need to be produced and assessed. 2f the two main methods to obtain fine horizontal resolution projections of climate, the dynamical (rather than statistical) downscaling is used here.Several previous dynamical downscaling studies of future urban climate have used regional climate models (RCMs) with a coupled urban canopy model (UCM) (hereafter RCM/UCM).For example, Kusaka et al. 3 predicted the 2070s summertime climate for three Japanese megacities using the Weather Research and Forecasting (WRF) model 4 with the single-layer UCM, 5,6 assuming pseudo-global warming (PGW), 7,8 at a 3-km resolution.This future climate projection (dynamical downscaling) method uses reanalysis data with the climatic change component from global climate models (GCMs; see 'Methods' section).Similar methods have been used for cities in Asia, [9][10][11] Europe, 12,13 North America 14,15 and Oceania. 16In these RCM/UCM studies, anthropogenic heat (Q F ) is an important term in the urban surface heat balance but is assumed to have a constant diurnal pattern (e.g.Fig. 3 of Kusaka et al. 17 ).Thus Q F impact on future urban temperatures is kept the same as current climate if all other factors (e.g.urban structure and human activities) do not change.This assumption may be reasonable in cities with little air-conditioning (AC) use (e.g.northern European cities); however, where AC use is already common (e.g.Asian cities), it may not.
If cities experience positive feedback from interactions between urban warming and AC use, [18][19][20][21] this assumption would not hold.When air temperatures increase, energy consumption associated with AC use surges (Fig. 1a, Path 1) (i.e. increase in Q F from AC use (Q F, AC )).In turn, the energy released outdoors will enhance urban warming (Fig. 1a, Path 2). 22Using a simple one-dimensional mixing-layer model, Ohashi et al. 23 estimated the increase in summertime air temperature from AC waste heat in Osaka City to be +0.36 to +0.72 °C.Previous studies focus on only one side of this process.9][30][31][32][33] However, the amount of 'additional warming' from urban warming-AC feedback has yet to be assessed.Hence, the impact of this feedback process on future urban climate remains unknown.
In this study, we define the temperature difference caused by AC use as δT AC (Fig. 1b, orange): where δT AC ≠ FB is the temperature difference caused by AC use but without AC feedback (Fig. 1b purple, without Path 2 Fig. 1a [no-Q F,AC simulation], see 'Methods' [Model settings]).The 'additional temperature difference' caused by the feedback process is δT AC → FB (Fig. 1b, green, with Path 2 Fig. 1a).In this study, we use model simulations to determine the different components (Fig. 1b and 'Methods' provide more details).The Salamanca et al.'s 22 WRF/BEP+BEM (building effect parameterisation and building energy model) [34][35][36] simulations, with and without Q F,AC release, allows δT AC to be estimated (Fig. 1b, reference (R) Salamanca et al., 22 R S14 ).The 10-day Phoenix heat wave study obtained an extreme δT AC of 1.0 °C at night.However, as they did not separate δT AC , δT AC → FB remains unknown.Given that δT AC → FB may enhance electricity consumption (hereafter δEC AC → FB ), we do not have a basis for understanding these changes with climate change (Fig. 1b, red).
With increases in both global temperatures 37 and AC demand, 38 it is important to explore this positive feedback phenomenon and evaluate its impact on the urban climate.Projected increases in AC demand and urban expansion may increase the feedback effect enhancing future projections uncertainties, whereas improved AC performance may offset demand from a warmer ambient environment.These could be as important as other uncertainties, such as the selection of emission scenarios, RCMs, urban planning scenarios and so on.Thus the feedback process consequently affects the impact assessment results and policy decisions associated with the Paris Agreement. 39irect feedbacks can be calculated from climate projections for cities with AC usage when the RCM/UCM is coupled to a BEM.Electricity demand and summer thermal comfort projections for Nagoya in the 2070s 40 used WRF with coupled multi-layer UCM (CM) 41 and BEM 19 (WRF/CM+BEM 42 ) (Fig. 1b, R T15 ).Impacts of urban expansion and global warming on future urban temperatures and cooling demand for the Phoenix and Tucson metropolitan areas 43 explored with WRF/BEP+BEM (Fig. 1b, R T17 ) did not assess the 'additional warming' from positive feedbacks (δT AC → FB ) (Fig. 1b, green) but did determine the amount of urban warming from current and future climate (ΔT) and the future electricity demand.This study addresses: (1) What is the magnitude of the impact (and associated uncertainty) associated with such a feedback (δT AC → FB and δEC AC → FB ) on future urban climate (Fig. 1b, green)?(2) How is δT AC → FB and δEC AC → FB changed by climate change (i.e.increase/decrease, linear/nonlinear) (Fig. 1b, red)?
To address these, objectives WRF/BEP+BEM is used ('Methods').From the results, we propose a simple parameterisation to account for δT AC → FB and δEC AC → FB in urban climate studies.We focus on Japan's second-largest megacity, Osaka (Fig. 2, Supplementary Fig. 2), as it experiences the warmest summertime mean temperatures in Japan. 44Osaka's humid climate (normal annual precipitation (1981-2010) is 1279 mm) results in a more significant daytime urban heat island intensity than other dry climate cities. 45

RESULTS
Changes in Q F,AC and urban air temperature (ΔT) The numerical model (see 'Methods') is verified (see Supplementary Note) for wooden detached dwellings, fireproof apartments and commercial and office buildings (hereafter simply office).Given the similarity of the results between the two dwelling types, only the fireproof apartments (hereafter simply residential) results are presented.
Fig. 1 Processes simulated through numerical experiments.a The interaction between the indoor energy use by air-conditioning (AC) can feedback (FB) on the outdoor environment (e.g.enhancing the temperature T, Path 2), which in turn enhances AC use (Path 1).b To evaluate these impacts, a series of simulations are undertaken: (i) control case (AC → FB) and (ii) no-Q F, AC case (AC ≠ FB), for current and future climates.From analysing the simulations, the trends (arrows) caused by global warming (ΔT GW ) (grey), urban warming calculated by AC → FB (ΔT AC → FB ) (red) and AC ≠ FB (ΔT AC ≠ FB ) (blue) and urban warming if Q F,AC or Q F in the future are the same as in the current climate (ΔT const.QF ) (black) are shown.Temperature differences are caused by AC use (orange), anthropogenic heat emitted by AC use without feedback (purple) and 'additional' impact of Q F,AC on temperature (δT AC → FB ), which are impacted by the 'additional' urban warming difference between the AC → FB and AC ≠ FB (ΔT AC → FB and ΔT AC ≠ FB ) (green).Climate simulations were undertaken for eleven August periods for current and future climates by AC → FB (red circles) and AC ≠ FB (blue circles).Approaches adopted by R S14 , 22 R T15 , 40 R T17 43 and R K12 3 (asterisks) are indicated.Inset shows feedback process caused by the interaction between outdoor weather change (urban warming due to climate change) and air-conditioning (AC) use.Here, Δ (e.g.ΔT) is used to distinguish the difference between current and future climate and between cases with the same ΔT GW , δ is used (e.g.δT is temperature difference between AC → FB and AC ≠ FB at current climate or a future climate) The mean Q F,AC at 14:00 local time (LT) (i.e.close to a daily maximum of air temperature) in August is larger for office than residential areas for all seven control simulations (hereafter AC → FB) (Fig. 2a-g and 'Methods').Q F,AC increases from current to future climates (ΔQ F,AC ) as ΔT GW increases, which is associated with additional AC use (Fig. 2a-g) causing a close to linear increase in both the office (ΔQ F,AC /ΔT GW = 1.76 W m −2 °C−1 ) and residential areas (3.32 W m -2 °C-1 ).When ΔT GW is +3.0 °C, ΔQ F,AC is 1.37 (office) to 1.81 (residential) times larger than for current climate conditions (Supplementary Table 1).
The mean August 2 m air temperature (24 h, 11 years) in Osaka is warmer than the surrounding areas in both current and future climates (Fig. 2h-n).The temperature difference between the AC → FB and the no-Q F, AC case (AC ≠ FB) simulations increases along with ΔT GW (Fig. 3), indicating that urban heat will increase over time.This is caused by ΔQ F,AC with larger differences in the residential than in the office areas.
The differences between simulations with (ΔT AC → FB ) and without (ΔT AC ≠ FB ) anthropogenic heat fluxes allows the δT AC → FB to be determined (Figs.1a 2o-t and 3).The δT AC → FB (05:00, 11 August months) and ΔT GW increase and the δT AC → FB in Osaka's surrounding areas are higher than those of the centre of Osaka in both current and future climates (Fig. 2o-t).The diurnal variation of the δT AC → FB in residential areas (Fig. 3b) has two main features: (1) δT AC → FB for current and future climates are relatively small during daytime (notably 09:00-18:00) but larger at night (notably 18:00-06:00), with a peak around 05:00; and (2) δT AC → FB increases gradually as ΔT GW increases (+0.5 °C to +3.0 °C) during the night, indicating that urban warmth will increase overnight.
The relatively small daytime δT AC → FB , despite the relatively large daytime Q F,AC , is because of the relatively high daytime mixed layer.As Q F,AC is mixed in a large volume, the impact of Q F, AC on surface air temperature is reduced.At night, Q F,AC is smaller, but the mixed layer height is lower.Similar results in previous work 31 support this proposed mechanism.
The relation between ΔT GW and downscaled urban warming (ΔT AC → FB and ΔT AC ≠ FB ) has two main features for office areas (Fig. 3d): (1) δT AC → FB tends to increase linearly with ΔT GW .This increase is from a feedback: ΔT GW modifies ΔQ F,AC , contributing to additional urban warming (Fig. 1a, Path 1 → Path 2) This is the first study to estimate the 'additional' positive feedback impact on future climate associated with AC.We conclude that ΔT AC → FB has an approximately linear trend (Fig. 3c-f, red).The resulting AC → FB case slope (ΔT AC → FB /ΔT GW ) is 1.18 °C °C-1 for the office areas.This is larger than for the AC ≠ FB case (1.13 °C °C-1 ).These results suggest that it is relatively easy to estimate (parameterise) the impact of the feedback.
The daily average of normalised δT AC → FB (see 'Methods') is roughly constant (3-5%) and not dependent on ΔT GW for the office area (Supplementary Fig. 1a).The average of normalised δT AC → FB at 05:00 increases from 3% (+0.5 °C) to 10% (+3.0 °C) as ΔT GW increases (Supplementary Fig. 1a, Supplementary Table 1).Results for the residential area are similar to the office area, but feature (1) is more evident, especially at 05:00 (Fig. 3e, f).The ΔT AC → FB /ΔT GW (1.38 °C °C-1 ) is larger than the AC ≠ FB case (1.16 °C °C-1 ).As AC is used 24 h per day in residential areas, this cumulatively contributes more to warming than in office areas.The daily average of normalised δT AC → FB increases slightly from +0.5 to +3.0 °C, reaching 8% (at +3.0 °C).Thus about 0.25 °C additional urban warming is caused by the feedback when ΔT GW is +3.0 °C in the residential area.The normalised δT AC → FB at 05:00 increases from 8% (at +0.5 °C) to 20% (at +3.0 °C).This means that about 0.6 °C additional urban warming is caused by the feedback when ΔT GW is +3.0 °C (Fig. 3e).

DISCUSSION
Importance of the feedback in adapting to future urban climate change The choice of emission scenario causes uncertainty in projections.For example, the 2070s August mean urban temperatures in Osaka in the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios with the A2 scenario projected 0.3 °C higher than A1b. 46Doan and Kusaka's 11 urban warming projections (current to 2050s) for greater Ho Chi Minh City based on IPCC representative concentration pathway 4.5 (RCP4.5)and RCP8.5 scenarios found different scenarios to cause urban warming differences of 0.5 °C.The feedback uncertainty (~0.25 °C for residential areas when ΔT GW = + 3.0 °C) is comparable to these emission scenario differences of 0.3 and 0.5 °C, respectively (Table 1a).Further uncertainty arises from the RCM choice.Kusaka et al.'s 10 August 2050s urban temperature projections in central Tokyo, using WRF and NHRCM (nonhydrostatic model of the Japan Meteorological Agency 47 ), had a 10-year mean temperature increase uncertainty between RCMs of ~0.2 °C (i.e.like the feedback reported here, see Table 1a).
Our work shows δT AC → FB could reach 0.25 °C in residential areas.This is comparable to the 0.3-0.5 °C century −1 estimated global warming in Japan. 48This feedback impact of 0.25 °C in residential areas is comparable to urban planning scenarios (Adachi et al., 49 Kusaka et al. 10 ) for 'dispersed' (~0.34 and 0.1 °C, respectively) and 'compact' (~0.1 and 0.4 °C, respectively) cities and future urbanisation on temperature of 0.5 °C11 (Table 1a).
The feedback impact is small in the daytime but large nocturnally, increasing with ΔT GW (Fig. 3).The normalised δT AC → FB at 05:00 increases with ΔT GW : 10% in the office areas and 20% in the residential areas, when ΔT GW is +3.0 °C (Supplementary Fig. 1).These results suggest that mitigating the feedback process before the impact becomes large may be an effective strategy to adapt to climate change in cities as the feedback effect could be comparable to those of past global warming, affecting urban planning scenarios.For example (Supplementary Discussion), the feedback uncertainty for a ΔT GW of +2.0 °C is ~10 years (Fig. 4b, Table 1b), if we can stop the feedback technologically (e.g.improved coefficient of performance (COP), geothermal energy use) we could postpone a +2.0 °C world by about 10 years.
Application: δT AC → FB and δEC AC → FB estimates for the recent heat waves in Japan As shown ('Results'), ΔT AC → FB and δT AC → FB have near linear trends with ΔT GW .Also δT AC → FB and δEC AC → FB are close to linear with ΔT AC → FB .This suggests that to estimate (parameterise) δT AC → FB and δEC AC → FB for other future urban climates or for specific events (e.g.heat waves) is straightforward.Table 2 shows gradients (a) of linear regression (y = ax) between ΔT AC → FB (x) and δT AC → FB (δEC AC → FB ) (y).Here we estimate the δT AC → FB and δEC AC → FB for the recent heat waves of July 2018 and August 2013 in Japan using our proposed linear relations to illustrate an application.
The 2013 summer was the warmest on record for Japan (statistics since 1946 50 ), with 41.0 °C recorded in western Japan (Shimanto City).In Osaka City, the August monthly mean (30.0 °C) was 0.99 °C higher than the 11-year mean (2000-2010).July 2018 was one of the hottest Julys in Japan since 1946, with Osaka's monthly mean temperature (29.5 °C) 1.63 °C higher than the 11year mean (2000-2010).Therefore, August 2013 and July 2018 roughly correspond to situations when ΔT AC → FB = +1.0°C and +1.5 °C, respectively.Figure 3a, b shows a diurnal variation of the δT AC → FB for the two heat waves estimated by linear relations (Table 2a).The δT AC → FB for the two heat waves during the night to morning are higher than that for the daytime, especially in the residential areas.The 24-h mean δT AC → FB for July 2018 is 0.11 °C in the residential areas, which is higher than the 0.07 °C of August 2013.Similarly, the 24-h mean δEC AC → FB for July 2018 is 0.05 W floor-m −2 , which is higher than 0.03 W floor-m −2 of August 2013.Estimates of δT AC → FB and δEC AC → FB could be easily assessed for other heat waves and/or future urban climates in hot cities where significant AC use is common.

Summary
Here a positive feedback from the interaction between urban warming and AC use is identified and quantified.Analysis of simulations for current and six future climate scenarios (global warming: ΔT GW ) is undertaken.For the latter, CMIP5 GCMs simulations with the highest IPCC GHG emissions scenario (RCP8.5)were used.The megacity of Osaka is analysed for August when AC use is at its greatest.From this, it is concluded that: (i) Anthropogenic heat emissions from AC use (Q F,AC ) are predicted to increase linearly with ΔT GW from current to future climates.Monthly total Q F,AC in commercial & office and residential areas are projected to be 1.37

METHODS
To distinguish the difference between current and future climate Δ (e.g.ΔT) is used, and between cases with the same ΔT GW , δ is used (e.g.δT is a temperature difference between AC → FB and AC ≠ FB (explained later) at current climate or a future climate).First, the impacts of the feedback (δT AC → FB ) on future urban climate from global temperature scenarios (hereafter ΔT GW ) are evaluated.Second, we explore how δT AC → FB changes in relation to ΔT GW .The urban air temperature and EC calculated by WRF/BEP+BEM were verified using detailed observational data for a year. 51scription of the numerical model Previously, WRF/BEP+BEM simulated urban air temperature and EC for April 2013 to March 2014 in Osaka were verified. 51The modified model reproduced current diurnal and horizontal variations of 2 m urban air temperature and EC at a 1-h temporal resolution for 12 electric power substations within Osaka City (about 1 km 2 area).
At each time step, Q F,AC is calculated for each grid as: 19,51 where H out and E out are the sensible and latent heat supplied from the AC system for cooling (see section 2.5 of Salamanca et al. 35 ), respectively, and COP is the coefficient of performance.The H out and E out are calculated from the total sensible, H in , and latent, E in , heat loads per floor. 35The H in considers four physical processes: (1) solar radiation through the window and the heat exchange between the windows and the indoor air, (2) heat conduction through the walls and heat exchange between the wall, ceiling and pavement and the indoor air, (3) sensible heat exchange through ventilation, and (4) the internal sensible heat generation from equipment and occupants (Fig. 1a).Outdoor weather directly affects the H in through processes (1-3), and the affected H in contributes to an increase in Q F,AC and EC through the H out (Fig. 1a Path 1).Here Q F,AC is split into sensible heat Q F,AC,S and latent heat, Q F,AC,L , as commercial and office building areas are considered 52 The WRF/BEP+ BEM model system assumptions include: 51 (1) individual AC units exist (using a simple parameterisation); (2) a constant COP (implications discussed in Supplemental information) occurs; (3) no Q F from traffic occurs; and (4) only typical weekdays occur.Therefore, the results of this study are for these conditions (i.e.weekdays).
Although Q F from traffic (Q F,traffic ) is important 53,54 in some cities, in Osaka 42,52 Q F,AC is at least two times larger.As traffic is already common in Asian megacities, while AC demand is dramatically increasing 38 (e.g. in Vietnam penetration rates of AC, motorcar and motorbike are 10.8%, 77.4%, and 81.8%, respectively 55 ), Q F,AC has the potential to be a big problem in the future (cf.Q F,traffic ) but will benefit from COP increases.Future studies will consider Q F,AC and Q F,traffic using models such as WRF/ CM+BEM. 42However, the conclusions of this report are not impacted by Q F,traffic as it is assumed to be constant and uses the single-layer UCM (SLUCM) with a static Q F 5,6 profile (traffic and AC) giving parallel results to urban warming without Q F emissions (see next section 'Model settings').

Model settings
Here the Advanced Research WRF model (ver.3.5.1) 4 was used with the same model parameters (Supplementary Table 3) and physics as previously described. 515][36] As previously indicated, these models can accurately reproduce the diurnal variation and horizontal distribution of summertime surface air temperature and EC due to AC use. 51he model domain (Supplementary Fig. 2a) covers western Japan with 126 x and y grids in both d01 and d02 domains (horizontal resolution 5 and 1 km, respectively, and two-way nesting).As the reanalysis captures the summer synoptic-scale features (high-pressure system) around Japan, 63 a larger coarse grid is not needed.The model top is 50 hPa, with 35 vertical sigma levels.The vertical resolution close to the ground (WRF atmospheric first layer height) is nearly 50 m.Therefore, 10 building layers below 50 m (5 m resolution) in BEP/BEM is used, as in previous studies, 42,51 due to high model reproducibility for 2 m air temperature and EC in Osaka with the first layer at 50 m.With a mean building height of about 25 m in the office area, 50 m corresponds to the constant flux layer.The innermost domain urban grid classifications are based on the dominant building type, using land use and land cover and topographic data sets from the Geospatial Information Authority of Japan data and Osaka City GIS polygon (building footprint) data, with building use (type), building height and total floor area for each building (Supplementary Fig. 2c, d) grouped into: (i) commercial and business grids (i.e.mainly commercial and office buildings; hereafter C), (ii) residential grids with predominantly fireproof apartments (Rr), and (iii) residential grids with predominately wooden detached dwellings (Rw).The outer and inner domain urban grids outside Osaka City are classified as Rw.
For current climate simulation, initial and boundary conditions are derived from the National Center for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) reanalysis data 64 and merged satellite and in situ global daily sea surface temperature (MGDSST) data. 65For future climate simulations, we also use a modified the NCEP-NCAR and MGDSST as initial and boundary conditions (see details in the next subsection).As previously, 3 time integration is conducted from 00:00 UTC July 27 to September 1 (i.e.August with 5 days spin-up) for 11 years 2000-2010 to allow climatological discussion of the simulation results.Other studies integrate for only 1 month, several weeks or days (current and future climate), which is insufficient for climatological discussion in general.In Osaka, August is the hottest month, with more clear skies (cf.June, July and September, rainy season) and the highest EC recorded.These model settings are termed the control simulation (case AC → FB) (Fig. 1b, red arrow) for all climates, including current and six future climates (see next subsection).
The no-Q F, AC simulation (case AC ≠ FB) differs from case AC → FB as Q F, AC = 0 W m −2 (Fig. 1b, blue arrow); i.e. the [AC ≠ FB] − [AC → FB] difference is caused by Q F, AC and it causes the feedback (δT AC → FB ).Current climate conditions (August 2000-2010) and six future climates are simulated.Here δT AC → FB is estimated by ΔT AC → FB − ΔT AC ≠ FB (Fig. 1b).In the current climate case, δT AC → FB is 0 °C as we assume there is no longterm climate change (decades) (i.e.no increase in forcing temperature and ΔT AC → FB and ΔT AC ≠ FB are 0 °C).Of course, δT AC → FB may not be 0 °C in the real current climate, but it is difficult to estimate δT AC → FB from current climate only, and our objectives are to clarify δT AC → FB in relation to future urban climate and how δT AC → FB is changed by climate change (see 'Introduction').Therefore, here 'δT AC → FB driven by decades of climate change' is the focus and is estimated by a time slice experiment, including current climate simulation and six future projections.In a similar way, δEC AC → FB is estimated by ΔEC AC → FB − ΔEC AC ≠ FB .To determine δT AC → FB and δEC AC → FB , it is assumed all conditions (e.g.urban structures, human activities and AC system technology) remain constant, except for background climate change.Although AC technology will improve, it is important to know the maximum (potential) impact of the feedback driven by climate change on future urban air temperature.A constant COP scenario is used as the RCP8.5, and Shared Socioeconomic Pathway 3 scenario assumes slow technological change in the energy sector. 66,67e assume that δT AC ≠ FB is constant with ΔT GW ; i.e. black dashed arrow (ΔT const.QF ) parallel to the blue arrow (ΔT AC ≠ FB ) in Fig. 1b.To confirm this, the impact of Q F on surface air temperature (i.e.like δT AC ≠ FB ), on future climate, was conducted with simulations using SLUCM, which has a static Q F 5,6 (includes traffic and AC).Urban warming, when Q F in the future was the same as for the current climate (ΔT const.QF : Fig. 1b, black dashed arrow), paralleled urban warming with a no-Q F emitted case, like ΔT AC ≠ FB , although different urban parameterisation was used.A reason why we used SLUCM (not BEP/BEM) to confirm black dashed arrow (ΔT const.QF )'s gradient (ΔT const.QF /ΔT GW ) (and then compared with that of blue arrow (ΔT AC ≠ FB ) by BEP/BEM) is that it is quite difficult to calculate the black dashed arrow by the BEP/BEM, because Q F, AC at every time step (6 s) in future climates has to be completely same as that of current climate.To confirm the black arrow's gradient by SLUCM is technically much easier and essentially same as that by the BEP/BEM.

Climate projections
Six future climates, with background temperature increases (global warming: ΔT GW : +0.5, +1.0, +1.5, +2.0, +2.5, and +3.0 °C) relative to the current climate, are simulated.These are based on the ensemble mean results from four GCMs used in the Climate Model Intercomparison Project (CMIP5) 68 : CCSM4, 69 CESM1 (CAM5), 70 GFDL-CM3 71 , and INM-CM4. 72owever, the ensemble mean of only three GCMs (CCSM4, CESM1 and GFDL-CM3) is used when ΔT GW > +3.0 °C, as INM-CM4 does not reach +3.0 °C at 2100.The simulations considered the highest IPCC GHG emissions scenario (RCP8.5)(Fig. 4).CESM1 (CAM5) and CCSM4 have good model climate performance index scores 73 and summertime synoptic pressure patterns around Japan perform better than other GCMs.GFDL-CM3 and INM-CM4 are selected to cover minimum to maximum uncertainties in GCM projections.By using several GCMs, the uncertainties in the GCM projections (horizontal colour bars, Fig. 4) are incorporated in the RCM.A similar GCMs selection method for future urban climate projection was used by Adachi et al. 9 For future projection experiments, climate difference components between current and future are estimated by the four individual GCMs.For example, for ΔT GW = +0.5 °C, it was identified that the August mean surface air temperature difference, averaged spatially around Japan, became nearly +0.5 °C for each GCM.The other climate variables (e.g.geopotential height, horizontal wind and sea surface temperature) are extracted for each GCM for the identified years, and an ensemble mean is calculated for all variables.For ΔT GW = +0.5 °C, the ensemble mean surface air temperature difference component is +0.4725 °C (Supplementary Fig. 3) (i.e.nearly +0.5 °C).The actual ensemble mean surface air temperature differences for ΔT GW = +0.5, +1.0, +1.5, +2.0, +2.5 and +3.0 °C are +0.4725,+0.9550, +1.4750, +1.9650, +2.4450 and +2.8825 °C, respectively (Supplementary Fig. 3).For the six ΔT GW cases, the climate difference for each variable (i.e.wind components, geopotential height and temperature) are added to the NCEP-NCAR and MGDSST data (Supplementary Fig. 3), therefore changing the atmospheric variables such as long-wave radiation.][76][77][78] The PGW method 7,8 was used to create regional climate projections 3,10,11,21,[79][80][81] and has been previously verified. 80,82,83For example, Kawase et al. 82 projected East Asia's 1960s climate by the PGW method using 1990s reanalysis data and PGW increment, calculated by the climate difference between the 1990s and 1960s.They compared the 1960s projection results with a hindcast simulation for the same period using reanalysis data.They showed that the 1960s projection could reproduce some decadal changes in the rainfall shown within the hindcast simulation, thus indicating that the PGW method could reproduce actual energy balance change using reanalysis, as precipitation is calculated in response to the atmospheric energy balance.
Yoshikane et al. 83 compared PGW and direct downscaling from GCM for East Asian future climate projections and found no significant differences in temperature and precipitation (i.e.atmospheric energy balances were similar).Snowfall 79,81,84 research, a process strongly influenced by atmospheric energy balance near the surface, and torrential rain 85 studies, caused by atmospheric buoyancy near the surface, have also successfully used PGW methods.Therefore, the PGW method is reasonable for use in future projection research.An advantage of using the PGW method is a reduction in GCM climate bias as the method uses modified objective analysis/re-analysis data rather than GCM output. 10A disadvantage is that perturbations in meteorological variables are not considered. 86As Adachi et al. 87 show, the perturbation component is equally important as the climatology component for summer precipitation, hence the PGW method is inappropriate for this.However, the perturbations do not affect nearsurface air temperatures change (current to future climates) in western Japan, including Osaka during summer. 87Thus the difference in climatology projected by a GCM has a larger influence than the difference in perturbation on the downscaled temperature change and its variability in the climate projections. 87Hence, as the main contributor to the temperatures change was climate change (PGW) component, their results indicate that PGW method is appropriate for our purpose.As we focus on August, a dry period in Japan, surface air temperature is little impacted by precipitation.The PGW method is used as it enabled the study to focus on internal physical (feedback) processes in the RCM and to remove bias from individual GCMs.It also allows for direct comparison with previous PGW studies that report projected urban temperature uncertainties with selections of RCM, GHG emission scenarios and urban planning scenarios (see Table 1 and main text).

Normalised feedback impact
Three urban land uses are selected for analysis (Supplementary Fig. 2d): commercial and office buildings, fireproof apartments, and wooden detached dwellings.These types of areas were evaluated for surface air temperature and EC model performance. 48o address the question (1), a normalised δT AC → FB (feedback impact) per future climate is determined: where ΔT AC → FB and ΔT AC ≠ FB are urban warming calculated by the AC → FB and AC ≠ FB, respectively, as explained in Fig. 1.
All times referred to are LT.Japan does not use summertime.
Published in partnership with CECCR at King Abdulaziz University npj Climate and Atmospheric Science (2019) 39

Fig. 3
Fig. 3 Impact of the feedback process on air temperature (δT AC → FB ).a, b Diurnal variations of the δT AC → FB for current climate (grey line, i.e. 0 °C), six future climates simulated via ΔT GW and for heat waves in July 2018 (purple circles) and August 2013 (orange squares) in Japan estimated by linear relation between ΔT AC → FB and δT AC → FB for each time (Table 2a) for a office and b residential areas.c-f Relation between ΔT GW and ΔT (δT AC → FB , green triangles) calculated from the AC → FB (red) and AC ≠ FB (blue) simulations for c, d office and e, f residential areas.c, e August 05:00 mean and d, f 24 h mean air temperatures.Regression lines (dotted).Error bars indicate 25th and 75th percentiles

Table 1 .
ΔT GW is +3.0 °C.The δT AC → FB is expected to increase EC through a feedback process:ΔT GW → ΔQ F, AC → δT AC → FB → δEC AC → FB ('additional' EC increase from the feedback).The EC increase calculated by case AC → FB (ΔEC AC → FB ) is greater than that calculated by case AC ≠ FB (ΔEC AC ≠ FB ).The difference between these (ΔEC AC → FB − ΔEC AC ≠ FB ): δEC AC → FB increased linearly with ΔT GW ; Uncertainties in summer urban climate simulation projections: a Temperatures and b time (years) for the indicated warming to occur with Using linear relations proposed in this study (between ΔT AC → FB and δT AC → FB ), the δT AC → FB during recent Japanese heat waves (July 2018 and August 2013) are estimated to be 0.11 and 0.07 °C, respectively.Such an approach can be used to estimate δT AC → FB for other heat waves and/or future urban climates in hot cities where significant AC use is common.
GW → ΔQ F,AC (Q F,AC increase) → δT AC → FB ('additional' temperature increase from the feedback) Urban warming calculated when anthropogenic heat fluxes are permitted [case AC → FB (ΔT AC → FB )] is greater than when they are not included [case AC ≠ FB (ΔT AC ≠ FB )].The difference between these two [ΔT AC → FB − ΔT AC ≠ FB ] provides δT AC → FB .This increases almost linearly with ΔT *Calculated from three GCMs (see Supplementary Discussion) Fig.4Changes in August mean surface air temperature for Japan.a Projected increase compared to 2000s (decadal running averages) for four global climate models (CCSM4, CESM1 (CAM5), GFDL-CM3 and INM-CM4; for references, see 'Methods') and the average ΔT GW (black line) for four GCMs assuming the RCP8.5 scenario (see text for more details).b Time series of ΔT GW (thick black line) and ΔT AC → FB (circles) with maximum uncertainties (colour lines) of four GCMs and linear approximation of ΔT AC → FB (dashed black line).Values are averaged for three urban categories: commercial and office buildings, fireproof apartments, and wooden detached dwellings (iii) Uncertainty introduced by the feedback (δT AC → FB ) is comparable to that introduced from selection of emission scenario, RCM and/or urban planning scenario.Thus this feedback should not be neglected in future urban climate projection, especially in hot cities where significant AC use is common; and (iv)

Table 2 .
Gradients of linear regression (sensitivities) between a ΔT AC → FB and δT AC → FB (°C °C−1 ) by time of day and location and between b ΔT AC → FB and δEC AC → FB (W floor-m −2 °C−1 ) *Mean when AC in use Published in partnership with CECCR at King Abdulaziz University npj Climate and Atmospheric Science (2019) 39