Mesoscale Observational Analysis of Isolated Convection Associated with the Interaction of the Sea Breeze Front and the Gust Front in the Context of the Urban Heat Humid Island Effect

: An isolated convection was unexpectedly initiated in the evening of 1 August 2019 around the Tianjin urban region (TUR), which happened at some distance from the shear line at lower level and the preexisting convection to the South, analyzed by using ERA5 reanalysis data and observations from surface weather stations, and a S-band radar. The results show that, 42 min before the initiation of the convection, the atmospheric thermodynamic conditions around TUR were favorable for the initiation of the isolated convection, although the southerly and vertical shear of the horizontal wind at the lower level was weak. A sea-breeze front approached the TUR and continued to move West, leading to the triggering of the isolated convection in the context of the urban humid heat island (UHHI) effect. Subsequently, the gust front, which was formed between the cold pool away from the TUR and the warm and humid air of the UHHI, moved northward, approached the convection, and collided with sea breeze front, resulting in ﬁve reﬂectivity centers of isolated convection being merged and the convection’s development. Finally, the isolated convection split into two convections that moved away from the TUR and disappeared at 20:36 Beijing Time. The isolated convection was initiated and developed by the interaction of the sea breeze front and gust front in the context of the UHHI effect. The sea breeze front triggered the isolated convection around TUR in the context of the UHHI effect, and the gust front produced by the early convective storms to the south played a vital role in the development of the isolated convection.


Introduction
Rainstorms in North China, especially in the West Bank of the Bohai Sea, are characterized by their short duration, rapid generation and high intensity and being strongly localized. Most of local rainstorms are caused by small-scale or meso-scale weather systems, which are very difficult to predict. The research has showed that the convergence caused by the local terrain provides better dynamic conditions for the initiation of local storms [1,2]. The interaction between the mesoscale cold pool and easterly winds [3,4], and the triggering effect of small-scale cyclones [3] on local heavy rainfall clouds has also been studied. In addition, scholars have also studied the role of the gust front [3], dew point front [5] and sea breeze front [6] in the triggering and development of local convections. Undoubtedly, these studies have significantly advanced the understanding of rainfall-generated mesoscale convective storms. Despite the many advances in the understanding of the structure and life cycle of convections, the accurate forecasting of convective initiation remains one of the greatest challenges in meteorology [7][8][9][10][11], which is especially true for the convective initiation (CI) of isolated convections [12]. From 2018 to 2021, the initiation and development of isolated convection occurred many times around the Tianjin urban region (TUR) on the West Bank of the Bohai Sea, such as the convections in the afternoon of 22 July 2018, the afternoon of 18 July 2019, the evening of 1 August 2019 and the evening of 7 July 2021. The multi-scale mechanism of isolated convection around the TUR in the afternoon of 18 July 2019 was analyzed [13], and the results showed that the initiation of the isolated convection was closely related to the urban heat island (UHI) phenomenon. The resulting UHI and urban canopy effects may lead to changes in the local and regional atmospheric circulation and precipitation [14][15][16][17][18][19][20]. The analysis of 113 locally developed non-typhoon cyclone extreme rainfall events in the summer of 2011-2018 also shows the contribution of the UHI effect to the occurrence of more hourly extreme rainfall [21].
As everyone knows, the UHI effect is the thermal mechanism of isolated convection around urban regions; however, the exact initial time and location of isolated convective systems remain difficult to predict [12]. Previous studies have shown that the convective cloud is mainly initiated from the convergence zone in the boundary layer, such as the zone near the front and the dry line [5,6,[22][23][24][25][26][27][28][29][30][31]. In addition, the sea breeze front is also one of the significant mesoscale weather systems triggering convective storms in Tianjin, as a coastal megacity in North China, and it usually causes localized storms as it moves further inland [32,33]. A lot of researches have been performed on convections caused by sea breeze fronts [6,[34][35][36][37][38][39]. A three-dimensional, cloud-resolving model was used to investigate the interaction between the sea-breeze circulation and boundary layer roll convection [40,41], and the results provide further evidence that horizontal convective rolls can play an important role in the initiation and modulation of convection along a sea-breeze front. Large-eddy simulations were performed to investigate the process of deep-convection initiation over a peninsula, and the result showed that convection is produced through the collision of the sea-breeze fronts or the collision between the gust front and sea breeze front [42]. In addition, the results of a lot of researches show that sea breeze front can not only interact with dynamic systems, such as sea breeze front [42], gust front [6,35] and horizontal convective roll [40,41], but also interact with thermal systems, such as the UHI effect [13,[43][44][45][46][47].
Nevertheless, few studies have been performed to reveal the interaction of the sea breeze front and gust front in the context of the urban humid heat island (UHHI) effect, which led to the generation and development of the isolated convection around the TUR. Thus, the present study focuses on the isolated convection around the TUR on 1 August 2019. This paper aims to better understand the role of the sea breeze front, gust front and the UHHI effect in the initiation and development of isolated convection. Moreover, this research also aims to acquire a better understanding about the physical mechanisms of this type of weather case, thus improving the ability of forecasting and early warning.

Dataset
The fifth generation of the European Center for Medium-Range Weather Forecasts (ECMWF) reanalysis database (ERA5) was used, which provides atmospheric variables at 1 h intervals. The spatial resolution of this dataset is 0.25 • × 0.25 • . This dataset has 37 levels in vertical direction [48]. The ERA5 data are available at https://cds.climate.copernicus. eu/#!/search?text=ERA5&type=dataset (accessed on 8 April 2020).
A total of 289 surface weather stations were used in this study, including 11 national and 278 regional observation stations, whose distribution is shown in Figure 1. These surface stations are relatively evenly distributed with an average distance of about 5 km. The meteorological elements observed by surface stations include precipitation, wind, temperature and dew point temperature at 1 h intervals. The data were interpolated into 0.1 • × 0.1 • grids by the inverse distance weight method in order to plot the contour map and the color map. The meteorological elements observed by surface stations include precipitation, wind, temperature and dew point temperature at 1 h intervals. The data were interpolated into 0.1° × 0.1° grids by the inverse distance weight method in order to plot the contour map and the color map.
Additionally, the Tanggu Doppler S-band Radar was also used in this paper. The location of this radar (117.72° E, 39.04° N) is shown in Figure 1.

Definition of Urban Humid Heat Island
As it is well known, the "urban heat island" effect refers to the "high temperature" of the city caused by a large number of artificial heating, high heat storage bodies, such as buildings and roads, and the reduction in green spaces. The temperature in the urban region is significantly higher than that in the rural region. However, another phenomenon was identified in our research, which is also worthy of our attention, that is, that the city is located in a high temperature zone and high humidity zone (high dew point temperature or high specific humidity zone), and is neither located in the center of high temperature nor in the center of high humidity, but is located in the center of high pseudo-equivalent potential temperature, indicating that the urban region has the characteristics of high temperature and high humidity compared with the surrounding rural area. In this paper, this phenomenon was defined as the phenomenon of urban humid heat island (UHHI).

Selection of Standard Time
The time of Eastern Eighth Time Zone was defined Beijing Time (BT). BT is 8 h earlier than universal time, that is, BT equals universal time plus 8 h. In this paper, Beijing time was taken as the standard time to describe the convective events.

Methodology of Creating Sounding
Due to their higher horizontal and temporal resolutions, reanalysis soundings have been used for climatological storm environment studies by a growing number of studies [49][50][51][52][53][54]. Therefore, the present study created proximity soundings for this case, which included creating reanalysis soundings using the ERA5 reanalysis dataset and correcting Additionally, the Tanggu Doppler S-band Radar was also used in this paper. The location of this radar (117.72 • E, 39.04 • N) is shown in Figure 1.

Definition of Urban Humid Heat Island
As it is well known, the "urban heat island" effect refers to the "high temperature" of the city caused by a large number of artificial heating, high heat storage bodies, such as buildings and roads, and the reduction in green spaces. The temperature in the urban region is significantly higher than that in the rural region. However, another phenomenon was identified in our research, which is also worthy of our attention, that is, that the city is located in a high temperature zone and high humidity zone (high dew point temperature or high specific humidity zone), and is neither located in the center of high temperature nor in the center of high humidity, but is located in the center of high pseudo-equivalent potential temperature, indicating that the urban region has the characteristics of high temperature and high humidity compared with the surrounding rural area. In this paper, this phenomenon was defined as the phenomenon of urban humid heat island (UHHI).

Selection of Standard Time
The time of Eastern Eighth Time Zone was defined Beijing Time (BT). BT is 8 h earlier than universal time, that is, BT equals universal time plus 8 h. In this paper, Beijing time was taken as the standard time to describe the convective events.

Methodology of Creating Sounding
Due to their higher horizontal and temporal resolutions, reanalysis soundings have been used for climatological storm environment studies by a growing number of studies [49][50][51][52][53][54]. Therefore, the present study created proximity soundings for this case, which included creating reanalysis soundings using the ERA5 reanalysis dataset and correcting the surface reanalysis soundings using 1 h surface observations. The construction of the sounding was as follows. First, 27 isobaric layers from 1000 hPa to 100 hPa of pressure, altitude, temperature, dew-point temperature, and horizontal wind at the certain grid (39 • N, 117 • E) and time (14:00 BT and 17:00 BT) were extracted from the ERA5 data. Next, the observation data of the surface station in the TUR were used to take the place of the lowest level data of ERA5 sounding.  Figure 2 displays the distribution of the accumulated rainfall amounts during this weather event, which was from 17:00 to 20:00 BT on 1 August 2019. The short-time strong rainfall event of interest occurred around the TUR on 1 August 2019. The 3 h (17:00 BT-20:00 BT) accumulated rainfall exhibited an isolated precipitation (the maximum rainfall intensity of 36.1 mm·h −1 located in Tiaoyuanli in the TUR during 18:00-19:00 BT) area, characterized by rapid generation, short duration, high strength and strong locality. Rainfall amounts in 9 stations exceeded 25 mm, while the maximum rainfall amounts of 43.8 mm appeared in the Erweilu weather station of the TUR. This type of local rainstorm in Tianjin in summer is not rare (there were four isolated convection events around the TUR from 2018 to 2021, and another three occurred at noon on 22 July 2018, in the afternoon on 18 July 2019, and in the evening on 7 July 2021), which also brings challenges to early and accurate warning. the surface reanalysis soundings using 1 h surface observations. The construction of the sounding was as follows. First, 27 isobaric layers from 1000 hPa to 100 hPa of pressure, altitude, temperature, dew-point temperature, and horizontal wind at the certain grid (39° N, 117° E) and time (14:00 BT and 17:00 BT) were extracted from the ERA5 data. Next, the observation data of the surface station in the TUR were used to take the place of the lowest level data of ERA5 sounding. Figure 2 displays the distribution of the accumulated rainfall amounts during this weather event, which was from 17:00 to 20:00 BT on 1 August 2019. The short-time strong rainfall event of interest occurred around the TUR on 1 August 2019. The 3 h (17:00 BT-20:00 BT) accumulated rainfall exhibited an isolated precipitation (the maximum rainfall intensity of 36.1 mm•h −1 located in Tiaoyuanli in the TUR during 18:00-19:00 BT) area, characterized by rapid generation, short duration, high strength and strong locality. Rainfall amounts in 9 stations exceeded 25 mm, while the maximum rainfall amounts of 43.8 mm appeared in the Erweilu weather station of the TUR. This type of local rainstorm in Tianjin in summer is not rare (there were four isolated convection events around the TUR from 2018 to 2021, and another three occurred at noon on 22 July 2018, in the afternoon on 18 July 2019, and in the evening on 7 July 2021), which also brings challenges to early and accurate warning.   Figure 3, the accumulated precipitation around the TUR was contributed by the locally developed, isolated convection during 17:40-20:00 BT on 1 August. At 17:30 BT on 1 August (see Figure 3a), a mature convective storm moving from East to West was located to the south of Tianjin, with a gust front in front of its movement  Figure 3 displays the evolution of the rainy storm observed by the Tanggu S-band Doppler Radar. As shown in Figure 3, the accumulated precipitation around the TUR was contributed by the locally developed, isolated convection during 17:40-20:00 BT on 1 August. At 17:30 BT on 1 August (see Figure 3a), a mature convective storm moving from East to West was located to the south of Tianjin, with a gust front in front of its movement and there was no convection around the TUR. The isolated convection was initiated at the junction of the TUR, XQD and JND at 17:42 BT (see Figure 3b), and the center of reflectivity was greater than 35 dBZ, reaching levels above 40 dBZ at 17:48, and located in the XQD (see Figure 3c). The continuous convection was triggered and the subsequent growth of the convections led to the convection having five centers of reflectivity, which were located in the south of the TUR, but the highest reflectivity of the convections was still 40 dBZ at 18:00 BT (see Figure 3d). grow into its mature stage, whose highest reflectivity reached levels of more than 55 dBz, locating at the junction of the TUR and XQD. The isolated convection moved northward from 18:12 to 18:30 BT, and the main body of the isolated convection influenced the TUR at 18:30 BT (see Figure 3g). After 19:00 BT (see Figure 3h), the isolated convection around the TUR rapidly weakened, dissipated, moved northward, and was gradually split into two convections, one (convection A) in the North of the TUR and the other (convection B) in the southeast of the TUR. After 20:00 BT (see Figure 3i), convection A moved northward, and convection B moved northeast by East. Then, the two convections gradually weakened and disappeared at 20:36 BT. and there was no convection around the TUR. The isolated convection was initiated at the junction of the TUR, XQD and JND at 17:42 BT (see Figure 3b), and the center of reflectivity was greater than 35 dBZ, reaching levels above 40 dBZ at 17:48, and located in the XQD (see Figure 3c). The continuous convection was triggered and the subsequent growth of the convections led to the convection having five centers of reflectivity, which were located in the south of the TUR, but the highest reflectivity of the convections was still 40 dBZ at 18:00 BT (see Figure 3d). From 18:06 BT (see Figure 3e) to 18:12 BT (see Figure 3f), the merging of the five reflectivity centers in the South of the TUR made the isolated convection grow into its mature stage, whose highest reflectivity reached levels of more than 55 dBz, locating at the junction of the TUR and XQD. The isolated convection moved northward from 18:12 to 18:30 BT, and the main body of the isolated convection influenced the TUR at 18:30 BT (see Figure 3g). After 19:00 BT (see Figure 3h), the isolated convection around the TUR rapidly weakened, dissipated, moved northward, and was gradually split into two convections, one (convection A) in the North of the TUR and the other (convection B) in the southeast of the TUR. After 20:00 BT (see Figure 3i), convection A moved northward, and convection B moved northeast by East. Then, the two convections gradually weakened and disappeared at 20:36 BT.

Synoptic Background
The major features of the synoptic conditions in the lower and middle troposphere (from 1000 hPa to 500 hPa) before the isolated convection was formed are displayed in Figure 4. The yellow pentagram in Figure 4 is used to indicate the position of the TUR. At 14:00 BT on 1 August, the TUR was located in the area of high pseudo-equivalent potential temperature with a weak southerly flow and about 250 km away from a 925 hPa convergence line to the north. The convergence line (red dashed line in Figure 4a,b) remained at about 43 • N~44 • N (see Figure 4a). Due to the imbalance of wind fields on the north and south sides of the western part of the convergence line at this time, 3 h later (42 min before the initiation of the convection, 17:00 BT), the western part of the convergence line was obviously moving southward (from 43 • N to 41 • N), but the eastern part was moving very slowly, so the distance between the convergence line and the TUR had no significant change, which suggests that the boundary of the low-level convergence line or baroclinic zone never intruded the TUR, while the isolated convection was initiated. In addition, as shown in Figure 4c,d, 3 h before the initiation of the isolated convection, the atmosphere of the middle level and low level were always in a convective unstable condition (pseudo-equivalent potential temperature decreased with height), and 42 min before the initiation of the isolated convection, there was a high pseudo-equivalent potential temperature center greater than 360 K in the lower layer of the TUR, increasing the vertical gradient of the pseudo-equivalent potential temperature over the TUR, which was conducive to the initiation and development of the convection. potential temperature with a weak southerly flow and about 250 km away from a 92 hPa convergence line to the north. The convergence line (red dashed line in Figure 4a, remained at about 43° N~44° N (see Figure 4a). Due to the imbalance of wind fields o the north and south sides of the western part of the convergence line at this time, 3 h lat (42 min before the initiation of the convection, 17:00 BT), the western part of th convergence line was obviously moving southward (from 43° N to 41° N), but th eastern part was moving very slowly, so the distance between the convergence line an the TUR had no significant change, which suggests that the boundary of the low-lev convergence line or baroclinic zone never intruded the TUR, while the isolate convection was initiated. In addition, as shown in Figure 4c,d, 3 h before the initiation the isolated convection, the atmosphere of the middle level and low level were always a convective unstable condition (pseudo-equivalent potential temperature decrease with height), and 42 min before the initiation of the isolated convection, there was a hig pseudo-equivalent potential temperature center greater than 360 K in the lower layer the TUR, increasing the vertical gradient of the pseudo-equivalent potential temperatu over the TUR, which was conducive to the initiation and development of the convection  There is no sounding station in Tianjin City. When analyzing and researching the strong convection case, Beijing or Laoting sounding closing to Tianjin are usually used. However, in fact, these two soundings are more than 100 km away from Tianjin, which is difficult to really reflect the atmospheric structure over Tianjin. Therefore, the hourly temperature, wind, humidity and pressure of surface meteorology observation stations in the TUR were combined with the ERA5 data of the adjacent grid points (117 • E, 39 • N) above the boundary layer to form the hourly sounding data (see Figure 5), which was analyzed to further investigate the dynamic and thermodynamic conditions for the initiation and development of the isolated convection in the TUR. Three hours before the initiation of the isolated convection (see Figure 5a), the temperature lapse rate in the lower troposphere over the TUR was higher than the dry adiabatic lapse rate, which means that the atmosphere in the TUR was absolutely unstable. At the same time, the convective inhibition (CIN) was 0 J/kg, the lifting condensation level (LCL) was only 827.5 hPa, and the convective available potential energy (CAPE) reached 2375 J/kg, which suggests that the atmospheric conditions around the TUR were greatly beneficial to convective triggering with very little lifting force required. The precipitable water was 43.76 mm, which means that the water vapor condition over the TUR was not bad. During the initiation period of the isolated convection (see Figure 5b), the atmosphere in the TUR was still absolutely unstable (CIN = 0 J/kg). In addition, the CAPE was further increased to 3771 J/kg, and the precipitable water increased to 47.11 mm, indicating that the atmospheric conditions were quite favorable to convective triggering and maintenance. Significantly, it can be seen from Figure 5 that the vertical shear of wind was relatively weak at a lower troposphere, barely reaching 4 m·s −1 between 1000 hPa to 850 hPa before the initiation of the isolated convection.
difficult to really reflect the atmospheric structure over Tianjin. Therefore, the hour temperature, wind, humidity and pressure of surface meteorology observation station in the TUR were combined with the ERA5 data of the adjacent grid points (117° E, 39° N above the boundary layer to form the hourly sounding data (see Figure 5), which w analyzed to further investigate the dynamic and thermodynamic conditions for the init ation and development of the isolated convection in the TUR. Three hours before the in tiation of the isolated convection (see Figure 5a), the temperature lapse rate in the low troposphere over the TUR was higher than the dry adiabatic lapse rate, which means th the atmosphere in the TUR was absolutely unstable. At the same time, the convectiv inhibition (CIN) was 0 J/kg, the lifting condensation level (LCL) was only 827.5 hPa, an the convective available potential energy (CAPE) reached 2375 J/kg, which suggests th the atmospheric conditions around the TUR were greatly beneficial to convective tri gering with very little lifting force required. The precipitable water was 43.76 mm, whic means that the water vapor condition over the TUR was not bad. During the initiatio period of the isolated convection (see Figure 5b), the atmosphere in the TUR was st absolutely unstable (CIN = 0 J/kg). In addition, the CAPE was further increased to 377 J/kg, and the precipitable water increased to 47.11 mm, indicating that the atmospher conditions were quite favorable to convective triggering and maintenance. Significantl it can be seen from Figure 5 that the vertical shear of wind was relatively weak at a low troposphere, barely reaching 4 m•s −1 between 1000 hPa to 850 hPa before the initiation the isolated convection. In summary, although the wind speed of airflow and the vertical shear of horizont wind were relatively weak, the atmospheric thermodynamic conditions promoted th start-up of the isolated convection in the TUR. However, because of the absence of a si nificant synoptic scale lifting mechanism in the TUR, and the initiation of the isolate convection cannot be related to the dynamic system, we hypothesize that some sma or mesoscale weather system might be responsible for the precise location and timing the isolated convection occurrences. Based on the analysis of the surface meteorologic elements seen by surface observations and the reflectivity detected by the Doppler R dar, we demonstrate in the following section that the sea breeze front and gust front we In summary, although the wind speed of airflow and the vertical shear of horizontal wind were relatively weak, the atmospheric thermodynamic conditions promoted the start-up of the isolated convection in the TUR. However, because of the absence of a significant synoptic scale lifting mechanism in the TUR, and the initiation of the isolated convection cannot be related to the dynamic system, we hypothesize that some small or mesoscale weather system might be responsible for the precise location and timing of the isolated convection occurrences. Based on the analysis of the surface meteorological elements seen by surface observations and the reflectivity detected by the Doppler Radar, we demonstrate in the following section that the sea breeze front and gust front were key to the occurrence and development of the isolated convection around the TUR in the context of the UHHI effect.  Figure 6b shows an enlarged view of the red dashed box in Figure 6a. As shown in Figure 6a, at 17:00 BT, the gust front (blue dashed line) was formed and its formation was related to the previous convections, which was 20-30 km away from the TUR. Meanwhile, there was velocity convergence in easterly (sea-breeze) front around the TUR at 17:00 BT, where the isolated convection was observed at 17:42 (see Figure 6b). That is to say, the sea-breeze front was the key mesoscale system, which was conductive to the isolated convection's initiation there. Figure 6c shows the composite reflectivity and surface wind at 18:00 BT on 1 August 2019, and Figure 6d shows an enlarged view of the blue dashed box in Figure 6c. As is shown in Figure 6c,d, at 18:00 BT, the gust front moved northward, and the sea breeze front moved westward. As a result, the gust front influenced the isolated convection and collided with the sea breeze front around the TUR, making five reflectivity centers of convection merge and strengthen rapidly. At 18:00 BT, there were five Atmosphere 2022, 13, 603 8 of 15 echo centers greater than 40 dBz, and by 18:30, they had merged into an isolated convection with an echo center greater than 55 dBz (see Figure 3g). Figure 6b). That is to say, the sea-breeze front was the key mesoscale system, which wa conductive to the isolated convection's initiation there. Figure 6c shows the composi reflectivity and surface wind at 18:00 BT on 1 August 2019, and Figure 6d shows an en larged view of the blue dashed box in Figure 6c. As is shown in Figure 6c,d, at 18:00 BT the gust front moved northward, and the sea breeze front moved westward. As a resul the gust front influenced the isolated convection and collided with the sea breeze fron around the TUR, making five reflectivity centers of convection merge and strengthe rapidly. At 18:00 BT, there were five echo centers greater than 40 dBz, and by 18:30, the had merged into an isolated convection with an echo center greater than 55 dBz (se Figure 3g). Figure 7a-f are the distributions of surface temperature, dew-point temperature and pseudo-equivalent potential temperature. As shown in Figure 7a-c, at 17:00 BT, due to the evaporative cooling caused by previous convections and clouds shielding from solar heating, the temperature fields at the surface exhibited a quasi-circular area of lower temperature (the center below 26 • C) to the South of Tianjin and there were also low value centers of dew point temperature (see Figure 7b) and pseudo-equivalent temperature (see Figure 7c) in close proximity to the same location. At the same time, there was a high temperature area extending from West to East at nearly 39 • N, with a high value center of more than 34 degrees, which is located to the west of the TUR, and the TUR is located at the east edge of the high temperature area (see Figure 7a). As shown in Figure 7b, there was a large dew point temperature area extending from southeast to northwest, and the TUR was also located at northwest edge of the high dew point temperature area. As a result, the TUR was located at the intersection of the high temperature area and high dew point temperature area. That is to say, the thermal characteristics of the TUR were high temperature and high humidity. As shown in Figure 7c, the pseudo-equivalent potential temperature center was located around the TUR, with a central value of 383 K, which can be defined as the UHHI effect, while the low value center of the pseudo-equivalent potential temperature was formed by early precipitation to the south of Tianjin, with a central value lower than 350 K, forming a strong thermal gradient and providing better thermodynamic conditions for the initiation of the isolated convection. As shown in Figure 7d-f, the thermal Atmosphere 2022, 13, 603 9 of 15 conditions at 18:00 BT were basically the same as those at 17:00 BT, while the cold pool to the south of Tianjin was more powerful. Corresponding to the further northward advance of the gust front, the high center of pseudo-equivalent potential temperature moved slightly northward, but it was still around the TUR, providing the thermal conditions for the further development and strengthening of the isolated convection. This suggests that, in addition to the gust front at the boundary layer, the baroclinic zone also moved northward during the development of the isolated convection. As shown in Figure 4c,d, the pseudo-equivalent potential temperature decreased with height, which means the atmosphere of the middle level and low level were always in a convective unstable condition. That is to say, the UHHI effect, corresponding to the high center of pseudo-equivalent potential temperature at the surface, strengthened the local instability around the TUR. So, the atmospheric thermodynamic conditions around the TUR were more conducive to the initiation of convection than those in rural areas.

Initiation of the Isolated Convection
high temperature and high humidity. As shown in Figure 7c, the pseudo-equivalent potential temperature center was located around the TUR, with a central value of 383 K, which can be defined as the UHHI effect, while the low value center of the pseudo-equivalent potential temperature was formed by early precipitation to the south of Tianjin, with a central value lower than 350 K, forming a strong thermal gradient and providing better thermodynamic conditions for the initiation of the isolated convection. As shown in Figure 7d-f, the thermal conditions at 18:00 BT were basically the same as those at 17:00 BT, while the cold pool to the south of Tianjin was more powerful. Corresponding to the further northward advance of the gust front, the high center of pseudo-equivalent potential temperature moved slightly northward, but it was still around the TUR, providing the thermal conditions for the further development and strengthening of the isolated convection. This suggests that, in addition to the gust front at the boundary layer, the baroclinic zone also moved northward during the development of the isolated convection. As shown in Figure 4c,d, the pseudo-equivalent potential temperature decreased with height, which means the atmosphere of the middle level and low level were always in a convective unstable condition. That is to say, the UHHI effect, corresponding to the high center of pseudo-equivalent potential temperature at the surface, strengthened the local instability around the TUR. So, the atmospheric thermodynamic conditions around the TUR were more conducive to the initiation of convection than those in rural areas. Thus, it can be seen that the isolated convection around the TUR that happened in the evening of 1 August 2019 was initiated and developed due to the interaction of the sea breeze front and gust front in the context of the UHHI effect. The sea breeze front approached the TUR and triggered the convection around the TUR in the context of the UHHI effect. The gust front produced by the early heavy precipitation played a vital role in the merging of the convection's five echo centers and the development of the isolated convection. Figure 8 is the composite reflectivity and surface wind at 19:00 BT on 1 August 2019. As shown in Figure 8, at 19:00 BT, the reflectivity center of the isolated convection (convection A) was located in the west of the TUR, and the value of reflectivity reached more than 55 dBz. Corresponding to the center of reflectivity, there was a divergent center of the surface wind field, and its eastward divergent air flow had shear with the environmental air flow (see Figure 8b), leading to a new convection (convection B) being triggered in the east of the TUR at 19:06 BT (see Figure 9a), and the reflectivity center reaching more than 50 dBz. As shown in Figure 9b Thus, it can be seen that the isolated convection around the TUR that happened in the evening of 1 August 2019 was initiated and developed due to the interaction of the sea breeze front and gust front in the context of the UHHI effect. The sea breeze front approached the TUR and triggered the convection around the TUR in the context of the UHHI effect. The gust front produced by the early heavy precipitation played a vital role in the merging of the convection's five echo centers and the development of the isolated convection. Figure 8 is the composite reflectivity and surface wind at 19:00 BT on 1 August 2019. As shown in Figure 8, at 19:00 BT, the reflectivity center of the isolated convection (convection A) was located in the west of the TUR, and the value of reflectivity reached more than 55 dBz. Corresponding to the center of reflectivity, there was a divergent center of the surface wind field, and its eastward divergent air flow had shear with the environmental air flow (see Figure 8b), leading to a new convection (convection B) being triggered in the east of the TUR at 19:06 BT (see Figure 9a), and the reflectivity center reaching more than 50 dBz. As shown in Figure 9b-d, between 19:12 BT and 19:24 BT, the intensity and area of convection B increased continuously, while convection A moved to the northeast, whose intensity and area decreased gradually. convection. Figure 8 is the composite reflectivity and surface wind at 19:00 BT on 1 August 201 As shown in Figure 8, at 19:00 BT, the reflectivity center of the isolated convectio (convection A) was located in the west of the TUR, and the value of reflectivity reache more than 55 dBz. Corresponding to the center of reflectivity, there was a diverge center of the surface wind field, and its eastward divergent air flow had shear with th environmental air flow (see Figure 8b), leading to a new convection (convection B) bein triggered in the east of the TUR at 19:06 BT (see Figure 9a), and the reflectivity cent reaching more than 50 dBz. As shown in Figure 9b-d, between 19:12 BT and 19:24 BT, th intensity and area of convection B increased continuously, while convection A moved the northeast, whose intensity and area decreased gradually.   Figure 10b), convection B continued to strengthen and the echo center was at the height of 1 km (the lowest height at which the echo can be detected by the Tanggu Radar at this position). After that (see Figure 10c,d), convection B continued to strengthen, the echo center was still at the height of 1 km, and reflectivity of the echo center increased to more than 50 dBz. However, convection A continued to weaken. 9a-d. As shown in Figure 10, at 19:06 BT (see Figure 10a), convection B was formed. The echo center was at the height of 3-6 km. At 19:12 BT (see Figure 10b), convection B continued to strengthen and the echo center was at the height of 1 km (the lowest height at which the echo can be detected by the Tanggu Radar at this position). After that (see Figure 10c,d), convection B continued to strengthen, the echo center was still at the height of 1 km, and reflectivity of the echo center increased to more than 50 dBz. However, convection A continued to weaken.

Discussion
In the current research, the triggering and development of the isolated convection around the TUR had a close relationship with the interaction of the sea breeze front and gust front in the context of the UHHI effect. Previous research results show that the convections were more likely to form at the crossing point of the collision between the gust front and sea breeze front, and the convections occurred after the collision of the sea breeze front and gust front [6,39]. The literature [6] further revealed the physical characteristics of convections at the collision intersection through numerical simulation, that is, there are obvious vertical motion and convergence characteristics, and the distribution characteristics of the generalized Richardson number are also significant, whose thickness is about 1.0 KM, and the CAPE value increased significantly before the initiation of convections. However, another phenomenon was identified in our research, and we present some interesting findings. The isolated convection was first triggered by the sea breeze front, and then the gust front advanced to the convection area and interacted with the sea breeze front, resulting in the merging of multiple convective cores in the isolated convection and the continuous strengthening of the isolated convection.
The UHI effect refers to the temperature in the urban region being significantly higher than that in the rural region. The investigation of a UHI event showed that the interactions between the UHI effect and the sea breeze front indeed existed [44,55]. Previous research shows that the high-temperature area of the UHI effect could provide favorable thermodynamic conditions for the isolated convections [13]. However, another phenomenon was identified in our research, which is also worthy of our attention, that is, that the urban region is not located in the center of high temperature, but is located around the center of high pseudo-equivalent potential temperature, indicating that the urban region has the characteristics of high temperature and high humidity compared with the surrounding rural area. As described in Section 4 of this paper, this meteorological condition is also conducive to forming unstable conditions and providing favorable thermal conditions for the initiation and development of strong convective weather, which played the same role as the UHI effect. However, this phenomenon does not meet the definition of the UHI effect. In order to better explain the physical mechanism of the isolated convection in urban regions, this phenomenon was defined as the phenomenon of urban humid heat island (UHHI) in this paper.
As is well known, it is very challenging for the most advanced numerical weather (c) (d)   Figure 10a), convection B was formed. The echo center was at the height of 3-6 km. At 19:12 BT (see Figure 10b), convection B continued to strengthen and the echo center was at the height of 1 km (the lowest height at which the echo can be detected by the Tanggu Radar at this position). After that (see Figure 10c,d), convection B continued to strengthen, the echo center was still at the height of 1 km, and reflectivity of the echo center increased to more than 50 dBz. However, convection A continued to weaken.

Discussion
In the current research, the triggering and development of the isolated convection around the TUR had a close relationship with the interaction of the sea breeze front and gust front in the context of the UHHI effect. Previous research results show that the convections were more likely to form at the crossing point of the collision between the gust front and sea breeze front, and the convections occurred after the collision of the sea breeze front and gust front [6,39]. The literature [6] further revealed the physical characteristics of convections at the collision intersection through numerical simulation, that is, there are obvious vertical motion and convergence characteristics, and the distribution characteristics of the generalized Richardson number are also significant, whose thickness is about 1.0 KM, and the CAPE value increased significantly before the initiation of convections. However, another phenomenon was identified in our research, and we present some interesting findings. The isolated convection was first triggered by the sea breeze front, and then the gust front advanced to the convection area and interacted with the sea breeze front, resulting in the merging of multiple convective cores in the isolated convection and the continuous strengthening of the isolated convection.
The UHI effect refers to the temperature in the urban region being significantly higher than that in the rural region. The investigation of a UHI event showed that the interactions between the UHI effect and the sea breeze front indeed existed [44,55]. Previous research shows that the high-temperature area of the UHI effect could provide favorable thermodynamic conditions for the isolated convections [13]. However, another phenomenon was identified in our research, which is also worthy of our attention, that is, that the urban region is not located in the center of high temperature, but is located around the center of high pseudo-equivalent potential temperature, indicating that the urban region has the characteristics of high temperature and high humidity compared with the surrounding rural area. As described in Section 4 of this paper, this meteorological condition is also conducive to forming unstable conditions and providing favorable thermal conditions for the initiation and development of strong convective weather, which played the same role as the UHI effect. However, this phenomenon does not meet the definition of the UHI effect. In order to better explain the physical mechanism of the isolated convection in urban regions, this phenomenon was defined as the phenomenon of urban humid heat island (UHHI) in this paper.
As is well known, it is very challenging for the most advanced numerical weather prediction (NWP) models to forecast the accurate position and time of isolated convection around urban regions. Real-time surface observations, radar measurements and other observations can help to reveal the mechanism of isolated convection. In our forthcoming study, more cases of isolated convection will be collected and analyzed to further examine the roles of the UHI (or UHHI) effect, mesoscale boundary and convective interaction in determining the CI timing and location of the isolated convection around the TUR.

Conclusions
The isolated convection was observed in the TUR on 1 August 2019. The present study analyzed the ERA5 reanalysis data, the surface meteorological elements data observed by surface observations, and the reflectivity detected by the Tanggu S-band Radar, and discussed the mechanisms that led to the initiation, maintenance, and development of the isolated convection. The main conclusions are summarized as follows: (1) This isolated convection event occurred about 250 km away from a low-level convergence line to the north and in the warm and humid area dominated by southerly flow at a low level. The vertical shear of wind velocity at a low level was not obvious (mostly less than 4 m·s −1 between 1000 hPa and 850 hPa), but the atmospheric thermodynamic conditions (i.e., the CAPE of more than 3500 J/kg and CIN of nearly zero) were beneficial for the triggering and development of the isolated convection.
(2) Moreover, 42 min before the initiation of the isolated convection, the sea-breeze front started to influence the TUR and continued to move west. At the same time, the TUR was located at the intersection of the high temperature band and high dew point temperature band, with characteristics of high temperature and high humidity. The surface pseudo-equivalent potential temperature around the TUR was more than 383 K, with 2-3 K contrast of the surrounding area, forming the UHHI effect. That is to say, the sea breeze front triggered the isolated convection around the TUR in the context of the UHHI effect.
(3) The gust front was formed between the cold pool and the warm and humid air of UHHI at 17:00 BT on 1 August 2019. The formation of cold pool was caused by the previous convections, which was 20-30 km away from the TUR. After the initiation of the isolated convection (18:00 BT), the gust front moved northward, started to influence the convection and collided with the sea breeze front, which made the many convections merge into an isolated convection. In the following 30 min (from 18:00 to 18:30 BT), the convection developed continuously and grew into its mature stage with the interaction of the sea breeze front and gust front. Maybe the gust front was not the direct reason of the isolated convection's initiation, but rather the gust front moving northward, approaching the convection and interacting with the sea breeze front, resulted in convections merging into an isolated convection and the convection's development.
(4) Due to the heavy precipitation produced by the isolated convection, a new cold pool was formed around the TUR. The westerly cold outflow subsequently met with easterly winds, leading to a new convection being triggered in the east of the TUR. Subsequently, the old convection continued to weaken, and the new convection continued to strengthen, but it also weakened and disappeared at 20:36 BT.