Assimilation of Doppler Radar Observations with a Regional 3DVAR System: Impact of Doppler Velocities on Forecasts of a Heavy Rainfall Case

2005 ◽  
Vol 44 (6) ◽  
pp. 768-788 ◽  
Author(s):  
Qingnong Xiao ◽  
Ying-Hwa Kuo ◽  
Juanzhen Sun ◽  
Wen-Chau Lee ◽  
Eunha Lim ◽  
...  

Abstract In this paper, the impact of Doppler radar radial velocity on the prediction of a heavy rainfall event is examined. The three-dimensional variational data assimilation (3DVAR) system for use with the fifth-generation Pennsylvania State University–NCAR Mesoscale Model (MM5) is further developed to enable the assimilation of radial velocity observations. Doppler velocities from the Korean Jindo radar are assimilated into MM5 using the 3DVAR system for a heavy rainfall case that occurred on 10 June 2002. The results show that the assimilation of Doppler velocities has a positive impact on the short-range prediction of heavy rainfall. The dynamic balance between atmospheric wind and thermodynamic fields, based on the Richardson equation, is introduced to the 3DVAR system. Vertical velocity (w) increments are included in the 3DVAR system to enable the assimilation of the vertical velocity component of the Doppler radial velocity observation. The forecast of the hydrometeor variables of cloud water (qc) and rainwater (qr) is used in the 3DVAR background fields. The observation operator for Doppler radial velocity is developed and implemented within the 3DVAR system. A series of experiments, assimilating the Korean Jindo radar data for the 10 June 2002 heavy rainfall case, indicates that the scheme for Doppler velocity assimilation is stable and robust in a cycling mode making use of high-frequency radar data. The 3DVAR with assimilation of Doppler radial velocities is shown to improve the prediction of the rainband movement and intensity change. As a result, an improved skill for the short-range heavy rainfall forecast is obtained. The forecasts of other quantities, for example, winds, are also improved. Continuous assimilation with 3-h update cycles is important in producing an improved heavy rainfall forecast. Assimilation of Doppler radar radial velocities using the 3DVAR background fields from a cycling procedure produces skillful rainfall forecasts when verified against observations.

2010 ◽  
Vol 138 (4) ◽  
pp. 1250-1272 ◽  
Author(s):  
David J. Stensrud ◽  
Jidong Gao

Abstract The assimilation of operational Doppler radar observations into convection-resolving numerical weather prediction models for very short-range forecasting represents a significant scientific and technological challenge. Numerical experiments over the past few years indicate that convective-scale forecasts are sensitive to the details of the data assimilation methodology, the quality of the radar data, the parameterized microphysics, and the storm environment. In this study, the importance of horizontal environmental variability to very short-range (0–1 h) convective-scale ensemble forecasts initialized using Doppler radar observations is investigated for the 4–5 May 2007 Greensburg, Kansas, tornadic thunderstorm event. Radar observations of reflectivity and radial velocity from the operational Doppler radar network at 0230 UTC 5 May 2007, during the time of the first large tornado, are assimilated into each ensemble member using a three-dimensional variational data assimilation system (3DVAR) developed at the Center for Analysis and Prediction of Storms (CAPS). Very short-range forecasts are made using the nonhydrostatic Advanced Regional Prediction System (ARPS) model from each ensemble member and the results are compared with the observations. Explicit three-dimensional environmental variability information is provided to the convective-scale ensemble using analyses from a 30-km mesoscale ensemble data assimilation system. Comparisons between convective-scale ensembles with initial conditions produced by 3DVAR using 1) background fields that are horizontally homogeneous but vertically inhomogeneous (i.e., have different vertical environmental profiles) and 2) background fields that are horizontally and vertically inhomogeneous are undertaken. Results show that the ensemble with horizontally and vertically inhomogeneous background fields provides improved predictions of thunderstorm structure, mesocyclone track, and low-level circulation track than the ensemble with horizontally homogeneous background fields. This suggests that knowledge of horizontal environmental variability is important to successful convective-scale ensemble predictions and needs to be included in real-data experiments.


2005 ◽  
Vol 62 (1) ◽  
pp. 220-230 ◽  
Author(s):  
Robert Nissen ◽  
Roland List ◽  
David Hudak ◽  
Greg M. McFarquhar ◽  
R. Paul Lawson ◽  
...  

Abstract For nonconvective, steady light rain with rain rates <5 mm h−1 the mean Doppler velocity of raindrop spectra was found to be constant below the melting band, when the drop-free fall speed was adjusted for pressure. The Doppler radar–weighted raindrop diameters varied from case to case from 1.5 to 2.5 mm while rain rates changed from 1.2 to 2.9 mm h−1. Significant changes of advected velocity moments were observed over periods of 4 min. These findings were corroborated by three independent systems: a Doppler radar for establishing vertical air speed and mean terminal drop speeds [using extended Velocity Azimuth Display (EVAD) analyses], a Joss–Waldvogel disdrometer at the ground, and a Particle Measuring System (PMS) 2-DP probe flown on an aircraft. These measurements were supported by data from upper-air soundings. The reason why inferred raindrop spectra do not change with height is the negligible interaction rate between raindrops at low rain rates. At low rain rates, numerical box models of drop collisions strongly support this interpretation. It was found that increasing characteristic drop diameters are correlated with increasing rain rates.


Author(s):  
VINCENT T. WOOD ◽  
ROBERT P. DAVIES-JONES ◽  
ALAN SHAPIRO

AbstractSingle-Doppler radar data are often missing in important regions of a severe storm due to low return power, low signal-to-noise ratio, ground clutter associated with normal and anomalous propagation, and missing radials associated with partial or total beam blockage. Missing data impact the ability of WSR-88D algorithms to detect severe weather. To aid the algorithms, we develop a variational technique that fills in Doppler velocity data voids smoothly by minimizing Doppler velocity gradients while not modifying good data. This method provides estimates of the analysed variable in data voids without creating extrema.Actual single-Doppler radar data of four tornadoes are used to demonstrate the variational algorithm. In two cases, data are missing in the original data, and in the other two, data are voided artificially. The filled-in data match the voided data well in smoothly varying Doppler velocity fields. Near singularities such as tornadic vortex signatures, the match is poor as anticipated. The algorithm does not create any velocity peaks in the former data voids, thus preventing false triggering of tornado warnings. Doppler circulation is used herein as a far-field tornado detection and advance-warning parameter. In almost all cases, the measured circulation is quite insensitive to the data that have been voided and then filled. The tornado threat is still apparent.


2017 ◽  
Vol 145 (10) ◽  
pp. 4187-4203 ◽  
Author(s):  
Feng Chen ◽  
Xudong Liang ◽  
Hao Ma

An improved Doppler radar radial velocity assimilation observation operator is proposed based on the integrating velocity–azimuth process (IVAP) method. This improved operator can ingest both radial wind and its spatial distribution characteristics to deduce the two components of the mean wind within a given area. With this operator, the system can be used to assimilate information from tangential wind and radial wind. On the other hand, because the improved observation operator is defined within a given area, which can be uniformly chosen in both the observation and analysis coordinate systems, it has a thinning function. The traditional observation operator and the improved observation operator, along with their corresponding data processing modules, were implemented in the community Gridpoint Statistical Interpolation analysis system (GSI) to demonstrate the superiority of the improved operator. The results of single analysis unit experiments revealed that the two operators are comparable when the analysis unit is small. When the analysis unit becomes larger, the analysis results of the improved operator are better than those of the traditional operator because the former can ingest more wind information than the latter. The results of a typhoon case study indicated that both operators effectively ingested radial wind information and produced more reasonable typhoon structures than those in the background fields. The tangential velocity relative to the radar was retrieved by the improved operator through ingesting tangential wind information from the spatial distribution characteristics of radial wind. Because of the improved vortex intensity and structure, obvious improvements were seen in both track and intensity predictions when the improved operator was used.


2012 ◽  
Vol 140 (5) ◽  
pp. 1603-1619 ◽  
Author(s):  
Yu-Chieng Liou ◽  
Shao-Fan Chang ◽  
Juanzhen Sun

This study develops an extension of a variational-based multiple-Doppler radar synthesis method to construct the three-dimensional wind field over complex topography. The immersed boundary method (IBM) is implemented to take into account the influence imposed by a nonflat surface. The IBM has the merit of providing realistic topographic forcing without the need to change the Cartesian grid configuration into a terrain-following coordinate system. Both Dirichlet and Neumann boundary conditions for the wind fields can be incorporated. The wind fields above the terrain are obtained by variationally adjusting the solutions to satisfy a series of weak constraints, which include the multiple-radar radial velocity observations, anelastic continuity equation, vertical vorticity equation, background wind, and spatial smoothness terms. Experiments using model-simulated data reveal that the flow structures over complex orography can be successfully retrieved using radial velocity measurements from multiple Doppler radars. The primary advantages of the original synthesis method are still maintained, that is, the winds along and near the radar baseline are well retrieved, and the resulting three-dimensional flow fields can be used directly for vorticity budget diagnosis. If compared with the traditional wind synthesis algorithm, this method is able to merge data from different sources, and utilize data from any number of radars. This provides more flexibility in designing various scanning strategies, so that the atmosphere may be probed more efficiently using a multiple-radar network. This method is also tested using the radar data collected during the Southwest Monsoon Experiment (SoWMEX), which was conducted in Taiwan from May to June 2008 with reasonable results being obtained.


2020 ◽  
Author(s):  
Hongli Li ◽  
Yang Hu ◽  
Zhimin Zhou

<p>During the Meiyu period, floods are prone to occur in the middle and lower reaches of the Yangtze River due to the highly concentrated and heavy rainfall, which caused huge life and economic losses. Based on numerical simulation by assimilating Doppler radar, radiosonde, and surface meteorological observations, the evolution mechanism for the initiation, development and decaying of a Meiyu frontal rainstorm that occurred from 4th to 5th July 2014 is analyzed in this study. Results show that the numerical experiment can well reproduce the temporal variability of heavy precipitation and successfully simulate accumulative precipitation and its evolution over the key rainstorm area. The simulated “rainbelt training” is consistent with observed “echo training” on both spatial structure and temporal evolution. The convective cells in the mesoscale convective belt propagated from southwest to northeast across the key rainstorm area, leading to large accumulative precipitation and rainstorm in this area. There existed convective instability in lower levels above the key rainstorm area, while strong ascending motion developed during period of heavy rainfall. Combined with abundant water vapor supply, the above condition was favorable for the formation and development of heavy rainfall. The Low level jet (LLJ) provided sufficient energy for the rainstorm system, and the low-level convergence intensified, which was an important reason for the maintenance of precipitation system and its eventual intensification to rainstorm. At its mature stage, the rainstorm system demonstrated vertically tilted structure with strong ascending motion in the key rainstorm area, which was favorable for the occurrence of heavy rainfall. In the decaying stage, unstable energy decreased, and the rainstorm no longer had sufficient energy to sustain. The rapid weakening of LLJ resulted in smaller energy supply to the convective system, and the stratification tended to be stable in the middle and lower levels. The ascending motion weakened correspondingly, which made it hard for the convective system to maintain.</p>


2010 ◽  
Vol 27 (7) ◽  
pp. 1140-1152 ◽  
Author(s):  
Eunha Lim ◽  
Juanzhen Sun

Abstract A Doppler velocity dealiasing algorithm is developed within the storm-scale four-dimensional radar data assimilation system known as the Variational Doppler Radar Analysis System (VDRAS). The innovative aspect of the algorithm is that it dealiases Doppler velocity at each grid point independently by using three-dimensional wind fields obtained either from an objective analysis using conventional observations and mesoscale model output or from a rapidly updated analysis of VDRAS that assimilates radar data. This algorithm consists of three steps: preserving horizontal shear, global dealiasing using reference wind from the objective analysis or the VDRAS analysis, and local dealiasing. It is automated and intended to be used operationally for radar data assimilation using numerical weather prediction models. The algorithm was tested with 384 volumes of radar data observed from the Next Generation Weather Radar (NEXRAD) for a severe thunderstorm that occurred during 15 June 2002. It showed that the algorithm was effective in dealiasing large areas of aliased velocities when the wind from the objective analysis was used as the reference and that more accurate dealiasing was achieved by using the continuously cycled VDRAS analysis.


2010 ◽  
Vol 2010 ◽  
pp. 1-14 ◽  
Author(s):  
Qin Xu ◽  
Li Wei ◽  
Wei Gu ◽  
Jiandong Gong ◽  
Qingyun Zhao

A 3.5-dimensional variational method is developed for Doppler radar data assimilation. In this method, incremental analyses are performed in three steps to update the model state upon the background state provided by the model prediction. First, radar radial-velocity observations from three consecutive volume scans are analyzed on the model grid. The analyzed radial-velocity fields are then used in step 2 to produce incremental analyses for the vector velocity fields at two time levels between the three volume scans. The analyzed vector velocity fields are used in step 3 to produce incremental analyses for the thermodynamic fields at the central time level accompanied by the adjustments in water vapor and hydrometeor mixing ratios based on radar reflectivity observations. The finite element B-spline representations and recursive filter are used to reduce the dimension of the analysis space and enhance the computational efficiency. The method is applied to a squall line case observed by the phased-array radar with rapid volume scans at the National Weather Radar Testbed and is shown to be effective in assimilating the phased-array radar observations and improve the prediction of the subsequent evolution of the squall line.


2017 ◽  
Vol 56 (12) ◽  
pp. 3263-3283 ◽  
Author(s):  
J. Rémillard ◽  
A. M. Fridlind ◽  
A. S. Ackerman ◽  
G. Tselioudis ◽  
P. Kollias ◽  
...  

AbstractA case study of persistent stratocumulus over the Azores is simulated using two independent large-eddy simulation (LES) models with bin microphysics, and forward-simulated cloud radar Doppler moments and spectra are compared with observations. Neither model is able to reproduce the monotonic increase of downward mean Doppler velocity with increasing reflectivity that is observed under a variety of conditions, but for differing reasons. To a varying degree, both models also exhibit a tendency to produce too many of the largest droplets, leading to excessive skewness in Doppler velocity distributions, especially below cloud base. Excessive skewness appears to be associated with an insufficiently sharp reduction in droplet number concentration at diameters larger than ~200 μm, where a pronounced shoulder is found for in situ observations and a sharp reduction in reflectivity size distribution is associated with relatively narrow observed Doppler spectra. Effectively using LES with bin microphysics to study drizzle formation and evolution in cloud Doppler radar data evidently requires reducing numerical diffusivity in the treatment of the stochastic collection equation; if that is accomplished sufficiently to reproduce typical spectra, progress toward understanding drizzle processes is likely.


2007 ◽  
Vol 135 (10) ◽  
pp. 3381-3404 ◽  
Author(s):  
Qingnong Xiao ◽  
Juanzhen Sun

Abstract The impact of multiple–Doppler radar data assimilation on quantitative precipitation forecasting (QPF) is examined in this study. The newly developed Weather Research and Forecasting (WRF) model Advanced Research WRF (ARW) and its three-dimensional variational data assimilation system (WRF 3DVAR) are used. In this study, multiple–Doppler radar data assimilation is applied in WRF 3DVAR cycling mode to initialize a squall-line convective system on 13 June 2002 during the International H2O Project (IHOP_2002) and the ARW QPF skills are evaluated for the case. Numerical experiments demonstrate that WRF 3DVAR can successfully assimilate Doppler radial velocity and reflectivity from multiple radar sites and extract useful information from the radar data to initiate the squall-line convective system. Assimilation of both radial velocity and reflectivity results in sound analyses that show adjustments in both the dynamical and thermodynamical fields that are consistent with the WRF 3DVAR balance constraint and background error correlation. The cycling of the Doppler radar data from the 12 radar sites at 2100 UTC 12 June and 0000 UTC 13 June produces a more detailed mesoscale structure of the squall-line convection in the model initial conditions at 0000 UTC 13 June. Evaluations of the ARW QPF skills with initialization via Doppler radar data assimilation demonstrate that the more radar data in the temporal and spatial dimensions are assimilated, the more positive is the impact on the QPF skill. Assimilation of both radial velocity and reflectivity has more positive impact on the QPF skill than does assimilation of either radial velocity or reflectivity only. The improvement of the QPF skill with multiple-radar data assimilation is more clearly observed in heavy rainfall than in light rainfall. In addition to the improvement of the QPF skill, the simulated structure of the squall line is also enhanced by the multiple–Doppler radar data assimilation in the WRF 3DVAR cycling experiment. The vertical airflow pattern shows typical characteristics of squall-line convection. The cold pool and its related squall-line convection triggering process are better initiated in the WRF 3DVAR analysis and simulated in the ARW forecast when multiple–Doppler radar data are assimilated.


Sign in / Sign up

Export Citation Format

Share Document