Snowfall retrieval at X, Ka and W bands : consistency of backscattering and microphysical properties using BAECC ground-based measurements

Radar-based snowfall intensity retrieval is investigated at centimeter and millimeter wavelengths using co-located ground-based multi-frequency radar and video-disdrometer observations. Using data from four snowfall events, recorded during the Biogenic Aerosols Effects on Clouds and Climate (BAECC) campaign in Finland, measurements of liquid-water-equivalent snowfall rate S are correlated to radar equivalent reflectivity factors Z(e), measured by the Atmospheric Radiation Measurement (ARM) cloud radars operating at X, Ka and W frequency bands. From these combined observations, power-law Z(e)-S relationships are derived for all three frequencies considering the influence of riming Using microwave radiometer observations of liquid water path, the measured precipitation is divided into lightly, moderately and heavily rimed snow. Interestingly lightly rimed snow events show a spectrally distinct signature of Z(e)-S with respect to moderately or heavily rimed snow cases. In order to understand the connection between snowflake microphysical and multi-frequency backscattering properties, numerical simulations are performed by using the particle size distribution provided by the in situ video disdrometer and retrieved ice particle masses. The latter are carried out by using both the T-matrix method (TMM) applied to soft-spheroid particle models with different aspect ratios and exploiting a pre-computed discrete dipole approximation (DDA) database for rimed aggregates. Based on the presented results, it is concluded that the soft-spheroid approximation can be adopted to explain the observed multifrequency Z(e)-S relations if a proper spheroid aspect ratio is selected. The latter may depend on the degree of riming in snowfall. A further analysis of the backscattering simulations reveals that TMM cross sections are higher than the DDA ones for small ice particles, but lower for larger particles. The differences of computed cross sections for larger and smaller particles are compensating for each other. This may explain why the soft-spheroid approximation is satisfactory for radar reflectivity simulations under study.Peer reviewe

Evaluation of convective boundary layer height estimates using radars operating at different frequency bands

Knowledge of the atmospheric boundary layer state and evolution is important for understanding air pollution and low-level cloud development, among other things. There are a number of instruments and methods that are currently used to estimate boundary layer height (BLH). However, no single instrument is capable of providing BLH measurements in all weather conditions. We proposed a method to derive a daytime convective BLH using clear air echoes in radar observations and investigated the consistency of these retrievals between different radar frequencies. We utilized data from three vertically pointing radars that are available at the SMEAR II station in Finland, i.e. the C band (5 GHz), Ka band (35 GHz) and W band (94 GHz). The Ka- or W-band cloud radars are an integral part of cloud profiling stations of pan-European Aerosol, Clouds and Trace Gases Research Infrastructure (ACTRIS). Our method will be utilized at ACTRIS stations to serve as an additional estimate of the BLH during summer months. During this period, insects and Bragg scatter are often responsible for clear air echoes recorded by weather and cloud radars. To retrieve a BLH, we suggested a mechanism to separate passive and independently flying insects that works for all analysed frequency bands. At the lower frequency (the C band) insect scattering has been separated from Bragg scattering using a combination of the radar reflectivity factor and linear depolarization ratio. Retrieved values of the BLH from all radars are in a good agreement when compared to the BLH obtained with the co-located HALO Doppler lidar and ERA5 reanalysis data set. Our method showed some underestimation of the BLH after nighttime heavy precipitation yet demonstrated a potential to serve as a reliable method to obtain a BLH during clear-sky days. Additionally, the entrainment zone was observed by the C-band radar above the CBL in the form of a Bragg scatter layer. Aircraft observations of vertical profiles of potential temperature and water vapour concentration, collected in the vicinity of the radar, demonstrated some agreement with the Bragg scatter layer.Peer reviewe

Estimation of extreme precipitation events in Estonia and Italy using dual-polarization weather radar quantitative precipitation estimations

Evaluating extreme rainfall for a certain location is commonly considered when designing stormwater management systems. Rain gauge data are widely used to estimate rainfall intensities for a given return period. However, the poor spatial and temporal resolution of operational gauges is the main limiting factor. Several studies have used rainfall estimates based on weather radar horizontal reflectivity (Zh), but they come with a great caveat: while proven reliable for low or moderate rainfall rates, they are subject to major errors in extreme rainfall and convective cases. It is widely known that C-band weather radar can underestimate precipitation intensity due to signal attenuation or overestimate it due to hail and clutter contamination. From the late 1990s, dual-polarization weather radar started to become operational in the national surveillance radar network in Europe, providing innovative quantitative precipitation estimation (QPE) based on polarimetric variables. This study circumvents Zh shortcomings by using specific differential-phase (Kdp) data from operational dual-polarization C-band weather radars. The rain intensity estimates based on a specific differential-phase data are immune to attenuation and less affected by hail contamination. In this study, for the first time, QPEs based on polarimetric observations by operational C-band weather radars and without any rain gauge adjustments are analyzed. The purpose is to estimate return periods for 1 h rainfall total computed from polarimetric weather radar data using non-adjusted QPEs based on R(Zh,Kdp) data and to compare the results with those derived using R(Zh) and rain gauge data. Only the warm period during the year is considered here, as most of the extreme precipitation events for such a duration occur for both places studied (Italy and Estonia) at this time. Limiting the dataset to warm periods also allows us to use the radar-based rainfall quantitative precipitation estimations, which are more reliable than the snowfall ones. Data from operational dual polarimetric C-band weather radar sites are used from both Italy and Estonia. Given climatologically homogeneous regions, this study demonstrates that polarimetric weather radar observations can provide reliable QPEs compared to single-polarization estimates with respect to rain gauges and that they can provide a reliable estimation of return periods of 1 h rainfall total, even for relatively short time series.</p

Introducing the Video In Situ Snowfall Sensor (VISSS)

The open-source Video In Situ Snowfall Sensor (VISSS) is introduced as a novel instrument for the characterization of particle shape and size in snowfall. The VISSS consists of two cameras with LED backlights and telecentric lenses that allow accurate sizing and combine a large observation volume with relatively high pixel resolution and a design that limits wind disturbance. VISSS data products include various particle properties such as maximum extent, cross-sectional area, perimeter, complexity, and sedimentation velocity. Initial analysis shows that the VISSS provides robust statistics based on up to 10 000 unique particle observations per minute. Comparison of the VISSS with the collocated PIP (Precipitation Imaging Package) and Parsivel instruments at Hyytiälä, Finland, shows excellent agreement with the Parsivel but reveals some differences for the PIP that are likely related to PIP data processing and limitations of the PIP with respect to observing smaller particles. The open-source nature of the VISSS hardware plans, data acquisition software, and data processing libraries invites the community to contribute to the development of the instrument, which has many potential applications in atmospheric science and beyond.</p

Dual-wavelength radar technique development for snow rate estimation: a case study from GCPEx

quantitative precipitation estimation (QPE) of snowfall has generally been expressed in power-law form between equivalent radar reflectivity factor (Ze) and liquid equivalent snow rate (SR). It is known that there is large variability in the prefactor of the power law due to changes in particle size distribution (PSD), density, and fall velocity, whereas the variability of the exponent is considerably smaller. The dual-wavelength radar reflectivity ratio (DWR) technique can improve SR accuracy by estimating one of the PSD parameters (characteristic diameter), thus reducing the variability due to the prefactor. The two frequencies commonly used in dual-wavelength techniques are Ku- and Ka-bands. The basic idea of DWR is that the snow particle size-to-wavelength ratio is falls in the Rayleigh region at Ku-band but in the Mie region at Ka-band. We propose a method for snow rate estimation by using NASA D3R radar DWR and Ka-band reflectivity observations collected during a long-duration synoptic snow event on 30–31 January 2012 during the GCPEx (GPM Cold-season Precipitation Experiment). Since the particle mass can be estimated using 2-D video disdrometer (2DVD) fall speed data and hydrodynamic theory, we simulate the DWR and compare it directly with D3R radar measurements. We also use the 2DVD-based mass to compute the 2DVD-based SR. Using three different mass estimation methods, we arrive at three respective sets of Z–SR and SR(Zh, DWR) relationships. We then use these relationships with D3R measurements to compute radar-based SR. Finally, we validate our method by comparing the D3R radar-retrieved SR with accumulated SR directly measured by a well-shielded Pluvio gauge for the entire synoptic event.</p

Translating the Physics of Snowfall to Radar-Based Validation of GPM

No abstract availabl

The importance of particle size distribution and internal structure for triple-frequency radar retrievals of the morphology of snow

The accurate representation of ice particles is essential for both remotely sensed estimates of clouds and precipitation and numerical models of the atmosphere. As it is typical in radar retrievals to assume that all snow is composed of aggregate snowflakes, both denser rimed snow and the mixed-phase cloud in which riming occurs may be under-diagnosed in retrievals and therefore difficult to evaluate in weather and climate models. Recent experimental and numerical studies have yielded methods for using triple-frequency radar measurements to interrogate the internal structure of aggregate snowflakes and to distinguish more dense and homogeneous rimed particles from aggregates. In this study we investigate which parameters of the morphology and size distribution of ice particles most affect the triple-frequency radar signature and must therefore be accounted for in order to carry out triple-frequency radar retrievals of snow. A range of ice particle morphologies are represented, using a fractal representation for the internal structure of aggregate snowflakes and homogeneous spheroids to represent graupel-like particles; the mass-size and area-size relations are modulated by a density factor. We find that the particle size distribution (PSD) shape parameter and the parameters controlling the internal structure of aggregate snowflakes both have significant influences on triple-frequency radar signature and are at least as important as that of the density factor. We explore how these parameters may be allowed to vary in order to prevent triple-frequency radar retrievals of snow from being over-constrained, using two case studies from the Biogenic Aerosols - Effects of Clouds and Climate (BAECC) 2014 field campaign at Hyytiala, Finland. In a case including heavily rimed snow followed by large aggregate snowflakes, we show that triple-frequency radar measurements provide a strong constraint on the PSD shape parameter, which can be estimated from an ensemble of retrievals; however, resolving variations in the PSD shape parameter has a limited impact on estimates of snowfall rate from radar. Particle density is more effectively constrained by the Doppler velocity than triple-frequency radar measurements, due to the strong dependence of particle fall speed on density. Due to the characteristic signatures of aggregate snowflakes, a third radar frequency is essential for effectively constraining the size of large aggregates. In a case featuring rime splintering, differences in the internal structures of aggregate snowflakes are revealed in the triple-frequency radar measurements. We compare retrievals assuming different aggregate snowflake models against in situ measurements at the surface and show significant uncertainties in radar retrievals of snow rate due to changes in the internal structure of aggregates. The importance of the PSD shape parameter and snowflake internal structure to triple-frequency radar retrievals of snow highlights that the processes by which ice particles interact may need to be better understood and parameterized before triple-frequency radar measurements can be used to constrain retrievals of ice particle morphology.Peer reviewe

Observing wind, aerosol particles, clouds and precipitation: Finland's new ground-based remote-sensing network

The Finnish Meteorological Institute, in collaboration with the University of Helsinki, has established a new ground-based remote-sensing network in Finland. The network consists of five topographically, ecologically and climatically different sites distributed from southern to northern Finland. The main goal of the network is to monitor air pollution and boundary layer properties in near real time, with a Doppler lidar and ceilometer at each site. In addition to these operational tasks, two sites are members of the Aerosols, Clouds and Trace gases Research InfraStructure Network (ACTRIS); a Ka band cloud radar at Sodankylä will provide cloud retrievals within CloudNet, and a multi-wavelength Raman lidar, PollyXT (POrtabLe Lidar sYstem eXTended), in Kuopio provides optical and microphysical aerosol properties through EARLINET (the European Aerosol Research Lidar Network). Three C-band weather radars are located in the Helsinki metropolitan area and are deployed for operational and research applications. We performed two inter-comparison campaigns to investigate the Doppler lidar performance, compare the backscatter signal and wind profiles, and to optimize the lidar sensitivity through adjusting the telescope focus length and data-integration time to ensure sufficient signal-to-noise ratio (SNR) in low-aerosol-content environments. In terms of statistical characterization, the wind-profile comparison showed good agreement between different lidars. Initially, there was a discrepancy in the SNR and attenuated backscatter coefficient profiles which arose from an incorrectly reported telescope focus setting from one instrument, together with the need to calibrate. After diagnosing the true telescope focus length, calculating a new attenuated backscatter coefficient profile with the new telescope function and taking into account calibration, the resulting attenuated backscatter profiles all showed good agreement with each other. It was thought that harsh Finnish winters could pose problems, but, due to the built-in heating systems, low ambient temperatures had no, or only a minor, impact on the lidar operation – including scanning-head motion. However, accumulation of snow and ice on the lens has been observed, which can lead to the formation of a water/ice layer thus attenuating the signal inconsistently. Thus, care must be taken to ensure continuous snow removal

Single hadron response measurement and calorimeter jet energy scale uncertainty with the ATLAS detector at the LHC

The uncertainty on the calorimeter energy response to jets of particles is derived for the ATLAS experiment at the Large Hadron Collider (LHC). First, the calorimeter response to single isolated charged hadrons is measured and compared to the Monte Carlo simulation using proton-proton collisions at centre-of-mass energies of sqrt(s) = 900 GeV and 7 TeV collected during 2009 and 2010. Then, using the decay of K_s and Lambda particles, the calorimeter response to specific types of particles (positively and negatively charged pions, protons, and anti-protons) is measured and compared to the Monte Carlo predictions. Finally, the jet energy scale uncertainty is determined by propagating the response uncertainty for single charged and neutral particles to jets. The response uncertainty is 2-5% for central isolated hadrons and 1-3% for the final calorimeter jet energy scale.Comment: 24 pages plus author list (36 pages total), 23 figures, 1 table, submitted to European Physical Journal

Standalone vertex finding in the ATLAS muon spectrometer

A dedicated reconstruction algorithm to find decay vertices in the ATLAS muon spectrometer is presented. The algorithm searches the region just upstream of or inside the muon spectrometer volume for multi-particle vertices that originate from the decay of particles with long decay paths. The performance of the algorithm is evaluated using both a sample of simulated Higgs boson events, in which the Higgs boson decays to long-lived neutral particles that in turn decay to bbar b final states, and pp collision data at √s = 7 TeV collected with the ATLAS detector at the LHC during 2011