Next Article in Journal
A Community-Based Sensor Network for Monitoring the Air Quality in Urban Romania
Next Article in Special Issue
Individual and Combined Effects of 3D Buildings and Green Spaces on the Urban Thermal Environment: A Case Study in Jinan, China
Previous Article in Journal
Microplastic Pollution Research Based on the VOS Viewer Software: Research Trends, Ecological Effects, and Testing Methods
Previous Article in Special Issue
Evolution over Time of Urban Thermal Conditions of a City Immersed in a Basin Geography and Mitigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Urban Heat Island Effect, Heat Stress and Public Health in Colombo, Sri Lanka and Shenzhen, China

1
State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
State Environmental Protection Scientific Observation and Research Station for Ecology and Environment of Rapid Urbanization Region, Shenzhen Environmental Monitoring Center, Shenzhen 518049, China
4
School of Life Sciences, Division of Life Sciences and Medicine, University of Science and Technology of China, Hefei 230027, China
*
Author to whom correspondence should be addressed.
Submission received: 25 March 2023 / Revised: 26 April 2023 / Accepted: 28 April 2023 / Published: 8 May 2023

Abstract

:
Human health, energy and comfort are determined by the climate that remains in the physical environment. Regarding urban climate, few studies assess the urban heat island effect, heat stress, and public health as geographical representations. This study seeks to fill this gap by selecting Colombo, Sri Lanka, and Shenzhen, China, comparatively, two coastal cities with different climate conditions. We quantified and compared the effects of heat waves and their impacts on public health and the effect of urbanization on urban heat islands (UHI). Heat-related public health issues have been calculated using the Wet-Bulb Globe Temperature (WBGT) index. The Urban Heat Island (UHI) effect was analyzed using Land Surface Temperature (LST), created based on Landsat images obtained in 1997, 2009 and 2019. A rapid increase in temperature and humidity creates an uncomfortable environment in both cities, but apparent differences can be observed in climatic phenomena. During the summer (June to August), the prevailing atmospheric condition in Shenzhen makes a “Very severe stress” with Heatstroke highly likely. Nevertheless, seven months (November to April) are found as “Comfortable” without having any heat-related health injuries. However, Colombo has never been classified as “Comfortable” throughout the year. Out of twelve, five months (April to August) are found as “Very severe stress” with Heatstroke highly likely. When considering the urban expansion and UHI, a fast expansion can be observed in Colombo than in Shenzhen. Consequently, with the more severe heat-related public health and rapid urban heat island expansion, Colombo makes it more stressful than Shenzhen city. Our findings highlight the comparison between heat-related public health and urban heat island between two coastal cities with different climate conditions and under rapid urbanization processes. Therefore, it is imperative to assess these risks and respond effectively.

1. Introduction

The climate is one of the physical sources that is extremely affected not only by thinking of human beings but also by human behaviour and activities [1]. Human health, energy and comfort are determined by climate remaining in the physical environment [2,3,4]. Climate change refers to significant, long-term changes in the global climate [5]. It directly influences global warming and is projected to affect human health, with primarily negative consequences of increasing the number of excess deaths and hospital admission worldwide [6,7,8]. Urbanization results in remarked removal and replacement of vegetation cover by various built-up structures, which leads to the urban heat island (UHI) Effect [9,10]. Climate change and urbanization combine to create a perfect disaster of global health risks, such as general discomfort, respiratory issues, heat cramps and strokes, exhaustion, and even heat-related mortality [11,12,13]. Most of these adverse environmental effects can be minimized by identifying problems and implementing proper urban planning systems with sustainable solutions [14,15,16]. This study examines the heat waves and their impacts on public health and the effect of urbanization on UHI by selecting the highly urbanized cities of Colombo, Sri Lanka, and Shenzhen, China—comparatively, two coastal cities with different climate conditions.
China is the largest developing country and has a faster increase in surface air temperature than the global average [17,18]. Although the Chinese Government has made increasing efforts to address the impact of climate change, it appears that the health implications have received less attention compared with the energy, economic and agriculture sector [19]. People may adapt to heat stress through modifications in activities, increased use of air conditioning, and alternative building designs. But a vulnerable population may amplify future heat-related health impacts [20]. Shenzhen is a national demonstration zone for sustainable development, with innovation driving China’s long-term growth of megacities. Although it is the country’s fourth most populous city proper, driven as a global center in technology, research, manufacturing, finance, and transportation, the advancements and sustainable megacity reform initiatives in Shenzhen will be invaluable for other areas [21]. But the trend of increased LST, losses of vegetation, and heat-related health issues from the extreme heat of summer temperatures in the city are increasing, and limited quantified studies have been conducted to analyze them.
Specifically, this research was conducted as a comparative study on the urban heat island effect, heat stress and public health issues in Shenzhen, China, and Colombo, Sri Lanka. Both cities are at high risk of increasing urban heat due to rapid urbanization. Sri Lanka is an island surrounded by the Indian Ocean. The hot weather is dangerously growing [22]. Today Sri Lanka also faced a highly stressful situation with heat stress [23,24,25]. The trend of mean annual temperatures has been increasing in all meteorological stations in the country. Out of 15 stations, 13 displayed statistically significant (p < 0.0001) growing trends [26]. Even traditionally cooler areas are experiencing warmer weather conditions. Also, there have been many hospital admissions because of heat stroke, heat rash (skin), and dehydration [27]. Colombo is the commercial capital and the largest city of Sri Lanka by population [28]. And also it is the administrative capital of the Western Province and the district capital of Colombo District [29,30]. Because of its central location in the Indian Ocean, Colombo port has been one of the most famous in South Asia since ancient times. Due to the high rate of urbanisation, vegetation cover is rapidly diminishing and being replaced by buildings, roads, parking lots, pavements, and other structures.
Furthermore, the region’s high vehicle density contributes to the emission of waste heat and pollutant gases. Therefore, previous research has identified Colombo as the most polluted city on the island [31]. As a result of these factors, there is a significant potential for the formation of UHIs within the urban core of Colombo [29].
Both cities of Shenzhen and Colombo are busy with various development projects, business activities and trading. Additionally, the rapid population growth and expansion of urban areas toward the suburban and rural sides address the number of risk cases directly and indirectly combined with urban heat island effects and human health. Although different climatic classifications and research can be seen to identify the various climatic patterns in the cities, attempts to measure heat stress and the urban heat island effect by combining heat-related illnesses are rare and limited to a certain framework. So far, comparing studies addressing all the above issues in spatial and temporal aspects is hard. It seems to be a significant shortcoming in meteorological and urban studies.
This study identifies these issues and their relationships by comparing two coastal cities in tropical and subtropical regions. The results can be used for future urban planning projects to understand the negative impacts of heat stress, UHI effects, and heat-related health issues. At the same time, this research would benefit future studies addressing issues like quantifying the negative heat effect of urban expansion, identifying the population at risk, and identifying the vulnerable sites. And it will also help plan and make policies for various development sectors in both cities.

Conceptual Framework

According to the theoretical framework (Figure 1), a man’s thermal comfort mainly depends on four factors. Namely physical activities (indoor and outdoor), climate factors (temperature, humidity, wind, sunshine hours, and Rainfall), clothing pattern (exposure and cover), and internal factors (sweating, metabolic rate, personal adjustment). The specific environmental factors causing heat stress are less movement of air, high humidity, high temperature and radiant heat. Physical works, such as indoor and outdoor activities, also contribute to the total heat stress of a job. In another way, thermal characteristics and the type of clothing worn are also important. Because some clothes are helpful to heat exchange between the skin and the air, but some are stressful for the remaining atmospheric condition.
The rapid growth of urbanization leads to multiple diseases, mental and physical stress, and so on. As a result of the development or expansion of cities and towns, significant vegetation loss can be seen. In addition, urban surfaces are covered with buildings, resulting in less shade and less moisture to keep urban areas cool. That caused an increase in the pattern of urban heat patches. Compared to rural surroundings, these conceptual issues warm the urban areas and generate heat-related health issues.

2. Materials and Methods

2.1. Description of the Study Area

Colombo metropolitan area in Sri Lanka and Shenzhen city in China have been selected as study areas for the research (Figure 2). Shenzhen is located in the southeast coastal region of China between 113°46′ E and 114°37′ E and between 2°27′ N and 22°52′ N, with a low altitude in a subtropical region. The city is on the east bank of the Pearl River estuary on the central coast of Guangdong’s southern province, with Margin Hong Kong to the south, Dongguan to the north, and Huizhou to the northeast [32]. Its total continental area is 1952 km2. Shenzhen’s climate is highly influenced by the monsoon, with very mild winters and hot, muggy, and rainy summers [33]. In recent decades, the city has evolved from a small border town to a highly urbanized regional economic center. Therefore, it has a solid foundation for exploring the sustainable development of megacities. Therefore, Shenzhen should follow the national sustainable development strategy, formulate distinctive urban sustainable development strategies, and decompose the national sustainable development strategy goals into local goals with specific implementation plans and operability [34,35]. But Shenzhen’s urban expansion has nearly reached saturation due to the rapid transformation of urban ecological land into hard urban construction lands. That creates several problems in the city center.
Colombo is the commercial capital and the most populous city of Sri Lanka. Due to the location of Sri Lanka, within the tropics between 5°55′ N to 9°51′ N and between 79°42′ E to 81°53′ E, the climate of the island could be characterized as tropical [36]. The city has a warm, humid climate without significant seasonal variations in temperature and humidity [37]. The city has become highly urbanized in the South Asian region due to modern urban development and planning projects. Among the cities of South Asia, the Colombo Metropolitan Area (CMA) is the only one that has seen rapid urban growth and fast changes in surface cover in the last two decades [38]. These massive changes in the urban environment also directly influence the urban heat island hazard and heat-related health problems.
Furthermore, the city is mostly unplanned, and one of the major sources of UHIs is buildings with large footprints roofed with low albedo roofing material such as asbestos or concrete [29]. That creates several heat hazards inside the city. Considering all these issues, this research tried to analyse the urban heat island effect, heat stress and public health issues in Colombo, Sri Lanka and Shenzhen, China, comparatively.

2.2. Heat-Related Illnesses Based on the WBGT Index

They analyzed heat-related illnesses based on the Wet Bulb Global Temperature Index (WBGT) from 2008 to 2018. The Australian Bureau of Meteorology (ABM) has published a simple formula to calculate the WBGT index, and that requires as input only water vapour pressure (ρ) and air temperature (Ta). The website (http://www.bom.gov.au, accessed on 1 March 2023) gives a standard physical science formula to calculate ρ from RH and Ta. The study was based on secondary data from two Meteorology Departments in Colombo, Sri Lanka and Shenzhen, China. Thermal sensation and heat-related health impacts in humans have been classified using the index values of WBGT [39]. The equation to calculate the WBGT index is as follows [40].
WBGT   ° C = 0.567   Ta   + 0.393   ρ   + 3.94
ρ   hPa =   RH / 100 × 6.105 exp 17.27 Ta / 237.7 + Ta  
where Ta—air temperature; ρ—water vapour pressure, RH—Relative Humidity.
The classification was applied for both cities of Colombo and Shenzhen to identify the population at risk monthly and seasonally (Table 1). The WBGT was initially formulated by the United State Navy [41,42] and was developed as a screening tool for analyzing environmental heat stress during industrial, military, sports, and occupational activities. The index has been recommended by the International Organization for Standardization (ISO) certification [40,43]. According to that, WBGT can be taken as a user-friendly and reliable measure of heat stress and health impact assessment. SPSS software and Mann Kendall test [44,45,46] have been applied for statistical analysis.

2.3. Method of Land Use Classification

There is an interaction between less vegetation, expansion of built-up surfaces and Land Surface Temperature (LST). These combinations, directly and indirectly, lead to the development of an uncomfortable urban environment and making several health hazards. Landsat-5 TM images and Landsat- 8 OLI/TIRS images obtained in 1997, 2009 and 2019 were used for land-use classification and LST calculation in this study (Table A1). Spatial analysis was entirely based on the satellite images downloaded from https://earthexplorer.usgs.gov (accessed on 1 March 2023), and ArcMap 10.7 software was used for the map preparations.
The Normalized Difference Vegetation Index (NDVI) is developed for estimating vegetation cover from the reflective bands of satellite data [47]. NDVI has been one of the most commonly used vegetation indices in remote sensing since its introduction in the 1970s. Moreover, created NDVI images could be used to identify the pattern of land use changes between two different dates [48]. Using Landsat imageries, NDVI values are calculated using the following formula. But the band number of the images changes with the images of Landsat 05 (NDVI = (Band 4 − Band 3)/(Band 4 + Band 3)) and Landsat 08 (NDVI = (Band 5 − Band 4)/(Band 5 + Band 4)) (NDVI, FAQ, 2019). The equation to calculate the NDVI index is as follows [49].
NDVI   = NIR     RED / NIR   +   RED  
where, NDVI = Normalized Deference Vegetation Index, NIR = Near Infrared Band, RED = Red Bad

2.4. LST Estimation

Landsat data has been downloaded from the USGS ’Earth Explorer’ website. The multispectral bands of the Landsat-5 TM data have a 30 m spatial resolution, while the thermal band (band 6) has a 120 m spatial resolution (https://landsat.usgs.gov, accessed on 1 March 2023). For the Landsat 8 OLI/TIRS data, the multispectral bands also have 30 m spatial resolution, and its thermal bands (band 10 and 11) have 100 m spatial resolution, which has also been resampled to 30 m by the USGS [50]. During the selection of image data, cloud-free images (<10%) have been selected. In addition, radiometric calibration and atmospheric correction have been done for all data sets. The equations to calculate the LST by Landsat Images are as follows [50,51,52].
Landsat-5 TM data
  • Step 1. Thermal Image
Band   6 . TIF
  • Step 2. Conversion DN to Radience
L λ   = LMAX λ LMIN λ QCALMAX QCALMIN × BAND 6 CALMIN +   LMIN λ
  • Step 3. Convert Radiance in to BT (in Kelvin)
T = K 2 ln K 1 L λ + 1
  • Step 4. Convert Degree Kelvin in to Degree Celsius
C = K 273.15
Landsat-8 OLI/TIRS data
  • Step 1. Thermal Image
Band   10 . TIF
  • Step 2. Convert DN to Radience
L λ = M L ×   Q cal + A L
  • Step 3. Brightness temperature
T 10 = K 2 ln K 1 L λ + 1   273.15
  • Step 4. Calculate NDVI
NDVI = Band   5     Band   4 Band   5 +   Band   4
  • Step 5. Proportion of Vegetation (Pv)
Pv = Square   NDVI     NDVImin / NDVImax     NDVImin
  • Step 6. Land Surface Emissivity
ε   = 0.004 ×   Pv   + 0.986
  • Step 7. Land Surface Temperature (LST)
LST   = BT / 1 + 0.00115 ×   BT / 1.4388 ×   Ln ε
The LST from raw Landsat datasets requires the conversion of the DN values of the thermal bands into radiance values, then which are converted to the satellite brightness temperatures (BT) (in Kelvin) [53]. Finally, convert Kelvin (K) into the degree of Celsius (°C) (Equations above have described the way of LST map preparation by using Landsat 5 TM (Thermal Band 06) and Landsat 8 OLI/TIRS data (Thermal Band 10)). Then, NDVI, Proportion of Vegetation (Pv) and Land Surface Emissivity (ε) must be calculated to obtain the final result of LST in Landsat 8 OLI/TIRS data. Here the brightness temperature is the temperature of a blackbody that would emit an identical amount of radiation at a definite wavelength [54]. Thermal Infrared (TIR) pixel values are firstly converted into radiance from Digital Number (DN) values. Then, radiances for the TIR band of Landsat 5 TM, the LMIN, LMAX, QCALMIN, and QCALMAX values are obtained from the metadata file of the Landsat 5 TM data set.
Their T values are referred as a satellite brightness temperature in Kelvin, K1 (Watts/(m2_srad__m)) and K2 (Kelvin) are the calibration constants, and L is the spectral radiance ((Landsat 5 (Band6); K1 = 607.76 and K2 = 1260.56) (Landsat 8 (Band 10); K1 = 774.89 and K2 = 1321.08)) [50,51]. For Landsat 8 data, calculation of NDVI, Proportion of Vegetation (Pv), and Land Surface Emissivity (LSE) (ε) is essential. In this study, we used m = 0.004 and n = 0.986. Pv is the proportion of vegetation extracted from the NDVI Equation. The NDVImin and NDVImax are the minimum and maximum values of the NDVI, respectively. The Calculation has been done step by step to achieve the final result of LST separately for Landsat-5 and Landsat-8 data sets.

3. Results

3.1. The Temperature Fluctuations of Both Cities

When comparing both cities, apparent differences can be observed in the pattern of climatic phenomena (Table 2). Except from June to September, the temperature in Colombo is higher than in Shenzhen throughout the year. Both cities are located in the coastal belt, and the influence of the sea breeze makes a somewhat comfortable environment. Temperature fluctuation also can be observed as a clear pattern due to the monsoon and seasonal behavior of the cities. This research mainly compares the two cities located in two different climatic zones (Figure 3). Therefore, the temperature transitions from subtropical (Shenzhen, China) to tropical (Colombo, Sri Lanka) can be observed.
Especially the fluctuation of minimum and maximum temperatures of both cities is making a significant pattern throughout the year. During the winter, the maximum temperature of Shenzhen city is below the minimum temperature of Colombo. Except for the period from June to October, throughout the year, the temperature of Shenzhen city is lower than the Colombo city. From January to December, Colombo city’s Minimum and Maximum temperatures flow as an equal line. But the curve shape of the pattern in temperature fluctuation can be observed in Shenzhen city due to the seasonal variation (Figure 4).
Only in the summer season, the line of temperature in Shenzhen exceeds rather than the temperature in Colombo. In Shenzhen, the temperature fluctuates from 15 °C to more than 29 °C. As a result, human heat stress and city comfortability also change with the fluctuation of temperature in the city. But significant changes cannot be identified in the pattern of fluctuation of temperature in Colombo. Throughout the year, the temperature of Colombo has varied from 27 °C to 28 °C. Therefore, the whole year, hot and humid climatic conditions can be observed in the atmosphere. It leads to making stressful conditions inside the city. But other factors, such as wind speeds, amount of rainfall in the year and seasonal behavior are affected to minimize the uncomfortable condition in both cities.
According to the annual rainfall pattern, Colombo receives more than 250 mm of rainfall between April to May and September to November. Shenzhen also receives more than 250 mm from June to September. Therefore, the peak temperatures of both cities are directly related to the highest rainfall months of the year.

3.2. Heat-Related Health Impact Based on the WBGT Index

Increasing humidity and temperature are having a negative effect on human body comfort simultaneously. Therefore, it is clear that there is a mutual relationship between temperature, humidity, and human body comfort, as a comparison study here mainly collected temperature and RH% data for both cities of Shenzhen and Colombo to compare and contrast the thermal sensation and heat-related health effects on the public, based on the WBGT index. A rapid increase in temperature and humidity creates an uncomfortable environment and unfavourable effects on the body’s metabolic rate. Therefore, there is a strong interconnection between temperature and RH%. On the other hand, normal air temperature between 25–26 °C and atmospheric RH concentration between 75% and 80% create a comfortable environment (Table 3). Thereby heat stress is minimized, and it leads to a convenient environment for the existence of human beings.
Seasonal influences are directly affected to determine the pattern of comfortable and uncomfortable thermal sensations in the human body. Especially during the summer, the temperature is higher than 29 °C in Shenzhen city. At the same time, atmospheric RH is also higher than 80%. This particular combination is supposed to increase the values of the WBGT index. Therefore, “Very severe stress” (WBGT: ≥32.2) with “Heatstroke highly likely” condition is created by the above factors in the city of Shenzhen during the period from June to August. Most heat-related health impacts can also be seen in this period. Autumn is “Partial discomfort” (WBGT: 27.8–29.4). But during this period, severe heat-related health impacts cannot be seen. The most common condition is “Fatigue”, which everyone feels. But in the winter (December to February), the above process goes the opposite direction. During this season, the temperature is <20 °C and RH is ≤75%. As a result of that, the index value of WBGT is recorded as <27.7. Therefore, the entire season is “comfortable”, and heat-related health effects cannot be observed. The same comfortable condition can also be observed in the spring (March to May) without any health effects related to heat stress (Table 4).
When Colombo and Shenzhen cities are compared, huge differences can be observed in the pattern of monsoons. Throughout the year, Colombo does not come under the “Comfortable” category (WBGT: ≤27.7). Because during the whole year, the temperature does not come down below 27 °C, and RH is also fluctuating between 75% and 84%. As a result, throughout the year, the value of the WBGT index is ≤30.0. Especially in the period of the First Inter Monsoon (FIM) and South West Monsoon (SWM), Colombo city remains under the “Very severe stress” category with “Heatstroke highly likely”. But During the Second Inter Monsoon (SIM) period, the “Severe stress” condition is recorded with a WBGT index score of 31.6, and health effects such as “Heat cramps or heat exhaustion likely and heatstroke are possible”. North East Monsoon (NEM) is better than the other three monsoons in the year. Still, according to the classification, NEM comes under the “Discomfort” category: "Heat cramps and heat exhaustion are possible”. Because of the low temperature and RH values recorded during this period, NEM is fairly more comfortable than other monsoons (Table 4).
Apparent differences in the monthly pattern of WBGT can be observed between Colombo and Shenzhen. “Very severe stress” months in Shenzhen are June, July, and August (3 months). But Colombo city comes under “Very severe stress” with “Heatstroke highly likely” during the five months from April to August. “Comfortable” and “Partial discomfort” conditions cannot be observed in Colombo even in one month of the year. But January, February and December are somewhat better than the other months of the year. Generally, throughout the year, Colombo city is in discomfort and stressful condition. But the “Comfortable” condition (WBGT ≤ 27.7) is maintained continuously in Shenzhen city for the six months from November to April. The reason for that is the temperature of Shenzhen city goes down under 23 °C, and RH fluctuates between 64–81% during the above months.
According to the summary statistic of the WBGT index (Table 5), three months (July, September, and October) recorded a significant positive increasing trend of the WBGT index in Colombo city. November is also very close to the significant level (p = 0.05). Throughout the year, all the mean values are recorded as WBGT ≥ 30. When considering the city of Shenzhen, five months (April, May, June, July, and August) recorded a significant positive increasing trend (p < 0.05) in the WBGT index, and many of them are included in the summer season. On the other hand, January and November are very close to a significant level (p = 0.05) (Figure 5). From May to September, all the mean values are recorded as WBGT ≥ 30. But December to February Mean values is recorded as WBGT < 20.
Both cities of Colombo and Shenzhen have a positive co-relationship between Temperature–WBGT and RH–WBGT. But it is clear that when Colombo and Shenzhen cities are compared, different distribution ranges of values in temperature, RH and WBGT can be observed. When Colombo city is concerned, no huge difference in values between temperature and WBGT index can be observed. Values of the upper level are gathered in one place when Figure 6 is considered. The main reason for that is temperature, and WBGT index are reached up to 27 °C in Colombo city throughout the year. This condition directly affected the upper distribution of co-relationship in WBGT and the Temperature of the city of Colombo. When Shenzhen city is concerned, different ranges of temperature and WBGT values distribution can be observed throughout the year. The main reason for this distribution is temperature, and WBGT values fluctuate between ≤17 and ≤33.8 due to seasonal behavior.
Nevertheless, the same pattern of positive co-relationship can be observed between RH and WBGT in Colombo and Shenzhen. When Colombo city is concerned, RH values are established between 75≤ and ≤84 in the entire year. In the meantime, values of the WBGT index are moved up to higher ranges. But in Shenzhen city, RH values fluctuate between 65≤ and ≤84, and values of WBGT fluctuate between 17.4≤ and ≤33.8. As a result of this pattern, the co-relationship of RH and WBGT has been distributed, as indicated in Figure 6.

3.3. Land Use Pattern and LST in Colombo Metropolitan Area

When Colombo and Shenzhen cities are concerned, the expansion of built-up areas has been influenced by reducing vegetation cover. Pertaining to that pattern, land surface temperature also can be seen as an increasing trend. When the study was carried on in the Colombo Metropolitan area, extreme expansion of Built-up areas and LST could be observed in the city between 1997 and 2019. Due to the urbanization process, special changes can be visible in areas A, B, C, and D (Figure 7). Fort City and Colombo Harbour are situated in “A” locations, and relatively in 1997 and 2009, Colombo fort city was expanded towards the sea 2019. Colombo municipal council has been symbolized by “D”, Biyagama Free Trade Zone is represented by “C” location, and Ratmalana Airport and its surrounding areas are symbolized by “B”. All these locations are under the rapid urbanization process.
In 1997, the Thick Vegetation cover was 18 km2 and was limited to 4 km2 by 2019. In the same way, Agriculture lands/light vegetation cover also decreased drastically. In 1997 it was 163 km2, but in 2009 it reduced to 145 km2 and was limited to 102 km2 by 2019. But Built-up areas were broadly and rapidly expanded in the city center. In 1997 Built-up areas were 39 km2, and in 2009 it, was increased to 54 km2; by 2019, it has been further increased to 82 km2.
With the expansion of Built-up areas, the land surface temperature is also rapidly expanded in Colombo. The areas related to construction works are appeared to be covered with increased LST and decreased LST values can be observed in vegetation lands. The lands belonging to LST: ≥31 °C category was 2 km2 in 1997. By 2009 and 2019, it gradually increased from 10 km2 to 21 km2. This pattern can be observed clearly through the maps (Table 6). From 1997 to 2019, most of the land areas belong to LST: 27–29 °C category. The expansion of the built-up regions on the coastal side can be seen spatially and temporally. Parallel to that, a rapid increase of LST in coastal lines can also be observed.

3.4. Land Use Pattern and LST in Shenzhen

Shenzhen city has a vast and clear difference in comparison to Colombo city. A considerable diffusion in Thick Vegetation cover and Agriculture Lands can be observed.
Built-up lands have been rapidly expanded in the A, B, C, and D locations. For example, Shenzhen Bao’an international airport is represented by “A”, Mawan Port by “B”, Yantian district is symbolized by “C”, and the city of Shenzhen is represented by “D” (Figure 8). In 1997, the Thick Vegetation cover was 287 km2. Gradually it was reduced to 165 km2 in 2009. Further, it has been reduced to 101 km2 in 2019. As a result, a vast increase in Built-up areas can be seen. In 1997, Built-up areas were expanded up to 426 km2. In 2009 it increased to 592 km2, and by 2019 it increased to 710 km2 (Table 7).
The broad urbanization process is indicated by these special locations, and urbanization leads to the expansion of LST in the city. The special thing in Shenzhen city is that LST distribution has been slightly controlled owing to broad vegetation cover and agricultural Lands. The pattern is indicated in Figure 8. The lands belonging to LST: ≥31 °C category was 4 km2 in 1997. By 2019 it had been expanded up to 16 km2. On the other hand, most of the lands belong to LST: 25–27 °C category. A clear increasing trend has been shown by LST: 29–31 °C category. Under this LST category, 5 km2 of land was recorded in 1997. In 2009, it was increased to 11 km2; by 2019, it was expanded to 24 km2 (Table 7).
Especially in the city of Colombo, minimum LST values are increasing year by year (Figure 9). In 1997 it was recorded as 17 °C. But in 2019, the minimum LST value has grown to 24 °C. When comparing the mean temperature of LST, Colombo consistently outperforms Shenzhen. In both cities, maximum LST values fluctuate between 33 °C to 36 °C. And also, minimum, maximum, and mean LST values are gradually increasing in 2019 compared to 1997 and 2009.

4. Discussion

Colombo has recorded more WBGT index values, air temperature, and LST annually. When considering the WBGT, from April to June, the index values reached more than 32.5. Several heat-related problems are arising in this period. This is consistent with earlier research found that March and May, when the solar zenith angle is at its lowest, are the most uncomfortable months in Sri Lanka. During this period country was right under the sun—high relative humidity results from low wind speed during ITCZ passage and pre-monsoonal showers [55]. From 1997 to 2019, the LST in the city center also expanded rapidly.
The main reason is the fast land-use changes covered by the unplanned urbanization process. Especially in the west coastline, including the port city area, is identified as the highest-risk area in the city for heat stress and heat-related illnesses. Previous research also showed the same LST pattern in Colombo city. Urban growth in the Colombo Metropolitan Area follows a diffusion pattern, particularly along the western coastal line. Because urban surface materials have higher radiation temperatures, the LST pattern also shows a similar gradient compared with the landscape metrics in the city [31]. There have been considerable and rapidly accelerating changes to the land cover since 2010, with roughly 30% of all new construction occurring between 1990 and 2015 in wetlands and agricultural areas [56].
In the same way, the maximum and minimum LST in the city also shows the same positive increasing trend. The minimum LST values are increasing faster than the maximum LST values. Previous studies also found that the nighttime LST trends show the highest statistically significant trend. It is increasing faster than the daytime trends [25].
When considering the city of Shenzhen, Summer is the most vulnerable period with extreme heat-related injuries. Many researchers found that Summer is the hottest season in China. The daily maximum temperature considerably impacted the number of deaths during this period. The daily death toll slowly rises as the maximum temperature exceeds 35 °C [57]. The city is vulnerable to several heat-related injuries from June to August, and the WBGT index values also exceed 32.1. But all other seasons are in Comfort and Partial discomfort categories considering the summer. The western part of the city shows the fast changes in the land use cover compared to the other areas. Previous researchers also found that the west and central parts had worse climates, which were less comfortable. An ideal climate region with a high wind speed and low temperature is southeast of Shenzhen [58]. This research also observed the same spatial pattern of the LST in the city.
Heat Stress and Heat Related Public health depend on many factors, such as climatic factors, physical activities, clothing patterns and internal factors of the body. But this research only discussed the climatic factors and urban heat environment related to the heat-related health impacts on the human body. This is a certain limitation of this study. In addition, we used land surface temperature data rather than air temperature data to quantify the urban heat island and their relations with urbanization. However, air temperature might be a better variable to characterize the actual thermal conditions, especially to address human perception and comfort issues with heat. However, the limitation of air temperature is that it is generally monitored at only a few meteorological stations. Thus, it is difficult to fully represent the spatial heterogeneity of the urban thermal environment. In our previous studies conducted in Shenzhen, a significant positive relationship between LST and air temperature from the nearest meteorological station was found [59]. We thus used LST as a proxy of the urban thermal environment because it is challenging to derive air temperature with very high spatial resolution. Still, it is also a point to be considered in the future.

5. Conclusions

Urban heat island effect, heat stress and heat-related public health issues can be observed in Colombo, Sri Lanka and Shenzhen, China. When comparing both cities, apparent differences can be observed in the temperature pattern. Colombo experiences higher air temperature throughout the year than Shenzhen City, except in the summer. Human heat stress, thermal sensation and city comfortability are also changed with air temperature fluctuation and relative humidity.
(1)
Throughout the year, hot and humid climatic conditions can be observed in Colombo. As a result, Colombo does not come under the “Comfortable” category of the WBGT index. Moreover, nearly half of the year (5 months), the city belongs to the “Very severe stress” category with “Heatstroke highly likely”.
(2)
But seasonal influences are directly affected to determine the pattern of comfortable and uncomfortable conditions in Shenzhen city. Summer is the only season that recorded the “Very severe stress” condition. Spring and winter create a “Comfortable” environment in the city, and WBGT index values also come down. Both cities of Colombo and Shenzhen have a positive co-relationship between the Temperature–WBGT and RH–WBGT index.
(3)
With the rapid growth in Built-up lands, air temperature and LST are also expanded relatively in the city centers. Expansion of LST is a more rapid process in Colombo city than in Shenzhen. Colombo port city and its surrounding areas in the North-Western part along the coast are identified as the most vulnerable locations for LST expansion. In Shenzhen, the urban land along the western coastal line and main roads are expanding faster while the LST patches are increasing.
By comparison of all the factors, this research concludes that Colombo city is more stressful than Shenzhen city. It does not seem like a big issue today. But in the near future, the increasing trend of heat stress, urban heat island effect and heat-related health issues will cause several complications for all nations worldwide. The use of cool and green roofs and the proper distribution of vegetation cover can be applied to the highly vulnerable areas identified in this research. Therefore, it is necessary to assess these risks and respond effectively.

Author Contributions

Conceptualization, S.N. and W.Z.; methodology, S.N. and W.Z.; formal analysis, S.N., W.Z. and J.C.; Resources, W.W.; writing—original draft preparation, S.N. and W.Z.; writing—review and editing, S.N., W.W., J.C., J.W. and W.Z.; supervision, W.Z.; funding acquisition, W.Z. and J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 32001160) and the National Key Research and Development Program (Grant No. 2022YFF1301104).

Data Availability Statement

The data presented in this study are available on request from the author.

Acknowledgments

The help provided by Shenzhen Ecology and Environment National Observation and Research Station for this study was gratefully acknowledged. Thanks to anonymous experts for their suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Details of the Landsat data acquired.
Table A1. Details of the Landsat data acquired.
Date AcquiredSatelliteWRS PathWRS RowDay/Night IndicatorSensor Identifier
Shenzhen29 October 1997LANDSAT_5122044DAYTM
17 October 2009LANDSAT_5122044DAYTM
14 November 2019LANDSAT_8122044DAYOLI_TIRS
Colombo7 February 1997LANDSAT_5141055DAYTM
8 February 2009LANDSAT_5141055DAYTM

References

  1. Semple, E.C. Influences of Geographic Environment: On the Basis of Ratzel’s System of Anthropo-Geography. The Project Gutenberg. 2005. Available online: https://www.gutenberg.org/ebooks/15293 (accessed on 1 March 2023).
  2. Howard, J.C. General Climatology; Prentice Hall of India Ltd.: New Delhi, India, 1995. [Google Scholar]
  3. Becker, J.A.; Stewart, L.K. Heat-Related Illness. Natl. Libr. Med. Cent. Biotechnol. Inf. 2011, 83, 1325–1330. Available online: https://pubmed.ncbi.nlm.nih.gov/21661715/ (accessed on 1 March 2023).
  4. Hughes, L.; Hanna, E.; Jacqui, F. The Silent Killer: Climate Change and the Health Impacts of Extreme Heat. Climate Council of Australia Limited. 2016. Available online: https://www.climatecouncil.org.au/resources/silentkillerreport (accessed on 1 March 2023).
  5. Umair, S. Global Warming: Causes, Effects and Solutions. Durreesamin J. 2015, 1, 1–7. [Google Scholar]
  6. Wang, Y.; Wang, A.; Zhai, J.; Tao, H.; Jiang, T.; Su, B.; Yang, J.; Wang, G.; Liu, Q.; Gao, C.; et al. Tens of thousands additional deaths annually in cities of China between 1.5 °C and 2.0 °C warming. Nat. Commun. 2019, 10, 3376. [Google Scholar] [CrossRef] [PubMed]
  7. Na, W.; Jae-Yeon, J.; Kyung, E.L.; Hyunyoung, K.; Byungyool, J.; Kwon, J.W.; Soo-Nam, J. The Effects of Temperature on Heat-related Illness According to the Characteristics of Patients During the Summer of 2012 in the Republic of Korea. J. Prev. Med. Public Health 2013, 46, 19. Available online: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pmc/articles/PMC3567322 (accessed on 1 March 2023). [CrossRef]
  8. CDC. Heat Related Deaths—United States—1999–2003. Morbidity and Mortality Weekly Report. 2016. Available online: https://www.cdc.gov/mmwr/preview/mmwrhtml/mm5529a2.htm (accessed on 1 March 2023).
  9. Zhou, W.; Yu, W.; Qian, Y.; Han, L.; Pickett, S.T.A.; Wang, J.; Li, W.; Ouyang, Z. Beyond city expansion: Multi-scale enviromental impacts of urban megaregion formation in China. Natl. Sci. Rev. 2022, 9, nwab107. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, J.; Zhou, W.; Wang, J. Time-Series Analysis Reveals Intensified Urban Heat Island Effects but without Significant Urban Warming. Remote Sens. 2019, 11, 1129. [Google Scholar] [CrossRef]
  11. Milica, L.; Jelena, M. UTCI Based Assessment of Urban Outdoor Thermal Comfort in Belgrade: Serbia. International Scientific Conference on Information Technology and Data Related Research: Sinteza. 2020. Available online: https://www.researchgate.net/publication/344942716_UTCI_Based_Assessment_of_Urban_Outdoor_Thermal_Comfort_in_Belgrade_Serbia (accessed on 1 March 2023).
  12. Elspeth, O.; Matt, B.; Lisa, L.; James, A.S.; Alan, C.; Kerstin, Z. Heat, health, and humidity in Australia’s monsoon tropics: A critical review of the problematization of ‘heat’ in a changing climate. WIREs Clim. Chang. Wiley Period. 2017, 8, e468. [Google Scholar]
  13. David, E.P. Urban Heat Island Effects on Estimates of Observed Climate Change; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2010; Volume 1, pp. 123–133. [Google Scholar] [CrossRef]
  14. Li, X.; Zhou, Y.; Yu, S.; Jia, G.; Li, H.; Li, W. Urban heat island impacts on building energy consumption: A review of approaches and findings. Energy 2019, 174, 407–419. [Google Scholar] [CrossRef]
  15. Zhou, W.; Huang, G.; Pickett, S.T.A.; Wang, J.; Cadenasso, M.L.; McPhearson, T.; Grove, J.M.; Wang, J. Urban tree canopy has greater cooling effects in socially vulnerable communities in the US. One Earth 2021, 4, 1764–1775. [Google Scholar] [CrossRef]
  16. Wang, J.; Zhou, W.; Jiao, M. Location matters: Planting urban trees in the right places improves cooling. Front. Ecol. Environ. 2022, 20, 147–151. [Google Scholar] [CrossRef]
  17. Li, Y.; Ren, T.; Kinney, P.L.; Joyner, A.; Zhang, W. Projecting future climate change impacts on heat-related mortality in large urban areas in China. Environ. Res. 2018, 163, 171–185. [Google Scholar] [CrossRef] [PubMed]
  18. Sun, X.; Sun, Q.; Zhou, X.; Li, X.; Yang, M.; Yu, A.; Geng, F. Heat wave impact on mortality in Pudong New Area, China in 2013. Sci. Total. Environ. 2014, 493, 789–794. [Google Scholar] [CrossRef] [PubMed]
  19. Bai, L.; Ding, G.; Gu, S.; Bi, P.; Su, B.; Qin, D.; Liu, Q. The effects of summer tem-perature and heat waves on heat-related illness in a coastal city of China: 2011–2013. Environ. Res. 2014, 132, 212–219. [Google Scholar] [CrossRef]
  20. He, Y.L.; Deng, S.Z.; Ho, H.C.; Wang, H.B.; Chen, Y.; Hajat, S.; Huang, C.R. The half-degree matters for heat-related health impacts under the 1.5 °C and 2 °C warming scenarios: Evidence from ambulance data in Shenzhen, China. Adv. Clim. Change Res. 2021, 12, 628–637. [Google Scholar] [CrossRef]
  21. Anna Bulanda, J. Evolution of the Metropolitan Area of Shenzhen, Analysis: From Theory to Selected Example. In IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2019; Volume 471, p. 112058. [Google Scholar]
  22. Disaster Services. Dangerous Heat Stress Away from the Hills of Sri Lanka Since April. WINNER, Environ. Green Energy. 2019. Available online: https://disaster.lk/category/heat-stress (accessed on 1 March 2023).
  23. Rekha, K.W.G.N. Recent trends of climate change in Sri Lanka. J. Geo Environ. Bangladesh 2009, 3, 37–49. [Google Scholar]
  24. Nanayakkara, N.W.S.; Nianthi, K.W.G.R. Analysis of Human Heat Stress in Sri Lanka: Using Temperature Humidity Index (THI). Int. J. Sci. Res. Publ. 2018, 8, 624–630. [Google Scholar] [CrossRef]
  25. Emmanuel, R. Thermal comfort implications of urbanization in a warm-humid city: The Colombo Metropolitan Region (CMR), Sri Lanka. Build. Environ. 2005, 40, 1591–1601. [Google Scholar] [CrossRef]
  26. Meegahakotuwa, U.S.; Rekha Nianthi, K.W.G. Spatial and Temporal Variation of Temperature Trends in Last Century of Sri Lanka. In Proceedings of the 2nd International Research Symposium—IRSUWC, Uva Wellassa University, Badulla, Sri Lanka, 1–2 February 2018. [Google Scholar]
  27. Sathisraja, A. Don’t Get Burned by the Fire in the Sky. Sunday Times. 2016. Available online: https://www.sundaytimes.lk/160508/news/dont-get-burned-by-the-fire-in-the-sky-193262.html (accessed on 1 March 2023).
  28. Kumarage, A.S. Impacts of Transportation Infrastructure and Services on Urban Poverty and Land Development in Colombo, Sri Lanka. Glob. Urban Dev. 2007, 3, 1–15. [Google Scholar]
  29. Manjula, R.; Takehiro, M.; Matamyo, S.; Yuji, M. Spatial Analysis of Urbanization Patterns in Four Rapidly Growing South Asian Cities Using Sentinel-2 Data. Remote Sens. 2021, 13, 1531. [Google Scholar] [CrossRef]
  30. Liyanage, J. Air Quality Management, City Report-Colombo, Sri Lanka. In Proceedings of the Kitakyushu Initiative Seminar on Urban Air Quality Management, Bangkok, Thailand, 20–21 February 2003; Available online: https://www.iges.or.jp/en/pub/kitakyushu-initiative-seminar-urban-air/en (accessed on 1 March 2023).
  31. Senanayake, I.; Welivitiya, W.; Nadeeka, P. Remote sensing based analysis of urban heat islands with vegetation cover in Colombo city, Sri Lanka using Landsat-7 ETM+ data. Urban Clim. 2013, 5, 19–35. [Google Scholar] [CrossRef]
  32. Shen, J. Urban Growth and Sustainable Development in Shenzhen City 1980–2006. Open Environ. J. 2008, 2, 71–79. [Google Scholar] [CrossRef]
  33. Wang, W.; Liu, K.; Tang, R.; Wang, S. Remote sensing image-based analysis of the urban heat island effect in Shenzhen, China. Phys. Chem. Earth Parts A/B/C 2019, 110, 168–175. [Google Scholar] [CrossRef]
  34. Wang, Y.; Lu, Y.; He, G.; Wang, C.; Yuan, J.; Cao, X. Spatial Variability of sustainable development goals in China: A pro-vinciallevel evaluation. Environ. Dev. 2020, 35, 100483. [Google Scholar] [CrossRef]
  35. Jonathan, Z. Shenzhen Ranked World’s 2nd Best Tourist Destination for 2019. That’s. 2018. Available online: https://www.thatsmags.com/shenzhen/post/25456/shenzhen-ranked-among-top-10-cities-to-visit-in-2019 (accessed on 1 March 2023).
  36. Perera, N.G.R.; Emmanue, R.A. “Local Climate Zone” based approach to urban planning in Colombo, Sri Lanka. Urban Clim. 2018, 23, 188–203. [Google Scholar] [CrossRef]
  37. Senanayake, I.P.; Welivitiya, W.D.D.P.; Nadeeka, P.M. Urban green spaces analysis for development planning in Colombo, Sri Lanka, utilizing THEOS satellite imagery—A remote sensing and GIS approach. Urban For. Urban Green. 2013, 12, 307–314. [Google Scholar] [CrossRef]
  38. Subasinghe, S.; Nianthi, R.; Rajapaksha, G.; Gamage, I. Monitoring the Impacts of Urbanisation on Environmental Sustainability Using Geospatial Techniques: A Case Study in Colombo District, Sri Lanka. J. Geospat. Surv. 2021, 1, 1–13. [Google Scholar] [CrossRef]
  39. Md Motakabber, A.; Gourab, B.; Arkajit; Rina, B. A Journey through heat stresses and its impact on population. Techno. Int. J. Health Eng. Manag. Sci. 2018, 2, 25–30. Available online: https://www.academia.edu/37211400/A_Journey_through_heat_stresses_and_its_impact_on_population (accessed on 1 March 2023).
  40. Bruno, L.; Tord, K. Calculating Workplace WBGT from Meteorological Data: A Tool for Climate Change Assessment. Ind. Health 2012, 50, 267–278. [Google Scholar] [CrossRef]
  41. Yaglou, C.P.; Minard, D. Control of heat casualties at military training centers. Arch. Ind. Health 1957, 16, 302–316. Available online: https://pubmed.ncbi.nlm.nih.gov/13457450 (accessed on 1 March 2023).
  42. Maia, P.A.; Ruas, C.; Bitencourt, D.P. Wet-bulb globe temperature index estimation using meteorological data from São Paulo State, Brazil. Int. J. Biometeorol. 2015, 59, 1395–1403. [Google Scholar] [CrossRef] [PubMed]
  43. Franck, B.; Gregoire, P.M. Is the Wet-Bulb Globe Temperature (WBGT) Index Relevant for Exercise in the Heat? Sport. Med. 2015, 45, 1619–1621. [Google Scholar] [CrossRef]
  44. Arun, M.; Sananda, K.; Anirban, M. Rainfall Trend Analysis by Mann-Kendall Test: A Case Study of north-Eastern Part of Cuttack District, Orissa. Int. J. Geol. Earth Environ. Sci. 2012, 2, 70–78. [Google Scholar]
  45. Kendall, M.G. Rank Correlation Methods, 4th ed.; Griffin: London, UK, 1975. [Google Scholar]
  46. Mann, H.B. Non-parametric test against trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  47. Xu, D. Compare NDVI Extracted from Landsat 8 Imagery with that from Landsat 7 Imagery. Am. J. Remote Sens. 2014, 2, 10–14. [Google Scholar] [CrossRef]
  48. Ehsan, S.; Dashtekian, K. Analysis of land use-land covers changes using normalized difference vegetation index (NDVI) differencing and classification methods. African J. Agric. Res. 2013, 8, 4614–4622. [Google Scholar] [CrossRef]
  49. Ozyavuz, M.; Bilgili, B.C.; Salici, A. Determination of Vegetation Changes with NDVI Method. J. Environ. Prot. Ecol. 2015, 16, 264–273. [Google Scholar]
  50. Manjula, R.; Ronald, C.E. An Urban Heat Island Study of the Colombo Metropolitan Area, Sri Lanka, Based on Landsat Data (1997–2017). Int. J. Geo-Inf. ISPRS 2017, 6, 189. [Google Scholar]
  51. Yasir, M.; Hui, S.; Rahman, S.U.; Ilyas, M.; Zafar, A.; Mehmood, A. Estimation of Land Surface Temperature using LANDSAT-8 Data—A Case Study of District Malakand, Khyber Pakhtunkhwa, Pakistan. J. Lib. Arts Humanit. 2020, 1, 140–148. [Google Scholar]
  52. Dozier, J.; Warren, S.G. Effect of viewing angle on the infrared brightness temperature of snow. Water Resour. Res. 1982, 18, 1424–1434. [Google Scholar] [CrossRef]
  53. Estoque, R.C.; Murayama, Y.; Myint, S.W. Effects of landscape composition and pattern on land surface temperature: An urban heat island study in the megacities of Southeast Asia. Sci. Total Environ. 2017, 577, 349–359. [Google Scholar] [CrossRef] [PubMed]
  54. Sekertekin, A.; Bonafoni, S. Land Surface Temperature Retrieval from Landsat 5, 7, and 8 over Rural Areas: Assessment of Different Retrieval Algorithms and Emissivity Models and Toolbox Implementation. Remote Sens. 2020, 12, 294. [Google Scholar] [CrossRef]
  55. Simath, S.; Emmanuel, R. Urban thermal comfort trends in Sri Lanka: The increasing overheating problem and its potential mitigation. Int. J. Biometeorol. 2022, 66, 1865–1876. [Google Scholar] [CrossRef] [PubMed]
  56. Saparamadu, S.; Yi, Z.; Zongping, Z. Temporal Changes of Land Use Land Cover and Environmental Impacts: A Case Study in Colombo, Sri Lanka. Int. J. Earth Environ. Sci. 2018, 3, 1–13. [Google Scholar] [CrossRef] [PubMed]
  57. Liu, X.; Tian, Z.; Sun, L.; Liu, J.; Wu, W.; Xu, H.; Sun, L.; Wang, C. Mitigating heat-related mortality risk in Shanghai, China: System dynamics modeling simulations. Environ. Geochem. Health 2020, 42, 3171–3184. [Google Scholar] [CrossRef] [PubMed]
  58. Zhang, S.; Fang, X.; Cheng, C.; Chen, L.; Zhang, L.; Yu, Y.; Li, L.; Luo, H. Research on the Planning Method and Strategy of Urban Wind and Heat Environment Optimization—Taking Shenzhen, a Sub-Tropical Megacity in Southern China, as an Example. Atmosphere 2022, 13, 1395. [Google Scholar] [CrossRef]
  59. Cao, J.; Zhou, W.; Zheng, Z.; Ren, T.; Wang, W. Within-city spatial and temporal heterogeneity of air temperature and its relationship with land surface temperature. Landsc. Urban Plan. 2020, 206, 103979. [Google Scholar] [CrossRef]
Figure 1. Conceptual framework.
Figure 1. Conceptual framework.
Atmosphere 14 00839 g001
Figure 2. Location of the study area.
Figure 2. Location of the study area.
Atmosphere 14 00839 g002
Figure 3. Pattern of climatic variabilities in Colombo and Shenzhen (2008–2018).
Figure 3. Pattern of climatic variabilities in Colombo and Shenzhen (2008–2018).
Atmosphere 14 00839 g003
Figure 4. Minimum & Maximum temperature fluctuations.
Figure 4. Minimum & Maximum temperature fluctuations.
Atmosphere 14 00839 g004
Figure 5. The monthly pattern of p values distribution from 2008 to 2018.
Figure 5. The monthly pattern of p values distribution from 2008 to 2018.
Atmosphere 14 00839 g005
Figure 6. Co-relationship between Temperature−WBGT and RH%−WBGT.
Figure 6. Co-relationship between Temperature−WBGT and RH%−WBGT.
Atmosphere 14 00839 g006
Figure 7. LST and NDVI changes in Colombo.
Figure 7. LST and NDVI changes in Colombo.
Atmosphere 14 00839 g007
Figure 8. LST and NDVI changes in Shenzhen.
Figure 8. LST and NDVI changes in Shenzhen.
Atmosphere 14 00839 g008
Figure 9. Histograms of LST diffusion in Colombo and Shenzhen.
Figure 9. Histograms of LST diffusion in Colombo and Shenzhen.
Atmosphere 14 00839 g009
Table 1. Heat exposure and associated health effects with WBGT.
Table 1. Heat exposure and associated health effects with WBGT.
WBGTThermal SensationHealth Effects
≤27.7Comfortable---
27.8–29.4Partial discomfortFatigue
29.5–31.0DiscomfortHeat cramps and heat exhaustion are possible
31.1–32.1Severe stressHeat cramps or heat exhaustion likely, and heatstroke is possible
≥32.2Very severe stressHeatstroke highly likely
Source: [39].
Table 2. Summary statistics of temperature and RH (%) fluctuations (2008–2018).
Table 2. Summary statistics of temperature and RH (%) fluctuations (2008–2018).
MonthsCityMinimumMaximumMeanStd. Deviationp-ValueSen’s Slope
TRHTRHTRHTRHTRHTRH
JanSH12.056.818.284.015.271.11.69.10.270.080.251.72
Co26.672.528.379.527.275.30.42.00.390.520.05−0.18
FebSH12.665.421.181.716.675.12.34.10.750.43−0.11−0.14
Co25.773.028.780.027.476.80.82.10.430.260.05−0.25
MarSH17.665.520.987.819.578.31.06.90.750.01−0.031.33
Co27.477.529.582.528.379.20.51.50.230.470.08−0.16
AprSH21.671.524.394.823.281.90.76.50.350.030.051.24
Co27.577.529.984.028.681.50.72.10.130.040.14−0.41
MaySH25.873.427.993.726.683.00.65.90.060.040.081.33
Co28.377.029.285.028.981.80.32.31.000.300.000.30
JunSH26.978.929.892.928.685.20.84.50.040.080.160.98
Co27.881.029.084.028.482.70.41.10.290.810.060.00
JulSH28.777.330.289.429.383.20.44.20.430.010.021.04
Co27.479.529.083.528.381.40.51.30.020.050.13−0.31
AugSH28.875.330.290.229.581.50.44.80.870.00−0.001.25
Co27.578.529.183.528.281.30.41.40.580.330.05−0.16
SepSH28.270.629.682.728.977.40.54.80.750.020.011.18
Co27.779.028.686.028.081.70.22.10.130.510.030.07
OctSH25.261.027.081.326.269.70.66.40.750.270.040.85
Co27.280.028.486.027.783.20.32.00.690.110.010.28
NovSH19.558.523.978.821.971.11.16.70.270.030.061.52
Co27.081.527.685.527.383.80.21.50.750.510.000.00
DecSH15.052.119.179.117.063.31.18.40.871.000.010.10
Co26.576.527.986.027.380.10.42.90.180.930.050.00
Note: SH: Shenzhen, Co: Colombo, T: Temperature (°C), RH: Relative Humidity (%). p < 0.05 shows a significant trend in the distribution of data.
Table 3. Relationship of Temperature, RH% and WBGT index (2008–2018).
Table 3. Relationship of Temperature, RH% and WBGT index (2008–2018).
CityJanFebMarAprMayJunJulAugSepOctNovDecData
Shenzhen707578818385838177697164RH
15.316.719.623.226.628.729.329.528.926.322.017.0MT
17.419.022.026.230.333.233.833.832.128.123.718.4WBGT
Colombo30.030.532.032.833.132.732.332.232.031.931.430.8WBGT
27.227.528.428.728.928.528.328.228.127.827.427.3MT
757779828283818182838480RH
Note: RH: Relative Humidity, MT: Mean Temperature, WBGT: Wet Bulb Global Temperature. Atmosphere 14 00839 i001 ≤27.7, Atmosphere 14 00839 i002 27.8–29.4, Atmosphere 14 00839 i003 29.5–31.0, Atmosphere 14 00839 i004 31.1–32.1, Atmosphere 14 00839 i005 ≥32.2.
Table 4. Monthly and seasonal variation of WBGT index. The classification of WBGT, Thermal Sensation and Heat Effects were referred from the previous research [39].
Table 4. Monthly and seasonal variation of WBGT index. The classification of WBGT, Thermal Sensation and Heat Effects were referred from the previous research [39].
WBGTThermal SensationsHealth Effects ColomboShenzhen
MonthMonsoonMonthSeason
≤27.7Comfortable------January, February, March, April, November, DecemberWinter, Spring
27.8–29.4Partial discomfortFatigue----OctoberAutumn
29.5–31.0DiscomfortHeat cramps and heat exhaustion are possibleJanuary, February, DecemberNEMMay--
31.1–32.1Severe StressHeat cramps or heat exhaustion likely, and heatstroke is possibleMarch, September, October, NovemberSIMSeptember--
≥32.2Very severe stressHeatstroke highly likelyApril, May, June, July, AugustFIM, SWMJune, July, AugustSummer
Note: FIM: First Inter Monsoon, SWM: Southwest Monsoon, SIM: Second Inter Monsoon, NEM: North East Monsoon.
Table 5. Summary statistics of WBGT in Shenzhen and Colombo (2008–2018).
Table 5. Summary statistics of WBGT in Shenzhen and Colombo (2008–2018).
VariableMinimumMaximumMeanStd. Deviationp-Value
CoShCoShCoShCoShCoSh
Jan29.414.031.420.630.017.40.61.70.870.06
Fe28.614.832.123.430.519.01.02.30.750.87
Mar31.319.133.323.932.022.00.51.20.100.53
Apr31.624.233.929.032.826.20.71.40.110.04
May32.528.233.532.833.130.40.31.40.530.03
Jun31.831.033.335.832.733.30.51.50.310.00
Jul31.332.732.935.732.333.80.50.90.020.01
Aug31.532.333.034.832.233.80.50.80.430.00
Sep31.430.832.434.032.032.40.30.90.040.21
Oct31.226.132.330.431.928.10.31.20.020.21
Nov31.020.932.026.431.423.70.31.50.050.06
Dec30.116.231.920.530.818.40.51.20.350.76
Note: A common threshold of the p-value is 0.05. It is a probability score used in the statistical test to identify the statistical significance of an observed effect; Co: Colombo, Sh: Shenzhen.
Table 6. Land use pattern changes and LST expansion in Colombo (Area Km2).
Table 6. Land use pattern changes and LST expansion in Colombo (Area Km2).
YearsNDVI
Water BodiesBuild-up AreasHome Gardens/Sub-urbanAgricultural Lands/Light VegetationThick Vegetation Cover
1997104398216318
200910354971457
201998821201024
LST (°C)
≤2525–2727–2929–3131≥
1997109137128302
2009102681458110
2019924114910321
Table 7. Land use pattern changes and LST expansion in Shenzhen (area km2).
Table 7. Land use pattern changes and LST expansion in Shenzhen (area km2).
YearsNDVI
Water BodiesBuild-up AreaHome Gardens/Sub-urbanAgricultural Lands/Light VegetationThick Vegetation Cover
199735426371395287
200925592440344165
201922710408325101
LST (°C)
≤2121–2323–2525–2727–2929–3131≥
19972255693843205254
200999304521503120118
2019782984105681722416
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nanayakkara, S.; Wang, W.; Cao, J.; Wang, J.; Zhou, W. Analysis of Urban Heat Island Effect, Heat Stress and Public Health in Colombo, Sri Lanka and Shenzhen, China. Atmosphere 2023, 14, 839. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos14050839

AMA Style

Nanayakkara S, Wang W, Cao J, Wang J, Zhou W. Analysis of Urban Heat Island Effect, Heat Stress and Public Health in Colombo, Sri Lanka and Shenzhen, China. Atmosphere. 2023; 14(5):839. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos14050839

Chicago/Turabian Style

Nanayakkara, Srimalee, Weimin Wang, Jie Cao, Jia Wang, and Weiqi Zhou. 2023. "Analysis of Urban Heat Island Effect, Heat Stress and Public Health in Colombo, Sri Lanka and Shenzhen, China" Atmosphere 14, no. 5: 839. https://0-doi-org.brum.beds.ac.uk/10.3390/atmos14050839

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop