Search

The Royal Netherlands Meteorological Institute (KNMI) operates two dual-polarization C-band weather radars in simultaneous transmission and reception (STAR; i.e., horizontally and vertically polarized pulses are transmitted simultaneously) mode, providing 2D radar rainfall products. Despite the application of Doppler and speckle filtering, remaining nonmeteorological echoes (especially sea clutter) mainly due to anomalous propagation still pose a problem. This calls for additional filtering algorithms, which can be realized by means of polarimetry. Here we explore the effectiveness of the open-source wradlib fuzzy echo classification and clutter identification based on polarimetric moments. Based on our study, this has recently been extended with the depolarization ratio and clutter phase alignment as new decision variables. Optimal values for weights of the different membership functions and threshold are determined employing a 4-h calibration dataset from one radar. The method is applied to a full year of volumetric data from the two radars in the Dutch temperate climate. The verification focuses on the presence of remaining nonmeteorological echoes by mapping the number of exceedances of radar reflectivity factors for given thresholds. Moreover, accumulated rainfall maps are obtained to detect unrealistically large rainfall depths. The results are compared to those for which no further filtering has been applied. Verification against rain gauge data reveals that only a little precipitation is removed. Because the fuzzy logic algorithm removes many nonmeteorological echoes, the practice to composite data from both radars in logarithmic space to hide these echoes is abandoned and replaced by linearly averaging reflectivities. [ABSTRACT FROM AUTHOR]

Anovel way to process polarimetric radar data collected via plan position indicator (PPI) scans and display those data in a time-height format is introduced. The columnar vertical profile (CVP) methodology uses radar data collected via multiple elevation scans, limited to data within a set region in range and azimuth relative to the radar, to create vertical profiles of polarimetric radar data representative of that limited region in space. This technique is compared to others existing in the literature, and various applications are discussed. Polarimetric ice microphysical retrievals are performed on CVPs created within the stratiform rain region of two mesoscale convective systems sampled during two field campaigns, where CVPs follow the track of research aircraft. Aircraft in situ data are collocated to microphysical retrieval data, and the accuracy of these retrievals is tested against other retrieval techniques in the literature. [ABSTRACT FROM AUTHOR]

Sampling error stems from the use of ensemble-based data assimilation (EDA) with a limited ensemble size and can result in spurious background error covariances, leading to false analysis corrections. The WRF-LETKF radar assimilation system (WLRAS) is performed separately with 256 and 40 members to investigate the characteristics of convective-scale sampling errors in the EDA and its impact on precipitation prediction based on a heavy rainfall event on 16 June 2008. The results suggest that the sampling errors for this event are sensitive to the relationships between the simulated observations and model variables, the intensity of reflectivity, and how the prevailing wind projects to the radial wind in the areas that the radar cannot resolve U or V wind. The sampling errors lead to an underprediction of heavy rainfall when the horizontal localization radius is inadequately large, but this can be mitigated when a more accurate moisture condition is provided. In addition, being able to use a larger vertical localization plays an important role in providing necessary adjustments for representing the vertical thermodynamical structure of convection, which further improves precipitation prediction. A strategy mitigating the impact of sampling errors associated with the limitation of radial wind measurement by inflating the observation error over sensitive areas can bring benefits to precipitation prediction. [ABSTRACT FROM AUTHOR]

Weather radar provides real-time, spatially distributed precipitation estimates, whereas traditional gauge data are restricted in space. The use of radar quantitative precipitation estimates (QPEs) as an input of hydrological models for hydrometeorological applications has increased with the development of weather radar worldwide. New dual-polarization technology and algorithms are showing improvements to radar QPEs. This study evaluates radar QPEs from C-band radar at King City, Canada (WKR), and NEXRAD S-band radar at Buffalo, New York (KBUF), to verify the reliability and accuracy for operational use in the Humber River (semiurban) and Don River (urban) watersheds in the Greater Toronto Area (GTA), Canada. Twenty rainfall events that occurred from 2011 to 2017 were determined from hourly gauge measurements and compared with nine radar QPEs. Rain rates were estimated with different algorithms using three dualpolarized reflectivity values: horizontal reflectivity Z, differential reflectivity ZDR, and specific differential phase KDP. The correlation coefficient, bias, detection, and root-mean-square error were calculated and averaged over all events for each gauge station to show the spatial distribution and in a similar pattern to represent the variation by the event. The quality of the results in terms of accuracy and reliability indicates that the radar QPEs from KBUF S-band and WKR C-band multiparameter rain rate estimators can be effectively used as precipitation forcing of hydrological models for hydrometeorological applications. The high spatiotemporal resolution, long-term data archive, and good percent detection of radar QPEs can facilitate hydrometeorological applications by providing a continuous time series for hydrological models. [ABSTRACT FROM AUTHOR]

Thunderstorms that produce hail accumulations at the surface can impact residents by obstructing roadways, closing airports, and causing localized flooding from hail-clogged drainages. These storms have recently gained an increased interest within the scientific community. However, differences that are observable in real time between these storms and storms that produce nonimpactful hail accumulations have yet to be documented. Similarly, the characteristics within a single storm that are useful to quantify or predict hail accumulations are not fully understood. This study uses lightning and dual-polarization radar data to characterize hail accumulations from three storms that occurred on the same day along the Colorado–Wyoming Front Range. Each storm’s characteristics are verified against radar-derived hail accumulation maps and in situ observations. The storms differed in maximum accumulation, either producing 22 cm, 7 cm, or no accumulation. The magnitude of surface hail accumulations is found to be dependent on a combination of in-cloud hail production, storm translation speed, and hailstone melting. The optimal combination for substantial hail accumulations is enhanced in-cloud hail production, slow storm speed, and limited hailstone melting. However, during periods of similar in-cloud hail production, lesser accumulations are derived when storm speed and/or hailstone melting, identified by radar presentation, is sufficiently large. These results will aid forecasters in identifying when hail accumulations are occurring in real time. [ABSTRACT FROM AUTHOR]

Specification of suitable initial conditions to accurately forecast high-impact weather events associated with intense thunderstorms still poses a significant challenge for convective-scale forecasting. Radar data assimilation has been showing encouraging results to produce an accurate estimate of the state of the atmosphere at the mesoscale, as it combines high-spatiotemporal-resolution observations with convection-permitting numerical weather prediction models. However, many open questions remain regarding the configuration of state-of-the-art data assimilation systems at the mesoscale and their potential impact upon short-range weather forecasts. In this work, several observing system simulation experiments of a mesoscale convective system were performed to assess the sensitivity of the local ensemble transform Kalman filter to both relaxation-to-prior spread (RTPS) inflation and horizontal localization of the error covariance matrix. Realistic large-scale forcing and model errors have been taken into account in the simulation of reflectivity and Doppler velocity observations. Overall, the most accurate analyses in terms of RMSE were produced with a relatively small horizontal localization cutoff radius (;3.6–7.3 km) and large RTPS inflation parameter (;0.9–0.95). Additionally, the impact of horizontal localization on short-range ensemble forecast was larger compared to inflation, almost doubling the lead times up to which the effect of using a more accurate state to initialize the forecast persisted. [ABSTRACT FROM AUTHOR]

The impacts of the differential phase of incident radar waves (ψi) on measured differential reflectivity (ZDR), differential phase, and correlation coefficient from ice cloud particles are presented for radars employing simultaneous transmission and reception of orthogonally polarized waves (SHV radar design). The maximal values of ZDR and the differential phase upon scattering (δ) from ice particles are obtained as functions of ψi. It is shown that SHV δ from ice particles can exceed a dozen degrees whereas the intrinsic δ is of a few hundredths of a degree. In melting layers, the δ values from particles obeying the Rayleigh scattering law can be several degrees depending on ψi so that, to explain such δ values, an assumption of resonance scattering is not necessary. The phase δ affects the estimation of specific differential phase (KDP) in icy media and, therefore, the phase δ should be measured. The radar differential phase upon transmission ψt is a part of ψi and, therefore, affects the δ values. A radar capability to alter ψi by varying ψt could deliver additional information about scattering media. [ABSTRACT FROM AUTHOR]

Two Ocean State Monitoring and Analyzing Radar (OSMAR071) (7.8 MHz) high-frequency (HF) radars and four moored ADCPs were operated concurrently in the southwestern Taiwan Strait during January-March 2013. Qualitative and quantitative comparisons of surface currents were conducted between the HF radars and the ADCPs. Except for a location probably affected by shallow water and sand waves on the Taiwan Banks, the HF-radar-derived radial currents (radials) showed good agreement with the ADCP measured results (correlation coefficient: 0.89-0.98; rms difference: 0.07-0.13 m s-1). To provide further insight into the geophysical processes involved, the performance of the HF-radar-derived radials was further evaluated under different sea states (sea states: 2-6). It was found that both the data returns of the radar-derived radials and the differences between the radar-derived radials and the ADCP-derived radials varied with sea state. The HF radar performed best at sea state 4 in terms of data returns. The spatial coverage increased rapidly as the waves increased from sea state 2 to 4. However, it decreased slowly from sea state 4 to 6. Second, the radial differences were relatively high under lower sea states (2 and 3) at the location where the best agreement was obtained between the radar and ADCP radials, whereas the differences increased as the sea states increased at the other three locations. The differences between the radials measured by the HF radars and the ADCPs could be attributed to wave-induced Stokes drift and spatial sampling differences. [ABSTRACT FROM AUTHOR]

On 10 November 2018, during the RELAMPAGO field campaign in Argentina, South America, a thunderstorm with supercell characteristics was observed by an array of mobile observing instruments, including three Doppler on Wheels radars. In contrast to the archetypal supercell described in the Glossary of Meteorology, the updraft rotation in this storm was rather short lived (;25 min), causing some initial doubt as to whether this indeed was a supercell. However, retrieved 3D winds from dual-Doppler radar scans were used to document a high spatial correspondence between midlevel vertical velocity and vertical vorticity in this storm, thus providing evidence to support the supercell categorization. Additional data collected within the RELAMPAGO domain revealed other storms with this behavior, which appears to be attributable in part to effects of the local terrain. Specifically, the IOP4 supercell and other short-duration supercell cases presented had storm motions that were nearly perpendicular to the long axis of the Sierras de Córdoba Mountains; a long-duration supercell case, on the other hand, had a storm motion nearly parallel to these mountains. Sounding observations as well as model simulations indicate that a mountain-perpendicular storm motion results in a relatively short storm residence time within the narrow zone of terrain-enhanced vertical wind shear. Such a motion and short residence time would limit the upward tilting, by the left-moving supercell updraft, of the storm-relative, antistreamwise horizontal vorticity associated with anabatic flow near complex terrain. [ABSTRACT FROM AUTHOR]

The attenuation-based rainfall estimator is less sensitive to the variability of raindrop size distributions (DSDs) than conventional radar rainfall estimators. For the attenuation-based quantitative precipitation estimation (QPE), the key is to accurately estimate the horizontal specific attenuation AH, which requires a good estimate of the ray-averaged ratio between AH and specific differential phase KDP, also known as the coefficient a. In this study, a variational approach is proposed to optimize the coefficient a for better estimates of AH and rainfall. The performance of the variational approach is illustrated using observations from an S-band operational weather radar with rigorous quality control in south China, by comparing against the α optimization approach using a slope of differential reflectivity ZDR dependence on horizontal reflectivity factor ZH. Similar to the ZDR-slope approach, the variational approach can obtain the optimized a consistent with the DSD properties of precipitation on a sweep-to-sweep basis. The attenuation-based hourly rainfall estimates using the sweep-averaged a values from these two approaches show comparable accuracy when verified against the gauge measurements. One advantage of the variational approach is its feasibility to optimize a for each radar ray, which mitigates the impact of the azimuthal DSD variabilities on rainfall estimation. It is found that, based on the optimized a for radar rays, the hourly rainfall amounts derived from the variational approach are consistent with gauge measurements, showing lower bias (1.0%), higher correlation coefficient (0.92), and lower root-mean-square error (2.35 mm) than the results based on the sweepaveraged α. [ABSTRACT FROM AUTHOR]

Weather radars and gauge observations are the primary observations to determine the coverage and magnitude of precipitation; however, radar and gauge networks have significant coverage gaps, which can underrepresent or even miss the occurrence of precipitation. This is especially noticeable in mountainous regions and in shallow precipitation regimes. The following study presents a methodology to improve spatial representations of precipitation by seamlessly blending multiple precipitation sources within the Multi-Radar Multi-Sensor (MRMS) system. A high spatiotemporal resolution multisensor merged quantitative precipitation estimation (QPE) product (MSQPE) is generated by using gauge-corrected radar QPE as a primary precipitation source with a combination of hourly gauge observations, monthly precipitation climatologies, numerical weather prediction short-term precipitation forecasts, and satellite observations to use in areas of insufficient radar coverage. The merging of the precipitation sources is dependent upon radar coverage based on an updated MRMS radar quality index, surface and atmospheric conditions, topography, gauge locations, and precipitation values. Evaluations of the MSQPE product over the western United States resulted in improved statistical measures over its individual input precipitation sources, particularly the locally gauge-corrected radar QPE. TheMSQPE scheme demonstrated its ability to sufficiently fill in areas where radar alone failed to detect precipitation due to significant beam blockage or poor coverage while minimizing the generation of false precipitation and underestimation biases that resulted from radar overshooting precipitation. [ABSTRACT FROM AUTHOR]

This work describes a phase-resolving wave-forecasting algorithm that is based on the assimilation of marine radar image time series. The algorithm is tested against synthetic data and field observations. The algorithm couples X-band marine radar observations with a phase-resolving wave model that uses the linear mild slope equations for reconstruction of water surface elevations over a large domain of O(km) and a prescribed time window of O(min). The reconstruction also enables wave-by-wave forecasting through forward propagation in space and time. Marine radar image time series provide the input wave observations through a previously given relationship between backscatter intensity and the radial component of the sea surface slope. The algorithm assimilates the wave slope information into the model via a best-fit wave source function at the boundary that minimizes the slope reconstruction error over an annular region at the outer ranges of the radar images. The wave model is then able to propagate the waves across a polar domain to a location of interest at nearer ranges. The constraints on the method for achieving real-time forecasting are identified, and the algorithm is verified against synthetic data and tested using field observations. [ABSTRACT FROM AUTHOR]

The structure of vortex Rossby waves (VRWs) and their role in the development of a secondary eyewall in Hurricane Matthew (2016) is examined from observations taken during the NOAA Sensing Hazards with Operational Unmanned Technology (SHOUT) field experiment. Radar measurements from ground-based and airborne systems, with a focus on the NASA High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) instrument on the Global Hawk aircraft, revealed the presence of; 12-15-km-wavelength spiral bands breaking from the inner-core eyewall in the downshear-right quadrant. The vorticity characteristics and calculations of the intrinsic phase speeds of the bands are shown to be consistent with sheared VRWs. A new angular momentum budget methodology is presented that allows an understanding of the secondary eyewall development process with narrow-swath radar measurements. Filtering of the governing equations enables explicit insight into the nonlinear dynamics of scale interactions and the role of the VRWs in the storm structure change. The results indicate that the large-scale (scales .15 km) vertical flux convergence of angular momentum associated with the VRWs dominates the time tendency with smaller effects from the radial flux term. The small-scale (scales # 15 km) vertical term produces weak, but nonnegligible nonlinear forcing of the large scales primarily through the Reynolds and cross-stress components. The projection of the wave kinematics onto the low-wavenumber (0 and 1) fields appears to be the more significant dynamic process. Flight-level observations show secondary peaks in tangential winds in the radial region where the VRW forcing signatures are active, connecting them with the secondary eyewall formation process. [ABSTRACT FROM AUTHOR]

The sensitivity of microwave brightness temperatures (TBs) to hydrometeors at frequencies between 89 and 190 GHz is investigated by comparing Fengyun-3C(FY-3C)MicrowaveHumidity Sounder-2 (MWHS-2) measurements with radar reflectivity profiles and retrieved products from the Global Precipitation Measurement mission's Dual-Frequency Precipitation Radar (DPR). Scattering-induced TB depressions (DTBs), calculated by subtracting simulated cloud-free TBs from bias-corrected observed TBs for each channel, are compared withDPR-retrieved hydrometeor water path (HWP) and vertically integrated radar reflectivity ZINT. We also account for the number of hydrometeors actually visible in each MWHS-2 channel by weightingHWP with the channel's cloud-free gas transmission profile and the observation slant path. We denote these transmission-weighted, slant-path-integrated quantities with a superscript asterisk (e.g., HWP*). The so-derived linear sensitivity of DTB with respect to HWP* increases with frequency roughly to the power of 1.78. A retrieved HWP* of 1 kgm22 at 89GHz on average corresponds to a decrease in observed TB, relative to a cloud-free background, of 11 K. At 183 GHz, the decrease is about 34-53K. We perform a similar analysis using the vertically integrated, transmission-weighted slant-path radar reflectivity Z*INT and find that DTB also decreases approximately linearly with (Z*INT)0:58. The exponent of 0.58 corresponds to the one we find in the purely DPR-retrieval-based ZINT-HWP relation. The observed sensitivities of DTB with respect to Z*INT and HWP* allow for the validation of hydrometeor scattering models. [ABSTRACT FROM AUTHOR]

Although radar is our most useful tool for monitoring severe weather, the benefits of assimilating its data are often short lived. To understand why, we documented the assimilation requirements, the data characteristics, and the common practices that could hinder optimum data assimilation by traditional approaches. Within storms, radars provide dense measurements of a few highly variable storm outcomes (precipitation and wind) in atmospherically unstable conditions. However, statistical relationships between errors of observed and unobserved quantities often become nonlinear because the errors in these areas tend to become large rapidly. Beyond precipitating areas lie large regions for which radars provide limited new information, yet whose properties will soon shape the outcome of future storms. For those areas, any innovation must consequently be projected from sometimes distant precipitating areas. Thus, radar data assimilation must contend with a double need at odds with many traditional assimilation implementations: correcting in-storm properties with complex errors while projecting information at unusually far distances outside precipitating areas. To further complicate the issue, other data properties and practices, such as assimilating reflectivity in logarithmic units, are not optimal to correct all state variables. Therefore, many characteristics of radar measurements and common practices of their assimilation are incompatible with necessary conditions for successful data assimilation. Facing these dataset-specific challenges may force us to consider new approaches that use the available information differently. [ABSTRACT FROM AUTHOR]

Raindrop size distribution (DSD) observed using a disdrometer can be represented by a constrained-gamma (C-G) DSD model based on the empirical relationship between shape (µ) and slope (Λ). The C-G DSD model can be used to retrieve DSDs and rain microphysical parameters from dual-polarization radar measurements of reflectivity (ZH) and differential reflectivity (ZDR). This study presents a new µ-Λ relationship to characterize rain microphysics in South Korea using a two-dimensional video disdrometer (2DVD) and Yong-in S-band dual-polarization radar. To minimize sampling errors from the 2DVD and radar measurements, measured size distributions are truncated by particle size and velocity-based filtering and compared with rain gauge measurement. The calibration biases of radar ZH and ZDR were calculated using the self-consistency constraint and vertical pointing measurements. The derived µ-Λ relationship was verified using the mass-weighted mean diameter (Dm) and standard deviation of the size distribution (σm), calculated from the 2DVD, for comparison with existing µ-Λ relationships for Florida and Oklahoma. The Dm-σm relationship derived from the 2DVD corresponded well with the µ-Λ relationship. The µ-Λ relationship derived for the Korean Peninsula was similar to Florida, and both generally had larger µ values than Oklahoma for the same Λ. The derived µ-Λ relationship was applied to retrieve DSD parameters from polarimetric radar data, and the retrieved DSDs and derived physical parameters were evaluated and compared with the 2DVD measurements. The polarization radar-based C-G DSD model characterized rain microphysics more accurately than the exponential DSD model. The C-G DSD model based on the newly derived µ-Λ relationship performed the best at retrieving rain microphysical parameters. [ABSTRACT FROM AUTHOR]

This article illustrates how multifrequency radar observations can refine the mass-size parameterization of frozen hydrometeors in scattering models and improve the correlation between the radar observations and in situ measurements of microphysical properties of ice and snow. The data presented in this article were collected during the GPM Cold Season Precipitation Experiment (GCPEx) (2012) and Olympic Mountain Experiment (OLYMPEx) (2015) field campaigns, where the true mass-size relationship was not measured. Starting from size and shape distributions of ice particles measured in situ, scattering models are used to simulate an ensemble of reflectivity factors for various assumed mass-size parameterizations (MSP) of the power-law type. This ensemble is then collocated to airborne and ground-based radar observations, and the MSPs are refined by retaining only those that reproduce the radar observations to a prescribed level of accuracy. A versatile "retrieval dashboard" is built to jointly analyze the optimal MSPs and associated retrievals. The analysis shows that the optimality of an MSP depends on the physical assumptions made in the scattering simulators. This work confirms also the existence of a relationship between parameters of the optimal MSPs. Through the MSP optimization, the retrievals of ice water content M and mean diameter Dm seem robust to the change in meteorological regime (between GCPEx andOLYMPEx); whereas the retrieval of the diameter spread Sm seems more campaign dependent. [ABSTRACT FROM AUTHOR]

Severe thunderstorms and their associated tornadoes pose significant threats to life and property, and using radar data to accurately measure the rotational velocity of circulations in thunderstorms is essential for appropriate, timely warnings. One key factor in accurately measuring circulation velocity is the azimuthal spacing between radar data points, which is referred to as the azimuthal sampling interval. Previous studies have shown that reducing the azimuthal sampling interval can aid in measuring circulation velocity; however, this comes at the price of increased computational complexity. Thus, choosing the best compromise requires knowledge of the relationship between the radar azimuthal sampling interval and the accuracy of the circulation strength as measured from the radar data. In this work, we use simulations to quantify the impact of azimuthal sampling on the strength of radar-observed circulations and show that the improvements get progressively smaller as the azimuthal sampling interval decreases. Thus, improved characterization of circulations can be achieved without using the finest possible sampling grid. We use real data to validate the results of the simulations, which can be used to inform the selection of an appropriate azimuthal sampling interval that balances the accuracy of the radar-observed circulations and computational complexity. [ABSTRACT FROM AUTHOR]

Typhoon Bolaven caused significant damage with severe rainfall all over South Korea, including Cheju Island, which received more than 250 mm in 2 days in August 2012. It was regarded as the most powerful storm to strike the Korean Peninsula in nearly a decade. The rainfall-rate datasets were obtained from S-band radar operated by the Korea Meteorological Administration to be analyzed and compared with the mesoscale Cloud Resolving Storm Simulator (CReSS) model simulation. Multifractal analysis was conducted to understand the structure of the rainfall rate with height in the typhoon system. The radar rainfall data presented with strong intermittency across scales at lower altitudes (1 and 2 km) and a more homogeneous rainfall field at high altitude (5 km) with two parameters (fractal codimension and multifractality index). The statistical scaling moment function and maximal singularities show clear significant differences between radar and the CReSS model. [ABSTRACT FROM AUTHOR]

In a 2018 paper by Bukovčić et al., polarimetric bivariate power-law relations for estimating snowfall rate S and ice water content (IWC), S ⁡ (K DP , Z) = γ K DP α Z β and IWC ⁡ (K DP , Z) = γ 2 K DP α 2 Z β 2 , were developed utilizing 2D video disdrometer snow measurements in Oklahoma. Herein, these disdrometer-based relations are generalized for the range of particle aspect ratios from 0.5 to 0.8 and the width of the canting angle distribution from 0° to 40° and are validated via analytical/theoretical derivations and simulations. In addition, a novel S(KDP, Zdr) polarimetric relation utilizing the ratio between specific differential phase KDP and differential reflectivity Zdr, K DP / ⁡ (1 − Z dr − 1 ) , is derived. Both KDP and (1 − Z dr − 1 ) are proportionally affected by the ice particles' aspect ratio and width of the canting angle distribution; therefore, the variables' ratio tends to be almost invariant to the changes in these parameters. The S(KDP, Z) and S(KDP, Zdr) relations are applied to the polarimetric S-band WSR-88D data obtained from three geographical locations in Virginia, Oklahoma, and Colorado, and their performance is compared with estimations from the standard S(Z) relations and ground snow measurements. The polarimetric estimates of snow accumulations from the three cases exhibit smaller bias in comparison with the S(Z), indicating good potential for more reliable radar snow measurements. [ABSTRACT FROM AUTHOR]

Improving our ability to predict future weather and climate conditions is strongly linked to achieving significant advancements in our understanding of cloud and precipitation processes. Observations are critical to making these advancements because they both improve our understanding of these processes and provide constraints on numerical models. Historically, instruments for observing cloud properties have limited cloud–aerosol investigations to a small subset of cloud-process interactions. To address these challenges, the last decade has seen the U.S. DOE ARM facility significantly upgrade and expand its surveillance radar capabilities toward providing holistic and multiscale observations of clouds and precipitation. These upgrades include radars that operate at four frequency bands covering a wide range of scattering regimes, improving upon the information contained in earlier ARM observations. The traditional ARM emphasis on the vertical column is maintained, providing more comprehensive, calibrated, and multiparametric measurements of clouds and precipitation. In addition, the ARM radar network now features multiple scanning dual-polarization Doppler radars to exploit polarimetric and multi-Doppler capabilities that provide a wealth of information on storm microphysics and dynamics under a wide range of conditions. Although the diversity in wavelengths and detection capabilities are unprecedented, there is still considerable work ahead before the full potential of these radar advancements is realized. This includes synergy with other observations, improved forward and inverse modeling methods, and well-designed data–model integration methods. The overarching goal is to provide a comprehensive characterization of a complete volume of the cloudy atmosphere and to act as a natural laboratory for the study of cloud processes. [ABSTRACT FROM AUTHOR]

Radar and surface rainfall observations are two sources of operational data crucial for heavy rainfall prediction. Their individual values on improving convective forecasting through data assimilation have been examined in the past using convection-permitting numerical models. However, the benefit of their simultaneous assimilations has not yet been evaluated. The objective of this study is to demonstrate that, using a 4D-Var data assimilation system with a microphysical scheme, these two data sources can be assimilated simultaneously and the combined assimilation of radar data and estimated rainfall data from radar reflectivity and surface network can lead to improved short-term heavy rainfall prediction. In our study, a combined data assimilation experiment is compared with a rainfall-only and a radar-only (with or without reflectivity) experiments for a heavy rainfall event occurring in Taiwan during the passage of a mei-yu system. These experiments are conducted by applying the Weather Research and Forecasting (WRF) 4D-Var data assimilation system with a 20-min time window aiming to improve 6-h convective heavy rainfall prediction. Our results indicate that the rainfall data assimilation contributes significantly to the analyses of humidity and temperature whereas the radar data assimilation plays a crucial role in wind analysis, and further, combining the two data sources results in reasonable analyses of all three fields by eliminating large, unphysical analysis increments from the experiments of assimilating individual data only. The results also show that the combined assimilation improves forecasts of heavy rainfall location and intensity of 6-h accumulated rainfall for the case studied. [ABSTRACT FROM AUTHOR]

A detailed damage survey is combined with high-resolution mobile, rapid-scanning X-band polarimetric radar data collected on the Shawnee, Oklahoma, tornado of 19 May 2013. The focus of this study is the radar data collected during a period when the tornado was producing damage rated EF3. Vertical profiles of mobile radar data, centered on the tornado, revealed that the radar reflectivity was approximately uniform with height and increased in magnitude as more debris was lofted. There was a large decrease in both the cross-correlation coefficient (ρhv) and differential radar reflectivity (ZDR) immediately after the tornado exited the damaged area rated EF3. Low ρhv and ZDR occurred near the surface where debris loading was the greatest. The 10th percentile of ρhv decreased markedly after large amounts of debris were lofted after the tornado leveled a number of structures. Subsequently, ρhv quickly recovered to higher values. This recovery suggests that the largest debris had been centrifuged or fallen out whereas light debris remained or continued to be lofted. Range–height profiles of the dual-Doppler analyses that were azimuthally averaged around the tornado revealed a zone of maximum radial convergence at a smaller radius relative to the leading edge of lofted debris. Low-level inflow into the tornado encountering a positive bias in the tornado-relative radial velocities could explain the existence of the zone. The vertical structure of the convergence zone was shown for the first time. [ABSTRACT FROM AUTHOR]

A well-designed deployment of well-maintained surface instruments as well as abundant rainfall provided an excellent dataset with which to evaluate the Micro Rain Radar (MRR) performance for estimating raindrop size distribution (DSD) and its integral rainfall parameters with respect to the consolidated devices during the Iowa Flood Studies (IFloodS) field campaign. The MRR was collocated with two-dimensional video disdrometer (2DVD) and Autonomous Parsivel2 Unit (APU) at three different sites located at 5–70-km distances from the National Aeronautics and Space Administration's S-band dual-polarization Doppler radar (NPOL). A comparative study between MRR, 2DVD, APU, and NPOL was conducted including all rainy minutes as well as minutes of stratiform rain and convective rain. Considering 2DVD as a primary reference, a good agreement was evident for reflectivity between MRR's lowest reliable height and 2DVD with an absolute bias of less than 2 dB even in convective rain except for one site. For rainfall rate, the percent absolute bias between MRR and 2DVD ranged between 25% and 35% in stratiform rain and about 10% higher in convective rain. Agreement for mean mass-weighted raindrop diameter was good (bias less than 0.1 mm), whereas MRR overestimated the normalized intercept parameter of the gamma DSD [mean bias among the three sites was −0.13 log(mm−1 m−3)]. The agreement between MRR and APU was slightly worse than the one between MRR and 2DVD. When the horizontal and differential reflectivities of NPOL were compared with the ones derived from the MRR DSD resampled within the radar volume, we found an absolute bias of approximately 3 and 0.4 dB, respectively. [ABSTRACT FROM AUTHOR]

Quasi-vertical profiles (QVPs) obtained from a database of U.S. WSR-88D data are used to document polarimetric characteristics of the melting layer (ML) in cold-season storms with high vertical resolution and accuracy. A polarimetric technique to define the top and bottom of the ML is first introduced. Using the QVPs, statistical relationships are developed to gain insight into the evolution of microphysical processes above, within, and below the ML, leading to a statistical polarimetric model of the ML that reveals characteristics that reflectivity data alone are not able to provide, particularly in regions of weak reflectivity factor at horizontal polarization ZH. QVP ML statistics are examined for two regimes in the ML data: ZH ≥ 20 dBZ and ZH < 20 dBZ. Regions of ZH ≥ 20 dBZ indicate locations of MLs collocated with enhanced differential reflectivity ZDR and reduced copolar correlation coefficient ρhv, while for ZH < 20 dBZ a well-defined ML is difficult to discern using ZH alone. Evidence of large ZDR up to 4 dB, backscatter differential phase δ up to 8°, and low ρhv down to 0.80 associated with lower ZH (from −10 to 20 dBZ) in the ML is observed when pristine, nonaggregated ice falls through it. Positive correlation is documented between maximum specific differential phase KDP and maximum ZH in the ML; these are the first QVP observations of KDP in MLs documented at S band. Negative correlation occurs between minimum ρhv in the ML and ML depth and between minimum ρhv in the ML and the corresponding enhancement of ZH (ΔZH = ZHmax − ZHrain). [ABSTRACT FROM AUTHOR]

To fulfill the evolving observational needs of the National Weather Service (NWS), future weather radar systems will have to meet demanding requirements. Designing such systems will likely involve trade-offs between system cost and operational performance. A potential cost driver for future weather radars that could cause significant data-quality impacts on forecasters is the required angular resolution and sidelobe performance, which are mainly dictated by the antenna radiation pattern. Typical antenna radiation patterns can be characterized by the width of the main lobe and their sidelobe levels, which are traditionally measured across the azimuthal and elevation dimensions. In this work, we study the impact of increasing sidelobe levels on NWS forecasters' data interpretation during warning operations. The resulting impact model can be used by decision-makers to better understand the cost–benefit trade-offs inherent in any radar system design. [ABSTRACT FROM AUTHOR]

Supercell thunderstorms produce a variety of hazards, including tornadoes. A supercell will often exist for some time prior to producing a tornado, while other supercells never become tornadic. In this study, a series of hypotheses is tested regarding the ability of S-band polarimetric radar fields to distinguish pretornadic from nontornadic supercell storms. Several quantified polarimetric radar metrics are examined that are related to storm inflow, updraft, and hailfall characteristics in samples of 19–30 pretornadic and 18–31 nontornadic supercells. The results indicate that pretornadic supercells are characterized by smaller hail extent and echo appendages with larger mean drop size. Additionally, differential reflectivity ZDR column size is larger and less variable in the pretornadic storms in the 25–30 min prior to initial tornadogenesis. Many of the results indicate relatively small polarimetric differences that will likely be difficult to translate to operational use. Hail extent and ZDR column size, however, may exhibit operationally useful differences between pretornadic and nontornadic supercells. [ABSTRACT FROM AUTHOR]

Radar reflectivity (Z) data are either directly assimilated using 3DVar, 4DVar, or ensemble Kalman filter, or indirectly assimilated using, for example, cloud analysis that preretrieves hydrometeors from Z. When directly assimilating radar data variationally, issues related to the highly nonlinear Z operator arise that can cause nonconvergence and bad analyses. To alleviate the issues, treatments are proposed in this study and their performances are examined via observing system simulation experiments. They include the following: 1) When using hydrometeor mixing ratios as control variables (CVq), small background Z can cause extremely large cost function gradient. Lower limits are imposed on the mixing ratios (qLim treatment) or the equivalent reflectivity (ZeLim treatment) in Z observation operator. ZeLim is found to work better than qLim in terms of analysis accuracy and convergence speed. 2) With CVq, the assimilation of radial velocity (Vr) is ineffective when assimilated together with Z data due to the much smaller cost function gradient associated with Vr. A procedure (VrPass) that assimilates Vr data in a separate pass is found very helpful. 3) Using logarithmic hydrometeor mixing ratios as control variables (CVlogq) can also avoid extremely large cost function gradient, and has much faster convergence. However, spurious analysis increments can be created when transforming the analysis increments back to mixing ratios. A background smoothing and a lower limit are applied to the background mixing ratios, and are shown to be effective. Using CVlogq with associated treatments produces better reflectivity analysis that is much closer to the observation without resorting to multiple analysis passes, and the cost function minimization also converges faster. CVlogq is therefore recommended for variational radar data assimilation. [ABSTRACT FROM AUTHOR]

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Zdr and larger specific differential phase Kdp above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Zdr and Kdp observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Zdr and Kdp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics. [ABSTRACT FROM AUTHOR]

The ground-based-radar-derived raindrop size distribution (DSD) parameters—mass-weighted drop diameter Dmass and normalized intercept parameter NW—are the sole resource for direct validation of the National Aeronautics and Space Administration (NASA) Global Precipitation Measurement (GPM) mission Core Observatory satellite-based retrieved DSD. Both Dmass and NW are obtained from radar-measured reflectivity ZH and differential reflectivity ZDR through empirical relationships. This study uses existing relationships that were determined for the GPM ground validation (GV) program and directly compares the NASA S-band polarimetric radar (NPOL) observables of ZH and ZDR and derived Dmass and NW with those calculated by two-dimensional video disdrometer (2DVD). The joint NPOL and 2DVD datasets were acquired during three GPM GV field campaigns conducted in eastern Iowa, southern Appalachia, and western Washington State. The comparative study quantifies the level of agreement for ZH, ZDR, Dmass, and log(NW) at an optimum distance (15–40 km) from the radar as well as at distances greater than 60 km from radar and over mountainous terrain. Interestingly, roughly 10%–15% of the NPOL ZH–ZDR pairs were well outside the envelope of 2DVD-estimated ZH–ZDR pairs. The exclusion of these pairs improved the comparisons noticeably. [ABSTRACT FROM AUTHOR]

The potential future installation of a multifunction phased-array radar (MPAR) network will provide capabilities of case-specific adaptive scanning. Knowing the impacts adaptive scanning may have on short-term forecasts will influence scanning strategy decision-making in hopes to produce the most optimal ensemble forecast while also benefiting human severe weather warning decision-making. An ensemble-based targeted observation algorithm is applied to an observing system simulation experiment (OSSE) where the impacts of synthetic idealized supercell radial velocity observations are estimated before the observations are "collected" and assimilated. The forecast metric of interest is the low-level rotation forecast metric (0–1-km updraft helicity), a surrogate for tornado prediction. It is found that the ensemble-based targeted observation approach can reasonably estimate the true error variance reduction when an effective method that treats sampling error is applied, the period of model forecast is associated with less degrees of nonlinearity, and the observation information content relative to the background forecast is larger. In some scenarios, a subset of a full-volume scan assimilation produces better forecasts than all observations within the full volume. Assimilating the full-volume scan increases the number of potential spurious correlations arising between the forecast metric and radial velocity observation induced state perturbations, which may degrade the forecast metric accuracy. [ABSTRACT FROM AUTHOR]

Severe weather nowcasting is a crucial mission of atmospheric science for the betterment of society to save life, limb, and property. In this study, composite radar data from the Central Weather Bureau of 16 typhoons are collected to examine the statistical performance of the McGill Algorithm for Precipitation nowcasting using Lagrangian Extrapolation (MAPLE) over Taiwan, an extrapolation algorithm that predicts future precipitation based on current radar echoes. In addition, instead of mixing the precipitation between radar extrapolation and numerical model forecast as in previous studies, a blending system is formed by synthesizing the wind information from model forecast with the echo extrapolation motion field via a variational algorithm to improve the nowcasting system. The statistical results of the radar echo extrapolation for 16 typhoon cases show that while the quantitative precipitation nowcasting skill can persist for up to 2 h, significant distortion for the rotational system is found after 2 h. On the other hand, the blending system helps to capture and maintain the rotation of typhoon rainband structures. The blending system extends the nowcasting skill by 1 h to a total of 3 h. Furthermore, the blending scheme performs especially well after the typhoon makes landfall in Taiwan. For disaster prevention and mitigation, this blending nowcasting technique may provide effective weather information immediately. [ABSTRACT FROM AUTHOR]

The performance of a numerical weather prediction model using convective-scale ensemble data assimilation with ground-based global navigation satellite systems-zenith total delay (ZTD) and radar data is investigated on a heavy rainfall event that occurred in Taiwan on 10 June 2012. The assimilation of ZTD and/or radar data is performed using the framework of the WRF local ensemble transform Kalman filter with a model grid spacing of 2 km. Assimilating radar data is beneficial for predicting the rainfall intensity of this local event but produces overprediction in southern Taiwan and underprediction in central Taiwan during the first 3 h. Both errors are largely overcome by assimilating ZTD data to improve mesoconvective-scale moisture analyses. Consequently, assimilating both the ZTD and radar data show advantages in terms of the location and intensity of the heavy rainfall. Sensitivity experiments involving this event indicate that the impact of ZTD data is improved by using a broader horizontal localization scale than the convective scale used for radar data assimilation. This optimization is necessary in order to consider more fully the network density of the ZTD observations and the horizontal scale of the moisture transport by the southwesterly flow in this case. [ABSTRACT FROM AUTHOR]

In radar observations of hydrometeors, the 0°C isotherm in the atmosphere (i.e., the freezing level) usually appears as a region of enhanced reflectivity. This region is known as the bright band (BB). In this study, observations over 12 months from a vertically pointing 35-GHz radar and a collocated disdrometer at the Natural Environment Research Council (NERC) Facility for Atmospheric and Radio Research (NFARR) are used to identify and compare microphysical differences between BB and non-brightband (NBB) periods. From these observations, the relationship between radar reflectivity Z and rainfall intensity R is found to be Z = 772R0.57 for BB periods and Z = 108R0.99 for NBB periods. Additionally, the brightband strength (BBS) was calculated using a novel method derived from the Michelson contrast equation in an attempt to explain the observed variability in BB precipitation. A series of Z–R relationships are computed with respect to BBS. The coefficients increase with increasing BBS from 227 to 926, while the exponents decrease with increasing BBS from 0.85 to 0.38. The results also indicate that NBB periods identified in the presence of a 0°C isotherm in other studies may be misclassified due to their inability to identify weak brightband periods. As such, it is hypothesized that NBB periods are solely due to warm rain processes. [ABSTRACT FROM AUTHOR]

A potential replacement candidate for the aging operational WSR-88D infrastructure currently in place is the phased array radar (PAR) system. The current WSR-88Ds take ~5 min to produce a full volumetric scan of the atmosphere, whereas PAR technology allows for full volumetric scanning of the same atmosphere every ~1 min. How this increase in temporal frequency of radar observations might affect the National Severe Storms Laboratory's (NSSL) Warn-on-Forecast system (WoFS), which is a storm-scale ensemble data assimilation and forecast system for severe convective weather, is unclear. Since radar data assimilation is critical for the WoFS, this study explores the optimal temporal frequency of PAR observations for storm-scale data assimilation using the 31 May 2013 El Reno, Oklahoma, tornadic supercell event. The National Severe Storms Laboratory's National Weather Radar Testbed PAR in Norman, Oklahoma, began scanning this event more than an hour before the first (and strongest) tornado developed near El Reno, and scanned most of the tornadic supercell's evolution. Several experiments using various cycling and data frequencies to synchronously and asynchronously assimilate these PAR observations are conducted to produce analyses and very short-term forecasts of the El Reno supercell. Forecasts of low-level reflectivity and midlevel updraft helicity are subjectively evaluated and objectively verified using spatial and object-based techniques. Results indicate that assimilating more frequent PAR observations can lead to more accurate analyses and probabilistic forecasts of the El Reno supercell at longer lead times. Hence, PAR is a promising radar platform for WoFS. [ABSTRACT FROM AUTHOR]

Radar retrievals of drop size distribution (DSD) parameters are developed and evaluated over the mountainous Olympic Peninsula of Washington State. The observations used to develop retrievals were collected during the 2015/16 Olympic Mountain Experiment (OLYMPEX) and included the NASA S-band dual-polarimetric (NPOL) radar and a collection of second-generation Particle Size and Velocity (PARSIVEL2) disdrometers over the windward slopes of the barrier. Nonlinear and random forest regressions are applied to the PARSIVEL2 data to develop retrievals for median volume diameter, liquid water content, and rain rate. Improvement in DSD retrieval accuracy, defined by the mean error of the retrieval relative to PARSIVEL2 observations, was achieved when using the random forest model when compared with nonlinear regression. Evaluation of disdrometer observations and the retrievals from NPOL indicate that the radar retrievals can accurately reproduce observed DSDs in this region, including the common wintertime regime of small but numerous raindrops that is important there. NPOL retrievals during the OLYMPEX period are further evaluated using two-dimensional video disdrometers (2DVD) and vertically pointing Micro Rain Radars. Results indicate that radar retrievals using random forests may be skillful in capturing DSD characteristics in the lowest portions of the atmosphere. [ABSTRACT FROM AUTHOR]

This study utilizes data collected by the University of Oklahoma Advanced Radar Research Center's Polarimetric Radar for Innovations in Meteorology and Engineering (OU-PRIME) C-band radar as well as the federal KTLX and KOUN WSR-88D S-band radars to study a supercell that simultaneously produced a long-track EF-4 tornado and an EF-2 landspout tornado (EF indicates the enhanced Fujita scale) near Norman, Oklahoma, on 10 May 2010. Contrasting polarimetric characteristics of two tornadoes over similar land cover but with different intensities are documented. Also, the storm-scale sedimentation of debris within the supercell is investigated, which includes observations of rotation and elongation of a tornadic debris signature with height. A dual-wavelength comparison of debris at S and C bands is performed. These analyses indicate that lofted debris within the tornado was larger than debris located outside the damage path of the tornado and that debris size outside the tornado increased with height, likely as the result of centrifuging. Profiles of polarimetric variables were observed to become more vertically homogeneous with time. [ABSTRACT FROM AUTHOR]

A dataset of convection initiation (CI) is of great value in studying the triggering mechanisms of deep moist convection and evaluating the performances of numerical models. In recent years, the data quality of the operationally generated radar mosaics over China has been greatly improved, which provides an opportunity to retrieve a CI dataset from that region. In this work, an attempt is made to reveal the potential of applying a simple framework of objective CI detection for the study of CI climatology in China. The framework was tested using radar mosaic maps in South China that were accessible online. The identified CI events were validated in both direct and indirect ways. On the basis of a direct manual check, nearly all of the identified CI cells had an organized motion. The precipitation echoes of the cells had a median duration of approximately 2.5 h. The CI occurrences were further compared with rainfall estimates to ensure physical consistency. The diurnal cycle of CI occurrence exhibits three major modes: a late-night-to-morning peak at the windward coasts and offshore, a noon-to-late-afternoon peak on the coastal land, and an evening-to-early-morning peak over the northwestern highland. These spatial modes agree well with those of rainfall, indirectly suggesting the reliability of the CI statistics. By processing radar mosaic maps, such a framework could be applied for studying CI climatology over China and other regions. [ABSTRACT FROM AUTHOR]

This case study analyzes a nocturnal mesoscale convective system (MCS) that was observed on 25-26 June 2015 in northeastern Kansas during the Plains Elevated Convection At Night (PECAN) project. Over the course of the observational period, a broken line of elevated nocturnal convective cells initiated around 0230UTCon the cool side of a stationary front and subsequently merged to form a quasi-linearMCSthat later developed strong, surface-based outflow and a trailing stratiform region. This study combines radar observations with mobile and fixed mesonet and sounding data taken during PECAN to analyze the kinematics and thermodynamics of the MCS from 0300 to 0630 UTC. This study is unique in that 38 consecutive multi- Doppler wind analyses are examined over the 3.5 h observation period, facilitating a long-duration analysis of the kinematic evolution of the nocturnal MCS. Radar analyses reveal that the initial convective cells and linear MCS are elevated and sustained by an elevated residual layer formed via weak ascent over the stationary front. During upscale growth, individual convective cells develop storm-scale cold pools due to pockets of descending rear-to-front flow that are measured by mobile mesonets. By 0500 UTC, kinematic analysis and mesonet observations show that the MCS has a surface-based cold pool and that convective line updrafts are ingesting parcels from below the stable layer. In this environment, the elevated system has become surface based since the cold pool lifting is sufficient for surface-based parcels to overcome the CIN associated with the frontal stable layer. [ABSTRACT FROM AUTHOR]

We investigate the practical predictability limits of deep convection in a state-of-the-art, high-resolution, limited-area ensemble prediction system. A combination of sophisticated predictability measures, namely, believable and decorrelation scale, are applied to determine the predictable scales of short-term forecasts in a hierarchy of model configurations. First, we consider an idealized perfect model setup that includes both small-scale and synoptic-scale perturbations. We find increased predictability in the presence of orography and a strongly beneficial impact of radar data assimilation, which extends the forecast horizon by up to 6 h. Second, we examine realistic COSMO-KENDAsimulations, including assimilation of radar and conventional data and a representation of model errors, for a convectively active two-week summer period over Germany. The results confirm increased predictability in orographic regions. We find that both latent heat nudging and ensemble Kalman filter assimilation of radar data lead to increased forecast skill, but the impact is smaller than in the idealized experiments. This highlights the need to assimilate spatially and temporally dense data, but also indicates room for further improvement. Finally, the examination of operational COSMO-DE-EPS ensemble forecasts for three summer periods confirms the beneficial impact of orography in a statistical sense and also reveals increased predictability in weather regimes controlled by synoptic forcing, as defined by the convective adjustment time scale. [ABSTRACT FROM AUTHOR]

The statistical properties of the radar echoes from biological, precipitation, and ground targets observed with the McGill S-band dual-polarization radar have been used to devise a polarimetric and a nonpolarimetric fuzzy logic algorithm for pixel-by-pixel target identification. Radar observations of migrating birds show distinctly different polarimetric features during their relative approach and departure from the radar site illustrating the dependency of radar parameters on the canting angle and scattering cross section. The devised algorithms have been tested with two independent events, each consisting of 2 h of radar observations with a 5-min temporal resolution. One event consisted of precipitation without birds while the other contained only birds. The misclassifications were 10.12% and 9.6%, respectively, for the two cases for the nonpolarimetric algorithm, and 1.99% and 0.92% for the polarimetric algorithm. The results indicate that even though nonpolarimetric radar membership functions may be considered adequate for separating radar echo returns from birds, precipitation, and ground targets, they are not sufficiently skilled if a greater accuracy is required. Target identification without polarimetric variables especially fails in the region of zero isodop and in precipitation with an echo top below 4 km. [ABSTRACT FROM AUTHOR]

Hurricane Harvey hit the Texas Gulf Coast as a major hurricane on 25 August 2017 before exiting the state as a tropical storm on 29 August 2017. Left in itswake was historic flooding, with some locations measuringmore than 60 in. (150 cm) of rain over a 5-day period. TheWSR-88D radar (KHGX)maintained operations for the entirety of the event. Rain gauge data fromtheHarris County Flood Warning System (HCFWS) was used for validation with the full radar dataset to retrieve daily and event-total precipitation estimates for the period 25-29 August 2017. The KHGX precipitation estimates were then compared with the HCFWS gauges. Three different hybrid polarimetric rainfall retrievals were used, alongwith attenuation-based retrieval that employs the radar-observed differential propagation. An advantage of using a attenuation-based retrieval is its immunity to partial beam blockage and calibration errors in reflectivity and differential reflectivity. All of the retrievals are susceptible to changes in the observed drop size distribution (DSD). No in situDSDdata were available over the study area, so changes in the DSD were interpreted by examining the observed radar data.We examined the parameter space of two key values in the attenuation retrieval to test the sensitivity of the rain retrieval. Selecting a value of a 5 0.015 and b 5 0.600 provided the best overall results, relative to the gauges, but more work needs to be done to develop an automated technique to account for changes in the ambient DSD. [ABSTRACT FROM AUTHOR]

A variational retrieval of rain microphysics from polarimetric radar data (PRD) has been developed through the use of S-band parameterized polarimetric observation operators. Polarimetric observations allow for the optimal retrieval of cloud and precipitation microphysics for weather quantification and data assimilation for convective-scale numerical weather prediction (NWP) by linking PRD to physical parameters. Rain polarimetric observation operators for reflectivity ZH, differential reflectivity ZDR, and specific differential phase KDP were derived for S-band PRD using T-matrix scattering amplitudes. These observation operators link the PRD to the physical parameters of water contentWand mass-/volume-weighted diameter Dm for rain, which can be used to calculate other microphysical information. The S-band observation operators were tested using a 1D variational retrieval that uses the (nonlinear) Gauss--Newton method to iteratively minimize the cost function to find an optimal estimate of Dm and W separately for each azimuth of radar data, which can be applied to a plan position indicator (PPI) radar scan (i.e., a single elevation). Experiments on two-dimensional video disdrometer (2DVD) data demonstrated the advantages of including FDP observations and using the nonlinear solution rather than the (linear) optimal interpolation (OI) solution. PRDcollected by the Norman, Oklahoma(KOUN)WSR-88D on 15 June 2011 were used to successfully test the retrieval method on radar data. The successful variational retrieval from the2DVDand the radar data demonstrate the utility of the proposed method. [ABSTRACT FROM AUTHOR]

A Bayesian Markov chain Monte Carlo (MCMC) algorithm is utilized to compare the skill of an A-Trainlike observing system with a cloud, convection, and precipitation (CCP) observing system like that contemplated for the 2020s by the 2017 National Academy of Sciences Decadal Survey. The main objective is to demonstrate a framework for observational trade space studies. This initial work focuses on weakly precipitating warm shallow cumulus constructed from in situ data. Radiative computations are based on Mie theory with spherical assumptions. Simulated measurements in the CCP configuration consist of W- and Kaband radar reflectivity and path-integrated attenuation, 31 and 94 GHz brightness temperatures (Tb), and visible and near-infrared reflectances. The collection of measurements in the CloudSat configuration is identical, but includes a single 94 GHz radar frequency, and the uncertainty in the 94 GHz microwave brightness temperature is increased to mimic the CloudSat Tb product. The experiments demonstrate that it remains a challenge to diagnose cloud properties in the presence of light rain because of the tendency of microwave remote sensing to respond to the higher moments of the hydrometeor populations. Rain properties are significantly better constrained than cloud properties, even in the optimal CCP configuration. The addition of Ka-band measurements places substantial constraints on the precipitation rain effective radius and rain rates. The Tb offers important information regarding the column-integrated condensate mass, the measurement accuracy of which appears more likely to affect the retrievals of clouds with low liquid water path. The constraints provided by reflectances are largely restricted to regions near the cloud top, particularly in the raining cases. [ABSTRACT FROM AUTHOR]

Radar and rawinsonde data from four ground-based observing stations in the tropical Indo-Pacific warm pool were used to identify possible associations of environmental state variables and their vertical profiles with radar-derived rain rate inside a mesoscale radar domain when the column-integrated relative humidity (CRH) exceeds 80%. At CRH exceeding 80%, a wide range--from near 0 to ;50mmday21--in rain rate is observed; therefore, tropospheric moisture was a necessary but insufficient condition for deep convection. This study seeks to identify possible factors that inhibit rainfall when the atmosphere is sufficiently moist to support large precipitation rates. The domain-mean rain rate was highly sensitive to the areal coverage of intense, convective rainfall that occurs. There were two fundamentally different instances in which convective area was low. One was when the radar domain is primarily occupied by weakly precipitating, stratiform echoes. The other was when the radar domain contained almost no precipitating echoes of any type. While the former was dependent upon the stage of the convective life cycle seen by radar, the latter was probably dependent upon the convective environment. Areal coverage of convective echoes was largely determined by the number of individual convective echoes rather than their sizes, so changes in the clear-air environment of updrafts might have governed how many updrafts grew into deep cumulonimbi. The most likely environmental influence on convective rainfall identified using rawinsonde data was 900--700-hPa lapse rate; however, processes occurring on spatial scales smaller than a radar domain were probably also important but not investigated. [ABSTRACT FROM AUTHOR]

Remote sensing observations, especially those from ground-based radars, have been used extensively to discriminate between severe and nonsevere storms. Recent upgrades to operational remote sensing networks in the United States have provided unprecedented spatial and temporal sampling to study such storms. These networks help forecasters subjectively identify storms capable of producing severe weather at the ground; however, uncertainties remain in how to objectively identify severe thunderstorms using the same data. Here, three large-area datasets (geostationary satellite, ground-based radar, and ground-based lightning detection) are used over 28 recent events in an attempt to objectively discriminate between severe and nonsevere storms, with an additional focus on severe storms that produce tornadoes. Among these datasets, radar observations, specifically those at mid- and upper levels (altitudes at and above 4 km), are shown to provide the greatest objective discrimination. Physical and kinematic storm characteristics from all analyzed datasets imply that significantly severe [≥2-in. (5.08 cm) hail and/or ≥65-kt (33.4 m s-1) straight-line winds] and tornadic storms have stronger upward motion and rotation than nonsevere and less severe storms. In addition, these metrics are greatest in tornadic storms during the time in which tornadoes occur. [ABSTRACT FROM AUTHOR]

Descriptions of the experimental design and research highlights obtained from a series of four multiagency field projects held near Cape Canaveral, Florida, are presented. The experiments featured a 3 MW, dual-polarization, C-band Doppler radar that serves in a dual capacity as both a precipitation and cloud radar. This duality stems from a combination of the radar's high sensitivity and extremely small-resolution volumes produced by the narrow 0.22° beamwidth and the 0.543 m along-range resolution. Experimental highlights focus on the radar's real-time aircraft tracking capability as well as the finescale reflectivity and eddy structure of a thin nonprecipitating stratus layer. Examples of precipitating storm systems focus on the analysis of the distinctive and nearly linear radar reflectivity signatures (referred to as "streaks") that are caused as individual hydrometeors traverse the narrow radar beam. Each streak leaves a unique radar reflectivity signature that is analyzed with regard to estimating the underlying particle properties such as size, fall speed, and oscillation characteristics. The observed along-streak reflectivity oscillations are complex and discussed in terms of diameter-dependent drop dynamics (oscillation frequency and viscous damping time scales) as well as radar-dependent factors governing the near-field Fresnel radiation pattern and inferred drop–drop interference. [ABSTRACT FROM AUTHOR]

This study explores the benefits of assimilating infrared (IR) brightness temperature (BT) observations from geostationary satellites jointly with radial velocity (Vr) and reflectivity (Z) observations from Doppler weather radars within an ensemble Kalman filter (EnKF) data assimilation system to the convection-allowing ensemble analysis and prediction of a tornadic supercell thunderstorm event on 12 June 2017 across Wyoming and Nebraska. While radar observations sample the three-dimensional storm structures with high fidelity, BT observations provide information about clouds prior to the formation of precipitation particles when instorm radar observations are not yet available and also provide information on the environment outside the thunderstorms. To better understand the strengths and limitations of each observation type, the satellite and Doppler radar observations are assimilated separately and jointly, and the ensemble analyses and forecasts are compared with available observations. Results show that assimilating BT observations has the potential to increase the forecast and warning lead times of severe weather events compared with radar observations and may also potentially complement the sparse surface observations in some regions as revealed by the probabilistic prediction of mesocyclone tracks initialized from EnKF analyses as various times. Additionally, the assimilation of both BT and Vr observations yields the best ensemble forecasts, providing higher confidence, improved accuracy, and longer lead times on the probabilistic prediction of midlevel mesocyclones. [ABSTRACT FROM AUTHOR]

During spring 2016 and spring 2017, a vertically pointing, S-band Frequency Modulated Continuous Wave radar (UMass FMCW) was deployed in northern Alabama under the auspices of the Verification of the Origins of Rotation in Tornadoes Experiment (VORTEX)-Southeast. In total, ~14 weeks of data were collected, in conditions ranging from quiescent clear skies to severe thunderstorms. The principal objective of these deployments was to characterize the boundary layer evolution near the VORTEX-Southeast domain. In this paper, we describe intermediate results in service of this objective. Specifically, we describe updates to the UMass FMCW system, document its deployments for VORTEX-Southeast, and apply four automated algorithms: 1) a dealiasing algorithm to the Doppler velocities, 2) a fuzzy logic scatterer classification scheme to separate precipitation from nonprecipitation observations, 3) a brightband/melting-layer identification algorithm for stratiform precipitation, and 4) an extended Kalman filter--based convective boundary layer depth (mixing height) measurement algorithm for nonprecipitation observations. Results from the latter two applications are qualitatively verified against retrieved soundings from a collocated thermodynamic profiling system. [ABSTRACT FROM AUTHOR]

This study uses a dense rain gauge network, radar observations, and an upslope model to document the detailed aspects of precipitation over Da-Tun Mountain (DT) of northern Taiwan under prevailing northeasterly monsoonal flow. DT is one of the most concentrated areas of heavy rainfall in Taiwan, with a size of ~15 km and a terrain peak of ~1 km (MSL), and it constitutes a concave ridge with two windward ridge arms that encompass a funnel-shaped valley. Twenty-one rainfall events identified over DT from January 2011 to February 2015 are chosen for analysis. More than half of the studied cases exhibit two local maxima of rainfall over the ridge arms, and asymmetric characteristics between these two maxima are evident. The other frequent rainfall pattern, observed as upstream winds are more from the east-northeast, is characterized by a local maximum of precipitation inside the valley. Analyses from the upslope model confirm that the occurrence of the windward maxima of rainfall is primarily caused by upslope lifting, and their asymmetric characteristics are closely related to the difference in the azimuthal variations of slope steepness between the two ridge arms. The observed characteristics of rainfall intensities inside the valley, however, cannot be well described by the upslope model. It is found that the lateral flow confluence induced by the deflected flows over the ridge arms may play an essential role in intensifying upslope-forced precipitation within the valley. This effect emerges as upstream winds are roughly parallel to the central axis of funneling regions located between the ridge arms. [ABSTRACT FROM AUTHOR]