On Thunderstorm Microphysics under Urban Heat Island, Sea Breeze, and Cold Front Effects in the Metropolitan Area of São Paulo, Brazil

This work discusses issues related to the impact of urbanization on the microphysical processes of precipitating systems associated with synoptic, mesoscale, and local scale systems. Among the issues addressed is the impact of urban heat island (UHI) in São Paulo city center and urban densification (UD) in the Metropolitan Area of São Paulo (MASP) on the microphysical, dynamic, and thermodynamic properties and distribution of precipitation and heavier rainfall from sea-breeze (SB) e cold-front (CF) combined during their space-time evolutions. For this purpose, it used four components: classification of hydrometeors with fuzzy logic, calculation of the raindrop diameters, an estimate of liquid water mass and ice mass from polarimetric-variables measured with dual-polarization X-band meteorological radar. The results indicated that urban densification (UD) and heat island (UHI) of the São Paulo city center impact the formation of precipitation, liquid water mass, and ice mass, depth, and duration of a thunderstorm. It was also observed the asymmetric configuration of the thunderstorm is induced by the strong convergence in the São Paulo city center, and this strong convergence is induced by the intense heat island (UHI) in the São Paulo city center. Was also observed that this event that is formed in the Metropolitan Area of São Paulo (MASP) depends on microphysical processes of mixed-phase of the cloud (water and ice) above the 0 ̊C isotherm for the production of intense rain and cold pool at the surface. These important microphysical processes within long-lasting secondary convective cells over the São Paulo city center should be taken into account on convective parameterization schemes as well as the associated cold pool dynamics.


Introduction
The Metropolitan Area of São Paulo (MASP) is among the five largest megacities of the planet. The diurnal cycle of surface diabatic heating, induced local circulation in an aerosol-rich environment is the main ingredients of deep thunderstorms in MASP [1]. The MASP is recurrently affected by flash floods, mudslides, heavy winds, lightning, and hail that cause significant social and economic impacts [2].
The geometric center of MASP is different from suburban surroundings in both material and geometry, with multiple intersecting paved roads, impervious surfaces and denser irregular constructions in urban areas causing the heat island (UHI) of the São Paulo city [1]. About 70% of all severe weather events within MASP were associated with local circulation induced by topographic and UHI effects [3]. The SB is the most common local circulation. In midafternoon hours while the SB front pushes continent ward from southeast against prevailing northwest winds (approaching CF), injecting a deeper layer of moisture over the MASP warmer and drier urban boundary layer, often causing major flood episodes [4]. The UHI is caused by large sensible heat storage within the urban canopy in the early hours of the day and released during the afternoon and night over the center of the MASP [1].
The SB fronts tend to increase environmental shear and associated vorticity and moist [5] and MASP urban aerosol convergence [2]. When the SB front and gust fronts of the primary cells collide there is a reduction of moist static stability and new convective cells initiate over the center of the MASP [6]. According to [7] the UHI forms a strong convergence zone in the center of the MASP and thereby accelerates the SB front toward the center of the MASP. After it arrives in the urban region, the SB front stalls over the MASP for 2 h. Subsequently, the sea breeze progresses beyond the city when the heat island dissipates. Thereafter, the SB front propagates beyond the urban area at a decelerated rate. The effects of the UHI on convective cell initiation and thunderstorm, resulting in the temperature contrast as a condition for storm initiation over the MASP [8]. The UHI effects amplify the thunderstorm activity by local factors (sea and country breeze, synoptic events, and terrain) and the higher flash rates throughout the urban region, which influences the lightning density in the MASP [8].
Under these local circulations and aerosol effects, microphysical processes tend to play an important role in the dynamics and thermodynamics of thunderstorms in MASP [9] [10] [11]. [12] analyzed measurements of drop size distribution in two thunderstorms in MASP. According to [12], the effects of the urban aerosol load were the factor affecting the deep convection of two thunderstorms (lower aerosol load and higher aerosol load) of the SB events. The thunderstorms associate with the higher aerosol load (on January 11, 2010) longer-lasting convection and larger raindrop concentration than the one observed in the Amazon Forest during LBA [11] [12]. [11] reported about the thunderstorm (on January 11, 2010) caused by synoptic conditions and typical deep SB circulation over MASP. The CF moved through over South Brazil which intensified northwest circulation and convergence against the SB front midafternoon. This scenario triggered thunderstorms under a moist unstable environment and strong convergence over MASP.
In the most intense urban heat island conditions on the polluted MASP occurs deep storms that depend on the humidity provided by the sea breeze and the prefrontal flows over the MASP, which caused flooding and other severe climatic characteristics on January 11, 2010. The amount of rain recorded was 177.3 mm, 68.7% of the historical average for January which is 258 mm [11]. The state of São Paulo was with the twenty-one (21) cities in an emergency and two (2) in a public calamity situation, caused by a sequence of 11 days of heavy rains in January 2010. In MASP it rained almost 70% of the total precipitation of January between days 1 and 12. This SB event was associated with the interaction of MASP heat island and CF and caused a high accumulation of rain. However, observational studies, modeling, surface data, and remote sensing estimates by meteorological radars and satellites were conducted to better understand the dynamics of SB and CF interaction with local circulations induced by the presence of MASP. The processes by which convection starts and is organized were understood to adequately represent the daytime precipitation cycle in numerical models [5] [7] [11] [13]. But the evolution of microphysical processes must be studied to provide a physical basis for estimating precipitation by remote sensing [14].
Hydrometeor classification using polarimetric radar data provides detailed information to validate and improve microphysical and convective parameterizations.
Microphysical information is known to significantly affect local and mesoscale numerical models [15] [16]. Therefore, microphysical effects are of extreme relevance for understanding the behavior of convective episodes [17].
Microphysical knowledge of thunderstorms is vital because, in recent decades, damage caused by these events seems to have increased because of urban microclimatic changes due to anthropogenic actions. The impact of severe convention episodes can be minimized by accurately predicting, saving lives, and preventing economic losses [18].
The MXPOL weather radar measurements detected with high spatial-temporal resolution the penetration of the sea breeze at MASP where the heat island was already intense. The satellite image sequences showed the arrival of the cold front over MASP at the same time with penetration of the sea breeze which probably caused hail precipitation and microbursts at one end of the storm. Thus, this work was motivated to analyze the microphysics of a rare event induced by the combination of systems of different scales on January 11, 2010, that, eventually, the microphysical inference can provide information for parameterizations of cloud and mesoscale models.

Study Area, Dataset, Instrumentation, and Microphysics Analyses
The study was performed over the Metropolitan Area of São Paulo-MASP ( Figure  1). The MASP is a megacity environment composed of 39 cities [19]. It has been increasing in vertical growth with very limited vegetated areas in the São Paulo city. The MXPOL was strategically placed to the west of São Paulo city to monitor the propagation of sea breeze fronts. The Joss-Waldvogel disdrometer (JWD) was placed approximately 30 km away from the MXPOL (Figure 1(a)). Figure  1(b) shows a visible satellite image of MASP and its surroundings. It is quite remarkable the color contrast between the MASP surface (brownish) and its ocean (blueish) and rural (greenish) surroundings (Figure 1(b)). The land cover change (brownish) caused by anthropogenic action through the last century has resulted in drastic changes in its energy and water budgets in MASP ( Figure 1(b)). About 70% of all severe weather events within MASP were associated with local circulation induced by topographic and UHI effects [2].

Technical Features of the MXPOL
MXPOL is a mobile X-band (9.37 GHz), dual-polarization Doppler weather radar system switchable to single polarization mode with a maximum of 601 gates at 150 km range. In horizontal-only transmitted mode with a longer pulse length, the sensitivity of the radar is better than 3 dB when compared to a dual-polarization mode of operation. In the horizontal-only transmitted mode, the dualchannel receiver is capable of measuring both co-and cross-polarized backscattered signals. Table 1 shows the technical characteristics of the MXPOL antenna, pedestal, transmitter, and receiver reflector. The volumetric scanning strategy for 11 January 2010 was a clockwise azimuth scan low to high elevation angle at a 5-minute time interval (0.6˚, 1.2˚, 2˚, 3˚, 4.5˚, 6˚, 9˚, 12˚, 16˚, 20˚, and 25˚). The collection of successive PPIs (Plan Position Indicator) of the radar is a  volumetric scanning. Radar is one of the best tools available for monitoring and forecasting as it provides real-time monitoring with height spatial resolution and allows a threedimensional view of the data [20]. The radar consists of a transmitter, a receiver, an antenna, and a decoding and processing system ( Table 1). Radar is a remote sensing tool, that is, through electromagnetic radiation, collects data without coming into direct contact with the study target.
Data collection occurs when the radar transmits a beam of electromagnetic radiation and captures the energy reflected by the particles present in the atmosphere. This energy is amplified and, given the time interval between emission and echo return, it is possible to determine the distance from the target to the radar and to describe the hydrometeors types.
The following variables, obtained from MXPOL volume data were used in this study: Reflectivity (Z). The Z is the measure of a target's efficiency in intercepting and returning the energy emitted by the radar. Thus, reflectivity is dependent on the sizes, formats, quantity, and dielectric properties of the targets [21]. High reflectivity values (up to 45 dBZ) are generally associated with heavy rainfall, if these values are between 5 and 12 km of altitude, they are associated with severe storms [22].
Radial velocity (V). The V is how fast the detected target is moving away or approaching the radar. Using the Doppler effect, the radar estimates the velocity at which the target moves toward or moves away from radar [21]. By convention, blue tonality (negative values) represents particles coming toward the radar, and red tonality (positive values), particles moving away from radar [22]. Differential Reflectivity (Z DR ). The Z DR is the difference between Z HH (reflec-tivity with the signal emitted and captured horizontally) and Z VV (reflectivity with the signal emitted and captured vertically) [21]. The variable Z DR is widely used to differentiate rain, hail, snow, and even non-meteorological targets [23].
Phase Differential (PHIDP). The PHIDP is the difference between the phase emitted and received horizontally and the phase emitted and received vertically in a pulse [22]. The variable PHIDP has radial variations, since cumulative changes in the phase difference for the complete pulse journey.
Specific Phase Differential (K DP ). The K DP variable is the spatial derivative of PHIDP and is an excellent estimator for precipitation, as presented by [24], values above 2˚/km indicate a significant amount of liquid or oblate water. Copolar Correlation (RHOHV). The RHOH represents the correlation between the horizontal and vertical polarized Z signals at a given point in space [20].
The values of RHOHV vary between 0 and 1.0 and can be seen as a measure of the similarity between horizontal and vertical polarization signals [23]. This variable describes the physical characteristics of the target, being closer to 1.0 the more uniform the target since the vertical and horizontal signals tend to be the same.
The data collected by the MXPOL are in spherical coordinates, where the position of each data is described in terms of (θ, φ, r), where θ is the angle of antenna elevation, φ is the angle of clockwise rotation of the antenna relative to the geographic North, called the azimuth angle, and r is the distance from the radar to target.

The MXPOL and JW Disdrometer Datasets
The 11 January 2010 thunderstorm microphysics was analyzed through Dual Doppler polarization (DDP) band X radar and JW disdrometer data. Five hours of polarimetric data namely, the reflectivity of horizontally polarized waves Z H ; differential reflectivity Z DR ; specific differential phase K DP ; and correlation coefficient between horizontally and vertically polarized signals ρ HV were used in the microphysical analyses through the supervised fuzzy logic method during the spatial-temporal evolution of the storm. Three and half hours of JW Disdrometer (JWD) data were used for droplet spectrum analysis during the spatial-temporal evolution of the storm.
Hydrometeor types were mapped on PPIs to analyze the thunderstorm mixedphase microphysics. Datasets were selected under MASP heat island, polluted conditions, SB, and CF effects. The above MXPOL polarimetric variables and respective vertical temperature profiles were input to the hydrometeors classification (HC) method and then visually verified against typical hydrometeor polarimetric signatures.

Hydrometeor Classification (HC)
The fuzzy logic method is a technique for synthesizing information from polarimetric variables into a hydrometeor classification algorithm [25]. We used the fuzzy logic hydrometeor identification algorithm for XPOL radar data developed by [26]. Four polarimetric variables (Z H , Z DR , K DP , and ρ HV ) from MXPOL radar and local altitude (temperature) were used as input data for hydrometeor classification of precipitating systems using the fuzzy logic method. Temperature (˚C) and height (m) data from the radiosonde of Campo de Marte aerodrome were interpolated for each grid point at each PPI. The interrelations of polarimetric variables and altitude generated information for automatic identification of the hydrometeors represented at each PPI grid point.
• Python ARM Radar Toolkit (Py-ART) The Python ARM Radar Toolkit (Py-ART), an open source library for working with weather radar data [27]. Py-ART was used to read and visualize the polarimetric MXPOL weather radar data and perform the hydrometeor classification using the fuzzy logic method, calculate the raindrop diameter (D 0 ), the solid water mass and the liquid water mass at each grid point in all the IPPs and generate graphs with the algorithms available in the CSU Radar Tools platform developed by [25]. The algorithms are implemented in modules such as Scientific Python, including NumPy, SciPy and matplotlib and Python to realize the interface with weather radar polarimetric variables with radiosonde data, etc, and computationally accelerate their execution.
• The fuzzy logic method The HC in the thunderstorm was performed using the Fuzzy logic method using python programming language [27]. This method uses smoother transitions within the limits of the observed polarimetric variables, for each type of hydrometeor and simple rules rather than analytical equations to describe the system in question [25]. This method is used to decide on the right output class-based on 4 input polarimetric variables and temperature ( Figure 2). HC was carried out discriminates several types of hydrometeors for X-pol based on weather radar measurements, T-matrix simulations and Mueller-matrix scattering models aiming at defining membership beta functions (MBFs) implemented by [28]. Figure 2 shows the implemented structure of the fuzzy logic algorithm consisting of four steps [29] and the typical values used for dual-pol variables are shown in Table 2. The algorithm uses the associated predefined functions to automatically classify the multidimensional input data. The classification was applied at each gate of the weather radars such that a complete field of hydrometeor type was produced consistent mainly with the polarimetric data of the SPOL radar. The hydrometeor types were mapped as a PPI for microphysical analysis during the spatial-temporal evolution. Table 2 contains the hydrometeors types corresponding to the thresholds of the Z H , Z DR , K DP , ρ HV , and temperature.
The method is robust enough that its performance cannot be affected due to typical measurement error in some of the input variables [30] [31]. This procedure is reliable for automatic classification and is most suitable for hydrometeor classification [28].
The fuzzy-logic method consists of four major components ( Figure 2): fuzzification, Table 2. Hydrometeor types corresponding to the thresholds of the ZH (dBZ) reflectivity, ZDR (dB) differential reflectivity, KDP (km −1 ) specific differential reflectivity, ρHV correlation coefficient, and altitude data. Source: [29].  Figure 2. Illustration of a fuzzy logic hydrometeor classification algorithm developed by [28]. Adapted from [29]. inference rule, aggregation, and defuzzification [28]. Fuzzification: fuzzy logic schemes transform discrete data into Fuzzy sets with their corresponding weight. A specific input variable may belong to different Fuzzy sets, but also with different weights. The most important component of this step is the "member function" [28], which is used to describe the relationship between input data and the fuzzy set. During the fuzzification of each variable, all hydrometeors receive a value between 0 and 1 through the member function (MBFs), which is the degree to which that hydrometeor is related to the fuzzy set it belongs.
Inference: in a Fuzzy logic system, rules are used to describe the relationship between input and output fuzzy variables in Boolean form, that is, if the input data is x, then the output data must be y [28]. The degree of importance's rules within the Boolean form can be determined by a series of methods, such as minimum correlation, product correlation, and minimum-maximum relationship [28]. The result of fuzzification is multiplied by a predetermined weight, a value that is determined by its usefulness in classifying and measuring the accuracy of input variables.
Aggregation: inference methods can be used to calculate the weight of each rule, so the aggregation method is used to determine the entire Fuzzy space to be analyzed. The values obtained by the inference rule are then aggregated to form a single value for each Fuzzy set, that is, hydrometeor or hydrometeor set.
Defuzzification: the output of the aggregation process is a fuzzy value. In many applications, such as hydrometeor classification, it is necessary to find a discrete output value that best represents the fuzzy output value, so the data is transformed into dimensional space. The average of the maximums that "defuzzify" the data and determine the strength of each fuzzy set, defined by this value the class of hydrometeors that it belongs to. In addition, it generates an output field, where the hydrometeors are classified into 10 categories (Figure 2).

• Membership functions (MBF)
The MBF determination is an important step of the fuzzy-logic-based HC method. [31] Developed the parameters and the membership functions for X band polarimetric variables for each hydrometeor that was used in this work, by the following equation: where x = polarimetric variables; m = center; a = width; b = slope. Table 3 contains the parameters set of the membership functions for four polarimetric variables and temperature for each hydrometeor type used in Equation (1).
• Estimates of mass of liquid water, the mass of solid water, and mean droplet diameter. The methodology of [32] was applied for the estimates of liquid water mass and solid water mass [26]. The mass of ice and water was calculated at specified height levels, given Z H and Z DR at each grid point [28] from the following Equations: where 10 10 DR Z H = horizontal reflectivity of the hydrometeor classified as rain (mm 6 •m −3 ); Z V = vertical reflectivity of the hydrometeor classified as rain (mm 6 •m −3 ); Z DR = differential reflectivity of hydrometeor classified as rain (dB).
To calculate the mean raindrop diameter, the following relationship obtained by [32]:

Synoptic Analysis
The IR satellite images (Figure 3) show the displacement of a CF towards the southeastern region. Where it is observed the most intense cloud band extremity over the Santa Catarina and Paraná states at 1600 UTC (Figure 3(a)) and the eastern part of the São Paulo state is covered with warm cloud of prefrontal system. However, Figure 3(b) shows an intense cloud-band moving towards MASP (1730 UTC, highlighted area), that might be associated to the sea breeze front. Figure 3(c) shows that part of the eastern region of the São Paulo state is cloud-covered, but except MASP (2000 UTC), already with the SB front over the center of MASP. The cloud-band is asymmetrically configured because of the blocking of the SB-front by intense MASP heat island and northwest winds.
This SB-front is temporarily stagnant over the center of MASP caused by the intense MASP UHI. At 2100 UTC, while the SB-front moves outside MASP push continent ward from southeast against prevailing northwest winds, injecting a deeper layer of moisture over the MASP warmer and drier urban boundary layer, causing a deep cloud trigger (highlighted area).

Thermodynamic and Dynamic Analysis
Figure 4(a) shows the temperature, humidity, and wind profile measured with the radiosounding at Campo de Marte Airport (23˚30'31"S, 46˚38'15"W) at 1200 UTC on January 11, 2010. It shows a very high mixing ratio between 700 -800 hPa in almost all troposphere. It is observed a strong entanglement of dry and hot air mass between 450 hPa and 700 hPa of medium atmosphere, which favors the deep mixing of the boundary layer, causing strong updrafts that carry moisture to the tropopause. These thermodynamic conditions strengthen convective cells [33].
The CAPE indicated 1491 J•kg −1 and Lifted index of −4˚C (Figure 4(a)). The superficial heating during the day and the injection of moisture in the afternoon ( Figure  4(a)), increased the instability indexes in an environment of weak vertical shear. Figure 4(b) shows the vertical moisture profile cross-section at longitude 47˚ 60'W from latitudes 22˚25'S to 24˚95'S of MASP at 1800 UTC on January 11, 2010. A moisture layer is observed varying between 60 -90 to tropopause over the MASP during the SB-front passage. A strong dry air mass is observed between 400 and 300 hPa. The surface heating during the day and the ejection of humidity by the SB-front and the northwest wind canalized by the frontal system during the afternoon (Figure 4(b)). This situation favored the increase of instability indices in an environment of weak vertical shear at 1800 UTC. Figure 4(c) shows a mesoscale signature in the radial velocity field at 25 km away from the radar. A convergence region induced by SB-front (blue tone), and the northwest wind (red tone), is observed over the center of the MASP at 2025 UTC, highlighted area. The radial velocity measured of the flows is ~8 m•s −1 . The air mass flow associated with northwest wind associated with the jet (red tone), propagated against the SB-front (blue tone), triggering the convective cells over the center of São Paulo city. A very fine line of reflectivity between 0 dBZ and 10 dBZ at the 40-km range almost perpendicular to 150˚ azimuth (inside the rectangle) was moving towards MXPOL, indicated by a black arrow (Figure 5(a)). It is the SB-front and its effect on convection over MASP. The SB-front intensifies as it spreads to the center of MASP. It shows PPI 0.6˚ elevation of unfiltered reflectivity fields at 1753 UTC, 1850 UTC, 1940 UTC, 1955 UTC, 2015 UTC, and 2025 UTC, when the SB-front reaches the geometrical center of the RMSP, it presents an asymmetric configuration induced by the intense UHI in the region (Figure 5(d)). The low elevations PPIs of Figure 4(b) is passing through the large sensible heat storage within the urban canopy, where the first cell convection trigger at the eastern extreme of São Paulo city (1850UTC), which later moves out of São Paulo city (1940 UTC). At 1955UTC, the incoming SB-front interacts with the MASP UHI over the center of the São Paulo city and it stagnates temporarily over the center of São Paulo center, triggering convective cells in its rear ( Figure 5(d)). Figure 5(e) shows the convective cells with larger dimensions located towards northeast and east MXPOL at 2015UTC. These cells produced gust fronts, and these moved along the SB-front and triggered secondary convective cell (highlighted in the circle) intensifying and interconnecting the old cells at 2025UTC (Figure 5(f)). It is the beginning of the production of very deep convection with heavy rainfall, gusty winds, and lightning over the MASP geometrical center. Figure 5(h) shows the unique cell thunderstorm, resulting from the aggregation of several cells. It observed the thunderstorm with asymmetric configuration induced by the intense UHI in the MASP geometrical center and also shows that the south end of the thunderstorm interacts directly with the SB-front. The fine line of reflectivity between 20 dBZ and 25 dBZ is the gusts of wind. These gusts of wind resulted from microburst which causing the cold pool to spread in all directions on the surface at 2050UTC and 2105UTC ( Figure 5(h), Figure  5(i)).

Temporal Thunderstorm Evolution
The SB-front was already outside of the study area ( Figure 5(i)). The convective cell continues to produce gust front witch moving towards MXPOL, but without triggering new cells due to the absence of the SB-front on-site at 2105 UTC and 2115 UTC ( Figure 5(j)). It is also observed that the thunderstorm loses the asymmetrical configuration. Therefore, the thunderstorm began to weaken, with a predominance of a stratiform system from 2135 UTC (Figure note show) to 2155 UTC ( Figure 5(k)).
In the highlighted area of Figure 6(f), we have the Zdr values, where green tones (Zdr ~ 3.2 dB) are noted, which indicates the presence of large drops  between 1 mm and 5 mm in the two secondary cells incorporated. In the highlighted area, Figure 5(f), the reflectivity field presents values > 55 dBZ. Figure 7 shows the correlation coefficient that indicates all the cloud boundary areas felling the convective cells, and the storm cloud resulting from the incorporation of the  secondary convective cells. Figure 8 shows the specific differential reflectivity, where the most positives values indicate more intense rainfall areas with a predominance of oblate raindrops.

Hydrometeor Classification (HC), Liquid Water Mass, Liquid Ice mass, and Drop Spectrum
• High hail concentration in the secondary cell Figure 9 show the high elevation PPIs (above the 0˚C isotherm) of hydrometeor I. M. T. J. Ihadua, A. J. P. Filho classification, reflectivity, differential reflectivity, specific differential reflectivity, liquid water mass, liquid ice mass, and drop diameter through the intense convective cell resulting from the feedback between the gust fronts of the first convective cells with SB-front, which caused the secondary cell trigger which caused greater hail production than other thunderstorm over the São Paulo city center. The convective cell exhibited reflectivity values above 55 dBZ at freezing temperatures (Figure 9(b), in the highlighted area). In the highlight of Figure 9(c), Figure 9(d), and Figure 9(g), we have the differential reflectivity, specific differential reflectivity, and drop diameter values, where red tones (Z DR between −0.5 and 2 dB, and K DP between 0.5˚ and 6˚ km −1 ) and drop diameters between 1 and 2 mm are noted, which indicates the presence of large drops and hail.
The highlight in Figure 9(e) indicates liquid water mass in the convective cell and ice mass with very high values at this altitude for liquid water (between 2.5 Figure 9. PPIs at 9.0˚ of (a) hydrometeor classification (HC), (b) reflectivity (ZH), (c) differential reflectivity (ZDR), (d) specific differential reflectivity (KDP), (e) liquid water mass (f) liquid ice mass and (g) drop diameter at 2032 UTC on 11 January 2010. In the highlighted area, is the intense convective cell resulting from the integration of two secondary convective cells. MXPOL centered circumferences at ∆r = 25 km. Scale, Geographic contours, and UTC times are indicated. Hydrometeor types, reflectivity (dBZ), differential reflectivity (dB), specific differential reflectivity (km −1 ), water, and ice masses (g•m −3 ) in color scales. . Figure 9(a) indicates the hydrometeor classification and in the highlighted area is the region with the highest hail concentration of the precipitating system which is in agreement with the values of the polarimetric variables and the reflectivity.
The dynamics of strong convergence over the center of the São Paulo City was the cause of these peculiar microphysical characteristics. When hot-humid air associated with urban pollution over the center of São Paulo City is loaded to the cloud at negative temperatures, water vapor condenses producing a high amount of droplets and latent heat performed in the cloud. This increased instability and intensified the convergence of air masses of different properties (Figure 4(c)), which carried a high amount of water drops at negative temperatures.
• Strong downdrafts Figure 10 shows the lower elevation (0.6˚) PPIs of hydrometeor classification, reflectivity, differential reflectivity, specific differential reflectivity, correlation coefficient, liquid water mass, liquid ice mass, and drop diameter of the storm's mature stage resulting from the integration of secondary convective cells over MASP.
Note that the storm presents an asymmetric configuration similar to the SB front caused by the intense convergence at low levels induced by the strong UHI in the São Paulo city center. The rectangular highlight in Figure 10  The hydrometeors classified in the PPI 0.6˚ elevation in Figure 10(a), are the same as those observed in the PPIs 3.0˚ elevation in Figure 11(a), characteristic of the collapse of the liquid water mass and solid water mass from the upper part from the cloud to the surface as observed highlighted areas. This aspect is characteristic of a microburst (Figure 10(a)), highlighted area. It is also observed reflectivity between 20 and 25 dBZ within concentric associated with the cold Figure 10. PPIs at 0.6˚ of (a) hydrometeor classification (HC), (b) reflectivity (ZH), (c) differential reflectivity (ZDR), (d) specific differential reflectivity (KDP), (e) correlation coefficient (ρHV), (f) liquid water mass (g) liquid ice mass and (h) drop diameter at 2055 UTC on 11 January 2010. In the highlighted area, is observed mixture of rain, large drops, and hail is observed due to the probable occurrence of a microburst. MXPOL centered circumferences at ∆r = 25 km. Scale, Geographic contours, and UTC times are indicated. Hydrometeor types, reflectivity (dBZ), differential reflectivity (dB), specific differential reflectivity (km −1 ), correlation coefficient (without unity), water and ice masses (g•m −3 ) in color scales. pool expelled by the gust front resulting from the strong downdrafts ( Figure   10(b)).
The diameter of the raindrops (Figure 10(h) and Figure 11(h)) shows the concentration of drops of almost the same size at two different levels at the moment of collapse of liquid water mass and ice mass from cold thunderstorm region to surface at 2055 UTC. These results corroborate with DSD measured on the surface (Figure 12). Probably the DSD modification was caused by the collision-coalescence of big drops on cloud droplets, breakup by the collision between big drops, and breakup by aerodynamic force combined (Figure 10(h) and Figure 11(h)). This mechanism acted together to increase a large number of raindrops diameter from 0.551 to 1.116 mm (highlights 1 of Figure 12 Figure 12(b)) and from 1.7 to 2.5 mm (highlights 2 of Figure 12(a) and Figure   12(b)). Figure 12 shows the Time evolution of the numbers and size of raindrops of the storm in the mature and decay phases and Averaged drop size distributions estimated with JWD measurements of 76 rainfall storm from 1830 to 2100 on 11 JAN 2010. The sea breeze had just was over the JWD at 1753 UTC (Figure 4(a)). A first convective cell developed right above de JWD location at 1955 UTC ( Figure 5(d)). JWD measurements of the main Storm with asymmetric and completely integrated configuration resulting from secondary convective cells started at 2050 UTC ( Figure 5(h)). By 2155 UTC, the dominant rainfall was of the stratiform type ( Figure 5(k)).
The time evolution of the drop spectra indicates the convective (in highlights 1 and 2), transition, and stratiform (in highlight 3) phases of the storm presenting raindrop diameters in all classes. It is interesting to note the high number of droplets estimated and the concentration of raindrops in highlights 1 and 2 of Figure 12(a) and Figure 12(b), and decrease in the number of larger diameter  [19]. raindrops and concentration of raindrops. It suggested that the injection of a deeper layer of moisture, as well as boundary layer shear, induces a higher drop concentration at all diameters. Furthermore, urban densification has impacted the increase of raindrop concentration.

Discussion
The fuzzy logic applied to MXPOL radar data as reflectivity, differential reflectivity, specific differential reflectivity, and the correlation coefficient was successfully performed to obtain the hydrometeors classification present in the thunderstorm. The attenuation was not corrected to the spatiotemporal evolution of the sea breeze front over MASP. However, this thunderstorm was intense and caused signal attenuation, and probably influenced the hydrometeors classification as graupel and hail. For the X-band radar still a case study under development for the hydrometeor classification consistency from an intense storm. Another aspect that may also have influenced hydrometeors classification is the altitude and temperature data. The sounding data from Campo de Marte Airport of São Paulo city at 1200 UTC were used the nearest to the event's occurrence. This event began at 1700 UTC, five (5) hours before the start of the event. Was also calculated the liquid water mass and ice mass, and raindrop diameter using the reflectivity and differential reflectivity.
The reflectivity fields showed primary and secondary cell triggers at the rear of the SB front as it propagated towards the geometric center of the MASP. These convective cells were aggregated when the SB front stagnated over the geometric center of MASP, forming a single convective cell (thunderstorm) with a north-south orientation over the São Paulo city.
The SB fronts tend to increase environmental shear and associated vorticity and moist and MASP urban aerosol convergence [2]. When the SB front and gust fronts of the primary cells collide there is a reduction of moist static stability and new convective cells initiate over MASP [6].
The UHI caused by urban densification at the geometric center of MASP induced strong convergence of southeast and northwest flow over São Paulo city center which blocked the advance of the SB front for almost2hours.TheMASP UHI is caused by large sensible heat storage within the urban canopy in the early hours of the day and released during the afternoon and night over the center of São Paulo [5]. This effect induced an asymmetric configuration of the SB front and the thunderstorm. Part of the southern end of the storm interacted directly with the stagnant SB front. Subsequently, the sea breeze progresses beyond the São Paulo city when the heat island probably dissipated, and the thunderstorm loses its asymmetrical configuration and sometime later begins to dissipate.
Environments with high aerosol concentration may favor the development of strong updrafts reaching high altitudes to the base of the tropopause, allowing the cloud to develop to this level [34]. The radial velocity fields show strong convergence at low levels over the São Paulo city center. Several authors [5] [6] [11] have claimed that the SB front moves towards RMSP, where surface rugosity and UHI increase humidity convergence and upward movement, respectively.
A strong downdraft probably was also occurred at the southern end of the thunderstorm where the microphysical characteristics of the hydrometeors (big drops and hail) found after the base of the 0 o C isotherm were the same found in the lower part of the thunderstorm. As observed by [35], observations of airplane sonice crystals in the cloud, suggest that the mechanism of production in the secondary droplet freezing process in ice production is in the strong updrafts that carry these droplets and operate for the formation of freezing nuclei. The consequence of this is the production of small ice particles colliding with supercooled droplets, which produce the hailstones within the strong updrafts [35], as observed in the results.
According to [20], as the ice mass precipitates below the 0˚C isotherm, part of the solid water mass thaws, providing the formation of big drops. According to some authors, for the thawing process to occur, energy is needed from the surrounding air with more hot temperatures [20]. The thawing process cools the air and causes cold air to accelerate downdrafts.
The results showed a high concentration of liquid water mass and ice mass of 12.5 gm −3 above the 0˚C isotherm in the area where the peculiar concentration of hail occurred. A high concentration of liquid water mass of ~17.5 gm −3 was also observed after the base of the 0˚C isotherm and in the lower part of the thunderstorm where the strong downdrafts occurred as well as the probable occurrence of microburst. But the concentration of ice mass varied, it was observed high concentration of the 5 gm −3 and a low concentration of 2.5 gm −3 after the base of 0˚C isotherm and a high concentration of 4 gm −3 in the lower part of the thunderstorm. According to [36] and [37], latent cooling performed or latent heat absorbed during freezing, evaporation, and sublimation negatively increases buoyancy, intensifying downdrafts. These processes are the cause of strong downdrafts for the occurrence of microburst [38]. As the volume of liquid water mass and solid water mass precipitates from the upper region of the cloud (as observed in the results), it drags the air together causing the cold pool in the surface.
The concentration of raindrop diameter also varied in the strong downdrafts probably occurred, there was a high concentration of >2.5 mm diameter raindrops in the lower part of the storm and a high concentration of >2 mm raindrops diameter after the base of the 0˚C isotherm. These results corroborated with those observed by JWD in the surface, where raindrops of 2,259 mm diameter were the fourth peak of the highest concentration of raindrops in DSD.
Almost all surface rain during convective events comes from the melting of hail and graupel [39]. leading to a strong convergence zone in the center of the city [7]. When its arrival in the urban region, the SB front stalls over the center of the city for about 2h. Posteriorly, the SB progresses beyond the city when the UHI dissipates [7].

Conclusions
The thunderstorm resulted from the simultaneous interaction of a cold front, low-level jets (northwest flow), sea breeze circulation and large urban area (with heat island) over the MASP. The urban densification (and heat island) in São Paulo city centre induced strong convergence of the air masses of the sea breeze circulation and of the low level jet (northwest), and prevents the advance of the SB front for almost 2 hover the region.
The interaction of the sea breeze front, gust fronts, and the urban densification of São Paulo city center caused triggers of the most intense secondary cell in the rear of the breeze front and resulted in a peculiar microphysics with higher hail production above the 0 o C isotherm than other areas of the thunderstorm.
The mesoscale system (sea breeze circulation) and synoptic scale (frontal system) would not be enough to cause an intense thunderstorm over the region. Anthropogenic action (urbanized areas and heat island) is the main factor for the formation of violent storms over MASP.
These important microphysical processes within long-lasting secondary convective cells should be taken into account on convective parameterization schemes as well as the associated cold pool dynamics.