| Home | E-Submission | Sitemap | Contact Us |  
Environ Eng Res > Volume 26(5); 2021 > Article
Kachenchart, Kamlangkla, Puttanapong, and Limsakul: Urbanization effects on surface air temperature trends in Thailand during 1970–2019


Continued urban expansion undergone in the last decades has converted many weather stations in Thailand into suburban and urban setting. Based on homogenized data during 1970–2019, therefore, this study examines urbanization effects on mean surface air temperature (Tmean) trends in Thailand. Analysis shows that urban-type stations register the strongest warming trends while rural-type stations exhibit the smallest trends. Across Thailand, annual urban-warming contribution exhibits a wide range (< 5% to 77%), probably manifesting the Urban Heat Island (UHI) differences from city to city resulting from the varied urban characteristics and climatic background. Country-wide average urban warming contribution shows a significant increasing trend of 0.15 oC per decade, accounting for 40.5% of the overall warming. This evidence indicates that urban expansion has great influence on surface warming, and the urban-warming bias contributes large fraction of rising temperature trends in Thailand. The increasing trend of annual Tmean for Thailand as a whole after adjusting urban-warming bias is brought down to the same rate as the annual global mean temperature trend, reflecting a national baseline signal driven by large-scale anthropogenic-induced climate change. Our results provide a scientific reference for policy makers and urban planners to mitigate substantial fraction of the UHI warming.

1. Introduction

Over the past several decades, the world has witnessed rapid urbanization which results in the expansion of built-up areas in and around cities, changes in artificial landscape and generation of several environmental problems [1, 2]. Urbanization strongly affects local land surface and substantially alters the urban climate, leading to the well-known UHI effects [35]. The anthropogenically driven UHI represents the differences in albedo, roughness, heat and moisture flux exchange of the land surface. It is increasingly recognized that localized UHI effects interacted with global climate change have considerable social-economic consequences and high potential for health impacts [69]. Moreover, it usually introduces a warming bias in the observed surface air temperature series, similar to that of the increasing concentration of atmospheric greenhouse gases [3, 1012]. Therefore, study of the urbanization effects on climate becomes an important component of urban climate change.
A number of studies have demonstrated that urbanization has contributed to the observed warming trends [1016]. At a global scale, evidence shows that the UHI has a trivial influence on global temperature averages due to the urban areas only comprising about 1% of the Earth’s surface [11, 1719]. However, the urbanization-induced warming becomes larger in local or national temperature observations which some locations have a similar magnitude or greater than that from greenhouse gas forced climate change [11, 12, 2022]. Over the United States, Kalnay and Cai [23] found that, in nearly one century, the annual mean surface air temperature increased by 0.27 °C owing to land-use changes and urbanization. Zhang et al. [24] estimated the UHI contribution for the whole of China to be 27% during 1961–2004. At a sub-regional scale for southeastern, eastern, and northern China, the studies showed a wide range (10% to 40%) of urbanization contribution to observed warming for different periods during 1951–2010 [10, 13, 25, 26]. Based on various approaches, Tang et al. [27] and Wang and Ge [28] showed urbanization contribution during 1961–2009 of 35%–58% for large cities. The recent study based on homogenized observations and using improved methods has showed that the urbanization effects could account for about 20% of the observed warming trend at the Beijing Observatory during the last few decades [12]. Large fraction of the urbanization contribution to the observed warming is also reported in other big cities such as Colombo, Delhi, London, Vietnam, San Salvador and Wuhan [15, 2932].
The urban minus rural (UMR) method is a direct technique widely used to examine the urbanization effects on observed surface temperature trends [11, 12, 19, 33]. It involves making comparison between the temperature trends at urban stations with those at nearby rural stations. A number of approaches which include socioeconomic indices (population, urban built-up area and growth development product) and nighttime light intensity and remote-sensed land-cover data are commonly employed to identify the urban and rural stations [11, 12, 19, 33]. In addition to the techniques used, the data quality is a comparably important factor for reliable estimation of the urbanization effects. Inhomogeneity caused by a number of non-climate factors such as changes in instruments, station moves, or different processing practices in the observed surface temperature data is often found as a major source of uncertainty in the estimation of urbanization-related warming [1113, 19].
Based on homogenized data during 1970–2019 and the UMR method, this study investigates the urbanization effects on Tmean trends for Thailand where has recently experienced fast urbanization but the scientific information on the degree to which urbanization influences local and national surface climate is not well documented. Most of the previous studies focus on UHI analysis at an individual city scale especially Bangkok Metropolis [3438]. Therefore, the main objectives are to analyze the differences in the rates of Tmean changes between the urban and suburban stations and rural stations, and to estimate the relative contribution of the urbanization effects to the trends of both individual station and country-wide average.

2. Data and Methods

This study has main steps of workflow consisting of compilation and selection of Tmean data and their quality checks, collection and analysis of land use and population data and classification of urban-type and non urban-type stations, and analysis of Tmean trends for different station types and the urbanization effects on them. The details of each main step of workflow are described in the following subsections.

2.1. Data Used and Quality Control

Monthly Tmean data from a dense observation network of Thai Meteorological Department (TMD) distributed across Thailand are used [39]. The TMD weather network was expanded in the early 1960s with the existing network having been more or less complete by the early 1970s. Data quality assurance is routinely performed, both real-time and post-processing [40, 41]. Air temperature is measured using thermometers placed in a Stevenson screen, a standardized well-ventilated shaded shelter at a height of approximately 1.2 m above the ground [40]. On the basis of record length and completeness, a set of 69 records with 50-year period (1970–2019) is initially selected to optimize data availability and spatial coverage. Overall, each of the selected records is at least 99% complete. Previous studies have demonstrated that temperature data from the TMD database are relatively good quality for climate studies [42, 43]. Nevertheless, all records are subjected to further statistical quality control checks. The objective approaches including tests of outliers, data missing interpolation and homogeneity test are employed to evaluate the quality of the selected data [44, 45]. For each record, outliers are identified by comparing their values to adjacent days and to the same day at nearby stations before they are edited or removed. Because of high correlations observed among Tmean data from nearby stations, the regression-based method of Feng et al. [44] is then applied to fill in the missing values by using data from the nearest neighboring stations.
The inhomogeneity of the data series is examined and adjusted, based on the penalized maximal t (PMT) test and the penalized maximal F (PMF) test [45, 46] using the RHtestsV4 software [47]. This technique is capable of detecting multiple step changes in time series by identifying the changepoints in a two-phase regression model [46]. In this study, inhomogeneity examination and adjustment are conducted on the monthly series. Analysis shows that 62 (89.9%) inhomogeneity records of Tmean data are detected. It is evident that most of original Tmean series underwent a systematic jump around 2006–2007. Further investigation reveals that the procedural calculation of daily Tmean was changed during that time. An absolute method which is not based on the nearby homogeneous stations as the reference series is then employed to adjust such inhomogeneity [48]. The recent period from 2007 onwards is used as the reference base to adjust the shift in the previous period, and the PMT test which assumes that the time series being tested has zero-trend and Gaussian errors [46] is applied to each independent inhomogeneous station. However, a relative technique which the candidate series is adjusted with the nearby homogeneous reference stations is applied for a few inhomogeneity records with multiple and random changes [48]. The homogeneous nearby stations well correlated with the inhomogeneous records are chosen as the reference series. All of the nearby stations used to construct the reference series have high correlations (r > 0.7) with the candidate series. The data after adjustments become more continuous and homogeneous, making they reasonably suitable for examining the urbanization effects on their trends. With this exercise, a set of 69 homogenized Tmean data spanning from 1970 to 2019 are obtained for further analysis (Fig. 1).

2.2. Classification of Stations

A key issue in examining the urbanization effects on temperature trend with the UMR method lies in the classification of urban-type and non-urban-type stations and the choice of reference stations which should be able to reflect macro-scale background climate change. In this study, we use two approaches to categorize urban-type and non-urban-type stations. Following Yang et al. [49], the differences in underlying surface and surrounding building environment of meteorological stations are first adopted to classify urban-type and rural-type stations. It is done by analyzing land use types within a radius of 2.5 km centered at each station. The analysis is based on aerial photograph-based land use data in 2012 from Land Development Department of Thailand with a scale of 1:25,000 [50]. The proportion of construction land-use type (urban and built-up) area to green area suggested by Yang et al. [49] is utilized to classify urban, suburban or rural stations (Table 1). If the stations with the proportion of construction land-use surroundings in 2012 greater than 50% with a radius of 2.5 km, they are then defined as urban-type weather stations. On the other hand, the proportion of construction land-use type area in 2012 less than 25% in 2.5 km radius buffer zone around the station is the criteria for selecting rural stations. We further analyze the LANDSAT 5/8 satellite images in 1988, 1998, 2008 and 2015 for some selected stations of each group in order to check validity of the above classification of land-use type (Figs. S1, S2 and S3). The results indicate that the proportion of urban or built-up land in 2000s for each station group is well met the criteria shown in Table 1. Based on this check, it provides additional confidence that the classification of urban, suburban, and rural stations based on land use data only in 2012 reflects the recent situations of surrounding building environment of each station group. Secondly, classification of stations on the basis of land use types is subsequently re-checked with population data. The population data in 2014 with more updated and complete version at sub-district level where each meteorological station is located are extracted from the database of Department of Provincial Administration [51]. We apply number and density of population (Table 1) as the criteria, since they are officially used to define different city sizes in Thailand as specified in the Municipality Act B.E. 2496. It is noted that the approaches used here are consistent well with field observations. Based on the above-described criteria, number of urban stations are largest (30), followed by rural stations (22) and suburban stations (17) in order (Fig. 1).

2.3. Analysis of Urbanization Effects on Observed Tmean

Following the studies of Ren et al. [13] and Yang et al. [25, 49], a pair of stations consisting of an urban/suburban station and a rural station with their distance ranging from 5 to 150 km is defined. For a given urban/suburban station, there may be more than one rural station within a radius of 50 km, and in that case the average is calculated. The urbanization effects (Tue) for an individual pair can be calculated as follows [13, 49, 52]:
where Tu/s is the observed temperature of the urban or substation station and Tr is the observed temperature of the rural station. Similarly, ΔTtre, defined as the difference in observational temperature trends (°C per decade) between urban or suburban station and rural station [13, 49, 52], can be expressed as:
Tu/stre and Trtre are the linear trends (°C per decade) for the observed temperatures of the urban or suburban stations and rural stations, respectively. The contribution of urban-related warming, Cu, to total warming is then calculated [13, 49] as:
In the trend analysis, a Mann-Kendall test which has been widely used to detect trends in meteorological and hydrological time series is applied for each station and for Thailand as a whole [5356]. The Mann-Kendal trend test is non-parametric ranked based procedure which the correlation between the rank order of the observed values and their order in time is considered [53, 56]. It is good for use with non-normally distributed variables and highly resistant to the effects of outliers [53, 56]. In this study, we used a modified version proposed by Hamed and Rao [56] which can account for time series autocorrelation by modifying the variance of the Mann-Kendal statistic S. A trend is taken to be statistically significant when the probability is at 95% level (p = 0.05). Country-wide time series for each station group is also created by averaging the anomalies of all stations to the base period from 1971 to 2000.

3. Results and Discussion

3.1. Trends of annual Tmean in Thailand

Station-by-station trends estimated based on the newly homogenized data reveal a spatially coherent and widespread significant increase of annual Tmean in Thailand during 1970–2019, consistent with the earlier study of Limjirakan and Limsakul [42] and general warming observed at regional and global scales [5762]. Analysis shows that almost all stations (98.6%) record statistically significant increase ranging from 0.05 to 0.92°C per decade (Fig. 2(a)). One of prominent features emerged from Fig. 2(a) is that the stations located at the urban and suburban areas exhibit higher increasing trends in annual Tmean compared with those situated at the rural area. These obvious differences may reflect varying degree of land use change especially urban expansion in addition to large-scale background changes. An exceptionally large increase in annual Tmean is observed at the highly urbanized Don Muang Station located in Bangkok Metropolis. Based on LANDSAT 5/8 thermal infrared data, Sanecharoen et al. [37] also found that land surface temperature in inner Bangkok city during 2008–2014 was greater 3°C higher than suburban area. Averaging all stations to present Thailand as a whole, annual Tmean has increased significantly by 0.30°C per decade during 1970–2019 (1.50°C per 50 years) (Fig. 2(b)). Using the same observed data but shorter period for the period 1970–2009, Limjirakan and Limsakul [42] reported an increase of annual Tmean in the whole of Thailand by 0.23°C per decade (0.92°C per 40 years). In comparison with the earlier results, the increasing rate of annual Tmean derived from this study is larger, probably indicating either improved inhomogeneity adjustment or a greater increase in the recent years detected based on longer quality controlled and homogenized data, or combination of both. However, our country-wide average value is in range of those values reported in Vietnam [59], Mainland China [58], and Central Asia [62].

3.2. Urban Warming Trends and Their Contribution to Overall Warming

Seasonal feature of Tmean trends in Thailand averaged each three-month season -winter (DJF), spring (MAM), summer (JJA), and autumn (SON) for three station groups, which is weak and insignificant correlations with the seasonal pattern of rainfall, is consistent with the previous studies [13, 19, 25, 30], with the largest warming appearing in wintertime (Table 2). Greater greenhouse effect resulting from larger temperature differential between the atmosphere and the sea surface in winter than other seasons may be one of possibilities causing such observed strongest warming [63]. The second largest warming occurs in autumn, followed by summertime with the least warming seen in springtime (Table 2). In the four seasons, each of the station groups experiences different warming which the strongest warming trends in Tmean all appear at urban stations, whereas rural stations show the smallest trends (Table 2). These results suggest that each of the station groups in Thailand exhibits different degree of warming and seasonally dependent, agreement with the evidence observed in China [12, 13, 19, 25, 30].
Relative contribution of urban warming to the overall warming in Thailand shows different seasonal pattern with large spatial variations (Table 2). The results show that wintertime contributes lower urban warming than springtime and summertime even though it registers the largest overall warming (Table 2). Such small contribution in winter indicates that the UHI is not a major contributor to the large wintertime warming in Thailand, in line with the results of Ren et al. [13] and Yang et al. [25] for north and east China. The weaker wintertime UHI observed in Thailand may result from the combined effect of lower solar radiation and decreased anthropogenic heat release during winter season [64, 65] thereby reducing urban to rural contrasts in energy balance [66]. On the other hand, the seasonal means of urban warming contribution generally increase from winter to spring and gradually decline thereafter (Table 2). The spring peak of urban warming contribution in Thailand is consistent with the results of Jauregui et al. [67], Ren et al. [30] and Yang et al. [25] that the UHI intensity in Mexico City, Beijing and east China is largest during the dry season. Previous studies have also shown that seasonal variations of the UHI in the tropical cities are evidently more related to urban-rural surface moisture characteristics with higher intensities during the dry season [68, 69]. Our results also agree well with the studies carried out in Thailand based on both station observations and satellite-measured data showing that UHI intensity in Bangkok, Chiang Mai and Songkhla is strongest in dry season [3438]. The maximum UHI intensity of around 6–7°C in Bangkok is found during dry season [35]. Anthropogenically generated heat in the urban area which can produce a shallow boundary layer in late-winter to mid-spring when local atmospheric circulation is relatively weak and climate conditions are clam may be possible explanation [30]. These climate conditions are developed well over Thailand during dry and hot boreal spring season which can subsequently induces enhanced UHI intensity with its impact being especially strong as observed in springtime. In addition, the studies of Roth [68] and Cui and de Foy [69] pointed out that the signal of UHI in the tropical cities was relatively weak during the rainy season. They showed that UHI intensities during this period fell dramatically as a result of strong wind and cloudiness and high rainfall amount. We believe that gradual decline of urban warming contribution observed during the rainy season (from spring to summer) (Table 2) may be caused by the effects of clouds and moisture availability, because rainfall amount and cloudiness in Thailand are substantially higher during this period.
Spatial map of annual urban warming to the overall warming in Thailand is marked by a wide range of contribution (less than 5% up to 77%) (Fig. 3). Such contribution probably reflects the UHI differences from city to city as a result of the varied urban characteristics and climatic background. In case of Bangkok Metropolis, for example, recent studies have shown that city structures, population and weather especially southwesterly wind direction and persistence are important factors governing and modulating the UHI variations [35, 38]. There are 39 (83%) out of 47 urban and suburban stations recorded significant urban warming contribution. Across Thailand, large urban warming contribution is geographically concentrated in central and eastern parts of the country (Fig. 3) where rapid urbanization and industrial and economic development have been taken place over the last few decades [70]. The first three stations with the largest urban-warming contribution are observed at Don Muang, Chonburi and Chantaburi (Fig. 3). Noticeably substantial urban-warming contribution can also be seen at some stations located in the North and Northeast (Fig. 3). When considering contribution separately for each station group, it is evident that urban stations contribute nearly two times more warming than suburban stations do, but the warming contribution from suburban stations exhibits larger variation (Fig. 3).
To further illustrate the overall picture of warming trends and warming contribution in Thailand for each station group, analysis based on country-wide time series of annual average Tmean shows different significant positive trends for all groups (Fig. 4). The largest increase occurs at the urban station group with a linear trend reaching 0.37 °C per decade; Next comes to the suburban group station with a linear trend of 0.28 °C per decade, while the smallest rise which is only 0.20°C per decade is observed at the rural station group (Fig. 4). As a consequence, the difference between urban and suburban stations combined and rural stations as urban warming contribution shows a significant upward trend of 0.15°C per decade for the last 50 years, accounting for 40.5% of the overall warming (Fig. 4). Urban warming contribution observed in Thailand as a whole is comparable to some previous studies carried out in China and the United States [10, 13, 24, 25, 27, 28, 30, 71]. Our results indicate that urban expansion has great influence on surface warming, and the urban-warming bias contributes large fraction of observed rising annual Tmean in Thailand.

3.3. Relationships between Annual Tmean Trends and Some Urbanization-Proxy Variables

Correlations between observed annual Tmean trends and some selected urbanization-proxy variables which are electricity use, Normalized Difference Vegetation Index (NDVI) and Nighttime Light (NTL) are analyzed in order to provide additional evidence illustrating possible temperature trend-urbanization relationships in Thailand. The electricity use data downloaded from the website of National Statistical Office of Thailand [72] are available only for provincial-level values which are then averaged for the period 2006–2015. Whereas, the NDVI and NTL data are district-level averaged values for the shorter period 2017–2019 obtained from the project ‘Spatial analysis of inequality in Thailand: Applications of satellite data and Geographic Information System (GIS)’ which both of data are derived from Moderate Resolution Imaging Spectroradiometer (MODIS) and Visible Near-Infrared (NVIR) bands of the Defense Meteorological Satellite Program (DMSP) Operational Linescan System (OLS) [73]. The details of how these data have been processed and calculated can be found in Puttanapong et al. [73].
Analysis shows that annual Tmean trends in Thailand have good positive correlations with provincial electricity use, NTL and NTL density (Fig. 5). Their correlations are statistical significance at the 95% confidence level. It can be seen that the relationships between annual Tmean trends and provincial electricity use and NTL are overwhelmingly influenced by the exceptionally extreme Bangkok values (Figs. 5(a) and (d). After taking the Bangkok data out, however, correlations still hold significant (Figs. 5(b) and (e)). Note that the NTL density normalized by area shows significant positive correlation (r = 0.41, p < 0.01, n = 61) with annual Tmean trends (Fig. 5(f)), while negative but insignificant correlation is observed for NDVI (Fig. 5(c)). It is known that NDVI as the most commonly used vegetation index relates to an urbanized environment especially urban land cover types [7477]. Higher NDVI values typically indicate a larger fraction of vegetation cover which in turn determines land surface temperature (LST) by the latent heat flux from surface to atmosphere through evapotranspiration [74, 75]. A negative correlation between LST and NDVI has been demonstrated by a number of studies [7477]. In other words, the NDVI values tend to be low or even negative where LST is high and vice-versa. However, NDVI measurements are subject to seasonal variations and the LST-NDVI relationship is known to be nonlinear with predominantly bare soil locations experiencing a wider variation in surface radiant temperature than densely vegetated locations [74, 78, 79]. These factors in combination with shorter period data used here may cause insignificant correlation between annual Tmean trends and NDVI as shown in Fig. 5c. Our correlation analysis generally supports the results in the previous section, highlighting that the rising trends of annual Tmean in Thailand, to a certain extent, are proportional to increasing degree of urbanization of each city in recent years as indicated by electricity use and NTL. Our results are also consistent with the recent work showing that significant positive relationships between UHI magnitude and increasing nightlight luminosity are detected worldwide [16].

3.4. Urban-warming Adjusted Annual Tmean in Thailand

Our results demonstrate that the urbanization effects on station Tmean trends in Thailand have significantly increased the country-wide average Tmean trend. Thus, trend of national baseline Tmean and detection of regional climate change and attribution to large-scale natural and anthropogenic-induced factors can better analyzed, only when the urban-warming effect is removed. As the main weather stations for climate change study in Thailand are currently situated near and within cities, further attempt is to adjust the UHI effects on seasonal and annual Tmean data.
After the adjustment of the urbanization effects, a spatially coherent and widespread increase in annual Tmean with the number of stations showing significant trends as compared to the unadjusted data remains the same (Figs. 2(a) and 6(a)). However, the warming magnitudes of the adjusted annual Tmean reduce from 0.05–0.92°C per decade to 0.04–0.35°C per decade. As a result of substantial fraction of warming caused by the UHI effects, the warming trend for Thailand as a whole is brought down from the original 0.30°C per decade to 0.19°C per decade (Fig. 6(b)). The baseline annual Tmean warming is therefore merely 0.95°C over the 50 years. It is interesting to note that the increasing trend of annual Tmean in Thailand after adjusting urban-warming bias has the same rate as annual global mean temperature trend (Fig. 6(b)), reflecting a real country-level signal driven by large-scale anthropogenic-induced climate change. Country-wide average warming in all seasons after the UHI adjustment is also significant except in spring. Winter warming decreases from 0.37°C per decade to 0.27°C per decade while spring, summer and autumn increases drop from 0.21, 0.23 and 0.33°C per decade to 0.10, 0.14 and 0.26°C per decade, respectively. Therefore, the post-adjustment seasonal temperature trends indicate that the national baseline warming magnitudes in Thailand during the 50-year periods are 1.35°C for wintertime, 0.50°C for springtime, 0.70°C for summertime and 1.30°C for autumn, respectively.

4. Conclusions

Continued alternation in land surface resulting primarily from urban expansion undergone in the last decades has converted many weather stations in Thailand into suburban or urban setting. On the basis of the newly homogenized data and the UMR method, therefore, the urbanization effects on Tmean trends and their contribution to the overall warming in Thailand during 1970–2019 are analyzed. The study highlights that widespread significant warming observed in Thailand during 1970–2019 has been affected by urban land-use expansion. From a seasonal perspective, urban-warming contribution to the overall warming in Thailand is marked by the spring peak. This seasonal feature may be due to a shallow boundary layer produced by anthropogenic-generated heat and weaker local atmospheric circulations which are well developed during dry and hot boreal spring season. Across Thailand, annual urban warming exhibits a wide range of contribution to the overall warming, probably manifesting the UHI differences from city to city as a result of the varied urban characteristics and climatic background. For Thailand as a whole, urban warming contribution shows a significant increasing trend for the last 50 years, accounting for 40.5% of the overall warming. This evidence indicates that land use change especially urban expansion has great influence on surface warming, and the urban-warming bias contributes large fraction of observed rising temperature trends in Thailand. It is also evident that, for some large cities, the significant urban-warming contribution to temperature change is comparable to that of greenhouse gas concentration, suggesting that land-surface processes can play a vital role in shaping present and future climate change in Thailand. If such trends continue, certain large city areas may experience a rate of warming well beyond the range projected by regional climate models. The increasing trend of annual Tmean for the whole of Thailand after adjusting urban-warming bias is brought down to the same rate as the annual global mean temperature trend, reflecting a national baseline signal driven by large-scale anthropogenic-induced climate change.
Our results provide a scientific reference for policy makers and urban planners to mitigate substantial fraction of the UHI warming. Nevertheless, the UHI effects in Thailand still needs further research attention. The research should be focused at least on the mechanisms leading to the evidence observed in this study, the multiple drivers of UHI and their interactions within urban domains and heat stress and its impact. Therefore, extensive remote sensing observations will be needed in future work.

Supplementary Information


We would like to thank Thai Meteorological Department, Land Development Department, Department of Provincial Administration and United States Geological Survey for kindly providing valuable temperature, land use and population data and LANDSAT imagery dataset. We also extend our thanks to the editor and three anonymous reviewers for their insightful critiques and constructive comments to substantially improve the quality of the earlier version of the manuscript.


Author Contributions

B.K. (Ph.D) analyzed temperature data and supported writing the initial version of the manuscript. C.K. (M.S. student) analyzed land use data and supported wringing the initial version of the manuscript. N.P. (Assistant Professor) analyzed Normalized Difference Vegetation Index and Nighttime Light data. A.L. (Ph.D) interpreted the results and wrote the initial version of the manuscript. All authors discussed the results and contributed to the final version of the manuscript.


1. Farrell K. The rapid urban growth triad: a new conceptual framework for examining the urban transition in developing countries. Sustainability. 2017;9:1407

2. Haase D, Güneralp B, Dahiya B, Bai X, Elmqvist T. Global urbanization. Griffith C, Maddox D, Simon D, editorsUrban Planet: Knowledge towards Sustainable Cities. Cambridge: Cambridge University Press; 2018. p. 19–44.

3. Parker DE. Urban heat island effects on estimates of observed climate change. WIREs Climate Change. 2010;1:123–133.

4. Ren G. Urbanization as a major driver of urban climate change. Adv Clim Chang Res. 2015;6:1–6.

5. Chapman S, Watson JEM, Salazar A, Thatcher M, McAlpine CA. The impact of urbanization and climate change on urban temperatures: a systematic review. Landsc Ecol. 2017;32:1921–1935.

6. Masson V. Urban surface modeling and the meso-scale impact of cities. Theor Appl Climatol. 2006;84:35–45.

7. Schlünzen KH, Bohnenstengel SI. Socio-economic Impacts-Urban Climate. Quante M, Colijn F, editorsNorth Sea Region Climate Change Assessment. Cham: Springer International Publishing; 2016. p. 417–429.

8. Leal Filho W, Echevarria Icaza L, Emanche VO, Quasem Al-Amin A. An evidence-based review of impacts, strategies and tools to mitigate urban heat islands. Int J Environ Res Public Health. 2017;14:1600

9. Heaviside C, Macintyre H, Vardoulakis S. The urban heat island: implications for health in a changing environment. Curr Environ health Rep. 2017;4:296–305.

10. Jones PD, Lister DH, Li Q. Urbanization effects in large-scale temperature records, with an emphasis on China. J Geophys Res. 2008;113:D16122

11. Wang J, Yan ZW. Urbanization-related warming in local temperature records: a review. Atmos OceanSci Lett. 2016;9:129–138.

12. Yan Z-W, Wang J, Xia J-J, Feng J-M. Review of recent studies of the climatic effects of urbanization in China. Adv Clim Chang Res. 2016;7:154–168.

13. Ren G, Zhou Y, Chu Z, et al. Urbanization effects on observed surface air temperature trends in North China. J Clim. 2008;21:1333–1348.

14. Almazroui M, Islam MN, Jones PD. Urbanization effects on the air temperature rise in Saudi Arabia. Clim Change. 2013;120:109–122.

15. Mohan M, Kandya A, Battiprolu A. Urban heat island effect over national capital region of India: a study using the temperature trends. J Environ Prot. 2011;2:465–472.

16. Paranunzio R, Ceola S, Laio F, Montanari A. Evaluating the effects of urbanization evolution on air temperature trends using nightlight satellite data. Atmosphere. 2019;10:117

17. Hansen J, Ruedy R, Sato M, Lo K. Global surface temperature change. Rev Geophys. 2010;48:RG4004

18. Hartmann DL, Klein Tank AMG, Rusticucci M, et al. Observations: atmosphere and surface. Stocker TF, Qin D, Plattner G-K, editorsClimate change 2013 the physical science basis: Working group I contribution to the fifth assessment report of the intergovernmental panel on climate change. Cambridge University Press; 2013. p. 159–254.

19. Wang F, Ge Q, Wang S, Li Q, Jones PD. A new estimation of urbanization’s contribution to the warming trend in China. J Clim. 2015;28:8923–8938.

20. Stone B. Urban and rural temperature trends in proximity to large US cities: 1951–2000. Int J Climatol. 2007;27:1801–1807.

21. Fujibe F. Detection of urban warming in recent temperature trends in Japan. Int J Climatol. 2009;29:1811–1822.

22. McCarthy MP, Best MJ, Betts RA. Climate change in cities due to global warming and urban effects. Geophys Res Letts. 2010;37:L09705

23. Kalnay E, Cai M. Impact of urbanization and land-use change on climate. Nature. 2003;423:528–531.

24. Zhang AY, Ren G, Zhou JX, Chu ZY, Ren YY, Tang GL. Urbanization effect on surface air temperature trends over China. Acta Meteor Sin. 2010;68:957–966.

25. Yang X, Hou Y, Chen B. Observed surface warming induced by urbanization in East China. J Geophys Res. 2011;116:D14113

26. Zhao P, Jones P, Cao L, et al. Trend of surface air temperature in Eastern China and associated large-scale climate variability over the last 100 years. J Clim. 2014;27:4693–4703.

27. Tang GL, Ren G, Zhou JX. Change of urban heat island intensity and its effect on surface mean air temperature records in Southwest China. J Appl Meteor Sci. 2008;19:722–730.

28. Wang F, Ge Q. Estimation of urbanization bias in observed surface temperature change in China from 1980 to 2009 using satellite land-use data. Chin Sci Bull. 2012;57:1708–1715.

29. Jones PD, Groisman PY, Coughlan M, Plummer N, Wang WC, Karl TR. Assessment of urbanization effects in time series of surface air temperature over land. Nature. 1990;347:169–172.

30. Ren GY, Chu ZY, Chen ZH, Ren YY. Implications of temporal change in urban heat island intensity observed at Beijing and Wuhan stations. Geophys Res Letts. 2007;34:L05711

31. Thanh Son N, Farn CC, Chen CR, Chen SH. Characterizing urbanization effects on land surface temperature in the capital city of El Salvador. In : EGU General Assembly Conference Abstracts; 8–13 April 2018; Vienna. p. 5664

32. Fonseka HPU, Zhang H, Sun Y, Su H, Lin H, Lin Y. Urbanization and its impacts on land surface temperature in colombo metro-politan area, Sri Lanka, from 1988 to 2016. Remote Sens. 2019;11:957–974.

33. Lamptey B. An analytical framework for estimating the urban effect on climate. Int J Climatol. 2009;30:72–88.

34. Jongtanom Y, Bualert S. Temporal variations of urban heat island intensity in three major cities, Thailand. Mod Appl Sci. 2011;5:105–110.

35. Arifwidodo S, Tanaka T. The characteristics of urban heat island in Bangkok, Thailand. Procedia Soc Behav Sci. 2015;195:423–428.

36. Ruthirako P, Darnsawasdi R, Chatupote W. Intensity and pattern of land surface temperature in Hat Yai City, Thailand. Walailak J Sci Tech. 2015;12:83–94.

37. Sanecharoen W, Nakhapakorn K, Mutchimwong A, Jirakajohnkool S, Onchang R. Assessment of urban heat island patterns in Bangkok Metropolitan area using time-series of LANDSAT thermal infrared data. Environ Nat Resour J. 2019;17:87–102.

38. Kamma J, Manomaiphiboon K, Aman N, Thongkamdee T, Chuangchote S, Bonnet S. Urban heat island analysis for Bangkok: multi-scale temporal variation, associated factors, directional dependence, and cool island condition. ScienceAsia. 2020;46:213–223.

39. Thai Meteorological Department. Mereological data [Internet]. Meteorological Telecommunication and Information Bureau; c2020. [cited 20 February 2020]. Available from: https://www.tmd.go.th/services/services.php

40. Ouprasitwong N. Indices and trends in extreme rainfall and temperature in Thailand [Technical document]. Bangkok: Meteorological Department of Thailand; 2002. p. 1–44.

41. Agsorn S. In : JMA / WMO workshop on quality management in surface, climate and upper-air observations in RA II (Asia); Tokyo: Japan Meteorological Agency; 2010; p. 1–197.

42. Limjirakan S, Limsakul A. Observed trends in surface air temperatures and their extremes in Thailand from 1970 to 2009. J Meteor Soc Japan. 2012;90:647–662.

43. Limsakul A, Goes JI. Empirical evidence for interannual and longer period variability in Thailand surface air temperatures. Atmos Res. 2008;87:89–102.

44. Feng S, Hu Q, Qian W. Quality control of daily meteorological data in China, 1951–2000: a new dataset. Int J Climatol. 2004;24:853–870.

45. Wang XL. Penalized maximal F test for detecting undocumented mean shift without trend change. J Atmos Ocean Technol. 2008;25:368–384.

46. Wang XL. Accounting for autocorrelation in detecting mean shifts in climate data series using the penalized maximal t or F test. J Appl Meteorol Climatol. 2008;47:2423–2444.

47. Zhang X. RHtestsV4 software [Internet]. Expert Team on Climate Change Detection and Indices (ETCCDI). c2020. [cited 20 January 2020]. Available from: http://etccdi.pacificclimate.org/software.shtml

48. Li-Juan C, Zhong-Wei Y. Progress in research on homogenization of climate data. Adv Clim Chang Res. 2012;3:59–67.

49. Yang Y-J, Wu B-W, Shi C-e, et al. Impacts of urbanization and station-relocation on surface air temperature series in Anhui Province, China. Pure Appl Geophys. 2013;170:1969–1983.

50. Land Development Department. Land use data in Thailand [Internet]. Information Service of the Land Development Department (LDD); c2019. [cited 10 July 2019]. Available from: http://www.ldd.go.th/ldd_en/en-US/information-service-of-the-land-development-department--ldd-/

51. Department of Provincial Administration. Population data [Internet]. Official statistics registration systems. c2019. [cited 17 November 2019]. Available from: https://stat.bora.dopa.go.th/stat/statnew/statMenu/newStat/home.php

52. Zhou Y, Ren G. Change in extreme temperature event frequency over mainland China, 1961–2008. Clim Res. 2011;50:125–139.

53. Sen PK. Estimates of the regression coefficient based on Kendall’s tau. J Am Stat Assoc. 1968;63:1379–1389.

54. Choi G, Collins D, Ren G, et al. Changes in means and extreme events of temperature and precipitation in the Asia-Pacific Network region, 1955–2007. Int J Climatol. 2009;29:1906–1925.

55. Donat MG, Alexander LV, Yang H, et al. Updated analyses of temperature and precipitation extreme indices since the beginning of the twentieth century: The HadEX2 dataset. J Geophys Res. 2013;118:2098–2118.

56. Hamed KH, Rao AR. A modified Mann-Kendall trend test for autocorrelated data. J Hydrol. 1998;204:182–196.

57. Jones PD, Lister DH, Osborn TJ, Harpham C, Salmon M, Morice CP. Hemispheric and large-scale land-surface air temperature variations: an extensive revision and an update to 2010. J Geophys Res. 2012;117:D05127

58. Ren G, Ding Y, Tang G. An overview of mainland China temperature change research. J Meteorol Res. 2017;31:3–16.

59. Nguyen DQ, Renwick J, McGregor J. Variations of surface temperature and rainfall in Vietnam from 1971 to 2010. Int J Climatol. 2014;34:249–264.

60. Babar ZA, Zhi X, Ge F, et al. Assessment of Southwest Asia surface temperature Changes: CMIP5 20th and 21st century simulations. Pakistan J Meteorol. 2016;13

61. Simmons AJ, Berrisford P, Dee DP, Hersbach H, Hirahara S, Thépaut JN. A reassessment of temperature variations and trends from global reanalyses and monthly surface climato-logical datasets. Q J R Meteorol Soc. 2017;143:101–119.

62. Peng D, Zhou T, Zhang L, Zou L. Detecting human influence on the temperature changes in Central Asia. Clim Dyn. 2019;53:4553–4568.

63. Webb M, Slingol A, Stephens G. Seasonal variations of the clear-sky greenhouse effect: the role of changes in atmospheric temperatures and humidities. Clim Dyn. 1993;9:117–129.

64. Janjai S, Masiri I, Pattarapanitchai S, Laksanaboonsong J. Mapping global solar radiation from long-term satellite data in the Tropics using an improved model. Int J Photoenergy. 2013;2013:210159

65. Poolsawat K, Tachajapong W, Prasitwattanaseree S, Wongsapai W. Electricity consumption characteristics in Thailand residential sector and its saving potential. Energy Reports. 2020;6:337–343.

66. Oleson K. Contrasts between urban and rural climate in CCSM4 CMIP5 climate change scenarios. J Clim. 2012;25:1390–1412.

67. Jauregui E, Godinez L, Cruz F. Aspects of heat-island development in Guadalajara, Mexico. Atmos Environ. 1992;26:391–396.

68. Roth M. Review of urban climate research in (sub)tropical regions. Int J Climatol. 2007;27:1859–1873.

69. Cui YY, de Foy B. Seasonal variations of the urban heat island at the surface and the near-surface and reductions due to urban vegetation in Mexico City. J Appl Meteor Climatol. 2012;51:855–868.

70. Nitivattananon V, Srinonil S. Enhancing coastal areas governance for sustainable tourism in the context of urbanization and climate change in eastern Thailand. Adv Clim Chang Res. 2019;10:47–58.

71. Hua LJ, Ma ZG, Guo WD. The impact of urbanization on air temperature across China. Theor Appl Climatol. 2008;93:179–194.

72. Office NS. Statistics for electricity use during 2006–2015 [Internet]. Tape water and electricity. c2020. [cited 12 March 2020]. Available from: http://service.nso.go.th/nso/web/statseries/statseries18.html

73. Puttanapong N, Luenam A, Jongwattanakul P. Spatial analysis of inequality in Thailand: Applications of satellite data and geographic information system (GIS). Thailand Research Fund; 2020.

74. Yuan F, Bauer ME. Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in Landsat imagery. Remote Sens Environ. 2007;106:375–386.

75. Yue W, Xu J, Tan W, Xu L. The relationship between land surface temperature and NDVI with remote sensing: application to Shanghai Landsat 7 ETM+ data. Int J Remote Sens. 2007;28:3205–3226.

76. Molnár G. Analysis of land surface temperature and NDVI distribution for Budapest using Land 7 ETM+ data. Acta Climatologica et Chorologica. 2016;49–50:49–61.

77. Guha S, Govil H. An assessment on the relationship between land surface temperature and normalized difference vegetation index. Environ Dev Sustain. 2020; https://doi.org/10.1007/s10668-020-00657-6

78. Guha S, Govil H, Diwan P. Analytical study of seasonal variability in land surface temperature with normalized difference vegetation index, normalized difference water index, normalized difference built-up index, and normalized multiband drought index. J Appl Remote Sens. 2019;13:024518

79. Owen TW, Carlson TN, Gillies RR. An assessment of satellite remotely-sensed land cover parameters in quantitatively describing the climatic effect of urbanization. Int J Remote Sens. 1998;19:1663–1681.

Fig. 1
Geographical distribution of quality controlled and homogenized monthly Tmean in Thailand during 1970–2019. The major city names where most of stations are located are also provided.
Fig. 2
Station-by-station annual trends of quality controlled and homogenized Tmean in Thailand during 1970–2019 (a) and country-wide average time series and its linear trend (b).
Fig. 3
Station-by-station annual warming contributions of urban and suburban stations to the overall warning in Thailand. Insertion is box plot of warming contribution between urban and suburban stations with central line, box edges and whiskers representing mean, 25% and 75% percentile and inter quartile range (IQR), respectively.
Fig. 4
All-Thailand averaged series and its linear trends of homogenized annual Tmean for (a) rural station group, (b) suburban station group, (c) urban station group and (d) difference between urban and suburban stations combined and rural station.
Fig. 5
Correlations between annual Tmean trends and provincial electricity use with Bangkok data (a) and without Bangkok data (b), NDVI (c), nighttime light (NTL) with Bangkok data (d) and without Bangkok data (e) and NTL density (f). Annual Tmean trends averaged at provincial level are plotted in Figures a and b, while those at individual station and district-level averages of NDVI, NTL and NTL density are used for the rest of Figures. Bangkok data are averaged values.
Fig. 6
Station-by-station trends of urban-warming adjusted annual Thailand’s Tmean during 1970–2019 (a) and its country-wide average time series (blue dash line) and annual global mean temperature (GMT) time series (black solid line) (b). The annual GMT data are global land and ocean combined temperature anomalies relative to the 1951–1980 average obtained from Goddard Institute for Space Studies (GISS) temperature analysis (GISTEMP), National Aeronautics and Space Administration (NASA). (https://climate.nasa.gov/vital-signs/global-temperature)
Table 1
Criteria used to Classify Urban, Suburban and Rural Stations
Type of stations Land use-based Population-based

Proportion of construction land-use type area to green area Number Density per area Size of municipality
Urban > 50% > 20,000 > 3,000 Large municipality
Suburban 25–50% 10,000–20,000 1,000 – 3,000 Medium municipality
Rural < 25% < 10,000 < 1,000 Small municipality
Table 2
Averages and Standard Deviations of seasonal and Annual Thailand’s Tmean Trends for Each Station Group and Urban Warming Contribution during 1970–2019. The Number in the Parentheses Denotes Stations with Significant at 5% Level
Annual Winter Spring Summer Autumn
Urban (°C per decade) 0.37 ± 0.16 (30) 0.47 ± 0.17 (29) 0.31 ± 0.18 (25) 0.32 ± 0.17 (29) 0.42 ± 0.16 (30)
Suburban (°C per decade) 0.28 ± 0.09 (17) 0.34 ± 0.19 (12) 0.15 ± 0.13 (12) 0.19 ± 0.12 (14) 0.29 ± 0.18 (15)
Rural (°C per decade) 0.20 ± 0.11 (21) 0.27 ± 0.17 (19) 0.12 ± 0.16 (12) 0.16 ± 0.15 (18) 0.25 ± 0.16 (20)
Urban warming contribution (%) 39.3 ± 22.8 (38) 36.7 ± 22.0 (41) 60.8 ± 30.7 (38) 43.8 ± 29.9 (35) 31.3 ± 26.0 (36)
PDF Links  PDF Links
PubReader  PubReader
Full text via DOI  Full text via DOI
Download Citation  Download Citation
Supplement  Supplement
Editorial Office
464 Cheongpa-ro, #726, Jung-gu, Seoul 04510, Republic of Korea
TEL : +82-2-383-9697   FAX : +82-2-383-9654   E-mail : eer@kosenv.or.kr

Copyright© Korean Society of Environmental Engineers. All rights reserved.        Developed in M2community
About |  Browse Articles |  Current Issue |  For Authors and Reviewers