Wind Structure of a Subtropical Squall Line in China : Results from Dual-Doppler Radar Data

A subtropical squall line moved from Guangxi to Guangdong province in South China on 23-24 April 2007, which resulted in gale and heavy precipitation.The three-dimensional (3D) wind field of the squall line in its mature period was retrieved by GuangzhouShenzhen dual-Doppler data.The 3D conceptual model of this squall line was proposed. On the horizontal plane, the storm-relative front-to-rear inflow prevailed at the lower altitudes of the leading edge. The rear-to-front cold inflow in the stratiform region was observed below 3 km height, which enhanced the convergence in the convective region. At the middle altitudes of the squall line, the front-to-rear horizontal flow prevailed. Strong updrafts were observed at the lower and middle altitudes of the leading edge. Some convergence centers were located at the lower altitudes of the convective region. Furthermore, the storm-relative flow in the vertical cross-section perpendicular to the squall line was revealed. The front-to-rear warm flow extended from the surface to 7.5 km altitude at the leading edge. Above it, part of the front-to-rear inflow blew upward and then forward, and the other part of the inflow blew backward. The descending rear-to-front cold flow was only seen below 3 km height in the stratiform region.

Squall lines often form in the tropical area.The tropical squall line structure has been investigated by dual-radar retrieval wind [2,3,10,[13][14][15]38]. Chong et al. [3] have found that the front-to-rear flow is located at the leading edge, and the rear-to-front flow is located in a layer of 3 km thick in the stratiform region.The mesoscale downdraft in the reflectivity trough [3] (i.e., transition zone) is located in a deep layer from the surface to approximately 6 km, while the downdraft in the stratiform region occurs below the anvil-base altitudes (0-4 km).
The structural characteristics of the midlatitude squall line have been reported in many studies [4,12,18,22,[39][40][41][42].Kessinger et al. [4] have pointed out the differences between a tropical squall line and a midlatitude squall line.The latter showed intense convective cloud along the leading edge, a less-pronounced bright band, and a reversal in the direction of storm-relative flow at the upper altitudes ahead of the squall line.Bright band is a distinct characteristic observed by a weather radar that denotes the freezing level of the atmosphere.The freezing level in a cloud contains ice particles that are coated with liquid water.These particles reflect significantly more radiation than the portions of the cloud above and below the freezing layer.Those characteristics of the bright bands in tropical, subtropical, and midlatitude squall lines are consistent.Houze et al. [9] have proposed the conceptual model of a midlatitude squall line following a dual-Doppler analysis.The storm-relative front-to-rear inflow is observed from the lower to the upper altitudes, and the storm-relative subsiding rear inflow is found at the lower and middle altitudes.
Dual-Doppler radar retrieval data have been applied to investigate the kinematic structure of subtropical squall lines in the Taiwan Strait in Taiwan Area Mesoscale Experiment (TAMEX) [43][44][45][46] and a squall-line-like tropical cyclone rainband [47].Wang et al. [43] have pointed out that the front-to-rear flow prevails at all altitudes at the front of the squall line.The rear-to-front flow prevailing at the rear of the squall line is shallow compared with that of the midlatitude squall line [9].However, the average updraft is slightly stronger than that of a tropical squall line.
By contrast, the study on the three-dimensional (3D) fine structure of the squall line on Mainland China is insufficient.In the past decade, the China Meteorology Administrator (CMA) has set up more than 170 China New Generation Doppler Weather Radars (CINRAD).The CINRAD include S-band and C-band Doppler weather radar.In some areas, some dual-Doppler weather radar networks have been established.On 24 April 2007, a subtropical squall line affected Guangdong province in South China.It caused gale, thunderstorm, and heavy precipitation.In its mature stage, the squall line passed through Guangzhou-Shenzhen dual-Doppler observation domain.It is a good opportunity to study the 3D kinematic structure of this subtropical squall line.What were the structural characteristics of this squall line?What were the new structures compared with the previously studied squall lines in the tropical, subtropical, and midlatitude areas?
This paper focuses on these questions.We begin by describing the data and the dual-Doppler wind retrieval method in Section 2. In Section 3, we first investigate the environment characteristics associated with the squall line and analyse the mesoscale system evolution using the reflectivity mosaic data and the mesonet data simultaneously.Then, we study the horizontal and vertical structural features of the squall line in its mature period using dual-Doppler retrieval wind.Section 4 discusses its characteristics compared with those of other squall lines in the subtropical areas.The 3D conceptual model of the squall line is summarized in Section 5.

Data and Method
The dual-Doppler radar network includes CINRAD/SA (Sband, A-type) Doppler radars in Guangzhou and Shenzhen, which is in the central of Pearl River Delta of China.Figure 1 illustrates the locations of the dual-Doppler radar network, the automatic weather station (AWS) network, and the sounding stations.The two radars are S-band radars with 10-cm wavelength and approximately 1 ∘ beamwidth.The radar parameters are similar to those of Weather Surveillance Radar-1988s (WSR-88D), so are the operation modes.In this case, the two radars operated in Volume Coverage Pattern (VCP) 21 mode.A volume scan is composed of nine constant elevations which vary from 0.5 ∘ to 19.5 ∘ , and the temporal resolution is approximately 5.5 min.
First, the raw radar data quality control was processed by the software system [48].The data quality control included removing noise, clutter, and other incorrect values of radar reflectivity and radial velocity with the methods similar to [49].The automated 2D multipass radial velocity dealiasing scheme has been applied to correct velocity fields [50].Then, the reflectivity mosaic data was generated by the mosaic system [51].The reflectivity fields from individual radars in the radar spherical coordinate were remapped into a common Cartesian grid by bilinear interpolation firstly.The remapped reflectivity fields from multiple radars were then merged to produce a unified 3D reflectivity grid.If the grid cell was covered by one radar only, the analysis value from that radar was assigned to the grid cell.If the grid cell was covered by two or more than two radars, the multiple-radar analysis values were merged to the grid cell by weighting method.The weight was based on the distance between an individual grid cell and the radar location as described in [52].The 3D wind fields were retrieved by the multiple-Doppler synthesis and continuity adjustment technique (MUSCAT) from the dualradar data.Bousquet and Chong proposed the MUSCAT to retrieve 3D wind field using an airborne Doppler radar [53].The comparative study with the extended overdetermined dual-Doppler (EODD) applications shows that MUSCAT provides the more regular description of the airflow.Numerical tests definitely reveal the real improvements of MUSCAT in Doppler radar data retrieval.The statistics of the errors between the retrieved and simulated 3D wind components are small.The mean values are less than 0.2 m/s, while standard deviations remain within 0.7 m/s below 12 km altitude [53].Chong and Bousquet have optimized this algorithm for ground-based dual-radar retrieval and introduced the complex terrain technology [54].Now, it is an operational multiple-Doppler wind retrieval system in France [55].Zhou and Zhang [56] have applied the MUSCAT for groundbased dual-radar wind retrieval and have tested the MUSCAT accuracy using the numerical simulation.The result was the same with the former one [53].Because the 3D wind field deduced from the MUSCAT in the least squares sense does not rigorously satisfy the air mass continuity equation, in order to improve the vertical velocity estimation, a refinement is processed.The refinement is a posteriori integration of the continuity equation to refine the first estimation of the vertical velocity.It consists of adjusting the boundary condition in order to optimize the vertical velocity, followed by upward integration of the continuity equation [57].Previously, we investigated the 3D wind field of a heavy precipitation in South China [58], which is not associated with a squall line, and the 3D structure of a rainfall event along the Meiyu (Baiu, Changma) front [59], using the MUSCAT dualradar retrieval wind.We have revealed that the meso--scale convergence line at the lower and middle altitudes is a cause of the heavy precipitation in both cases.The horizontal scale of this meso--scale convergence line was 40-80 km, and the vertical scale was about 1.5-4.5 km in altitude.
The origin of the wind retrieval coordinate system is at 112.47 ∘ E, 22.37 ∘ N, indicated by "o" in Figure 1.The axis is parallel to the dual-radar baseline at the direction of 127.3 ∘ from the north, the -axis is perpendicular to the baseline, and the -axis points to the direction opposite to the gravity vector.The grid resolutions are 1 km in the horizontal direction and 250 m in the vertical direction.The surface mesonet consists of the AWSs established by the CMA.Two rawinsonde stations are located near the dual-radar network.One is in Qingyuan and the other in Hong Kong.23 April 2007, a temperature trough was showed along Chongqing, Guizhou, and Guangxi provinces at 850 hPa.The region of the squall initiation after 8 hours was ahead of the temperature trough.A pressure trough was over Guangxi, Hunan, and Jiangxi provinces, and the squall line initiation region was on the south of it.The warm, moist southwest wind from South China Sea prevailed in a widespread area ahead of the pressure trough.At 700 hPa, a low-level jet (LLJ) with a maximum speed of 18 m/s prevailed over Guangxi and Guangdong provinces in front of the pressure trough.Behind the trough, the winds were from the northwest with low wind speed.Abundant moisture was available in a wide area ahead of the trough with dewpoint depression of more than 4 ∘ C. At 500 hPa, a shallow trough was shown along Chongqing and Sichuan provinces.The west wind prevailed above the LLJ of 700 hPa.The wind speed was more than 20 m/s over Guangxi and Guangdong provinces.At 200 hPa, a region of divergence occurred over Guangdong province.The contribution of divergence in the upper troposphere contributed to the upward motion in the warm region [43].At 0200 LST 24 April, a quasi-stationary front was extended from China west to Japan with large latitudinal gradient of potential temperature and equivalent potential temperature at 850 hPa.The south wind was on the south of the front and north wind on the north of it.The squall line area was on the south side of the front.A transversal trough laid at 500 hPa in north China (line AB in Figure 2(a)).The shallow trough was still shown along Chongqing and Sichuan (line CD in Figure 2(a)).By 0800 LST, a pressure trough was in the south of Guangxi and Guangdong at 850 hPa.A convergence line was located on the south part of Guangxi, Guangdong, and Fujian provinces.The southwesterly wind prevailed on the south side of the convergence line, and the northwesterly wind was observed on the north of it.The convergence line was also shown at 700 hPa.The LLJ prevailed with the maximum speed of 30 m/s over Guangxi and Guangdong.At 500 hPa, the two troughs moved southward (Figure 2(b)).The wind over Guangxi and Guangdong provinces was from west.This configuration was conducive to convective system development.

Presquall and Postsquall
Conditions.The vertical structure of the atmospheric environment was deduced from the sounding data.Qingyuan and Hong Kong soundings at 0800 LST 24 April were chosen.The former is in the trail of the squall line and the latter in the front-to-rear inflow region of the squall line.Figure 3 shows the skew T-lnP diagram of Hong Kong sounding.The storm-relative wind hodographs of Qingyuan and Hong Kong sounding are illustrated in Figure 4. To calculate the storm-relative motion, it separately averaged the  and V components of the wind at all levels.Then, it added the mean  and V vectors.As shown in Figure 3, the lifting condensation level (LCL) was at 986 hPa and the convective condensation level (CCL) at 960 hPa.The level of free convection (LFC) was low with 946 hPa height.The equilibrium temperature level (ETL) was at 218 hPa.Horizontal winds veered with height below 700 hPa, which indicates warm air advection in the lower troposphere.The southwesterly flow from the South China Sea prevailed at the lower levels with adequate water vapour, while the northwesterly wind from the mainland with cold and dry air prevailed at the upper levels.The vertical shear was about 18 × 10 −3 s −1 between the 0 and 3 km altitudes and 23 × 10 −3 s −1 between the 0 and 5 km altitudes.The moisture was abundant at the lower levels, with the relative humidity below 800 hPa being approximately 5880 5860 5840 5820 5800 5780 5760 85%.There was a convective instability below 3 km.The warm south air contributed to the convective instability.The value of the convective available potential energy (CAPE) was approximately 1676 J/kg.The large CAPE denotes that the presquall environment was convectively unstable and was favourable for squall line development.The water vapour mixing ratio (q) was large below 823 hPa, with a value of 13.87 g/kg at this level.It was 6.16 g/kg at 716 hPa and 0.72 g/kg at 546 hPa (about 5 km altitude).The profile of equivalent potential temperature (e) in the presquall line environment showed convective instability below 625 hPa.The equivalent potential temperature had a maximum value of 352 K at the surface and had a minimum value of 327 K at 625 hPa.The difference was 25 K, which was larger than that of the other subtropical squall line [43].At the same time, Qingyuan is in the stratiform cloud region.Its CAPE value was only 11 J/kg.The LFC in Qingyuan was at 609 hPa and the ETL at 579 hPa.The relative humidity in Qingyuan was small compared with that in Hong Kong.

Surface Observations.
By 0200 LST 24 April, the surface map showed a mesolow in the northwest region of Guangdong that the squall line was initiation.A weak temperature trough was indicated in the vicinity of Hunan and Guangxi.The temperature ahead of the trough was warmer.A convergence line was shown at the northwest region of Guangdong  and moved southeastward.The squall line developed along the convergence line.As indicated in Figure 5, a surface convergence line was presented along the leading edge.The southerly wind from the South China Sea was warm and humid.It was dominant ahead of the convergence line.The northerly flow prevailing in the rear of the convergence line was cold.A weak presquall mesoscale low (1007 hPa) was located on the southwest side of the leading edge.Two thunderstorm mesohighs (1010 hPa), resulting from the cool convective outflow, can be seen along the convective line with heavy precipitation.The weak mesolows were at the back edge of the stratiform precipitation region.The surface backward flow and the environment flow converged, resulting in the convergence line (the short dashed-line in Figure 5) at the front of the weak low center.The horizontal temperature gradient reached a high value across the squall front.
Figure 6 demonstrates the time series of the surface observations at Zhuhai station from 0800 to 1400 LST 24 April 2007.Before the squall line arrival, the temperature and the dew point decreased by approximately 1 ∘ C. The relative humidity dropped by 3%.The wind speed steadied around 2.3 m/s, and the wind direction remained between 180 ∘ and 200 ∘ .The pressure slowly increased when the squall line approached.The squall front arrived in Zhuhai at 1015 LST.In a few minutes, the temperature and dew point dropped by approximately 3 ∘ C and 2 ∘ C, respectively.The relative humidity dropped by 3% rapidly, and the pressure increased by approximately 1.4 mb.The wind direction rotated from south to northeast, with a wind gust of 6 m/s.The precipitation rate reached a peak value of 31 mm/h between 1000 and 1100 LST.The strongest precipitation was 8.7 mm from 1024 to 1030 LST.The accumulative precipitation between 1000 and 1200 LST was 56.8 mm.After 1200 LST, the meteorological parameters were changed as follows.The temperature tended to be uniform at 21 ∘ C over a period of 3 h, and the dew point was at 20.5 ∘ C. The pressure decreased to 1010.8 hPa on 1215 LST and increased to 1011.5 hPa on 1300 LST.After 30 minutes, it decreased to 1010.3 hPa and increased to 1011 hPa on 1400 LST.The wind speed kept an approximate direction of 120 ∘ with a mean value of 4 m/s between 1200 and 1300 LST and then decreased slowly.At the same time, the precipitation rate decreased markedly with a rate close to 1.6 mm/h at 1300 LST.The relative humidity at the surface increased to 95% at 1200 LST as a result of the precipitation.This sequence of the meteorological events was the typical feature at the passage of a squall line [1,2,9,43,46].Figure 6(e) illustrates the time-height cross-section reflectivity at Zhuhai station observed by Guangzhou radar.In this study, the convective cell was defined as the cell with  radar reflectivity of over 40 dBZ.The stratiform region was defined as the wide region with the moderate radar reflectivity (30 to 40 dBZ) [3].Two convective cells passed through Zhuhai station one after another.The first convective cell was more than 4.5 km thick.The core with the reflectivity of over 45 dBZ was about 3 km thick.It lasted for approximately half an hour and produced the first precipitation peak (a peak at 10:20-10:45).The second cell was about 6 km in height, and the core was close to 3 km thick.It produced the second precipitation peak (a peak at 11:10-11:45).The figure shows clearly that the two convective cells resulted in the intense precipitation.
The squall line produced the intense precipitation in a short time in the south part of Guangdong province.The AWS observations showed clearly that the accumulative precipitation in 6 h was more than 20 mm over a large region (Figure 7).There were two heavy precipitation centers with the precipitation of more than 70 mm around the Pearl River Delta.The precipitation in Zhuhai city was 59 mm in 6 h.The precipitation from 1000 to 1200 LST was 56.8 mm in Zhuhai.It also revealed that the precipitation by the squall line was intense and short-lived.

Evolution of the Mesoscale System.
The reflectivity mosaic data with five Doppler weather radars is used for investigating the evolution of the squall line.These radars are all CINRAD/SA, located in Guangzhou (GZ), Shaoguan (SG), Shantou (ST), Shenzhen (SZ), and Yangjiang (YJ). Figure 8 illustrates the mosaic data at 2.5 km and AWSs wind from 0400 to 1300 LST 24 April 2007.2.5 km is a representative altitude that is free from terrain effect.
At 0400 LST, a southwest-northeast echo band with 35 dBZ developed in the adjacent region between Guangxi and Guangdong provinces (Figure 8(a)).The strong cells of over 50 dBZ were observed in the echo band.Some Advances in Meteorology intense cells developed around Shaoguan radar station.Both the strong cells and the echo band developed quickly and moved southeastward (Figure 8(b)).The maximum reflectivity embedded in the convective cloud was over 55 dBZ.After one hour (Figure 8(c)), the leading edge arrived at Qingyuan station.The surface wind speed at Qingyuan station was more than 12 m/s.Some intense cells with 50 dBZ lined up along the leading edge.The squall line continued to move southeastward.At 0700 LST (Figure 8(d)), the squall line was more than 250 km in length in the west-east direction.The rear of the squall line was the wide stratiform region with the moderate radar reflectivity (30 to 40 dBZ), while the leading edge was a narrow convective region of over 30 dBZ.Between the stratiform region and the convective region, a narrow weak reflectivity band under 30 dBZ formed.It was similar to the reflectivity trough [3].A convergence line was located at the surface.At the front part of the squall line, the southerly wind prevailed.On the other hand, the northwesterly wind was observed behind the squall front.The convective cloud was developing quickly in the convergence area.The squall line moved southeastward and passed through Heyuan sounding station and Guangzhou radar station (Figure 8(e)).The surface convergence line moved southeastward, too.The wind speeds at the two stations exceeded 10 m/s and 6 m/s, respectively.The squall line developed continually.It was more than 500 km in length in the west-east direction at 0900 LST (Figure 8(f)).Seven intense cells of over 50 dBZ were shown at its front.The rear stratiform region broadened.At 1000 LST (Figure 8(g)), the squall line moved into the dual-Doppler radar wind retrieval domain and continued to develop.Its length was approximately 600 km in the westeast direction.The asymmetric cloud structure of the squall line was clear.The stratiform cloud continued to broaden in the east part of the squall line.On the other hand, the Advances stratiform cloud in the west part of the squall line (the region around Qiangyuan station) was relatively weak.At 1100 LST (Figure 8(h)), the squall front moved to the south part of the dual-radar retrieval domain.As illustrated in Figures 8(i) and 8(j), the squall line moved to the South China Sea with weak reflectivity compared with the former one finally.

Dual-Doppler Analysis.
The 3D wind fields were retrieved by the MUSCAT using Guangzhou-Shenzhen dual-Doppler synchronous volume scan data.At the time of the analysis (1000 to 1041 LST), the squall line was moving at a nearly constant speed of 12 m/s in the direction of 112.5 ∘ .The mean motion was deduced from the reflectivity data in the successive volume scans by Guangzhou radar.The storm-relative horizontal wind was subsequently calculated by subtracting the squall line propagation velocity from the dual-Doppler retrieval wind.

3.3.1.
Wind Field on the Horizontal Plane. Figure 9 presents the radar reflectivity and storm-relative wind field at 1005 LST 24 April 2007.At the lower altitudes, the leading edge was oriented in the southwest-northeast direction.Eight intense cells of over 45 dBZ were located along the leading edge at the altitude of 1.75 km (Figure 9(a)).The maximum reflectivity was over 50 dBZ in the northeast part of the leading edge.
The storm-relative front-to-rear horizontal inflow and the storm-relative rear-to-front cold inflow were illustrated in Figure 9(a).The storm-relative front-to-rear horizontal inflow was in the right of the convergence line (line AB in Figure 9(a)) and prevailed at the leading edge.The region of large presquall line wind vectors was located at  < 60 km.Wind vectors were nearly paralleling the line at a small angle intercepting the squall line.The strongest inflow was about 17 m/s.The reflectivity gradient was intense at the leading edge.The storm-relative rear-to-front cold inflow was in the left of the convergence line.The two flows eventually converged at the convective line (line AB in Figure 9(a)), resulting in the updraft (Figure 9(c)) and convergence band (Figure 9(e)).Some convergence centers were observed in the convergence band (Figure 9(e)).The absolute value of the strongest convergence was more than 30 × 10 −4 s −1 .Above 3 km, the storm-relative front-to-rear horizontal inflow was stronger than that of the low levels at the leading edge.
The surface observation at 1000 LST (Figure 8(g)) indicated that the southerly flow prevailed on the front part of the squall line.On the other hand, the northwesterly flow was observed in the stratiform region.On the whole, the dual-Doppler retrieval at the lowest attitude is consistent with the surface wind patterns.
Qingyuan sounding station was in the stratiform region by LST.The sounding showed that the northwest flow prevailed below 3 km altitude.Hong Kong sounding station was in the convective region.The sounding indicated that southeast wind and southwest wind adequate water vapour were shown below 3 km height (Figure 3).These wind patterns were in accordance with dual-radar retrieval on the whole.
The VWP (VAD Wind Profile) reveals the horizontal wind vertical profile near the radar station region.Guangzhou radar was in the stratiform region at 1005 LST.The north wind prevails below 0.9 km altitude, with speed between 4 and 12 m/s.It changed to west direction anticlockwise from 0.9 to 3.7 km altitude.The west wind was 22 m/s at 5.2 km altitude.Above it, the southwest prevailed.The wind direction was consistent with dual-radar retrieval.
Figure 10 shows the reflectivity and the storm-relative radial velocity observed by Guangzhou radar at 1005 LST.At the elevation of 0.5 ∘ (Figure 10(b)), the rear-to-front inflow prevailed in the large region of dual-radar analysis domain.At the leading edge, the front-to-rear flow was showed.The convergence line was at 1.2 km altitude (line AB).At the layer of 1.5 ∘ (Figure 10(c)), the convergence line (line AB), which was corresponding to the convection line (Figure 10(a)), was at 2.4 km height.At the elevation of 2.4 ∘ (Figure 10(d)), the convergence line (line AB) was at 3 km altitude in the dualradar domain.The convergence line was at 3.5 km altitude at the layer of 3.4 ∘ (Figure 10(e)).At the elevation of 4.3 ∘ (Figure 10(f)), the divergence was at 6.8 km altitude.The maximum front outflow was 14 m/s, and the maximum rear outflow was about −20 m/s.Compared with the reflectivity at the low layers, part of the front outflow was in the overhang echo.The storm-relative radial velocity was associated with dual-radar analysis.
Updrafts were observed in the convective region (Figure 9(c)), corresponding to the convergence region and intense reflectivity.The leading updrafts were associated with the storm-relative horizontal southeasterly inflow and spread to the middle and upper altitudes along the squall front.Many intense updraft cores were located along the leading edge.The maximum convective updraft was about 7 m/s (Figure 9(d)).On the other hand, the downdraft zones were located in the stratiform region.The downdrafts were weaker than the updrafts at the same altitude.The updraft and downdraft were weaker than those of a tropical squall line [3].
The reflectivity at 5 km (Figure 9(b)) was weaker than that at the lower altitudes.The intense reflectivity band over 45 dBZ was narrower than that at the lower altitudes, too.However, the maximum reflectivity was still over 50 dBZ.The storm-relative front-to-rear horizontal wind prevailed in both convective cloud and tailing stratiform region.The wind speed was larger than that at the lower altitudes.The maximum speed was about 29 m/s in the southwest part of the domain ( < 30 km).Some strong convergence centers less than −20×10 −4 s −1 were observed in the convective region (Figure 9(f)).
Figures 11 and 12 present the reflectivity and stormrelative wind field at 1023 and 1041 LST, respectively.The leading edge preserved its shape below 4 km altitude.Although the echo band of over 30 dBZ split into two parts at 5 km altitude at 1023 LST (Figure 11(b)), the reflectivity structure (Figures 11(a

Wind Field Structure along Vertical Cross-Section
(1) Radar Reflectivity Structure.Figure 13(a) shows the reflectivity and storm-relative horizontal flow contours in the vertical cross-section perpendicular to the squall line.The vertical distribution of the radar echo had two contrastive structures, that is, the convective cloud in the front of the system and the stratiform cloud in the rear.The overhanging radar echo was observed above 7.5 km ahead of the leading edge.This pattern is also observed in the subtropical [60] and midlatitude [9] squall line but is not in some tropical squall lines [1,3].The storm-relative front-to-rear wind pattern, which was different from those in the literatures, formed the overhanging radar echo as discussed later in Section 3.3.2(4).The reflectivity in the convective region was intense.The maximum reflectivity was over 50 dBZ, but it was still weaker than that of a tropical squall line [3].In the latter, the air is moister and more unstable at the low and middle levels.The water vapour mixing ratio profiles in the two cases were different.In the latter, it is more than 8 g/kg at 3 km altitude and 5 k/kg at 5 km altitude.The ETL is 13 km altitude, which is 1 km higher than that of our case as shown in Section 3.1.2.As a result, the reflectivity of over 50 dBZ can extend from the surface to the altitude of 5 km in the latter.
Three convective cells over 40 dBZ lined up along the squall line from left to right.Two cells on the right were higher than the left one.The first cell ( = 80km) was approximately 6.5 km in height, and the reflectivity of over 50 dBZ was located at the altitude of 2 km.These cells were vertical in general, and only the cell near the stratiform region was backward slightly.The vertical structure is dependent on the development stage [61] and the environment kinematic condition.The upright structure was similar to that of midlatitude [9,18] and other subtropical [60] squall lines.
The CAPE, which is a measure of the amount of buoyant energy available in the environment to drive updrafts, is 2300 J/kg [62] in a midlatitude squall line [9].Its environmental instability is more intense than that of our case.This updraft derived by the environmental instability transports moister air to higher region.As a result, the reflectivity with 30 dBZ is dominant compared to 13.4 km altitude in the midlatitude case [9].In our case, it was approximately 9.5 km in height, which was relatively shallow.On the contrary, the convective cell in our case was stronger than that of the TAMEX subtropical squall line [43] but weaker than that of Taiwan other squall lines reported by Chang et al. [60].
The environment instability and moisture contribute to the reflectivity structure.
(2) Storm-Relative Horizontal Airflow Structure.As shown in Figure 13(a), the storm-relative front-to-rear horizontal flow entered the system from the surface to the altitude of 7.5 km at the leading edge, which was associated with intense upward motion.The intense upward motion was defined as the value over 5 m/s and the intense downward motion as the value under −2 m/s.Above 7.5 km of the squall front, part of the flow slopped forward and outward and the other part of the flow rearward.The front-to-rear flow prevailed in the stratiform region and the rear part of the convective region above 3 km.The intense flow was observed at the middle altitudes (3.5-8 km) mainly, with the speed of more than 10 m/s.This rearward flow speed attained the maximum value of more than 20 m/s at the altitude of 3.5 km ( = 73 km), while the rearward flow approaches the maximum of 27 m/s in a midlatitude squall line [9].This inflow was important for the horizontal transports of momentum and moisture into the stratiform region.
The rear-to-font horizontal dry cold inflow entered the stratiform region below 3 km height and descended from the stratiform region to the convective region.It reduced the hydrometeor mass by evaporation and sublimation [5].This inflow played an important role in feeding the cold air.Moreover, it reached the maximum speed of approximately 10 m/s at the lower altitudes of the second convective cell and finally reached the rear of the third convective cell at the lower altitudes ( = 80 km).As the descending air approaching the surface, it collided with the front-to-rear low-level environmental inflow and formed a gust front.The interaction between the front-to-rear environmental low-level inflow and the rear-to-front environmental midtropospheric flow was responsible for maintaining the updraft-downdraft in the convective region [43].The depth of the rear-to-front inflow depends on the environmental vertical wind shear.Although the depth in our case was similar to that of the subtropical squall line in the TAMEX [43], the inflow speed was stronger in our case.It was shallower than that of the tropical [1] and midlatitude [9] squall lines.In the midlatitude squall line with higher shear, the inflow can also reach the height of 7.5 km, and the wind speed is stronger than that of our case.cells and the formation of the heavy precipitation.The interaction between the convective updraft and downdraft played an important role in maintaining the 3D circulation within the squall line [43].As shown in the sounding, the LLJ in our case was more intense than that of other subtropical squall line [43].The CAPE is 1369 J/kg [43] which is smaller than that of our case.The environmental instability (/) in our case from surface to midlevel (600 hPa) was 5 K/(100 hPa), which was more unstable than that of the other subtropical squall line [60].As higher vertical wind shear and higher environmental instability, the updraft was stronger than those of other subtropical squall lines [43,60].In the midlatitude squall line [9], the CAPE and the environmental instability shown in Section 3.3.2(1)are more intense than that of our case.Furthermore, the lifting in the vicinity of surface front helps to initiate the squall line.As a result, the updraft in the convective region is 13.4 km altitude, which is stronger and deeper than that of our case.The cores of the updrafts decreased in intensity, which was similar to the result reported in [9].The layer of front-to-rear inflow in the convective cell was dominated by upward motion, which is also found in the midlatitude [9] and other subtropical [43,60] squall lines.There was a shallow and weak downdraft with the speed of 2 m/s at the lower altitudes behind the second convective cell in our case.The updraft and downdraft were observed alternately in the stratiform region.The front-to-rear flow transported the unstable and moist presquall line air into the convective cells.The updraft further transferred the water vapour into the upper levels of the troposphere.Convergence was prominent at the low and middle levels of the convective region.Flow above 7.5 km was divergent in the zone of convective updraft.As the result of divergence from the convective updrafts, part of the flow carrying the particles sloped gradually into upstream and formed trailing stratiform precipitation [60].Other part of this flow diverted forward toward the leading.It transported the particles downstream and played an important role in forming and maintaining the overhanging radar echo at the middle and upper altitudes of the leading edge.The overhanging wind structure depended on the vertical wind structure of the environment relative to the orientation of the squall line.This flow distribution was different from that of the tropical [1,3] and other subtropical [46] squall lines.In the previous studies, the storm-relative front-to-rear wind directs into the squall lines from the front at all levels [1,3,46], with a relative minimum appearing at 700-600 mb generally [3].
The storm-relative rear-to-front cold flow prevailed below the altitude of 3 km only in the trailing stratiform region.It entered the back edge of the convective cloud at the lower altitudes.The midtropospheric air was cooler than the surrounding air.Its buoyancy was further reduced due to the precipitation loading in the high reflectivity region behind the leading edge.The rear-to-front flow reached the convective line and enhanced the convergence at the lower altitudes.An enclosed circulation was observed at the lower and middle altitudes of the squall line.The interaction between the frontto-rear warm flow and the rear-to-front cold flow at the lower altitudes is conducive to maintain the updraft-downdraft structure in the convective region.

Discussions
As depicted above, many structural characteristics of this squall line were similar to those of the other subtropical squall line [43].The squall line was composed of a well-organized convective cloud region, a weak reflectivity transition region under 30 dBZ, and a broad stratiform region.The stormrelative front-to-rear wind entered the squall line at the lower and middle altitudes of the leading edge, and the stormrelative rear-to-front wind entered the system at the lower altitudes of the stratiform region.
Our results, however, had presented some important kinematic characteristics different from the reported squall lines.The echo top height is defined as the altitude of 20 dBZ [3,18].First, the echo top in our case was approximately 10 km, and it is 8 km height in some squall lines in the Taiwan Strait [43,46].In the other squall line of the Taiwan Strait [60], the echo top is more than 13 km.In our case, the overhanging radar echo was observed above 7.5 km at the leading edge.This echo feature is not seen in the other subtropical squall line [43].The vertical wind shear in the environment contributes to the echo structure.
Second, the storm-relative front-to-rear horizontal inflow entered the squall line from the surface to the altitude of 7.5 km at the leading edge and then converted to two outward flows.One was forward and the other rearward.In the other subtropical squall lines [43,46], the front-to-rear wind prevailed at all altitudes.
Third, the rear-to-font cold inflow entered the squall line from the stratiform region below 3 km only.On the other hand, the inflow prevails at all altitudes in the stratiform region in the other subtropical squall line [60].The speed was stronger in our case, although the wind pattern was similar to that of a subtropical line in the same latitude zone [43].The vertical wind structure of the environment contributes to the flow pattern.
Moreover, as shown in Figure 13(b), the intense updrafts prevailed in the squall front and extended to 10 km in height.The updraft and downdraft were stronger than those of the other subtropical squall line [43].

Conclusions
The goal of this study is to reveal the evolution features of a subtropical squall line in South China the mature period.The retrieval results demonstrated the possibility to perform the dual-Doppler retrieval analysis from operational radars, too.We have presented some new structures as demonstrated in Section 4.
A conceptual model of this squall line in the cross-section perpendicular to the squall line is summarized in Figure 14.The squall line included a convective region, a weak reflectivity transition region under 30 dBZ, and a trailing stratiform region.The convective region was intense.The storm-relative front-to-rear warm flow entered the squall line from the surface to the altitude of 7.5 km at the leading edge, which was associated with intense upward motions.Above 7.5 km height, part of the inflow sloped gradually into the trailing stratiform cloud; part of the front-to-rear inflow sloped into the leading convective cloud and then moved forward.Below the altitude of 3 km in the wide trailing stratiform region, the descending shallow rear-to-front cold inflow entered the squall line, playing an important role in feeding the cold air.The rear-to-front inflow blew forward, and finally part of the cold inflow converged with the frontto-rear deep warm inflow at the convective line.It resulted in intense updraft and strong convergence at the lower altitudes.New cells formed at the squall front, due to the front-to-rear inflow at the lower altitudes of the gust front.It should be of interest and importance to know how the dynamics within the squall line connected with microphysics.This topic will be explored by numerical simulation in future.The squall line studied in the present was only a particular case.In order to summarize the general features of the squall line over South China, more cases should be studied.Climate (CNR-ISAC) for his insightful comments and constructive suggestions to improving this manuscript.

Figure 1 :
Figure 1: Map of the dual-Doppler radar observation network and the mesonet weather stations.The colour shading shows the composite reflectivity (dBZ) on 1017 LST 24 April 2007 and the gray shading for the altitude (m).The two large circles represent Guangzhou (GZ) and Shenzhen (SZ) radars scan domains, respectively.Rectangle A shows the dual-Doppler wind retrieval domain, and "o" indicates the origin of the wind retrieval coordinate system.ZH represents Zhuhai weather station. Feature

Figure 2 :Figure 3 :
Figure 2: National Centers for Environment Prediction (NCEP) geopotential height analyses at 500 hPa from 0200 LST (a) and 0800 LST (b) 24 April 2007.The blue rectangles indicate the position of the squall line and the thick solid lines are for the trough.

Figure 4 :
Figure 4: Storm-relative wind hodograph from the upstream sounding (Qingyuan: magenta line) and downstream sounding (Hong Kong: green line) at 0800 LST 24 April 2007.Every ring presents 20 knots of increment in speed, and the numbers indicate the pressure levels (hpa).

Figure 5 :
Figure 5: Surface chat at 0800 LST 24 April 2007.The colour shading is radar reflectivity (dBZ) at 2 km height.Solid contours are pressure field (hPa).Dashed lines indicate the two surface convergence lines.Temperature and dew point ( ∘ C) are given at individual weather station.

Figure 9 :
Figure 9: (a-b) The storm-relative horizontal wind fields at 1.75 and 5 km, respectively, at 1005 LST 24 April 2007.(c-d) The vertical velocities at 1.75 and 5 km, respectively.The solid lines indicate the updraft and the dashed lines for the downdraft.(e-f) The divergence field (10 −4 s −1 ) at 1.75 and 5 km, respectively.Solid lines are for positive values (divergence) and dashed lines for negative value (convergence); the contour interval is 5 × 10 −4 s −1 .The -axis and -axis are, respectively, parallel and perpendicular to the radar baseline.The colour shading is radar reflectivity (dBZ).The squall line movement direction (C) and the north (N) are indicated in (a).Line AB in (a) denotes the convergence line.

Figure 11 :
Figure 11: (a-b) The storm-relative horizontal wind fields at 1.75 and 5 km, respectively, at 1023 LST 24 April 2007.(c) The vertical velocity at 1.75 km.The solid lines are the updraft, and the dashed lines are the downdraft.(d) The divergence field (10 −4 s −1 ) at 1.75 km.Solid contours are for positive values (divergence) and dashed contours for negative value (convergence); the contour interval is 5 × 10 −4 s −1 .The colour shading is radar reflectivity (dBZ).

( 3 )Figure 13 :
Figure 13: Reflectivity, the storm-relative horizontal air flow (a), the vertical velocity (b), and the storm-relative wind field (c) in the vertical cross-section perpendicular to the squall line (along AB in Figure 11(a)).

( 4 )
Storm-Relative Wind Field Structure.As shown in Figure 13(c), the storm-relative front-to-rear deep ascending wind began at the lower altitudes of the leading edge and extended upward through the convective region where the strong updraft over 10 m/s was located (Figure 13(b)).The updraft was forced by the cold rear inflow at the low levels.

FFigure 14 :
Figure 14: A 3D conceptual model of the squall line.White isosurfaces present the radar reflectivity of over 40 dBZ.The horizontal colour shading presents the radar reflectivity at 1.75 km.The vertical colour shading is the reflectivity in the vertical cross-section perpendicular to the squall line.The white thin solid arrows in the vertical cross-section are the storm-relative wind field.The arrow indicates the north (N).The 3D gray arrow presents the storm movement direction.