Search

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. [ABSTRACT FROM AUTHOR]

Fine-scale characteristics of summer precipitation over Cang Mountain, a long and narrow mountain with a quasi-north–south orientation in Southwest China, are studied using station and radar data. Three kinds of rainfall processes are classified according to the initial stations of regional rainfall events (RREs) by utilizing minute-scale rain gauge data. RREs initiating in the western part of Cang Mountain exhibit eastward evolution and tend to reach their maximum rainfall intensity on the mountaintop. The results indicate differences in the precipitation evolution characteristics between short-duration (1–3 h) and long-duration (at least 6 h) events. Short-duration events begin farther from the mountaintop and then propagate eastward, whereas long-duration events remain longer around the mountaintop. RREs that initiate from the eastern part of Cang Mountain display westward propagation and frequently reach their maximum rainfall intensity over the eastern slope of the mountain. Among them, short-duration events tend to propagate farther west of Cang Mountain at high speeds, but the westward evolution of long-duration events is mainly confined to the eastern part of Cang Mountain. For mountaintop-originated RREs, precipitation quickly reaches its maximum intensity after it starts and then continues for a long time around the mountaintop during the period from late afternoon to early morning. These findings provide references for the fine-scale prediction of precipitation evolution in small-scale mountainous areas. [ABSTRACT FROM AUTHOR]

The Integrated Multisatellite Retrievals for Global Precipitation Measurement Mission (IMERG) is a global precipitation product that uses precipitation retrievals from the virtual constellation of satellites with passive microwave (PMW) sensors, as available. In the absence of PMW observations, IMERG uses a Kalman filter scheme to morph precipitation from one PMW observation to the next. In this study, an analysis of convective systems observed during the Convective Process Experiment (CPEX) suggests that IMERG precipitation depends more strongly on the availability of PMW observations than previously suspected. Following this evidence, we explore systematic biases in IMERG through bulk statistics. In two CPEX case studies, cloud photographs, pilot's radar, and infrared imagery suggest that IMERG represents the spatial extent of precipitation relatively well when there is a PMW observation but sometimes produces spurious precipitation areas in the absence of PMW observations. Also, considering an observed convective system as a precipitation object in IMERG, the maximum rain rate peaked during PMW overpasses, with lower values between them. Bulk statistics reveal that these biases occur throughout IMERG Version 06.Wefind that locations and times without PMW observations have a higher frequency of light precipitation rates and a lower frequency of heavy precipitation rates due to retrieval artifacts. These results reveal deficiencies in the IMERG Kalman filter scheme, which have led to the development of the Scheme for Histogram Adjustment with Ranked Precipitation Estimates in the Neighborhood (SHARPEN; described in a companion paper) that will be applied in the next version of IMERG. [ABSTRACT FROM AUTHOR]

Flood forecasting is essential to minimize the impacts and costs of floods, especially in urbanized watersheds. Radar rainfall estimates are becoming increasingly popular in flood forecasting because they provide the much-needed realtime spatially distributed precipitation information. The current study evaluates the use of radar quantitative precipitation estimates (QPEs) in hydrological model calibration for streamflow simulation and flood mapping in an urban setting. First, S-band and C-band radar QPEs were integrated into event-based hydrological models to improve the calibration of model parameters. Second, rain gauge and radar precipitation estimates' performances were compared for hydrological modeling in an urban watershed to assess radar QPE's effects on streamflow simulation accuracy. Third, flood extent maps were produced using coupled hydrological--hydraulic models integrated within the Hydrologic Engineering Center Real-Time Simulation (HEC-RTS) framework. It is shown that the bias correction of radar QPEs can enhance the hydrological model calibration. The radar--gauge merging obtained a Kling--Gupta efficiency, modified peak flow criterion, Nash--Sutcliffe efficiency, and volume error improvement of about +0.42, +0.12, +0.78, and -0.23, respectively, for S-band and +0.64, +0.36, +1.12, and -0.34, respectively, for C-band radar QPEs. Merged radar QPEs are also helpful in running hydrological models calibrated using gauge data. The HEC-RTS framework can be used to produce flood forecast maps using the bias-corrected radar QPEs. Therefore, radar rainfall estimates could be efficiently used to forecast floods in urbanized areas for effective flood management and mitigation. Canadian flood forecasting systems could be efficiently updated by integrating bias-corrected radar QPEs to simulate streamflow and produce flood inundation maps. [ABSTRACT FROM AUTHOR]

Decades of research have investigated processes that contribute to downburst development, as well as identified precursor radar signatures that can accompany these events. These advancements have increased downburst predictability, but downbursts still pose a significant forecast challenge, especially in low-shear environments that produce short-lived single and multicell thunderstorms. Additional information provided by dual-polarization radar data may prove useful in anticipating downburst development. One such radar signature is the KDP core (where KDP is the specific differential phase), which can indicate processes such as melting and precipitation loading that increase negative buoyancy and can result in downburst development. Therefore, KDP cores associated with 81 different downbursts across 10 states are examined to explore if this signature could provide forecasters with a reliable and useable downburst precursor signature. The KDP core characteristics near the environmental melting layer, vertical gradients of KDP, and environmental conditions were quantified to identify any differences between downbursts of varying intensities. The analysis shows that 1) KDP cores near the environmental melting layer are a reliable signal that a downburst will develop; 2) while using KDP cores to anticipate an impending downburst's intensity is challenging, larger KDP near the melting layer and larger vertical gradients of KDP are more commonly associated with strong downbursts than weak ones; 3) downbursts occurring in environments with less favorable conditions for downbursts are associated with higher magnitude KDP cores, and 4) KDP cores evolve relatively slowly (typically longer than 15 min), which makes them easily observable with the 5-min volumetric updates currently provided by operational radars. [ABSTRACT FROM AUTHOR]

Range-oversampling processing is a technique that can be used to lower the variance of radar-variable estimates, reduce radar update times, or a mixture of both. There are two main assumptions for using range-oversampling processing: accurate knowledge of the range correlation and uniform reflectivity in the radar resolution volume. The first assumption has been addressed in previous research; this work focuses on the uniform reflectivity assumption. Earlier research shows that significant reflectivity gradients can occur in storms; we utilized those results to develop realistic simulations of radar returns that include effects of reflectivity gradients in range. An important consideration when using range-oversampling processing is the resulting change in the range weighting function. The range weighting function can change for different types of range-oversampling processing, and some techniques, such as adaptive pseudowhitening, can lead to different range weighting functions at each range gate. To quantify the possible effects of differing range weighting functions in the presence of reflectivity gradients, we developed simulations to examine varying types of range-oversampling processing with two receiver filters: a matched receiver filter and a wider-bandwidth receiver filter (as recommended for use with range oversampling). Simulation results show that differences in range weighting functions are the only contributor to differences in radar reflectivity measurements. Results from real weather data demonstrate that the reflectivity gradients that occur in typical severe storms do not cause significant changes in reflectivity measurements and that the benefits from range-oversampling processing outweigh the possible isolated effects from large reflectivity gradients. [ABSTRACT FROM AUTHOR]

The retrieval of the mass-weighted mean diameter (Dm) is a fundamental component of spaceborne precipitation retrievals. The Dual-Frequency Precipitation Radar (DPR) on the Global Precipitation Measurement (GPM) satellite is the first satellite to use dual-wavelength ratio measurements—the quotient of radar reflectivity factors (Z) measured at Ku and Ka wavelengths—to retrieve Dm. While it is established that DWR, being theoretically insensitive to changes in ice crystal mass and concentration, can provide a superior retrieval of Dm compared to Z-based retrievals, the benefits of this retrieval have yet to be directly observed or quantified. In this study, DWR–Dm and Z–Dm relationships are empirically generated from collocated airborne radar and in situ cloud particle probe measurements. Data are collected during nine intensive observation periods (IOPs) from three experiments representing different locations and times of year. Across IOPs with varying ice crystal concentrations, cloud temperatures, and storm types, Z–Dm relationships vary considerably while the DWR–Dm relationship remains consistent. This study confirms that a DWR–Dm relationship can provide a more accurate and consistent Dm retrieval than a Z–Dm relationship, quantified by a reduced overall RMSE (0.19 and 0.25 mm, respectively) and a reduced range of biases between experiments (0.11 and 0.32 mm, respectively). [ABSTRACT FROM AUTHOR]

The Flexible Array of Radars and Mesonets (FARM) Facility is an extensive mobile/quickly deployable (MQD) multiple-Doppler radar and in situ instrumentation network. The FARM includes four radars: two 3-cm dual polarization, dual frequency (DPDF), Doppler on Wheels (DOW6/DOW7), the Rapid-Scan DOW (RSDOW), and a quickly deployable (QD) DPDF 5-cm C band on Wheels (COW). The FARM includes three mobile mesonet (MM) vehicles with 3.5-m masts, an array of rugged QD weather stations (PODNET), QD weather stations deployed on infrastructure such as light/power poles (POLENET), four disdrometers, six MQD upper-air sounding systems and a Mobile Operations and Repair Center (MORC). The FARM serves a wide variety of research/educational uses. Components have deployed to >30 projects during 1995–2020 in the United States, Europe, and South America, obtaining pioneering observations of a myriad of small spatial- and temporal-scale phenomena including tornadoes, hurricanes, lake-effect snow storms, aircraft-affecting turbulence, convection initiation, microbursts, intense precipitation, boundary layer structures and evolution, airborne hazardous substances, coastal storms, wildfires and wildfire suppression efforts, weather modification effects, and mountain/alpine winds and precipitation. The radars and other FARM systems support innovative educational efforts, deploying >40 times to universities/colleges, providing hands-on access to cutting-edge instrumentation for their students. The FARM provides integrated multiple radar, mesonet, sounding, and related capabilities enabling diverse and robust coordinated sampling of three-dimensional vector winds, precipitation, and thermodynamics increasingly central to a wide range of mesoscale research. Planned innovations include S-band on Wheels Network (SOWNET) and Bistatic Adaptable Radar Network (BARN), offering more qualitative improvements to the field project observational paradigm, providing broad, flexible, and inexpensive 10-cm radar coverage and vector wind field measurements. [ABSTRACT FROM AUTHOR]

A convective-scale static background-error covariance (BEC) matrix is further developed to include the capability of direct reflectivity assimilation and evaluated within the GSI-based three-dimensional variational (3DVar) and hybrid ensemble–variational (EnVar) methods. Specific developments are summarized as follows: 1) Control variables (CVs) are extended to include reflectivity, vertical velocity, and all hydrometeor types. Various horizontal momentum and moisture CV options are included. 2) Cross correlations between all CVs are established. 3) A storm intensity-dependent binning method is adopted to separately calculate static error matrices for clear-air and storms with varying intensities. The resultant static BEC matrices are simultaneously applied at proper locations guided by the observed reflectivity. 4) The EnVar is extended to adaptively incorporate static BECs based on the quality of ensemble covariances. Evaluation and examination of the new static BECs are first performed on the 8 May 2003 Oklahoma City supercell. Detailed diagnostics and 3DVar examinations suggest zonal/meridional winds and pseudo–relative humidity are selected as horizontal momentum and moisture CVs for direct reflectivity assimilation, respectively; inclusion of cross correlations favors spin up and maintains the analyzed storms; application of binning improves characteristics and persistence of the simulated storm. Relative to an experiment using the full ensemble BECs (Exp-PureEnVar), incorporating static BECs in hybrid EnVar reduces spinup time and better analyzes reflectivity distributions while the background ensemble is deficient in sampling errors. Compared to both pure 3DVar and Exp-PureEnVar, hybrid EnVar better predicts reflectivity distributions and better maintains a strong mesocyclone. Further examination through the 20 May 2013 Oklahoma supercells confirms these results and additionally demonstrates the effectiveness of adaptive hybridization. [ABSTRACT FROM AUTHOR]

Applications like drought monitoring and forecasting can profit from the global and near-real-time availability of satellite-based precipitation estimates once their related uncertainties and challenges are identified and treated. To this end, this study evaluates the IMERG V06B Late Run precipitation product from the Global Precipitation Measurement mission (GPM), a multisatellite product that combines space-based radar, passive microwave (PMW), and infrared (IR) data into gridded precipitation estimates. The evaluation is performed on the spatiotemporal resolution of IMERG (0.1° × 0.1°, 30 min) over the Netherlands over a 5-yr period. A gauge-adjusted radar precipitation product from the Royal Netherlands Meteorological Institute (KNMI) is used as reference, against which IMERG shows a large positive bias. To find the origin of this systematic overestimation, the data are divided into seasons, rainfall intensity ranges, echo top height (ETH) ranges, and categories based on the relative contributions of IR, morphing, and PMW data to the IMERG estimates. Furthermore, the specific radiometer is identified for each PMW-based estimate. IMERG's detection performance improves with higher ETH and rainfall intensity, but the associated error and relative bias increase as well. Severe overestimation occurs during low-intensity rainfall events and is especially linked to PMW observations. All individual PMW instruments show the same pattern: overestimation of low-intensity events and underestimation of high-intensity events. IMERG misses a large fraction of shallow rainfall events, which is amplified when IR data are included. Space-based retrieval of shallow and low-intensity precipitation events should improve before IMERG can become accurate over the middle and high latitudes. [ABSTRACT FROM AUTHOR]

It remains a challenge to provide accurate and timely flood warnings in many parts of the western United States. As part of the Advanced Quantitative Precipitation Information (AQPI) project, this study explores the potential of using the AQPI gap-filling radar network for streamflow simulation of selected storm events in the San Francisco Bay Area under a WRF-Hydro modeling system. Two types of watersheds including natural and human-affected among the most flood-prone region of the Bay Area are investigated. Based on the high-resolution AQPI X-band radar rainfall estimates, three basic routing configurations, including Grid, Reach, and National Water Model (NWM), are used to quantify the impact of different model physics options on the simulated streamflow. It is found that the NWM performs better in terms of reproducing streamflow volumes and hydrograph shapes than the other routing configurations when reservoirs exist in the watershed. Additionally, the AQPI X-band radar rainfall estimates (without gauge correction) provide reasonable streamflow simulations, and they show better performance in reproducing the hydrograph peaks compared with the gauge-corrected rainfall estimates based on the operational S-band Next Generation Weather Radar network. Also, a sensitivity test reveals that surficial conditions have a significant influence on the streamflow simulation during the storm: the discharge increases to a higher level as the infiltration factor (REFKDT) decreases, and its peak goes down and lags as surface roughness coefficient (Mann) increases. The time delay analysis of precipitation input on the streamflow at the two outfalls of the surveyed watersheds further demonstrates the link between AQPI gap-filling radar observations and streamflow changes in this urban region. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

NASA Precipitation Measurement Mission observations are used to evaluate the diurnal cycle of precipitation from three CMIP6 models (NCAR-CESM2, CNRM-CM6.1, CNRM-ESM2.1) and the ERA5 reanalysis. NASA's global-gridded IMERG product, which combines spaceborne microwave radiometer, infrared sensor, and ground-based gauge measurements, provides high-spatiotemporal-resolution (0.1° and half-hourly) estimates that are suitable for evaluating the diurnal cycle in models, as determined against the ground-based radar network over the conterminous United States. IMERG estimates are coarsened to the spatial and hourly resolution of the state-of-the-art CMIP6 and ERA5 products, and their diurnal cycles are compared across multiple decades of June–August in the 60°N–60°S domain (IMERG and ERA5: 2000–19; NCAR and CNRM: 1979–2008). Low-precipitation regions (and weak-amplitude regions when analyzing the diurnal phase) are excluded from analyses so as to assess only robust diurnal signals. Observations identify greater diurnal amplitudes over land (26%–134% of the precipitation mean; 5th–95th percentile) than over ocean (14%–66%). ERA5, NCAR, and CNRM underestimate amplitudes over ocean, and ERA5 overestimates over land. IMERG observes a distinct diurnal cycle only in certain regions, with precipitation peaking broadly between 1400 and 2100 LST over land (2100–0600 LST over mountainous and varying-terrain regions) and 0000 and 1200 LST over ocean. The simulated diurnal cycle is unrealistically early when compared with observations, particularly over land (NCAR-CESM2 AMIP: −1 h; ERA5: −2 h; CNRM-CM6.1 AMIP: −4 h on average) with nocturnal maxima not well represented over mountainous regions. Furthermore, ERA5's representation of the diurnal cycle is too simplified, with less interannual variability in the time of maximum relative to observations over many regions. SIGNIFICANCE STATEMENT: Identifying and addressing climate model errors in representing the diurnal cycle of precipitation are critical to improving their accuracy. This study provides an update on the diurnal cycle performance of state-of-the-art climate models and reanalysis against state-of-the-art satellite observations. The models and reanalysis have varying biases in diurnal amplitude over land, where amplitudes are stronger, and they underestimate amplitudes over ocean. They also simulate precipitation over land to peak too early in the day, from −1 to −4 h on average depending on the model. Nocturnal maxima in mountainous regions are not well simulated, although the reanalysis outperforms the models in this case. Future work can use these findings to improve realism in the next generation of climate models. [ABSTRACT FROM AUTHOR]

The Method of Object-based Diagnostic Evaluation (MODE) is used to perform an object-based verification of approximately 1400 forecasts of composite reflectivity from the operational HRRR during April-September 2019. In this study, MODE is configured to prioritize deep, moist convective storm cells typical of those that produce severe weather across the central and eastern United States during the warm season. In particular, attributes related to distance and size are given the greatest attribute weights for computing interest in MODE. HRRR tends to overforecast all objects, but substantially overforecasts both small objects at low-reflectivity thresholds and large objects at high-reflectivity thresholds. HRRR tends to either underforecast objects in the southern and central plains or has a correct frequency bias there, whereas it overforecasts objects across the southern and eastern United States. Attribute comparisons reveal the inability of the HRRR to fully resolve convective-scale features and the impact of data assimilation and loss of skill during the initial hours of the forecasts. Scalar metrics are defined and computed based on MODE output, chiefly relying on the interest value. The object-based threat score (OTS), in particular, reveals similar performance of HRRR forecasts as does the Heidke skill score, but with differing magnitudes, suggesting value in adopting an object-based approach to forecast verification. The typical distance between centroids of objects is also analyzed and shows gradual degradation with increasing forecast length. SIGNIFICANCE STATEMENT: Improving weather forecast models requires determining where the model does well and where it does not. Gridpoint-based methods for assessing model forecasts have known shortfalls when applied to high-resolution models that can forecast individual thunderstorms. We present an object-based verification procedure that focuses on identifying actual meteorological features such as thunderstorms instead of gridpoint-by-gridpoint comparison between forecasts and verifying truth. This article reveals some of the information ascertained from this assessment and illustrates the enhancement of information obtained from object-based verification to gridpoint-based assessment. [ABSTRACT FROM AUTHOR]

The Royal Netherlands Meteorological Institute (KNMI) operates two operational dual-polarization C-band weather radars providing 2D radar rainfall products. Attenuation can result in severe underestimation of rainfall amounts, particularly in convective situations that are known to have high impact on society. To improve the radar-based precipitation estimates, two attenuation correction methods are evaluated and compared: 1) modified Kraemer (MK) method, i.e., Hitschfeld–Bordan where parameters of the power-law Zh–kh relation are adjusted such that reflectivities in the entire dataset do not exceed 59 dBZh and attenuation correction is limited to 10 dB; and 2) using two-way path-integrated attenuation computed from the dual-polarization moment specific differential phase Kdp (Kdp method). In both cases the open-source Python library wradlib is employed for the actual attenuation correction. A radar voxel only contributes to the computed path-integrated attenuation if its height is below the forecasted freezing-level height from the numerical weather prediction model HARMONIE-AROME. The methods are effective in improving hourly and daily quantitative precipitation estimation (QPE), where the Kdp method performs best. The verification against rain gauge data shows that the underestimation diminishes from 55% to 37% for hourly rainfall for the Kdp method when the gauge indicates more than 10 mm of rain in that hour. The improvement for the MK method is less pronounced, with a resulting underestimation of 40%. The stability of the MK method holds a promise for application to data archives from single-polarization radars. [ABSTRACT FROM AUTHOR]

In this work, we present a new quantitative precipitation estimation (QPE) quality-control (QC) algorithm for the U.K. weather radar network. The real-time adaptive algorithm uses a neural network (NN) to select data from the lowest useable elevation scan to optimize the combined performance of two other radar data correction algorithms: ground-clutter mitigation [using Clutter Environment Analysis using Adaptive Processing (CLEAN-AP)] and vertical profile of reflectivity (VPR) correction. The NN is trained using 3D tiles of observed uncontaminated weather signals that are systematically combined with ground-clutter signals collected under dry weather conditions. This approach provides a way to simulate radar signals with a wide range of clutter contamination conditions and with realistic spatial structures while providing the uncontaminated "truth" with respect to which the performance of the QC algorithm can be measured. An evaluation of QPE products obtained with the proposed QC algorithm demonstrates superior performance as compared to those obtained with the QC algorithm currently used in operations. Similar improvements are also illustrated using radar observations from two periods of prolonged precipitation, showing a better balance between overestimation errors from using clutter-contaminated low-elevation radar data and VPR-induced errors from using high-elevation radar data. [ABSTRACT FROM AUTHOR]

Using 17-yr spaceborne precipitation radar measurements, this study investigates how diurnal cycles of rainfall and convective characteristics over the South China Sea region are modulated by the boreal summer intraseasonal oscillation (BSISO). Generally, diurnal cycles change significantly between suppressed and active BSISO periods. Over the Philippines and Indochina, where the low-level monsoon flows impinge on coast lines, diurnal cycles of rainfall and many convective properties are enhanced during suppressed periods. During active periods, diurnal variation of convection is still significant over land but diminishes over water. Also, afternoon peaks of rainfall and MCS populations over land are obviously extended in active periods, mainly through the enhancement of stratiform precipitation. Over Borneo, where the prevailing low-level winds are parallel to coasts, diurnal cycles (both onshore and offshore) are actually stronger during active periods. Radar profiles also demonstrate a pronounced nocturnal offshore propagation of deep convection over western Borneo in active periods. During suppressed periods, coastal afternoon convection over Borneo is reduced, and peak convection occurs over the mountains until the convective suppression is overcome in the late afternoon or evening. A major portion (>70%) of the total precipitation over the Philippines and Indochina during suppressed periods falls from afternoon isolated to medium-sized systems (<10 000 km²), but more than 70% of the active BSISO rainfall is contributed by nocturnal (after 1800 LT) broad precipitation systems (>10 000 km²). However, offshore total precipitation is dominated by large precipitation systems (>10 000 km²) regardless of BSISO phases and regions. [ABSTRACT FROM AUTHOR]

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 (08C) 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 08, 2108, and above 2108C; 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. [ABSTRACT FROM AUTHOR]

The NASA Goddard Space Flight Center's (GSFC's) W-band (94 GHz) Cloud Radar System (CRS) has been comprehensively updated to modern solid-state and digital technology. This W-band (94 GHz) radar flies in nadir-pointing mode on the NASA ER-2 high-altitude aircraft, providing polarimetric reflectivity and Doppler measurements of clouds and precipitation. This paper describes the design and signal processing of the upgraded CRS. It includes details on the hardware upgrades [solid-state power amplifier (SSPA) transmitter, antenna, and digital receiver] including a new reflectarray antenna and solid-state transmitter. It also includes algorithms, including internal loop-back calibration, external calibration using a direct relationship between volume reflectivity and the range-integrated backscatter of the ocean, and a modified staggered–pulse repetition frequency (PRF) Doppler algorithm that is highly resistant to unfolding errors. Data samples obtained by upgraded CRS through recent NASA airborne science missions are provided. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In this study, we propose a new way to obtain motion vectors using the integrating velocity–azimuth process (IVAP) method for extrapolation nowcasting. Traditional tracking methods rely on tracking radar echoes of a few time slices. In contrast, the IVAP method does not depend on the past variation of radar echoes; it only needs the radar echo and radial velocity observations at the latest time. To demonstrate it is practical to use IVAP-retrieved winds to extrapolate radar echoes, we carried out nowcasting experiments using the IVAP method, and compared these results with the results using a traditional method, namely, the tracking radar echoes by correlation (TREC) method. Comparison based on a series of large-scale mature rainfall cases showed that the IVAP method has similar accuracy to that of the TREC method. In addition, the IVAP method provides the vertical wind profile that can be used to anticipate storm type and motion deviations. [ABSTRACT FROM AUTHOR]

A method based on machine learning and image processing techniques has been developed to track the surface expression of internal waves in near–real time. X-band radar scans are first preprocessed and averaged to suppress surface wave clutter and enhance the signal-to-noise ratio of persistent backscatter features driven by gradients in surface currents. A machine learning algorithm utilizing a support vector machine (SVM) model is then used to classify whether or not the image contains an internal solitary wave (ISW) or internal tide bore (bore). The use of machine learning is found to allow rapid assessment of the large dataset, and provides insight on characterizing optimal environmental conditions to allow for radar illumination and detection of ISWs and bores. Radon transforms and local maxima detections are used to locate these features within images that are determined to contain an ISW or bore. The resulting time series of locations is used to create a map of propagation speed and direction that captures the spatiotemporal variability of the ISW or bore in the coastal environment. This technique is applied to 2 months of data collected near Point Sal, California, and captures ISW and bore propagation speed and direction information that currently cannot be measured with instruments such as moorings and synthetic aperture radar (SAR). [ABSTRACT FROM AUTHOR]

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, σ e ⁡ (K DP , Z) = γ K DP α Z β . 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. [ABSTRACT FROM AUTHOR]

A robust and computationally efficient object tracking algorithm is developed by incorporating various tracking techniques. Physical properties of the objects, such as brightness temperature or reflectivity, are not considered. Therefore, the algorithm is adaptable for tracking convection-like features in simulated data and remotely sensed two-dimensional images. In this algorithm, a first guess of the motion, estimated using the Fourier phase shift, is used to predict the candidates for matching. A disparity score is computed for each target–candidate pair. The disparity also incorporates overlapping criteria in the case of large objects. Then the Hungarian method is applied to identify the best pairs by minimizing the global disparity. The high-disparity pairs are unmatched, and their target and candidate are declared expired and newly initiated objects, respectively. They are tested for merger and split on the basis of their size and overlap with the other objects. The sensitivity of track duration is shown for different disparity and size thresholds. The paper highlights the algorithm's ability to study convective life cycles using radar and simulated data over Darwin, Australia. The algorithm skillfully tracks individual convective cells (a few pixels in size) and large convective systems. The duration of tracks and cell size are found to be lognormally distributed over Darwin. The evolution of size and precipitation types of isolated convective cells is presented in the Lagrangian perspective. This algorithm is part of a vision for a modular platform [viz., TINT is not TITAN (TINT) and Tracking and Object-Based Analysis of Clouds (tobac)] that will evolve into a sustainable choice to analyze atmospheric features. [ABSTRACT FROM AUTHOR]

Assessments of spatiotemporal severe hailfall characteristics using hail reports are plagued by serious limitations in report databases, including biases in reported sizes, occurrence time, and location. Multiple studies have used Next Generation Weather Radar (NEXRAD) network observations or environmental hail proxies from reanalyses. Previous work has specifically utilized the single-polarization radar parameter maximum expected size of hail (MESH). In addition to previous work being temporally limited, updates are needed to include recent improvements that have been made to MESH. This study aims to quantify severe hailfall characteristics during a 23-yr period, markedly longer than previous studies, using both radar observations and reanalysis data. First, the improved MESH configuration is applied to the full archive of gridded hourly radar observations known as GridRad (1995–2017). Next, environmental constraints from the Modern-Era Retrospective Analysis for Research and Applications, version 2, are applied to the MESH distributions to produce a corrected hailfall climatology that accounts for the reduced likelihood of hail reaching the ground. Spatial, diurnal, and seasonal patterns show that in contrast to the report climatology indicating one high-frequency hail maximum centered on the Great Plains, the MESH-only method characterizes two regions: the Great Plains and the Gulf Coast. The environmentally filtered MESH climatology reveals improved agreement between report characteristics (frequency, location, and timing) and the recently improved MESH calculation methods, and it reveals an overall increase in diagnosed hail days and westward broadening in the spatial maximum in the Great Plains than that seen in reports. [ABSTRACT FROM AUTHOR]

Cloud microphysics significantly impact tropical cyclone precipitation. A prior polarimetric radar observational study by the authors revealed the ice-phase microphysical processes as the dominant microphysics mechanisms responsible for the heavy precipitation in the outer rainband of Typhoon Nida (2016). To assess the model performance regarding microphysics, three double-moment microphysics schemes (i.e., Thompson, Morrison, and WDM6) are evaluated by performing a set of simulations of the same case. While these simulations capture the outer rainband's general structure, microphysics in the outer rainbands are strikingly different from the observations. This discrepancy is primarily attributed to different microphysics parameterizations in these schemes, rather than the differences in large-scale environments due to cloud–environment interactions. An interesting finding in this study is that the surface rain rate or liquid water content is inversely proportional to the simulated mean raindrop sizes. The mass-weighted raindrop diameters are overestimated in the Morrison and Thompson schemes and underestimated in the WDM6 scheme, while the former two schemes produce lower liquid water content than WDM6. Compared with the observed ice water content based on a new polarimetric radar retrieval method, the ice water content above the environmental 0°C level in all simulations is highly underestimated, especially at heights above 12 km MSL where large concentrations of small ice particles are typically prevalent. This finding suggests that the improper treatment of ice-phase processes is potentially an important error source in these microphysics schemes. Another error source identified in the WDM6 scheme is overactive warm-rain processes that produce excessive concentrations of smaller raindrops. [ABSTRACT FROM AUTHOR]

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 thanDPR_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 1mmh-1 and underestimation when gauge estimates are larger than 1mmh-1. The DPR_HS and DPR_MS agree better with gauge estimates below and above 2.5mmh-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. [ABSTRACT FROM AUTHOR]

Conventionally, Micro Rain Radars (MRRs) have been used as a tool to calibrate reflectivity from weather radars, estimate the relation between rainfall rate and reflectivity, and study microphysical processes in precipitation. However, limited attention has been given to the reliability of the retrieved drop size distributions (DSDs) from MRRs. This study sheds more light on this aspect by examining the sensitivity of retrieved DSDs to the assumptions made to map Doppler spectra into size distributions, and investigates the capability of an MRR to assess polarimetric observations from operational weather radars. For that, an MRR was installed near the Cabauw observatory in the Netherlands, between the International Research Center for Telecommunications and Radar (IRCTR) Drizzle Radar (IDRA) X-band radar and the Herwijnen operational C-band radar. The measurements of the MRR from November 2018 to February 2019 were used to retrieve DSDs and simulate horizontal reflectivity Ze, differential reflectivity ZDR, and specific differential phase KDP in rain. Attention is given to the impact of aliased spectra and right-hand-side truncation on the simulation of polarimetric variables. From a quantitative assessment, the correlations of Ze and ZDR between theMRRand Herwijnen radar were 0.93 and 0.70, respectively, while those between the MRR and IDRA were 0.91 and 0.69. However, Ze and ZDR from the Herwijnen radar showed slight biases of 1.07 and 0.25 dB. For IDRA, the corresponding biases were 2.67 and20.93 dB. Our results show that MRR measurements are advantageous to inspect the calibration of scanning radars and validate polarimetric estimates in rain, provided that the DSDs are correctly retrieved and controlled for quality assurance. [ABSTRACT FROM AUTHOR]

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 operationalWSR-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 performvery 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. [ABSTRACT FROM AUTHOR]

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 D m l and the liquid equivalent normalized intercept parameter N w l . 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 D m l , N w l , 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 D m l , N w l , 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. [ABSTRACT FROM AUTHOR]

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

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. [ABSTRACT FROM AUTHOR]

Two power-law relations linking equivalent radar reflectivity factor Ze and snowfall rate S are derived for a K-band Micro Rain Radar (MRR) and for a W-band cloud radar. For the development of these Ze–S relationships, a dataset of calculated and measured variables is used. Surface-based video-disdrometer measurements were collected during snowfall events over five winters at the high-latitude site in Hyytiälä, Finland. The data from 2014 to 2018 include particle size distributions (PSD) and their fall velocities, from which snowflake masses were derived. The K- and W-band Ze values are computed using these surface-based observations and snowflake scattering properties as provided by T-matrix and single-particle scattering tables, respectively. The uncertainty analysis shows that the K-band snowfall-rate estimation is significantly improved by including the intercept parameter N0 of the PSD calculated from concurrent disdrometer measurements. If N0 is used to adjust the prefactor of the Ze–S relationship, the RMSE of the snowfall-rate estimate can be reduced from 0.37 to around 0.11 mm h−1. For W-band radar, a Ze–S relationship with constant parameters for all available snow events shows a similar uncertainty when compared with the method that includes the PSD intercept parameter. To demonstrate the performance of the proposed Ze–S relationships, they are applied to measurements of the MRR and the W-band microwave radar for Arctic clouds at the Arctic research base operated by the German Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI) and the French Polar Institute Paul Emile Victor (IPEV) (AWIPEV) in Ny-Ålesund, Svalbard, Norway. The resulting snowfall-rate estimates show good agreement with in situ snowfall observations while other Ze–S relationships from literature reveal larger differences. [ABSTRACT FROM AUTHOR]

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 CVlogq for the analyses and forecasts of five convective storm cases from the spring of 2017. Results show that CVpq with p 5 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 CVlogq. Compared to CVq and CVlogq, 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. [ABSTRACT FROM AUTHOR]

High-resolution datasets offer the potential to improve our understanding of spatial and temporal precipitation patterns and storm structures. The goal of this study is to evaluate the similarities and differences of object-based storm characteristics as observed using space- or land-based sensors. The Method of Object-based Diagnostic Evaluation (MODE) Time Domain (MTD) is used to identify and track storm objects in two high-resolution merged datasets: the Integrated Multisatellite Retrievals for Global Precipitation Measurement (IMERG) final product V06B and gauge-corrected ground-radar-based Multi-Radar Multi-Sensor (MRMS) quantitative precipitation estimations. Characteristics associated with landfalling hurricanes were also examined as a separate category of storm. The results reveal that IMERG and MRMS agree reasonably well across many object-based storm characteristics. However, there are some discrepancies that are statistically significant. MRMS storms are more concentrated, with smaller areas and higher peak intensities, which implies higher flash flood risks associated with the storms. On the other hand, IMERG storms can travel longer distances with a higher volume of precipitation, which implies higher risk of riverine flooding. Agreement between the datasets is higher for faster-moving hurricanes in terms of the averaged intensity. Finally, MRMS indicates a higher average precipitation intensity during the hurricane's lifetime. However, in non-hurricanes, the opposite result was observed. This is likely related to MRMS having higher resolution; monitoring the hurricanes from many viewing angles, leading to different signal saturation properties compared to IMERG; and/or the dominance of droplet aggregation effects over evaporation effects at lower altitudes. [ABSTRACT FROM AUTHOR]

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). [ABSTRACT FROM AUTHOR]

Techniques to mitigate analysis errors arising from the nonsimultaneity of data collections typically use advection-correction procedures based on the hypothesis (frozen turbulence) that the analyzed field can be represented as a pattern of unchanging form in horizontal translation. It is more difficult to advection correct the radial velocity than the reflectivity because even if the vector velocity field satisfies this hypothesis, its radial component does not—but that component does satisfy a second-derivative condition. We treat the advection correction of the radial velocity (υr) as a variational problem in which errors in that second-derivative condition are minimized subject to smoothness constraints on spatially variable pattern-translation components (U, V). The Euler–Lagrange equations are derived, and an iterative trajectory-based solution is developed in which U, V, and υr are analyzed together. The analysis code is first verified using analytical data, and then tested using Atmospheric Imaging Radar (AIR) data from a band of heavy rainfall on 4 September 2018 near El Reno, Oklahoma, and a decaying tornado on 27 May 2015 near Canadian, Texas. In both cases, the analyzed υr field has smaller root-mean-square errors and larger correlation coefficients than in analyses based on persistence, linear time interpolation, or advection correction using constant U and V. As some experimentation is needed to obtain appropriate parameter values, the procedure is more suitable for non-real-time applications than use in an operational setting. In particular, the degree of spatial variability in U and V, and the associated errors in the analyzed υr field are strongly dependent on a smoothness parameter. [ABSTRACT FROM AUTHOR]

Polarimetric radar observations from the NEXRAD WSR-88D operational radar network in the contiguous United States, routinely available since 2013, are used to reveal three prominent microphysical signatures in landfalling tropical cyclones: 1) hydrometeor size sorting within the eyewall convection, 2) vertical displacement of the melting layer within the inner core, and 3) dendritic growth layers within stratiform regions of the inner core. Size sorting signatures within eyewall convection are observed with greater frequency and prominence in more intense hurricanes, and are observed predominantly within the deep-layer environmental wind shear vector-relative quadrants that harbor the greatest frequency of deep convection (i.e., downshear and left-of-shear). Melting-layer displacements are shown that exceed 1 km in altitude compared to melting-layer altitudes in outer rainbands and are complemented by analyses of archived dropsonde data. Dendritic growth and attendant snow aggregation signatures in the inner core are found to occur more often when echo-top altitudes are low (≤10 km MSL), nearer the −15°C isotherm commonly associated with dendritic growth. These signatures, uniquely observed by polarimetric radar, provide greater insight into the physical structure and thermodynamic characteristics of tropical cyclones, which are important for improving rainfall estimation and the representation of tropical cyclones in numerical models. [ABSTRACT FROM AUTHOR]

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 B 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% B 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. [ABSTRACT FROM AUTHOR]

Super resolution involves synthetically increasing the resolution of gridded data beyond their native resolution. Typically, this is done using interpolation schemes, which estimate sub-grid-scale values from neighboring data, and perform the same operation everywhere regardless of the large-scale context, or by requiring a network of radars with overlapping fields of view. Recently, significant progress has been made in single-image super resolution using convolutional neural networks. Conceptually, a neural network may be able to learn relations between large-scale precipitation features and the associated sub-pixel-scale variability and outperform interpolation schemes. Here, we use a deep convolutional neural network to artificially enhance the resolution of NEXRAD PPI scans. The model is trained on 6 months of reflectivity observations from the Langley Hill, Washington, radar (KLGX), and we find that it substantially outperforms common interpolation schemes for 4× and 8× resolution increases based on several objective error and perceptual quality metrics. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

An improved understanding of common differences between tornadic and nontornadic supercells is sought using a large set of observations from the operational NEXRAD WSR-88D polarimetric radar network in the contiguous United States. In particular, data from 478 nontornadic and 294 tornadic supercells during a 7-yr period (2011–17) are used to produce probability-matched composite means of microphysical and kinematic variables. Means, which are centered on echo-top maxima and in a horizontal coordinate system rotated such that storm motion points in the positive x dimension, are created in altitude relative to ground level at times of peak echo-top altitude and peak midlevel rotation for nontornadic supercells and times at and prior to the first tornado in tornadic supercells. Robust differences between supercell types are found, with consistent characteristics at and preceding tornadogenesis in tornadic storms. In particular, the mesocyclone is found to be vertically aligned in tornadic supercells and misaligned in nontornadic supercells. Microphysical differences found include a low-level radar reflectivity hook echo aligned with and ~10 km right of storm center in tornadic supercells and displaced 5–10 km down-motion in nontornadic supercells, a low-to-midlevel differential radar reflectivity dipole that is oriented more parallel to storm motion in tornadic supercells and more perpendicular in nontornadic supercells, and a separation between enhanced differential radar reflectivity and specific differential phase (with unique displacement-relative correlation coefficient reductions) at low levels that is more perpendicular to storm motion in tornadic supercells and more parallel in nontornadic supercells. [ABSTRACT FROM AUTHOR]

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 radarQPE 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. [ABSTRACT FROM AUTHOR]

An assumption related to clouds is one of uncertain factors in precipitation retrievals by the Dual-Frequency Precipitation Radar (DPR) on board the Global Precipitation Measurement (GPM) Core Observatory. While an attenuation due to cloud ice is negligibly small for Ku and Ka bands, attenuation by cloud liquid water is larger in the Ka band and estimating precipitation intensity with high accuracy from Ka-band observations can require developing a method to estimate the attenuation due to cloud liquid water content (CLWC). This paper describes a CLWC database used in the DPR level-2 algorithm for the GPM V06A product. In the algorithm, the CLWC value is assumed using the database with inputs of precipitation-related variables, temperature, and geolocation information. A calculation of the database was made using the 3.5-km-mesh global atmospheric simulation derived from the Nonhydrostatic Icosahedral Atmospheric Model (NICAM) global cloud-system-resolving model. Impacts of current CLWC assumptions for surface precipitation estimates were evaluated by comparisons of precipitation retrieval results between default values and 0 mg m−3 of the CLWC. The impacts were quantified by the normalized mean absolute difference (NMAD) and the NMAD values showed 2.3% for the Ku, 9.9% for the Ka, and 6.5% for the dual-frequency algorithms in global averages, while they were larger in the tropics than in high latitudes. Effects of the precipitation estimates from the CLWC assumption were examined further in terms of retrieval processes affected by the CLWC assumption. This study emphasizes the CLWC assumption provided more effects on the precipitation estimates through estimating path-integrated attenuation due to rain. [ABSTRACT FROM AUTHOR]

High-resolution boundary layer water vapor profile observations are essential for understanding the interplay between shallow convection, cloudiness, and climate in the trade wind atmosphere. As current observation techniques can be limited by low spatial or temporal resolution, the synergistic benefit of combining ground-based microwave radiometer (MWR) and dual-frequency radar is investigated by analyzing the retrieval information content and uncertainty. Synthetic MWR brightness temperatures, as well as simulated dual-wavelength ratios of two radar frequencies are generated for a combination of Ka and W band (KaW), as well as differential absorption radar (DAR) G-band frequencies (167 and 174.8 GHz, G2). The synergy analysis is based on an optimal estimation scheme by varying the configuration of the observation vector. Combining MWR and KaW only marginally increases the retrieval information content. The synergy of MWR with G2 radar is more beneficial due to increasing degrees of freedom (4.5), decreasing retrieval errors, and a more realistic retrieved profile within the cloud layer. The information and profile below and within the cloud is driven by the radar observations, whereas the synergistic benefit is largest above the cloud layer, where information content is enhanced compared to an MWR-only or DAR-only setup. For full synergistic benefits, however, G-band radar sensitivities need to allow full-cloud profiling; in this case, the results suggest that a combined retrieval of MWR and G-band DAR can help close the observational gap of current techniques. [ABSTRACT FROM AUTHOR]

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 reflectivity Z between −30 and 0 dBZ and all values of lidar backscatter β observations. No partitioning of cloud or drizzle is required on the basis of an arbitrary threshold of Z as 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. [ABSTRACT FROM AUTHOR]

Observations of the air vertical velocities (wair) in supercell updrafts are presented, including uncertainty estimates, from radiosonde GPS measurements in two supercells. These in situ observations were collected during the Colorado State University Convective Cloud Outflows and Updrafts Experiment (C³LOUD-Ex) in moderately unstable environments in Colorado and Wyoming. Based on the radiosonde accelerations, instances when the radiosonde balloon likely bursts within the updraft are determined, and adjustments are made to account for the subsequent reduction in radiosonde buoyancy. Before and after these adjustments, the maximum estimated wair values are 36.2 and 49.9ms-1, respectively. Radar data are used to contextualize the in situ observations and suggest that most of the radiosonde observations were located several kilometers away from the most intense vertical motions. Therefore, the radiosonde-based wair values presented likely underestimate the maximum values within these storms due to these sampling biases, as well as the impacts from hydrometeors, which are not accounted for. When possible, radiosonde-based wair values were compared to estimates from dual-Doppler methods and from parcel theory. When the radiosondes observed their highest wair values, dual-Doppler methods generally produced 15–20ms-1 lower wair for the same location, which could be related to the differences in the observing systems’ resolutions. In situ observations within supercell updrafts, which have been limited in recent decades, can be used to improve our understanding and modeling of storm dynamics. This study provides new in situ observations, as well as methods and lessons that could be applied to future field campaigns. [ABSTRACT FROM AUTHOR]

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. SIGNIFICANCE STATEMENT The official nationwide radar network in the United States is capable of detecting some tornadoes, but this capability generally decreases with increasing distance from the radar. In this study, we examined two peculiar radar signatures observed in a tornado-producing thunderstorm in Oklahoma in May 2016. These signatures looked similar to those seen in anticyclonic (i.e., clockwise-rotating in the Northern Hemisphere) tornadoes, but the location within the storm would have been highly unusual for a tornado. Data from a rapid-scan mobile radar closer to the storm revealed that the signatures were not from anticyclonic tornadoes but, rather, may have been quasi-horizontal vortices centered;1 km above ground. This case serves as an example of the benefits expected from a denser radar network of next-generation, rapid-scan weather radars. [ABSTRACT FROM AUTHOR]

Fully polarimetric scanning and vertically pointing Doppler spectral data from the state-of-the-art Stony Brook University Ka-band Scanning Polarimetric Radar (KASPR) are analyzed for a long-duration case of ice pellets over central Long Island in New York from 12 February 2019. Throughout the period of ice pellets, a classic refreezing signature was present, consisting of a secondary enhancement of differential reflectivityZDR beneath the melting layer within a region of decreasing reflectivity factor at horizontal polarization ZH and reduced copolar correlation coefficient rhv. The KASPR radar data allow for evaluation of previously proposed hypotheses to explain the refreezing signature. It is found that, upon entering a layer of locally generated columnar ice crystals and undergoing contact nucleation, smaller raindrops preferentially refreeze into ice pellets prior to the complete freezing of larger drops. Refreezing particles exhibit deformations in shape during freezing, leading to reduced rhv, reduced co-to-cross-polar correlation coefficient rxh, and enhanced linear depolarization ratio, but these shape changes do not explain the ZDR signature. The presence of columnar ice crystals, though apparently crucial for instigating the refreezing process, does not contribute enough backscattered power to affect the ZDR signature, either. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]