scholarly journals From model to radar variables: a new forward polarimetric radar operator for COSMO

2018 ◽  
Vol 11 (7) ◽  
pp. 3883-3916 ◽  
Author(s):  
Daniel Wolfensberger ◽  
Alexis Berne

Abstract. In this work, a new forward polarimetric radar operator for the COSMO numerical weather prediction (NWP) model is proposed. This operator is able to simulate measurements of radar reflectivity at horizontal polarization, differential reflectivity as well as specific differential phase shift and Doppler variables for ground based or spaceborne radar scans from atmospheric conditions simulated by COSMO. The operator includes a new Doppler scheme, which allows estimation of the full Doppler spectrum, as well a melting scheme which allows representing the very specific polarimetric signature of melting hydrometeors. In addition, the operator is adapted to both the operational one-moment microphysical scheme of COSMO and its more advanced two-moment scheme. The parameters of the relationships between the microphysical and scattering properties of the various hydrometeors are derived either from the literature or, in the case of graupel and aggregates, from observations collected in Switzerland. The operator is evaluated by comparing the simulated fields of radar observables with observations from the Swiss operational radar network, from a high resolution X-band research radar and from the dual-frequency precipitation radar of the Global Precipitation Measurement satellite (GPM-DPR). This evaluation shows that the operator is able to simulate an accurate Doppler spectrum and accurate radial velocities as well as realistic distributions of polarimetric variables in the liquid phase. In the solid phase, the simulated reflectivities agree relatively well with radar observations, but the simulated differential reflectivity and specific differential phase shift upon propagation tend to be underestimated. This radar operator makes it possible to compare directly radar observations from various sources with COSMO simulations and as such is a valuable tool to evaluate and test the microphysical parameterizations of the model.

2017 ◽  
Author(s):  
Daniel Wolfensberger ◽  
Alexis Berne

Abstract. In this work, a new forward polarimetric radar operator for the COSMO numerical weather prediction (NWP) model is proposed. This operator is able to simulate measurements of radar reflectivity at horizontal polarization, differential reflectivity as well as specific differential phase shift and Doppler variables for ground based or spaceborne radar scans from atmospheric conditions simulated by COSMO. The operator includes a new Doppler scheme, which allows to estimate the full Doppler spectrum, as well a melting scheme which allows to represent the very specific polarimetric signature of melting hydrometeors. In addition, the operator is adapted to both the operational one-moment microphysical scheme of COSMO and its more advanced two-moment scheme. The parameters of the relationships between the microphysical and scattering properties of the various hydrometeors are derived either from the literature or, in the case of graupel and aggregates, from observations collected in Switzerland. The operator is evaluated by comparing the simulated fields of radar observables with observations from the Swiss operational radar network, from a high resolution X-band research radar and from the dual-frequency precipitation radar of the Global Precipitation Measurement satellite (GPM-DPR). This evaluation shows that the operator is able to simulate an accurate Doppler spectrum and accurate radial velocities as well as realistic distributions of polarimetric variables in the liquid phase. In the solid phase, the simulated reflectivities agree relatively well with radar observations, but the simulated differential reflectivity and specific differential phase shift upon propagation tend to be underestimated. This radar operator makes it possible to compare directly radar observations from various sources with COSMO simulations and as such is a valuable tool to evaluate and test the microphysical parameterizations of the model.


2015 ◽  
Vol 8 (11) ◽  
pp. 4681-4698 ◽  
Author(s):  
G. Vulpiani ◽  
L. Baldini ◽  
N. Roberto

Abstract. This work documents the effective use of X-band radar observations for monitoring severe storms in an operational framework. Two severe hail-bearing Mediterranean storms that occurred in 2013 in southern Italy, flooding two important Sicilian cities, are described in terms of their polarimetric radar signatures and retrieved rainfall fields. The X-band dual-polarization radar operating inside the Catania airport (Sicily, Italy), managed by the Italian Department of Civil Protection, is considered here. A suitable processing is applied to X-band radar measurements. The crucial procedural step relies on the differential phase processing, being preparatory for attenuation correction and rainfall estimation. It is based on an iterative approach that uses a very short-length (1 km) moving window, allowing proper capture of the observed high radial gradients of the differential phase. The parameterization of the attenuation correction algorithm, which uses the reconstructed differential phase shift, is derived from electromagnetic simulations based on 3 years of drop size distribution (DSD) observations collected in Rome (Italy). A fuzzy logic hydrometeor classification algorithm was also adopted to support the analysis of the storm characteristics. The precipitation field amounts were reconstructed using a combined polarimetric rainfall algorithm based on reflectivity and specific differential phase. The first storm was observed on 21 February when a winter convective system that originated in the Tyrrhenian Sea, marginally hit the central-eastern coastline of Sicily, causing a flash flood in Catania. Due to an optimal location (the system is located a few kilometers from the city center), it was possible to retrieve the storm characteristics fairly well, including the amount of rainfall field at the ground. Extemporaneous signal extinction, caused by close-range hail core causing significant differential phase shift in a very short-range path, is documented. The second storm, on 21 August 2013, was a summer mesoscale convective system that originated from a Mediterranean low pressure system lasting a few hours that eventually flooded the city of Syracuse. The undergoing physical process, including the storm dynamics, is inferred by analyzing the vertical sections of the polarimetric radar measurements. The high registered amount of precipitation was fairly well reconstructed, although with a trend toward underestimation at increasing distances. Several episodes of signal extinction were clearly manifested during the mature stage of the observed supercells.


2021 ◽  
Author(s):  
Gregor Möller ◽  
Florian Ewald ◽  
Silke Groß ◽  
Martin Hagen ◽  
Christoph Knote ◽  
...  

<p>The representation of microphysical processes in numerical weather prediction models remains a main source of uncertainty. To tackle this issue, we exploit the synergy of two polarimetric radars to provide novel observations of model microphysics parameterizations. In the framework of the IcePolCKa project (Investigation of the initiation of Convection and the Evolution of Precipitationusing simulatiOns and poLarimetric radar observations at C- and Ka-band) we use these observations to study the initiation of convection as well as the evolution of precipitation. At a distance of 23 km between the C-band PoldiRad radar of the German Aerospace Center (DLR) in Oberpfaffenhofen and the Ka-band Mira35 radar of the Ludwig-Maximilians-University of Munich (LMU), the two radar systems allow targeted observations and coordinated scan patterns. A second C-band radar located in Isen and operated by the German Weather Service (DWD) provides area coverage and larger spatial context. By tracking the precipitation movement, the dual-frequency and polarimetric radar observations allow us to characterize important microphysical parameters, such as predominant hydrometeor class or conversion rates between these classes over a significant fraction of the life time of a convective cell. A WRF (Weather Research and Forecasting Model) simulation setup has been established including a Europe-, a nested Germany- and a nested Munich- domain. The Munich domain covers the overlap area of our two radars Mira35 and Poldirad with a horizontal resolution of 400 m. For each of our measurement days we conduct a WRF hindcast simulation with differing microphysics schemes. To allow for a comparison between model world and observation space, we make use of the radar forward-simulator CR-SIM. The measurements so far include 240 coordinated scans of 36 different convective cells over 10 measurement days between end of April and mid July 2019 as well as 40 days of general dual-frequency volume scans between mid April and early October 2020. The performance of each microphysics scheme is analyzed through a comparison to our radar measurements on a statistical basis over all our measurements.</p>


2014 ◽  
Vol 2014 ◽  
pp. 1-10 ◽  
Author(s):  
Cheol-Hwan You ◽  
Dong-In Lee ◽  
Mi-Young Kang

To assess the performance of rainfall estimation using specific differential phase observed by Bislsan radar, the first polarimetric radar in Korea, three rainfall cases occurring in 2011 were selected, each caused by different conditions: the first is the Changma front and typhoon, the second is only the Changma front, and the third is only a typhoon. For quantitative use of specific differential phase (KDP), a data quality algorithm was developed for differential phase shift (ΦDP), composed of two steps; the first involves removal of scattered noise and the second is unfolding ofΦDP. This order of the algorithm is necessary so as not to remove unfolded areas, which are the real meteorological target. All noise was removed and the foldedΦDPwere unfolded successfully for this study.RKDPrelations for S-band radar were calculated for 84,754 samples of observed drop size distribution (DSD) using different drop shape assumptions. The relation for the Bringi drop shape showed the best statistics: 0.28 for normalized error, and 6.7 mm for root mean square error for rainfall heavier than 10 mm h-1. Because the drop shape assumption affects the accuracy of rainfall estimation differently for different rainfall types, such characteristics should be taken into account to estimate rainfall more accurately using polarimetric variables.


2005 ◽  
Vol 6 (3) ◽  
pp. 248-262 ◽  
Author(s):  
Sergey Y. Matrosov ◽  
David E. Kingsmill ◽  
Brooks E. Martner ◽  
F. Martin Ralph

Abstract The utility of X-band polarimetric radar for quantitative retrievals of rainfall parameters is analyzed using observations collected along the U.S. west coast near the mouth of the Russian River during the Hydrometeorological Testbed project conducted by NOAA’s Environmental Technology and National Severe Storms Laboratories in December 2003 through March 2004. It is demonstrated that the rain attenuation effects in measurements of reflectivity (Ze) and differential attenuation effects in measurements of differential reflectivity (ZDR) can be efficiently corrected in near–real time using differential phase shift data. A scheme for correcting gaseous attenuation effects that are important at longer ranges is introduced. The use of polarimetric rainfall estimators that utilize specific differential phase and differential reflectivity data often provides results that are superior to estimators that use fixed reflectivity-based relations, even if these relations were derived from the ensemble of drop size distributions collected in a given geographical region. Comparisons of polarimetrically derived rainfall accumulations with data from the high-resolution rain gauges located along the coast indicated deviation between radar and gauge estimates of about 25%. The ZDR measurements corrected for differential attenuation were also used to retrieve median raindrop sizes, D0. Because of uncertainties in differential reflectivity measurements, these retrievals are typically performed only for D0 > 0.75 mm. The D0 estimates from an impact disdrometer located at 25 km from the radar were in good agreement with the radar retrievals. The experience of operating the transportable polarimetric X-band radar in the coastal area that does not have good coverage by the National Weather Service radar network showed the value of such radar in filling the gaps in the network coverage. The NOAA X-band radar was effective in covering an area up to 40–50 km in radius offshore adjacent to a region that is prone to flooding during wintertime landfalling Pacific storms.


2015 ◽  
Vol 8 (7) ◽  
pp. 7201-7237
Author(s):  
G. Vulpiani ◽  
L. Baldini ◽  
N. Roberto

Abstract. This work documents the fruitul use of X-band radar observations for the monitoring of severe storms in an operational framework. More specifically, a couple of severe hail-bearing Mediterranean storms occurred in 2013 in southern Italy, flooding two important cities of Sicily, are described in terms of their polarimetric radar signatures and retrieved rainfall fields. It is used the X-band dual-polarization radar operating inside the Catania airport (Sicily, Italy), managed by the Italian Department of Civil Protection. A suitable processing is applied to X-band radar measurements. The crucial procedural step relies on the differential phase processing based on an iterative approach that uses a very short-length (1 km) moving window allowing to properly catch the observed high radial gradients of the differential phase. The parameterization of the attenuation correction algorithm, which use the reconstructed differential phase shift, is derived from electromagnetic simulations based on 3 years of DSD observations collected in Rome (Italy). A Fuzzy Logic hydrometeor classification algorithm was also adopted to support the analysis of the storm characteristics. The precipitation fields amount were reconstructed using a combined polarimetric rainfall algorithm based on reflectivity and specific differential phase. The first considered storm was observed on the 21 February, when a winter convective system, originated in the Tyrrhenian sea, hit only marginally the central-eastern coastline of Sicily causing the flash-flood of Catania. Due to the optimal radar location (the system is located at just few kilometers from the city center), it was possible to well retrieve the storm characteristics, including the amount of rainfall field at ground. Extemporaneous signal extinction, caused by close-range hail core causing significant differential phase shift in very short range path, is documented. The second storm, occurred on 21 August 2013, is a summer mesoscale convective system originated by the temperature gradient between sea and land surface, lasted a few hours and eventually flooded the city of Siracusa. The undergoing physical process, including the storm dynamics, is inferred by analysing the vertical sections of the polarimetric radar measurements. The high registered precipitation amount was fairly well reconstructed even though with a trend to underestimation at increasing distances. Several episodes of signal extinction clearly manifested during the mature stage of the observed supercell.


2021 ◽  
Vol 13 (16) ◽  
pp. 3060
Author(s):  
Muyun Du ◽  
Jidong Gao ◽  
Guifu Zhang ◽  
Yunheng Wang ◽  
Pamela L. Heiselman ◽  
...  

Polarimetric radar data (PRD) have potential to be used in numerical weather prediction (NWP) models to improve convective-scale weather forecasts. However, thus far only a few studies have been undertaken in this research direction. To assimilate PRD in NWP models, a forward operator, also called a PRD simulator, is needed to establish the relation between model physics parameters and polarimetric radar variables. Such a forward operator needs to be accurate enough to make quantitative comparisons between radar observations and model output feasible, and to be computationally efficient so that these observations can be easily incorporated into a data assimilation (DA) scheme. To address this concern, a set of parameterized PRD simulators for the horizontal reflectivity, differential reflectivity, specific differential phase, and cross-correlation coefficient were developed. In this study, we have tested the performance of these new operators in a variational DA system. Firstly, the tangent linear and adjoint (TL/AD) models for these PRD simulators have been developed and checked for the validity. Then, both the forward operator and its adjoint model have been built into the three-dimensional variational (3DVAR) system. Finally, some preliminary DA experiments have been performed with an idealized supercell storm. It is found that the assimilation of PRD, including differential reflectivity and specific differential phase, in addition to radar radial velocity and horizontal reflectivity, can enhance the accuracy of both initial conditions for model hydrometer state variables and ensuing model forecasts. The usefulness of the cross-correlation coefficient is very limited in terms of improving convective-scale data analysis and NWP.


2009 ◽  
Vol 11 (4) ◽  
pp. 045020 ◽  
Author(s):  
Lijun Ma ◽  
S Nam ◽  
Hai Xu ◽  
B Baek ◽  
Tiejun Chang ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document