删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Recent Progress in Dual-Polarization Radar Research and Applications in China

本站小编 Free考研考试/2022-01-02

Kun ZHAO1,2,*,,,
Hao HUANG1,2,
Mingjun WANG1,2,
Wen-Chau LEE3,
Gang CHEN1,
Long WEN1,
Jing WEN1,
Guifu ZHANG1,4,
Ming XUE1,4,
Zhengwei YANG1,
Liping LIU2,
Chong WU2,
Zhiqun HU2,
Sheng CHEN5

Corresponding author: Kun ZHAO,zhaokun@nju.edu.cn;
1.Key Laboratory of Mesoscale Severe Weather of Ministry of Education and School of Atmospheric Sciences, Nanjing University, 163 Xianlin Road, Nanjing 210023, China
2.State Key Laboratory of Severe Weather and Joint Center for Atmospheric Radar Research of China Meteorological Administration and Nanjing University, Chinese Academy of Meteorological Sciences, Beijing 100081, China
3.National Center for Atmospheric Research, Boulder, Colorado 80301, USA
4.School of Meteorology and Advanced Radar Research Center, University of Oklahoma, Norman, Oklahoma 73019, USA
5.School of Atmospheric Sciences, and Guangdong Province Key Laboratory for Climate Change and Natural Disaster Studies, Sun Yat-sen University, Guangzhou 510006, China
Manuscript received: 2019-03-21
Manuscript revised: 2019-06-07
Manuscript accepted: 2019-06-11
Abstract:Dual-polarization (dual-pol) radar can measure additional parameters that provide more microphysical information of precipitation systems than those provided by conventional Doppler radar. The dual-pol parameters have been successfully utilized to investigate precipitation microphysics and improve radar quantitative precipitation estimation (QPE). The recent progress in dual-pol radar research and applications in China is summarized in four aspects. Firstly, the characteristics of several representative dual-pol radars are reviewed. Various approaches have been developed for radar data quality control, including calibration, attenuation correction, calculation of specific differential phase shift, and identification and removal of non-meteorological echoes. Using dual-pol radar measurements, the microphysical characteristics derived from raindrop size distribution retrieval, hydrometeor classification, and QPE is better understood in China. The limited number of studies in China that have sought to use dual-pol radar data to validate the microphysical parameterization and initialization of numerical models and assimilate dual-pol data into numerical models are summarized. The challenges of applying dual-pol data in numerical models and emerging technologies that may make significant impacts on the field of radar meteorology are discussed.
Keywords: dual-polarization radar,
quantitative precipitation estimation,
precipitation microphysics,
drop size distribution,
numerical model
摘要:同常规多普勒雷达相比, 双偏振雷达可测量更多反映降水系统微物理信息的参数,因此被广泛用于研究降水微物理特征和改进雷达定量降水估测. 本文总结了我国近期双偏振雷达研究和应用的进展. 首先, 回顾了我国一些代表性的双偏振雷达特性和雷达数据质量控制方法, 包括雷达标定、衰减订正、比差分传播相移的计算,以及非气象回波识别和去除. 基于双偏振雷达的雨滴谱反演、水凝物相态分类和降雨估测产品, 揭示了我国典型降水系统内部的微物理特征和过程. 同时, 总结了双偏振雷达观测在我国数值模式微物理参数化方案评估、资料同化和模式初始场改进中的应用. 最后, 讨论了利用双偏振雷达观测改进数值模式面临的挑战和天气雷达技术发展的趋势, 及其对雷达气象学领域研究的影响.
关键词:双偏振雷达,
定量降水估测,
降水微物理,
雨滴谱,
数值模式





--> --> -->
1. Introduction
The first dual-polarization (dual-pol) research radar was developed in the late 1970s in the United States (US) (Seliga and Bringi, 1976, 1978). By adding a vertically (V) polarized channel in addition to the horizontally (H) polarized channel to a conventional radar, dual-pol radars can measure more parameters, including differential reflectivity (ZDR), co-polar cross-correlation coefficient (ρhv), differential phase ( ΦDP), and specific differential phase (KDP), which can provide more microphysical information of precipitation systems (shape, phase, and type of hydrometeors) (e.g., Doviak and Zrni?, 1993; Bringi and Chandrasekar, 2001; Zhang, 2016). In practice, the majority of dual-pol radars transmit in simultaneous H and V mode, for which only the co-polar signals are measured. For several research dual-pol radars, the cross-polar signals can be measured by transmitting alternating H and V modes with two receivers. The US operational radar network, weather surveillance radar (WSR)-88D, completed the upgrade to dual-pol radars in 2013 (Kumjian, 2013). Other regions, including Europe, Canada, Japan, Korea, and China have begun to upgrade their operational radar networks in recent years. The theory behind dual-pol radars and the applications of dual-pol radar data in understanding microphysical processes and quantitative precipitation estimation (QPE) can be found in textbooks (e.g., Doviak and Zrni?, 1993; Bringi and Chandrasekar, 2001; Zhang, 2016) and review articles (e.g., Bluestein et al., 2014; Hubbert et al., 2018; Zhang et al., 2019).
China is one of the countries in the world that suffers severe damage from high-impact weather [e.g., mesoscale convective systems (MCSs), tropical cyclones etc.] and the accompanying flooding and mudslides. To improve the ability in monitoring and nowcasting these high-impact weather events, China has deployed a nationwide radar network since 1999 composed of more than 200 China's New Generation Doppler Weather Radars (CINRAD 98D) in both S- and C-bands. In recent years, with the advent of dual-pol radar technologies in China, several X- and C-band mobile dual-pol radars have been developed by universities and research institutes, including Nanjing University's C-band radar (NJU CPOL) and the Institute of Atmospheric Physics' X-band radar (IAP XPOL). These radars have been used to observe severe weather in field campaigns such as the Observation, Prediction and Analysis of Severe Convection of China (OPACC) (Xue, 2016), the Southern China Monsoon Rainfall Experiment (SCMREX) (Luo et al., 2017), the Understanding and Prediction of Rainfall Associated with Landfalling Tropical Cyclones (UPDRAFT) (Wang et al., 2018a), and the Third Tibetan Plateau Atmospheric Scientific Experiment (TIPEX-III) (Zhao et al., 2018). In 2013, the first operational dual-pol radar (CINRAD 98DP) in China was in service in Zhuhai, Guangdong Province. Since then, more than 10 S-band weather radars in Shanghai City and the provinces of Guangdong, Fujian, and Anhui have been upgraded to dual-pol radars, and more than 100 dual-pol radars will be built or upgraded in China by 2020. To date, the data collected from these dual-pol radars in China have been used and analyzed to improve understanding of precipitation microphysics (e.g., Wang et al., 2016b; Wen et al., 2017) as well as radar QPE and quantitative precipitation forecast (QPF) in China (e.g., Chen et al., 2017; Huang et al., 2018a).
The purpose of this paper is to review the progress and status of dual-pol radars, research results and applications in China, including (1) the characteristics of dual-pol radars and data quality control procedures, (2) rainfall estimation and microphysical retrieval methods from dual-pol radars, (3) precipitation and microphysical characteristics in severe weather deduced from dual-pol radars, and (4) applications of dual-pol radar data in numerical models.

2. Characteristics of dual-pol radars and data quality control in China
The characteristics of CINRAD-98D with polarimetric capacity (CIRAD-98DP) and several C- and X-band mobile dual-pol radars in China are summarized in Table 1. Three polarimetric variables are important for radar applications in terms of precipitation and microphysical characteristics: ZDR, which is a function of drop shape and a good measure of the median drop diameter; KDP, which is more linearly related to rain rate (R) than equivalent reflectivity factor at H polarization (ZH) and is immune from radar calibration, attenuation, and partial beam blockage; and ρhv, which is a measure of the diversity of particles (e.g., types, shapes and/or orientations) within a radar sampling volume. Meteorological targets typically possess ρhv>0.8, where ρhv for pure rain and snow is close to 1.
For QPE or hydrometeor classification (HC), the desired accuracies for ZH and ZDR are 1 dB and 0.1-0.2 dB, respectively (Bringi and Chandrasekar, 2001; Chandrasekar et al., 2015). Common procedures to calibrate ZH and ZDR include using a test signal, sun scans, and backscatter signal from a metallic sphere (Hubbert et al., 2003; Ryzhkov et al., 2005b; Zrnic et al., 2006). ZDR can also be calibrated using signals from drizzle and dry aggregated snow (Ryzhkov et al., 2005b). The accuracy of ZH can also be benchmarked against the values calculated from ZDR and KDP according to the self-consistency of radar variables (Gorgucci et al., 1992; Vivekanandan et al., 2003). These methods have been widely utilized in recent studies in China (Du et al., 2012, Du et al., 2013; Hu et al., 2014; Huang et al., 2017). Generally, different ZDR calibration methods result in similar results. However, the quality of rotary joints can degrade over time and use, and influence the consistency of the dual channels of dual-pol radars. Thus, the bias of ZDR can vary with azimuth, which should also be considered in ZDR calibration (Chen et al., 2018; Hu et al., 2018).
S-band dual-pol radars are less subject to attenuation in heavy rainfall and are primarily used for the national weather surveillance radar network in China, similar to those used in the US. However, X- and C-band dual-pol radars are also widely used in China (e.g., citywide radar networks and research radars), which require attenuation correction for quantitative or even qualitative applications when using algorithms developed for S-band radars (Carey et al., 2000). ΦDP and KDP are mostly unaffected by attenuation despite it being able to reduce the signal-to-noise ratio (SNR) and degrade the quality of ΦDP. ΦDP can be utilized to improve the accuracy of attenuation correction (Bringi and Chandrasekar, 2001). Two approaches for attenuation correction are commonly used in China: (1) specific attenuation (AH) and specific differential attenuation (ADP) are linearly related to KDP (Bringi et al., 1990; Lei, 2014; Wu and Huang, 2014; Huang et al., 2018b), and (2) the "ZPHI" rain-profiling algorithm, where AH is derived from attenuated reflectivity factor under the constraints of total path-integrated attenuation (Testud et al., 2000). To mitigate the impacts of raindrop size distribution (DSD) variabilities, the ratio α used in the ZPHI algorithm can be adaptively adjusted using the self-consistency of polarimetric variables.
KDP is widely utilized in QPE because the relationship between KDP and R is less affected by the DSD variabilities than Z-R relationships (Bringi and Chandrasekar, 2001). However, KDP is not directly measured; instead, it is estimated from the range derivative of filtered ΦDP (Hubbert et al., 1993; Hubbert and Bringi, 1995) to avoid random errors in Φ DP propagating to KDP and producing erroneous negative values for rain or for rain mixtures. Thus, Φ DP is usually first processed/filtered using running average method (Wei et al., 2014), median average method (Wei et al., 2014; Wu et al., 2017), finite-impulse response filter (Hubbert et al., 1993; Hubbert and Bringi, 1995), Kalman filter (Wang and Chandrasekar, 2009), wavelet analysis (Hu and Liu, 2014), or the linear fitting and recurrence method (Sun et al., 2015). (Huang et al., 2017) proposed a hybrid method by combining the physical constraints of KDP (calculated from ZH and ZDR using self-consistent relationship, Vivekanandan et al., 2013) and the linear programming algorithm (Giangrande et al., 2013) to improve the estimates of KDP and rainfall.

3. Precipitation characteristics retrieved from dual-pol radars
2
3.1. Hydrometeor Classification
--> Dual-pol radar is capable of identifying the primary hydrometeor type in a radar sampling volume because hydrometeor size, shape, orientation, phase, and bulk density affect dual-pol radar observables to different degrees. The relationships between the distributions of dual-pol radar measurements and the hydrometeor types overlap and are not well defined. Therefore, a fuzzy logic approach, which assigns membership functions for each radar observable to account for the overlapping and soft boundaries, has been widely used in HC for both research and operations (Vivekanandan et al., 1999; Park et al., 2009; Dolan et al., 2013). These membership functions should be tuned for different weather regimes, geographical regions, and type of radars. In China, most HC studies adjusted the membership functions of existing methods, (e.g., Park et al., 2009), for different regions (Wu et al., 2018a) and different radar frequencies (Gu et al., 2015; Ran et al., 2017; Feng et al., 2018). In (Wu et al., 2018a), the method of (Park et al., 2009) was used for radars in South China. Due to the differences in polarimetric characteristics of some hydrometeors between China and the US, applying US HC membership functions in South China results in insufficient discrimination of aggregate values. Discontinuities are found in hail, graupel, wet snow, and heavy rainfall areas (within the dotted line in Fig. 1b). By tuning (statistics-based optimization) the membership functions, the HC results in Fig. 1c are more coherent.
Figure1. Vertical structure of a squall line in South China observed by Zhuhai dual-pol radar, 10 May 2014: (a) horizontal reflectivity; (b) HC based on fuzzy logic; (c) optimized HC. The colors in (b, c) represent different classes of scatterers, including ground clutter or anomalous propagation, biological scatterers, dry snow, wet snow, crystal, graupel, big drops, rain, moderate light and moderate rain, heavy rain, and hail or the mixture of rain and hail.


Considering the limitation of the fuzzy logic-based HC method, statistical decision theories, e.g., the maximum likelihood and Bayesian theory, have also been applied for HC in China in recent years (Marzano et al., 2008), where the hydrometeor types are determined using a posteriori probability. The statistical information can also be used to constrain HC by using the a priori distribution. In the work of (Wen et al., 2015) and (Wen et al., 2016), the conditional probability distribution of the polarimetric variables and ambient temperature corresponding to different hydrometeor types were derived by applying clustering techniques, and were successfully used for the HC of hailstorms and shallow Arctic mixed-phase clouds. (Yang et al., 2017) recently proposed a Bayesian-based HC algorithm, in which the conditional probability functions of polarimetric variables are constructed for seven different hydrometeor types. Since the method is statistically trained using radar observations in China, it has been proven to produce more reasonable hydrometeor types than the fuzzy logic method for a squall line event that occurred on 30 July 2014 in eastern China. It could be a promising way to achieve HC for dual-pol radar measurements.

2
3.2. QPE
--> Many dual-pol radar rainfall estimators have been developed, including R(ZH,ZDR), R(KDP), and R(KDP,ZDR) (Ryzhkov and Zrni?, 1995; Gorgucci et al., 2001; Ryzhkov et al., 2005a; Lee, 2006; Bringi et al., 2011), and they have yielded better rainfall estimation than the conventional Z-R relation, R(ZH), particularly for moderate and heavy rain. Dual-pol radar rainfall estimators mainly suffer from uncertainty in two aspects: the model errors caused by DSD variabilities and the measurement errors. To make rainfall estimators more consistent with microphysical climatology in China, DSDs derived from disdrometers have been used to tune dual-pol radar rainfall estimators, and these estimators have been widely applied and evaluated in China (Gao et al., 2014; Zheng et al., 2014; Wei et al., 2016; Chen et al., 2017; Zhang et al., 2017c). Among them, R(KDP) provided the best rainfall estimation for X- and C-band radars that are susceptible to severe attenuation in heavy precipitation (Wei et al., 2016; Chen et al., 2017).
For light rain, the advantage of polarimetric rainfall estimators over the conventional Z-R relationship diminishes, because measurement errors carry a greater weight than the useful information contained in ZDR and KDP. To improve rainfall estimation, (Chen et al., 2017) proposed a new composite rainfall estimator, R(ZH,KDP,ZDR), which is constructed by combining R(ZH), R(ZH, ZDR) and R(KDP), based on the statistical QPE error in the ZH-ZDR space, and was proven to outperform any single rainfall estimator in typical heavy rainfall events (e.g., mei-yu, typhoon rainbands and squall lines) in East China. However, the composite estimator is sometimes discontinuous owing to the hard thresholds for switching among different rainfall estimators. To overcome this drawback, (Huang et al., 2018a) proposed using a variational approach for QPE, which statistically combines the information provided by radar measurements (ZH and ΦDP) and applies spatial continuity of rainfall in a unified framework. In this method, the R-KDP relationship, tuned using DSD observations in South China, is used for the construction of the forward operator; the tuned R(ZH) is used as the a priori, with its error covariance matrix statistically determined, which can help to reduce the effect of measurement errors in Φ DP. It is found that the variational approach produces better rainfall estimation than the traditional rainfall estimator R(KDP) or composite algorithm in multiple rainfall cases, showing higher correlation coefficients and lower normalized absolute errors (Fig. 2).
Figure2. Hourly rainfall comparisons at rain gauge sites for (a) the variational approach of (Huang et al., 2018a) with R(ZH) as the a priori and (b) the conventional KDP-based approach. The places where the rain gauges were deployed are shown as circles, wherein the size of the circles represents the correlation coefficient between the time series of the radar-derived accumulated rainfalls (AR) and the time series of the gauge-derived AR, and the color represents the normalized absolute error (NE) between them. The NE is defined as $\rm NE=\frac1\rm N\sum_i=1^N|R_\rm e(i)-R_\rm g(i)|/\overline{R_\rm g}$, where N is the total sampling number at each gauge site, Rg(Re) is the hourly rainfall from gauge measurements (radar estimation), and $\overline{R_\rm g}$ is the corresponding mean value. [Reprinted from (Huang et al., 2018a). ? American Meteorological Society. Used with permission.]



2
3.3. Retrieval of DSD
--> DSD is a fundamental characteristic of rain microphysics, which can be used to represent all rain physical parameters. Since a DSD contains numerous unknowns, the exponential distribution (Blanchard, 1953; Seliga and Bringi, 1978) and the gamma distribution (Ulbrich, 1983) have been proposed to approximate natural DSDs. It is well-known that retrieving DSDs from polarimetric data using a three-parameter gamma distribution is ill-posed (e.g., Huang et al., 2019). An extra physical constraint for the gamma distribution model, e.g., the statistical relation between the parameters μ and $\Lambda$ (the slope term) or a fixed value for μ, helps to improve the accuracy of DSD retrieval from ZH and ZDR (Seliga and Bringi, 1978; Zhang et al., 2001). As revealed by the result in (Huang et al., 2019) (Fig. 3), when the μ and $\Lambda$ relation is utilized in the retrieval, the radar-derived precipitation parameters (R, mass-weighted mean diameter Dm, and total number concentration Nt) are generally consistent with the measurements from disdrometers; when a three-parameter gamma distribution is used as the model for DSD retrieval, the correlation coefficients of R, Dm and Nt between the estimates and measurements decrease to 0.58, 0.48 and 0.04 (not shown), respectively. Since DSDs can vary with different climate regions and different geophysical locations, μ-$\Lambda$ relationships need to be refined in different locations of China. (Li et al., 2015), (Wen et al., 2018) and (Liu et al., 2018) have constructed and applied μ-$\Lambda$ relationships to DSD retrievals in Northeast, East and South China, respectively.
Figure3. Comparisons of radar-retrieved (a) R, (b) Dm, and (c) total number concentration Nt with those calculated from 2DVD data (black lines). The red dots and green circles represent the results from error minimization analysis (EMA)-based retrieval using a constrained-gamma distribution (CG) and the three-parameter gamma distribution (GM) with KDP measurements included. [Reprinted from (Huang et al., 2019). ? American Meteorological Society. Used with permission.]


To reduce the impact of measurement errors on the retrievals, (Huang, 2018) proposed using variational analysis for DSD retrieval. In this optimization, the attenuation effects are considered in observation operators, which help to avoid the error propagation from attenuation correction to DSD retrieval. The measurement errors are also mitigated by an azimuthal Kalman filter and a radial B-spline filter. Verification using C- and S-band radar observations shows satisfactory performance of the variational approach.

4. Precipitation microphysics and processes in China revealed by dual-pol radar
Precipitation microphysics is one of the key factors determining the behaviors of convective systems owing to the nonlinear interactions between microphysics and dynamics through latent heat release or absorption in microphysical processes. The 3D predominant hydrometeor and DSD distribution can be obtained by using HC algorithms and DSD retrieval methods, as described in section 3, to portray microphysical characteristics with high temporal and spatial resolution and infer the dominant microphysical process (Kumjian and Ryzhkov, 2010, Kumjian and Ryzhkov, 2012; Kumjian and Prat, 2014; Barnes and Houze, 2016; Wang et al., 2018b).
DSDs of convective rain in different climatic regimes exhibit two clusters within the framework of generalized intercept and median volume diameter, known as the so-called "maritime-like" and "continental-like" DSD characteristics (Bringi et al., 2003). In general, maritime-type convective precipitation possesses a higher number concentration of small/medium sized raindrops than continental-type convective precipitation. Recently, (Dolan et al., 2018) revealed that the variation of environmental conditions, in addition to geographic locations, also affects the microphysical properties of convective systems, based on twelve sets of disdrometer observations across three latitudinal bands across the globe. Over the past five years, precipitation microphysics has been investigated extensively in different convective systems in China by combining dual-pol radar observations and disdrometer observations or numerical model results.

2
4.1. Polarimetric signatures and precipitation microphysics in convective systems
--> The microphysical characteristics of MCSs in the US have been well documented, especially for supercell thunderstorms. Typical polarimetric signatures in supercells were reviewed and summarized by (Kumjian and Ryzhkov, 2008), including the ZDR arc, KDP foot, ZDR column/ring, KDP column, ρhv ring, large hail signature, and so on.
Similar signatures have also been observed in different convective systems in China, such as hailstorms (Chen et al., 2014), supercell storms (Zhang et al., 2017a), and MCSs (Zhang et al., 2017b). A supercell case that occurred in Qingyuan was studied by (Zhang et al., 2017a) using data collected by an S-band dual-pol radar. In that case, a similar large hail signature was presented as large ZH, reduced ρhv, and near-zero ZDR. (Wang et al., 2018a) presented a ZDR column within updrafts and ZDR arc near the forward-flank downdraft from a supercell (Fig. 4). (Zhang et al., 2018) used an X-band dual-pol radar to identify the tornadic debris signature from the Foshan tornado within an outer rainband of Typhoon Mujigae.
Figure4. PPI of Zhuhai S-band dual-pol radar at 0.5° elevation at 0909 UTC 20 April 2015: (a) ZH; (b) ZDR; (c) ρhv; (d) KDP. [Reprinted from (Wang et al., 2018a)]


Based on dual-pol radar observations from OPACC, SCMREX and TIPEX-III, microphysical processes in different regions of China have been investigated (e.g., Gao et al., 2016; Luo et al., 2017; Wen et al., 2017; Wang et al., 2019). (Wen et al., 2017) investigated the variations of microphysical characteristics within the convective region during the formative, intensifying, and mature stages of a subtropical squall line in summer using the NJU-CPOL observations during OPACC in eastern China. The radar-derived DSD in the convection region of a squall line evolved from more continental-like to maritime-like characteristics when the system developed from the formative stage to the mature stage (Fig. 5), which is different from previous studies where the DSD characteristics of a convective line mostly depend on the geographical location rather than within the life cycle of a squall line (Petersen and Rutledge, 2001). The dual-pol radar-derived liquid water content below the freezing level in the convective region was three times higher than the ice water content above the freezing level, indicating the dominance of the warm rain process within this squall line. (Luo et al., 2017) showed RHI scans of two MCSs over Guangdong collected by a C-band dual-pol radar in the SCMREX field campaign (Fig. 6). ZDR and KDP columns were identified within convective regions, indicating vigorous updrafts. The increases of ZH, ZDR and KDP toward the ground provided clear signatures of rainwater growth through warm-rain processes. Raindrop breakup was also noticed below the altitude of 2 km, which was characterized as KDP and ZDR decreasing toward the ground. Contrary to MCSs in East and South China where warm-rain processes are dominant owing to the influence of the East Asian summer monsoon, MCSs over the Tibetan Plateau develop much deeper with more distinct ice processes (Mei et al., 2018).
Figure5. (a-c) The CAPPI of Z km above ground level from the NJU C-POL radar at 2157 LST (formative stage), 2217 LST (developing stage), and 2251 LST (mature stage), respectively, 30 July 2014. The convective region is enclosed by the black solid lines. (d, e) Frequency distribution of Dm and lgNw retrieved using the constrained-gamma model from the C-POL radar data of convective regions, at 1-km elevation only, for the three stages (a-c). The outermost gray line represents 5% contours. The mean Nw and Dm values for all convective regions are represented by the black plus signs. The two gray rectangles correspond to the maritime and continental convective clusters reported by (Bringi et al., 2003). In (f), the square signs represent mean values for the convective center (CC), and the triangle signs represent those for the convective edge (CE) combined. [Reprinted from (Wen et al., 2017).]


Figure6. Vertical cross section at about 1752 Local Standard Time (LST) 8 May 2014 of the Heshan C-POL measurements of (a) reflectivity ZH, (b) differential reflectivity ZDR, (c) specific differential phase KDP, and (d) correlation coefficient ρhv. The black dashed and solid lines represent the 0°C level (4.6 km) and -15°C level (7 km), respectively, according to sounding data. (e-h) As in (a-d), respectively, but at about 1604 LST 22 May 2014; the 0°C and -15°C levels are 5.2 and 8.1 km, respectively. [Reprinted from (Luo et al., 2017). ? American Meteorological Society. Used with permission.]



2
4.2. Precipitation microphysics of landfall typhoons
--> In China, the DSDs of landfalling TCs observed by 2-dimensional video disdrometers (2DVDs) mainly consist of very small drops and high number concentrations——more like maritime-type convection than those of TCs in Taiwan (Chang et al., 2009; Wen et al., 2018). The DSDs in the inner rainband of Typhoon Matmo (2014) observed by a 2DVD and retrieved from dual-pol radar measurements also show the characteristics of typical maritime-type convection (Fig.7) (Wang et al., 2016b; Wen et al., 2018). It is also found that warm-rain processes were predominant within the convective region of the inner rainband of Typhoon Matmo (2014).
Figure7. The (a) reflectivity and (b) differential reflectivity at 0.5° elevation observed by Lishui Radar (LSRD) at 1100 UTC 24 September 2014. (c) Frequency of occurrences (color shaded) of Dm (units: mm) and logarithmic Nw (units: mm-1 m-3) of the retrieved DSDs from LSRD. The gray crosses represent the Dm and Nw values calculated from 2DVD data. The dashed line indicates the rainfall rate of 10 mm h-1. The two outlined solid/dashed squares represent the maritime/continental types of convective systems. The gray square, black square, and black dot indicate the mean value of Dm and Nw from 2DVD, LSRD, and the study of (Chang et al., 2009) for rainfall rates over 10 mm h-1. [Reprinted from (Wang et al., 2016b).]


Based on HC, (Wang et al., 2018b) further found that heavy rainfall tends to locate in the updraft and downdraft regions affected by graupel. Within the updraft region, heavy rainfall was generally produced by the warm-rain processes of auto-conversion, accretion, and coalescence from 5 km to 0.5 km in altitude, while melting of graupel particles dominated in the downdraft region.
(Wu et al., 2018b) examined the microphysics of convective cells in an outer rainband of Typhoon Nida (2016) using an S-band dual-pol radar. Combining ZH, ZDR and KDP information suggested a layered microphysical structure with riming near the -5°C level, aggregation around the -15°C level, and deposition almost everywhere above the freezing level. Ice processes dominated the precipitation in outer rainbands, being characterized by a much higher ZH and ZDR (Fig. 8).
Figure8. Median profiles of (a) reflectivity, (b) ZDR, (c) ice water content, and (d) liquid water content at the convective center in the inner rainband (blue lines) and the mature stage of the outer rainband (red lines) of Typhoon Nida (2016). [Reprinted from (Wu et al., 2018b). ? American Meteorological Society. Used with permission.]



5. Use of dual-pol radar data for nowcasting and NWPs
Dual-pol radar observations have been used to validate numerical model outputs using forward operators (Vivekanandan et al., 1991; Ryzhkov et al., 2011), as well as parameterize and initialize numerical models by radar data assimilation (e.g., Jung et al., 2008; Posselt et al., 2015). As more dual-pol data have become available from several field experiments in China, these data have provided unique opportunities to validate numerical model results (Gao et al., 2016; Wang et al., 2016a; Wen, 2017). (Wang et al., 2016a) developed a simulator to transfer model outputs into S-band dual-polar radar parameters based on Rayleigh-Gans scattering theory. The simulator can calculate polarimetric variables by using cloud mixing ratios and number concentrations from microphysics schemes and the axis ratio, relative dielectric constant and canting angles of particles. It can reproduce typical polarimetric radar signatures of a mature 2D idealized squall line, including hail with high ZH and low ZDR, and a ZDR column in the convective updraft region. (Gao et al., 2016) found that ZDR was higher (lower) in the convective (stratiform) regions compared to observations for a plateau summertime rainfall event, and thus identified bias in modeled hydrometeor types. Recently, (Wen, 2017) simulated a squall line system using the ARPS model with three bulk schemes. They found that microphysics schemes with different moments had large impacts on the simulated polarimetric radar variables, and the three-moment scheme reproduced the best characteristics of the simulated squall line.
These limited number of studies focused mainly on summertime convective systems and compared simulations to dual-pol radar observations. As more dual-pol CINRAD-98Ds become available in the future, there will be more opportunities in China for model comparisons, forecast evaluations, and radar data assimilation.

6. Summary and future outlook
This paper has reviewed recent advancements in dual-pol radar in China, including radar data quality control, microphysical retrieval algorithms and microphysical characteristics of summer precipitation systems, to guide efforts to further utilize dual-pol radars in the understanding, warning and forecasting severe weather. Studies indicate the importance of environmental conditions rather than latitude-dependence in determining the dominant microphysical processes and DSDs suggested in (Dolan et al., 2018). The limited number of studies in China that have sought to use dual-pol radar data to validate the microphysical parameterization and initialization of numerical models and assimilate dual-pol data into numerical models have been summarized. This line of work remains a ongoing research topic, and thus the research community will face many technical challenges in the foreseeable future.
Microphysical parameterization is very important for numerical weather models to accurately simulate precipitation systems. To improve parameterization schemes, efforts are still needed to gain knowledge of precipitation microphysics using different types of instruments, especially dual-pol radars, in the future. In addition, developing techniques for initialization and assimilating dual-pol radar measurements is urgent for the numerical forecasting of mesoscale weather systems.
Multi-frequency and phased arrays are two emerging technologies that can provide additional microphysical information and reduce the radar sampling time. Phased array dual-pol radars can measure precipitation systems more rapidly, which also gains more dynamical and microphysical information. However, the H- and V-polarized beam matching off the broadside affecting the dual-pol measurement remains a challenge for radar engineers. In contrast, the system of multi-frequency dual-pol radars is based on the variations of the Mie scattering effect with respect to the frequency of electromagnetic waves, which also provides information on the sizes of hydrometeors. It is anticipated that these two types of technologies will mature in the next decade and make significant impacts in the future of radar meteorology.

相关话题/Recent Progress Polarization