MODIS Summer SUHI Cross-sections Anomalies over the Megacities of the Monsoon Asia Region and Global Trends

RESEARCH ARTICLE MODIS Summer SUHI Cross-sections Anomalies over the Megacities of the Monsoon Asia Region and Global Trends Hofit Itzhak-Ben-Shalom, Pinhas Alpert, Oded Potchter and Rana Samuels The Porter School of Environmental Studies, Tel Aviv University, 69978 Tel Aviv, Israel Department of Geosciences, Tel-Aviv University, P.O. Box 39040, Tel Aviv, Israel Department of Geography and the Human Environment, Tel Aviv University, Tel Aviv, Israel Department of Geography, Beit Berl College, Beit Berl, Israel


INTRODUCTION
50% of the world's population now lives in cities, with this percentage projected to increase [1].Urban areas with populations exceeding 10 million inhabitants are known as megacities; and such large concentrated populations have been identified as a major environmental issue, especially in developing countries [2].Although cities occupy a small fraction of land, amounting to only 2.8 per cent of the Earth's land area [3], they are considered to have indirectly triggered global climate change due to their insatiable appetite for energy and materials [4].Alpert et al. [5], for instance, found that in megacities, which tend to have high levels of aerosol emissions, the Aerosol Optical Depth (AOD) indicates significantly increasing levels of anthropogenic pollution.Furthermore, rapidly expanding cities, as major consumers of energy and materials, have deteriorating air quality, as population, traffic, energy use, and industrialization grow [2,6,7].Fossil fuels are consumed due to large vehicular populations.Given their very high population densities and high-rise buildings, cities also absorb a greater amount of solar radiation.Seven of the 10 most populous countries in the world − China, India, Indonesia, Russia, Pakistan, Bangladesh and Japan − are located in Asia [8].
The Asian region is marked by continuing high-level emissions of greenhouse gases from human activities [9].IPCC Global Climate Models (GCMs) predict that the increase in annual mean warming could reach about 3°C by the 2050s and about 5°C by the 2080s over selected regions in Asia [4].A recent study using a Regional Climate Model (RCM) found that by the end of this century Southwest Asia will be uninhabitable [10,11].
In light of the burgeoning population, social environment scientists and policy makers have recently begun to take special notice of urban climate change [12,13].Various studies have also argued that significant collaboration between city planners and climate change scientists is required if these areas are to continue to maintain their livability and reduce their effect on global climate change [14 -19].
Until recently, research on urban warming and its links to regional as well as global warming has been limited to observational measurements taken in mature cities, the only sites providing good sources of the data required for interpreting trends in urban climate change [20].These in situ data have the advantage of long-term records but the disadvantage of poor spatial resolution [21 -24].More recently, observation from satellites, one of the most novel and effective monitoring tools developed for recording temperature change phenomena [25] over broad geographic areas has come into use due to its methodological efficacy in recording local information [4].
Land surface temperature (LST) plays a significant feature of global temperature [20].The field of remote sensing lends itself well to urban climatology studies due to its greatly improved spectral and spatial resolution, which facilitates tracking the development of the respective phenomena [26].The advantage it provides in measuring LST as a feature of global temperature has been noted in several studies that have focused on individual urban areas in China [27 -30], India [31], and Dubai [32].
Recent studies have linked change in LST with atmospheric circulation; see for example, Zhuo [33] in China and Zhou and Wang [34] over global deserts.Lensky and Dayan [35] attempted to quantify the impact of several atmospheric circulations on LST patterns over the Eastern Mediterranean for 2000-2012 MODIS data.They found that the differences between the LST data retrieved from satellites and 2-m air temperature data as reported by land-based meteorological stations generally increase with stronger sun insolation but decline with the turbulence caused by strong winds.These parameters also vary according to synoptic-scale circulation, which is affected by horizontal pressure gradients and cloud cover.Lensky and Dayan [35] have shown that synoptic circulation induces LST patterns, with variations in LST related, among others, to vegetation cover, clouds, water vapor, surface heterogeneity, and topography [36].In another study, Lensky and Dayan [37] found that the high spatial variability of surface temperatures induces similar variability in 2-m air temperatures, hampering representation of these temperatures in numerical models [38].As LST differs from 2-m air temperature in terms of physical magnitude, the derivation of 2-m air temperatures from LSTs is not straightforward [39].
In consequence, methods applying statistical approaches and spatial smoothing are employed when predicting air temperatures from LSTs even though the outcomes are applicable only to the domain for which they were developed [40,41] and are sensitive to station density [42].Although Blandford [43] adopted a physical approach using surface data at different altitudes to retrieve the thermal profile of the ~1500 m boundary layer, this method is unable to detect near-surface (2-m) temperatures.
In the last 2-3 decades, the number of comparative studies using satellite data has grown, with the technology allowing coverage of large geographic areas, thereby creating new opportunities for climate change research [20].Imhoff [44], for example, examined the LST data obtained from MODIS for an average of three years (2003)(2004)(2005) in the 38 most populous cities in the United States.The results demonstrated that Urban Heat Island (UHI) amplitude increases with city size and is seasonally asymmetric at a large number of sites.Peng [45] selected 428 of the world's largest cities when exploring the global drivers of the Surface Urban Heat Island (SUHI).Jin and Dickinson [39], analyzed the MODIS skin temperature observations during 9 years (2000-2009) for diurnal, seasonal, and inter-annual variations at a 0.05• latitude/longitude grid over the global land surface.In 2010, Clinton and Gong [14] performed a composite Aqua and Terra -MODIS analysis (with 5-and 10-km buffers) of the factors influencing urban vs. rural SUHI temperatures, followed by its global analysis and ranking of the urban areas surveyed.Both studies mapped urban areas while here the focus is given to the individual megacities in Asia.Also, extending time series longer-term LST trends (2002-2012) as they appeared in all the world's megacities are investigated here, in an effort to explain the intermegacity variations.Tran et al. [20] also employed a 1-km MODIS resolution for 2001-2003 in their study of the SUHI's spatial patterns for eight Asian megacities as well as its relationship to land-cover properties such as vegetation, population and city area.While here, the 1-km MODIS resolution was employed when investigating longer-term LST trends (2003-2014) over the world's 29 largest megacities (>10 million inhabitants).It should be noticed that there are many more studies of the UHI of individual cities that are not cited here.
The overall objective of this study is to provide local comparative cross-sectional data of the SUHI for the top 10 megacities in the Monsoon Asian Region, the world's most rapidly growing geographic area.A further goal, for global comparison, was to provide an overview of trends in global temperature changes for 29 megacities based on MODIS LST data retrievals.In addition, comparisons between LST trends with different MODIS resolutions, i.e., 1 deg and 1 km, available in the Monsoon Asian Region, were also conducted.

Population Data
For the purpose of this study, we defined global megacities as such with populations exceeding 10 million, based on Brinkhoff's recommendations [46].Agglomerations include a central city and neighboring communities, linked to the central city by either continuous built-up areas or a predetermined number of commuters.Some agglomerations may have two or more central cities (e.g., The Ruhr) although the majority have only one.Countries, cities, population and latitude-longitude coordinates appear in (Table 1 and Fig. 1).

MODIS Data
MODIS, the Moderate Resolution Imaging Spectroradiometer, is a satellite programmed to collect data for the spatial and temporal characteristics of the global LST during daytime and nighttime [47].The monitor has 36 channels, spanning the spectral range of 0.41 to 15 µm.MODIS, with its 2330-km viewing swath, provides almost daily global coverage.We employed Collection 5 (MOD08_M3.050) of Terra -MODIS and Collection 5.1 of Aqua -MODIS (MYD08_M3.051)level-3 monthly data with a global 1°×1° grid for the summer months -June to August -for the tenyear period 2002-2012.
An Aqua and Terra -MODIS database was constructed to enable estimation of LST trends.Global and local overviews in specific megacities were obtained using the two aforementioned sensors, positioned aboard two NASA satellites: Terra (launched in December 1999) and Aqua (launched in May 2002).With cameras focused directly (straight down) at the earth's surface, MODIS technology records a broad range of global temperature and other data.The LST trends themselves were calculated for the summer months, when data from the two sensors were available [48 -51].

Cloudiness Effect
The MODIS expanse becomes limited with increasing cloud cover [52,53] due to the aerosol effect.This means that the satellite's data retrieval under overcast conditions is less accurate than when cloud cover is nonexistent or thin due to aerosols.According to Remer et al. [52,53], percent cloud cover exceeding 0.8, may lead to overestimated satellite retrieval of AOD because cloud droplets can become mixed with coarse mode particles.In order to minimize retrieval uncertainty, the data used was collected only during months characterized by percent cloud cover below 0.7.The effect of cloudiness on land surface temperature may be significantly smaller (personal communication with MODIS people) still the present analysis screened high-cloudiness months.It is assumed that this screening did not have a significant influence on our SUHI results.

The MODIS 8-Day 1-km LST
Products in the Giovanni cover the Monsoon Asian Integrated Regional Study (MAIRS) region [50].The 8-day product's observation time is a simple average of daily observation time.Terra's daily local equator-crossing times are 10:30 AM and 10:30 PM, whereas Aqua's daily passes occur at local 1:30 AM and 1:30 PM.Orbit width is 2,330 km; hence, local time at its edge indicates a difference of approximately 1.5 hours.For the SUHI investigation we focused only on Aqua -MODIS data to ensure that the satellite sensors could differentiate between LST trends over megacities and the surrounding rural areas at 1:30 PM, the hottest hour of the day.The subsequent data file, which covers the entire Asian Monsoon Region, was comprised of averaged daytime and nighttime LSTs under cloudless conditions, with quality flags, during the summer months June to August, 2003-2014.
It should be clarified that the Terra and Aqua -MODIS readings obtained to extract trends and mean LST data have two different resolutions; first, 1°×1° (summer 2002-2012), for the world's 233 largest cities and second employing the finer 1 km horizontal grid (summer 2002-2014), for the 10 megacities in the Asian Monsoon Region.In addition, it was assumed here that the detectable temperature differences and/or absolute temperature accuracies of the MODIS sensor had a small effect on the calculated trends.This, because changes in the MODIS sensor accuracy during the 11-12 years period were assumed to be small enough.

Temperature Analysis
Our approach to estimating the effect of urbanization on temperature change across the world's top ten megacities was based on an analysis of long-term variations in LST trends for the chosen sites in the selected time period.The LST parameters estimated in three stages.First, we investigated the latitudinal and longitudinal spread of each city center and confirmed our results with Google Earth.In order to compare mean LST distributions between city centers, mean LST values, labeled "running mean", acted as a filter for the calculated vectors as follows in Equation ( 1): (1) Y (i) = sum(X (j)) / (2*M+1), for j = (i-M): (i+M), and i=1: length(X) Where Y = RUNMEAN (X, M) calculates a running mean on the elements of vector X, which in this study is the LST Average, using a window of 2*M+1 data points.M is defined as the length of each city's vector for each city's latitudinal and longitude distribution.
In the second stage, we created a global map at 1°×1°, day and night, covering the world's 233 largest cities (i.e., population exceeding 2 million), including the world's 29 megacities, based on the Terra and Aqua -MODIS readings obtained during the study period (summer 2002-2012).To ensure that the satellite sensors would differentiate between LST trends in megacities and the surrounding rural areas, we used Aqua -MODIS (1 km resolution) data exclusively.
The 1 km resolution is better for understanding the impact of urbanization on local climate.The Aqua data are available for 10 megacities in the Asian Monsoon Region (counted among the world's 20 largest cities) during June-August, 2003-2014.It should note that cloudiness impedes data collection in tropical climates during the rainy season.Hence, no data is available for India among other sites.In addition, megacities for which more than six years of data were missing were deleted from the sample in order to minimize retrieval uncertainty.

Study Area
The Asian Monsoon Region extends from 0 o to 60 o N, 60 o to 150 o E; it includes eastern China, the Japanese islands, the Korean peninsula, Mongolia, Taiwan (Republic of China), Pakistan and India.Geographically, the region is located in the eastern segment of the Eurasian continent -the world's largest continent -and borders the Pacific, the world's largest ocean [4]; (Fig. 1 and Table 2).

The Surface Urban Heat Island 1 km Resolution Cross-Sections
First, the cross-sections of LST for 10 selected megacities according to the latitude (E-W) as shown in Figs.(2, 4) and longitude (N-S) Figs.(3,5) of their SUHI averaged for 2003-2014 period, were investigated.The exact location of the city center was taken from Brinkhoff [46] and confirmed by Google Earth.(2, 3) compare the Lat and Lon distributions of the SUHI Aqua -MODIS LSTs in the proximity of the city centers.They show that the LSTs are on the average decreasing with distance with an approximated bell-shaped distribution for both day and night, with a pronounced SUHI maximum over the city center.However, it should be noted that our figures are averages over a period of 12-y during which the city fringes have probably undergone significant LST changes.Hence the LST decreases with distance from city center represent a 12-y average picture which misses sub-period variations.
It is also interesting to note that during daytime, pronounced maximum LST are found especially over the megacities of Tokyo, Seoul and Beijing with a larger decrease for latitudinal (E-W) variation as shown in Fig. (2a) as compared with the longitudinal (N-S) change in (Fig. 3a).
During nighttime as shown in Figs.(2b, 3b), the LST slopes were found larger in the coldest cities, i.e., Tokyo,

Individual LST Cross-Sections
The "Urban minus Rural temperature" differences, ΔT u-r LST Aqua -MODIS data is calculated between the city center and the nearby urban area for each megacity, in terms of latitudinal as in Fig.In all the four graphs Figs.
(2, 5), it can be identified the location where the SUHI begins, ends and sometimes starts again.For example, Figs. (2,3) shows the SUHI appears very clearly in Tokyo, Beijing and Osaka during the daytime, while at nighttime the SUHI appears clearly in the tropical megacities of Kolkata and Seoul.
The next section presents maps of the day and night LST trends with time for all global megacities during the periods 2002-2012 and 2003-2014 for Terra and Aqua, respectively.

Maps of Global Megacities LST Trends with Time
In order to obtain a global map of day and night LST trends based on Terra and Aqua -MODIS temperatures that are covering the largest cities similar to the methodology employed by Alpert [5] for aerosol trends, were employed.However, here the trends were calculated by the LST time-slopes.The maps in Figs.Based on Terra -MODIS data from the 29 megacities (population over 10 million) in the sample, in Fig. (6a) shows a significant LST increase in Asian megacities, i.e., Cairo, Jakarta, Istanbul, Seoul, Tehran, Osaka and especially in Moscow.Increasing LST is also found in Mexico City and New York.In contrast, London and Paris, Los Angeles, São Paulo, Rio de Janeiro and Buenos Aires, and some Asian areas (Tokyo, Beijing, Tianjin and Karachi), exhibit declining LSTs.Out of the 29 megacities Terra and Aqua -MODIS (1°×1°) show similar increasing trends in Cairo, Istanbul, Osaka and Moscow.Similar declining LSTs were found for both Terra/Aqua in London, Paris, Los Angeles, São Paulo, Rio de Janeiro, Buenos Aires and Tokyo (Fig. 6 and Table 2).The 10 megacities LST trends in Fig. (8) and Table (3) (based on Aqua -MODIS) are more pronounced due to the 1 km resolution as compared with the previous maps at 1°×1° in Figs.(6,7).The Aqua -MODIS results for the Monsoon Asian Region clearly show dominant increasing LSTs during daytime.Furthermore, declining nighttime LST was indicated in the three largest megacities, i.e., Tokyo, Guangzhou and Shanghai.The Daily Temperature Range (DTR) was found to increase except for Osaka.In addition, increases in the LST throughout the day and night were found in Seoul, Karachi, Osaka and Beijing.Our statistical estimates for the slopes during daytime hours were found to be significant at the 95% confidence level (p<0.05)only for Karachi (p<0.03).Note that the Karachi data are for 11 years only, ending in 2011; (Table 3).

DISCUSSION
According to the report issued by the UN's Department of Economic and Social Affairs [1], the world's human population is expected to exceed 10 billion in the near future.Countries with high population densities will consequently witness rising pressure on the environment.The pressure in the Asian region is likely to increase even more when compared to other countries due to the rising demand for grains and economic development [9].
Since the beginning of the 21st century, most studies have examined climate change with in-situ data, usually collected in one or two cities [20,51].Here, in contrast, the use of latitudinal and longitudinal cross-sectional Aqua -MODIS LST data enabled mapping of 10 Asian megacities during June-August.Based on these data, we found that the sites exhibiting maximum SUHI intensity roughly coincided with the heavily built-up areas located in major city centers.Hence, cities with extremely urbanized areas, characterized by low albedo and dense construction, exhibited significantly higher daytime surface temperatures when compared to the surrounding rural area comprised of relatively moister vegetated areas [20].With respect to the LST magnitudes shown in Figs. ( [20] study, which produced a similar ordering of cities by temperature gradient, i.e., Tokyo, followed by Beijing, Seoul and Shanghai.Temperatures were also found to decrease outward from the city center, distributed according to an approximate bell-shaped curve during the day, a finding that can be explained by the UHI phenomenon.During the night, megacities such as Guangzhou and Shanghai showed far gentler SUHI, which can be explained by the effect of peripheral cities or industrial centers, a condition that influences LST.Another factor could be the sun insolation that plays an important contribution to the increase of daily SUHI.Tran et al. [21] also suggested that SUHI diurnal variations are due to differences in urban-rural characteristics such as moisture availability and surface cooling rates, further affected by topography, rural surroundings and the geographic location of cities relative to the sea [20,52]. In contrast to in-situ UHI studies, in which minimum daily temperatures tend to increase at a higher rate than do daily temperatures (the latter inducing a decreasing long-term diurnal temperature range; see for example [52,53]), the satellite-basedF SUHIs temperatures observed tending to be higher during the day, with peak temperatures reached in areas of large buildings or paved surfaces, and lowest during the night, due to difference in surface properties and cooling rates [20,51].Therefore, in the current study, in order to monitor the seasonal variations in SUHIs (at a 1 km resolution), we focused only on Aqua -MODIS rather than Terra data, thus ensuring that the satellite sensors could differentiate between LST trends observed in megacities and those observed in the surrounding rural areas, both at 13:30 local time (LT), frequently observed to be the hottest hour of the day (as opposed to the 10:30 LT Terra measurement).Here, the DTR increased during 2002-2014 Table (3), which fit the earlier findings reported by Roth [51] and Tran [20], among others.Throughout the research period, a negative DTR was found only in Osaka.Furthermore, according to the IPCC's Atmosphere-Ocean Global Climate Model (OGCM) simulations the DTR will increase over central Asia and be significantly higher there than in other regions during the summer months of the 2050s and 2080s [54 -56].
We should stress that, as Roth [57] has already noted, the limited range of urban morphologies characterizing the cities, as well as the role of building density and materials, prevented detailed exploration of temperature trends.Hence, further study is required in order to explore the contribution of these and other parameters, such as population density, energy consumption, and transportation volume, to urban warming.

CONCLUSION
This article has stressed the advantages provided by satellite remote sensing monitoring techniques when performing quantitative climatic observations worldwide, an ability that has significantly enhanced our ability to identify long-term trends in climate change [41].These techniques have made reliable high-quality global coverage at high spatial resolution available, thereby creating opportunities for researchers to easily obtain improved global meteorological data and thus providing significant advantage relative to ground-based measurements.We should recall, however, that unlike ground-based measurements, satellite data are limited regarding the parameters to be measured; they likewise cannot provide accurate information on climate variables, which are important for the construction of human comfort indices.
The 1 km LST data from Aqua -MODIS can distinguish urban from surrounding rural areas, information supporting statistical analyses by making it possible to create bell-shaped data distributions (for latitude as well longitude magnitudes), which can indicate changes of trends in LST.By doing so, satellites provide us with the opportunity to compare global-temperature trends in different cities on the basis of data obtained from the same source and with the same measurement technology, irrespective of city size.
The same Aqua -MODIS (1 km) satellite source indicates that for Asia's three largest megacities --Tokyo, Guangzhou and Shanghai --a recent decrease in nighttime LST trends (2002-2014) has been observed.In contrast, the LST slope for both day and night was found to increase in Seoul, Karachi, Osaka and Beijing.Further research on the physical, economical, and social environment is required in order to better understand and evaluate climate changeinduced vulnerabilities in the Asian Monsoon Region.These data are also relevant for assessing the feasibility of introducing adaptation measures in other regions.We should also anticipate that given the vital roles played by surface water and groundwater resources in forestry, agriculture, fisheries, livestock production and industrial activity in the Asian countries, the agricultural sector will continue to be highly sensitive to climate change impacts in that region [4].
As to the Terra -MODIS and Aqua -MODIS (1°×1°) LST trends, Alpert et al. [5] have suggested that in cases where the different sensors show similar trends, the results can be considered reliable.This conclusion is supported by our findings for daytime hours, which were observed as increasing in the Mid-East and the Asian Monsoon region and in the southeastern US during 2002-2012.In contrast, trends in declining LST appeared to be dominant in Europe sites, Saudi Arabia, parts of Northwest America, and the Eastern segment part of South America.One suggestion is that those cities are aware of climate change, therefore implement climate change issues into their policies.
In general, we concluded that the MODIS satellite can be better-utilized if we were to adapt our current technologies and devise new methodologies.While most climatological studies conducted to date have employed in-situ data, the current study incorporated satellite data having global spatial coverage despite being limited in the variables to be observed.Furthermore, despite the methodological limitations of the current study, we suggest that the optimal direction of future research to involve: (a) Obtaining data with higher resolution MODIS 8-day 1 km LST for additional regions; (b) Expanding the list of variables, e.g., air humidity and wind speed, as tracked from land-based meteorological stations; and (c) Applying further combinations of high resolution satellite data with Global Climate Model (GCM) data for evaluation.

ETHICS APPROVAL AND CONSENT TO PARTICIPATE
Not applicable.

HUMAN AND ANIMAL RIGHTS
No Animals/Humans were used for studies that are base of this research.

CONSENT FOR PUBLICATION
Not applicable.

Fig. ( 1 ).
Fig. (1).Global distribution of population (2010), over the world's 233 largest cities (population exceeding 2 million).The largest cities in the Monsoon Asian Integrated Region are in bold.Population magnitude is designated by circles of different diameters and colors, as elaborated in the bottom panel.

Fig. ( 2 ).
Fig. (2).Latitudinal distribution of 12-year (2003-2014) mean LST for (a) daytime and (b) nighttime, based on 1 km Aqua -MODIS.The 10 Asian megacities are listed on the panel to the right.Data are limited in time-span for Karachi (up to 2011), Dhaka (up to 2013), Delhi (up to 2008) and Kolkata (nighttime only).

Fig ( 4 )
Fig (4).The change in Urban-Rural of the longitudinal distribution of 12-year mean (a) daytime and (b) nighttime LST based on Aqua -MODIS data sets for 1 km for each megacity.

Figs.
Figs.(2,3) compare the Lat and Lon distributions of the SUHI Aqua -MODIS LSTs in the proximity of the city centers.They show that the LSTs are on the average decreasing with distance with an approximated bell-shaped distribution for both day and night, with a pronounced SUHI maximum over the city center.However, it should be noted that our figures are averages over a period of 12-y during which the city fringes have probably undergone significant LST changes.Hence the LST decreases with distance from city center represent a 12-y average picture which misses sub-period variations.
Latitudinal distribution of 12-year mean during daytime (b) Latitudinal distribution of 12-year mean during nighttime Seoul, Osaka and Beijing (the cities with lower temperatures).In addition, a relatively flat LST distribution along both latitude and longitude was found in Guangzhou, Shanghai, Delhi and Dhaka; it should be noted, however, that SUHI can still reach ~2-4 o C (as particularly noticed in the individual cities' cross-sections in Figs.(4, 5).In Shanghai, the flat distribution at night is contrasted during the daytime by a pronounced bell-shaped distribution in both latitude/ longitude of ~8-11 o C.

Fig ( 5 )
Fig (5).The change in Urban-Rural of the longitudinal distribution of 12-year mean (a) daytime and (b) nighttime LST based on Aqua -MODIS data sets for 1 km for each megacity.

( 4 )
and longitudinal in Fig. (5) for (a) day and (b) night.In Fig. (4), a pronounced distribution was found in Tokyo, Seoul, Osaka and Beijing (arranged in the decreasing order) with a maximum of 15 o C, following with Shanghai (12 o C); Guangzhou (8.8 o C) and Karachi with 5.1 o C SUHI.Fig. (4b) represents the latitude nighttime, the maximum SUHI was found with smaller values (compared to the daytime), in Beijing (7.6 o C); Tokyo (6.4 o C), Osaka and Seoul (~5.5 o C); Guangzhou and Dhaka (~3.5 o C).In the other megacities, ΔT u-r reaches maximum values of 2.5 o C SUHI only.According to Fig. (5a), for the longitude daytime data, the maximum SUHI found in two cities Osaka and Beijing in the order of 15 o C. Following are Seoul, Shanghai and Guangzhou with ~11 o C; and Karachi with 6.5 o C SUHI.Fig. (5b) represents the longitude nighttime data; the maximum temperature differences were also found to be more moderate compared to the daytime as in Fig.(4).For example, the maximum SUHI was found in Tokyo(6

( 6 , 7 )
represent the world's 233 largest cities having populations exceeding 2 million [46].Circles of various diameters and colors indicate the extent and sign of the LST tendencies (see bottom panel).The circles radii represent the value of the LST trends, while the color indicates the sign of the trend, blue/red indicates decrease/increase.

Table 2 . The world's top 29 megacities, listed by number of inhabitants in decreasing order (see Table 1). The 2002-2012 LST trends ( o C/11y) are for the two MODIS sensors, Terra and Aqua, at 1°×1°; NA = not available. Statistically significant slopes at a 95% (i.e., P < 0.05). In parentheses, the (R 2 ; P -value) are also indicated.
The top 10 megacities in the Monsoon Asian Integrated Region are indicated in bold.

Table 3 . The 10 megacities in the Asian Monsoon Region; LST slope (°C/12y) during June-August 2003-2014, based on 1 km Aqua -MODIS. Statistically significant slopes at a 95% (i.e. P < 0.05) are in bold. NA = not available.
Data for Karachi, Dhaka, Delhi and Kolkata are limited as indicated earlier.