Infuence of urbanization on the thermal environment of meteorological station:Satellite-observed evidence

,YANG Yun-Jin,*

aKey Laboratory of Atmospheric Sciences and Satellite Remote Sensing of Anhui Province,Anhui Institute of Meteorological Sciences,Hefei 230031,China

bWuhu Meteorological Bureau,Wuhu,241000,China

cNanjing University of Information Science and Technology,Nanjing 210044,China

Infuence of urbanization on the thermal environment of meteorological station:Satellite-observed evidence

SHI Taoa,b,HUANG Yonga,WANG Hongc,SHI Chun-Ea,YANG Yuan-Jiana,*

aKey Laboratory of Atmospheric Sciences and Satellite Remote Sensing of Anhui Province,Anhui Institute of Meteorological Sciences,Hefei 230031,China

bWuhu Meteorological Bureau,Wuhu,241000,China

cNanjing University of Information Science and Technology,Nanjing 210044,China

In this paper,fve national meteorological stations in Anhui province are taken as typical examples to explore the effects of local urbanization on their thermal environment by using Landsat data from 1990 to 2010.Satellite-based land use/land cover(LULC),land surface temperature (LST),normalized difference vegetation index(NDVI)are used to investigate the effects.The study shows that LULC around meteorological stations changed signifcantly due to urban expansion.Fast urbanization is the main factor that affects the spatial-temporal distribution of thermal environment around meteorological stations.Moreover,the normalized LSTand NDVI exhibit strong inverse correlations around meteorological stations,so the variability of LST can be monitored through evaluating the variability of NDVI.In addition,station-relocation plays an important role in improving representativeness of thermal environment.Notably,the environment representativeness was improved,but when using the data from the station to study climate change,the relocation-induced inhomogeneous data should be considered and adjusted.Consequently, controlling the scale and layout of the urban buildings and constructions around meteorological stations is an effective method to ameliorate observational thermal environment and to improve regional representativeness of station observation.The present work provides observational evidences that high resolution Landsat images can be used to evaluate the thermal environment of meteorological stations.

Urbanization;Thermal environment;Representativeness;Land surface temperature;Normalized difference vegetation index(NDVI)

1.Introduction

Driven by the recent thirty-year economic booming,China has undergone rapid development and urbanization.Many meteorological stations used to be in rural area now are in urban area and their observational environments have changeddramatically(Ren et al.,2010;Zhang et al.,2010;Ren and Ren,2011;Shao et al.,2011;Yang et al.,2011,2013;Li et al.,2015).The changed thermal environment around meteorological stations signifcantly infuences the observations(Yan et al.,2010;Li et al.,2012;Shao et al.,2011;Ren and Ren,2011),which will further disturb local weather and climate analysis,such as the evaluation of heat island effect, one of the main features of a modern city(Kalnay and Cai, 2003;Li et al.,2004;Zhou et al.,2004;Chen et al.,2007; Ren et al.,2007,2008;Shi et al.,2011a;Zhao et al.,2013; Yang et al.,2013).For these reasons,objectively quantifying changes in thermal environment is crucial to evaluate the representativeness of meteorological stations.

Land surface temperature(LST)is an important index to represent the thermal environment around meteorologicalstations.Satellite remote sensing,as a recently developed technique,provides a unique opportunity to monitor and study macroscopical and dynamic-continuous LST in different spatial scales.Satellite-derived LST images have been widely used in land use classifcation,urban heat island research, thermal environment and hydrological investigation in an urban or even larger scales(Yang,2000;Weng,2001;Zhang et al.,2005,2011;Hung et al.,2006;Shi et al.,2011b, 2013).However,only in recent years have satellite-derived LST data been used specifcally for evaluating thermal environment and observational stations'representativeness in China.For example,by using 1-km resolution MODIS LST dataset,Wang et al.(2011)evaluated the representativeness of 142 weather stations and investigated the relationship between the representativeness of a station and its surrounding conditions.Ren et al.(2010)and Ren and Ren(2011)proposed to combine remote sensing images from Google Earth and MODIS LST data to evaluate the representativeness of observational stations and investigate the effect of urban heat islands.Using Landsat remote sensing data,Li et al.(2015) calculated the normalized LST and the heat effect contribution index(HECI)of different land use/land cover(LULC) type to classify the stations'observational environment into three types(urban,sub-urban and rural),and these two indexes could be used conveniently,effectively and quantitatively to choose a reference station when analyzing observational data in weather and climate research.However,thermal environment and representativeness around many Chinese meteorological stations are still not clear,especially the ones located in the southeastern China where has undergone rapid urbanization.Therefore,it is essential to use satellite-derived LST to monitor spatial-temporal variations of thermal environment and representativeness of these meteorological stations.

In this study,taking fve typical meteorological stations as samples in Anhui province,a southeastern province in China, where the tremendous growth of urban sprawl,population, vehicles and economy have occurred since the 1990s(Shi et al.,2008;Li et al.,2012;Yang et al.,2011,2013),spatiotemporal variations of LULC and thermal environment (LST)around meteorological stations are systematically explored by using Landsat remote sensing data.Moreover,the effects of urbanization on the thermal environment of meteorological stations are quantitatively evaluated by using LULC change around these meteorological stations.Finally,the relationship between LST and normalized difference vegetation index(NDVI)are quantitatively investigated.

2.Data and method

2.1.Data

The Landsat-5 remote sensing data used in this study are obtained from the Open Spatial Data Sharing Project,which was launched by the Institute of Remote Sensing and Digital Earth of the Chinese Academy of Sciences(http://ids.ceode. ac.cn/en/).The TM(ETM)sensor which is on board of the Landsat-5 satellite has seven bands,and the sixth band(TM6, with the band wavelength 10.40-12.50 μm)is selected here for LST retrieval(Li et al.,2015).In addition,the fourth band (TM4,with the near infrared band wavelength 0.62-0.69 μm) and the third band(TM3,with the red band wavelength 0.76-0.96 μm)are also selected here for NDVI retrieval(Shi et al.,2011b).In order to better capture the local vegetation information in a similar season,only data during the vegetation growth periods(i.e.,May to September in China)in 1990, 2000 and 2010 are selected.

To investigate the temporal-spatial variability of the meteorological stations'surrounding environment,the LULC types in the buffer zone of the station are categorized into three types,vegetation(including farmland,forest and grass land), water(including lakes,rivers and pools),and construction (including buildings and roads),which are derived exactly through supervised classifcation together with visual interpretation(Li et al.,2015).

Five stations under rapid urbanization are selected as typical cases(Fig.1),i.e.,Suzhou(SZ),Hefei(HF),Chuzhou (CZ),Anqing(AQ)and Wuhu(WH);and fve reference stations in rural areas are also selected corresponding to the fve urban stations,i.e.,Lingbi(LB),Feixi(FX),Quanjiao(QJ), Huaining(HN)and Nanling(NL).In addition,the underlying surfaces around the fve reference stations experienced little changes during 1990-2010.The annual mean air temperatures for the period 1980-2010 were recorded at the above mentioned 10 meteorological stations in Anhui province.

2.2.Method

2.2.1.Critical radius for the buffer zone around station

The buffer zone here is defned as the area within 2 km from the station based on two following reasons:1)In Anhui province,through feld survey it has been found that most of the meteorological stations were initially built 2 km away from the edge of cities.Therefore,by studying the 2 km radius area of a station,the infuence of city expansion on the stations'observational environment can be captured.2)According to recent studies(Ren et al.,2010;Yang et al.,2013;Li et al.,2015),the infuence of urbanization on a meteorological station in China can be refected by the observational environment within its 2-km radius.Therefore,2 km is chosen as the critical radius for the buffer zone.

2.2.2.Retrieval of LST

Brightness temperature is usually retrieved from thermal infrared band data.First,digital number values(QDN)of TM6 images are converted to thermal radiation intensity(R)(Weng, 2003),following,

Here,QDNis the digital number of a pixel in the TM6 image. RMAXand RMINare the maximum and minimum radiation intensity received by TM6,which are set as 1.896 and 0.1534 mW cm-2sr-1mw-1μm-1,respectively.The maximum and minimum digital number QMAXand QMINare 255 and 0,respectively.As such,Equation(1)can be written as:

Brightness temperature TB(units:K)for a black body is then calculated following:

Here,K1=60.776 and K2=1260.56.However,the Earth is not a black body and has different underlying surface conditions at difference places.LST,which is different from TB,is then calculated using(Artis and Carnahan,1982;Weng,2001)

Here,λ(≈11.5 μm)is the effective wavelength of TM6. ρ=hc/σ=1.438×10-2mK,where σ is the Boltzmann constant(1.38× 10-23J K-1),h is the Planck constant (6.626 × 10-34J s) and c is the lightvelocity (2.998×108m s-1).ε is emissivity whose value is 0.9925, 0.95 and 0.92 for water,vegetation and construction surfaces, respectively(Nichol,1994).

2.2.3.Normalized LST

To reveal the effects of changing underlying surface conditions on the stations'surrounding thermal environment during the recent 20 years,it is important to exclude the infuence of background weather conditions.To do so,LST is normalized as Iifollowing Li et al.(2015):

Where Tiis the LSTof pixel i.Tmaxand Tminare the maximum and minimum LST in the considering area.Iiis the normalized LST index for pixel i,which ranges 0-1.With the data of all the pixels in the 2 km buffer zones of the 5 stations,I is categorized into 5 zones(i.e.,high,sub-high,medium,sublow and low temperature zones)using the Natural Breaks method in ArcGIS(Table 1).

In addition,linear regression was applied to calculate the relationship between normalized LST and NDVI around the fve stations in 1990,2000,and 2010.Due to 60-m resolution of TM6,about 3000 pixels can be obtained in the 2 km buffer zones.After qualities controlled,about 2600 pixels in the 2 km buffer zones per panel are selected as samples for linear regression at each station.All correlation coeffcients were tested using a t-distribution.A p value associated with this test was applied.

2.2.4.Heat effect contribution index(HECI)

In the surrounding area of an observational station,various land use types contribute differently to the thermal environment.To quantitatively evaluate the heat contribution from each type of land use,an index(HECI)is introduced following Li et al.(2015):

Where T is the average LST of the whole area.N is the total pixel number within the area.Tijis the LST of type i LULC in its jth pixel.Niis the pixel number of the area covered by type i LULC.HECIirepresents the heat contribution LULC from type i,and it varies from 0 to 100%.Larger HECIimeans larger contribution of type i LULC to the thermal environment.The sum of HECIioverall types LULC equals 100%.

Table 1 Five zones of the thermal environment categorized by normalized LST.

3.Results

3.1.Impacts of urbanization on LULC around meteorological stations

Fig.2 shows the temporal variation of LULC in the 2 km buffer zone of the fve urban stations.Generally,it can be seen that the LULC changed signifcantly primarily due to urban expansion.Taking CZ as an example,in 1990,the construction area covered approximately 30%of the buffer zone(Fig.2a), which then increased to approximately 45%in 2000(Fig.2b), while to approximately 70%in 2010(Fig.2c).Similarly,it is clear that other four stations have been in urban areas due to the tremendous urban construction sprawl since the 1990.In fact,they were in rural areas before the 1980s.Ten years later, the construction areas of the buffer zone increased continuously,accompanying with vegetation and water reductions around these four stations(Fig.2b).Until 2010,SZ and AQ had entered into city passively,except HF and WH(Fig.2c).It is evident that the measurements can be strongly affected by surrounding buildings.Therefore,in order to avoid disturbances from urban developments,HF and WH moved to rural areas in 2004 and 2006,respectively.After their relocations, vegetation areas increased signifcantly,while construction decreased drastically.

3.2.Impacts of LULC on thermal environment around observational stations

Fig.3 shows the fve-level thermal environment based on the normalized LST in the 2 km buffer zones of the fve stations.For SZ/AQ/CZ,the high and sub-high temperature zones with the construction LULC type kept increasing signifcantly, while the area of low,sub-low and medium temperature zones, as well as the vegetation and water LULC type kept decreasing from 1990 to 2010(Figs.2 and 3).As a result,the HECI of construction for SZ/AQ/CZ also kept increasing,e.g.,65.47%/ 54.20%/48.50%in 1990,80.96%/71.37%/50.66%in 2000, and 89.81%/80.17%/61.32%in 2010,respectively(Table 2), while that of vegetation and water also kept decreasing. Similarly,increase/decrease trends for thermal environment of construction/vegetation and water are found for HF and WH during 1990-2000(Fig.3,Table 2).However,after their relocation,the high and sub-high temperature zone and the construction LULC type in the 2 km buffer zone were drastically reduced,and the low,sub-low and medium temperature zone,as well as the vegetation and water LULC type were increased.The HECI of construction for HF and WH were 69.72%and 79.04%in 2000,while decreased to 21.90%and 29.97%in 2010,respectively(Table 2).The above results clearly indicate that the observational environment of ameteorological station can be greatly infuenced by LULC under local fast urbanization.

Overall,in low and sub-low temperature zones,vegetation and water are the main contributors,while in sub-high and high temperature zones,construction becomes dominant.In the medium temperature zone,vegetation and construction are mixed.Consequently,controlling the scale and layout of the urban buildings around the stations is an effective method to ameliorate observation thermal environment and to improve region representativeness of station observation.

3.3.Relationship between normalized LST and NDVI around meteorological stations

NDVI,as an important indicator of LULC,can be derived by visible/infrared(VIR)sensors aboard most satellites with the opportunity of high sampling frequency.Therefore,it is very important to investigate the relationship between LST and NDVI to study the thermal environment around meteorological stations,which can be further used to effectively improve and populate the technique that using VIR measurements to monitor the thermal environment around meteorological stations.Previous studies indicated that LST and NDVI were signifcantly negative correlation in urban areas(Shi et al.,2011b,2013;Zhang et al.,2011).To assess quantitative relationship between NDVI and LST around meteorological stations,the linear regressions of NDVI on the normalized LST in the 2 km buffer zone of each station are shown for the year 1990,2000,and 2010,respectively (Fig.4).In general,the normalized LST and NDVI exhibit greatnegative correlations,with the normalized LST increasing while NDVI decreasing around the stations (Fig.4).The correlation coeffcients are statistically signifcant at the 99%confdence level.These results indicate that NDVI can refect the changes of thermal environment conditions in respond to LULC changes around meteorologicalstations affected by urbanization.Therefore,the variability of thermal environment could also be monitored through evaluating the variability of NDVI.

Table 2 The HECI of different LULC type around stations(unit:%).

4.Discussion

Different remote-sensing products have different advantages and disadvantages.For example,although some historical images are missing,images provided by Google Earth usually have high resolution(approximately 15 m)and can be used to evaluate the representativeness of observational stations in recent years(Ren et al.,2010;Ren and Ren,2011). Large quantities of MODIS LST products have been gathered for more than 14 years since the launch of MODIS,and they provide a unique opportunity to study the long-term representativeness of observational stations(Wang et al.,2011;Ren and Ren,2011).However,the resolution of MODIS LST products is 1 km,which is too coarse for investigating the surrounding conditions within a station's 1 km buffer zone. Moreover,because of the low spatial resolution of these products,it is diffcult to precisely construct the relationship between LULC and LST.In contrast,remote-sensing images provided by the Landsat project have both high spatial resolution(approximately 30 or 60 m)and long-term coverage(from 1979 to present).Landsat images have been widely used for land use classifcation and thermal environmental assessment around stations on a 1-km or smaller scale(Shao et al., 2011;Yang et al.,2013;Li et al.,2015).In the present study, Landsat images were reliable for constructing a relationship between LULC(or NDVI)and LST.Consequently,we found advantages in the combined LULC and LST method for station-type classifcation.

To evaluate urbanization-related warming,Fig.5 shows the time series of the differences between the observed temperatures at fve urban stations and at reference stations between 1980 and 2010.These differences can be treated as the intensities of UHI for urban stations(Liu et al.,2007;Ren et al., 2008).The fgure clearly shows that UHI values for the fve urban stations increased with time over the past 20 years.In particular,the UHI values for the SZ,CZ,and AQ stations continued to increase signifcantly after 2000;however,this was not observed for the HF and WH stations because of their relocation.This implies that urbanization-related warming cannot be ignored(Ren and Zhou,2014;Ren et al.,2014). Wang et al.(2010)also indicated that the urban heat island effect and the decadal variation in ocean thermohaline circulation were responsible for the continuing warming in China, although a global warming hiatus occurred after 2000.In addition,Stewart and Oke(2012)pointed out that the effects of terrain and local wind patterns would infuence the representativeness of the stations and the local UHI.Therefore,the exact mechanisms underlying warming in China are complex and can vary from place to place.The analysis and separation of these mechanisms are beyond the scope of this paper but will be addressed in future research.

Abrupt changes(i.e.,signifcant drops)in UHI values occurred at the HF(Fig.5b)and WH stations(Fig.5e)in 2004 and 2006,respectively,as indicated by the arrows.These detectible inhomogeneities were mainly associated with the station relocations in 2004 and 2006,which were also consistent with the LULC and thermal environment changes (Table 1,Figs.2 and 3).After relocation,the environmental representativeness improved,but the inhomogeneities in the data must be taken into account.There are methods available for adjusting some of these inhomogeneities(Yan et al.,2010; Li and Yan,2010;Cao and Yan,2012;Yang et al.,2013;Yang and Li,2014).

We selected LST rather than in situ air temperatures for the present study because of the fact that LST can be monitored using macroscopic,large-scale,dynamic,continuous satellite remote-sensing images.In contrast,in situ air temperatures,as observed by meteorological stations,are taken from a single space point and therefore have weak spatial representativeness,particularly over heterogeneous urban surfaces.In fact,air temperature and LST are closely related,and air temperature can be estimated from satellitebased LST using the temperature/vegetation index(TVX) method and a statistical regression approach(Prihodko and Goward,1997;Florio et al.,2004;Stathopoulou et al., 2006;Shen and Leptoukh,2011).Therefore,the use of LST rather than in situ air temperatures was more advantageous in the present study.

5.Summary

We used Landsat-based LULC,LST,NDVI values to study the effects of local urbanization on the thermal environments of fve national meteorological stations in the Anhui province from 1990 to 2010.The results obtained in this study are summarized as follows:

First,LULC around the observational stations changed signifcantly because of rapid local urban expansion between 1990 and 2010.Construction areas increased continuously, accompanied by reductions in vegetation and water bodies.Toavoid disturbances from these urban developments,the HF and WH stations were moved to new locations in 2004 and 2006, respectively.Therefore,in the areas surrounding these two stations,vegetation increased and construction decreased after the relocation,accompanied by improved environmental representativeness and inhomogeneity in the data.

Second,the observational environment of a meteorological station can be greatly infuenced by LULC changes.In particular,vegetation and water are the main contributors in low and sub-low temperature zones,while construction is dominant in sub-high and high temperature zones.In the medium temperature zone,the infuences of vegetation and construction are mixed.

Finally,the normalized LST and NDVI values signifcantly exhibit negative correlations,with LST increasing and NDVI decreasing around the stations.Therefore,it is reasonable to use NDVI as another indicator of thermal environmental conditions around meteorological stations.By evaluating the variability of NDVI values,the variability of the thermal environment can also be monitored.The scale and layout of construction around a station can be changed by its relocation, and this is an effective way to improve the observational thermal environment and to improve regional representativeness of station observations.In addition,when choosing a new site for a station,this study provides a method for examining the thermal environment around the intended location with high-resolution Landsat images.

Acknowledgements

This study was supported by the National Natural Science Foundation of China(41205126 and 41475085),Anhui Provincial Natural Science Foundation(1408085MKL60 and 1508085MD64)and Meteorological Research Fund of Anhui Meteorological Bureau(KM201520).We also appreciated the constructive comments and suggestions by the editors and two anonymous reviewers.

Artis,D.A.,Carnahan,W.H.,1982.Survey of emissivity variability in thermography of urban areas.Remote Sens.Environ.12(4),313-329.

Cao,L.-J.,Yan,Z.-W.,2012.Progress in research on homogenization of climate data.Adv.Clim.Change Res.3(2),59-67.

Chen,Z.-H.,Wang,H.-J.,Ren,G.-Y.,2007.Asymmetrical change of urban heat island intensity in Wuhan,China.Adv.Clim.Change Res.3(5), 282-286(in Chinese).

Florio,E.N.,Lele,S.R.,Chang,C.Y.,et al.,2004.Integrating AVHRR satellite data and NOAA ground observations to predict surface air temperature:a statistical approach.Int.J.Remote Sens.25(15),2797-2994.

Hung,T.,Uchihama,D.,Ochi,S.,2006.Assessment with satellite data of the urban heat island effects in Asian mega cities.Int.J.Appl.Earth Obs.8 (1),34-38.

Kalnay,E.,Cai,M.,2003.Impact of urbanization and land use change on climate.Nature 423(6939),528-531.

Li,Q.,Zhang,H.,Liu,X.,et al.,2004.Urban heat island effect on annual mean temperature during the last 50 years in China.Theor.Appl.Climatol. 79(3-4),165-174.

Li,Y.-B.,Shi,T.,Yang,Y.-J.,et al.,2015.Satellite-based investigation and evaluation of the observational environment of meteorological stations in Anhui province,China.Pure Appl.Geophys.172(6),1717-1733.http:// dx.doi.org/10.1007/s00024-014-1011-8.

Li,Z.,Yan,Z.,2010.Application of multiple analysis of series for homogenization(MASH)to Beijing daily temperature series 1960-2006. Adv.Atmos.Sci.27(4),777-787.http://dx.doi.org/10.1007/s00376-009-9052-0.

Li,Z.-H.,Yang,J.,Shi,C.-E.,et al.,2012.Urbanization effects on fog in China:feld research and modeling.Pure Appl.Geophys.169(5-6), 927-939.

Liu,W.,Ji,C.,Zhong,J.,et al.,2007.Temporal characteristics of the Beijing urban heat island.Theor.Appl.Clim.87(1-4),213-221.http:// dx.doi.org/10.1007/s00704-005-0192-6.

Nichol,J.E.,1994.A GIS-based approach to microclimate monitoring in Singapore's high-rise housing estates.Photogramm.Eng.Rem.S.60(10), 1225-1232.

Prihodko,L.,Goward,S.N.,1997.Estimation of air temperature from remotely sensed surface observations.Remote Sens.Environ.60(3),335-346.

Ren,G.-Y.,Zhou,Y.-Q.,2014.Urbanization effect on trends of extreme temperature indices of national stations over mainland China,1961-2008. J.Clim.27(6),2340-2360.

Ren,G.-Y.,Chu,Z.-Y.,Chen,Z.-H.,et al.,2007.Implications of temporal change in urban heat island intensity observed at Beijing and Wuhan stations.Geophys.Res.Lett.34 http://dx.doi.org/10.1029/2006GL027927. Ren,G.-Y.,Zhou,Y.-Q.,Chu,Z.-Y.,et al.,2008.Urbanization effects on observed surface air temperature trends in North China.J.Clim.21(6), 1333-1348.

Ren,G.-Y.,Zhang,A.-Y.,Chu,Z.-Y.,et al.,2010.Principles and procedure for selecting reference surface air temperature stations in China.Meteor.Sci. Technol.38(1),78-85(in Chinese).

Ren,G.-Y.,Ren,Y.-Y.,Li,Q.-X.,et al.,2014.An overview on global land surface air temperature change.Adv.Earth Sci.29(8),934-946(in Chinese).

Ren,Y.-Y.,Ren,G.-Y.,2011.A remote-sensing method of selecting reference stations for evaluating urbanization effect on surface air temperature trends.J.Clim.24(13),3179-3189.

Shao,Q.-Q.,Sun,C.-Y.,Liu,J.-Y.,et al.,2011.Impact of urban expansion on meteorological observation data and overestimation to regional air temperature in China.J.Geogr.Sci.21(6),994-1006.

Shen,S.-H.,Leptoukh,G.-G.,2011.Estimation of surface air temperature over central and eastern Eurasia from MODIS land surface temperature.Environ.Res.Lett.6(4),http://dx.doi.org/10.1088/1748-9326/6/4/045206.

Shi,C.,Roth,M.,Zhang,H.,et al.,2008.Impacts of urbanization on longterm variation of fog in Anhui province,China.Atmos.Environ.42 (11),8484-8492.

Shi,T.,Yang,Y.-J.,Jiang,Y.-L.,et al.,2011a.Impact of the variation of urban heat island intensity on temperature series in Anhui province.Clim.Environ.Res.16 (6),779-788.http://dx.doi.org/10.3878/j.issn.1006-9585.06.13(in Chinese).

Shi,T.,Yang,Y.-J.,Zhang,A.-M.,et al.,2011b.Study of thermal environment of Hefei city based on TM and GIS.Remote Sens.Technol.Appl.26(2), 156-162(in Chinese).

Shi,T.,Yang,Y.-J.,Ma,J.,et al.,2013.Study on spatial-temporal characteristics of urban heat island in Anhui province based on MODIS:typical cases study.J.Appl.Meteror.Sci.24(4),484-494(in Chinese).

Stathopoulou,M.,Cartalis,C.,Chrysoulakis,N.,2006.Using midday surface temperature to estimate cooling degree-days from NOAA-AVHRR thermal infrared data:an application for Athens,Greece.Sol.Energy 80(4), 414-422.

Stewart,I.D.,Oke,T.R.,2012.Local climate zones for urban temperature studies.Bull.Amer.Meteorol.Soc.93(12),1879-1900.

Wang,S.-W.,Wen,X.-Y.,Luo,Y.,et al.,2010.Does the global warming pause in the last decade:1999-2008?Adv.Clim.Change Res.1(1),49-54.

Wang,Y.-Y.,Li,G.-C.,Zhang,Y.,2011.Regional representativeness analysis of national reference climatological stations based on MODIS/LST product.J.Appl.Meteror.Sci.22(2),214-220(in Chinese).

Weng,Q.,2001.A remote sensing-GIS evaluation of urban expansion and its impact on surface temperature in the Zhujiang Delta,China.Int.J.Remote Sens.22(10),1999-2014.

Weng,Q.,2003.Fractal analysis of satellite-detected urban heat island effect. Photogramm.Eng.Rem.S.69(5),555-565.

Yan,Z.,Li,Z.,Li,Q.,et al.,2010.Effects of site-change and urbanisation in the Beijing temperature series 1977-2006.Int.J.Climatol.30(8), 1226-1234.

Yang,L.,2000.Integration of a numerical model and remotely sensed data to study urban/rural land surface climate processes.Comput.Geosci.26(4), 451-468.

Yang,S.,Li,Q.-X.,2014.Improvement in homogeneity analysis method and update of China precipitation data.Prog.Inquisitiones De.Mutat.Clim.10 (4),276-281(in Chinese).

Yang,Y.-J.,Shi,T.,Tang,W.-A.,et al.,2011.Study of observational environment of meteorological station based on remote sensing:cases in six stations of Anhui province.Remote Sens.Technol.Appl.26(6),791-797 (in Chinese).

Yang,Y.-J.,Wu,B.-W.,Shi,C.-E.,et al.,2013.Impacts of urbanization and station-relocation on surface air temperature series in Anhui province, China.Pure Appl.Geophys.170(11),1969-1983.

Zhang,A.-Y.,Ren,G.-Y.,Zhou,J.-X.,et al.,2010.On the urbanization effect on surface air temperature trends over China.Acta Meteor.Sin.68(6), 957-966(in Chinese).

Zhang,H.-Q.,Yang,Y.-J.,Xun,S.-P.,et al.,2011.Seasonal and spatial variability of vegetation and land surface temperature in Anhui province.J. Appl.Meteror.Sci.22(2),232-240(in Chinese).

Zhang,J.-H.,Hou,Y.-Y.,Li,G.-C.,et al.,2005.The diurnal and seasonal characteristics of urban heat island variation in Beijing city and surrounding areas and impact factors based on remote sensing satellite data. Sci.China Earth Sci.48,S220-S229.

Zhao,Z.-C.,Luo,Y.,Huang,J.-B.,2013.Are there impacts of urban heat island on future climate change?Adv.Clim.Change Res.4(2), 133-136.

Zhou,L.M.,Robert,E.D.,Tian,Y.H.,et al.,2004.Evidence for a signifcant urbanization effect on climate in China.Proc.Natl.Acad.Sci.U.S.A.101 (26),9540-9544.

Received 9 February 2015;revised 17 April 2015;accepted 17 July 2015

Available online 28 July 2015

*Corresponding author.

E-mail address:yyj1985@mail.ustc.edu.cn(YANG Y.-J.).

Peer review under responsibility of National Climate Center(China Meteorological Administration).

Production and Hosting by Elsevier on behalf of KeAi

http://dx.doi.org/10.1016/j.accre.2015.07.001

1674-9278/Copyright?2015,National Climate Center(China Meteorological Administration).Production and hosting by Elsevier B.V.on behalf of KeAi. This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).