Search

In this study, we examined the variations of precipitation morphology and rainfall in relation to the simultaneous passages of a Madden–Julian oscillation (MJO) event and convectively coupled equatorial waves (CCEWs) observed during the Years of the Maritime Continent pilot study. We utilized globally merged infrared brightness temperature data and the radiosonde and radar data observed aboard the Research Vessel Mirai at 4°4′S, 101°54′E. As well as the observed MJO event, equatorial Rossby waves (ERWs), Kelvin waves (KWs), and mixed Rossby–gravity waves (MRGWs) were identified. The radar data exhibited high-frequency variation, mainly caused by KWs and MRGWs, and low-frequency variation, mainly caused by the MJO and ERWs. The MRGWs predominantly modulated convective echo areas and both convective and stratiform volumetric rainfall. In contrast, the MJO event had little influence on the variance of convective echoes. Moreover, stratiform echo areas and volumetric rainfall were more strongly modulated by the combined effects of the MJO, ERWs, KWs, and MRGWs than their convective counterparts. The intense development of stratiform echo areas and volumetric rainfall was coherent with the superimposition of the active phases of the MJO event and all the analyzed CCEWs. The strongest development and a significant reduction of convective echo-top heights before and after the peak MJO date, respectively, were coherent with the passages of ERWs and MRGWs, which were the dominant wave types in modulating echo-top heights. Thus, it appears that the superimposition of the CCEWs on the MJO event exerted complex modulations on the convective activities within the MJO event.

Comparisons of precipitation between general circulation models (GCMs) and observations are often confounded by a mismatch between model output and instrument measurements, including variable type and temporal and spatial resolutions. To mitigate these differences, the radar-simulator Quickbeam within the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) simulates reflectivity from model variables at the subgrid scale. This work adapts Quickbeam to the dual-frequency precipitation radar (DPR) on board the Global Precipitation Measurement (GPM) satellite. The longer wavelength of the DPR is used to evaluate moderate to heavy precipitation in GCMs, which is missed when Quickbeam is used as a cloud radar simulator. Latitudinal and land–ocean comparisons are made between COSP output from the Community Atmospheric Model version 5 (CAM5) and DPR data. Additionally, this work improves the COSP subgrid algorithm by applying a more realistic, nondeterministic approach to assigning GCM gridbox convective cloud cover when convective cloud is not provided as a model output. Instead of assuming a static 5% convective cloud coverage, DPR convective precipitation coverage is used as a proxy for convective cloud coverage. For example, DPR observations show that convective rain typically only covers about 1% of a 2° grid box, but that the median convective rain area increases to over 10% in heavy rain cases. In our CAM5 tests, the updated subgrid algorithm improved the comparison between reflectivity distributions when the convective cloud cover is provided versus the default 5% convective cloud-cover assumption.

This article summarizes research and risk reduction that will inform acquisition decisions regarding NOAA’s future national operational weather radar network. A key alternative being evaluated is polarimetric phased-array radar (PAR). Research indicates PAR can plausibly achieve fast, adaptive volumetric scanning, with associated benefits for severe-weather warning performance. We assess these benefits using storm observations and analyses, observing system simulation experiments, and real radar-data assimilation studies. Changes in the number and/or locations of radars in the future network could improve coverage at low altitude. Analysis of benefits that might be so realized indicates the possibility for additional improvement in severe-weather and flash-flood warning performance, with associated reduction in casualties. Simulations are used to evaluate techniques for rapid volumetric scanning and assess data quality characteristics of PAR. Finally, we describe progress in developing methods to compensate for polarimetric variable estimate biases introduced by electronic beam-steering. A research-to-operations (R2O) strategy for the PAR alternative for the WSR-88D replacement network is presented.

Hydrological hazards usually occur after heavy precipitation, especially during strong convection. Therefore, accurately identifying convective precipitation is very helpful for hydrological warning and forecasting. However, separating the convective, bright band (BB), and stratiform precipitation is found to be challenging when the convection is adjacent to or within the BB region. A new convection/BB/stratiform precipitation segregation algorithm is proposed in this study to resolve this challenging issue. This algorithm is applicable for a single radar volume scan data in native (polar) coordinates and consists of four processes: 1) check the freezing (0°C) level to roughly assess whether convection is occurring or not; 2) identify the convective cores through analyzing composite reflectivity (maximum reflectivity for a given range gate among all the sweeps), vertically integrated liquid water (VIL), VIL horizontal gradient, and reflectivity at the levels of 0°, −10°, and above −10°C; 3) delineate the whole convective region through the seeded region growing method by taking account of the microphysical differences between the BB and convective regions; and 4) delineate BB features in the stratiform region. The proposed algorithm utilizes the physical characteristics of different precipitation types for precisely segregating the convective, BB, and stratiform precipitation. This algorithm has been tested with radar data of different precipitation events and evaluated with three months of rain gauge data. The results show that the proposed algorithm performs consistently well for challenging precipitation events with the convection adjacent to or within a strong BB. Furthermore, the proposed algorithm could be used to improve the vertical profile of reflectivity (VPR) correction and reduce the overestimation of rainfall in the BB precipitation region.

With the launch of GOES-16 in November 2016, effective utilization of its data in convective-scale numerical weather prediction (NWP) has the potential to improve high-impact weather (HIWeather) forecasts. In this study, the impact of satellite-derived layered precipitable water (LPW) and cloud water path (CWP) in addition to NEXRAD observations on short-term convective-scale NWP forecasts are examined using three severe weather cases that occurred in May 2017. In each case, satellite-derived CWP and LPW products and radar observations are assimilated into the Advanced Research Weather Research and Forecasting (WRF-ARW) Model using the NSSL hybrid Warn-on-Forecast (WoF) analysis and forecast system. The system includes two components: the GSI-EnKF system and a deterministic 3DEnVAR system. This study examines deterministic 0–6-h forecasts launched from the hybrid 3DEnVAR analyses for the three severe weather events. Three types of experiments are conducted and compared: (i) the control experiment (CTRL) without assimilating any data, (ii) the radar experiment (RAD) with the assimilation of radar and surface observations, and (iii) the satellite experiment (RADSAT) with the assimilation of all observations including surface-, radar-, and satellite-derived CWP and LPW. The results show that assimilating additional GOES products improves short-range forecasts by providing more accurate initial conditions, especially for moisture and temperature variables.

Remote sensing techniques have been used to study and track wildfire smoke plume structure and evolution; however, knowledge gaps remain because of the limited availability of observational datasets aimed at understanding fine-scale fire–atmosphere interactions and plume microphysics. Meteorological radars have been used to investigate the evolution of plume rise in time and space, but highly resolved plume observations are limited. In this study, we present a new mobile millimeter-wave (Ka band) Doppler radar system acquired to sample the fine-scale kinematics and microphysical properties of active wildfire smoke plumes from both wildfires and large prescribed fires. Four field deployments were conducted in autumn of 2019 during two wildfires in California and one prescribed burn in Utah. Radar parameters investigated in this study include reflectivity, radial velocity, Doppler spectrum width, differential reflectivity ZDR, and copolarized correlation coefficient ρHV. Observed radar reflectivity ranged between −15 and 20 dBZ in plume, and radial velocity ranged from 0 to 16 m s−1. Dual-polarimetric observations revealed that scattering sources within wildfire plumes are primarily nonspherical and oblate-shaped targets as indicated by ZDR values measuring above 0 and ρHV values below 0.8 within the plume. Doppler spectrum width maxima were located near the updraft core region and were associated with radar reflectivity maxima.

The intrinsic uncertainty of radar-based retrievals in snow originates from a large diversity of snow growth habits, densities, and particle size distributions, all of which can make interpreting radar measurements of snow very challenging. The application of polarimetric radar for snow measurements can mitigate some of these issues. In this study, a novel polarimetric method for quantification of the extinction coefficient and visibility in snow, based on the joint use of radar reflectivity at horizontal polarization Z and specific differential phase KDP, is introduced. A large 2D-video-disdrometer snow dataset from central Oklahoma is used to derive a polarimetric bivariate power-law relation for the extinction coefficient, . The relation is derived for particle aspect ratios ranging from 0.5 to 0.8 and the width of the canting angle distribution ranging from 0° to 40°, values typical of aggregated snow, and validated via theoretical and analytical derivations/simulations. The multiplier of the relation is sensitive to variations in particles’ densities, the width of the canting angle distribution, and particles’ aspect ratios, whereas the relation’s exponents are practically invariant to changes in the latter two parameters. This novel approach is applied to polarimetric S-band WSR-88D data and verified against previous studies and in situ measurements of the extinction coefficient for four snow events in the eastern United States. The polarimetric radar estimates of the extinction coefficient exhibit smaller biases in comparison to previous studies concerning the ground measurements. The results indicate that there is good potential for reliable radar estimates of visibility from polarimetric weather radars, a parameter inversely proportional to the extinction coefficient.

When directly assimilating radar data within a variational framework using hydrometeor mixing ratios (q) as control variables (CVq), the gradient of the cost function becomes extremely large when background mixing ratio is close to zero. This significantly slows down minimization convergence and makes the assimilation of radial velocity and other observations ineffective because of the dominance of the reflectivity observation term in the cost function gradient. Using logarithmic hydrometeor mixing ratios as control variables (CV logq) can alleviate the problem but the high nonlinearity of logarithmic transformation can introduce spurious analysis increments into mixing ratios. In this study, power transform of hydrometeors is proposed to form new control variables (CVpq) where the nonlinearity of transformation can be adjusted by a tuning exponent or power parameter p. The performance of assimilating radar data using CVpq is compared with those using CVq and CV logq for the analyses and forecasts of five convective storm cases from the spring of 2017. Results show that CVpq with p = 0.4 (CVpq0.4) gives the best reflectivity forecasts in terms of root-mean-square error and equitable threat score. Furthermore, CVpq0.4 has faster convergence of cost function minimization than CVq and produces less spurious analysis increment than CV logq. Compared to CVq and CV logq, CVpq0.4 has better skills of 0–3-h composite reflectivity forecasts, and the updraft helicity tracks for the 16 May 2017 Texas and Oklahoma tornado outbreak case are more consistent with observations when using CVpq0.4.

Hurricane Harvey in 2017 generated one of the most catastrophic rainfall events in United States history. Numerous gauge observations in Texas exceeded 1200 mm, and the record accumulations resulted in 65 direct fatalities from rainfall-induced flooding. This was followed by Hurricane Florence in 2018, where multiple regions in North Carolina received over 750 mm of rainfall. The Multi-Radar Multi-Sensor (MRMS) system provides the unique perspective of applying fully automated seamless radar mosaics and locally gauge-corrected products for these two historical tropical cyclone rainfall events. This study investigates the performance of various MRMS quantitative precipitation estimation (QPE) products as it pertains to rare extreme tropical cyclone rainfall events. Various biases were identified in the radar-only approaches, which were mitigated in a new dual-polarimetric synthetic radar QPE approach. A local gauge correction of radar-derived QPE provided statistical improvements over the radar-only products but introduced consistent underestimation biases attributed to undercatch from tropical cyclone winds. This study then introduces a conceptual methodology to bulk correct for gauge wind undercatch across the numerous gauge networks ingested by the MRMS system. Adjusting the hourly gauge observations for wind undercatch resulted in increased storm-total accumulations for both tropical cyclones that better matched independent gauge observations, yet its application across large network collections highlighted the challenges of applying a singular wind undercatch correction scheme for significant wind events (e.g., tropical cyclones) while recognizing the need for increased metadata on gauge characteristics.

With the launch of the Global Precipitation Measurement Dual-Frequency Precipitation Radar (GPM-DPR) in 2014, renewed interest in retrievals of snowfall in the atmospheric column has occurred. The current operational GPM-DPR retrieval largely underestimates surface snowfall accumulation. Here, a neural network (NN) trained on data that are synthetically derived from state-of-the-art ice particle scattering models and measured in situ particle size distributions (PSDs) is used to retrieve two parameters of the PSD: liquid equivalent mass-weighted mean diameter and the liquid equivalent normalized intercept parameter . Evaluations against a test dataset showed statistically significantly improved ice water content (IWC) retrievals relative to a standard power-law approach and an estimate of the current GPM-DPR algorithm. Furthermore, estimated median percent errors (MPE) on the test dataset were −0.7%, +2.6%, and +1% for , , and IWC, respectively. An evaluation on three case studies with collocated radar observations and in situ microphysical data shows that the NN retrieval has MPE of −13%, +120%, and +10% for , , and IWC, respectively. The NN retrieval applied directly to GPM-DPR data provides improved snowfall retrievals relative to the default algorithm, removing the default algorithm’s ray-to-ray instabilities and recreating the high-resolution radar retrieval results to within 15% MPE. Future work should aim to improve the retrieval by including PSD data collected in more diverse conditions and rimed particles. Furthermore, different desired outputs such as the PSD shape parameter and snowfall rate could be included in future iterations.

The variational method formulated in Part I for analyzing vortex flow (VF), called VF-Var, is tested with simulated radar radial-velocity observations from idealized and pseudo-operational Doppler scans of analytically formulated benchmark vortices with spiral-band structures to resemble VFs in observed tornadic mesocyclones. The idealized Doppler scans are unidirectional in parallel along horizontal grid lines of a coarse-resolution grid, so they measure only the horizontal components of three-dimensional velocities in the analysis domain. The pseudo-operational Doppler scans mimic a scan mode used by operational WSR-88Ds for severe storms. Paired numerical experiments are designed and performed to test the two-step analysis versus single-step analysis formulated in VF-Var. Both analyses perform very well with dual-Doppler scans and reasonably well with single-Doppler scans. Errors in the analyzed velocities from single-Doppler scans are mainly in the unobserved velocity components and only in fractions of the benchmark velocities. When the vortex is upright or slanted in the direction perpendicular to idealized single-Doppler scans, the two-step analysis slightly outperforms the single-step analysis for idealized Doppler scans and pseudo-operational dual-Doppler scans. When the vortex becomes slanted in the direction largely along or against Doppler scans, both analyses become less (more) accurate in analyzing the horizontal (slantwise vertical) velocity, and the single-step analysis outperforms the two-step analysis especially for single-Doppler scans. By considering the projections of analyzed velocity on radar beams in the original Cartesian coordinates, useful insights are gained for understanding why and how the analysis accuracies are affected by vortex slanting.

A variational method is formulated with theoretical considerations for analyzing vortex flows in Doppler radar–scanned tornadic mesocyclones. The method has the following features. (i) The vortex center axis (estimated as a continuous function of time and height in the four-dimensional space) is used as the vertical coordinate, so the coordinate system used for the analysis is slantwise curvilinear and nonorthogonal in general. (ii) The vortex flow (VF), defined by the three-dimensional vector wind minus the horizontal moving velocity of vortex center axis, is expressed in terms of the covariant basis vectors (tangent to the coordinate curves), so its axisymmetric part can be properly defined in that slantwise-curvilinear coordinate system. (iii) To satisfy the mass continuity automatically, the axisymmetric part is expressed by the scalar fields of azimuthally averaged tangential velocity and cylindrical streamfunction and the remaining asymmetric part is expressed by the scalar fields of streamfunction and vertically integrated velocity potential. (iv) VF-dependent covariance functions are formulated for these scalar variables and then deconvoluted to construct the square root of background error covariance matrix analytically with the latter used to transform the control vector to precondition the cost function. (v) The deconvoluted covariance functions and their transformed control variables satisfy two required boundary conditions (i.e., zero vertical velocity at the lower rigid boundary and zero cross-axis velocity along the vortex center axis), so the analyzed VF satisfies not only the mass continuity but also the two boundary conditions automatically.

Several studies in the literature have developed approaches to diagnose hail storms from satellite-borne passive microwave imagery and build nearly global climatologies of hail. This paper uses spaceborne Ku-band radar measurements from the Global Precipitation Measurement (GPM) mission Dual-Frequency Precipitation Radar (DPR) to validate several passive microwave approaches. We assess the retrievals on the basis of how tightly they constrain the radar reflectivity at −20°C and how this measured radar reflectivity aloft varies geographically. The algorithm that combines minimum 19-GHz polarization corrected temperature (PCT) with a 37-GHz PCT depression normalized by tropopause height constrains the radar reflectivity most tightly and gives the least appearance of regional biases. A retrieval that is based on a 19-GHz PCT threshold of 261K also produces tightly clustered profiles of radar reflectivity, with little regional bias. An approach using regionally adjusted minimum 37-GHz PCT performs relatively well, but our results indicate it may overestimate hail in some subtropical and midlatitude regions. A threshold applied to the minimum 37-GHz PCT (≤230 K), without any scaling by region or probability of hail, overestimates hail in the tropics and underestimates beyond the tropics. For all retrieval approaches, storms identified as having hail tended to have radar reflectivity profiles that are consistent with general expectations for hailstorms (reflectivity > 50 dBZ below the 0°C level, and > 40 dBZ extending far above 0°C). Profiles from oceanic regions tended to have more rapidly decreasing reflectivity with height than profiles from other regions. Subtropical, high-latitude, and high-terrain land profiles had the slowest decreases of reflectivity with height.

The Dual-Frequency Precipitation Radar (DPR) on board the Global Precipitation Measurement (GPM) mission core satellite provides the new-generation global observation of rain since 2014. The main objective of this paper is to evaluate the suitability and limitation of GPM-DPR level-2 products over China. The DPR rain rate products are compared with rain gauge data during the summers of 5 years (2014–18). The ground observation network is composed of more than 50 000 rain gauges. The DPR precipitation products for all scans (DPR_NS, DPR_MS, and DPR_HS) generally underestimate rain rates. However, DPR_MS agrees better with gauge estimates than DPR_NS and DPR_HS, yielding the lowest mean error, systematic deviation, and highest Pearson correlation coefficient. In addition, all three swath types show obvious overestimation over gauge estimates between 0.5 and 1 mm h−1 and underestimation when gauge estimates are larger than 1 mm h−1. The DPR_HS and DPR_MS agree better with gauge estimates below and above 2.5 mm h−1, respectively. A deeper investigation was carried out to analyze the variation of DPR_MS’s performance with respect to terrains over China. An obvious underestimation, relative to gauge estimates, occurs in Tibetan Plateau while a slight overestimation occurs in the North China Plain. Furthermore, our comprehensive analysis suggests that in Sichuan Basin, the DPR_MS exhibit the best agreement with gauge estimates.

Over the last two decades, the Central Weather Bureau of Taiwan and the U.S. National Severe Storms Laboratory have been involved in a research and development collaboration to improve the monitoring and prediction of river flooding, flash floods, debris flows, and severe storms for Taiwan. The collaboration resulted in the Quantitative Precipitation Estimation and Segregation Using Multiple Sensors (QPESUMS) system. The QPESUMS system integrates observations from multiple mixed-band weather radars, rain gauges, and numerical weather prediction model fields to produce high-resolution (1 km) and rapid-update (10 min) rainfall and severe storm monitoring and prediction products. The rainfall products are widely used by government agencies and emergency managers in Taiwan for flood and mudslide warnings as well as for water resource management. The 3D reflectivity mosaic and QPE products are also used in high-resolution radar data assimilation and for the verification of numerical weather prediction model forecasts. The system facilitated collaborations with academic communities for research and development of radar applications, including quantitative precipitation estimation and nowcasting. This paper provides an overview of the operational QPE capabilities in the Taiwan QPESUMS system.

The Radar Applications in Northern Scotland (RAiNS) experiment took place from February to August 2016 near Inverness, Scotland. The campaign was motivated by the need to provide enhanced weather radar observations for hydrological applications for the Inverness region. Here we describe the campaign in detail and observations over the summer period of the campaign that show the improvements that high-resolution polarimetric radar observations may have on quantitative precipitation estimates in this region compared to concurrently generated operational radar quantitative precipitation estimates (QPEs). We further provide suggestions of methods for generating QPE using dual-polarization X-band radars in similar regions.

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. The combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

Accurate, physically based precipitation retrieval over global land surfaces is an important goal of the NASA/JAXA Global Precipitation Measurement Mission (GPM). This is a difficult problem for the passive microwave constellation, as the signal over radiometrically warm land surfaces in the microwave frequencies means that the measurements used are indirect and typically require inferring some type of relationship between an observed scattering signal and precipitation at the surface. GPM, with collocated radiometer and dual-frequency radar, is an excellent tool for tackling this problem and improving global retrievals. In the years following the launch of the GPM Core Observatory satellite, physically based passive microwave retrieval of precipitation over land continues to be challenging. Validation efforts suggest that the operational GPM passive microwave algorithm, the Goddard profiling algorithm (GPROF), tends to overestimate precipitation at the low (−1) end of the distribution over land. In this work, retrieval sensitivities to dynamic surface conditions are explored through enhancement of the algorithm with dynamic, retrieved information from a GPM-derived optimal estimation scheme. The retrieved parameters describing surface and background characteristics replace current static or ancillary GPROF information including emissivity, water vapor, and snow cover. Results show that adding this information decreases probability of false detection by 50% and, most importantly, the enhancements with retrieved parameters move the retrieval away from dependence on ancillary datasets and lead to improved physical consistency.

The objective of this research is to develop techniques for assimilating GOES-R series observations in precipitating scenes for the purpose of improving short-term convective-scale forecasts of high-impact weather hazards. Whereas one approach is radiance assimilation, the information content of GOES-R radiances from its Advanced Baseline Imager saturates in precipitating scenes, and radiance assimilation does not make use of lightning observations from the GOES Lightning Mapper. Here, a convolutional neural network (CNN) is developed to transform GOES-R radiances and lightning into synthetic radar reflectivity fields to make use of existing radar assimilation techniques. We find that the ability of CNNs to utilize spatial context is essential for this application and offers breakthrough improvement in skill compared to traditional pixel-by-pixel based approaches. To understand the improved performance, we use a novel analysis method that combines several techniques, each providing different insights into the network’s reasoning. Channel-withholding experiments and spatial information–withholding experiments are used to show that the CNN achieves skill at high reflectivity values from the information content in radiance gradients and the presence of lightning. The attribution method, layerwise relevance propagation, demonstrates that the CNN uses radiance and lightning information synergistically, where lightning helps the CNN focus on which neighboring locations are most important. Synthetic inputs are used to quantify the sensitivity to radiance gradients, showing that sharper gradients produce a stronger response in predicted reflectivity. Lightning observations are found to be uniquely valuable for their ability to pinpoint locations of strong radar echoes.

In this study, a hybrid En3DVar data assimilation (DA) scheme is compared with 3DVar, EnKF, and pure En3DVar for the assimilation of radar data in a real tornadic storm case. Results using hydrometeor mixing ratios (CVq) or logarithmic mixing ratios (CVlogq) as the control variables are compared in the variational DA framework. To address the lack of radial velocity impact issues when using CVq, a procedure that assimilates reflectivity and radial velocity data in two separate analysis passes is adopted. Comparisons are made in terms of the root-mean-square innovations (RMSIs) as well as the intensity and structure of the analyzed and forecast storms. For pure En3DVar that uses 100% ensemble covariance, CVlogq and CVq have similar RMSIs in the velocity analyses, but errors grow faster during forecasts when using CVlogq. Introducing static background error covariance at 5% in hybrid En3DVar (with CVlogq) significantly reduces the forecast error growth. Pure En3DVar produces more intense reflectivity analyses than EnKF that more closely match the observations. Hybrid En3DVar with 50% outperforms other weights in terms of the RMSIs and forecasts of updraft helicity and is thus used in the final comparison with 3DVar and EnKF. The hybrid En3DVar is found to outperform EnKF in better capturing the intensity and structure of the analyzed and forecast storms and outperform 3DVAR in better capturing the intensity and evolution of the rotating updraft.

Partial beam blockage (PBB) correction is an indispensable step in weather radar data quality control and subsequent quantitative applications, such as precipitation estimation, especially in urban and/or complex terrain environments. This paper developed a novel PBB correction procedure based on the improved ZPHI method for attenuation correction and regional specific differential propagation phase (KDP)–reflectivity (ZH) relationship derived from in situ raindrop size distribution (DSD) measurements. The practical performance of this PBB correction technique was evaluated through comparing the spatial continuity of reflectivity measurements, the consistency between radar-measured and DSD-derived KDP and ZH relationships, as well as rainfall estimates based on R(ZH) and R(KDP). The results showed that through incorporating attenuation and PBB corrections (i) the spatial continuity of ZH measurements can effectively be enhanced; (ii) the distribution of radar-measured KDP versus ZH is more consistent with the DSD-derived KDP versus ZH; (iii) the measured ZH from a C-band radar in the PBB-affected area becomes more consistent with collocated S-band measurements, particularly in the rainstorm center area where ZH is larger than 30 dBZ; and (iv) rainfall estimates based on R(ZH) in the PBB-affected area are incrementally improved with better spatial continuity and the performance tends to be more comparable with R(KDP).

Dual-frequency millimeter-wavelength radar observations in snowfall are analyzed in order to evaluate differences in conventional polarimetric radar variables such as differential reflectivity (ZDR) specific differential phase shift (KDP) and linear depolarization ratio (LDR) at traditional cloud radar frequencies at Ka and W bands (~35 and ~94 GHz, correspondingly). Low radar beam elevation (~5°) measurements were performed at Oliktok Point, Alaska, with a scanning fully polarimetric radar operating in the horizontal–vertical polarization basis. This radar has the same gate spacing and very close beam widths at both frequencies, which largely alleviates uncertainties associated with spatial and temporal data matching. It is shown that observed Ka- and W-band ZDR differences are, on average, less than about 0.5 dB and do not have a pronounced trend as a function of snowfall reflectivity. The observed ZDR differences agree well with modeling results obtained using integration over nonspherical ice particle size distributions. For higher signal-to-noise ratios, KDP data derived from differential phase measurements are approximately scaled as reciprocals of corresponding radar frequencies indicating that the influence of non-Rayleigh scattering effects on this variable is rather limited. This result is also in satisfactory agreement with data obtained by modeling using realistic particle size distributions. Observed Ka- and W-band LDR differences are strongly affected by the radar hardware system polarization “leak” and are generally less than 4 dB. Smaller differences are observed for higher depolarizations, where the polarization “leak” is less pronounced. Realistic assumptions about particle canting and the system polarization isolation lead to modeling results that satisfactorily agree with observational dual-frequency LDR data.

This paper examines nowcasts of precipitation from the High-Resolution Rapid Refresh (HRRRv2) model from the summer of 2017 along the Colorado Front Range. It was found that model nowcasts (2 h or less) of precipitation amount were less skillful than extrapolation of the KFTG WSR-88-D data at a spatial scale of 120 km. It was also found that local-scale (mesoscale) influences on rainfall intensity and amount have a much greater impact on rainfall intensity than large-scale (synoptic) influences. Thus, large-scale trends are not useful for modifying extrapolation nowcasts on the local scale. Errors in the HRRR nowcasts are attributed to an inability of the model and data assimilation to resolve convergence along outflow boundaries and other terrain-influenced mesogamma-scale flows that contribute to storm formation and evolution. While the HRRRv2 1-h nowcasts were strongly correlated with observed precipitation events, the nowcast precipitation amounts were in error by more than a factor of 2 about 50% of the time, with half of the cases being overestimates and half being underestimates. A large fraction of the HRRRv2 overestimates were associated with stratiform rain events. It is speculated that this was a result of misinterpretation of the radar bright band as more intense precipitation aloft by the data assimilation scheme. A large fraction of the HRRRv2 underestimates occurred when the data assimilation and model were unable to fully resolve the low-level convergence along small-scale, narrow boundaries that led to new storm initiation and/or storm growth.

The method of neural networks (aka deep learning) has opened up many new opportunities to utilize remotely sensed images in meteorology. Common applications include image classification, e.g., to determine whether an image contains a tropical cyclone, and image-to-image translation, e.g., to emulate radar imagery for satellites that only have passive channels. However, there are yet many open questions regarding the use of neural networks for working with meteorological images, such as best practices for evaluation, tuning, and interpretation. This article highlights several strategies and practical considerations for neural network development that have not yet received much attention in the meteorological community, such as the concept of receptive fields, underutilized meteorological performance measures, and methods for neural network interpretation, such as synthetic experiments and layer-wise relevance propagation. We also consider the process of neural network interpretation as a whole, recognizing it as an iterative meteorologist-driven discovery process that builds on experimental design and hypothesis generation and testing. Finally, while most work on neural network interpretation in meteorology has so far focused on networks for image classification tasks, we expand the focus to also include networks for image-to-image translation.

Tor Bergeron was a key member of the Bergen School of Meteorology that developed some of the most influential contributions to synoptic analysis in the twentieth century: airmass analysis, polar-front theory, and the Norwegian cyclone model. However, the eventual success of these so-called Bergen methods of synoptic analysis was not guaranteed. Concerns and criticisms of the methods—in part from the lack of referencing to prior studies, overly simplified conceptual models, and lack of real data in papers by J. Bjerknes and Solberg—were inhibiting worldwide adoption. Bergeron’s research output in the 1920s was aimed at addressing these concerns. His doctoral thesis, written in German, was published as a journal article in Geofysiske Publikasjoner in 1928. Here, an accessible and annotated English translation is provided along with a succinct overview of this seminal study. Major interlaced themes of Bergeron’s study were the first comprehensive description of the Bergen methods: a vigorous defense of cyclogenesis as primarily a lower-tropospheric process as opposed to an upper-tropospheric–lower-stratospheric one; a nuanced explanation of the assertion that meteorology constituted a distinct and special scientific discipline; and, very understandably, a thorough account of Bergeron’s own contributions to the Bergen School. His contributions included identifying how deformation results in frontogenesis and frontolysis, classifying the influence of aerosols on visibility, and explaining the role of the ambient conditions in the onset of drizzle as opposed to rain showers—a distinction that led the formulation of the Wegener–Bergeron–Findeisen process.

Recently, the Central Weather Bureau of Taiwan developed a WRF- and WRF data assimilation (WRFDA)-based convective-scale data assimilation system to increase model predictability toward high-impact weather. In this study, we focus on afternoon thunderstorm (AT) prediction and investigate the following questions: 1) Is the designation of a rapid update cycle strategy with a blending scheme effective? 2) Can surface data assimilation contribute positively to AT prediction under the complex geography of Taiwan island? 3) What is the relative importance between radar and surface observation to AT prediction? 4) Can we increase the AT forecast lead time in the morning through data assimilation? Consecutive ATs from 30 June to 8 July 2017 are investigated. Five experiments, each having 240 continuous cycles, are designed. Results show that employing continuous cycles with a blending scheme mitigates model spinup compared with downscaled forecasts. Although there are few radar echoes before AT initiation, assimilating radar observations is still crucial since it largely corrects model errors in cycles. However, assimilating surface observations is more important compared with radar in terms of extending forecast lead time in the morning. Either radar or surface observations contribute positively, and assimilating both has the highest QPF score. Assimilating surface observations systematically improves surface wind and temperature predictions based on 240 cases. A case study demonstrates that the model can capture the AT initiation and development by assimilating surface and radar observations. Its cold pool and outflow boundary prediction are also improved. In this case, the assimilation of surface wind and water vapor in the morning contributes more compared with temperature and pressure.

This research explores the benefits of radar data assimilation for short-range weather forecasts in southeastern Brazil using the Weather Research and Forecasting (WRF) Model’s three-dimensional variational data assimilation (3DVAR) system. Different data assimilation options are explored, including the cycling frequency, the number of outer loops, and the use of null-echo assimilation. Initially, four microphysics parameterizations are evaluated (Thompson, Morrison, WSM6, and WDM6). The Thompson parameterization produces the best results, while the other parameterizations generally overestimate the precipitation forecast, especially WDSM6. Additionally, the Thompson scheme tends to overestimate snow, while the Morrison scheme overestimates graupel. Regarding the data assimilation options, the results deteriorate and more spurious convection occurs when using a higher cycling frequency (i.e., 30 min instead of 60 min). The use of two outer loops produces worse precipitation forecasts than the use of one outer loop, and the null-echo assimilation is shown to be an effective way to suppress spurious convection. However, in some cases, the null-echo assimilation also removes convective clouds that are not observed by the radar and/or are still not producing rain, but have the potential to grow into an intense convective cloud with heavy rainfall. Finally, a cloud convective mask was implemented using ancillary satellite data to prevent null-echo assimilation from removing potential convective clouds. The mask was demonstrated to be beneficial in some circumstances, but it needs to be carefully evaluated in more cases to have a more robust conclusion regarding its use.

形態: カラー図版あり, Physical characteristics: Original contains color illustrations, 資料番号: PA2110044000

A new dual-polarization (DP) radar synthetic quantitative precipitation estimation (QPE) product was developed using a combination of specific attenuation A, specific differential phase KDP, and reflectivity Z to calculate the precipitation rate R. Specific attenuation has advantages of being insensitive to systematic biases in Z and differential reflectivity ZDR due to partial beam blockage, attenuation, and calibration while more linearly related to R than other radar variables. However, the R(A) relationship is not applicable in areas containing ice. Therefore, the new DP QPE applies R(A) in areas where radar is observing pure rain, R(KDP) in regions potentially containing hail, and R(Z) elsewhere. Further, an evaporation correction was applied to minimize false light precipitation related to virga. The new DP QPE was evaluated in real time over the conterminous United States and showed significant improvements over previous radar QPE techniques that were based solely on R(Z) relationships. The improvements included reduced dry biases in heavy to extreme precipitation during the warm season. The new DP QPE also reduced errors and spatial discontinuities in regions impacted by partial beam blockage. Further, the new DP QPE reduced wet bias for scattered light precipitation in the southwest and north central United States where there is significant boundary layer evaporation.

This paper describes a technique for estimating the liquid water content (LWC) and a characteristic particle diameter in stratocumulus clouds using radar and lidar observations. The uncertainty in LWC estimate from radar and lidar measurements is significantly reduced once the characteristic particle diameter is known. The technique is independent of the drop size distribution. It is applicable for a broad range of W-band reflectivityZbetween −30 and 0 dBZand all values of lidar backscatterβobservations. No partitioning of cloud or drizzle is required on the basis of an arbitrary threshold ofZas in prior studies. A method for estimating droplet diameter and LWC was derived from the electromagnetic simulations of radar and lidar observations. In situ stratocumulus cloud and drizzle probe spectra were input to the electromagnetic simulation. The retrieved droplet diameter and LWC were validated using in situ measurements from the southeastern Pacific Ocean. The retrieval method was applied to radar and lidar measurements from the northeastern Pacific. Uncertainty in the retrieved droplet diameter and LWC that are due to the measurement errors in radar and lidar backscatter measurements are 7% and 14%, respectively. The retrieved LWC was validated using the concurrent G-band radiometer estimates of the liquid water path.

Tornadic supercells moved across parts of Oklahoma on the afternoon and evening of 9 May 2016. One such supercell, while producing a long-lived tornado, was observed by nearby WSR-88D radars to contain a strong anticyclonic velocity couplet on the lowest elevation angle. This couplet was located in a very atypical position relative to the ongoing cyclonic tornado and to the supercell’s updraft. A storm survey team identified damage near where this couplet occurred, and, in the absence of evidence refuting otherwise, the damage was thought to have been produced by an anticyclonic tornado. However, such a tornado was not seen in near-ground, high-resolution radar data from a much closer, rapid-scan, mobile radar. Rather, an elongated velocity couplet was observed only at higher elevation angles at altitudes similar to those at which the WSR-88D radars observed the strong couplet. This paper examines observations from two WSR-88D radars and a mobile radar from which it is argued that the anticyclonic couplet (and a similar one ~10 min later) were actually quasi-horizontal vortices centered ~1–1.5 km AGL. The benefits of having data from a radar much closer to the convective storm being sampled (e.g., better spatial resolution and near-ground data coverage) and providing more rapid volume updates are readily apparent. An analysis of these additional radar data provides strong, but not irrefutable, evidence that the anticyclonic tornado that may be inferred from WSR-88D data did not exist; consequently, upon discussions with the National Weather Service, it was not included in Storm Data.

A novel 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.

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.

The Multi-Radar-Multi-Sensor (MRMS) system was transitioned into operations at the National Centers for Environmental Prediction in the fall of 2014. It provides high-quality and high-resolution severe weather and precipitation products for meteorology, hydrology, and aviation applications. Among processing modules, the radar data quality control (QC) plays a critical role in effectively identifying and removing various nonhydrometeor radar echoes for accurate quantitative precipitation estimation (QPE). Since its initial implementation in 2014, the radar QC has undergone continuous refinements and enhancements to ensure its robust performance across seasons and all regions in the continental United States and southern Canada. These updates include 1) improved melting-layer delineation, 2) clearance of wind farm contamination, 3) mitigation of corrupt data impacts due to hardware issues, 4) mitigation of sun spikes, and 5) mitigation of residual ground/lake/sea clutter due to sidelobe effects and anomalous propagation. This paper provides an overview of the MRMS radar data QC enhancements since 2014.

A polarimetric radar–based method for retrieving atmospheric ice particle shapes is applied to snowfall measurements by a scanning Ka-band radar deployed at Oliktok Point, Alaska (70.495°N, 149.883°W). The mean aspect ratio, which is defined by the hydrometeor minor-to-major dimension ratio for a spheroidal particle model, is retrieved as a particle shape parameter. The radar variables used for aspect ratio profile retrievals include reflectivity, differential reflectivity, and the copolar correlation coefficient. The retrievals indicate that hydrometeors with mean aspect ratios below 0.2–0.3 are usually present in regions with air temperatures warmer than approximately from −17° to −15°C, corresponding to a regime that has been shown to be favorable for growth of pristine ice crystals of planar habits. Radar reflectivities corresponding to the lowest mean aspect ratios are generally between −10 and 10 dBZ. For colder temperatures, mean aspect ratios are typically in a range between 0.3 and 0.8. There is a tendency for hydrometeor aspect ratios to increase as particles transition from altitudes in the temperature range from −17° to −15°C toward the ground. This increase is believed to result from aggregation and riming processes that cause particles to become more spherical and is associated with areas demonstrating differential reflectivity decreases with increasing reflectivity. Aspect ratio retrievals at the lowest altitudes are consistent with in situ measurements obtained using a surface-based multiangle snowflake camera. Pronounced gradients in particle aspect ratio profiles are observed at altitudes at which there is a change in the dominant hydrometeor species, as inferred by spectral measurements from a vertically pointing Doppler radar.

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.

In this study, a 10-yr (2008–17) radar-based mesoscale convective system (MCS) and derecho climatology for Poland is presented. This is one of the first attempts of a European country to investigate morphological and precipitation archetypes of MCSs as prior studies were mostly based on satellite data. Despite its ubiquity and significance for society, economy, agriculture, and water availability, little is known about the climatological aspects of MCSs over central Europe. Our results indicate that MCSs are not rare in Poland as an annual mean of 77 MCSs and 49 days with MCS can be depicted for Poland. Their lifetime ranges typically from 3 to 6 h, with initiation time around the afternoon hours (1200–1400 UTC) and dissipation stage in the evening (1900–2000 UTC). The most frequent morphological type of MCSs is a broken line (58% of cases), then areal/cluster (25%), and then quasi-linear convective systems (QLCS; 17%), which are usually associated with a bow echo (72% of QLCS). QLCS are the feature with the longest life cycle. Among precipitation archetypes of linear MCSs, trailing stratiform (73%) and parallel stratiform (25%) are the most common. MCSs are usually observed from April to September, with a peak in mid-July. A majority of MCSs travels from the west, southwest, and south sectors. A total of 16 derecho events were identified (1.5% of all MCS and 9.1% of all QLCS); the majority of them were produced by a warm-season QLCS, whereas only 4 were produced by cold-season narrow cold-frontal rainbands. Warm-season derechos produced a bigger impact than did cold-season events, even though their damage paths were shorter.

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.

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 dual-polarized 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.

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.

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.

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. The MSQPE 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.

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 α. In this study, a variational approach is proposed to optimize the coefficient α 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 α consistent with the DSD properties of precipitation on a sweep-to-sweep basis. The attenuation-based hourly rainfall estimates using the sweep-averaged α 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 α for each radar ray, which mitigates the impact of the azimuthal DSD variabilities on rainfall estimation. It is found that, based on the optimized α 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 sweep-averaged α.

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.

This paper describes the development of convolutional neural networks (CNN), a type of deep-learning method, to predict next-hour tornado occurrence. Predictors are a storm-centered radar image and a proximity sounding from the Rapid Refresh model. Radar images come from the Multiyear Reanalysis of Remotely Sensed Storms (MYRORSS) and Gridded NEXRAD WSR-88D Radar dataset (GridRad), both of which are multiradar composites. We train separate CNNs on MYRORSS and GridRad data, present an experiment to optimize the CNN settings, and evaluate the chosen CNNs on independent testing data. Both models achieve an area under the receiver-operating-characteristic curve (AUC) well above 0.9, which is considered to be excellent performance. The GridRad model achieves a critical success index (CSI) of 0.31, and the MYRORSS model achieves a CSI of 0.17. The difference is due primarily to event frequency (percentage of storms that are tornadic in the next hour), which is 3.52% for GridRad but only 0.24% for MYRORSS. The best CNN predictions (true positives and negatives) occur for strongly rotating tornadic supercells and weak nontornadic cells in mesoscale convective systems, respectively. The worst predictions (false positives and negatives) occur for strongly rotating nontornadic supercells and tornadic cells in quasi-linear convective systems, respectively. The performance of our CNNs is comparable to an operational machine-learning system for severe weather prediction, which suggests that they would be useful for real-time forecasting.

Differential absorption radar (DAR) offers an active remote sensing solution to the problem of measuring humidity profiles with high vertical and horizontal resolution in hydrometeor layers. The Vapor In-Cloud Profiling Radar (VIPR) is a frequency-modulated continuous-wave (FMCW) G-band DAR tunable from 167 to 174.8 GHz being developed at the Jet Propulsion Laboratory (JPL). Here we describe ground-based measurements from VIPR performed at the Department of Energy’s Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site for humidity product validation. Two distinct measurement capabilities are investigated: 1) humidity profiles inside of cloudy volumes with 180 m vertical resolution, and 2) integrated water vapor (IWV) between the surface and cloud base. High radar sensitivity permits detection of upper-tropospheric clouds and retrieval of humidity profiles above 10 km in height. We develop an improved humidity retrieval algorithm based on a regularized least squares method that includes detailed accounting of measurement covariances and systematic error sources. This regularization mitigates high-spatial-frequency humidity biases that arise from frequency-dependent hydrometeor scattering, which is an important limitation for DAR systems. Through comparisons with over 20 coincident radiosondes, we find close agreement between in situ and remotely sensed humidity profiles, with a correlation coefficient of r = 0.96, root-mean-square error (RMSE) of 0.8 g m−3, and median retrieval precision of 0.5 g m−3. Using a merged radiosonde and Raman lidar product for surface-to-cloud-base IWV, we demonstrate precise column sounding capabilities with r = 1.00, RMSE of 1.2 mm, and median retrieval precision of 0.25 mm.

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 contentMand mean diameterDmseem robust to the change in meteorological regime (between GCPEx and OLYMPEx); whereas the retrieval of the diameter spreadSmseems more campaign dependent.

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.

This study demonstrates an implementation of the prototype quantitative precipitation R estimation algorithm using specific attenuation A for S-band polarimetric radar. The performance of R(A) algorithm is assessed, compared to the conventional algorithm using radar reflectivity Z, at multiple temporal scales. Because the factor α, defined as the net ratio of A to specific differential phase, is a key parameter of the algorithm characterized by drop size distributions (e.g., differential reflectivity Zdr dependence on Z), the estimation equations of α and a proper number of Zdr–Z samples required for a reliable α estimation are examined. Based on the dynamic estimation of α, the event-based evaluation using hourly rain gauge observations reveals that the performance of R(A) is superior to that of R(Z), with better agreement and lower variability. Despite its superiority, the study finds that R(A) leads to quite consistent overestimations of about 10%–30%. It is demonstrated that the application of uniform α over the entire radar domain yields the observed uncertainty because of the heterogeneity of precipitation in the domain. A climatological range-dependent feature of R(A) and R(Z) is inspected in the multiyear evaluation at yearly scale using rain totals for April–October. While R(Z) exposes a systematic shift and overestimation, each of which arise from the radar miscalibration and bright band effects, R(A) combining with multiple R(Z) values for solid/mixed precipitation shows relatively robust performance without those effects. The immunity of R(A) to partial beam blockage (PBB) based on both qualitative and quantitative analyses is also verified. However, the capability of R(A) regarding PBB is limited by the presence of the melting layer and its application requirement for the total span of differential phase (e.g., 3°), which is another challenge for light rain.