Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (54)

Search Parameters:
Keywords = raindrop size distribution (DSD)

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
17 pages, 3131 KiB  
Article
Microphysical Characteristics of Precipitation for Four Types of Typical Weather Systems on Hainan Island
by Wupeng Xiao, Yun Zhang, Hepeng Zheng, Zuhang Wu, Yanqiong Xie and Yanbin Huang
Remote Sens. 2024, 16(22), 4144; https://doi.org/10.3390/rs16224144 - 6 Nov 2024
Viewed by 493
Abstract
The microphysical characteristics of precipitation and their differences among four typical weather systems over Hainan Island were investigated via multi-source observations from 2019 to 2023. We find that the cold fronts (CFs) have the greatest concentration of small raindrops, with a more substantial [...] Read more.
The microphysical characteristics of precipitation and their differences among four typical weather systems over Hainan Island were investigated via multi-source observations from 2019 to 2023. We find that the cold fronts (CFs) have the greatest concentration of small raindrops, with a more substantial raindrop condensation process. The subtropical highs (SHs), with primarily deep convection and more prominent evaporation at low levels, lead to greater medium-to-large raindrops (diameters > 1 mm). Tropical cyclones (TCs) are characterized mainly by raindrop condensation and breakup, resulting in high concentrations of small raindrops and low concentrations of large raindrops. The trough of low pressures (TLPs) produces the lowest concentration of small raindrops because of evaporation processes. The convective clusters of the SHs are between maritime-like and continental-like convective clusters, and those of the other three types of weather systems are closer to maritime-like convective clusters. The relationships between the shape parameter (μ) and the slope parameter (Λ), as well as between the reflectivity factors (Z) and the rain rates (R), were established for the four weather systems. These results could improve the accuracy of radar quantitative precipitation estimation and the microphysical parameterizations of numerical models for Hainan Island. Full article
(This article belongs to the Special Issue Synergetic Remote Sensing of Clouds and Precipitation II)
Show Figures

Figure 1

Figure 1
<p>Distribution of OTTs and AWSs on Hainan Island.</p>
Full article ">Figure 2
<p>Circulation conditions of the four types of typical weather systems on Hainan Island. Composite of the 500 hPa geopotential height (black contours), the 850 hPa wind vector (blue arrows), and the 700 hPa specific humidity (g kg<sup>−1</sup>, shading). (<b>a</b>) CFs—the blue curve approximates the position of the cold front. (<b>b</b>) SHs—the black bold contours represent the 5880 gpm lines. (<b>c</b>) TCs—the yellow typhoon marker indicates the location of the center of the tropical cyclone Lionrock. (<b>d</b>) TLPs—the brown curve represents the location of the trough of low pressure.</p>
Full article ">Figure 3
<p>Average raindrop size distributions of the four types of weather systems, where the blue, red, purple, and green solid lines represent CFs, SHs, TCs, and TLPs, respectively.</p>
Full article ">Figure 4
<p>Average raindrop size distributions of the four types of weather systems in different rainfall rates (<b>a</b>–<b>d</b>) represent <math display="inline"><semantics> <mi>R</mi> </semantics></math> ≤ 10 mm h<sup>−1</sup>, 10 &lt; <math display="inline"><semantics> <mi>R</mi> </semantics></math> ≤ 20 mm h<sup>−1</sup>, 20 &lt;<math display="inline"><semantics> <mi>R</mi> </semantics></math> ≤ 50 mm h<sup>−1</sup>, and <math display="inline"><semantics> <mi>R</mi> </semantics></math> &gt; 50 mm h<sup>−1</sup>, where the blue, red, purple, and green solid lines represent CFs, SHs, TCs, and TLPs, respectively.</p>
Full article ">Figure 5
<p>Relative contributions of raindrops to (<b>a</b>) the rainfall rate <math display="inline"><semantics> <mi>R</mi> </semantics></math> (<b>b</b>) the total raindrop concentration <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </mrow> </semantics></math>, and (<b>c</b>) the reflectivity factor <span class="html-italic">Z</span> in different diameter bins, where the blue, red, purple, and green regions represent CFs, SHs, TCs, and TLPs, respectively.</p>
Full article ">Figure 6
<p>Distribution of <math display="inline"><semantics> <mrow> <msub> <mi>D</mi> <mi>m</mi> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>l</mi> <mi>o</mi> <msub> <mi>g</mi> <mrow> <mn>10</mn> </mrow> </msub> <mo stretchy="false">(</mo> <msub> <mi>N</mi> <mi mathvariant="normal">w</mi> </msub> <mo stretchy="false">)</mo> </mrow> </semantics></math> for convective precipitation and stratiform precipitation for the four types of weather systems. The black box represents the region of maritime and continental convective precipitation, as defined by [<a href="#B15-remotesensing-16-04144" class="html-bibr">15</a>], and the thick black dashed line represents the stratiform precipitation fitting line. The thin black dashed lines represent the contours of the rainfall rate. The dark gray crosses and light gray dots represent convective and stratiform precipitation, respectively. The circle, square and rhombus symbols in (<b>a</b>) indicate the Meiyu front in Central China [<a href="#B43-remotesensing-16-04144" class="html-bibr">43</a>] and East China [<a href="#B44-remotesensing-16-04144" class="html-bibr">44</a>,<a href="#B45-remotesensing-16-04144" class="html-bibr">45</a>]. Red and blue shading represent convective precipitation and stratiform precipitation, respectively. (<b>b</b>) Western (WWP), southern (SWP), and northern (NWP) of the western Pacific subtropical high [<a href="#B21-remotesensing-16-04144" class="html-bibr">21</a>]. (<b>c</b>) The circle, square and rhombus denote the convective precipitation of tropical cyclones that made landfall in East China and South China [<a href="#B42-remotesensing-16-04144" class="html-bibr">42</a>], Taiwan [<a href="#B26-remotesensing-16-04144" class="html-bibr">26</a>], and Hainan [<a href="#B27-remotesensing-16-04144" class="html-bibr">27</a>], and (<b>d</b>) The circle, square and rhombus indicate the pre-, mid-, and post-monsoon periods in the South China Sea, respectively [<a href="#B22-remotesensing-16-04144" class="html-bibr">22</a>].</p>
Full article ">Figure 7
<p>Scatterplot density distributions of <math display="inline"><semantics> <mrow> <msub> <mi>D</mi> <mi>m</mi> </msub> </mrow> </semantics></math> (<b>a</b>–<b>d</b>) and <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </mrow> </semantics></math> (<b>e</b>–<b>h</b>) with <span class="html-italic">R</span>. The red curves are the least-squares-fitted <math display="inline"><semantics> <mrow> <msub> <mi>D</mi> <mi>m</mi> </msub> </mrow> </semantics></math>-<span class="html-italic">R</span> and <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </mrow> </semantics></math>-<span class="html-italic">R</span> relationships, and the gray dashed line represents the 10 mm h<sup>−1</sup> contour.</p>
Full article ">Figure 8
<p>Fitted relationships for <math display="inline"><semantics> <mrow> <msub> <mi>D</mi> <mi>m</mi> </msub> </mrow> </semantics></math> (<b>a</b>) and <math display="inline"><semantics> <mrow> <msub> <mi>N</mi> <mi>t</mi> </msub> </mrow> </semantics></math> (<b>b</b>) of the four types of weather systems with rainfall rate <span class="html-italic">R</span>, with the gray dashed line representing the 10 mm h<sup>−1</sup> contour.</p>
Full article ">Figure 9
<p>(<b>a</b>) Boxplot of CAPE. The red solid line represents the median, the blue dashed line represents the mean, and the red dots represent the anomalies; (<b>b</b>) Distribution of the mean probability density of the TBB derived from FY-4A, where the dashed lines represent the dividing lines of stratiform and shallow convection (−10 °C), moderate convection (−32 °C), deep convection (−60 °C), and extreme convection (−75 °C). The blue, red, purple, and green lines represent cold fronts, subtropical highs, tropical cyclones, and low-pressure troughs, respectively. (<b>c</b>) Boxplots of the LCL, 0 °C level height, and CTH.</p>
Full article ">Figure 10
<p>Vertical profiles of the (<b>a</b>) temperature, (<b>b</b>) wind speed, (<b>c</b>) relative humidity, and (<b>d</b>) specific humidity for the four types of weather systems on Hainan Island.</p>
Full article ">Figure 11
<p><span class="html-italic">μ-</span>Λ and <span class="html-italic">Z-R</span> relationships for the four types of weather systems. (<b>a</b>–<b>d</b>) show the probability density distributions of the <span class="html-italic">μ-</span>Λ relationship and the quadratic polynomial fitting curves. The color bars on the right side represent the densities of the points in the scatterplot, where the data with precipitation rates of <span class="html-italic">R</span> &lt; 5 mm h<sup>−1</sup> are excluded, and the fitting curves are shown in (<b>e</b>). (<b>f</b>) shows the <span class="html-italic">Z-R</span> relationship for the corresponding system and the WSR-88D empirical relationship [<a href="#B48-remotesensing-16-04144" class="html-bibr">48</a>].</p>
Full article ">
29 pages, 9650 KiB  
Article
Seasonal Variations in the Rainfall Kinetic Energy Estimation and the Dual-Polarization Radar Quantitative Precipitation Estimation Under Different Rainfall Types in the Tianshan Mountains, China
by Yong Zeng, Lianmei Yang, Zepeng Tong, Yufei Jiang, Abuduwaili Abulikemu, Xinyu Lu and Xiaomeng Li
Remote Sens. 2024, 16(20), 3859; https://doi.org/10.3390/rs16203859 - 17 Oct 2024
Viewed by 655
Abstract
Raindrop size distribution (DSD) has an essential effect on rainfall kinetic energy estimation (RKEE) and dual-polarization radar quantitative precipitation estimation (QPE); DSD is a key factor for establishing a dual-polarization radar QPE scheme and RKEE scheme, particularly in mountainous areas. To improve the [...] Read more.
Raindrop size distribution (DSD) has an essential effect on rainfall kinetic energy estimation (RKEE) and dual-polarization radar quantitative precipitation estimation (QPE); DSD is a key factor for establishing a dual-polarization radar QPE scheme and RKEE scheme, particularly in mountainous areas. To improve the understanding of seasonal DSD-based RKEE, dual-polarization radar QPE, and the impact of rainfall types and classification methods, we investigated RKEE schemes and dual-polarimetric radar QPE algorithms across seasons and rainfall types based on two classic classification methods (BR09 and BR03) and DSD data from a disdrometer in the Tianshan Mountains during 2020–2022. Two RKEE schemes were established: the rainfall kinetic energy flux–rain rate (KEtimeR) and the rainfall kinetic energy content–mass-weighted mean diameter (KEmmDm). Both showed seasonal variation, whether it was stratiform rainfall or convective rainfall, under BR03 and BR09. Both schemes had excellent performance, especially the KEmmDm relationship across seasons and rainfall types. In addition, four QPE schemes for dual-polarimetric radar—R(Kdp), R(Zh), R(Kdp,Zdr), and R(Zh,Zdr)—were established, and exhibited characteristics that varied with season and rainfall type. Overall, the performance of the single-parameter algorithms was inferior to that of the double-parameter algorithms, and the performance of the R(Zh) algorithm was inferior to that of the R(Kdp) algorithm. The results of this study show that it is necessary to consider different rainfall types and seasons, as well as classification methods of rainfall types, when applying RKEE and dual-polarization radar QPE. In this process, choosing a suitable estimator—KEtime(R), KEmm(Dm), R(Kdp), R(Zh), R(Kdp,Zdr), or R(Zh,Zdr)—is key to improving the accuracy of estimating the rainfall KE and R. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Topography (m) and location of the Tianshan Mountains, and (<b>b</b>) locations of Zhaosu (red dot) and Xinyuan (black dot; Zeng et al. [<a href="#B55-remotesensing-16-03859" class="html-bibr">55</a>]).</p>
Full article ">Figure 1 Cont.
<p>(<b>a</b>) Topography (m) and location of the Tianshan Mountains, and (<b>b</b>) locations of Zhaosu (red dot) and Xinyuan (black dot; Zeng et al. [<a href="#B55-remotesensing-16-03859" class="html-bibr">55</a>]).</p>
Full article ">Figure 2
<p>Seasonal variations in the distributions of (<b>a</b>) <span class="html-italic">KE<sub>time</sub></span> and (<b>b</b>) <span class="html-italic">KE<sub>mm</sub></span> at Zhaosu.</p>
Full article ">Figure 3
<p>Scatterplots of <span class="html-italic">KE<sub>time</sub></span> vs. <span class="html-italic">R</span> for the entire data and the fitted <span class="html-italic">KE<sub>time</sub></span>–<span class="html-italic">R</span> relationship across seasons at Zhaosu. Dashed lines represent the <span class="html-italic">KE<sub>time</sub></span>–<span class="html-italic">R</span> relationship reported by Zeng et al. [<a href="#B55-remotesensing-16-03859" class="html-bibr">55</a>], Seela et al. [<a href="#B83-remotesensing-16-03859" class="html-bibr">83</a>], and Wu et al. [<a href="#B36-remotesensing-16-03859" class="html-bibr">36</a>].</p>
Full article ">Figure 4
<p>Scatterplots of <span class="html-italic">KE<sub>mm</sub></span> vs. <span class="html-italic">D<sub>m</sub></span> for the entire data and the seasonal variation in fitted <span class="html-italic">KE<sub>mm</sub></span>–<span class="html-italic">D<sub>m</sub></span> at Zhaosu. Dashed lines represent the <span class="html-italic">KE<sub>mm</sub></span>–<span class="html-italic">D<sub>m</sub></span> relationship reported by Zeng et al. [<a href="#B55-remotesensing-16-03859" class="html-bibr">55</a>] and Seela et al. [<a href="#B83-remotesensing-16-03859" class="html-bibr">83</a>].</p>
Full article ">Figure 5
<p>Scatterplot of estimated <span class="html-italic">KE<sub>time</sub></span> from RKEE schemes versus <span class="html-italic">KE<sub>time</sub></span> calculated from DSD for (<b>a</b>) the entire data, (<b>b</b>) spring, (<b>c</b>) summer, and (<b>d</b>) fall at Zhaosu. Scatterplot of estimated <span class="html-italic">KE<sub>mm</sub></span> from RKEE schemes versus the <span class="html-italic">KE<sub>mm</sub></span> calculated from DSD for (<b>e</b>) the entire data, (<b>f</b>) spring, (<b>g</b>) summer, and (<b>h</b>) fall at Zhaosu in Tianshan Mountains. Black dashed lines represent the 1:1 relationship.</p>
Full article ">Figure 6
<p>Violin plots of seasonal variations in <span class="html-italic">KE<sub>time</sub></span> under (<b>a</b>) BR09_S, (<b>c</b>) BR09_C, (<b>e</b>) BR03_S, and (<b>g</b>) BR03_C, and violin plots of seasonal variations in <span class="html-italic">KE<sub>mm</sub></span> under (<b>b</b>) BR09_S, (<b>d</b>) BR09_C, (<b>f</b>) BR03_S, and (<b>h</b>) BR03_C at Zhaosu.</p>
Full article ">Figure 7
<p>Scatterplots of <span class="html-italic">KE<sub>time</sub></span> vs. <span class="html-italic">R</span> for the entire data and the seasonal variation of the fitted <span class="html-italic">KE<sub>time</sub></span>–<span class="html-italic">R</span> relationship at Zhaosu under (<b>a</b>) BR09_S, (<b>b</b>) BR09_C, (<b>c</b>) BR03_S, and (<b>d</b>) BR03_C.</p>
Full article ">Figure 8
<p>Scatterplot of estimated <span class="html-italic">KE<sub>time</sub></span> from RKEE schemes versus <span class="html-italic">KE<sub>time</sub></span> calculated from DSD for (<b>a</b>) the entire data, (<b>e</b>) spring, (<b>i</b>) summer, and (<b>m</b>) fall under BR09_S; those for (<b>b</b>) the entire data, (<b>f</b>) spring, (<b>j</b>) summer, and (<b>n</b>) fall under BR09_C; those for (<b>c</b>) the entire data, (<b>g</b>) spring, (<b>k</b>) summer, and (<b>o</b>) fall under BR03_S; and those for (<b>d</b>) the entire data, (<b>h</b>) spring, and (<b>l</b>) summer under BR03_C at Zhaosu. Black dashed lines represent the 1:1 relationship.</p>
Full article ">Figure 9
<p>Scatterplots of <span class="html-italic">KE<sub>mm</sub></span> vs. <span class="html-italic">D<sub>m</sub></span> for the entire data and the fitted <span class="html-italic">KE<sub>mm</sub></span>–<span class="html-italic">D<sub>m</sub></span> relationship across seasons at Zhaosu under (<b>a</b>) BR09_S, (<b>b</b>) BR09_C, (<b>c</b>) BR03_S, and (<b>d</b>) BR03_C.</p>
Full article ">Figure 10
<p>Scatterplot of estimated <span class="html-italic">KE<sub>mm</sub></span> from RKEE schemes versus <span class="html-italic">KE<sub>mm</sub></span> calculated from DSD for (<b>a</b>) the entire data, (<b>e</b>) spring, (<b>i</b>) summer, and (<b>m</b>) fall under BR09_S; for (<b>b</b>) the entire data, (<b>f</b>) spring, (<b>j</b>) summer, and (<b>n</b>) fall under BR09_C; for (<b>c</b>) the entire data, (<b>g</b>) spring, (<b>k</b>) summer, and (<b>o</b>) fall under BR03_S; and for (<b>d</b>) the entire data, (<b>h</b>) spring, and (<b>l</b>) summer under BR03_C at Zhaosu. Black dashed lines represent the 1:1 relationship.</p>
Full article ">Figure 11
<p>Seasonal variations in the distributions of (<b>a</b>) <span class="html-italic">Z<sub>h</sub></span>, (<b>b</b>) <span class="html-italic">Z<sub>dr</sub></span>, and (<b>c</b>) <span class="html-italic">K<sub>dp</sub></span> at Zhaosu.</p>
Full article ">Figure 12
<p>Scatterplot of estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>) algorithm for (<b>a</b>) the entire data, (<b>e</b>) spring, (<b>i</b>) summer, and (<b>m</b>) fall; estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>) algorithm for (<b>b</b>) the entire data, (<b>f</b>) spring, (<b>j</b>) summer, and (<b>n</b>) fall; estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm for (<b>c</b>) the entire data, (<b>g</b>) spring, (<b>k</b>) summer, and (<b>o</b>) fall; and estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm for (<b>d</b>) the entire data, (<b>h</b>) spring, (<b>l</b>) summer, and (<b>p</b>) fall versus calculated <span class="html-italic">R</span> according to Equation (4) at Zhaosu. Black dashed lines represent the 1:1 relationship.</p>
Full article ">Figure 13
<p>Seasonal variations in the distributions of <span class="html-italic">Z<sub>h</sub></span> under (<b>a</b>) BR09_S, (<b>d</b>) BR09_C, (<b>g</b>) BR03_S, and (<b>j</b>) BR03_C; those of <span class="html-italic">Z<sub>dr</sub></span> under (<b>b</b>) BR09_S, (<b>e</b>) BR09_C, (<b>h</b>) BR03_S, and (<b>k</b>) BR03_C; and those of <span class="html-italic">K<sub>dp</sub></span> under (<b>c</b>) BR09_S, (<b>f</b>) BR09_C, (<b>i</b>) BR03_S, and (<b>l</b>) BR03_C at Zhaosu.</p>
Full article ">Figure 14
<p>Scatterplot of estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>) algorithm during BR09_S for (<b>a</b>) the entire data, (<b>e</b>) spring, (<b>i</b>) summer, and (<b>m</b>) fall; that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>) algorithm during BR09_S for (<b>b</b>) the entire data, (<b>f</b>) spring, (<b>j</b>) summer, and (<b>n</b>) fall; that based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR09_S for (<b>c</b>) the entire data, (<b>g</b>) spring, (<b>k</b>) summer, and (<b>o</b>) fall; and that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR09_S for (<b>d</b>) the entire data, (<b>h</b>) spring, (<b>l</b>) summer, and (<b>p</b>) fall versus calculated <span class="html-italic">R</span> according to Equation (4) at Zhaosu. Black dashed lines represent the 1:1 relationship.</p>
Full article ">Figure 15
<p>Scatterplot of estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>) algorithm during BR09_C for (<b>a</b>) the entire data, (<b>e</b>) spring, (<b>i</b>) summer, and (<b>m</b>) fall; that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>) algorithm during BR09_C for (<b>b</b>) the entire data, (<b>f</b>) spring, (<b>j</b>) summer, and (<b>n</b>) fall; that based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR09_C for (<b>c</b>) the entire data, (<b>g</b>) spring, (<b>k</b>) summer, and (<b>o</b>) fall; and that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR09_C for (<b>d</b>) the entire data, (<b>h</b>) spring, (<b>l</b>) summer, and (<b>p</b>) fall versus calculated <span class="html-italic">R</span> according to Equation (4) at Zhaosu. Black dashed lines represent the 1:1 relationship.</p>
Full article ">Figure 16
<p>Scatterplot of estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>) algorithm during BR03_S for (<b>a</b>) the entire data, (<b>e</b>) spring, (<b>i</b>) summer, and (<b>m</b>) fall; that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>) algorithm during BR03_S for (<b>b</b>) the entire data, (<b>f</b>) spring, (<b>j</b>) summer, and (<b>n</b>) fall; that based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR03_S for (<b>c</b>) the entire data, (<b>g</b>) spring, (<b>k</b>) summer, and (<b>o</b>) fall; and that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR03_S for (<b>d</b>) the entire data, (<b>h</b>) spring, (<b>l</b>) summer, and (<b>p</b>) fall versus calculated <span class="html-italic">R</span> according to Equation (4) at Zhaosu. Black dashed lines represent the 1:1 relationship.</p>
Full article ">Figure 17
<p>Scatterplot of estimated <span class="html-italic">R</span> based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>) algorithm during BR03_C for (<b>a</b>) the entire data, (<b>e</b>) spring, and (<b>i</b>) summer; that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>) algorithm during BR03_C for (<b>b</b>) the entire data, (<b>f</b>) spring, and (<b>j</b>) summer; that based on the <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR03_C for (<b>c</b>) the entire data, (<b>g</b>) spring, and (<b>k</b>) summer; and that based on the <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) algorithm during BR03_C for (<b>d</b>) the entire data, (<b>h</b>) spring, and (<b>l</b>) summer versus calculated <span class="html-italic">R</span> according to Equation (4) at Zhaosu. Black dashed lines represent the 1:1 relationship.</p>
Full article ">
12 pages, 4993 KiB  
Article
Raindrop Size Distribution Characteristics of the Precipitation Process of 2216 Typhoon “Noru” in the Xisha Region
by Guozhang Wang, Lei Li, Chaoying Huang and Lili Zhang
Water 2024, 16(18), 2630; https://doi.org/10.3390/w16182630 - 16 Sep 2024
Viewed by 676
Abstract
This study focuses on the comparative analysis and research of the raindrop size distribution (DSD) in the outer rainband and inner rainband of Typhoon “Noru” in 2022, using the OTT-Parsivel raindrop spectrometer deployed on Yongxing Island, Sansha City. The results indicate that precipitation [...] Read more.
This study focuses on the comparative analysis and research of the raindrop size distribution (DSD) in the outer rainband and inner rainband of Typhoon “Noru” in 2022, using the OTT-Parsivel raindrop spectrometer deployed on Yongxing Island, Sansha City. The results indicate that precipitation intensity is lower when composed mainly of small and medium raindrops and increases with the presence of larger raindrops. Stronger precipitation is associated with a higher number of large raindrops. Due to the interaction of cold and warm air masses, the raindrop concentration is higher, and the raindrop diameters are larger compared to Typhoons “LEKIMA” and “RUMBIA”. The entire process predominantly consists of numerous small- to medium-sized raindrops, characteristic of a tropical typhoon. The precipitation in the inner and outer rainbands exhibits consistent types, characterized by a unimodal raindrop size distribution with a narrow spectral width, typical of stratiform-mixed cloud precipitation, where stratiform precipitation constitutes a significant portion. Strong echo reflectivity factors are often associated with higher raindrop number concentrations and larger particle sizes. The Z-R relationship of the precipitation shows a smaller coefficient but a consistent exponent compared to the standard. The calculated shape parameter slope relationship is Λ=0.02μ2+0.696μ+1.539, providing a reference for localizing the Z-R relationship in the South China Sea. Full article
(This article belongs to the Section Hydrology)
Show Figures

Figure 1

Figure 1
<p>During the influence of Typhoon Noru from 25 to 29 September 2022, the hourly rainfall on Yongxing Island in Xisha was compared using data from the rain gauge (RG) and the OTT disdrometer. The correlation coefficient (CC), standard deviation (SD), absolute bias (ab.bia (%)) and root mean square error (RMSE) were used as metrics for comparison.</p>
Full article ">Figure 2
<p>(<b>a</b>) shows the reflectivity map of Typhoon Noru’s DRP (UTC: 20220927-S054505-E071737), with the RMW representing the diameter of the typhoon eye, approximately 60 km. (<b>b</b>) displays the typhoon reflectivity profile along the blue line. It can be observed that rainfall within the RMW is relatively low, with the main precipitation bands concentrated in the inner rainbands (&lt;3 RMW) and the outer rainbands (&gt;3 RMW).</p>
Full article ">Figure 3
<p>The 500 hPa geopotential height fields derived from ERA5 reanalysis data, depicting snapshots at various times: (<b>a</b>) 20:00 on 26 September, (<b>b</b>) 08:00 on 27 September, (<b>c</b>) 20:00 on 27 September, and (<b>d</b>) 08:00 and 20:00 on 28 September.</p>
Full article ">Figure 4
<p>Temporal evolution of (<b>a</b>) total raindrop concentration <math display="inline"><semantics> <mrow> <mi>Nt</mi> <mo>/</mo> <mfenced> <mrow> <msup> <mrow> <mi>mm</mi> </mrow> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> <mo>·</mo> <msup> <mi mathvariant="normal">m</mi> <mrow> <mo>−</mo> <mn>3</mn> </mrow> </msup> </mrow> </mfenced> </mrow> </semantics></math>, (<b>b</b>) rainfall intensity <math display="inline"><semantics> <mrow> <mi mathvariant="normal">R</mi> <mo>/</mo> <mo stretchy="false">(</mo> <mi>mm</mi> <mo>·</mo> </mrow> </semantics></math>h<sup>−1</sup>), (<b>c</b>) mass-weighted mean diameter <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">D</mi> <mi mathvariant="normal">m</mi> </msub> <mo>/</mo> <mfenced> <mrow> <mi>mm</mi> </mrow> </mfenced> </mrow> </semantics></math>, and (<b>d</b>) radar reflectivity factor Z (dBZ) during Typhoon “Noru” at Yongxing station from 25 to 29 September 2022. S1 indicates the phase when the typhoon center is approaching Yongxing Island but is more than 180 km away, representing the outer rainband (&gt;3 RMW); S2 indicates the phase when the typhoon center is less than 180 km away from Yongxing Island, representing the inner rainband (&lt;3 RMW); S3 indicates the phase when the typhoon center is moving away from Yongxing Island, more than 180 km away, representing the outer rainband (&lt;3 RMW). The red lines in panel (<b>a</b>) mark the phase divisions, with red numbers indicating their corresponding times, where 26.2030 means 20:30 on the 26th and 27.0933 means 09:33 on the 27th.</p>
Full article ">Figure 5
<p>Temperature advection at 850 hPa at 08:00 UTC on 27 September 2022. The Chinese in the picture represents the place name where the equipment is located, and its position in the picture represents the location where the equipment is installed.</p>
Full article ">Figure 6
<p>(<b>a</b>) Mean raindrop spectra distribution for S1, S2, and S3 (solid lines) and their Gamma function fits (dashed lines); (<b>b</b>) scatter plot of generalized raindrop number concentration <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">N</mi> <mi mathvariant="normal">w</mi> </msub> </mrow> </semantics></math> versus mass-weighted mean diameter <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="normal">D</mi> <mi mathvariant="normal">m</mi> </msub> </mrow> </semantics></math>, with the red solid line representing the regression line for stratiform rain by Bringi et al. [<a href="#B29-water-16-02630" class="html-bibr">29</a>].</p>
Full article ">Figure 7
<p>(<b>a</b>) Z-R relationship calculated from raindrop spectrum data during Typhoon “Noru” (The red line represents the Z-R relationship of Typhoon “Noru”, the green line is the classical Z-R relationship, cited from Fulton, and the green and blue dashed lines are cited from Zhang et al. [<a href="#B9-water-16-02630" class="html-bibr">9</a>]) and (<b>b</b>) <math display="inline"><semantics> <mi>μ</mi> </semantics></math>-<math display="inline"><semantics> <mi>Λ</mi> </semantics></math> relationship (The red line represents the results obtained from the study of Typhoon Oulu, the deep blue represents the results obtained by Zhang et al. [<a href="#B30-water-16-02630" class="html-bibr">30</a>] in statistical analysis of raindrop spectra that conform to gamma distribution, the green represents the results obtained by Chang et al. [<a href="#B31-water-16-02630" class="html-bibr">31</a>] in statistical analysis of typhoons in the western Pacific, and the sky blue represents the results obtained by Chen et al. [<a href="#B24-water-16-02630" class="html-bibr">24</a>] in analyzing the rainfall process of a single typhoon MORAKOT. From the results, it can be seen that the two studies on individual typhoons have similar results).</p>
Full article ">
17 pages, 3218 KiB  
Article
Raindrop Size Distribution Characteristics for Typhoons over the Coast in Eastern China
by Dongdong Wang, Sheng Chen, Yang Kong, Xiaoli Gu, Xiaoyu Li, Xuejing Nan, Sujia Yue and Huayu Shen
Atmosphere 2024, 15(8), 951; https://doi.org/10.3390/atmos15080951 - 9 Aug 2024
Viewed by 689
Abstract
This study investigates the characteristics of the raindrop size distribution (DSD) for five typhoons that made landfall or passed by Zhejiang on the eastern coast of China, from 2019 to 2022. Additionally, it examines the raindrop shape–slope (µ-Λ) relationship, as well as the [...] Read more.
This study investigates the characteristics of the raindrop size distribution (DSD) for five typhoons that made landfall or passed by Zhejiang on the eastern coast of China, from 2019 to 2022. Additionally, it examines the raindrop shape–slope (µ-Λ) relationship, as well as the local Z-R relationship for these typhoons. The DSD datasets were collected by the DSG1 disdrometer located in Ningbo, Zhejiang Province. Based on rainfall rate (R), the DSD can be categorized into convective and stratiform rainfall types. Some rainfall parameters can also be derived from the DSDs to further analyze the specific characteristics of rainfall. The histograms of the generalized intercept parameter (log10Nw) exhibit negative skewness in both convective and stratiform rainfall, whereas the histograms of the mass-weighted mean diameter (Dm) of raindrops display positive skewness. During typhoon periods on the eastern coast of China, the DSD characteristic was composed of a lower number concentration of small and midsize raindrops (3.42 for log10Nw, 1.43 mm for Dm in the whole dataset) as compared to Jiangsu in eastern China, Tokyo, in Japan, Miryang, in South Korea, and Thiruvananthapuram in south India, respectively. At the same time, the scatter plots of Dm and log10Nw indicate that the convective rain during typhoon periods exhibits characteristics that are intermediate between “maritime-like” and “continental-like” clusters. Additionally, the raindrop spectra of convective rainfall and midsize raindrops in stratiform rainfall are well-represented by a three-parameter gamma distribution. The µ-Λ relation in this region is similar to Taiwan and Fujian, located along the southeastern coast of China. The Z-R relationship for eastern coastal China during typhoons based on filtered disdrometer data is Z = 175.04R1.53. These results could offer deeper insights into the microphysical characteristics of different rainfall types along the eastern coast of China and potentially improve the accuracy of precipitation estimates from weather radar observations. Full article
(This article belongs to the Special Issue Tropical Cyclones: Observations and Prediction)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Five typhoons that affect the eastern coastal areas of China, along with their differently colored trajectories. (<b>b</b>) Geographical location of the eastern coast of China corresponds to the red box in <a href="#atmosphere-15-00951-f001" class="html-fig">Figure 1</a>a. The red-colored solid circle indicates the location of the DSG1. The white polygon represents the Ningbo area of Zhejiang Province and the shading indicates topographic elevation levels (the data are derived from ETOPO 2022).</p>
Full article ">Figure 2
<p>The time series (local standard time; LST) of the DSDs observed on the eastern coast of China during the five typhoons. The shading indicates the logarithmic number concentration (mm<sup>−1</sup> m<sup>−3</sup>) of the DSDs. The left <span class="html-italic">y</span>-axis denotes the diameter (D, mm) of raindrops and the right <span class="html-italic">y</span>-axis represents the rainfall rate (R, mm h<sup>−1</sup>) derived from the disdromete, illustrated by a pink curve.</p>
Full article ">Figure 3
<p>Histogram panels for D<sub>m</sub> and log<sub>10</sub>N<sub>w</sub>, with D<sub>m</sub> depicted in gray and log<sub>10</sub>N<sub>w</sub> in black. Mean values, standard deviation (SD), and skewness (SK) are also shown in the respective panel. (<b>a</b>) For the whole, (<b>b</b>) for the stratiform, and (<b>c</b>) for the convective.</p>
Full article ">Figure 4
<p>Scatter distributions of log<sub>10</sub>N<sub>w</sub>-D<sub>m</sub> for convective and stratiform rainfall types during the five typhoons over the eastern coast of China, with convective rainfall depicted in orange and stratiform in sky blue. The two outlined rectangles represent the “maritime-like” and “continental-like” clusters identified by Bringi et al. [<a href="#B24-atmosphere-15-00951" class="html-bibr">24</a>]. The solid circle illustrates the mean of log<sub>10</sub>N<sub>w</sub>-D<sub>m</sub> for the convective rainfall of five typhoons on the eastern coast of China. Mean of convective log<sub>10</sub>N<sub>w</sub>-D<sub>m</sub> for the typhoon rainfall in the South China Sea (cross symbol) from Zheng et al. [<a href="#B8-atmosphere-15-00951" class="html-bibr">8</a>], Guangdong (solid triangle) from Feng et al. [<a href="#B27-atmosphere-15-00951" class="html-bibr">27</a>], Jiangsu (solid square) from Wen et al. [<a href="#B13-atmosphere-15-00951" class="html-bibr">13</a>], and northern Taiwan (plus sigh) from Chang et al. [<a href="#B23-atmosphere-15-00951" class="html-bibr">23</a>] are given as well, respectively.</p>
Full article ">Figure 5
<p>Relationship between log<sub>10</sub>N<sub>w</sub> (<b>a</b>,<b>b</b>) and D<sub>m</sub> (<b>c</b>,<b>d</b>) with R of the two rainfall types. The red solid lines denote their fitting curves. (<b>a</b>,<b>c</b>) for stratiform rainfall (SR) and (<b>b</b>,<b>d</b>) for convective rainfall (CR).</p>
Full article ">Figure 6
<p>The average raindrop spectral distribution and gamma fitting distribution of the two types of rainfall during the five typhoons.</p>
Full article ">Figure 7
<p>The µ-Λ relationship for convective rainfall during typhoons with the number of raindrops exceeding 1000. The black, blue, red, and orange lines represent the findings from this research, alongside the studies of Chang et al. [<a href="#B23-atmosphere-15-00951" class="html-bibr">23</a>], Chen et al. [<a href="#B28-atmosphere-15-00951" class="html-bibr">28</a>], and Wen et al. [<a href="#B13-atmosphere-15-00951" class="html-bibr">13</a>], respectively. The gray lines depict the relationship ΛD<sub>m</sub> = 4 + µ, with D<sub>m</sub> values of 1.0, 1.5, and 2.5 mm.</p>
Full article ">Figure 8
<p>Scatter plot and fitted line of radar reflectivity factor (Z) versus rainfall rate (R) for whole rainfall based on filtered disdrometer data, along with fitted lines from previous studies by Fulton et al. [<a href="#B55-atmosphere-15-00951" class="html-bibr">55</a>], Chen et al. [<a href="#B28-atmosphere-15-00951" class="html-bibr">28</a>], Wen et al. [<a href="#B13-atmosphere-15-00951" class="html-bibr">13</a>], and Chang et al. [<a href="#B23-atmosphere-15-00951" class="html-bibr">23</a>].</p>
Full article ">
21 pages, 10017 KiB  
Article
Seasonal Variation in Vertical Structure for Stratiform Rain at Mêdog Site in Southeastern Tibetan Plateau
by Jiaqi Wen, Gaili Wang, Renran Zhou, Ran Li, Suolang Zhaxi and Maqiao Bai
Remote Sens. 2024, 16(7), 1230; https://doi.org/10.3390/rs16071230 - 30 Mar 2024
Viewed by 1069
Abstract
Mêdog is located at the entrance of the water vapor channel of the Yarlung Tsangpo Great Canyon on the southeastern Tibetan Plateau (TP). In this study, the seasonal variation in the microphysical vertical structure of stratiform precipitation at the Mêdog site in 2022 [...] Read more.
Mêdog is located at the entrance of the water vapor channel of the Yarlung Tsangpo Great Canyon on the southeastern Tibetan Plateau (TP). In this study, the seasonal variation in the microphysical vertical structure of stratiform precipitation at the Mêdog site in 2022 was investigated using micro rain radar (MRR) observations, as there is a lack of similar studies in this region. The average melting layer height is the lowest in February, after which it gradually increases, reaches its peak in August, and then gradually decreases. For lower rain categories, the vertical distribution of small drops remains uniform in winter below the melting layer. The medium-sized drops show slight increases, leading to negative gradients in the microphysical profiles. Slight or evident decreases in concentrations of small drops are observed with decreasing height in the premonsoon, monsoon, and postmonsoon seasons, likely due to significant evaporation. The radar reflectivity, rain rate, and liquid water content profiles decrease with decreasing height according to the decrease in concentrations of small drops. With increasing rain rate, the drop size distribution (DSD) displays significant variations in winter, and the fall velocity decreases rapidly with decreasing height. In the premonsoon, monsoon, and postmonsoon seasons, the concentrations of large drops significantly decrease below the melting layer because of the breakup mechanism, leading to the decreases in the fall velocity profiles with decreasing height during these seasons. Raindrops with sizes ranging from 0.3–0.5 mm are predominant in terms of the total drop number concentration in all seasons. Precipitation in winter and postmonsoon seasons is mainly characterized by small raindrops, while that in premonsoon and monsoon seasons mainly comprises medium-sized raindrops. Understanding the seasonal variation in the vertical structure of precipitation in Mêdog will improve the radar quantitative estimation and the use of microphysical parameterization schemes in numerical weather forecast models over the TP. Full article
(This article belongs to the Special Issue Advance of Radar Meteorology and Hydrology II)
Show Figures

Figure 1

Figure 1
<p>Location of the Mêdog National Climate Observatory (red solid dots) and topography (shaded, unit: m) of the Tibetan Plateau, which are superimposed with the mean vertical integral of the water vapor flux (unit: kg m<sup>−1</sup> s<sup>−1</sup>) in different seasons in 2022.</p>
Full article ">Figure 2
<p>Comparison of rain rates (5 min average) between the MRR at 150 m above ground and rain gauge data of the precipitation process from 0100 to 0800 (local standard time) on 14 May 2022.</p>
Full article ">Figure 3
<p>(<b>a</b>) The GFV time series calculated from the MRR based on precipitation events from 0000 LT to 0700 LST on 10 February 2022. (<b>b</b>) The corresponding vertical profile of reflectivity from the MRR, in which the solid red line (the solid black line) marks the BB bottom calculated by the maximum GFV (the BB top calculated by the gradient of Ze).</p>
Full article ">Figure 3 Cont.
<p>(<b>a</b>) The GFV time series calculated from the MRR based on precipitation events from 0000 LT to 0700 LST on 10 February 2022. (<b>b</b>) The corresponding vertical profile of reflectivity from the MRR, in which the solid red line (the solid black line) marks the BB bottom calculated by the maximum GFV (the BB top calculated by the gradient of Ze).</p>
Full article ">Figure 4
<p>Monthly variation in the BB bottom height based on the GFV method in the Mêdog region.</p>
Full article ">Figure 5
<p>Vertical profiles of the radar reflectivity factor for each category in the four seasons. R1: 0.1 ≤ R &lt; 1.0 mm h<sup>−1</sup>, R2: 1.0 ≤ R &lt; 2.0 mm h<sup>−1</sup>, R3: 2.0 ≤ R &lt; 5.0 mm h<sup>−1</sup>, and R4: 5.0 ≤ R &lt; 10.0 mm h<sup>−1</sup>: (<b>a</b>) winter; (<b>b</b>) premonsoon; (<b>c</b>) monsoon; and (<b>d</b>) postmonsoon.</p>
Full article ">Figure 6
<p>Vertical profiles of the fall velocity for each category in the four seasons. R1: 0.1 ≤ R &lt; 1.0 mm h<sup>−1</sup>, R2: 1.0 ≤ R &lt; 2.0 mm h<sup>−1</sup>, R3: 2.0 ≤ R &lt; 5.0 mm h<sup>−1</sup>, and R4: 5.0 ≤ R &lt; 10.0 mm h<sup>−1</sup>: (<b>a</b>) winter; (<b>b</b>) premonsoon; (<b>c</b>) monsoon; and (<b>d</b>) postmonsoon.</p>
Full article ">Figure 7
<p>Vertical profiles of the rain rate for each category in the four seasons. R1: 0.1 ≤ R &lt; 1.0 mm h<sup>−1</sup>, R2: 1.0 ≤ R &lt; 2.0 mm h<sup>−1</sup>, R3: 2.0 ≤ R &lt; 5.0 mm h<sup>−1</sup>, and R4: 5.0 ≤ R &lt; 10.0 mm h<sup>−1</sup>: (<b>a</b>) winter; (<b>b</b>) premonsoon; (<b>c</b>) monsoon; and (<b>d</b>) postmonsoon.</p>
Full article ">Figure 8
<p>Vertical profiles of the liquid water content for each category in the four seasons. R1: 0.1 ≤ R &lt; 1.0 mm h<sup>−1</sup>, R2: 1.0 ≤ R &lt; 2.0 mm h<sup>−1</sup>, R3: 2.0 ≤ R &lt; 5.0 mm h<sup>−1</sup>, and R4: 5.0 ≤ R &lt; 10.0 mm h<sup>−1</sup>: (<b>a</b>) winter; (<b>b</b>) premonsoon; (<b>c</b>) monsoon; and (<b>d</b>) postmonsoon.</p>
Full article ">Figure 8 Cont.
<p>Vertical profiles of the liquid water content for each category in the four seasons. R1: 0.1 ≤ R &lt; 1.0 mm h<sup>−1</sup>, R2: 1.0 ≤ R &lt; 2.0 mm h<sup>−1</sup>, R3: 2.0 ≤ R &lt; 5.0 mm h<sup>−1</sup>, and R4: 5.0 ≤ R &lt; 10.0 mm h<sup>−1</sup>: (<b>a</b>) winter; (<b>b</b>) premonsoon; (<b>c</b>) monsoon; and (<b>d</b>) postmonsoon.</p>
Full article ">Figure 9
<p>Vertical structure of the DSD in different rain categories in the winter (<b>a</b>–<b>d</b>), premonsoon (<b>e</b>–<b>h</b>), monsoon (<b>i</b>–<b>l</b>), and postmonsoon (<b>m</b>–<b>p</b>) seasons.</p>
Full article ">Figure 10
<p>Seasonal variation in the average DSD at several heights.</p>
Full article ">Figure 11
<p>Contributions of different diameter categories to the total number concentration <span class="html-italic">N<sub>t</sub></span> (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) and precipitation intensity <span class="html-italic">R</span> (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) for stratiform samples in each season.</p>
Full article ">
19 pages, 7705 KiB  
Article
Spatial Variability of Raindrop Size Distribution at Beijing City Scale and Its Implications for Polarimetric Radar QPE
by Zhe Zhang, Huiqi Li, Donghuan Li and Youcun Qi
Remote Sens. 2023, 15(16), 3964; https://doi.org/10.3390/rs15163964 - 10 Aug 2023
Cited by 5 | Viewed by 1389
Abstract
Understanding the characteristics of the raindrop size distribution (DSD) is crucial to improve our knowledge of the microphysical processes of precipitation and to improve the accuracy of radar quantitative precipitation estimation (QPE). In this study, the spatial variability of DSD in different regions [...] Read more.
Understanding the characteristics of the raindrop size distribution (DSD) is crucial to improve our knowledge of the microphysical processes of precipitation and to improve the accuracy of radar quantitative precipitation estimation (QPE). In this study, the spatial variability of DSD in different regions of Beijing and its influence on radar QPE are analyzed using 11 disdrometers. The DSD data are categorized into three regions: Urban, suburban, and mountainous according to their locations. The DSD exhibits evidently different characteristics in the urban, suburban, and mountain regions of Beijing. The average raindrop diameter is smaller in the urban region compared to the suburban region. The average rain rate and raindrop number concentration are lower in the mountainous region compared to both urban and suburban regions. The difference in DSD between urban and suburban regions is due to the difference in DSD for the same precipitation types, while the difference in DSD between mountain and plains (i.e., urban and suburban regions) is the combined effect of the convection/stratiform ratio and the difference of DSD for the same precipitation types. Three DSD-based polarimetric radar QPE estimators were retrieved and estimated. Among these three QPE estimators, R(ZH), R(Kdp), and R(Kdp, ZDR), R(Kdp, ZDR) performs best, followed by R(Kdp), and R(ZH) performs worst. R(Kdp) is more sensitive to the representative parameters, while R(ZH) and R(Kdp, ZDR) are more sensitive to observational error and systematic bias (i.e., calibration). Full article
(This article belongs to the Special Issue Processing and Application of Weather Radar Data)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Location of Beijing in China and (<b>b</b>) topography of Beijing and locations of the disdrometers used in this study. The thin black lines in (<b>b</b>) denote the 6th Ring Road of Beijing.</p>
Full article ">Figure 2
<p>(<b>a</b>) Scatter density plot for <span class="html-italic">R</span> versus <span class="html-italic">D<sub>m</sub></span>, superimposed with the power–law relationship obtained using the least-square fit method and (<b>b</b>) scatter plot for <span class="html-italic">D<sub>m</sub></span> versus <span class="html-italic">N<sub>w</sub></span>. Red (blue) dots represent convection (stratiform). The star and square symbols represent the mean values for convection and stratiform, respectively. The black line is the log10(<span class="html-italic">N<sub>w</sub></span>)–<span class="html-italic">D<sub>m</sub></span> relationship for stratiform in BR03. Two rectangles indicate the maritime and continental convective clusters in BR03.</p>
Full article ">Figure 3
<p>Average raindrop spectra for different areas of Beijing.</p>
Full article ">Figure 4
<p>The probability distribution functions (PDF) of (<b>a</b>) <span class="html-italic">D<sub>m</sub></span>, (<b>b</b>) <span class="html-italic">N<sub>t</sub></span>, and (<b>c</b>) <span class="html-italic">R</span> for different areas of Beijing.</p>
Full article ">Figure 5
<p>Scatter density plots of <span class="html-italic">R</span> from all 11 disdrometer observations and R estimated using estimators: (<b>a</b>) R(<span class="html-italic">Z</span><sub>H</sub>), (<b>b</b>) R(<span class="html-italic">K</span><sub>dp</sub>), and (<b>c</b>) R(<span class="html-italic">K</span><sub>dp</sub>, <span class="html-italic">Z</span><sub>DR</sub>). The black line in each panel is the perfect fit line (i.e., y = x). Statistical scores of CC, RMSE, and RMB are superimposed.</p>
Full article ">Figure 6
<p>Scatter density plots of <span class="html-italic">R</span> in the whole region of Beijing from 11 disdrometer observations and <span class="html-italic">R</span> estimated using estimator R(<span class="html-italic">Z</span><sub>H</sub>): (<b>a</b>) Control experiment, (<b>b</b>) DSD variability experiment, (<b>c</b>) measurement error experiment, and (<b>d</b>) systematic bias experiment as described in <a href="#remotesensing-15-03964-t004" class="html-table">Table 4</a>. The black line in each panel is the perfect fit line (i.e., y = x). Statistical scores of CC, RMSE, and RMB are superimposed.</p>
Full article ">Figure 7
<p>Scatter density plots of R in the whole region of Beijing from 11 disdrometer observations and R estimated using estimator R(<span class="html-italic">K</span><sub>dp</sub>): (<b>a</b>) Control experiment, (<b>b</b>) DSD variability experiment, (<b>c</b>) measurement error experiment 1, and (<b>d</b>) measurement error experiment 2 as described in <a href="#remotesensing-15-03964-t005" class="html-table">Table 5</a>. The black line in each panel is the perfect fit line (i.e., y = x). Statistical scores of CC, RMSE, and RMB are superimposed.</p>
Full article ">Figure 8
<p>Scatter density plots of R in the whole region of Beijing from 11 disdrometer observations and R estimated using estimator R(K<sub>dp</sub>,Z<sub>dr</sub>): (<b>a</b>) Control experiment, (<b>b</b>) DSD variability experiment, (<b>c</b>) measurement error experiment, and (<b>d</b>) systematic bias experiment as described in <a href="#remotesensing-15-03964-t006" class="html-table">Table 6</a>. The black line in each panel is the perfect fit line (i.e., y = x). Statistical scores of CC, RMSE, and RMB are superimposed.</p>
Full article ">
17 pages, 49393 KiB  
Article
An Improved S-Band Polarimetric Radar-Based QPE Algorithm for Typhoons over South China Using 2DVD Observations
by Zeyong Guo, Sheng Hu, Guangyu Zeng, Xingdeng Chen, Honghao Zhang, Feng Xia, Jiahui Zhuang, Min Chen and Yuwen Fan
Atmosphere 2023, 14(6), 935; https://doi.org/10.3390/atmos14060935 - 26 May 2023
Cited by 1 | Viewed by 1286
Abstract
Polarimetric radar data are an important tool for quantitative precipitation estimation (QPE), which is essential for monitoring and forecasting precipitation. Previous studies have shown that the drop size distribution (DSD) and polarimetric radar parameters of typhoon-induced precipitation differ significantly from those of other [...] Read more.
Polarimetric radar data are an important tool for quantitative precipitation estimation (QPE), which is essential for monitoring and forecasting precipitation. Previous studies have shown that the drop size distribution (DSD) and polarimetric radar parameters of typhoon-induced precipitation differ significantly from those of other types of rainfall. South China is a region that frequently experiences typhoons and heavy rainfall, which can cause serious disasters. Therefore, it is critical to develop a QPE algorithm that is suitable for typhoon precipitation over South China. In this study, we constructed four simple QPE estimators, R(ZH), R(ZH, ZDR), R(KDP) and R(KDP, ZDR) based on two-dimensional video disdrometer (2DVD) DSD observations of typhoon-induced precipitation over South China in 2017–2018. We analyzed the DSD characteristics and the estimation accuracy of these four QPE estimators in the reflectivity–differential reflectivity (ZH–ZDR) space, as well as the S-band polarimetric radar (S-POL) data of seven typhoon-induced precipitation events that affected South China in 2017–2019. We used these data to quantitatively determine the optimal ranges of the estimators and establish a typhoon precipitation QPE algorithm for typhoon-induced precipitation over South China (2DVD-Typhoon). The evaluation results showed that: (1) compared to R(ZH) and R(KDP), R(ZH, ZDR) and R(KDP, ZDR) had lower performance in estimating typhoon-induced rainfall after incorporating the polarimetric parameter ZDR, as strong crosswind of the typhoon caused some bias in the raindrop-induced ZDR; (2) the 2DVD-Typhoon algorithm utilizes the respective advantages of the individual estimators to generate the best QPE results; (3) the QPE performance of 2DVD-Typhoon and the Colorado State University–Hydrometeor Identification Rainfall Optimization (CSU-HIDRO) is used as a comparison for hourly rainfall, cumulative rainfall and different rainfall intensity. The comparison shows that 2DVD-Typhoon gives a better normalized error (NE), root mean square error (RMSE) and correlation coefficient (CC), indicating its strength in rainfall estimation for typhoons over South China. The above results provide theoretical support for improving typhoon-induced rainfall monitoring and numerical weather forecasting models in South China. Full article
(This article belongs to the Special Issue Monsoon and Typhoon Precipitation in Asia: Observation and Prediction)
Show Figures

Figure 1

Figure 1
<p>Trajectory and intensity change of typhoons (TD: tropical depression, TS: tropical storm, STS: severe tropical storm, TY: typhoon, STY: severe typhoon).</p>
Full article ">Figure 2
<p>The location map of the two radars, 2DVDs and surface rain gauge stations in Guangzhou and Yangjiang. (Red crosses represent the locations of the two S-band polarimetric radars (S-POL) in Guangzhou and Yangjiang, Guangdong Province; green dots indicate the locations of the rain gauges within the area covered by the radars (5–100 km radius); blue diamonds represent the 2DVDs of the Longmen Cloud Physics Field Experiment Base of the China Meteorological Administration (CMA)).</p>
Full article ">Figure 3
<p>Quality control process of rain gauge data.</p>
Full article ">Figure 4
<p>Technical procedures for establishing the South China typhoon QPE scheme (2DVD-Typhoon).</p>
Full article ">Figure 5
<p>Scatterplot of rain rate retrieved from (<b>a</b>) R(Z<sub>H</sub>), (<b>b</b>) R(Z<sub>H</sub>, Z<sub>DR</sub>), (<b>c</b>) R(K<sub>DP</sub>) and (<b>d</b>) R(K<sub>DP</sub>, Z<sub>DR</sub>) vs. the directly observed rain rate from the 2DVDs during typhoon precipitation from 2017 to 2018.</p>
Full article ">Figure 6
<p>The total sample number (<b>a</b>), hourly QPE normalized error of (<b>b</b>) R(Z<sub>H</sub>), (<b>c</b>) R(Z<sub>H</sub>, Z<sub>DR</sub>), (<b>d</b>) R(K<sub>DP</sub>) and (<b>e</b>) R(K<sub>DP</sub>, Z<sub>DR</sub>) for all seven typhoon precipitation events in Z<sub>H</sub>–Z<sub>DR</sub> space; red lines in (<b>b</b>,<b>c</b>) represent the optimal threshold ranges of the four single QPE estimators used to calculate the composite QPE algorithm (2DVD-Typhoon).</p>
Full article ">Figure 6 Cont.
<p>The total sample number (<b>a</b>), hourly QPE normalized error of (<b>b</b>) R(Z<sub>H</sub>), (<b>c</b>) R(Z<sub>H</sub>, Z<sub>DR</sub>), (<b>d</b>) R(K<sub>DP</sub>) and (<b>e</b>) R(K<sub>DP</sub>, Z<sub>DR</sub>) for all seven typhoon precipitation events in Z<sub>H</sub>–Z<sub>DR</sub> space; red lines in (<b>b</b>,<b>c</b>) represent the optimal threshold ranges of the four single QPE estimators used to calculate the composite QPE algorithm (2DVD-Typhoon).</p>
Full article ">Figure 7
<p>The ratio diagrams of estimated cumulative precipitation derived from (<b>a</b>) R(Z<sub>H</sub>, Z<sub>DR</sub>), (<b>b</b>) R(K<sub>DP,</sub>, Z<sub>DR</sub>), (<b>c</b>) R(Z<sub>H</sub>) and (<b>d</b>) R(K<sub>DP</sub>) to rain gauge cumulative precipitation (The stars in the figure are the locations of the two S-band polarimetric radars (S-POL) in Guangzhou and Yangjiang, Guangdong Province. The scatter points in the figure are the locations of the rain gauge stations. The size of the scatter point represents the cumulative rainfall of the station and the color of the scatter point represents the ratio of QPE to the rain gauge. The filled contours represent the average instantaneous wind speed).</p>
Full article ">Figure 7 Cont.
<p>The ratio diagrams of estimated cumulative precipitation derived from (<b>a</b>) R(Z<sub>H</sub>, Z<sub>DR</sub>), (<b>b</b>) R(K<sub>DP,</sub>, Z<sub>DR</sub>), (<b>c</b>) R(Z<sub>H</sub>) and (<b>d</b>) R(K<sub>DP</sub>) to rain gauge cumulative precipitation (The stars in the figure are the locations of the two S-band polarimetric radars (S-POL) in Guangzhou and Yangjiang, Guangdong Province. The scatter points in the figure are the locations of the rain gauge stations. The size of the scatter point represents the cumulative rainfall of the station and the color of the scatter point represents the ratio of QPE to the rain gauge. The filled contours represent the average instantaneous wind speed).</p>
Full article ">Figure 7 Cont.
<p>The ratio diagrams of estimated cumulative precipitation derived from (<b>a</b>) R(Z<sub>H</sub>, Z<sub>DR</sub>), (<b>b</b>) R(K<sub>DP,</sub>, Z<sub>DR</sub>), (<b>c</b>) R(Z<sub>H</sub>) and (<b>d</b>) R(K<sub>DP</sub>) to rain gauge cumulative precipitation (The stars in the figure are the locations of the two S-band polarimetric radars (S-POL) in Guangzhou and Yangjiang, Guangdong Province. The scatter points in the figure are the locations of the rain gauge stations. The size of the scatter point represents the cumulative rainfall of the station and the color of the scatter point represents the ratio of QPE to the rain gauge. The filled contours represent the average instantaneous wind speed).</p>
Full article ">Figure 8
<p>Scattered plots comparing the two algorithms: hourly rainfall (<b>a</b>) and accumulated rainfall (<b>b</b>) of 2DVD-Typhoon, hourly rainfall (<b>c</b>) and accumulated rainfall (<b>d</b>) of CSU-HIDRO.</p>
Full article ">Figure 9
<p>Evaluation of the two algorithms under different rainfall intensities.(Comparison of, (<b>a</b>) NE, (<b>b</b>) RMSE, (<b>c</b>) cc, between 2DVD-Typhoon and CSU-HIDRO).</p>
Full article ">
22 pages, 7128 KiB  
Article
Regional Variability of Raindrop Size Distribution from a Network of Disdrometers over Complex Terrain in Southern China
by Asi Zhang, Chao Chen and Lin Wu
Remote Sens. 2023, 15(10), 2678; https://doi.org/10.3390/rs15102678 - 21 May 2023
Cited by 2 | Viewed by 1619
Abstract
Raindrop size distribution (DSD) over the complex terrain of Guangdong Province, southern China, was studied using six disdrometers operated by the Guangdong Meteorology Service during the period 1 March 2018 to 30 August 2022 (~5 years). To analyze the long-term DSD characteristics over [...] Read more.
Raindrop size distribution (DSD) over the complex terrain of Guangdong Province, southern China, was studied using six disdrometers operated by the Guangdong Meteorology Service during the period 1 March 2018 to 30 August 2022 (~5 years). To analyze the long-term DSD characteristics over complex topography in southern China, three stations on the windward side, Haifeng, Enping and Qingyuan, and three stations on the leeward side, Meixian, Luoding and Xuwen, were utilized. The median mass-weighted diameter (Dm) value was higher on the windward than on the leeward side, and the windward-side stations also showed greater Dm variability. With regard to the median generalized intercept (log10Nw) value, the log10Nw values decreased from coastal to mountainous areas. Although there were some differences in Dm, log10Nw and liquid water content (LWC) frequency between the six stations, there were still some similarities, with the Dm, log10Nw and LWC frequency all showing a single-peak curve. In addition, the diurnal variation of the mean log10Nw had a negative relationship with Dm diurnal variation although the inverse relationship was not particularly evident at the Haifeng site. The diurnal mean rainfall rate also peaked in the afternoon and exceeded the maximum at night which indicated that strong land heating in the daytime significantly influenced the local DSD variation. What is more, the number concentration of drops, N(D), showed an exponential shape which decreased monotonically for all rainfall rate types at the six observation sites, and an increase in diameter caused by increases in the rainfall rate was also noticeable. As the rainfall rate increased, the N(D) for sites on the windward side (i.e., Haifeng, Enping and Qingyuan) were higher than for the sites on the leeward side (i.e., Meixian, Luoding and Xuwen), and the difference between them also became distinct. The abovementioned DSD characteristic differences also showed appreciable variability in convective precipitation between stations on the leeward side (i.e., Meixian, Luoding and Xuwen) and those on the windward side (Haifeng and Enping, but not Qingyuan). This study enhances the precision of numerical weather forecast models in predicting precipitation and verifies the accuracy of measuring precipitation through remote sensing instruments, including weather radars located on the ground. Full article
(This article belongs to the Special Issue Remote Sensing of Clouds and Precipitation at Multiple Scales II)
Show Figures

Figure 1

Figure 1
<p>Guangdong province and its digital elevation model (DEM). The red symbols are the disdrometer locations used in this study. The locations of Nanling, Lianhua, Tianlu, Yunwu and Wuzhi Mountains are also displayed.</p>
Full article ">Figure 2
<p>Scatter plot of hourly accumulated rainfall collected from the HY-P1000 disdrometers and rain gauges from automatic weather stations.</p>
Full article ">Figure 3
<p>Box−and−whisker plot of (<b>a</b>) <span class="html-italic">D<sub>m</sub></span> and (<b>b</b>) log<sub>10</sub><span class="html-italic">N<sub>w</sub></span> distribution over the six different sites. The box represents the data between the 25th and 75th percentiles and the whiskers show the maximum and minimum values. The horizontal red line within the box represents the median value of the distribution.</p>
Full article ">Figure 4
<p>Frequency of (<b>a</b>) <span class="html-italic">D<sub>m</sub></span>, (<b>b</b>) log<sub>10</sub><span class="html-italic">N<sub>w</sub></span> and (<b>c</b>) <span class="html-italic">LWC</span> over the Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan areas.</p>
Full article ">Figure 5
<p>Diurnal variation of the mean (<b>a</b>) <span class="html-italic">D<sub>m</sub></span> (<b>b</b>) log<sub>10</sub><span class="html-italic">N<sub>w</sub></span> and (<b>c</b>) rain rate over the Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan areas.</p>
Full article ">Figure 6
<p>The number variation for every hour over the Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan areas.</p>
Full article ">Figure 7
<p>Variation of the mean raindrop concentration and raindrop diameter (Diameter, mm) in Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan over five years.</p>
Full article ">Figure 8
<p>Average raindrop spectra of Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan for rainfall in 12 rainfall rate classes (<b>a</b>) C1: 0.1–0.2, (<b>b</b>) C2: 0.2–0.4, (<b>c</b>) C3: 0.4–0.7, (<b>d</b>) C4: 0.7–1.0, (<b>e</b>) C5: 1.0–2.0, (<b>f</b>) C6: 2.0–5.0, (<b>g</b>) C7: 5.0–8.0, (<b>h</b>) C8: 8.0–12.0, (<b>i</b>) C9: 12.0–18.0, (<b>j</b>) C10: 18.0–25.0, (<b>k</b>) C11: 25.0–40.0 and (<b>l</b>) C12: &gt;40 mm h<sup>−1</sup>.</p>
Full article ">Figure 8 Cont.
<p>Average raindrop spectra of Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan for rainfall in 12 rainfall rate classes (<b>a</b>) C1: 0.1–0.2, (<b>b</b>) C2: 0.2–0.4, (<b>c</b>) C3: 0.4–0.7, (<b>d</b>) C4: 0.7–1.0, (<b>e</b>) C5: 1.0–2.0, (<b>f</b>) C6: 2.0–5.0, (<b>g</b>) C7: 5.0–8.0, (<b>h</b>) C8: 8.0–12.0, (<b>i</b>) C9: 12.0–18.0, (<b>j</b>) C10: 18.0–25.0, (<b>k</b>) C11: 25.0–40.0 and (<b>l</b>) C12: &gt;40 mm h<sup>−1</sup>.</p>
Full article ">Figure 9
<p>Distribution of (<b>a</b>) <span class="html-italic">D<sub>m</sub></span> and (<b>b</b>) log<sub>10</sub><span class="html-italic">N<sub>w</sub></span> for rainfall at Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan in relation to rainfall rates (C1: 0.1–0.2, C2: 0.2–0.4, C3: 0.4–0.7, C4: 0.7–1.0, C5: 1.0–2.0, C6: 2.0–5.0, C7: 5.0–8.0, C8: 8.0–12.0, C9: 12.0–18.0, C10: 18.0–25.0, C11: 25.0–40.0 and C12: &gt;40 mm h<sup>−1</sup>).</p>
Full article ">Figure 10
<p>The scatter plots of D<sub>m</sub> vs. rain rate for (<b>a</b>) stratiform and (<b>b</b>) convective precipitation and log<sub>10</sub>N<sub>w</sub> vs. rain rate for (<b>c</b>) stratiform and (<b>d</b>) convective precipitation.</p>
Full article ">Figure 11
<p>Average raindrop spectra of stratiform and convective precipitation at Meixian, Haifeng, Luoding, Enping, Xuwen and Qingyuan stations. the dashed line represents the stratiform precipitation and the solid line represents the convective precipitation.</p>
Full article ">Figure 12
<p>Radar reflectivity <span class="html-italic">Z</span> and rainfall rate <span class="html-italic">R</span> relationship of stratiform and convective precipitation at Meixian, Haifeng, Luoding, Enping, Xuwen, Qingyuan stations and for all data.</p>
Full article ">
19 pages, 7143 KiB  
Article
Microphysical Characteristics of Raindrop Size Distribution and Implications for Dual-Polarization Radar Quantitative Precipitation Estimations in the Tianshan Mountains, China
by Yong Zeng, Jiangang Li, Lianmei Yang, Haoyang Li, Xiaomeng Li, Zepeng Tong, Yufei Jiang, Jing Liu, Jinru Zhang and Yushu Zhou
Remote Sens. 2023, 15(10), 2668; https://doi.org/10.3390/rs15102668 - 20 May 2023
Cited by 5 | Viewed by 3093
Abstract
In order to improve the understanding of the microphysical characteristics of raindrop size distribution (DSD) under different rainfall rates (R) classes, and broaden the knowledge of the impact of radar wavelengths and R classes on the QPE of dual-polarization radars in [...] Read more.
In order to improve the understanding of the microphysical characteristics of raindrop size distribution (DSD) under different rainfall rates (R) classes, and broaden the knowledge of the impact of radar wavelengths and R classes on the QPE of dual-polarization radars in the Tianshan Mountains, a typical arid area in China, we investigated the microphysical characteristics of DSD across R classes and dual-polarimetric radar QPE relationships across radar wavelengths and R classes, based on the DSD data from a PARSIVEL2 disdrometer at Zhaosu in the Tianshan Mountains during the summers of 2020 and 2021. As the R class increased, the DSD became wider and flatter. The mean value of the mass-weighted mean diameters (Dm) increased, while the mean value of logarithm normalized intercept parameters (log10 Nw) decreased after increasing from C1 to C3, as the R class increased. The largest contributions to R and the radar reflectivity factor from large raindrops (diameter > 3 mm) accounted for approximately 50% and 97%, respectively, while 84% of the total raindrops were small raindrops (diameter < 1 mm). Dual-polarization radars—horizontal polarization reflectivity (Zh), differential reflectivity (Zdr), and specific differential phase (Kdp)—were retrieved based on the DSD data using the T-matrix scattering method. The DSD-based polarimetric radar QPE relations of a single-parameter (R(Zh), R(Kdp)), and double-parameters (R(Zh,Zdr), R(Kdp,Zdr)) on the S-, C-, and X-bands were derived and evaluated. Overall, the performance of the R(Kdp) (R(Kdp,Zdr)) scheme was better than that of R(Zh) (R(Zh,Zdr)) for the QPE in the three bands. Furthermore, we have for the first time confirmed and quantified the performance differences in the QPE relationship of dual-polarization radars under different schemes, radar wavelengths, and R classes in typical arid areas of China. Therefore, selecting an appropriate dual-polarization radar band and QPE scheme for different R classes is necessary to improve the QPE ability compared with an independent scheme under all R classes. Full article
(This article belongs to the Special Issue Processing and Application of Weather Radar Data)
Show Figures

Figure 1

Figure 1
<p>Location of Zhaosu (the black dot), with shading representing the topography (m) of the Tianshan Mountains.</p>
Full article ">Figure 2
<p>Accumulated rain duration (red histogram) and rain amount (blue line) for the six <span class="html-italic">R</span> classes in Zhaosu.</p>
Full article ">Figure 3
<p>Mean DSD variations for different <span class="html-italic">R</span> classes (color lines) and all samples (black line) in Zhaosu. The two vertical dashed lines on the left and right distinguish the raindrop spectrum of small and medium-size raindrops, and the raindrop spectrum of medium-size and large raindrops, respectively.</p>
Full article ">Figure 4
<p>Variations of the <span class="html-italic">D<sub>m</sub></span> and the log<sub>10</sub><span class="html-italic">N<sub>w</sub></span> in Zhaosu for the six <span class="html-italic">R</span> classes. The line and dot of the box indicate the mean (black line) and median (black dot), respectively. The bottom (top) lines of the box indicate the 25th (75th) percentiles. The bottom (top) lines of the vertical lines out of the box indicate the 5th (95th) percentiles.</p>
Full article ">Figure 5
<p>The contribution of small, medium, and large drops to <span class="html-italic">R</span>, <span class="html-italic">Z</span>, <span class="html-italic">LWC</span>, and <span class="html-italic">N<sub>t</sub></span> in Zhaosu.</p>
Full article ">Figure 6
<p>Scatterplots of <span class="html-italic">Z<sub>dr</sub></span> versus <span class="html-italic">Z<sub>h</sub></span>, and the <span class="html-italic">Z<sub>dr</sub></span>–<span class="html-italic">Z<sub>h</sub></span> relations represented by black line and equation on the (<b>a</b>) S-band, (<b>c</b>) C-band, and (<b>e</b>) X-band. Scatterplots of <span class="html-italic">K<sub>dp</sub></span> versus <span class="html-italic">Z<sub>h</sub></span>, and the <span class="html-italic">K<sub>dp</sub></span>–<span class="html-italic">Z<sub>h</sub></span> relations represented by black line and equation on the (<b>b</b>) S-band, (<b>d</b>) C-band, and (<b>f</b>) X-band.</p>
Full article ">Figure 7
<p>Scatterplot of <span class="html-italic">R</span> calculated from (<b>a</b>) <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>), (<b>b</b>) <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>), (<b>c</b>) <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>), and (<b>d</b>) <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) relations versus the <span class="html-italic">R</span> computed from DSD for S-band in Zhaosu.</p>
Full article ">Figure 8
<p>Scatterplot of <span class="html-italic">R</span> calculated from (<b>a</b>) <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>), (<b>b</b>) <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>), (<b>c</b>) <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>), and (<b>d</b>) <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) relations versus the <span class="html-italic">R</span> computed from DSD for C-band in Zhaosu.</p>
Full article ">Figure 9
<p>Scatterplot of <span class="html-italic">R</span> calculated from (<b>a</b>) <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span>), (<b>b</b>) <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span>), (<b>c</b>) <span class="html-italic">R</span>(<span class="html-italic">Z<sub>h</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>), and (<b>d</b>) <span class="html-italic">R</span>(<span class="html-italic">K<sub>dp</sub></span><sub>,</sub><span class="html-italic">Z<sub>dr</sub></span>) relations versus the <span class="html-italic">R</span> computed from DSD for X-band in Zhaosu.</p>
Full article ">Figure 10
<p>Variations of (<b>a</b>) <span class="html-italic">Z<sub>h</sub></span>, (<b>b</b>) <span class="html-italic">Z<sub>dr</sub></span>, and (<b>c</b>) <span class="html-italic">K<sub>dp</sub></span> on the S- (red), C- (green), and X-bands (purple) for the six <span class="html-italic">R</span> classes. The line of the box indicates the mean. The bottom (top) lines of the box indicate the 25th (75th) percentiles. The bottom (top) lines of the vertical lines out of the box indicate the 5th (95th) percentiles.</p>
Full article ">Figure 11
<p>The (<b>a</b>,<b>d</b>,<b>g</b>) CC, (<b>b</b>,<b>e</b>,<b>h</b>) RMSE, and (<b>c</b>,<b>f</b>,<b>i</b>) NMAE of <span class="html-italic">R</span> estimated from the dual-polarization radar QPE estimators against <span class="html-italic">R</span> calculated from the DSD under different <span class="html-italic">R</span> classes and different radar bands, (<b>a</b>–<b>c</b>) for S-band, (<b>d</b>–<b>f</b>) for C-band, and (<b>g</b>–<b>i</b>) for X-band, respectively.</p>
Full article ">
22 pages, 12411 KiB  
Article
Evaluating Simulated Microphysics of Stratiform and Convective Precipitation in a Squall Line Event Using Polarimetric Radar Observations
by Yuting Sun, Zhimin Zhou, Qingjiu Gao, Hongli Li and Minghuan Wang
Remote Sens. 2023, 15(6), 1507; https://doi.org/10.3390/rs15061507 - 9 Mar 2023
Cited by 3 | Viewed by 1973
Abstract
Recent upgrades to China’s radar network now allow for polarimetric measurements of convective systems in central China, providing an effective data set with which to evaluate the microphysics schemes employed in local squall line simulations. We compared polarimetric radar variables derived by Weather [...] Read more.
Recent upgrades to China’s radar network now allow for polarimetric measurements of convective systems in central China, providing an effective data set with which to evaluate the microphysics schemes employed in local squall line simulations. We compared polarimetric radar variables derived by Weather Research and Forecasting (WRF) and radar forward models and the corresponding hydrometeor species with radar observations and retrievals for a severe squall line observed over central China on 16 March 2022. Two microphysics schemes were tested and were able to accurately depict the contrast between convective and stratiform regions in terms of the drop size distribution (DSD) and reproduce the classical polarimetric signatures of the observed differential reflectivity (ZDR) and specific differential phase (KDP) columns. However, for the convective region, the simulated DSDs in both schemes exhibited lower proportions of large drops and lower liquid water content; by contrast, for the stratiform region, the proportion of large drops was found to be too high in the Morrison (MORR) scheme. The underprediction of ice-phase processes in the convective region, particularly the riming processes associated with graupel and hail, was likely responsible for the bias toward large raindrops at low levels. In the stratiform region, raindrop evaporation in the WRF Double-Moment 6-Class (WDM6) scheme, which partially offsets the overestimation of ice-phase processes, produced ground DSDs that more closely matched the observational data, and did not exhibit the overly strong warm-rain collisional growth processes of MORR. Full article
(This article belongs to the Special Issue Processing and Application of Weather Radar Data)
Show Figures

Figure 1

Figure 1
<p>The model domains used. The outer region is the parent domain (9 km), and d02 is the inner domain (3 km). The innermost box, shaded green, denotes the area of the squall line system (approximately 27°–31°N, 107°–113°E), and the triangle at the bottom of the green area indicates the location of the Huaihua station.</p>
Full article ">Figure 2
<p>The ERA5 data for (<b>a</b>) 500 hPa geopotential height (solid black, contoured every 30 m) and wind barbs (m/s) superimposed on total column water (kg/m<sup>2</sup>; grayscale shading) and (<b>b</b>) 850 hPa equivalent potential temperature (K; color shading) and wind barbs (m/s) at 0000 UTC on 16 March 2022. The gray shading in (<b>b</b>) denotes terrain height &gt; 1.5 km. The white rectangles in (<b>a</b>,<b>b</b>) represent the area of the squall line system.</p>
Full article ">Figure 3
<p>Skew T-logP diagram for Huaihua (triangle in <a href="#remotesensing-15-01507-f001" class="html-fig">Figure 1</a>) at 0000 UTC on 16 March 2022. The red and blue lines show temperature and dewpoint temperature profiles, respectively. The gray curve represents the ascending path of the most unstable parcel. The yellow background line, sloped at a 45° angle, denotes temperature lines.</p>
Full article ">Figure 4
<p>Radar mosaic showing radar reflectivity (dBZ) at (<b>a</b>) 0900 UTC, (<b>b</b>) 1200 UTC, (<b>c</b>) 1400 UTC, and (<b>d</b>) 1530 UTC on 16 March 2022. The triangle denotes the location of the S-POL radar.</p>
Full article ">Figure 5
<p>Spatial distribution of hourly rainfall above 20 mm/h during the mature stage of the squall line. The time values are given in the legend.</p>
Full article ">Figure 6
<p>Temporal evolution of observed and simulated radar reflectivities (dBZ) for the 16 March 2022 squall line case at 1200 UTC (left column), 1300 UTC (middle column), and 1400 UTC (right column) at an elevation of 3 km. Black contours show the hourly temperature change, indicating the location of the surface cold pool. The triangle in each panel indicates where S-POL is situated. The dark gray circles in the first three panels indicate the 150 km range from the radar.</p>
Full article ">Figure 7
<p>(<b>a</b>,<b>d</b>,<b>g</b>) Radar reflectivity (<span class="html-italic">Z</span>; dBZ), (<b>b</b>,<b>e</b>,<b>h</b>) differential reflectivity (<span class="html-italic">Z</span><sub>DR</sub>; dB), and (<b>c</b>,<b>f</b>,<b>i</b>) specific differential phase (<span class="html-italic">K</span><sub>DP</sub>; °/km) from (<b>a</b>−<b>c</b>) the S-POL radar observations and the POLARRIS-f simulations converted from the WRF output using (<b>d</b>–<b>f</b>) MORR and (<b>g</b>–<b>i</b>) WDM6 microphysics schemes, at 1300 UTC on 16 March 2022 at an elevation of 3 km. The convective region was divided based on the criteria revised by Powell et al. [<a href="#B59-remotesensing-15-01507" class="html-bibr">59</a>] and is indicated by black dots. The triangle in each panel indicates where S-POL is situated. Lines A-B in (a), (d), and (g) indicate cross-section lines of Figure 9.</p>
Full article ">Figure 8
<p>Joint radar reflectivity–differential reflectivity (<span class="html-italic">Z</span>-<span class="html-italic">Z</span><sub>DR</sub>; top) and reflectivity– specific differential phase (<span class="html-italic">Z</span>-<span class="html-italic">K</span><sub>DP</sub>; bottom) frequency distributions, normalized by maximum frequency, for 1200–1500 UTC (percent, color shadings from 2% to 100% for the convective region; black contours at 5, 30, 70, and 90% for the stratiform region). (<b>a</b>,<b>d</b>) The S-POL radar observations (below 3 km elevation) and POLARRIS-f simulations converted from the WRF output using the (<b>b</b>,<b>e</b>) MORR and (<b>c</b>,<b>f</b>) WDM6 microphysics schemes. The statistics from the simulations were limited to below approximately 2.75 km. <span class="html-italic">Z</span> was binned from 0 to 60 dBZ every 1 dBZ and both <span class="html-italic">Z</span><sub>DR</sub> and <span class="html-italic">K</span><sub>DP</sub> were binned from −0.5 to 4.5 dB (or °/km) every 0.05 dB (or °/km).</p>
Full article ">Figure 9
<p>Cross section of (<b>a</b>,<b>d</b>,<b>g</b>) radar reflectivity (<span class="html-italic">Z</span>; dBZ), (<b>b</b>,<b>e</b>,<b>h</b>) differential reflectivity (<span class="html-italic">Z</span><sub>DR</sub>; dB), and (<b>c</b>,<b>f</b>,<b>i</b>) <span class="html-italic">K</span><sub>DP</sub> (°/km) from the (<b>a</b>–<b>c</b>) S-POL radar and (<b>d</b>–<b>f</b>) MORR and (<b>g</b>–<b>i</b>) WDM6 microphysics schemes. The dashed black line indicates the freezing level and the black arrows represent wind (m/s). The black lines in <a href="#remotesensing-15-01507-f007" class="html-fig">Figure 7</a> a,d,g indicate the locations of the cross-section plots.</p>
Full article ">Figure 10
<p>Contoured frequency by altitude diagrams (CFADs) for reflectivity over (top panels) convective and (bottom panels) stratiform regions, from (<b>a</b>,<b>d</b>) S-POL radar and the (<b>b</b>,<b>e</b>) MORR and (<b>c</b>,<b>f</b>) WDM6 microphysics schemes. The time period of the analysis is 1200–1500 UTC.</p>
Full article ">Figure 11
<p>Comparison of (<b>a</b>,<b>b</b>) <span class="html-italic">Z</span>, (<b>c</b>,<b>d</b>) differential reflectivity (<span class="html-italic">Z</span><sub>DR</sub>), and (<b>e</b>,<b>f</b>) specific differential phase (<span class="html-italic">K</span><sub>DP</sub>) median profiles over (top) convective and (bottom) stratiform regions for S-POL radar and the MORR and WDM6 microphysics schemes in the sampling period 1200–1500 UTC. Error bars denote the interquartile range for each altitude layer. The gray dashed lines indicate the 0 °C, −10 °C, and −20 °C levels as recorded by soundings, with the lowest value at the bottom and the highest at the top.</p>
Full article ">Figure 12
<p>Hydrometeor identification frequency by height for convective (left) and stratiform (right) regions from (<b>a</b>,<b>b</b>) S-POL radar and the (<b>c</b>,<b>d</b>) MORR and (<b>e</b>,<b>f</b>) WDM6 microphysics schemes during 1200–1500 UTC. DZ = drizzle, RN = rain, CR = ice crystals, AG = aggregates, WS = wet snow, VI = vertical ice, LDG = low-density graupel, HDG = high-density graupel, HA = hail, BD = big drops/melting hail. The 0, −10, and −20 °C levels are indicated by gray dashed lines, as in <a href="#remotesensing-15-01507-f011" class="html-fig">Figure 11</a>.</p>
Full article ">
19 pages, 964 KiB  
Article
Study on the Backscatter Differential Phase Characteristics of X-Band Dual-Polarization Radar and its Processing Methods
by Fei Geng and Liping Liu
Remote Sens. 2023, 15(5), 1334; https://doi.org/10.3390/rs15051334 - 27 Feb 2023
Viewed by 1719
Abstract
The differential propagation phase (ΦDP) of X-band dual-polarization weather radar (including X-band dual-polarization phased-array weather radar, X-PAR) is important for estimating precipitation and classifying hydrometeors. However, the measured differential propagation phase contains the backscatter differential phase (δ), which [...] Read more.
The differential propagation phase (ΦDP) of X-band dual-polarization weather radar (including X-band dual-polarization phased-array weather radar, X-PAR) is important for estimating precipitation and classifying hydrometeors. However, the measured differential propagation phase contains the backscatter differential phase (δ), which poses difficulties for the application of the differential propagation phase from X-band radars. This paper presents the following: (1) the simulation and characteristics analysis of the backscatter differential phase based on disdrometer DSD (raindrop size distribution) measurement data; (2) an improved method of the specific differential propagation phase (KDP) estimation based on linear programming and backscatter differential phase elimination; (3) the effect of backscatter differential phase elimination on the specific differential propagation phase estimation of X-PAR. The results show the following: (1) For X-band weather radar, the raindrop equivalent diameters D > 2 mm may cause a backscatter differential phase between 0 and 20°; in particular, the backscatter differential phase varies sharply with raindrop size between 3.2 and 4.5 mm. (2) Using linear programming or smoothing filters to process the differential propagation phase could suppress the backscatter differential phase, but it is hard to completely eliminate the effect of the backscatter differential phase. (3) Backscatter differential phase correction may improve the calculation accuracy of the specific differential propagation phase, and the optimization was verified by the improved self-consistency of polarimetric variables, correlation between specific differential propagation phase estimations from S- and X-band radar and the accuracy of quantitative precipitation estimation. The X-PAR deployed in Shenzhen showed good observation performance and the potential to be used in radar mosaics with S-band weather radar. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The relationship between <span class="html-italic">δ</span> and raindrops’ equivalent diameter for (<b>a</b>) S-, (<b>b</b>) C-, and (<b>c</b>) X-band at different temperatures.</p>
Full article ">Figure 2
<p>Scatterplots of the <span class="html-italic">δ</span> calculated from DSD measurement with <span class="html-italic">Z<sub>DR</sub></span> for the (<b>a</b>) C-band and (<b>b</b>) X- band, and of <span class="html-italic">δ</span> calculated from gamma assumptions of DSD with <span class="html-italic">Z<sub>DR</sub></span> for the (<b>c</b>) C-band and (<b>d</b>) X- band. The red line represents the polynomial-fitted <span class="html-italic">δ–Z<sub>DR</sub></span> relationship curve.</p>
Full article ">Figure 3
<p>Simulated (<b>a</b>) <span class="html-italic">Z<sub>H</sub></span> and Z<sub>DR</sub> curves with distance using DSD measurements from 12:12 to 15:58 on 10 June 2019, (<b>b</b>) <span class="html-italic">Φ<sub>DP</sub></span> and <span class="html-italic">K<sub>DP</sub></span> calculated from DSD with <span class="html-italic">δ</span> effects (red curve) and without <span class="html-italic">δ</span> effects (black curve), (<b>c</b>) Rain intensity calculated by <span class="html-italic">K<sub>DP</sub></span> without <span class="html-italic">δ</span> effects (R<sub>1</sub>, blue curve) and with <span class="html-italic">δ</span> effects (R<sub>2</sub>, red curve), and (<b>d</b>) Backscatter differential phase.</p>
Full article ">Figure 4
<p>Comparison of the actual attenuation and attenuation correction for (<b>a</b>) <span class="html-italic">Z<sub>H</sub></span> and (<b>b</b>) <span class="html-italic">Z<sub>DR</sub></span>. The black line shows <span class="html-italic">Z<sub>H</sub></span> or <span class="html-italic">Z<sub>DR</sub></span>, the red solid line shows the actual PIA for reflectivity and path integrated differential attenuation (PIDA) for differential reflectivity, and the blue line shows the attenuation correction (<span class="html-italic">PIA<sub>δ</sub></span> and <span class="html-italic">PIDA<sub>δ</sub></span>); the differences between the red and blue curves are the errors in the attenuation correction caused by <span class="html-italic">δ</span>.</p>
Full article ">Figure 5
<p>(<b>a</b>) <span class="html-italic">K<sub>DP</sub></span> (unaffected by <span class="html-italic">δ</span>) and <span class="html-italic">K<sub>DPδ</sub></span> calculated from <span class="html-italic">Φ<sub>DPR</sub></span> (affected by <span class="html-italic">δ</span>); (<b>b</b>) <span class="html-italic">K<sub>DP</sub></span><sub>1</sub> after the elimination of <span class="html-italic">δ</span> from <span class="html-italic">Φ<sub>DPR</sub></span> using the <span class="html-italic">Z<sub>DR</sub></span> without attenuation effect, and <span class="html-italic">K<sub>DP</sub></span><sub>2</sub> using the corrected <span class="html-italic">Z<sub>DR</sub>, K<sub>DP</sub></span><sub>1</sub>, and <span class="html-italic">K<sub>DP</sub></span><sub>2</sub> denote the corrected <span class="html-italic">K<sub>DP</sub></span> for the two tests, respectively; (<b>c</b>) rain intensity estimated from <span class="html-italic">K<sub>DP1</sub> (R<sub>1</sub>)</span> and <span class="html-italic">K<sub>DP2</sub> (R<sub>2</sub>)</span>; (<b>d</b>) the actual rainfall intensity.</p>
Full article ">Figure 6
<p>Colored density plot of the simulated <span class="html-italic">K<sub>DP</sub></span> versus (<b>a</b>) the <span class="html-italic">K<sub>DPδ</sub></span> affected by <span class="html-italic">δ</span> and (<b>b</b>) the corrected <span class="html-italic">K<sub>DP1</sub></span> based on the <span class="html-italic">δ</span>-elimination using <span class="html-italic">Z<sub>DR</sub></span>.</p>
Full article ">Figure 7
<p>Typical radar sounding of (<b>a</b>) <span class="html-italic">Φ<sub>DP</sub></span>, <span class="html-italic">Z<sub>DR</sub></span>, and <span class="html-italic">δ</span> obtained by <span class="html-italic">Z<sub>DR</sub></span>, and <span class="html-italic">Φ<sub>DP2</sub></span> obtained by <span class="html-italic">δ</span> correction; (<b>b</b>) LP-processed results obtained before (blue curve) and after (red curve) <span class="html-italic">δ</span> correction; (<b>c</b>) <span class="html-italic">Z<sub>H</sub></span> after attenuation correction (dBZ) and <span class="html-italic">K<sub>DP</sub></span> before and after <span class="html-italic">δ</span> correction using LS; and (<b>d</b>) <span class="html-italic">Z<sub>H</sub></span> after attenuation correction (dBZ) and <span class="html-italic">K<sub>DP</sub></span> before and after <span class="html-italic">δ</span> correction using LP.</p>
Full article ">Figure 8
<p>Radar PPI of (<b>a</b>) <span class="html-italic">Z<sub>H</sub></span> from S-POL, (<b>b</b>) <span class="html-italic">Z<sub>DR</sub></span> from S-POL, (<b>c</b>) <span class="html-italic">K<sub>DP</sub></span> from S-POL, (<b>d</b>) <span class="html-italic">Z<sub>H</sub></span> from X-PAR, (<b>e</b>) <span class="html-italic">Z<sub>DR</sub></span> from X-PAR, (<b>f</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp1 (LS without <span class="html-italic">δ</span>-elimination), (<b>g</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp2 (LS with <span class="html-italic">δ</span>-elimination), (<b>h</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp3 (LP without <span class="html-italic">δ</span>-elimination), (<b>i</b>) <span class="html-italic">K<sub>DP</sub></span> from X-PAR calculated by Exp4 (LP with <span class="html-italic">δ</span>-elimination).</p>
Full article ">Figure 9
<p>Based on the <span class="html-italic">δ</span>-elimination processing of <span class="html-italic">Φ<sub>DP</sub></span>, this figure shows the distribution of (<b>a</b>) <span class="html-italic">K<sub>DPS</sub></span> from S-POL vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LS, (<b>b</b>) <span class="html-italic">K<sub>DPS</sub></span> from S-POL vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LP, (<b>c</b>) <span class="html-italic">K<sub>DPM</sub></span> calculated using fitting relation of polarimetric variables vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LS, (<b>d</b>) <span class="html-italic">K<sub>DPM</sub></span> calculated using fitting relation of polarimetric variables vs. <span class="html-italic">K<sub>DPX</sub></span> calculated by LP.</p>
Full article ">Figure 10
<p>Colored density plot of 1 h rainfall by rain gauges vs. (<b>a</b>) QPE by EXP1 <span class="html-italic">K<sub>DP</sub></span> (LS without δ-elimination) and (<b>b</b>) QPE by EXP4 <span class="html-italic">K<sub>DP</sub></span> (LP with <span class="html-italic">δ</span>-elimination).</p>
Full article ">
22 pages, 5580 KiB  
Article
Study on Attenuation Correction for the Reflectivity of X-Band Dual-Polarization Phased-Array Weather Radar Based on a Network with S-Band Weather Radar
by Fei Geng and Liping Liu
Remote Sens. 2023, 15(5), 1333; https://doi.org/10.3390/rs15051333 - 27 Feb 2023
Cited by 1 | Viewed by 1776
Abstract
X-band dual-polarization phased-array weather radars (X-PARs) have been used in South China extensively. Eliminating the attenuation and system bias of X-band radar data is the key to utilizing the advantage of X-PAR networks. In this paper, the disdrometer raindrop-size distribution (DSD) measurements are [...] Read more.
X-band dual-polarization phased-array weather radars (X-PARs) have been used in South China extensively. Eliminating the attenuation and system bias of X-band radar data is the key to utilizing the advantage of X-PAR networks. In this paper, the disdrometer raindrop-size distribution (DSD) measurements are used to calculate the radar polarimetric variables and analyze the characteristics of precipitation attenuation. Furthermore, based on the network of S-band dual-polarization Doppler weather radar (S-POL) and X-PARs, an attenuation-correction method for X-PAR reflectivity is proposed with S-POL constraints in view of the radar-mosaic requirements of a multi-radar network. Linear programming is used to calculate the attenuation-correction parameters of different rainfall areas, which realizes the attenuation correction for X-PAR. The results show that the attenuation-correction parameters simulated based on the disdrometer DSD vary with different precipitation classification; the attenuation-corrected reflectivity of X-PARs is consistent with S-POL and can realize a more precise observation of the evolution of the convective system. Compared with previous attenuation-correction methods with constant correction parameters, the improved method can reduce the deviation between X-PAR reflectivity and that of S-POL in heavy rainfall areas and areas of strong attenuation. The method proposed in this paper is stable and effective. After effective quality control, it is found that the X-PAR network deployed in South China observes data accurately and is consistent with S-POL; thus, it is expected to achieve high temporal–spatial resolution within a radar mosaic. Full article
Show Figures

Figure 1

Figure 1
<p>Location and Digital Elevation Model (DEM) information of the S-POL and two X-PARs. The center of the red and black circles shows the sites of the S-POL and X-PARs, respectively. The distance of the black concentric circle is 10 km, and the distance of the red concentric circle is 50 km.</p>
Full article ">Figure 2
<p>Calculation of <span class="html-italic">A<sub>H</sub></span> for X-PAR attenuation correction based on networked S-band radar and precipitation classification.</p>
Full article ">Figure 3
<p>Eight-gate interpolation diagram, where <span class="html-italic">P</span> is the X-PAR observation gate, (<span class="html-italic">a</span>, <span class="html-italic">e</span>, <span class="html-italic">r</span>) is the S-POL radar coordinate of <span class="html-italic">P</span>, and <math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mn>1</mn> </msub> </mrow> </semantics></math>–<math display="inline"><semantics> <mrow> <msub> <mi>f</mi> <mn>8</mn> </msub> </mrow> </semantics></math> are the eight gates adjacent to <span class="html-italic">P</span> in S-POL observation.</p>
Full article ">Figure 4
<p>Exponential-fitting results of <span class="html-italic">Z<sub>H</sub></span> and <span class="html-italic">A<sub>H</sub></span> (<b>a</b>); linear relationship (blue line) and exponential relationship (red line) between <span class="html-italic">K<sub>DP</sub></span> and <span class="html-italic">A<sub>H</sub></span> (<b>b</b>); probability distribution of the <span class="html-italic">γ</span> value calculated according to the linear relationship for <span class="html-italic">Z<sub>H</sub></span> &lt; 45 dBZ (<b>c</b>) and for <span class="html-italic">Z<sub>H</sub></span> &gt; 45 dBZ (<b>d</b>).</p>
Full article ">Figure 5
<p>Simulated polarimetric variables and <span class="html-italic">A<sub>H</sub></span> based on DSD measurement in typical rainfall process.</p>
Full article ">Figure 6
<p>Simulated reflectivity distribution at X-band (<span class="html-italic">Z<sub>HX</sub></span>) and S-band (<span class="html-italic">Z<sub>HS</sub></span>) using DSD measurement data.</p>
Full article ">Figure 7
<p>Precipitation classification (<b>a</b>,<b>d</b>), <math display="inline"><semantics> <mrow> <msub> <mi>φ</mi> <mrow> <mi>D</mi> <mi>P</mi> </mrow> </msub> </mrow> </semantics></math> (<b>b</b>,<b>e</b>), and <span class="html-italic">PIA</span> (<b>c</b>,<b>f</b>) of X-PAR1 at 4.5° elevation-angle PPI from 01:44 to 01:50, where (<b>a</b>–<b>c</b>) are the calculation results at 01:44, (<b>d</b>–<b>f</b>) are the calculation results at 01:50, (<b>c</b>) is the <span class="html-italic">PIA</span><sub>0</sub> calculated by Equation (16), and (<b>f</b>) is the <span class="html-italic">PIA<sub>t</sub></span> calculated by attenuation-correction method proposed in this paper.</p>
Full article ">Figure 8
<p>Reflectivity measured by X-PAR and S-POL from 01:44 to 01:50 (UTC) on 5 June 2020. (<b>a</b>–<b>e</b>) are <span class="html-italic">Z<sub>H</sub></span> of X-PAR without attenuation correction at 01:44:32, 01:46: 04, 01:47:38, 01:49:10, and 01:50:42; (<b>f</b>–<b>j</b>) are <span class="html-italic">Z<sub>H</sub></span> of X-PAR after attenuation correction; and (<b>k</b>,<b>l</b>) are <span class="html-italic">Z<sub>SX</sub></span> (<span class="html-italic">Z<sub>H</sub></span> from S-POL and convert to X-PAR) at 01:42 and 01:48, respectively; A, B and C are three rapidly evolving cells in the precipitation system.</p>
Full article ">Figure 9
<p>Frequency plots of reflectivity from S-POL (after wavelength conversion and system-bias correction) and X-PAR before (<b>a</b>) and after (<b>b</b>) attenuation correction and frequency plots of the gates with <span class="html-italic">φ<sub>DP</sub></span> &gt; 40° before (<b>c</b>) and after (<b>d</b>) attenuation correction at 01:50. <span class="html-italic">n/N</span> is the frequency of the X-PAR and S-POL at the corresponding reflectivity.</p>
Full article ">Figure 10
<p>Frequency plots of Z<span class="html-italic"><sub>H</sub></span>, and <span class="html-italic">A<sub>H</sub></span> before (<b>a</b>) and after (<b>b</b>) attenuation correction at 01:50 and frequency plots of gates with <span class="html-italic">φ<sub>DP</sub></span> &gt; 40° before (<b>c</b>) and after (<b>d</b>) attenuation correction, where the gray curve is the fitting results of DSD.</p>
Full article ">
26 pages, 3976 KiB  
Review
The Retrieval of Drop Size Distribution Parameters Using a Dual-Polarimetric Radar
by GyuWon Lee, Viswanathan Bringi and Merhala Thurai
Remote Sens. 2023, 15(4), 1063; https://doi.org/10.3390/rs15041063 - 15 Feb 2023
Cited by 2 | Viewed by 2813
Abstract
The raindrop size distribution (DSD) is vital for applications such as quantitative precipitation estimation, understanding microphysical processes, and validation/improvement of two-moment bulk microphysical schemes. We trace the history of the DSD representation and its linkage to polarimetric radar observables from functional forms (exponential, [...] Read more.
The raindrop size distribution (DSD) is vital for applications such as quantitative precipitation estimation, understanding microphysical processes, and validation/improvement of two-moment bulk microphysical schemes. We trace the history of the DSD representation and its linkage to polarimetric radar observables from functional forms (exponential, gamma, and generalized gamma models) and its normalization (un-normalized, single/double-moment scaling normalized). The four-parameter generalized gamma model is a good candidate for the optimal representation of the DSD variability. A radar-based disdrometer was found to describe the five archetypical shapes (from Montreal, Canada) consisting of drizzle, the larger precipitation drops and the ‘S’-shaped curvature that occurs frequently in between the drizzle and the larger-sized precipitation. Similar ‘S’-shaped DSDs were reproduced by combining the disdrometric measurements of small-sized drops from an optical array probe and large-sized drops from 2DVD. A unified theory based on the double-moment scaling normalization is described. The theory assumes the multiple power law among moments and DSDs are scaling normalized by the two characteristic parameters which are expressed as a combination of any two moments. The normalized DSDs are remarkably stable. Thus, the mean underlying shape is fitted to the generalized gamma model from which the ‘optimized’ two shape parameters are obtained. The other moments of the distribution are obtained as the product of power laws of the reference moments M3 and M6 along with the two shape parameters. These reference moments can be from dual-polarimetric measurements: M6 from the attenuation-corrected reflectivity and M3 from attenuation-corrected differential reflectivity and the specific differential propagation phase. Thus, all the moments of the distribution can be calculated, and the microphysical evolution of the DSD can be inferred. This is one of the major findings of this article. Full article
(This article belongs to the Special Issue Radar-Based Studies of Precipitation Systems and Their Microphysics)
Show Figures

Figure 1

Figure 1
<p>60 one min observed DSDs by the Precipitation Occurrence Sensor System (POSS) on 5 May 1998 in Montreal. The x-axis is diameter in linear scale and the y-axis is number concentration in logarithmic scale.</p>
Full article ">Figure 2
<p>Average DSDs (dotted lines) for different rainfall intensities and exponential fits (solid lines) introduced by [<a href="#B11-remotesensing-15-01063" class="html-bibr">11</a>]. Published 1948 by the American Meteorological Society.</p>
Full article ">Figure 3
<p>(<b>a</b>,<b>b</b>) Sample DSDs of different shapes that are generated by the generalized gamma distribution. The liquid water content (LWC) and characteristic diameter (<span class="html-italic">D<sub>m</sub><sup>’</sup></span>) are fixed as 0.5 g m<sup>−3</sup> and 1.5 mm. Adapted from [<a href="#B31-remotesensing-15-01063" class="html-bibr">31</a>]. © American Meteorological Society. Used with permission.</p>
Full article ">Figure 4
<p>(<b>a</b>) Samples of observed one-minute DSDs from POSS in Montreal. (<b>b</b>) Their moment values with moment order. Five DSDs were highlighted by different colors. Adapted from [<a href="#B43-remotesensing-15-01063" class="html-bibr">43</a>]. © American Meteorological Society. Used with permission.</p>
Full article ">Figure 5
<p>Frequency distribution of normalized DSDs in (<b>a</b>) South Korea and (<b>b</b>) Oklahoma. (<b>c</b>) Average <span class="html-italic">h</span>(<span class="html-italic">x</span>)s and (<b>d</b>) fitted <span class="html-italic">h</span>(<span class="html-italic">x</span>) with the generalized gamma function. Adapted from [<a href="#B33-remotesensing-15-01063" class="html-bibr">33</a>]. © American Meteorological Society. Used with permission.</p>
Full article ">Figure 6
<p>Hourly DSDs from the 2DVD (plus signs) and 2DVD–MPS combined (open diamonds) adapted from Figure 11 of [<a href="#B52-remotesensing-15-01063" class="html-bibr">52</a>]. DSDs show the so-called “S-shape” with higher number concentration in smaller and medium sizes. © American Meteorological Society. Used with permission.</p>
Full article ">Figure 7
<p><span class="html-italic">D<sub>m</sub></span> from the full DSD spectra at the ground instrumentation site vs. (<b>a</b>) the measured <span class="html-italic">Z<sub>dr</sub></span> and (<b>b</b>) the measured <span class="html-italic">Z<sub>h</sub></span> from the CSU CHILL S-band radar over the disdrometer location, during a light rain event period on 17 April 2015 in Greeley, Colorado.</p>
Full article ">Figure 8
<p>(<b>a</b>) Reflectivity from the NPOL gridded data for the 30 April 2020 event at 21:06 UTC; (<b>b</b>) DSD-based rain type classification (orange: convective, cyan: stratiform, purple: mixed); (<b>c</b>) Texture-based classification (red: convective, cyan: stratiform); (<b>d</b>) matched and mismatched rain-types (see text for details), with the ground instrumentation site marked as WFF. Adapted from [<a href="#B93-remotesensing-15-01063" class="html-bibr">93</a>]. © American Meteorological Society. Used with permission.</p>
Full article ">Figure 9
<p>(<b>a</b>) Attenuation corrected <span class="html-italic">Z<sub>h</sub></span>, (<b>b</b>) attenuation corrected <span class="html-italic">Z<sub>dr</sub>,</span> and (<b>c</b>) A<sub>h</sub> for a warm rain event which occurred on 10 August 2020; (<b>d</b>) retrieved <span class="html-italic">M<sub>6</sub></span>; (<b>e</b>) retrieved <span class="html-italic">M<sub>3</sub></span>; (<b>f</b>) DSDs derived from the retrieved moments averaged over 17–18 km range interval and averaged over various height intervals centered around the following: (purple—3.9 km; blue—3.5 km; cyan—3.1 km; yellow–green: 2.84 km; orange: 2.6 km; red—2.2 km; green—1.97 km; light-green: 1.69 km; black—1.29 km).</p>
Full article ">Figure 10
<p>Block diagram outlining the DSD estimation procedure and comparing with radar data.</p>
Full article ">Figure 11
<p>Vertical profiles of radar data compared with those recalculated using DSDs derived from the retrieved moments: From the left to right, <span class="html-italic">Z<sub>h</sub></span>, <span class="html-italic">Z<sub>dr</sub></span>, <span class="html-italic">A<sub>h</sub>,</span> and <span class="html-italic">K<sub>dp</sub></span> in the 17.7 to 17.8 km radar-range.</p>
Full article ">
17 pages, 14993 KiB  
Article
Study on the Quantitative Precipitation Estimation of X-Band Dual-Polarization Phased Array Radar from Specific Differential Phase
by Guo Zhao, Hao Huang, Ye Yu, Kun Zhao, Zhengwei Yang, Gang Chen and Yu Zhang
Remote Sens. 2023, 15(2), 359; https://doi.org/10.3390/rs15020359 - 6 Jan 2023
Cited by 4 | Viewed by 2438
Abstract
In this study, the quantitative precipitation estimation (QPE) capability of three X-band dual-polarization phased array radars (PAR) in Guangzhou, South China, was demonstrated, with an S-band operational dual-polarization radar as the benchmark. Rainfall rate (R) estimators based on the specific differential [...] Read more.
In this study, the quantitative precipitation estimation (QPE) capability of three X-band dual-polarization phased array radars (PAR) in Guangzhou, South China, was demonstrated, with an S-band operational dual-polarization radar as the benchmark. Rainfall rate (R) estimators based on the specific differential phase (KDP) for summer precipitation for both X-band and S-band radars were derived from the raindrop size distributions (DSDs) observed by a 2-dimensional video disdrometer (2DVD) in South China. Rainfall estimates from the radars were evaluated with gauge observations in three events, including pre-summer rainfall, typhoon precipitation, and local severe convective precipitation. Observational results showed that radar echoes from the X-band PARs suffered much more severely from attenuation than those from the S-band radar. Compared to S-band observations, the X-band echoes can disappear when the signal-to-noise ratio drops to a certain level due to severe attenuation, resulting in different estimated rainfall areas for X- and S-band radars. The attenuation corrected by KDP had good consistency with S-band observations, but the accuracy of attenuation correction was affected by DSD uncertainty and may vary in different types of precipitation. The QPE results demonstrated that the R(KDP) estimator produced better rainfall accumulations from the X-band PAR observations compared to the S-band observations. For both the X-band and S-band radars, the estimates of hourly accumulated rainfall became more accurate in heavier rainfall, due to the decreases of both the DSD uncertainty and the impact of measurement errors. In the heavy precipitation area, the estimation accuracy of the X-band radar was high, and the overestimation of the S-band radar was obvious. Through the analysis of the ZH-ZDR distribution in the three weather events, it was found that the X-band PAR with the capability of high spatiotemporal observations can capture minute-level changes in the microphysical characteristics, which help improve the estimation accuracy of ground rainfall. Full article
Show Figures

Figure 1

Figure 1
<p>Location of X-band dual-polarization phased array radar and S-band dual-polarization radar at Guangzhou on the digital terrain elevation map.</p>
Full article ">Figure 2
<p><span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) relationship for X-band (<b>a</b>) and S-band (<b>b</b>) radars, derived by the least squares method of the observations.</p>
Full article ">Figure 3
<p>A comparison between the <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) and the intrinsic rainfall rate <span class="html-italic">R</span><sub>cal</sub> directly observed by the 2DVD for X-band (<b>a</b>) and S-band (<b>b</b>).</p>
Full article ">Figure 4
<p>The <span class="html-italic">K</span><sub>DP</sub>-<span class="html-italic">A<sub>H</sub></span> (<b>a</b>) and <span class="html-italic">K</span><sub>DP</sub>-<span class="html-italic">A</span><sub>DP</sub> (<b>b</b>) relationships for the X-band PAR fitted by DSD data.</p>
Full article ">Figure 5
<p>Measurements of <span class="html-italic">Z<sub>H</sub></span> from the Panyu PAR (<b>a</b>) without attenuation correction and (<b>b</b>) after attenuation correction at the 6.3 deg elevation at 19:30 on 2 May 2018; (<b>c</b>) Measurements of <span class="html-italic">Z<sub>H</sub></span> from the Guangzhou S-band radar at the 6 deg elevation at 19:30 on 2 May 2018.</p>
Full article ">Figure 6
<p>The <span class="html-italic">Z<sub>H</sub></span> of the X-band PAR before (blue) and after (red) attenuation correction at an azimuth of 5 deg and an elevation of 6.3 deg corresponding to <a href="#remotesensing-15-00359-f005" class="html-fig">Figure 5</a>. The black line is from the S-band radar, also at an azimuth of 5 deg and an elevation of 6 deg.</p>
Full article ">Figure 7
<p>The NDFs of estimated hourly accumulated rainfall from (<b>a</b>) X-band <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) and (<b>b</b>) S-band <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) versus the gauge observations for all precipitation cases.</p>
Full article ">Figure 8
<p><span class="html-italic">Z<sub>H</sub></span>, <span class="html-italic">Z</span><sub>DR</sub>, and <span class="html-italic">K</span><sub>DP</sub> observed by the X-band PAR radar at the 2.7-deg-elevation in the (<b>a</b>–<b>c</b>) pre-summer rainband at 10:30 on June 7, the (<b>d</b>–<b>f</b>) typhoon rainband at 13:30 on August 19, and the (<b>g</b>–<b>i</b>) severe convective rainband at 17:30 on September 4.</p>
Full article ">Figure 9
<p>Comparison of hourly accumulated rainfall obtained using <span class="html-italic">R</span>(<span class="html-italic">K</span><sub>DP</sub>) with gauge observations in the pre-summer rainfall between 10:00–11:00 on June 7, the typhoon precipitation between 13:00–14:00 on August 19, and the severe convective precipitation rainband between 17:00–18:00 on September 4. (<b>a</b>,<b>c</b>,<b>e</b>) for X-band band radar and (<b>b</b>,<b>d</b>,<b>f</b>) for S-band band radar. The color of the bubble charts shows bias ratios (<math display="inline"><semantics> <mrow> <msub> <mi>R</mi> <mrow> <msub> <mi>K</mi> <mrow> <mi>DP</mi> </mrow> </msub> </mrow> </msub> </mrow> </semantics></math>/<span class="html-italic">R</span><sub>gauge</sub>) between the hourly accumulated rainfall estimates and independent gauge observations.</p>
Full article ">Figure 10
<p>The normalized density of <span class="html-italic">Z<sub>H</sub></span>-<span class="html-italic">Z</span><sub>DR.</sub> in the pre-summer rainband observed by the Huadu X-band PAR at (<b>a</b>) 10:48, (<b>b</b>) 10:50, (<b>c</b>) 10:52, and (<b>d</b>) the S-band radar at 10:54 on 7 June 2020.</p>
Full article ">Figure 11
<p>The normalized density of <span class="html-italic">Z<sub>H</sub></span>-<span class="html-italic">Z</span><sub>DR.</sub> in the Typhoon Higos rainband observed by the Nansha X-band PAR at (<b>a</b>) 13:00, (<b>b</b>) 13:02, (<b>c</b>) 13:04, and (<b>d</b>) the S-band radar at 13:00 on 19 August 2020.</p>
Full article ">Figure 12
<p>The normalized density of <span class="html-italic">Z<sub>H</sub></span>-<span class="html-italic">Z</span><sub>DR</sub> in the local heavy precipitation rainband observed by the Panyu X-band PAR at (<b>a</b>) 17:30, (<b>b</b>) 17:31, (<b>c</b>) 17:33, and (<b>d</b>) the S-band radar at 17:30 on 4 September 2020.</p>
Full article ">
19 pages, 3141 KiB  
Article
The Characteristics of Raindrop Size Distributions in Different Climatological Regions in South Korea
by Cheol-Hwan You, Hyeon-Joon Kim, Sung-Ho Suh, Woonseon Jung and Mi-Young Kang
Remote Sens. 2022, 14(20), 5137; https://doi.org/10.3390/rs14205137 - 14 Oct 2022
Viewed by 1621
Abstract
To understand the microphysical characteristics of rainfall in four different climatological regions (called BOS, BUS, CPO, and JIN) in South Korea, DSDs and their variables, including the mass-weighted mean diameter (Dm) and normalized number concentration (logNw), were examined. To [...] Read more.
To understand the microphysical characteristics of rainfall in four different climatological regions (called BOS, BUS, CPO, and JIN) in South Korea, DSDs and their variables, including the mass-weighted mean diameter (Dm) and normalized number concentration (logNw), were examined. To examine the characteristics of DSDs at four sites with different climatology and topography, data measured from Parsivel disdrometer and wind direction from Automatic Weather System (AWS) during rainy seasons from June to August for three years (2018 to 2020) were analyzed. The DSDs variables were calculated using Gamma distribution model. In the coastal area, larger raindrops with a lower number concentration occurred, whereas smaller raindrops with a higher number concentration dominated in the middle land and mountain region. The mountain area of CPO and middle land area of JIN had a larger contribution to the rain rate than that of the coastal area of BOS and JIN in the range of the smallest diameter. The contribution of the drop size to the total number concentration at the CPO and JIN sites was larger (smaller) than that at BOS and BUS in the smallest (larger) diameter. The average shape and slope parameter of gamma model were higher values at the mountain area than at other sites for both rain types, Z-R relation and polarimetric variables were also shown different values at the four studied sites. The intercept coefficient of Z-R relation showed higher values in the mountain area and middle land area than the coastal area. The slope values of Z-R relation were the smallest in the mountain area. The polarimetric variables of ZH and ZDR were shown highest (lowest) value at the coastal region of BOS (mountain area of CPO) site for both rain types. The Dm-rose, which shows the Dm distributions with the wind direction, was used in this study. In the coastal area (mountain and middle land area), the dominant wind was east–southeast (east) direction. The ratio of the smaller diameter to the middle size at BOS was much smaller than that at CPO. In the analysis of the hourly distribution of the Dm and logNw, there were two and four peaks of Dm at BUS and BOS, respectively. There was one peak of the Dm at the CPO and JIN sites. The time variation of the Dm was much higher than that of the logNw. Full article
(This article belongs to the Section Atmospheric Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The locations of the four PARSIVEL disdrometers (full circle), with topography. CPO site is located at the mountain area (849.15 m in altitude), JIN is a middle land area (111 m in altitude), BUS (1.41 m in altitude) and BOS (10 m in altitude) in a coastal area.</p>
Full article ">Figure 2
<p>The time series of annual rainfall (blue) and the trend (red) for 30 years at (<b>a</b>) CPO, (<b>b</b>) JIN, (<b>c</b>) BOS, and (<b>d</b>) BUS site. Cheonan is the nearest official station to JIN, Goheung is the nearest site to BOS, and Pusan is the nearest site to BUS.</p>
Full article ">Figure 3
<p>The workflow of this study.</p>
Full article ">Figure 4
<p>Average raindrop number concentration with diameter at four different sites. The rectangles in blue represent BOS, the triangles in red represent BUS, the circles in purple represent CPO, and the arrows in green represent the JIN site. The numbers in the legend are sampling numbers at each site.</p>
Full article ">Figure 5
<p>Average raindrop number concentration with diameter at four different sites with respect to rain rate (<b>a</b>) 0.1 &lt; R &lt; 5.0 mm h<sup>−1</sup>, (<b>b</b>) 5.0 &lt; R &lt; 10.0 mm h<sup>−1</sup>, and (<b>c</b>) 10.0 mm h<sup>−1</sup> &lt; R. The rectangles in blue represent BOS, the triangles in red represent BUS, the circles in purple represent CPO, and the arrows in green represent the JIN site. The numbers in the legend represent the percentage of each rain rate.</p>
Full article ">Figure 6
<p>The contribution ratios of the raindrop diameter to (<b>a</b>) the rain rate and (<b>b</b>) the total number concentration. The blue represents BOS site, the red represents BUS site, the purple represents CPO site, and the green represents JIN site.</p>
Full article ">Figure 7
<p>Average number concentrations of diameter with rain types (<b>a</b>) stratiform and (<b>b</b>) convective rain at the four studied sites. The number in the legend represent the sample number. The rectangles in blue represent BOS, the triangles in red represent BUS, the circles in purple represent CPO, and the arrows in green represent the JIN site.</p>
Full article ">Figure 8
<p>The occurrence frequency of shape parameters with rain type at (<b>a</b>) BOS, (<b>b</b>) BUS, (<b>c</b>) CPO, and (<b>d</b>) JIN site. The blue represents stratiform rain and the red represents convective rain. The legend shows the average (Avg.), standard deviation (Std.), and skewness (Skew) of shape parameter.</p>
Full article ">Figure 9
<p>The occurrence frequency of slope parameters with rain type at (<b>a</b>) BOS, (<b>b</b>) BUS, (<b>c</b>) CPO, and (<b>d</b>) JIN site. The blue represents stratiform rain and the red represents convective rain. The legend shows the average (Avg.), standard deviation (Std.), and skewness (Skew) of slope parameter.</p>
Full article ">Figure 10
<p>Boxplot of <span class="html-italic">Dm</span> and <span class="html-italic">logNw</span> values of convective and stratiform rain for each site. The line curve represents the separation line of rain types proposed by BR09. The section above (below) this line represents convective (stratiform) rain.</p>
Full article ">Figure 11
<p>The scatter plots of rain rate measured by DSD and R(Z) relations were obtained from DSDs at (<b>a</b>) BOS, (<b>b</b>) BUS, (<b>c</b>) CPO, and (<b>d</b>) JIN sites. The legend shows the Z-R relation, root mean square error (RMSE), and cross-correlation coefficient (CC).</p>
Full article ">Figure 12
<p>The occurrence frequency of reflectivity with rain types at (<b>a</b>) BOS, (<b>b</b>) BUS, (<b>c</b>) CPO, and (<b>d</b>) JIN site. The blue represents stratiform rain and the red represents convective rain. The legend shows the average (Avg.), standard deviation (Std.), and skewness (Skew) of reflectivity for both rain types.</p>
Full article ">Figure 13
<p>The occurrence frequency of differential reflectivity with rain types at (<b>a</b>) BOS, (<b>b</b>) BUS, (<b>c</b>) CPO, and (<b>d</b>) JIN site. The blue represents stratiform rain and the red represents convective rain. The legend shows the average (Avg.), standard deviation (Std.), and skewness (Skew) of differential reflectivity for both rain types.</p>
Full article ">Figure 14
<p>The occurrence frequency of specific differential phase with rain types at (<b>a</b>) BOS, (<b>b</b>) BUS, (<b>c</b>) CPO, and (<b>d</b>) JIN site. The blue represents stratiform rain and the red represents convective rain. The legend shows the average (Avg.), standard deviation (Std.), and skewness (Skew) of specific differential phase for both rain types.</p>
Full article ">Figure 15
<p>The <span class="html-italic">Dm</span> distribution with wind direction of the (<b>a</b>) BOS, (<b>b</b>) BUS, (<b>c</b>) CPO, and (<b>d</b>) JIN sites.</p>
Full article ">Figure 16
<p>Time series of average (<b>a</b>) <span class="html-italic">Dm</span> and (<b>b</b>) <span class="html-italic">logNw</span> at the four studied sites. The rectangles in blue represent BOS, the triangles in red represent BUS, the circles in purple represent CPO, and the arrows in green represent the JIN site.</p>
Full article ">
Back to TopTop