A New Model to Downscale Urban and Rural Surface and Air Temperatures Evaluated in Shanghai, China

Dongwei Liu Shanghai Institute of Meteorological Science, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Dongwei Liu in
Current site
Google Scholar
PubMed
Close
,
C. S. B. Grimmond Department of Meteorology, University of Reading, Reading, United Kingdom

Search for other papers by C. S. B. Grimmond in
Current site
Google Scholar
PubMed
Close
,
Jianguo Tan Shanghai Climate Center, Shanghai Meteorological Service, and Key Open Laboratory of Cities’ Mitigation and Adaptation to Climate Change in Shanghai, China Meteorological Administration, and Shanghai Key Laboratory of Meteorology and Health, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Jianguo Tan in
Current site
Google Scholar
PubMed
Close
,
Xiangyu Ao Shanghai Institute of Meteorological Science, and Shanghai Key Laboratory of Meteorology and Health, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Xiangyu Ao in
Current site
Google Scholar
PubMed
Close
,
Jie Peng Shanghai Institute of Meteorological Science, and Shanghai Key Laboratory of Meteorology and Health, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Jie Peng in
Current site
Google Scholar
PubMed
Close
,
Linli Cui Shanghai Institute of Meteorological Science, and Shanghai Key Laboratory of Meteorology and Health, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Linli Cui in
Current site
Google Scholar
PubMed
Close
,
Bingxin Ma Shanghai Institute of Meteorological Science, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Bingxin Ma in
Current site
Google Scholar
PubMed
Close
,
Yan Hu Shanghai Institute of Meteorological Science, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Yan Hu in
Current site
Google Scholar
PubMed
Close
, and
Mingbin Du Shanghai Institute of Meteorological Science, Shanghai Meteorological Service, Shanghai, China

Search for other papers by Mingbin Du in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

A simple model, the Surface Temperature and Near-Surface Air Temperature (at 2 m) Model (TsT2m), is developed to downscale numerical model output (such as from ECMWF) to obtain higher-temporal- and higher-spatial-resolution surface and near-surface air temperature. It is evaluated in Shanghai, China. Surface temperature (Ts) and near-surface air temperature (Ta) submodels account for variations in land cover and their different thermal properties, resulting in spatial variations of surface and air temperature. The net all-wave radiation parameterization (NARP) scheme is used to compute net wave radiation for the surface temperature submodel, the objective hysteresis model (OHM) is used to calculate the net storage heat fluxes, and the surface temperature is obtained by the force-restore method. The near-surface air temperature submodel considers the horizontal and vertical energy changes for a column of well-mixed air above the surface. Modeled surface temperatures reproduce the general pattern of MODIS images well, while providing more detailed patterns of the surface urban heat island. However, the simulated surface temperatures capture the warmer urban land cover and are 10.3°C warmer on average than those derived from the coarser MODIS data. For other land-cover types, values are more similar. Downscaled, higher-temporal- and higher-spatial-resolution air temperatures are compared to observations at 110 automatic weather stations across Shanghai. After downscaling with TsT2m, the average forecast accuracy of near-surface air temperature is improved by about 20%. The scheme developed has considerable potential for prediction and mitigation of urban climate conditions, particularly for weather and climate services related to heat stress.

© 2018 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Prof. Jianguo Tan, jianguot@21cn.com, tanjg@smb.gov.cn

Abstract

A simple model, the Surface Temperature and Near-Surface Air Temperature (at 2 m) Model (TsT2m), is developed to downscale numerical model output (such as from ECMWF) to obtain higher-temporal- and higher-spatial-resolution surface and near-surface air temperature. It is evaluated in Shanghai, China. Surface temperature (Ts) and near-surface air temperature (Ta) submodels account for variations in land cover and their different thermal properties, resulting in spatial variations of surface and air temperature. The net all-wave radiation parameterization (NARP) scheme is used to compute net wave radiation for the surface temperature submodel, the objective hysteresis model (OHM) is used to calculate the net storage heat fluxes, and the surface temperature is obtained by the force-restore method. The near-surface air temperature submodel considers the horizontal and vertical energy changes for a column of well-mixed air above the surface. Modeled surface temperatures reproduce the general pattern of MODIS images well, while providing more detailed patterns of the surface urban heat island. However, the simulated surface temperatures capture the warmer urban land cover and are 10.3°C warmer on average than those derived from the coarser MODIS data. For other land-cover types, values are more similar. Downscaled, higher-temporal- and higher-spatial-resolution air temperatures are compared to observations at 110 automatic weather stations across Shanghai. After downscaling with TsT2m, the average forecast accuracy of near-surface air temperature is improved by about 20%. The scheme developed has considerable potential for prediction and mitigation of urban climate conditions, particularly for weather and climate services related to heat stress.

© 2018 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Prof. Jianguo Tan, jianguot@21cn.com, tanjg@smb.gov.cn

1. Introduction

Rapid urbanization globally has resulted in land-use and land-cover changes that have significantly affected climate (Jin et al. 2011; Bechtel et al. 2012). One of the most well-known urban climatic effects is the enhanced surface and air temperatures in cities, the urban heat island (UHI; Oke 1973; Kim and Baik 2005; Yow 2007). Many aspects of urban life are impacted (Oke 1995; Stone and Rodgers 2001; Weng et al. 2004) including human health and well-being (Basu and Samet 2002), the operations of city infrastructure (O’Malley et al. 2014), and ecological (Raich and Schlesinger 1992; Trumbore et al. 1996) and hydrological (Grimm et al. 2008; Kalma et al. 2008) processes. Thus, the ability to simulate high-resolution surface or near-ground air temperatures in cities is of considerable interest for a wide range of applications.

Satellite remote sensing has many advantages for monitoring land surface temperatures. Past studies have shown that the UHI derived from remote sensing responds to spatial patterns of surface characteristics (Roth et al. 1989; Gallo et al. 1993; Nichol et al. 2009). However, for many applications the value of remote sensing is limited by cloud cover, low temporal resolution of imagery, and/or coarse spatial resolution, although downscaling methods exist to address the latter (Xu et al. 2008; Bechtel et al. 2012; Bonafoni 2016).

Numerical simulations provide an effective way to forecast surface and/or air temperatures. Models are available for different surface types, including crops (Luo et al. 1992; Mihailović and Eitzinger 2007) and roads (Shao and Lister 1996; Diefenderfer et al. 2006). Methods underpinning these include obtaining the ground surface temperature by solving the surface energy balance equation (e.g., Best 1998; Herb et al. 2008); statistical downscaling from coarser-resolution models assuming stationary empirical relations (Winkler et al. 1997); multiobjective genetic-programming-based method to downscaling near-surface temperature by high-resolution information on land surface properties (Zerenner et al. 2016); and dynamic downscaling using a mesoscale model, for example, the Weather Research and Forecasting (WRF) Model (Pan et al. 2012).

An extensive comparison of 33 urban land surface models (ULSM) conducted as part of an international urban energy balance models comparison project (Grimmond et al. 2010) found, in general, simpler models perform just as well as more complex models (Grimmond et al. 2011; Best and Grimmond 2015). This finding is attributed, in part, to the ability of individual models to use the available input parameters and the difficulties of determining some of the parameters employed in more complex schemes. Moreover, complex models may require more computing resources (e.g., time, power). However, advances are being made with simple, practical models for cities (e.g., Wang et al. 2017; Lindberg et al. 2018).

This paper presents a novel simple model [the Surface Temperature and Near-Surface Air Temperature (at 2 m) Model (TsT2m)] for simulating high-spatial-resolution surface and near-surface air temperatures (2-m height). The model accounts for differences in surface cover proportion within a city (i.e., beyond coarse urban and rural differences) while using standard meteorological station data and easily obtained parameters. Critically, it runs rapidly and so permits operational high-resolution air temperature forecasts that allow heat- and cold-related urban climate and weather services to be provided. It considers both the surface energy budget and the effects of horizontal airflow. TsT2m is evaluated in Shanghai, a city with greater than 23 million inhabitants, more than 2.6 million automobiles, and more than 32 000 tall buildings (>30 m tall) in the year 2010 (Tan et al. 2015). Results are compared with 10-day weather forecast datasets derived from ECMWF model output (Molteni et al. 1996) and automatic weather station (AWS) observations distributed across the city of Shanghai.

2. TsT2m description

a. Model structure

TsT2m uses the surface energy balance to calculate surface temperature Ts (Ts submodel) and a dynamic energy balance of a column of air to diagnose the near-surface air temperature at 2-m height Ta (T2m submodel).

b. Surface temperature model (Ts submodel)

For the surface temperature model, the surface energy balance for each land-cover type is described as (Oke 1988)
e1
where Q* is net all-wave radiation, and QF, QH, QE, and ΔQS are the anthropogenic heat, turbulent sensible heat, turbulent latent heat, and net storage heat fluxes, respectively. The net all-wave radiation parameterization (NARP) scheme (Offerle et al. 2003; Loridan et al. 2011) is used to calculate Q*. NARP requires routine meteorological data (near-surface air temperature Ta, humidity, pressure, cloud cover, wind, etc.) and surface physical properties (surface albedo α, broadband surface emissivity εs, etc.):
e2
e3
where εsky is broadband atmospheric emissivity, and σ is the Stefan–Boltzmann constant (Offerle et al. 2003):
e4
The value of K depends on solar geometry and the solar constant R0:
e5
where h is the solar zenith angle. The atmospheric transmittance Tr is calculated as follows (Shao and Lister 1996):
e6
where fe_s is a factor to account for deviation in mean Earth–sun distance based on day of year D:
e7
The absorption by water vapor tw is calculated as follows (McDonald 1960):
e8
from the precipitable water μ that requires an empirical constant G dependent upon time of year and latitude and dewpoint temperature Tdew (Smith 1966):
e9
and optical mass m (Shao and Lister 1996):
e10
The absorption and scattering of aerosols ta is calculated as follows (Meyers and Dale 1983):
e11
The Rayleigh scattering and absorption by the permanent gases tRtpg additionally requires pressure p (Atwater and Brown 1974):
e12
Cloud transmittance tc uses cloud fraction (CCi, 0–1) and is the cloud transmittance by cloud level (ci, where i = low, medium, and high cloud) following Atwater and Ball (1981):
e13
If only total cloud-cover fraction (CC1) and/or lower/middle (CC2) cloud cover are available, this becomes
e14
where c1 is cloud transmittance without low cloud, and c2 is cloud transmittance without high cloud. In this study, to evaluate c1 and c2, a sensitivity test (at 0.2 steps in each parameter, from 0.0 to 1.0) is performed using observations of cloud cover (CC1, CC2) and K at the Xujiahui (XJH) site in Shanghai (Fig. 1). The value of K is calculated [Eqs. (4)(14)] when K > 0 W m−2 and cloud observations are available [December 2012–August 2013 at 0800 and 1400 local time (LT) (UTC + 8 h), 712 thirty-minute data points]. The goodness-of-fit metrics used are the root-mean-square error (RMSE), standard deviation, and correlation coefficient. The best estimate of K occurred when c1 = 0.4 and c2 = 0.8 (Fig. 2, point A). These values are used in the remainder of this study.
Fig. 1.
Fig. 1.

Land use (IGBP classes: numbers and color) and location of AWSs in the Shanghai region (province and constituent districts: gray lines) based on 2010 MODIS imagery (see text). The AWS national station codes are as follows: BS: Baoshan; CM: Chongming; FX: Fengxian; JD: Jiading; JS: Jinshan; MH: Minghang; SJ: Songjiang; PD: Pudong; XJH: Xujiahui; and QP: Qingpu. The eddy covariance (EC) tower is also at XJH, about 50 m from the AWS site. Simulated area is the same extent as this map (30.5°–32.0°N, 120.7°–122.2°E).

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

Fig. 2.
Fig. 2.

Normalized Taylor (2001) diagram for simulated incoming shortwave radiation assessed with observations from the XJH flux tower (Ao et al. 2016a) with changing combinations of transmission coefficients c1 and c2 (0.2 increment steps) for the period December 2012–August 2013 (30 min, at 0800 and 1400 LT, 712 samples). The ideal model performance is point O on the x axis. Points (red) are located based on the normalized standard deviation, correlation coefficient, and normalized RMSE. The best fit (point A) occurs with c1 = 0.4 and c2 = 0.8.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

The total input of energy in an urban environment (QF*=Q* + QF) requires the anthropogenic heat flux QF to be determined. Here, the Large-Scale Urban Consumption of Energy Model (LUCY) model (Allen et al. 2011; Gabey et al. 2018) with three sources of QF (traffic, building, and human metabolic emissions) is used. The traffic data are based on International Road Federation (IFR) World Road Statistics 2005, with the average traffic speed set to 48 km h−1. The population data are from the 2000 Global Rural Urban Mapping Project (GRUMP). For the average daily air temperature, the XJH station (September 2012–August 2013) data are used.

We obtain ΔQS from the objective hysteresis model (OHM; Grimmond et al. 1991):
e15
where fi is the proportion of surface type i; a1i, a2i, and a3i are the OHM coefficients (Table 1); n is the number of surface types; and t is the time. This is a robust, well-tested method to calculate the heat flux of stored energy in urban areas (e.g., Arnfield and Grimmond 1998; Taha 1999; Roth et al. 1989; Grimmond and Oke 1999; Roberts et al. 2006).
Table 1.

Parameters required for OHM and values used for each land-use class (IGBP codes; Fig. 2).

Table 1.
Bhumralkar (1975) and Blackadar (1976) developed a two-layer force-restore scheme to describe surface temperature variations. This is applied widely in land–air schemes (e.g., Xue et al. 1991; Mihailović 1991, 1996). Using ΔQS, we calculate surface temperature TS by the force-restore method over different urban land-cover types (e.g., building, pavement, vegetation, trees, water, barren):
e16
where Td is a deep (z ~ ∞) soil temperature at the substrate heat reservoir (Hu and Islam 1995). The value of Td is usually set as a constant, but for a limited period such as 1 day, the deep (z ~ ∞) soil temperature is not modified and a shallower slowly varying daily temperature is appropriate. Here, observed soil temperature at a depth of 0.2 m is used. Also, C is the volumetric heat capacity, KT is the thermal conductivity, and C0 is the heat capacity per unit area, determined from
e17
where ω = 2π/τ and τ is 1 day (24 × 3600 s).

c. Near-surface air temperature model (T2m submodel)

Given a column of well-mixed air of height Zi above a surface, the temperature change depends on the internal energy change ΔQa, including the horizontal energy change ΔQa,h, and the vertical energy change ΔQL_V:
e18
The vertical energy change of the column of air ΔQL_V is set to equal the changes caused by surface longwave radiation heat ΔL↑0 and longwave radiation change at the top of the air column :
e19
where A is the grid area and Δt is the time span.
e20
where Ts,t and Ts,t−1 are the surface temperature at time t and the preceding time (t − 1), respectively.
e21
The air temperature high (850 hPa) within the expected mixed layer , should reasonably represents the synoptic-scale changes and be impacted minimally by ground-based heat sources.
The energy brought by horizontal airflow ΔQa,h is a function of the total mass of a column of air in the grid (ma), the volumetric heat capacity of air Ca and the change of air temperature (ΔTh = TinTout) caused by horizontal advection:
e22
where
e23
where ρ is the density of air. From the internal energy change ΔQa, the air temperature change in the column ΔT within time span Δt is obtained:
e24
then
e25
Precipitation modifies a number of important processes, including increases to the water vapor in the boundary layer, which affects the transmittance of the atmosphere air density and heat capacity; changes to the albedo of the surface (Ao et al. 2016a), which affects the shortwave radiation flux; increases to the surface latent heat flux, which changes the Bowen ratio (i.e., QH and QE partitioning) (Ao et al. 2016b); and soil moisture, which impacts the soil thermal conductivity and volumetric heat capacity. Although all of these processes impact surface temperature and air temperature, the model does not consider these hydrologic processes. Its performance on rainy days should be considered in this context. The influences of rainfall/surface moisture on temperature are complex and need to be explored further in future model versions.

3. Data and methods

a. Model parameter setting

The detailed land-cover classification used by the surface temperature model is derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD12Q1 dataset (http://modis.gsfc.nasa.gov/data/dataprod/mod12.php) with 500-m resolution, based on the International Geosphere–Biosphere Programme (IGBP) land-cover classification (Cohen et al. 2006). In the study area (30.5°–32.0°N, 120.7°–122.2°E), the original 16 IGBP classes (Fig. 1) are reclassified into six types for the OHM model (Table 1). The barren or sparsely vegetated category refers to land with no plants, building or water cover, but usually porous ground. Pavement is impermeable, such as concrete or asphalt road surfaces. As the 500-m MODIS data do not provide land-cover fractions (e.g., proportions of pavement, building, vegetation) detailed analyses of both residential and commercial areas in Shanghai use similar methods to Ao et al. (2016a) for XJH. From this the urban and built-up land-cover class is assumed to be have the following characteristics: paved areas = 58%, buildings = 25%, grass = 14%, and trees = 3%. To explore the model performance at the national AWS stations, the surface-cover metadata [based on Shanghai Institute of Surveying and Mapping data; see Table 4 in Tan et al. (2015)] within 1 km of each is analyzed.

The surface temperature submodel of TsT2m [Eqs. (1)(25)] is run with a spatial resolution of 0.015° and 1-h time step for the area of interest. It is initialized with a spatial dataset that includes surface temperature (bare soil), 2-m air temperature and humidity, 10-m wind direction and wind speed, and station pressure corrected to sea level. Data from the 10 AWS (Fig. 1) operated by the Shanghai Meteorological Service [see Tan et al. (2015) for more details about these stations] are interpolated by inverse distance weighting (IDW; Shepard 1968) to the model grid scale to provide model input. IDW is selected as it does not require other parameters (e.g., influence radius) and with the sites fixed in space the interpolation result is deterministic and results of adjacent grid points are continuous. Alternatives, such as nearest-neighbor (NN) interpolation may have adjacent grid values that change abruptly. IDW is advantageous for routine interpolation. It has smaller errors than kriging or NN methods and is recommended for local-scale patterns of air temperature (e.g., Honjo et al. 2015; Li et al. 2017).

To assess the impact of the IDW interpolation, a clear-sky case (24–25 July 2016; cloud cover around 10%) is analyzed with and without the XJH station. The initial surface temperature of the XJH station area is 31.2°C (30.2°C) with (without) XJH and the corresponding simulated surface temperature at 1400 LT are 59.0° (61.1°)C. The pattern of temperature distribution changes little (Fig. 3). For the investigation of the impact of the IDW interpolation method we have chosen to initiate the model with ECMWF run time after sunset at 2000 LT during which the temperature trends were not abruptly changing.

Fig. 3.
Fig. 3.

Interpolated 2000 LT 24 Jul surface temperature (a) without XJH, (b) with XJH, and the resulting simulated surface temperature at 1400 LT (c) without XJH and (d) with XJH. Gray lines delineate Shanghai and its districts.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

b. Evaluation method

The surface temperature submodel of TsT2m is evaluated for the period from 1 September 2012 to 31 August 2013. Simulated incoming shortwave and net all-wave radiation fluxes are compared against the XJH flux tower observations (Ao et al. 2016a) by season and weather type. Modeled surface temperatures are compared to twice daily (about 1100 and 1300 LT) MODIS land surface temperature products (MOD11C1), which have a resolution of 5 km (https://modis.gsfc.nasa.gov/data/dataprod/mod11.php).

The near-surface air temperature submodel of TsT2m is initialized with interpolated surface and air temperatures from the 10 AWS observations. The ECMWF forecast 10-m wind direction and wind speed, humidity, air pressure, total and lower cloud cover (resolution = 0.125°), and 850-hPa temperature (resolution = 0.25°) are used to force this submodel through interpolation. The new spatial resolution of 0.015° for the model grid is obtained by bilinear interpolation. The 1-h temporal resolution (from 3 h) is obtained by linear interpolation. Simulations are conducted from 1 May 2016 to 31 October 2016 using the next 24-h forecast (1200 UTC) with a 1-h time step. Simulated air temperatures are compared to the 3-h ECMWF forecast 2-m air temperature.

To evaluate the performance, coefficients of determination R2, mean bias error (MBE), mean absolute error (MAE), and RMSE are used:
e26
e27
e28
e29
where Mi indicates simulated and Oi observed value for grid i of the total N grids or time steps; and and the respective means. The forecast accuracy (FA) or hit rate is also used (Zhang et al. 2009):
e30
where NT is the total number of time steps evaluated, NG is the number of times that the forecast temperature has an absolute error ≤1°C.

4. Results and discussion

a. Surface temperatures

The simulated hourly incoming shortwave flux and net all-wave radiation flux for clear days and cloudy days are compared with observations from the XJH tower (Fig. 4). Performance for both incoming shortwave and net all-wave radiation is better on clear days (R2 = 0.94 and 0.98) than on cloudy days (R2 = 0.68 and 0.70).

Fig. 4.
Fig. 4.

Simulated hourly (a),(b) incoming shortwave and (c),(d) net all-wave radiation flux compared with observed data at the XJH tower for (a),(c) clear days (N = 664) and (b),(d) cloudy days (N = 1478) during December 2012–August 2013. Regression slope (forced through the origin) and coefficients of determination R2 are given.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

Violin plots allow both the distribution and frequency of the errors to be seen (Fig. 5). On clear days simulated incoming shortwave radiation is higher than observed in spring and summer, but lower in autumn and winter. On cloudy days generally, simulated incoming shortwave radiation is higher than the observed data for all seasons. The hourly incoming shortwave flux MAE are 35, 23, 34, and 24 W m−2 in spring, summer, autumn, and winter, respectively (RMSE = 110, 123, 132, and 100 W m−2, respectively). The simulated transmissivity Tr on clear days are 0.70, 0.68, 0.61, and 0.66 on average in spring, summer, autumn, and winter, respectively, which follows the observed monthly median (1100–1300 LT) of 0.6–0.7 for clear conditions documented by Ao et al. (2016a). On cloudy days, values are 0.65, 0.58, 0.53, and 0.55 on average (spring, summer, autumn, and winter, respectively), which are higher than those derived from observations (0.59, 0.49, 0.48, and 0.47, respectively) (Ao et al. 2016a).

Fig. 5.
Fig. 5.

Violin plot showing maximum, median, and minimum whiskers and the kernel density of simulated bias for (a),(b) incoming shortwave radiation flux (K↓mod – K↓obs) and (c),(d) net all-wave radiation flux by seasons: (a),(c) clear (number of hours per season N) and (b),(d) cloudy days.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

Following Leroyer et al. (2011) and Zhang et al. (2011), modeled surface temperatures are compared to the MODIS skin temperature. The simulation values reproduce the general pattern well and give more detailed patterns of the surface urban heat island than is evident from MODIS (Fig. 6). The high temperature in the center of city does not have the surrounding gradient evident in the MODIS skin temperature (Fig. 6a). Significantly lower temperatures near the Huangpu River in the city center are seen in the simulated temperatures (Fig. 6b). The Yangtze River and Chongming Island have similar temperature with the waters and land being hard to distinguish in Fig. 6a, but easy to distinguish in Fig. 6b. To evaluate the modeled results, the data are stratified by surface type (water, urban, short vegetation, trees and shrubs, and barren). The simulated surface temperature for urban areas (N = 16 600) are 10.3°C warmer on average than those derived from MODIS, while the others are more similar [water: −0.2°C (N = 36 936), short vegetation 3.1°C (N = 25 411), trees and shrubs 0.1°C (N = 12 203), and barren ground −0.3°C (N = 1158)]. The larger difference associated with the urban surfaces is not unexpected as the coarser MODIS data are unlikely to only have the warmer built surfaces but will typically include cooler surface covers (e.g., grass, trees) as well.

Fig. 6.
Fig. 6.

Surface temperature (°C) on 12 Aug 2013 derived from (a) MODIS (1040 LT) and (b) simulations (1100 LT). Black lines delimit Shanghai and its districts.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

Four cloud-free MODIS images are analyzed (1300 LT 3 June, 1300 LT 12 July, 1100 LT 11 August, and 1100 LT 12 August 2013) to investigate the performance of the surface temperature simulation for different land covers (Fig. 7). Observations of surface temperature are available at 10 national AWS sites. These are also used to evaluate the grids where they occur using the barren surface cover type. The barren surface temperatures’ MBE and RMSE are −2.8° and 4.8°C, while for the MODIS data they are −5.2° and 6.1°C, respectively (N = 36).

Fig. 7.
Fig. 7.

Variation in surface temperature on four clear days in summer (1300 LT 3 Jun, 1300 LT 12 Jul, 1100 LT 11 Aug, and 1100 LT 12 Aug 2013) determined from MODIS (MOD) and simulation (SIM) by surface type (see Table 1 for classes in domain). Boxplots show interquartile range (IQR), 1.5 × IQR whiskers, median (red line), and mean (black dashed line).

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

b. Near-surface air temperatures

To evaluate the near-surface air temperature of Shanghai, a day with high temperatures is selected (27 July 2016). A mix of cloud and sun occurred, with low wind speeds from inland (southwest) and near the sea from the southeast. The downscaled air temperature distribution is compared to air temperatures from 110 AWS (Fig. 8) and the ECMWF high-resolution 2-m temperature prediction with 3-h temporal resolution and 0.125° spatial resolution. The model is initialized at 2000 LT 26 July 2016 (i.e., the same time as the as ECMWF prediction) and temperatures simulated for 24 h. At 1400 LT 27 July 2016, the highest air temperature areas are in the downtown area of Shanghai for all three approaches: the observations are the warmest (~41°C), TsT2m is 39.8°C, and ECMWF (~37°C) the coolest. Both ECMWF and TsT2m simulate lower temperatures near the sea, because of a sea breeze flowing in from the cooler sea surface. ECMWF has a secondary hot spot in the southwest corner of Shanghai, but this is not evident in the AWS or the TsT2m results.

Fig. 8.
Fig. 8.

Air temperature at 2 m above ground level in Shanghai at 1400 LT 27 Jul 2016 at (a) AWS observations, (b) ECMWF (1200 UTC) forecast, and (c) TsT2m [initialized with 1200 UTC (2000 LT) 26 Jul 2016 ECMWF forecast].

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

For 1 May–31 October 2016, the near-surface air temperatures (Fig. 9) from ECMWF are mostly cooler than observed, except for some stations near the sea. In contrast, the TsT2m results are mostly warmer than observed, especially in the urban area. The RMSETsT2m are more spatially consistent than RMSEECMWF, as the latter are larger in the urban area than in the nonurban area. In August (Figs. 9c,d), the hottest month in Shanghai, the MBEECM is <−1.0°C for 77 stations, and for 15 stations RMSEECM is >2.0°C. Most of these are stations in the urban area. However, 96 stations have a MAETsT2m < 1.0°C. The four stations with RMSETsT2m > 2.0°C are located close to the coast.

Fig. 9.
Fig. 9.

Variation of MBE (color) and RMSE (dot size) for 2-m air temperature assessed with AWS data for (a),(b) May–October and (c),(d) August 2016 simulated by (a),(c) TsT2m and (b),(d) ECMWF.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

The average FA for the near-surface air temperature (Table 2) is better for FATsT2m than FAECM at almost all stations. In XJH and Chongming (CM) (Fig. 1), the full suite of TsT2m results are much better than those from the original ECMWF as both sites have a 20% higher FA. This also occurs in the MBE and RMSE (Fig. 10 shows XJH; others not shown). The FATsT2m is better than FAECM for every month, and much better in July (by 12.9%), August (24.8%) and September (11.4%). This indicates the TsT2m performance is much better on hot days (MBETsT2m < 1°C; RMSETsT2m < 3°C in each month) and better than the ECMWF forecast.

Table 2.

Air temperature (at 2 m) FA (%) of ECMWF (ECM) and the TsT2m (T2m) at each station (Fig. 1) by month (May–October 2016).

Table 2.
Fig. 10.
Fig. 10.

Boxplot [IQR, median (red line), mean (black dotted line), whiskers 1.5 × IQR] of hourly air temperature bias at XJH using TsT2m and ECMWF from May (5) to October (10).

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

To evaluate the effect of the surface cover, the fractions of building and pavement (%) for each station are determined (Fig. 11a, section 3a). FAECM has a negative correlation with impervious cover (Fig. 11b) (R2 = 0.45, t-test significance p = 0.002 for α = 0.05). However, there is little influence of urban fractional cover on FATsT2m. Thus, the ECMWF prediction of temperature in urban areas (e.g., Xujiahui, Minghang, Songjiang) is modified from that of nonurban area (e.g., Fengxian, Qingpu, Pudong). The new TsT2m 2-m air temperature, by accounting for urban land cover, results in much improved performance for areas with road and built-up cover.

Fig. 11.
Fig. 11.

Land-cover characteristics (a) within 1-km radius of the national stations (Fig. 1) and (b) urban percentage (paved + building) with FA of 2-m air temperature from May to October 2016.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

c. Application

A range of potential weather and climate service applications are possible with the new model. These include regular weather forecasts at higher spatial resolution and climate services mitigation applications.

Here the impact of changing land cover for heat mitigation is assessed by analyzing the hottest day in Shanghai in 2016 using the same forcing data as in section 4b. The air temperature is simulated with the original land cover (section 3a) to characterize urban grids (i.e., paved areas = 58%, buildings = 25%, grass = 14%, trees = 3%). This is used as a reference temperature. Two more sets of simulations are conducted that have decreasing paved proportions (58%–8%, 10% steps) that are replaced with (i) grass or (ii) trees. The maximum and minimum daily air temperature for each urban grid (T2m,max and T2m,min, respectively) are compared (simulation minus reference temperature simulation) to obtain the change (ΔT2m,max and ΔT2m,min, respectively) for each land-cover change (Fig. 12).

Fig. 12.
Fig. 12.

Impact of change in surface cover [section 4c; pavement replaced by grass (green) or tree (blue)] on simulated (by TsT2m) daily air temperature (27 Jul 2016) (a) maximum (ΔT2m,max) and (b) minimum (ΔT2m,min). Spatial difference boxplot for all urban grids with IQR, and 1.5 × IQR whiskers, median (red line), and mean (red diamond). Note y-scale differences.

Citation: Journal of Applied Meteorology and Climatology 57, 10; 10.1175/JAMC-D-17-0255.1

The mean T2m,max decreases with an increase in both grass and trees; therefore, ΔT2m,max becomes larger with an increase of grass or tree proportion. The reduction in T2m,max is on average greater for an increase in trees [0.45°C (10%)−1] than grass [0.37°C (10%)−1]. Similarly, the mean T2m,min decreases with an increase in both grass and trees (Fig. 12b). The decrease associated with the trees [0.063°C (10%)−1] is larger than with grass [0.026°C (10%)−1] but smaller than the impact on the maximum temperature (cf. Figs. 12a,b).

5. Discussion and conclusions

In this paper a new model (TsT2m) to downscale coarse NWP (e.g., ECMWF) forecast surface and air temperature based on underlying surface cover types is proposed. The approach provides higher temporal and spatial resolution outputs than traditional statistical downscaling. The model is assessed using AWS and MODIS data in Shanghai. For air temperature, the performance of the TsT2m is better than the coarser (spatially and temporally) ECMWF forecast, and reflects better the impact of the environment on the near-surface (2 m) temperature. In Shanghai, enhanced performance is most evident on hot days. From May to October, the ECMWF 2-m air temperature forecast in the urban area is systematically cooler, while the mean and median errors of TsT2m are both smaller. As TsT2m mostly uses conventional meteorological data (temperature, humidity, wind speed, air pressure, cloud cover, etc.), and requires limited computing resources, it has great potential to be translated and extended to other cities and applications.

TsT2m performs less well in winter than summer. This may result from the assumption that air near ground is well mixed, which it less likely under more stable conditions. Currently, TsT2m does not consider hydrological processes. As surface moisture processes are important, their inclusion is important in future improvements. Heterogeneity of urban morphology influences land surface temperature (Guo et al. 2016) and this also will be considered in future studies. As TsT2m is tested only in Shanghai, evaluation in other cities is important to establish the robustness of the approach and its broader utility.

Given a lack of high-resolution observations to initialize the model, it is necessary to interpolate data, which inevitably introduces biases in the model simulations. Here, IDW is chosen given the local scale and as it yields continuous results. Obviously, inclusion of more AWSs in the interpolation will result in a better spatial representation.

The storage heat flux is sometimes overestimated with the OHM mode (especially in the afternoon). This may lead to an overestimate of Ts. Given the simplicity of OHM. and the fact it does not account for impacts of surface moisture and wind speed variations (Grimmond and Oke 1999), the first-order results of this study suggest that assessing the new analytical OHM (AnOHM) model (Sun et al. 2017) is warranted to increase OHM coefficients in future simulations.

TsT2m has the potential to be used to study finer-scale urban processes and conditions, for example: thermodynamic conditions of local convection, the effects of which are too small scale to be captured by a NWP models. Thus, TsT2m provides a useful tool for evaluating the impact of urbanization on the thermal environment and to inform city planning and design. TsT2m is developed with the aim of improving climate and weather services for heat stress conditions, particularly for those most vulnerable. The improved spatial resolution, particularly in the densest urban areas, will enable improved predictions and interventions and aid long-term planning for extreme heat conditions.

Acknowledgments

This work is supported by The National Natural Science Foundation of China Grants 41775019 (J.G. Tan) and 41475040 (M.B. Du) and the UK–China Research and Innovation Partnership Fund through the Met Office Climate Science for Service Partnership (CSSP) China as part of the Newton Fund (Grimmond, D.W. Liu, and J. Peng UK visit). Useful discussions with Ting Sun (Reading) are appreciated.There are no conflicts of interest.

APPENDIX A

Symbols

Symbols used in this paper are listed in Table A1.

Table A1.

List of symbols.

Table A1.

REFERENCES

  • Allen, L., F. Lindberg, and C. S. B. Grimmond, 2011: Global to city scale urban anthropogenic heat flux: Model and variability. Int. J. Climatol., 31, 19902005, https://doi.org/10.1002/joc.2210.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Anandakumar, K., 1999: A study on the partition of net radiation into heat fluxes on a dry asphalt surface. Atmos. Environ., 33, 39113918, https://doi.org/10.1016/S1352-2310(99)00133-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ao, X., C. S. B. Grimmond, D. Liu, Z. Han, P. Hu, Y. Wang, X. Zhen, and J. Tan, 2016a: Radiation fluxes in a business district of Shanghai, China. J. Appl. Meteor.Climatol., 55, 24512468, https://doi.org/10.1175/JAMC-D-16-0082.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ao, X., C. S. B. Grimmond, Y. Chang, D. Liu, Y. Tang, and P. Hu, 2016b: Heat, water and carbon exchanges in the tall megacity of Shanghai: Challenges and results. Int. J. Climatol., 36, 46084624, https://doi.org/10.1002/joc.4657.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Arnfield, A. J., and C. S. B. Grimmond, 1998: An urban energy budget model and its application to urban storage heat flux modeling. Energy Build., 27, 6168, https://doi.org/10.1016/S0378-7788(97)00026-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Asaeda, T., and V. T. Ca, 1993: The subsurface transport of heat and moisture and its effect on the environment: A numerical model. Bound.-Layer Meteor., 65, 159179, https://doi.org/10.1007/BF00708822.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Atwater, M. A., and P. S. Brown, 1974: Numerical computations of the latitudinal variation of solar radiation for an atmosphere of varying opacity. J. Appl. Meteor., 13, 289297, https://doi.org/10.1175/1520-0450-13.2.289.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Atwater, M. A., and J. T. Ball, 1981: A surface solar radiation model for cloudy atmospheres. Mon. Wea. Rev., 109, 878888, https://doi.org/10.1175/1520-0493(1981)109<0878:ASSRMF>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Basu, R., and J. M. Samet, 2002: Relation between elevated ambient temperature and mortality: A review of the epidemiologic evidence. Epidemiol. Rev., 24, 190202, https://doi.org/10.1093/epirev/mxf007.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bechtel, B., Z. Klemen, and H. Gholamali, 2012: Downscaling land surface temperature in an urban area: A case study for Hamburg, Germany. Remote Sens., 4, 31843200, https://doi.org/10.3390/rs4103184.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Best, M. J., 1998: A model to predict surface temperatures. Bound.-Layer Meteor., 88, 279306, https://doi.org/10.1023/A:1001151927113.

  • Best, M. J., and C. S. B. Grimmond, 2015: Key conclusions of the First International Urban Land Surface Model Comparison Project. Bull. Amer. Meteor. Soc., 96, 805819, https://doi.org/10.1175/BAMS-D-14-00122.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bhumralkar, C. M., 1975: Numerical experiments on the computation of ground surface temperature in an atmospheric general circulation model. J. Appl. Meteor., 14, 12461258, https://doi.org/10.1175/1520-0450(1975)014<1246:NEOTCO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blackadar, A. K., 1976: Modeling the nocturnal boundary layer. Preprints, Third Symp. on Atmospheric Turbulence, Diffusion, and Air, Raleigh, NC, Amer. Meteor. Soc., 46–49.

  • Bonafoni, S., 2016: Downscaling of Landsat and MODIS land surface temperature over the heterogeneous urban area of Milan. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens., 9, 20192027, https://doi.org/10.1109/JSTARS.2016.2514367.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Chen, F., and J. Dudhia, 2001: Coupling an advanced land surface–hydrology model with the Penn State–NCAR MM5 modeling system. Part I: Model implementation and sensitivity. Mon. Wea. Rev., 129, 569585, https://doi.org/10.1175/1520-0493(2001)129<0569:CAALSH>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Cohen, W. B., and Coauthors, 2006: MODIS land cover and LAI collection 4 product quality across nine states in the western hemisphere. IEEE Trans. Geosci. Remote Sens., 44, 18431857, https://doi.org/10.1109/TGRS.2006.876026.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Diefenderfer, B. K., I. L. Al-Qadi, and S. D. Diefenderfer, 2006: Model to predict pavement temperature profile: Development and validation. J. Transp. Eng., 132, 162167, https://doi.org/10.1061/(ASCE)0733-947X(2006)132:2(162).

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Doll, D., J. K. S. Ching, and J. Kaneshiro, 1985: Parameterization of subsurface heating for soil and concrete using net radiation data. Bound.-Layer Meteor., 32, 351372, https://doi.org/10.1007/BF00122000.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gabey, A. M., C. S. B. Grimmond, and I. Capel-Timms, 2018: Anthropogenic heat flux: Advisable spatial resolutions when input data are scarce. Theor. Appl. Climatol., https://doi.org/10.1007/s00704-018-2367-y, in press.

    • Search Google Scholar
    • Export Citation
  • Gallo, P., A. L. McNab, T. R. Karl, J. F. Brown, J. J. Hood, and J. D. Tarpley, 1993: The use of NOAA AVHRR data for assessment of the urban heat island effect. J. Appl. Meteor., 32, 899908, https://doi.org/10.1175/1520-0450(1993)032<0899:TUONAD>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Grimm, N. B., S. H. Faeth, N. E. Golubiewski, C. L. Redman, J. G. Wu, X. M. Bai, and J. M. Briggs, 2008: Global change and the ecology of cities. Science, 319, 756760, https://doi.org/10.1126/science.1150195.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Grimmond, C. S. B., and T. R. Oke, 1999: Heat storage in urban areas: local-scale observations and evaluation of a simple model. J. Appl. Meteor., 38, 922940, https://doi.org/10.1175/1520-0450(1999)038<0922:HSIUAL>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Grimmond, C. S. B., H. A. Cleugh, and T. R. Oke, 1991: An objective urban heat storage model and its comparison with other schemes. Atmos. Environ., 25B, 311326, https://doi.org/10.1016/0957-1272(91)90003-W.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Grimmond, C. S. B., and Coauthors, 2010: The international Urban Energy Balance Models Comparison Project: First results from phase 1. J. Appl. Meteor. Climatol., 49, 12681292, https://doi.org/10.1175/2010JAMC2354.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Grimmond, C. S. B., and Coauthors, 2011: Initial results from phase 2 of the International Urban Energy Balance Model Comparison. Int. J. Climatol., 31, 244272, https://doi.org/10.1002/joc.2227.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Guo, G., X. Zhou, Z. Wu, R. Xiao, and Y. Chen, 2016: Characterizing the impact of urban morphology heterogeneity on land surface temperature in Guangzhou, China. Environ. Modell. Software, 84, 427439, https://doi.org/10.1016/j.envsoft.2016.06.021.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Herb, W. R., B. Janke, O. Mohseni, and H. G. Stefan, 2008: Ground surface temperature simulation for different land covers. J. Hydrol., 356, 327343, https://doi.org/10.1016/j.jhydrol.2008.04.020.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Honjo, T., H. Yamato, T. Mikami, and C. S. B. Grimmond, 2015: Network optimization for enhanced resilience of urban heat island measurements. Sustainable Cities Soc., 19, 319330, https://doi.org/10.1016/j.scs.2015.02.004.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hu, Z., and S. Islam, 1995: Prediction of ground surface temperature and soil moisture content by the force-restore method. Water Resour. Res., 31, 25312539, https://doi.org/10.1029/95WR01650.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Järvi, L., C. S. B. Grimmond, and M. Taka, 2014: Development of the Surface Urban Energy and Water Balance Scheme (SUEWS) for cold climate cities. Geosci. Model Dev., 7, 16911711, https://doi.org/10.5194/gmd-7-1691-2014.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Jin, M. S., W. Kessomkiat, and G. Pereira, 2011: Satellite-observed urbanization characters in Shanghai, China: Aerosols, urban heat island effect, and land–atmosphere interactions. Remote Sens., 3, 8399, https://doi.org/10.3390/rs3010083.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kalma, J. D., T. R. McVicar, and M. F. McCabe, 2008: Estimating land surface evaporation: A review of methods using remotely sensed surface temperature data. Surv. Geophys., 29, 421469, https://doi.org/10.1007/s10712-008-9037-z.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kim, Y. H., and J. J. Baik, 2005: Spatial and temporal structure of the urban heat island in Seoul. J. Appl. Meteor., 44, 591605, https://doi.org/10.1175/JAM2226.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Leroyer, S., S. Belair, and J. Mailhot, 2011: Microscale numerical prediction over Montreal with the Canadian External Urban Modeling System. J. Appl. Meteor. Climatol, 50, 24102428, https://doi.org/10.1175/JAMC-D-11-013.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Li, R., Y. Zhang, and W. U. Jie, 2017: Interpolation methods for temperature distribution of urban communities in southern China (in Chinese). Build. Sci., 33 (12), 2737.

    • Search Google Scholar
    • Export Citation
  • Lindberg, F., and Coauthors, 2018: Urban Multi-scale Environmental Predictor (UMEP): An integrated tool for city-based climate services. Environ. Modell. Software, 99, 7087, https://doi.org/10.1016/j.envsoft.2017.09.020.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Loridan, T., C. S. B. Grimmond, B. D. Offerle, D. T. Young, T. E. Smith, L. Järvi, and F. Lindberg, 2011: Local-Scale Urban Meteorological Parameterization Scheme (LUMPS): Longwave radiation parameterization and seasonality-related developments. J. Appl. Meteor. Climatol., 50, 185202, https://doi.org/10.1175/2010JAMC2474.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Luo, Y., R. S. Loomis, and T. C. Hsiao, 1992: Simulation of soil temperature in crops. Agric. For. Meteor., 61, 2338, https://doi.org/10.1016/0168-1923(92)90023-W.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McCaughey, J. H., 1985: Energy balance storage terms in a mature mixed forest at Petawawa, Ontario—A case study. Bound.-Layer Meteor., 31, 89101, https://doi.org/10.1007/BF00120036.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McDonald, J. E., 1960: Direct absorption of solar radiation by atmospheric water vapor. J. Meteor., 17, 319328, https://doi.org/10.1175/1520-0469(1960)017<0319:DAOSRB>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Meyers, T. P., and R. F. Dale, 1983: Predicting daily insolation with hourly cloud height and coverage. J. Climate Appl. Meteor., 22, 537545, https://doi.org/10.1175/1520-0450(1983)022<0537:PDIWHC>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Meyn, S. K., and T. R. Oke, 2009: Heat fluxes through roofs and their relevance to estimates of urban heat storage. Energy Build., 41, 745752, https://doi.org/10.1016/j.enbuild.2009.02.005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mihailović, D. T., 1991: A model for prediction of the soil temperature and the soil moisture content in three layers. Z. Meteor., 41, 2933.

    • Search Google Scholar
    • Export Citation
  • Mihailović, D. T., 1996: Description of a land-air parameterization scheme (LAPS). Global Planet. Change, 13, 207215, https://doi.org/10.1016/0921-8181(95)00048-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mihailović, D. T., and J. Eitzinger, 2007: Modelling temperatures of crop environment. Ecol. Modell., 202, 465475, https://doi.org/10.1016/j.ecolmodel.2006.11.009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Molteni, F., R. Buizza, T. N. Palmer, and T. Petroliagis, 1996: The ECMWF Ensemble Prediction System: Methodology and validation. Quart. J. Roy. Meteor. Soc., 122, 73119, https://doi.org/10.1002/qj.49712252905.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Narita, K., T. Sekine, and T. Tokuoka, 1984: Thermal properties of urban surface materials: Study on heat balance at asphalt pavement. Geogr. Rev. Japan, 57, 639651, https://doi.org/10.4157/grj1984a.57.9_639.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nichol, J. E., W. Y. Fung, K. Lam, and M. S. Wong, 2009: Urban heat island diagnosis using ASTER satellite images and in situ air temperature. Atmos. Res., 94, 276284, https://doi.org/10.1016/j.atmosres.2009.06.011.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Novak, M. D., 1981: The moisture and thermal regimes of a bare soil in the Lower Fraser Valley during spring. Ph.D. thesis, University of British Columbia, Vancouver, BC, Canada, 153 pp.

  • Offerle, B., C. S. B. Grimmond, and T. R. Oke, 2003: Parameterization of net all-wave radiation for urban areas. J. Appl. Meteor., 42, 11571173, https://doi.org/10.1175/1520-0450(2003)042<1157:PONARF>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Oke, T. R., 1973: City size and the urban heat island. Atmos. Environ., 7, 769779, https://doi.org/10.1016/0004-6981(73)90140-6.

  • Oke, T. R., 1988: The urban energy balance. Prog. Phys. Geogr., 12, 471508, https://doi.org/10.1177/030913338801200401.

  • Oke, T. R., 1995: The heat island of the urban boundary layer: Characteristics, causes and effects. Wind Climate in Cities, J. E. Cermak et al., Eds., NATO ASI Series E, Vol. 227, Kluwer Academic, 81–107.

    • Crossref
    • Export Citation
  • O’Malley, C., P. A. E. Piroozfarb, E. R.P. Farr, and J. Gates, 2014: An investigation into minimizing urban heat island (UHI) effects: A UK perspective. Energy Procedia, 62, 7280, https://doi.org/10.1016/j.egypro.2014.12.368.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Pan, X., X. Li, X. Shi, X. Han, and L. Luo, 2012: Dynamic downscaling of near-surface air temperature at the basin scale using WRF—A case study in the Heihe River basin, China. Front. Earth Sci., 6, 314323, https://doi.org/10.1007/s11707-012-0306-2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Raich, J. W., and W. H. Schlesinger, 1992: The global carbon dioxide flux in soil respiration and its relationship to vegetation and climate. Tellus, 44B, 8199, https://doi.org/10.3402/tellusb.v44i2.15428.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Roberts, S. M., T. R. Oke, C. S. B. Grimmond, and J. A. Voogt, 2006: Comparison of four methods to estimate urban heat storage. J. Appl. Meteor. Climatol., 45, 17661781, https://doi.org/10.1175/JAM2432.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Roth, M., T. R. Oke, and W. J. Emery, 1989: Satellite derived urban heat islands from three coastal cities and the utilization of such data in urban climatology. Int. J. Remote Sens., 10, 16991720, https://doi.org/10.1080/01431168908904002.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Shao, J., and P. J. Lister, 1996: An automated nowcasting model of road surface temperature and state for winter road maintenance. J. Appl. Meteor., 35, 13521361, https://doi.org/10.1175/1520-0450(1996)035<1352:AANMOR>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Shepard, D., 1968: A two-dimensional interpolation function for regularly spaced data. Proc. 23rd National Conf. of American Computing Machinery, Princeton, NJ, Association for Computing Machinery, 517–524.

    • Crossref
    • Export Citation
  • Smith, W. L., 1966: Note on the relationship between total precipitable water and surface dew point. J. Appl. Meteor., 5, 726727, https://doi.org/10.1175/1520-0450(1966)005<0726:NOTRBT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • South, C., C. S. B. Grimmond, and C. P. Wolfe, 1998: Evapotranspiration rates from wetlands with different disturbance histories: Indiana Dunes National Lakeshore. Wetlands, 18, 216229, https://doi.org/10.1007/BF03161657.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Stone, J. B., and M. O. Rodgers, 2001: Urban form and thermal efficiency: How the design of cities influences the urban heat island effect. J. Amer. Plann. Assoc., 67, 186198, https://doi.org/10.1080/01944360108976228.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Sun, T., Z. H. Wang, W. Oechel, and C. S. B. Grimmond, 2017: The Analytical Objective Hysteresis Model (AnOHM v1. 0): Methodology to determine bulk storage heat flux coefficients. Geosci. Model Dev., 10, 28752890, https://doi.org/10.5194/gmd-10-2875-2017.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taesler, R., 1980a: Studies of the development and thermal structure of the urban boundary layer in Uppsala. Part I: Experimental program. Meteorological Institution, Uppsala University Rep. 61, 63 pp.

  • Taesler, R., 1980b: Studies of the development and thermal structure of the urban boundary layer in Uppsala. Part II: Data analysis and results. Meteorological Institution, Uppsala University Rep. 61, 181 pp.

  • Taha, H., 1999: Modifying a mesoscale meteorological model to better incorporate urban heat storage: A bulk-parameterization approach. J. Appl. Meteor., 38, 466473, https://doi.org/10.1175/1520-0450(1999)038<0466:MAMMMT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Tan, J., and Coauthors, 2015: Urban integrated meteorological observations: Practice and experience in Shanghai, China. Bull. Amer. Meteor. Soc., 96, 85102, https://doi.org/10.1175/BAMS-D-13-00216.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Taylor, K. E., 2001: Summarizing multiple aspects of model performance in a single diagram. J. Geophys. Res., 106, 71837192, https://doi.org/10.1029/2000JD900719.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Trumbore, S., O. A. Chadwick, and R. Amundson, 1996: Rapid exchange between soil carbon and atmospheric carbon dioxide driven by temperature change. Science, 272, 393396, https://doi.org/10.1126/science.272.5260.393.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Van De Griend, A. A., and M. Owe, 1993: On the relationship between thermal emissivity and the normalized difference vegetation index for natural surfaces. Int. J. Remote Sens., 14, 11191131, https://doi.org/10.1080/01431169308904400.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wang, X., Y. Li, K. Wang, X. Yang, and P. W. Chan, 2017: A simple daily cycle temperature boundary condition for ground surfaces in CFD predictions of urban wind flows. J. Appl. Meteor. Climatol., 56, 29632980, https://doi.org/10.1175/JAMC-D-17-0095.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Weng, Q., D. Lu, and J. Schubring, 2004: Estimation of land surface temperature–vegetation abundance relationship for urban heat island studies. Remote Sens. Environ., 89, 467483, https://doi.org/10.1016/j.rse.2003.11.005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Winkler, J. A., J. P. Palutikof, J. A. Andresen, and C. M. Goodess, 1997: The simulation of daily temperature time series from GCM output. Part II: Sensitivity analysis of an empirical transfer function methodology. J. Climate, 10, 25142532, https://doi.org/10.1175/1520-0442(1997)010<2514:TSODTT>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Xu, W., M. J. Wooster, and C. S. B. Grimmond, 2008: Modelling of urban sensible heat flux at multiple spatial scales: A demonstration using airborne hyperspectral imagery of Shanghai and a temperature–emissivity separation approach. Remote Sens. Environ., 112, 34933510, https://doi.org/10.1016/j.rse.2008.04.009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Xue, Y., P. J. Sellers, J. L. Kinter, and J. Shukla, 1991: A simplified biosphere model for global climate studies. J. Climate, 4, 345364, https://doi.org/10.1175/1520-0442(1991)004<0345:ASBMFG>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Yoshida, A., K. Tominaga, and S. Watatani, 1990: Field measurements on energy balance of an urban canyon in the summer season. Energy Build., 15, 417423, https://doi.org/10.1016/0378-7788(90)90016-C.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Yow, D. M., 2007: Urban heat islands: Observations, impacts, and adaptation. Geogr. Compass, 1, 12271251, https://doi.org/10.1111/j.1749-8198.2007.00063.x.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zerenner, T., V. Venem, P. Friederichs, and C. Simmer, 2016: Downscaling near-surface atmospheric fields with multi-objective genetic programming. Environ. Modell. Software, 84, 8598, https://doi.org/10.1016/j.envsoft.2016.06.009.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zhang, D. L., Y. X. Shou, R. R. Dickerson, and F. Chen, 2011: Impact of upstream urbanization on the urban heat island effects along the Washington–Baltimore corridor. J. Appl. Meteor. Climatol., 50, 20122029, https://doi.org/10.1175/JAMC-D-10-05008.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zhang, Q., A. Y. Xiong, J. Y. Zhang, M. N. Feng, and B. M. Wang, 2009: Preliminary study on the scoring methods of cloud-free rainfall/snowfall and air temperature forecasts (in Chinese). J. Appl. Meteor. Sci., 20, 692698.

    • Search Google Scholar
    • Export Citation
Save