Contributions of sea–land breeze and local climate zones to daytime and nighttime heat island intensity

The acceleration of global urbanization has increased the frequency of the urban heat island (UHI) effect and heatwaves, which seriously endanger human health. We used Shenzhen as a case study to examine the daytime and nighttime differences in UHI intensity (UHII), considering different local climate zones (LCZs) and sea–land breezes. The diurnal UHII was >3 °C for 52% of the study period, whereas the nocturnal UHII was >3 °C for only 26% of the study period. The average diurnal and nocturnal building-type UHII values were 2.77 and 1.11 °C higher than those of the natural type, respectively. Sea breezes alleviated the UHI effect with a linear correlation coefficient of −0.68601 between them. Moreover, diurnal and nocturnal UHII showed differences across different gradients, which can help guide urban planning.


INTRODUCTION
Urbanization has been accelerating since the beginning of the 21st century [1][2][3] . In China, as the population increases annually 4 , the areas under impervious layers have also witnessed increases, causing surface temperatures to increase and exacerbate the urban heat island (UHI) effect, extreme weather, and heat weaves [5][6][7][8] , all of which affect human comfort and threaten human health [9][10][11][12][13] , and even cause death. The UHI phenomenon was first discovered by Howard 14 , who conducted temperature observations in London, UK, and urban thermal environments have been the subject of active research ever since.
Land surface temperature is an important parameter for urban thermal environment research 15 ; therefore, many international studies have investigated methods to determine it. Compared with traditional meteorological station observations, thermal infrared (TIR) image data can retrieve land surface temperature more effectively. Owing to the advantages of this method, such as its spatial continuity and easily accessible images, it has been widely used in research on UHIs [16][17][18][19][20][21][22] . Depending on the source of the image sensor, TIR images can include high, medium, and low spatial resolution data, such as Landsat/TIRS (100 m), ASTER (90 m), GF5/VIMS (40 m), HJ-1B/IRS (300 m), and MODIS (1000 m). Presently, studies on urban thermal environments commonly use Landsat TM/TIRS or MODIS data for land surface temperature inversion and UHI intensity (UHII) calculations, resulting in the development of heat island mitigation strategies 18,20,[23][24][25][26][27][28][29][30][31] . However, the former is difficult to acquire at night, whereas the latter has a low resolution and cannot accurately analyze thermal environments at the urban-block scale. Therefore, MODIS data are more suitable for large-scale research on regional thermal environments. ASTER data have a high spatial resolution and can provide suitable night images. These data have become the main source for daytime and nighttime UHI research, respectively; therefore, AST_08 data were used in this study to conduct land surface temperature inversion calculations.
UHIs are affected by many factors, including architecture, climate, and land map [32][33][34][35][36][37][38] . However, owing to rapid urbanization in China, available morphological and land-use data are inevitably unstandardized, difficult to obtain, or cannot be used for urban planning or environmental research. The local climate zone (LCZ) classification scheme constitutes a climate-related land cover classification system for urban structures, which was first proposed by Stewart et al. 39,40 . This classification scheme has been used worldwide as an international standard for climate-related research [41][42][43] and is widely used in UHI research [44][45][46][47][48][49] . Presently, LCZ classification methods can primarily be divided into two types. The first uses the World City Database Portal Tool (WUDAPT) [50][51][52][53] : the software required for this method is free, and the data are easily obtainable, but its spatial resolution is low; thus, it is only suitable for large-scale research. The second method uses a geographic information system (GIS) to calculate building parameters and vegetation indicators to divide LCZs into smaller study units [54][55][56][57] . Although obtaining data for this method is more difficult, the classification results correspond well with the actual field conditions. Therefore, the second classification method was used herein and building and land-use data were combined to classify the LCZs.
Recent studies have shown that urban ventilation can improve urban thermal comfort and alleviate the UHI effect [58][59][60][61] . Sea and land breezes, representative of typical local winds in coastal cities and which alternate during days and nights, refer to the mesoscale circulations formed via the temperature difference between the ocean and land. As they frequently occur in coastal cities, previous research has attempted to address the interactions between sea and land breezes and UHIs [62][63][64][65] . For instance, Sangobanwo 66 revealed that coastal cities are more susceptible to UHIs due to sea-land breezes and are generally more prone to extremely high temperatures than inland cities. Wang et al. 67 used the Weather Research and Forecasting (WRF) model to simulate the sea-land breeze and thermal environment, further analyzing its effect on the UHI. Shen and Yuan 68 used the large eddy simulation module in the WRF model to simulate five urban cases and one non-urban case. They analyzed the interactions between sea-land breezes and UHIs, and found that sea-land breeze circulation was stronger in urban areas; moreover, they revealed that UHIs were alleviated when the wind speed was high and sea breezes flooded city interior. Sea-land breezes can also significantly affect UHIs. However, most studies have only qualitatively analyzed the weakening effect of sea breeze on the strength of the heat island. Few studies analyze the correlation between sea and land wind speed and UHII. In addition, discussing the strength of the sea and land breeze and heat island under different gradients (distance from the coastline) is necessary. This has important implications for urban planning and a rational layout of the city. Therefore, in this study, the effects of sea-land breezes on the daytime and nighttime UHIIs were analyzed under different gradients based on meteorological data.
The main purpose of this study was to explore the internal differences in the urban thermal environment of Shenzhen from the perspective of LCZs and sea-land breezes. The main objectives were to: (1) analyze the UHII day and night spatial  differences based on the land surface temperature directly calculated from AST_08 images, (2) analyze the influence of different LCZ types on the intensity of heat islands based on Oke's classification standard for LCZs, and (3) analyze the daily variation in the UHIIs between typical and atypical sea and land breeze days based on weather station data, and investigate the influence of the sea and land wind speeds on the intensity of heat islands under different gradients using correlation analysis. This study can provide a reference for urban planning and for mitigating the UHI effect.

RESULTS
Daytime and nighttime UHII and LCZ spatial distributions Overall, the daytime and nighttime UHII spatial distributions ( Fig. 1) were highly consistent with the spatial distribution of the LCZ type (Fig. 2). Among these (Table 1), a large amount of building-type data were collected from the eastern Longhua District, and in the southeastern Bao'an, Nanshan, and Futian districts. The order of the LCZs of residential areas was LCZ5 > LCZA > LCZD > LCZF > LCZ4 > LCZ10 > LCZ3 > LCZE > LCZB > LCZ1 > LCZ2 > LCZ9 > LCZC > LCZ8 > LCZ6 > LCZ7. Moreover, the proportion of natural types in Bao'an, Longhua, Futian, and Nanshan Districts was higher than that of building types (Bao'an District: 0.240 > 0.233, Longhua District: 0.113 > 0.106, Futian District: 0.064 > 0.027, Nanshan District: 0.156 > 0.062). Among these, Bao'an and Longhua districts collected large-scale data on LCZ5, LCZA, LCZD, LCZE, and LCZF. In particular, the corresponding daytime and nighttime UHII values were relatively high for compact super high-rise (LCZ1), compact high-rise (LCZ2), compact mid-high-rise (LCZ3), and compact mid-rise (LCZ4) buildings. Meanwhile, a large amount of natural type data (LCZA-LCZF) and compact low-rise buildings (LCZ5) were present in western Bao'an District and western Longhua District, but the corresponding daytime and nighttime UHII values were low. Therefore, these building types (LCZ1-LCZ10) corresponded to higher UHII values, whereas the natural types (LCZA-LCZF) corresponded to lower UHII values. In addition, the statistics on the proportions of the daytime and nighttime UHII values revealed that the nighttime UHII values were lower than those during the daytime (Fig. 1). During the day, the UHII was >3°C for 52% of the study period, which was only true for 26% of the study period at night. For the negative heat island, the maximum intensity (<−3°C) also showed that day (29%) > night (19%); however, the UHII values in each interval, ranging from −1 to 3°C, show that the proportion at night is greater than that during the day (28% (night) > 7% (day),16% (night) > 6% (day), 11% (night) > 6% (day), and 13% (night) > 12%(day)). Table 2 lists the average daytime and nighttime UHII values corresponding to other LCZ types, except water bodies (LCZF). As AST_08 recorded the dynamic temperatures of the Earth's land surface, water bodies were not analyzed. LCZ5 and LCZA exhibited the lowest nighttime and daytime UHIIs (−0.84 and −5.31, respectively), whereas the highest nighttime and daytime UHIIs were observed for LCZ7(2.20) and LCZ4(3.14), respectively. Figure  3 shows the differences in the daytime and nighttime UHII among the various LCZ types. The building types (LCZ1-LCZ10) generally exhibited higher UHII values, whereas the natural types (LCZA-LCZE) showed lower values. The average nighttime and daytime UHII values of the building types were 1.11 and 2.77°C higher than those of the natural types, respectively. In addition, LCZ types with different vegetation coverage also exhibited differences in their daytime and nighttime UHIIs. During the daytime, LCZ types with high vegetation coverage (LCZ10, A, B, and C) exhibited lower UHIIs than LCZ types with low vegetation coverage (LCZ3, 4, and E). At night, the UHII values of LCZ7 and 8 were higher than those of LCZ1 and 2, which had lower levels of vegetation coverage. This is because the impervious surface of the city differs from the surfaces with vegetation, the reflectivity of solar radiation is small, the heat capacity is large, and the nearsurface layer is relatively stable at night, resulting in poor heat dissipation, high stored heat, and a significant heat island effect.

Influence of sea and land breeze on UHII
Taking the offshore G3643 station as an example, typical (no significant changes in wind speed and direction) and atypical (diurnal changes in wind speed and direction) sea-land breeze days were analyzed, as shown in Fig. 4 (the typical sea-land breeze day was November 1, 2019, and the atypical sea-land breeze day was November 5, 2019). The results showed that the wind speed and direction were evidently disturbed by the onset of the sea breeze. The wind direction changed from offshore at the last moment to onshore, with an increase of >90°, and the wind speed increased rapidly. However, as the sea breeze developed, the wind direction stabilized. In the onshore wind direction, the wind speed reached its daytime maximum at~14:00 (Beijing: GMT + 8). When the sea breeze ended, the wind direction suddenly shifted offshore, and the land breeze started developing. The wind speed reached its intraday maximum at~03:00 at night (Beijing: GMT + 8). For atypical sea and land breeze days, the dominant wind direction was southeast, and the wind speed remained at 1 m·s −1 .
To compare the influence of the presence or absence of sea-land breeze on the intensity of heat islands, the daily variation characteristics of the heat island intensity on typical and atypical sea-land breeze days were analyzed by taking the offshore station (G3643) as an example. As shown in Fig. 5, the diurnal variation characteristics in the heat island intensity corresponding to the sea-land breeze day exhibited a "V"-shaped distribution. When the sea breeze started, the UHII reduced sharply, reaching its intraday minimum at~14:00 (Beijing: GMT + 8). This was considerably lower than that of the atypical sea-land breeze day. The UHII value of the land breeze day subsequently began to increase, stabilized at 21:00 (Beijing: GMT + 8), and then continued to increase, reaching its intraday maximum at~03:00 (Beijing: GMT + 8). This was higher than the UHII value on the atypical sea-land breeze day, but not on the typical sea-land breeze day. The UHII increased at 18:00 (Beijing: GMT + 8), reached its maximum at 23:00 (Beijing: GMT + 8), and gradually decreased before stabilizing. It began to decrease at 07:00 (Beijing: GMT + 8) and then remained relatively stable and low until 17:00 (Beijing: GMT + 8); however, it was higher than that on the typical sea-land breeze day. This showed that the sea breeze could cool the city  and alleviate the heat island effect, whereas the land breeze could moderately increase the strength of the heat island effect. Therefore, compared with the atypical sea-land breeze day, the daily UHI variations were greater on the typical sea-land breeze day.
Both sea and land breezes had a certain effect on UHI. To further analyze the mechanism between the sea-land breezes and UHI, we combined the UHIIs during the start and end of the sea and land breezes to obtain the correlations between the average wind speed of the overall sea and land breezes and UHIIs in Shenzhen (Fig. 6); subsequently, both passed the confidence test with a significance level of α = 0.05. The sea breeze wind speed was negatively correlated with the UHII, with a correlation coefficient of −0.68601; i.e., the greater the sea breeze wind speed, the smaller the UHII, indicating that the sea breeze alleviated the heat island effect. However, land breeze wind speed and UHII were positively correlated, with the correlation coefficient of 0.81142; i.e., the greater the wind speed, the greater the UHII. The R 2 values between the sea and land wind speeds and UHII were 0.60862 and 0.26253, respectively, indicating that the impact of sea breeze was stronger than that of a land breeze. As the sea breeze was accompanied by damp and cold air currents, the UHI effect was greatly reduced.

DISCUSSION
Most of the currently used LCZ classification schemes employ the WUDAPT L0 method with a spatial resolution of 100 m 69 or remote sensing image data (spatial resolution of 30 m-1 km), combined with building data and other GIS analysis methods 70,71 . For example, Chen et al. 70 selected training samples based on MODIS and ASTER images, combined with Google Earth high-definition images. They then used the random forest method to perform LCZ classification to explore the spatial distribution of thermal environments in Guangzhou and Hong Kong, two of China's subtropical high-density cities; consequently, the two cities were found to be clearly spatially resolved. Among these methods, the WUDAPT L0 approach appears to be more suitable for performing LCZ classification over large-scale areas; however, analyzing the internal differences in small-and medium-sized cities requires higher-precision data. Based on the GIS analysis, this study used land-use and urban building data, with a spatial resolution of 10 m, to classify the LCZs. This approach resolved the low accuracy and yielded more comprehensive and accurate classification results. Additionally, this study referred to the LCZ classification scheme of Stewart & Oke 40 and combined the actual situation in Shenzhen to classify the buildings more finely, which can better reflect the effect of building height and density on the strength of the heat island.
The LCZ scheme was originally designed to use the difference in air temperature to quantify the intensity of the heat island 39 . However, the distribution of meteorological stations is uneven, and the acquired UHII cannot cover all LCZs. The UHII value of the study area ignores the differences within the city, and it is not convenient to study the effect of LCZs on UHII. This study also calculated the air temperature results ( Table 3) and found that the UHII of the building type during the daytime was higher than that of the natural type. The difference between the building types was not obvious. Therefore, this article analyzed from the perspective of LST to reduce the influence of unknown interference factors that cause interference.
When analyzing the effect of sea and land breeze on UHII, this study used the air temperature values acquired from the meteorological stations. In addition, to illustrate the difference between air temperature and land surface temperature, the relationship of UHII with sea and land breezes was also calculated using land surface temperature (Fig. 7). The results (all passed the confidence test at the level of α = 0.05) showed that sea and the land breeze had a positive and negative correlation with UHII, respectively, but both the R 2 values were small (3.44127E-6 and 0.03529, respectively), and the fitting effect was not ideal, Fig. 6 The relationship between sea-land breeze wind speed and heat island intensity. a Land breeze, b sea breeze. indicating that the effect of sea and land breeze on UHII was not ideal from the perspective of LST (for Shenzhen).
Most previous studies have only focused on the effects of LCZs or sea-land breeze on UHIs [72][73][74] . Furthermore, the influencing mechanism of the UHII has rarely been investigated from the comprehensive perspective of the sea-land breeze and LCZ. Therefore, Zhou et al. 75 and Martinelli et al. 76 . discussed the impacts of different LCZ types on the surface temperature from the perspective of LCZs by simultaneously qualitatively analyzing the impacts of sea and land breezes on urban temperatures in Sendai, Japan, and Bari, Italy, respectively. However, Sendai, Bari, and Shenzhen are characterized by different climate types, and differences in the LCZ distributions within these cities. In addition, they did not quantitatively analyze the effects of sea and land breezes on the UHII. Moreover, when analyzing the influences of sea and land breezes on UHII, they focused less on the differences between these two variables under different gradients. Therefore, based on the two directions of the sea and land breezes and LCZs, this study analyzed the impacts of different LCZ types on the UHII while also determining their correlation with UHII. The average daytime and nighttime UHII values of the building LCZ types were 1.11 and 2.77°C higher than those of the natural types, respectively. Sea breezes alleviated the UHI effect, whereas land wind moderately enhanced the UHI effect. The linear regression coefficients between the sea and land breeze wind speeds and UHII were −0.68601 and 0.81142, respectively.
To further analyze the variations in the sea and land wind speeds under different gradients and the influences of different LCZ types on the UHI strength, ArcMap 10.4 (Environmental Systems Research Institute, California, USA) was used to calculate the distances of various meteorological stations from the coastline and analyze the effects of the sea and land wind speeds and LCZ types on the UHII (UHII is obtained by air temperature) from different gradient angles (Fig. 8). Regardless of the type of LCZ at the meteorological station, the overall trends in the sea breeze wind speed and UHII values were inversely proportional, i.e., the greater the sea breeze wind speed, the smaller the UHII. Moreover, for the same LCZ type, areas farther from the coastline had lower sea breeze wind speeds and greater UHII values. For example, stations G3643, G3546, G3530, and G3634 were all LCZ4 type, and their distances to the coastline were 1.9, 2.6, 2.7, and 12.7 km, respectively. The sea breeze wind speeds were G3643 (3.55 m·s −1 ) Fig. 7 The relationship between sea-land breeze wind speed and heat island intensity. a Land breeze, b sea breeze. Fig. 8 The relationship between day and night UHII and sea-land wind speed, LCZ types, and distance to coastline at various weather stations. SID-LCZ indicates the station number and the type of LCZ to which it belongs, DIS stands for the distance between the station and the coastline, SWS and LWS represent wind speed of sea and land, respectively, DUHII and NUHII, respectively, represent the intensity of day and night heat island.
The urban spatial distribution of Shenzhen differs from that of Sendai (Japan), Tianjin, and other cities in China. Its urban center is near the coastline, and the suburbs are inland. Therefore, the conclusions drawn here may differ from those of other cities. The differences in the climate types may also cause differences in the results. Thus, local conditions should be taken into consideration when performing LCZ divisions and calculating the correlations between sea and land breezes and the UHII.
This study primarily analyzed the effects of different LCZ types, typical and atypical sea-land breeze days, and different gradient sea-land wind speeds on daytime and nighttime UHII. Our results provide an important reference for urban planning and government decision-making. However, this study also has certain limitations. First, as the AST_08 data do not contain water temperature information, the potential factors that can affect the strengths of UHIs were not taken into consideration. Second, due to time constraints associated with image acquisition, the daytime and nighttime UHIIs recorded in this study were not obtained for the same days. Although they were acquired under clear weather conditions within the same month, various daily conditions may have caused partial differences in the results. Finally, this study only considered the impacts of the sea and land breezes on the daytime and nighttime UHIIs only in November, and seasonal differences or the UHII mechanisms were not considered. Therefore, the seasonal and regional environmental differences in the sea and land breezes and daytime and nighttime UHIIs should be further explored.

Study area and data
Shenzhen (Fig. 9), also known as "Pengcheng," is located in southern Guangdong, China, on the eastern bank of the Pearl River Estuary, with Daya Bay and Dapeng Bay to the east and the Pearl River Estuary and Lingding Ocean to the west. The terrain is high in the southeast and low in the northwest. Most of the area comprises low hills with gentle terraces. The region has a subtropical oceanic climate. Due to the strong influence of monsoons, the dominant wind direction is easterly-tosoutheast, and southeasterly winds prevail in the summer. There are occasional monsoon lows and tropical cyclones. The northeast monsoon prevails during the remaining seasons. The weather is relatively dry, the climate is mild, and the annual average temperature is~22.4°C.
Data on land use, construction, administrative divisions, digital elevation models, Landsat remote sensing images, ASTER surface temperatures, and meteorological factors were used in this study (Table 4).

LCZ classification
An LCZ refers to a combination of similar thermal environment characteristics based on urban surface properties and morphologies. However, the LCZ classification system and its standards are not static. This study used the results of previous studies 39,40,[50][51][52] , along with the urban structure and architectural characteristics of Shenzhen, to establish its LCZ classification system. Thus, an LCZ system, comprising 16 categories, was constructed, which contained 10 buildings and six natural types (Table 5. Furthermore, the study area was divided into 30 × 30 m grids using the fishing net tool in ArcGIS 10.4 (Environmental Systems Research Institute, California, USA). The building data were then mapped to the fishing net. Moreover, the building-type areas were divided according to two morphological indicators: building height and building density. The classification results were labeled under LCZ1-LCZ10. Among these, "dense" referred to a building density of >0.4, while "open" referred to a building density of <0.4. In addition, the natural type areas were divided using land-use data. The corresponding classification results were labeled as LCZA-LCZF.

Sea and land wind speed
A sea-land breeze is a small-to-medium-scale thermal circulation caused by the temperature difference between the sea and land. It typically overlaps with the background wind field. Along the coastline of Shenzhen, the sea-land wind direction and monsoon wind direction overlap each other, which can easily lead to errors in assessing the sea-land winds. Therefore, effectively distinguishing the sea-land winds is crucial. This study referred to the sea and land wind distinction method proposed by Wei 77 . Thus, the wind was decomposed into a vector based on a trigonometric function, i.e., u component in the east-west direction and v component in the north-south direction. The corresponding formulas are as follows: where V is the measured wind speed at the weather station and D is the wind direction. The winds were divided into the measured, system, and local winds. Among these, the measured wind comprised hourly data recorded by various meteorological stations, system wind comprised large-scale background wind recorded daily at each station, and local wind comprised mesoscale circulation, such as sea-land and valley wind. After the vector decomposition of the measured wind values, the 24 h average values of the u and v components were calculated. The daily average wind was obtained after synthesis, i.e., the system wind of a particular day. The measured and system wind values were calculated as the vector difference, i.e., the local wind. As this study did not consider the influence of valley wind, the calculated local wind comprised the sea and land wind. The specific formula is as follows: and Further, in this study, data of November 2019 acquired from the meteorological stations in the Futian, Bao'an, Longhua, and Nanshan districts were analyzed (except for the adjacent mountains). Based on statistics, 02:00-09:00 (Beijing: GMT + 8) was selected as the duration for land winds, while 14:00-21:00 (Beijing: GMT + 8) time was selected as the duration for sea winds; subsequently, the average wind speeds for the sea and land winds during these times were calculated.

UHII
In addition to the LCZ theory, the definition of the heat island effect changed. When Stewart and Oke 40 proposed the LCZ theory, they also redefined the intensity of the heat island effect. This study adopted the redefinition of the intensity of the heat island effect proposed by Stewart and Oke 40 . The formula is as follows: where UHII LCZX represents the heat island effect intensity of LCZX and T LCZX and T LCZD are the surface temperature of the type X and D LCZs,  respectively. In other words, the intensity of the heat island effect is the temperature difference between each LCZ and the LCZD type (low vegetation).

DATA AVAILABILITY
The datasets aggregated and/or analyzed during the current study are available from the corresponding author on reasonable request.