Search

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.