Comprehensive T-Matrix Reference Database: A 2007-2009 Update

The T-matrix method is among the most versatile, efficient, and widely used theoretical techniques for the numerically exact computation of electromagnetic scattering by homogeneous and composite particles, clusters of particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper presents an update to the comprehensive database of T-matrix publications compiled by us previously and includes the publications that appeared since 2007. It also lists several earlier publications not included in the original database

Melting Layer Attenuation at Ka- and W-Bands as Derived From Multifrequency Radar Doppler Spectra Observations

Abstract The melting layer of precipitation has a major impact on remote sensing and telecommunications. However, there is a shortage of observational studies to validate and constrain the melting layer models especially for high-frequency radar bands. In this paper, we report how multifrequency radar Doppler spectra can be used to retrieve the melting layer attenuation at Ka- and W-bands. The presented analysis is based on identifying Rayleigh scattering regions in radar Doppler spectra measurements where dual-wavelength spectral ratios can be related to differential attenuation. We show that the estimated attenuation at Ka- and W-bands agrees reasonably well with previously reported studies, but there are indications of differences at higher rain rates. We advocate that this technique can be applied to long-term observations to advance our knowledge of the melting process. The parameterizations of melting layer attenuation as a function of rain rate and radar reflectivity are also presented.Peer reviewe

Opportunistic rain rate estimation from measurements of satellite downlink attenuation: A survey

Recent years have witnessed a growing interest in techniques and systems for rainfall surveillance on regional scale, with increasingly stringent requirements in terms of the following: (i) accuracy of rainfall rate measurements, (ii) adequate density of sensors over the territory, (iii) space‐time continuity and completeness of data and (iv) capability to elaborate rainfall maps in near real time. The devices deployed to monitor the precipitation fields are traditionally networks of rain gauges distributed throughout the territory, along with weather radars and satellite remote sensors operating in the optical or infrared band, none of which, however, are suitable for full compliance to all of the requirements cited above. More recently, a different approach to rain rate estimation techniques has been proposed and investigated, based on the measurement of the attenuation induced by rain on signals of pre‐existing radio networks either in terrestrial links, e.g., the backhaul connections in cellular networks, or in satellite‐to‐earth links and, among the latter, notably those between geostationary broadcast satellites and domestic subscriber terminals in the Ku and Ka bands. Knowledge of the above rain‐induced attenuation permits the retrieval of the corresponding rain intensity provided that a number of meteorological and geometric parameters are known and ultimately permits estimating the rain rate locally at the receiver site. In this survey paper, we specifically focus on such a type of “opportunistic” systems for rain field monitoring, which appear very promising in view of the wide diffusion over the territory of low‐cost domestic terminals for the reception of satellite signals, prospectively allowing for a considerable geographical capillarity in the distribution of sensors, at least in more densely populated areas. The purpose of the paper is to present a broad albeit synthetic overview of the numerous issues inherent in the above rain monitoring approach, along with a number of solutions and algorithms proposed in the literature in recent years, and ultimately to provide an exhaustive account of the current state of the art. Initially, the main relevant aspects of the satellite link are reviewed, including those related to satellite dynamics, frequency bands, signal formats, propagation channel and radio link geometry, all of which have a role in rainfall rate estimation algorithms. We discuss the impact of all these factors on rain estimation accuracy while also highlighting the substantial differences inherent in this approach in comparison with traditional rain monitoring techniques. We also review the basic formulas relating rain rate intensity to a variation of the received signal level or of the signal‐to-noise ratio. Furthermore, we present a comprehensive literature survey of the main research issues for the aforementioned scenario and provide a brief outline of the algorithms proposed for their solution, highlighting their points of strength and weakness. The paper includes an extensive list of bibliographic references from which the material presented herein was taken

Detection and characterization of the melting layer based on polarimetric radar scans

Stratiform rain situations are generally associated with the presence of a melting layer characterized by a strong signature in polarimetric radar variables. This layer is an important feature as it indicates the transition from solid to liquid precipitation. The melting layer remains poorly characterized, particularly from a polarimetric radar point of view. In this work a new algorithm to automatically detect the melting layer on polarimetric RHI radar scans using gradients of reflectivity and copolar correlation is first proposed. The algorithm was applied to high-resolution X-band polarimetric radar data and validated by comparing the height of the detected layer with freezing-level heights obtained from radiosoundings and was shown to give both small errors and bias. The algorithm was then used on a large selection of precipitation events (more than 4000 RHI scans) from different seasons and climatic regions (South of France, Swiss Alps and plateau, and Iowa, USA) to characterize the geometric and polarimetric signatures of the melting layer. The melting layer is shown to have a very similar geometry on average, independent of the topography and climatic conditions. Variations in the thickness of the melting layer during and between precipitation events was shown to be strongly related to the presence of rimed particles, to the vertical velocity of hydrometeors and to the intensity of the bright band

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

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.</p

A polarimetric radar operator to evaluate precipitation from the COSMO atmospheric model

Weather radars provide real-time measurements of precipitation at a high temporal and spatial resolution and over a large domain. A drawback, however, it that these measurements are indirect and require careful interpretation to yield relevant information about the mechanisms of precipitation. Radar observations are an invaluable asset for the numerical forecast of precipitation, both for data assimilation, parametrization of subscale phenomena and model verification. This thesis aims at investigating new uses for polarimetric radar data in numerical weather prediction. The first part of this work is devoted to the design of an algorithm able to automatically detect the location and extent of the melting layer of precipitation , an important feature of stratiform precipitation, from vertical radar scans. This algorithm is then used to provide a detailed characterization of the melting layer, in several climatological regions, providing thus relevant information for the parameterization of melting processes and the evaluation of simulated freezing level heights. The second part of this work uses a multi-scale approach based on the multifractal framework to evaluate precipitation fields simulated by the COSMO weather model with radar observations. A climatological analysis is first conducted to relate multifractal parameters to physical descriptors of precipitation. A short-term analysis, that focuses on three precipitation events over Switzerland, is then performed. The results indicate that the COSMO simulations exhibit spatial scaling breaks that are not present in the radar data. It is also shown that a more advanced microphysics parameterization generates larger extreme values, and more discontinuous precipitation fields, which agree better with radar observations. The last part of this thesis describes a new forward polarimetric radar operator, able to simulate realistic radar variables from outputs of the COSMO model, taking into account most physical aspects of beam propagation and scattering. An efficient numerical scheme is proposed to estimate the full Doppler spectrum, a type of measurement often performed by research radars, which provides rich information about the particle velocities and turbulence. The operator is evaluated with large datasets from various ground and spaceborne radars. This evaluation shows that the operator is able to simulate accurate Doppler variables and realistic distributions of polarimetric variables in the liquid phase. In the solid phase, the simulated reflectivities agree relatively well with radar observations, but the polarimetric variables tend to be underestimated. A detailed sensitivity analysis of the radar operator reveals that, in the liquid phase, the simulated radar variables depend very much on the hypothesis about drop geometry and drop size distributions. In the solid phase, the potential of more advanced scattering techniques is investigated, revealing that these methods could help to resolve the strong underestimation of polarimetric variables in snow and graupel

The impact of precipitation physical processes on the polarimetric radar variables

The ongoing upgrade of the National Weather Service WSR-88D radar network to polarimetric capabilities, as well as similar upgrades worldwide, will soon provide a wealth of data and information regarding storm precipitation physics. Fully understanding how a variety of microphysical processes are revealed in polarimetric data is necessary for the best use of these new data by operational and research meteorologists. The focus of this study is to quantify a number of these precipitation physics "fingerprints" in the polarimetric radar variables by using a synthesis of explicit microphysical modeling, electromagnetic scattering calculations, thought experiments, and polarimetric radar observations.The complete set of polarimetric variables available from linearly-orthogonal dual-polarization radars are derived from basic electromagnetic scattering principles. A detailed physical description of these variables is then provided for applications in precipitation and other atmospheric scatterers, as is a summary of common data artifacts. The impact of various precipitation physical processes on these radar variables is then quantified. Using explicit microphysical models, scattering calculations, observations, and thought experiments, the microphysical fingerprints are determined and quantified for raindrop thermal conduction, raindrop size sorting by sedimentation, updrafts, and vertical wind shear, evaporation of raindrops, coalescence of raindrops, and freezing of raindrops in deep convective storm updrafts. A catalogue of the qualitative fingerprints of a number of precipitation processes is summarized