Search for: All records

The assimilation of radar reflectivity requires an accurate and efficient forward operator that links the model state variables to radar observations. In this study, newly developed parameterized forward operators (PFO) for radar reflectivity with a new continuous melting model are implemented to assimilate observed radar data. To assess the impact of the novel parameterized reflectivity forward operators on convective storm analysis and forecasting, two distinct sets of cycled assimilation and forecast experiments are conducted. One set of experiments (ExpRFO) uses a conventional Rayleigh‐scattering‐approximation‐based forward operator (RFO) with hydrometeor classification, while the other uses the PFO (ExpPFO_New) for radar reflectivity with a new continuous melting model. Eight high‐impact severe convective weather events from the Hazardous Weather Testbed (HWT) 2019 Spring Experiments are selected for this study. The analysis and forecast results are first examined in detail for a classic tornadic supercell case on 24 May 2019, with the potential benefits provided by the PFO then evaluated for all eight cases. It is demonstrated that ExpPFO_New provides more robust results in terms of improving the short‐term severe weather forecasts. Compared to ExpRFO, ExpPFO_New better reproduces all observed supercells in the analysis field, yields a more continuous and reasonable reflectivity distribution near the melting layer, and improves the strength of the cold pool compared to observations. Overall, ExpPFO_New, initialized from the more accurate analysis fields, produces better forecasts of reflectivity and hourly precipitation with smaller biases, especially at heavy precipitation thresholds.

 

To improve short‐term severe weather forecasts through assimilation of polarimetric radar data (PRD), the use of accurate and efficient forward operators for polarimetric radar variables is required. In this study, a new melting model is proposed to estimate the mixing ratio and number concentration of melting hydrometeor species and incorporated in a set of parameterized polarimetric radar forward operators. The new melting model depends only on the mixing ratio and number concentration of rain and ice species and is characterized by its independence from ambient temperature and its simplicity and ease of linearization. To assess the impact of this newly proposed melting model on the simulated polarimetric radar variables, a real mesoscale convective system is simulated using three double‐moment microphysics schemes. Compared with the output of the original implementation of the parameterized forward operators (PFO_Old) that rely on an “old” melting model which only estimates the mixing ratio of the melting species, the updated implementation with the new melting model (PFO_New) that estimates both the mixing ratio and number concentration of melting species eliminates the very large mass/volume‐weighted mean diameter (D_m) at the bottom of the melting layer and produces more reasonable melting layer signatures for all three double‐moment microphysics schemes that more closely match the corresponding radar observations. This suggests that the new melting model has more reasonable implicit estimates of mixing ratios and number concentrations of melting hydrometeor species than the “old” melting model.

 

Abstract Quasi-vertical profiles (QVPs) of polarimetric radar data have emerged as a powerful tool for studying precipitation microphysics. Various studies have found enhancements in specific differential phase K dp in regions of suspected secondary ice production (SIP) due to rime splintering. Similar K dp enhancements have also been found in regions of sublimating snow, another proposed SIP process. This work explores these K dp signatures for two cases of sublimating snow using nearly collocated S- and Ka-band radars. The presence of the signature was inconsistent between the radars, prompting exploration of alternative causes. Idealized simulations are performed using a radar beam-broadening model to explore the impact of nonuniform beam filling (NBF) on the observed reflectivity Z and K dp within the sublimation layer. Rather than an intrinsic increase in ice concentration, the observed K dp enhancements can instead be explained by NBF in the presence of sharp vertical gradients of Z and K dp within the sublimation zone, which results in a K dp bias dipole. The severity of the bias is sensitive to the Z gradient and radar beamwidth and elevation angle, which explains its appearance at only one radar. In addition, differences in scanning strategies and range thresholds during QVP processing can constructively enhance these positive K dp biases by excluding the negative portion of the dipole. These results highlight the need to consider NBF effects in regions not traditionally considered (e.g., in pure snow) due to the increased K dp fidelity afforded by QVPs and the subsequent ramifications this has on the observability of sublimational SIP. Significance Statement Many different processes can cause snowflakes to break apart into numerous tiny pieces, including when they evaporate into dry air. Purported evidence of this phenomenon has been seen in data from some weather radars, but we noticed it was not seen in data from others. In this work we use case studies and models to show that this signature may actually be an artifact from the radar beam becoming too big and there being too much variability of the precipitation within it. While this breakup process may actually be occurring in reality, these results suggest we may have trouble observing it with typical weather radars. 

Abstract This study evaluates ice particle size distribution and aspect ratio φ Multi-Radar Multi-Sensor (MRMS) dual-polarization radar retrievals through a direct comparison with two legs of observational aircraft data obtained during a winter storm case from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. In situ cloud probes, satellite, and MRMS observations illustrate that the often-observed K dp and Z DR enhancement regions in the dendritic growth layer can either indicate a local number concentration increase of dry ice particles or the presence of ice particles mixed with a significant number of supercooled liquid droplets. Relative to in situ measurements, MRMS retrievals on average underestimated mean volume diameters by 50% and overestimated number concentrations by over 100%. IWC retrievals using Z DR and K dp within the dendritic growth layer were minimally biased relative to in situ calculations where retrievals yielded −2% median relative error for the entire aircraft leg. Incorporating φ retrievals decreased both the magnitude and spread of polarimetric retrievals below the dendritic growth layer. While φ radar retrievals suggest that observed dendritic growth layer particles were nonspherical (0.1 ≤ φ ≤ 0.2), in situ projected aspect ratios, idealized numerical simulations, and habit classifications from cloud probe images suggest that the population mean φ was generally much higher. Coordinated aircraft radar reflectivity with in situ observations suggests that the MRMS systematically underestimated reflectivity and could not resolve local peaks in mean volume diameter sizes. These results highlight the need to consider particle assumptions and radar limitations when performing retrievals. significance statement Developing snow is often detectable using weather radars. Meteorologists combine these radar measurements with mathematical equations to study how snow forms in order to determine how much snow will fall. This study evaluates current methods for estimating the total number and mass, sizes, and shapes of snowflakes from radar using images of individual snowflakes taken during two aircraft legs. Radar estimates of snowflake properties were most consistent with aircraft data inside regions with prominent radar signatures. However, radar estimates of snowflake shapes were not consistent with observed shapes estimated from the snowflake images. Although additional research is needed, these results bolster understanding of snow-growth physics and uncertainties between radar measurements and snow production that can improve future snowfall forecasting. 

The differential reflectivity (Z_DR) column is a notable polarimetric signature related to updrafts in deep moist convection. In this study, pseudo–water vapor (q_υ) observations are retrieved from observedZ_DRcolumns under the assumption that humidity is saturated within the convection whereZ_DRcolumns are detected, and are then assimilated within the 3DVar framework. The impacts of assimilating pseudo-q_υobservations fromZ_DRcolumns on short-term severe weather prediction are first evaluated for a squall-line case. Radar data analysis indicates that theZ_DRcolumns are mainly located on the inflow side of the high-reflectivity region. Assimilation of the pseudo-q_υobservations leads to an enhancement ofq_υwithin the convection, while concurrently reducing humidity in no-rain areas. Sensitivity experiments indicate that a tuned smaller observation error and a shorter horizontal decorrelation scale are optimal for a better assimilation of pseudo-q_υfromZ_DRcolumns, resulting in more stable rain rates during short-term forecasts. Additionally, a 15-min cycling assimilation frequency yields the best performance, providing the most accurate reflectivity forecast in terms of both location and intensity. Analysis of thermodynamic fields reveal that assimilatingZ_DRcolumns provides more favorable initial conditions for sustaining convection, including sustainable moisture condition, a strong cold pool, and divergent winds near the surface, consequently enhancing reflectivity and precipitation. With the optimal configuration determined from the sensitivity tests, a quantitative evaluation further demonstrates that assimilating the pseudo-q_υobservations fromZ_DRcolumns using the 3DVar method can improve the 0–3-h reflectivity and accumulated precipitation predictions of convective storms.

 

Abstract Snow sublimating in dry air is a forecasting challenge and can delay the onset of surface snowfall and affect storm-total accumulations. Despite this, it remains comparatively less studied than other microphysical processes. Herein, the characteristics of sublimating snow and the potential for nowcasting snowfall reaching the surface are explored through the use of dual-polarization radar. Twelve cases featuring prolific sublimation were analyzed using range-defined quasi-vertical profiles (RDQVPs) and compared with environmental model analyses. Overall, reflectivity Z significantly decreases, differential reflectivity Z DR slightly decreases, and copolar-correlation coefficient ρ hv remains nearly constant through the sublimation layer. Regions of enhanced specific differential phase K dp were frequently observed in the sublimation layer and are believed to be polarimetric evidence of secondary ice production via sublimation. A 1D bin model was initialized using particle size distributions retrieved from the RDQVPs using numerous novel polarimetric snowretrieval relations for a wide range of forecast lead times, with the model environment evolving in response to sublimation. It was found that the model was largely able to predict the snowfall start time up to six hours in advance, with a 6-h median bias of just -18.5 minutes. A more detailed case study of the 08 December 2013 snowstorm in the Philadelphia region was also performed, demonstrating good correspondence with observations and examples of model fields (e.g., cooling rate) hypothetically available from such a tool. The proof-of-concept results herein demonstrate the potential benefits of incorporating spatially averaged radar data in conjunction with simple 1D models into the nowcasting process. 

Abstract 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 K DP , 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. 

Raindrop size distributions (DSD) and rain rate have been estimated from polarimetric radar data using different approaches with the accuracy depending on the errors both in the radar measurements and the estimation methods. Herein, a deep neural network (DNN) technique was utilized to improve the estimation of the DSD and rain rate by mitigating these errors. The performance of this approach was evaluated using measurements from a two-dimensional video disdrometer (2DVD) at the Kessler Atmospheric and Ecological Field Station in Oklahoma as ground truth with the results compared against conventional estimation methods for the period 2006–17. Physical parameters (mass-/volume-weighted diameter and liquid water content), rain rate, and polarimetric radar variables (including radar reflectivity and differential reflectivity) were obtained from the DSD data. Three methods—physics-based inversion, empirical formula, and DNN—were applied to two different temporal domains (instantaneous and rain-event average) with three diverse error assumptions (fitting, measurement, and model errors). The DSD retrievals and rain estimates from 18 cases were evaluated by calculating the bias and root-mean-squared error (RMSE). DNN produced the best performance for most cases, with up to a 5% reduction in RMSE when model errors existed. DSD and rain estimated from a nearby polarimetric radar using the empirical and DNN methods were well correlated with the disdrometer observations; the rain-rate estimate bias of the DNN was significantly reduced (3.3% in DNN vs 50.1% in empirical). These results suggest that DNN has advantages over the physics-based and empirical methods in retrieving rain microphysics from radar observations.

 

Assimilating radar reflectivity into convective-scale NWP models remains a challenging topic in radar data assimilation. A primary reason is that the reflectivity forward observation operator is highly nonlinear. To address this challenge, a power transformation function is applied to the WRF Model’s hydrometeor and water vapor mixing ratio variables in this study. Three 3D variational data assimilation experiments are performed and compared for five high-impact weather events that occurred in 2019: (i) a control experiment that assimilates reflectivity using the original hydrometeor mixing ratios as control variables, (ii) an experiment that assimilates reflectivity using power-transformed hydrometeor mixing ratios as control variables, and (iii) an experiment that assimilates reflectivity and retrieved pseudo–water vapor observations using power-transformed hydrometeor and water vapor mixing ratios (q_υ) as control variables. Both qualitative and quantitative evaluations are performed for 0–3-h forecasts from the five cases. The analysis and forecast performance in the two experiments with power-transformed mixing ratios is better than the control experiment. Notably, the assimilation of pseudo–water vapor with power-transformedq_υas an additional control variable is found to improve the performance of the analysis and short-term forecasts for all cases. In addition, the convergence rate of the cost function minimization for the two experiments that use the power transformation is faster than that of the control experiments.

The effective use of radar reflectivity observations in any data assimilation scheme remains an important research topic because reflectivity observations explicitly include information about hydrometeors and also implicitly include information about the distribution of moisture within storms. However, it is difficult to assimilate reflectivity because the reflectivity forward observation operator is highly nonlinear. This study seeks to identify a more effective way to assimilate reflectivity into a convective-scale NWP model to improve the accuracy of predictions of high-impact weather events.