Search

In Australia, hailstorms present considerable public safety and economic risks, where they are considered the most damaging natural hazard in terms of annual insured losses. Despite these impacts, the current climatological distribution of hailfall across the continent is still comparatively poorly understood. This study aims to supplement previous national hail climatologies, such as those based on environmental proxies or satellite radiometer data, with more direct radar-based hail observations. The heterogeneous and incomplete nature of the Australian radar network complicates this task and prompts the introduction of some novel methodological elements. We introduce an empirical correction technique to account for hail reflectivity biases at C band, derived by comparing overlapping C- and S-band observations. Furthermore, we demonstrate how object-based hail swath analysis may be used to produce resolution-invariant hail frequencies, and describe an interpolation method used to create a spatially continuous hail climatology. The maximum estimated size of hail (MESH) parameter is then applied to a mixture of over 50 operational radars in the Australian radar archive, resulting in the first nationwide, radar-based hail climatology. The spatiotemporal distribution of hailstorms is examined, including their physical characteristics, seasonal and diurnal frequency, and regional variations of such properties across the continent. [ABSTRACT FROM AUTHOR]

Doppler radars measure Doppler velocity within the [−VN, VN] range, where VN is the Nyquist velocity. Doppler velocities outside this range are "folded" within this interval. All Doppler "unfolding" techniques use the folded velocities themselves. In this work, we investigate the potential of using velocities derived from optical flow techniques applied to the radar reflectivity field for that purpose. The analysis of wind speed errors using six months of multi-Doppler wind retrievals showed that 99.9% of all points are characterized by errors smaller than 26 m s−1 below 5-km height, corresponding to a failure rate of less than 0.1% if optical flow winds were used to unfold Doppler velocities for VN = 26 m s−1. These errors largely increase above 5-km height, indicating that vertical continuity tests should be included to reduce failure rates at higher elevations. Following these results, we have developed the Two-step Optical Flow Unfolding (TOFU) technique, with the specific objective to accurately unfold Doppler velocities with VN = 26 m s−1. The TOFU performance was assessed using challenging case studies, comparisons with an advanced Doppler unfolding technique using higher Nyquist velocities, and 6 months of high VN (47.2 m s−1) data artificially folded to 26 m s−1. TOFU failure rates were found to be very low. Three main situations contributed to these errors: high low-level wind shear, elevated cloud layers associated with high winds, and radar data artifacts. Our recommendation is to use these unfolded winds as the first step of advanced Doppler unfolding techniques. Significance Statement: The potential of using optical flow winds operationally to accurately unfold Doppler velocities is demonstrated in this work. The operational significance is that the Nyquist velocity can confidently be reduced to 26 m s−1, allowing for extended first trip radar maximum range and reduced contamination from dual pulse repetition frequency artifacts. [ABSTRACT FROM AUTHOR]

Three-dimensional wind retrievals from ground-based Doppler radars have played an important role in meteorological research and nowcasting over the past four decades. However, in recent years, the proliferation of open-source software and increased demands from applications such as convective parameterizations in numerical weather prediction models has led to a renewed interest in these analyses. In this study, we analyze how a major, yet often-overlooked, error source effects the quality of retrieved 3D wind fields. Namely, we investigate the effects of spatial interpolation, and show how the common practice of pregridding radial velocity data can degrade the accuracy of the results. Alternatively, we show that assimilating radar data directly at their observation locations improves the retrieval of important dynamic features such as the rear flank downdraft and mesocyclone within supercells, while also reducing errors in vertical vorticity, horizontal divergence, and all three velocity components. Significance Statement: We can attempt to estimate the wind speed and direction within a weather system when two weather radars measure it simultaneously. However, radars do not scan the whole atmosphere at once—instead, they measure along many cross sections, each at different heights. We show that a method commonly used to stitch the observations together degrades the accuracy of the winds. Additionally, we describe a way to feed the data directly into the analysis without stitching it together first, and show that this improves the wind retrievals considerably. We hope these improvements will help researchers better understand how various weather systems work, and help forecasters warn for dangerous weather such as tornadoes. [ABSTRACT FROM AUTHOR]

We present an integrated framework that leverages multiple weather radar calibration and monitoring techniques to provide real-time diagnostics on reflectivity calibration, antenna pointing, and dual-polarization moments. This framework uses a volume-matching technique to track the absolute calibration of radar reflectivity with respect to the Global Precipitation Measurement (GPM) spaceborne radar, the relative calibration adjustment (RCA) technique to track relative changes in the radar calibration constant, the solar calibration technique to track daily change in solar power and antenna pointing error, and techniques that track properties of light-rain medium to monitor the differential reflectivity and dual-polarization moments. This framework allows for an evaluation of various calibration and monitoring techniques. For example, we found that a change in the RCA is highly correlated to a change in absolute calibration, with respect to GPM, if a change in antenna pointing can first be ruled out. It is currently monitoring 67+ radars from the Australian radar network. Because of the diverse and evolving nature of the Australian radar network, flexibility and modularity are at the core of the calibration framework. The framework can tailor its diagnostics to the specific characteristics of a radar (band, beamwidth, etc.). Because of its modularity, it can be expanded with new techniques to provide additional diagnostics (e.g., monitoring of radar sensitivity). The results are presented in an interactive dashboard at different level of details for a wide and diverse audience (radar engineers, researchers, forecasters, and management), and it is operational at the Australian Bureau of Meteorology. Significance Statement: Weather radars, like all instruments, require maintenance and upgrades. Rainfall measurements are highly variable and sensitive to change, and this can lead to inconsistencies within a radar network. Calibration is the process to counteract those inconsistencies. Any calibration requires a fixed standard to which the changed/upgraded radar can be compared. The SCAR calibration framework presented herein makes use of several standards to retrieve a full set of diagnostics about the radar data. We apply these techniques over the entire Australian weather radar network and demonstrate that, by using this integrated approach, absolute calibration can be achieved to within 1 dBZ of reflectivity, antenna pointing can be monitored within 0.1°, and the various measurements of the radars can be quality controlled. [ABSTRACT FROM AUTHOR]

We use dual-polarization C-band data collected in the Southern Ocean to examine the properties of snow observed during a voyage in the austral summer of 2018. Using existing forward modeling formalisms based on an assumption of Rayleigh scattering by soft spheroids, an optimal estimation algorithm is implemented to infer snow properties from horizontally polarized radar reflectivity, the differential radar reflectivity, and the specific differential phase. From the dual-polarization observables, we estimate ice water content qi, the mass-mean particle size Dm, and the exponent of the mass–dimensional relationship bm that, with several assumptions, allow for evaluation of snow bulk density, and snow number concentration. Upon evaluating the uncertainties associated with measurement and forward model errors, we determine that the algorithm can retrieve qi, Dm, and bm within single-pixel uncertainties conservatively estimated in the range 120%, 60%, and 40%, respectively. Applying the algorithm to open-cellular convection in the Southern Ocean, we find evidence for secondary ice formation processes within multicellular complexes. In stratiform precipitation systems we find snow properties and infer processes that are distinctly different from the shallow convective systems with evidence for riming and aggregation being common. We also find that embedded convection within the frontal system produces precipitation properties consistent with graupel. Examining 5 weeks of data, we show that snow in open-cellular cumulus has higher overall bulk density than snow in stratiform precipitation systems with implications for interpreting measurements from space-based active remote sensors. [ABSTRACT FROM AUTHOR]

Observations made by weather radars play a central role in many aspects of meteorological research and forecasting. These applications commonly require that radar data be supplied on a Cartesian grid, necessitating a coordinate transformation and interpolation from the radar's native spherical geometry using a process known as gridding. In this study, we introduce a variational gridding method and, through a series of theoretical and real data experiments, show that it outperforms existing methods in terms of data resolution, noise filtering, spatial continuity, and more. Known problems with existing gridding methods (Cressman weighted average and nearest neighbor/linear interpolation) are also underscored, suggesting the potential for substantial improvements in many applications involving gridded radar data, including operational forecasting, hydrological retrievals, and three-dimensional wind retrievals. [ABSTRACT FROM AUTHOR]

Jakarta, a megacity in Indonesia, experiences recurrent floods associated with heavy rainfall. Characteristics of subdaily rainfall and the local factors influencing rainfall around Jakarta have not been thoroughly investigated, primarily because of data limitations. In this study, we examine the frequency and intensity of hourly and daily rain rate, including spatial characteristics and variations across time scales. We use 6-min C-band Doppler radar and 1-min in situ data during 2009–12 to resolve spatial rain-rate characteristics at higher resolution than previous studies. A reflectivity–rain rate (Z–R) relationship is derived (Z = 102.7R1.75) and applied to estimate hourly rain rate. Our results show that rain rate around Jakarta is spatially inhomogeneous. In the rainy season [December–February (DJF)], rain rate exhibits statistical properties markedly different from other seasons, with much higher frequency of rain, but, on average, less intense rain rate. In all seasons, there is a persistent higher hourly and daily mean rain rate found over mountainous areas, indicating the importance of local orographic effects. In contrast, for hourly rain-rate extremes, peaks are observed mostly over the coastal land and lowland areas. For the diurnal cycle of mean rain rate, a distinct afternoon peak is found developing earlier in DJF and later in the dry season. This study has implications for other analyses of mesoscale rain-rate extremes in areas of complex topography and suggests that coarse-grain products may miss major features of the rain-rate variability identified in our study. Significance Statement: For many years, Jakarta and its surrounding regions have been repeatedly inundated by flooding triggered by short-duration heavy rainfall or rainfall accumulated over multiple days. Little is known about the distribution of local rainfall and how it differs between seasons. In this study, we used high-resolution C-band Doppler radar during 2009–12 to understand the characteristics of rainfall over this complex topography. The results demonstrate that the rainfall features vary spatially and seasonally. In the wet season, rainfall is more frequent but, on average, lighter relative to other seasons. In all seasons, the highest hourly and daily mean rain rate persistently occurs over the mountains, indicating the vital role of topography in generating rainfall in the region. [ABSTRACT FROM AUTHOR]

For mechanically scanning weather radars, precise pointing of the antenna is a key factor in ensuring accurate observation of the atmosphere at far range. Since operational radars typically scan the atmosphere using a series of 360° sweeps at fixed elevation angles, level scanning during azimuthal rotation is also important but often is not actively monitored after installation. One method of gauging pointing accuracy of a radar is to use solar interference that occurs as the antenna sweeps past the sun. By comparing the observed position of the sun with its known position, an estimate of pointing error in both elevation and azimuth can be obtained. A basic model for this error assumes that the radar sweep is perfectly level and that biases in elevation are therefore independent of azimuth. We extend this model to allow for the possibility that the plane of rotation may not be level. Consequently, the direction and severity of tilt may be diagnosed in addition to any constant error in elevation and azimuth pointing. The extended model was applied to a subset of radars from the Australian weather radar network, resulting in the discovery of several out-of-level radars. One radar, Captains Flat near Canberra, showed a severe tilt of 0.81° prompting inspection by a technician. This revealed that mounting studs on the pedestal of the radar tower were badly worn and loose. Correction of this issue resolved the tilt component of the diagnosed elevation error and prevented further mechanical damage to the instrument. [ABSTRACT FROM AUTHOR]

Understanding the key dynamical and microphysical mechanisms driving precipitation in the Snowy Mountains region of southeast Australia, including the role of orography, can help improve precipitation forecasts, which is of great value for efficient water management. An intensive observation campaign was carried out during the 2018 austral winter, providing a comprehensive range of ground-based observations across the Snowy Mountains. We used data from three vertically pointing rain radars, cloud radar, a PARSIVEL disdrometer, and a network of 76 pluviometers. The observations reveal that all of the precipitation events were associated with cold front passages. About half accumulated during the frontal passage associated with deep, fully glaciated cloud tops, while the rest occurred in the postfrontal environment and were associated with clouds with supercooled liquid water (SLW) tops. About three-quarters of the accumulated precipitation was observed under blocked conditions, likely associated with blocked stratiform orographic enhancement. Specifically, more than a third of the precipitation resulted from moist cloudless air being lifted over stagnant air, upwind from the barrier, creating SLW-top clouds. These SLW clouds then produced stratiform precipitation mostly over the upwind slopes and mountain tops, with hydrometeors reaching the mountain tops mostly as rimed snow. Two precipitation events were studied in detail, which showed that during unblocked conditions, orographic convection invigoration and unblocked stratiform enhancement were the two main mechanisms driving the precipitation, with the latter being more prevalent after the frontal passage. During these events, ice particle growth was likely dominated by vapor deposition and aggregation during the frontal periods, while riming dominated during the postfrontal periods. [ABSTRACT FROM AUTHOR]

A spatial mismatch between radar-based hail swaths and surface hail reports is commonly noted in meteorological literature. The discrepancy is partly due to hailstone advection and melting between detection aloft and observation at the ground. This study aims to mitigate this problem by introducing a model named HailTrack, which estimates hailfall at the surface using radar observations. The model operates by detecting, tracking, and collating hailstone trajectories using dual-polarized, dual-Doppler radar retrievals. Notable improvements in hailfall forecasts were observed through the use of HailTrack, and initializing the model with radar retrievals of hail differential reflectivity HDR was found to produce the most accurate hailfall estimates. The analysis of a case study in Brisbane, Australia, demonstrated that trajectory modeling significantly improved the correlation between hail swaths and hail-related insurance losses, increasing Heidke skill scores from 0.48 to 0.58. The accumulated kinetic energy of hailstone impacts also showed some skill in identifying areas that were exposed to particularly severe hailfall. Other unique impact estimates are presented, such as hailstone advection information and hailstone impact angle statistics. The potential to run the model in real time and produce short-term (10–15 min) nowcasts is also introduced. Model applications include improving radar-based hail climatologies, validating hail detection techniques and insurance claims data, and providing real-time hail impact maps to improve public awareness of hail risk. [ABSTRACT FROM AUTHOR]

Conditions of high ice water content (HIWC; defined herein as at least 1.0 g m-3) are often found in the anvils of convective systems and can cause engine damage and/or failure in aircraft. We use ice water content (IWC) retrievals from satellite-borne radar and lidar (CloudSat and CALIOP) to provide the first analysis of global HIWC frequency using 11 years of data (2007-17). Results show that HIWC is generally present in 1%-2% of CloudSat and CALIOP IWC retrievals between flight level 270 (FL270; 27000 ft or 8.230 km) and FL420 (42000 ft or 12.801 km) in areas with frequent convection. Similar rates of HIWC are found over midlatitude oceans at relatively low altitudes (below FL270). Possible nonconvective mechanisms for the formation of this low-level HIWC are discussed, as are the uncertainties suggesting that the results at these low altitudes are an overestimation of the true threat of HIWC to aircraft engines. The satellite IWC retrievals are also used to validate an HIWC diagnostic tool that provides storm-scale statistics on HIWC over the contiguous United States (CONUS) during the summer convective season (May-August from 2012 to 2019). Results over the CONUS suggest that HIWC over the Great Plains is highest in June, when a point in the region is under HIWC conditions for approximately 25 h of 30 days on average. The mean area-equivalent diameters of HIWC conditions in some areas of the Great Plains exceed 350 km, and the conditions can persist for 4-5 h. [ABSTRACT FROM AUTHOR]

Unfold Radar Velocity (UNRAVEL) is an open-source modular Doppler velocity dealiasing algorithm for weather radars. UNRAVEL is an algorithm that does not need external reference velocity data, making it easily applicable. The proposed algorithm includes 11 core modules and 2 dealiasing strategies. UNRAVEL is an iterative algorithm. The goal is to build the dealiasing results starting with the strictest possible continuity tests in azimuth and range and, after each step, relaxing the parameters to include more results from a progressively growing number of reference points. UNRAVEL also has modules that perform 3D continuity checks. Thanks to this modular design, the number of dealiasing strategies can be expanded in order to optimize the dealiasing results. While the first driver dealiases Doppler velocity from each tilt independently from one another, the second driver also performs a three-dimensional continuity check of the velocity using successive elevations. The proposed dealiasing algorithm is tested using severe weather data from an S-band Doppler radar that have been aliased to mimic aliased radial velocity patterns that would be observed by a C-band Doppler radar. Artificially aliasing S-band data permits creation of a reference to which the performance of various dealiasing techniques can be compared. Comparisons show that UNRAVEL consistently outperforms other established dealiasing algorithms for the test period selected in this work. [ABSTRACT FROM AUTHOR]

Attenuation of the W-band (95 GHz) radar signal by atmospheric ice particles has long been neglected in cloud microphysics studies. In this work, 95-GHz airborne multibeam cloud radar observations in tropical stratiform ice anvils are used to estimate vertical profiles of 95-GHz attenuation. Two techniques are developed and compared, using very different assumptions. The first technique examines statistical reflectivity differences between repeated aircraft passes through the same cloud mass at different altitudes. The second technique exploits reflectivity differences between two different pathlengths through the same cloud, using the multibeam capabilities of the cloud radar. Using the first technique, the two-way attenuation coefficient produced by stratiform ice particles ranges between 1 and 1.6 dB km−1 for reflectivities between 13 and 18 dBZ, with an expected increase of attenuation with reflectivity. Using the second technique, the multibeam results confirm these high attenuation coefficient values and expand the reflectivity range, with typical attenuation coefficient values of up to 3–4 dB km−1 for reflectivities of 20 dBZ. The potential impact of attenuation on precipitating-ice-cloud microphysics retrievals is quantified using vertical profiles of the mean and the 99th percentile of ice water content derived from noncorrected and attenuation-corrected reflectivities. A large impact is found on the 99th percentile of ice water content, which increases by 0.3–0.4 g m−3 up to 11-km height. Finally, T-matrix calculations of attenuation constrained by measured particle size distributions, ice crystal mass–size, and projected area–size relationships are found to largely underestimate cloud radar attenuation estimates. [ABSTRACT FROM AUTHOR]

Interpolation of ground-based radar measurements is required when mapping data from their native spherical coordinates to a Cartesian grid. For reflectivity the question arises as to whether this processing should be performed in units of Z (mm6 m−3) or dBZ. This study addresses this question using one year of data from three radars, operating in diverse climates across Australia. For each radar, a subset of 800 volume scans is processed to identify "triads"—groups of three consecutive gates with valid data—in each of the three coordinate directions: range, azimuth, and elevation. For every triad, the reflectivity at the central gate is estimated by linearly interpolating between the outer two gates in both Z and dBZ. The resulting values are then compared with the true reflectivity at the central gate to quantify the interpolation errors. For all three sites and in all three coordinate directions, we find that interpolation in Z is more accurate on average, especially in regions of high reflectivity and strong reflectivity gradient (i.e., convective cores). However, interpolation in dBZ is better in regions of low and monotonically increasing/decreasing reflectivity. It is therefore recommended that reflectivities be converted from dBZ to Z prior to interpolation except when identifying echo-top height or other low-reflectivity boundaries. [ABSTRACT FROM AUTHOR]

The stability and accuracy of weather radar reflectivity calibration are imperative for quantitative applications, such as rainfall estimation, severe weather monitoring and nowcasting, and assimilation in numerical weather prediction models. Various radar calibration and monitoring techniques have been developed, but only recently have integrated approaches been proposed, that is, using different calibration techniques in combination. In this paper the following three techniques are used: 1) ground clutter monitoring, 2) comparisons with spaceborne radars, and 3) the self-consistency of polarimetric variables. These techniques are applied to a C-band polarimetric radar (CPOL) located in the Australian tropics since 1998. The ground clutter monitoring technique is applied to each radar volumetric scan and provides a means to reliably detect changes in calibration, relative to a baseline. It is remarkably stable to within a standard deviation of 0.1 dB. To obtain an absolute calibration value, CPOL observations are compared to spaceborne radars on board TRMM and GPM using a volume-matching technique. Using an iterative procedure and stable calibration periods identified by the ground echoes technique, we improve the accuracy of this technique to about 1 dB. Finally, we review the self-consistency technique and constrain its assumptions using results from the hybrid TRMM–GPM and ground echo technique. Small changes in the self-consistency parameterization can lead to 5 dB of variation in the reflectivity calibration. We find that the drop-shape model of Brandes et al. with a standard deviation of the canting angle of 12° best matches our dataset. [ABSTRACT FROM AUTHOR]

Current climate models cannot resolve individual convective clouds, and hence parameterizations are needed. The primary goal of convective parameterization is to represent the bulk impact of convection on the gridbox-scale variables. Spectral convective parameterizations also aim to represent the key features of the subgrid-scale convective cloud field such as cloud-top-height distribution and in-cloud vertical velocities in addition to precipitation rates. Ground-based radar retrievals of these quantities have been made available at Darwin, Australia, permitting direct comparisons of internal parameterization variables and providing new observational references for further model development. A spectral convective parameterization [the convective cloud field model (CCFM)] is discussed, and its internal equation of motion is improved. Results from the ECHAM–HAM model in single-column mode using the CCFM and the bulk mass flux Tiedtke–Nordeng scheme are compared with the radar retrievals at Darwin. The CCFM is found to outperform the Tiedtke–Nordeng scheme for cloud-top-height and precipitation-rate distributions. Radar observations are further used to propose a modified CCFM configuration with an aerodynamic drag and reduced entrainment parameter, further improving both the convective cloud-top-height distribution (important for large-scale impact of convection) and the in-cloud vertical velocities (important for aerosol activation). This study provides a new development in the CCFM, improving the representation of convective cloud spectrum characteristics observed in Darwin. This is a step toward an improved representation of convection and ultimately of aerosol effects on convection. It also shows how long-term radar observations of convective cloud properties can help constrain parameters of convective parameterization schemes. [ABSTRACT FROM AUTHOR]

Calibration error represents a significant source of uncertainty in quantitative applications of ground-based radar (GR) reflectivity data. Correcting it requires knowledge of the true reflectivity at well-defined locations and times during a volume scan. Previous work has demonstrated that observations from certain spaceborne radar (SR) platforms may be suitable for this purpose. Specifically, the Ku-band precipitation radars on board the Tropical Rainfall Measuring Mission (TRMM) satellite and its successor, the Global Precipitation Measurement (GPM) mission Core Observatory satellite together provide nearly two decades of wellcalibrated reflectivity measurements over low-latitude regions (±35°). However, when comparing SR andGR reflectivities, great care must be taken to account for differences in instrument sensitivity and frequency, and to ensure that the observations are spatially and temporally coincident. Here, a volume-matching method, developed as part of the ground validation network for GPM, is adapted and used to quantify historical calibration errors for three S-band radars in the vicinity of Sydney, Australia. Volume-matched GR-SR sample pairs are identified over a 7-yr period and carefully filtered to isolate reflectivity differences associated with GR calibration error. These are then used in combination with radar engineering work records to derive a piecewise-constant time series of calibration error for each site. The efficacy of this approach is verified through comparisons between GR reflectivities in regions of overlapping coverage, with improved agreement when the estimated errors are removed. [ABSTRACT FROM AUTHOR]

Understanding phase transitions in mixed-phase clouds is of great importance because the hydrometeor phase controls the lifetime and radiative effects of clouds. In high latitudes, these cloud radiative effects have a crucial impact on the surface energy budget and thus on the evolution of the ice cover. For a springtime low-level mixed-phase stratiform cloud case from Barrow, Alaska, a unique combination of instruments and retrieval methods is combined with multiple modeling perspectives to determine key processes that control cloud phase partitioning. The interplay of local cloud-scale versus large-scale processes is considered. Rapid changes in phase partitioning were found to be caused by several main factors. Major influences were the large-scale advection of different air masses with different aerosol concentrations and humidity content, cloud-scale processes such as a change in the thermodynamical coupling state, and local-scale dynamics influencing the residence time of ice particles. Other factors such as radiative shielding by a cirrus and the influence of the solar cycle were found to only play a minor role for the specific case study (11-12 March 2013). For an even better understanding of cloud phase transitions, observations of key aerosol parameters such as profiles of cloud condensation nucleus and ice nucleus concentration are desirable. [ABSTRACT FROM AUTHOR]

Doppler cloud radars are amazing tools to characterize cloud and fog properties and to improve their representation in models. However, commercially available cloud radars (35 and 95 GHz) are still very expensive, which hinders their widespread deployment. This study presents the development of a lower-cost semioperational 95-GHz Doppler cloud radar called the Bistatic Radar System for Atmospheric Studies (BASTA). To drastically reduce the cost of the instrument, a different approach is used compared to traditional pulsed radars: instead of transmitting a large amount of energy for a very short time period (as a pulse), a lower amount of energy is transmitted continuously. By using a specific signal processing technique, the radar can challenge expensive radars and provide high-quality measurements of cloud and fog. The latest version of the instrument has a sensitivity of about −50 dB Z at 1 km for 3-s integration and a vertical resolution of 25 m. The BASTA radar currently uses four successive modes for specific applications: the 12.5-m vertical resolution mode is dedicated to fog and low clouds, the 25-m mode is for liquid and ice midtropospheric clouds, and the 100- and 200-m modes are ideal for optically thin high-level ice clouds. The advantages of such a radar for calibration procedures and field operations are also highlighted. The radar comes with a set of products dedicated to cloud and fog studies. For instance, cloud mask, corrected Doppler velocity, and multimode products combining the high-sensitivity mode and high-resolution modes are provided. [ABSTRACT FROM AUTHOR]

Cumulus parameterizations in general circulation models (GCMs) frequently apply mass-flux schemes in their description of tropical convection. Mass flux constitutes the product of the fractional area covered by cumulus clouds in a model grid box and the vertical velocity within the cumulus clouds. The cumulus area fraction profiles can be derived from precipitating radar reflectivity volumes. However, the vertical velocities are difficult to observe, making the evaluation of mass-flux schemes difficult. In this paper, the authors develop and evaluate a parameterization of vertical velocity in convective (cumulus) clouds using only radar reflectivities collected by a C-band polarimetric research radar (CPOL), operating at Darwin, Australia. The parameterization is trained using vertical velocity retrievals from a dual-frequency wind profiler pair located within the field of view of CPOL. The parametric model uses two inputs derived from CPOL reflectivities: the 0-dB Z echo-top height (0-dB Z ETH) and a height-weighted column reflectivity index ( ZHWT). The 0-dB Z ETH determines the shape of the vertical velocity profile, while ZHWT determines its strength. The evaluation of these parameterized vertical velocities using (i) the training dataset, (ii) an independent wind-profiler-based dataset, and (iii) 1 month of dual-Doppler vertical velocity retrievals indicates that the statistical representation of vertical velocity is reasonably accurate up to the 75th percentile. However, the parametric model underestimates the extreme velocities. The method allows for the derivation of cumulus mass flux and its variability on current GCM scales based only on reflectivities from precipitating radar, which could be valuable to modelers. [ABSTRACT FROM AUTHOR]

Two new algorithms for hydrometeor classification using polarimetric radar observations are developed based on prototypes derived by applying clustering techniques (Part I of this two-part paper). Each prototype is defined as a probability distribution of the polarimetric variables and ambient temperature corresponding to a hydrometeor type. The first algorithm is a maximum prototype likelihood classifier that uses all prototypes attributed to the different hydrometeor types in Part I. The hydrometeor type is assigned as the prototype with the highest likelihood when comparing the polarimetric variables and temperature with each prototype. The second algorithm is a Bayesian classifier that uses the probability density functions (PDFs) as derived from the prototype set associated with the identical hydrometeor type. The posteriori probability in the Bayesian method is calculated from a combination of the PDFs and the prior probability, the maximum of which corresponds to the most likely hydrometeor type. The respective merits of the two techniques are discussed. The two classifiers are applied to CP-2 S-band radar observations of two hailstorms that occurred between 16 and 20 November 2008, including the so-called Gap storm, which produced a devastating microburst and large hail at the ground. Results from the classifiers are compared with those derived using the well-established National Center for Atmospheric Research fuzzy logic classifier. In general, good agreement is found, yielding overall confidence in the robustness of the new classifiers. However, large differences are found for the melting ice and ice crystal categories, which will need to be studied further. [ABSTRACT FROM AUTHOR]

Hydrometeor classification methods using polarimetric radar variables rely on probability density functions (PDFs) or membership functions derived empirically or by using electromagnetic scattering calculations. This paper describes an objective approach based on cluster analysis to deriving the PDFs. An iterative procedure with K-means clustering and expectation-maximization clustering based on Gaussian mixture models is developed to generate a series of prototypes for each hydrometeor type from several radar scans. The prototypes are then grouped together to produce a PDF for each hydrometeor type, which is modeled as a Gaussian mixture. The cluster-based method is applied to polarimetric radar data collected with the CP-2 S-band radar near Brisbane, Queensland, Australia. The results are illustrated and compared with theoretical classification boundaries in the literature. Some notable differences are found. Automated hydrometeor classification algorithms can be built using the PDFs of polarimetric variables associated with each hydrometeor type presented in this paper. [ABSTRACT FROM AUTHOR]

In this study, methods of convective/stratiform precipitation classification and surface rain-rate estimation based on the Atmospheric Radiation Measurement Program (ARM) cloud radar measurements were developed and evaluated. Simultaneous and collocated observations of the Ka-band ARM zenith radar (KAZR), two scanning precipitation radars [NCAR S-band/Ka-band Dual Polarization, Dual Wavelength Doppler Radar (S-PolKa) and Texas A&M University Shared Mobile Atmospheric Research and Teaching Radar (SMART-R)], and surface precipitation during the Dynamics of the Madden-Julian Oscillation/ARM MJO Investigation Experiment (DYNAMO/AMIE) field campaign were used. The motivation of this study is to apply the unique long-term ARM cloud radar observations without accompanying precipitation radars to the study of cloud life cycle and precipitation features under different weather and climate regimes. The resulting convective/stratiform classification from KAZR was evaluated against precipitation radars. Precipitation occurrence and classified convective/stratiform rain fractions from KAZR compared favorably to the collocated SMART-R and S-PolKa observations. Both KAZR and S-PolKa radars observed about 5% precipitation occurrence. The convective (stratiform) precipitation fraction is about 18% (82%). Collocated disdrometer observations of two days showed an increased number concentration of small and large raindrops in convective rain relative to dominant small raindrops in stratiform rain. The composite distributions of KAZR reflectivity and Doppler velocity also showed distinct structures for convective and stratiform rain. These evidences indicate that the method produces physically consistent results for the two types of rain. A new KAZR-based, two-parameter [the gradient of accumulative radar reflectivity Z e (GAZ) below 1 km and near-surface Z e] rain-rate estimation procedure was developed for both convective and stratiform rain. This estimate was compared with the exponential Z-R (reflectivity-rain rate) relation. The relative difference between the estimated and surface-measured rainfall rates showed that the two-parameter relation can improve rainfall estimation relative to the Z-R relation. [ABSTRACT FROM AUTHOR]

Clouds strongly affect the absorption and reflection of shortwave and longwave radiation in the atmosphere. A key bias in climate models is related to excess absorbed shortwave radiation in the high-latitude Southern Ocean. Model evaluation studies attribute these biases in part to midtopped clouds, and observations confirm significant midtopped clouds in the zone of interest. However, it is not yet clear what cloud properties can be attributed to the deficit in modeled clouds. Present approaches using observed cloud regimes do not sufficiently differentiate between potentially distinct types of midtopped clouds and their meteorological contexts. This study presents a refined set of midtopped cloud subregimes for the high-latitude Southern Ocean, which are distinct in their dynamical and thermodynamic background states. Active satellite observations from CloudSat and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations ( CALIPSO) are used to study the macrophysical structure and microphysical properties of the new cloud regimes. The subgrid-scale variability of cloud structure and microphysics is quantified within the cloud regimes by identifying representative physical cloud profiles at high resolution from the radar-lidar (DARDAR) cloud classification mask. The midtopped cloud subregimes distinguish between stratiform clouds under a high inversion and moderate subsidence; an optically thin cold-air advection cloud regime occurring under weak subsidence and including altostratus over low clouds; optically thick clouds with frequent deep structures under weak ascent and warm midlevel anomalies; and a midlevel convective cloud regime associated with strong ascent and warm advection. The new midtopped cloud regimes for the high-latitude Southern Ocean will provide a refined tool for model evaluation and the attribution of shortwave radiation biases to distinct cloud processes and properties. [ABSTRACT FROM AUTHOR]

C-band polarimetric radar measurements spanning two wet seasons are used to perform a critical evaluation of two algorithms for the classification of stratiform and convective precipitation. The first approach is based on the horizontal texture of the radar reflectivity field (two classes: stratiform, convective), and the second approach is based on the properties of the drop size distribution (DSD) parameters as derived from a set of polarimetric variables (three classes: stratiform, mixed, convective). To investigate how well those two methods compare quantitatively, probability density functions of reflectivity, rain rate, 5-dB Z echo top height, and DSD parameters (namely, the median volume diameter and the 'generalized' intercept parameter) are built. The study found that while the two methods agree well on the identification of stratiform precipitation, large differences are obtained for convective rainfall. The texture-based approach seems to classify too many points as being of convective nature compared to the DSD-based method. Among the points that are classified as convective by the texture-based approach, 25% correspond to low concentration of relatively small particles associated with rain rates below 10 mm h−1. This large proportion of unrealistically low convective rain rates is not produced by the DSD-based approach, which only classifies 4% of the convective points with rain rates below 10 mm h−1. These points were found to be mainly isolated points embedded within stratiform precipitation and associated with low cloud-top height, suggesting a misclassification of the texture-based approach. Thus, to improve the statistics of the convective class, three modified equations of the peakedness criterion used in the radar-based algorithm are proposed to decrease the number of misclassified points. [ABSTRACT FROM AUTHOR]

This study presents a summary of the properties of deep convective updraft and downdraft cores over the central plains of the United States, accomplished using a novel and now-standard Atmospheric Radiation Measurement Program (ARM) scanning mode for a commercial wind-profiler system. A unique profiler-based hydrometeor fall-speed correction method modeled for the convective environment was adopted. Accuracy of the velocity retrievals from this effort is expected to be within 2 m s−1, with minimal bias and base core resolution expected near 1 km. Updraft cores are found to behave with height in reasonable agreement with aircraft observations of previous continental convection, including those of the Thunderstorm Project. Intense updraft cores with magnitudes exceeding 15 m s−1 are routinely observed. Downdraft cores are less frequently observed, with weaker magnitudes than updrafts. Weak, positive correlations are found between updraft intensity (maximum) and updraft diameter length (coefficient r to 0.5 aloft). Negligible correlations are observed for downdraft core lengths and intensity. [ABSTRACT FROM AUTHOR]

C-band polarimetric radar measurements spanning two wet seasons are used to study the effects of the large-scale environment on the statistical properties of stratiform and convective rainfall around Darwin, Australia. The rainfall physical properties presented herein are the reflectivity fields, daily rainfall accumulations and raining area, rain rates, and drop size distribution (DSD) parameters (median volume diameter and 'normalized' intercept parameter). Each of these properties is then analyzed according to five different atmospheric regimes and further separated into stratiform or convective rain categories following a DSD-based approach. The regimes, objectively identified by radiosonde thermodynamic and wind measurements, represent typical wet-season atmospheric conditions: the active monsoon regime, the 'break' periods, the 'buildup' regime, the trade wind regime, and a mixture of inactive/break periods. The large-scale context is found to strongly modulate rainfall and cloud microphysical properties. For example, during the active monsoon regime, the daily rain accumulation is higher than in the other regimes, while this regime is associated with the lowest rain rates. Precipitation in this active monsoon regime is found to be widespread and mainly composed of small particles in high concentration compared to the other regimes. Vertical profiles of reflectivity and DSD parameters suggest that warm rain processes are dominant during this regime. In contrast, rainfall properties in the drier regimes (trade wind/buildup regimes) are mostly of continental origin, with rain rates higher than in the moister regimes. In these drier regimes, precipitation is mainly formed of large raindrops in relatively low concentration due to a larger contribution of the ice microphysical processes on the rainfall formation. [ABSTRACT FROM AUTHOR]

Comparisons between direct measurements and modeled values of vertical air motions in precipitating systems are complicated by differences in temporal and spatial scales. On one hand, vertically profiling radars more directly measure the vertical air motion but do not adequately capture full storm dynamics. On the other hand, vertical air motions retrieved from two or more scanning Doppler radars capture the full storm dynamics but require model constraints that may not capture all updraft features because of inadequate sampling, resolution, numerical constraints, and the fact that the storm is evolving as it is scanned by the radars. To investigate the veracity of radar-based retrievals, which can be used to verify numerically modeled vertical air motions, this article presents several case studies from storm events around Darwin, Northern Territory, Australia, in which measurements from a dual-frequency radar profiler system and volumetric radar-based wind retrievals are compared. While a direct comparison was not possible because of instrumentation location, an indirect comparison shows promising results, with volume retrievals comparing well to those obtained from the profiling system. This prompted a statistical analysis of an extended period of an active monsoon period during the Tropical Warm Pool International Cloud Experiment (TWP-ICE). Results show less vigorous deep convective cores with maximum updraft velocities occurring at lower heights than some cloud-resolving modeling studies suggest. [ABSTRACT FROM AUTHOR]

This study illustrates the high potential of RALI, the French airborne radar-lidar instrument, for studying cloud processes and evaluating satellite products when satellite overpasses are available. For an Arctic nimbostratus ice cloud collected on 1 April 2008 during the Polar Study using Aircraft, Remote Sensing, Surface Measurements and Models, of Climate, Chemistry, Aerosols, and Transport (POLARCAT) campaign, the capability of this synergistic instrument to retrieve cloud properties and to characterize the cloud phase at scales smaller than a kilometer, which is crucial for cloud process analysis, is demonstrated. A variational approach, which combines radar and lidar, is used to retrieve the ice-water content (IWC), extinction, and effective radius. The combination of radar and lidar is shown to provide better retrievals than do stand-alone methods and, in general, the radar overestimates and the lidar underestimates IWC. As the sampled ice cloud was simultaneously observed by CloudSat and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations ( CALIPSO) satellites, a new way to assess satellite cloud products by combining in situ and active remote sensing measurements is identified. It was then possible to compare RALI to three satellite ice cloud products: CloudSat, CALIPSO, and the Cloud-Aerosol-Water-Radiation Interactions (ICARE) center's radar-lidar project (DARDAR). [ABSTRACT FROM AUTHOR]

Doppler radar measurements at different frequencies (50 and 2835 MHz) are used to characterize the terminal fall speed of hydrometeors and the vertical air motion in tropical ice clouds and to evaluate statistical methods for retrieving these two parameters using a single vertically pointing cloud radar. For the observed vertical air motions, it is found that the mean vertical air velocity in ice clouds is small on average, as is assumed in terminal fall speed retrieval methods. The mean vertical air motions are slightly negative (downdraft) between the melting layer (5-km height) and 6.3-km height, and positive (updraft) above this altitude, with two peaks of 6 and 7 cm s−1 at 7.7- and 9.7-km height. For the retrieved hydrometeor terminal fall speeds, it is found that the variability of terminal fall speeds within narrow reflectivity ranges is typically within the acceptable uncertainties for using terminal fall speeds in ice cloud microphysical retrievals. This study also evaluates the performance of previously published statistical methods of separating terminal fall speed and vertical air velocity from vertically pointing Doppler radar measurements using the 50-/2835-MHz radar retrievals as a reference. It is found that the variability of the terminal fall speed-radar reflectivity relationship ( Vt- Ze) is large in ice clouds and cannot be parameterized accurately with a single relationship. A well-defined linear relationship is found between the two coefficients of a power-law Vt- Ze relationship, but a more accurate microphysical retrieval is obtained using Doppler velocity measurements to better constrain the Vt- Ze relationship for each cloud. When comparing the different statistical methods to the reference, the distribution of terminal fall speed residual is wide, with most residuals being in the ±30-40 cm s−1 range about the mean. The typical mean residual ranged from 15 to 20 cm s−1, with different methods having mean residuals of <10 cm s−1 at some heights, but not at the same heights for all methods. The so-called Vt- Ze technique was the most accurate above 9-km height, and the running-mean technique outperformed the other techniques below 9-km height. Sensitivity tests of the running-mean technique indicate that the 20-min average is the best trade-off for the type of ice clouds considered in this analysis. A new technique is proposed that incorporates simple averages of Doppler velocity for each ( Ze, H) couple in a given cloud. This technique, referred to as DOP- Ze- H, was found to outperform the three other methods at most heights, with a mean terminal fall residual of <10 cm s−1 at all heights. This error magnitude is compatible with the use of such retrieved terminal fall speeds for the retrieval of microphysical properties. [ABSTRACT FROM AUTHOR]

This article investigates the source and impact of artifacts produced by ordered linear interpolation techniques on variationally retrieved updraft intensities. Qualitative reasoning for the generation of periodic perturbations in gridded products is presented, and a simple analytical investigation into the impact of gridding artifacts on updraft retrieval is carried out. By projecting a nonconvergent flow typical of Darwin, Australia, onto the viewing geometry of a scanning radar, a numerical assessment of the impact of gridding artifacts is carried out. A simple enhancement to ordered linear interpolation, mixed-order linear interpolation, is proposed to reduce gridding artifacts. Radial velocity grids produced using both techniques are used to investigate the generation of spurious updrafts, with the simple ordered linear interpolation technique producing erroneous updrafts on the order of 2 m s−1. To investigate the impact on vertical velocities retrieved from a real weather event, radar-derived measurements taken during the active monsoon phase of Tropical Warm Pool International Cloud Experiment are gridded using both techniques, and vertical velocities are retrieved and contrasted. [ABSTRACT FROM AUTHOR]

A correction to the article on cloud thermodynamic phase for clustering algorithm is presented.