scholarly journals Assimilating FY-4A Lightning and Radar Data for Improving Short-Term Forecasts of a High-Impact Convective Event with a Dual-Resolution Hybrid 3DEnVAR Method

2021 ◽  
Vol 13 (16) ◽  
pp. 3090
Author(s):  
Peng Liu ◽  
Yi Yang ◽  
Anwei Lai ◽  
Yunheng Wang ◽  
Alexandre O. Fierro ◽  
...  

A dual-resolution, hybrid, three-dimensional ensemble-variational (3DEnVAR) data assimilation method combining static and ensemble background error covariances is used to assimilate radar data, and pseudo-water vapor observations to improve short-term severe weather forecasts with the Weather Research and Forecast (WRF) model. The higher-resolution deterministic forecast and the lower-resolution ensemble members have 3 and 9 km horizontal resolution, respectively. The water vapor pseudo-observations are derived from the combined use of total lightning data and cloud top height from the Fengyun-4A(FY-4A) geostationary satellite. First, a set of single-analysis experiments are conducted to provide a preliminary performance evaluation of the effectiveness of the hybrid method for assimilating multisource observations; second, a set of cycling analysis experiments are used to evaluate the forecast performance in convective-scale high-frequency analysis; finally, different hybrid coefficients are tested in both the single and cycling experiments. The single-analysis results show that the combined assimilation of radar data and water vapor pseudo-observations derived from the lightning data is able to generate reasonable vertical velocity, water vapor and hydrometeor adjustments, which help to trigger convection earlier in the forecast/analysis and reduce the spin-up time. The dual-resolution hybrid 3DEnVAR method is able to adjust the wind fields and hydrometeor variables with the assimilation of lightning data, which helps maintain the triggered convection longer and partially suppress spurious cells in the forecast compared with the three-dimensional variational (3DVAR) method. A cycling analysis that introduced a large number of observations with more frequent small adjustments is able to better resolve the observed convective events than a single-analysis approach. Different hybrid coefficients can affect the forecast results, either in the single deterministic or cycling analysis experiments. Overall, we found that a static coefficient of 0.4 and an ensemble coefficient of 0.6 yields the best forecast skill for this event.

2016 ◽  
Vol 144 (11) ◽  
pp. 4373-4393 ◽  
Author(s):  
Alexandre O. Fierro ◽  
Jidong Gao ◽  
Conrad L. Ziegler ◽  
Kristin M. Calhoun ◽  
Edward R. Mansell ◽  
...  

Abstract This work evaluates the performance of the assimilation of total lightning data within a three-dimensional variational (3DVAR) framework for the analysis and short-term forecast of the 24 May 2011 tornado outbreak using the Weather Research and Forecasting (WRF) Model at convection-allowing scales. Between the lifted condensation level and a fixed upper height, pseudo-observations for water vapor mass first are created based on either the flash extent densities derived from Oklahoma Lightning Mapping Array data or the lightning source densities derived from the Earth Networks pulse data, and then assimilated by the 3DVAR system. Assimilation of radar data with 3DVAR and a cloud analysis algorithm (RAD) also are performed as a baseline for comparison and in tandem with lightning to evaluate the added value of this lightning data assimilation (LDA) method. Given a scenario wherein the control experiment without radar or lightning data assimilation fails to accurately initiate and forecast the observed storms, the LDA and RAD yield comparable short-term forecast improvements. The RAD alone produces storms of similar strength to the observations during the first 30 min of forecast more rapidly than the LDA alone; however, the LDA is able to better depict individual supercellular features at 1-h forecast. When both the lightning and radar data are assimilated, the 30-min forecast showed noteworthy improvements over RAD in terms of the model’s ability to better resolve individual supercell structures and still maintained a 1-h forecast similar to that from the LDA. The results chiefly illustrate the potential value of assimilating total lightning data along with radar data.


2019 ◽  
Vol 12 (9) ◽  
pp. 4031-4051 ◽  
Author(s):  
Shizhang Wang ◽  
Zhiquan Liu

Abstract. A reflectivity forward operator and its associated tangent linear and adjoint operators (together named RadarVar) were developed for variational data assimilation (DA). RadarVar can analyze both rainwater and ice-phase species (snow and graupel) by directly assimilating radar reflectivity observations. The results of three-dimensional variational (3D-Var) DA experiments with a 3 km grid mesh setting of the Weather Research and Forecasting (WRF) model showed that RadarVar was effective at producing an analysis of reflectivity pattern and intensity similar to the observed data. Two to three outer loops with 50–100 iterations in each loop were needed to obtain a converged 3-D analysis of reflectivity, rainwater, snow, and graupel, including the melting layers with mixed-phase hydrometeors. It is shown that the deficiencies in the analysis using this operator, caused by the poor quality of the background fields and the use of the static background error covariance, can be partially resolved by using radar-retrieved hydrometeors in a preprocessing step and tuning the spatial correlation length scales of the background errors. The direct radar reflectivity assimilation using RadarVar also improved the short-term (2–5 h) precipitation forecasts compared to those of the experiment without DA.


2021 ◽  
Author(s):  
Alexandre Fierro ◽  
Junjun Hu ◽  
Yunheng Wang ◽  
Jidong Gao ◽  
Edward Mansell

<p>The GLM instruments aboard the GOES-16 and 17 satellites provides nearly uniform spatiotemporal coverage of total lightning over the Americas and adjacent vast oceanic regions of the western hemisphere. This work summarizes recent efforts from our group at CIMMS/NSSL geared towards the evaluation of the potential added value of assimilating GLM-observed total lightning data on short-term, convection-allowing scale (dx = 2-3 km) forecasts for higher impact weather events. Results using data assimilation (DA) approaches ranging from single deterministic three-dimensional variational (3DVAR) methods applied in real time to experimental ensemble-based VAR hybrid methods (3DEnVAR) will be highlighted. <br>The lightning data assimilation (DA) scheme in these frameworks follow the same core philosophy wherein background water vapor mass mixing ratio is adjusted (increased) locally at or around observed lightning locations, either throughout the entire atmospheric column or within a fixed, confined layer above the lifted condensation level. Toward a more systematic assimilation of real GLM data, emphasis will be directed toward: (i) sensitivity tests with deterministic 3DVAR experiments aimed at evaluating the impact of the horizontal decorrelation length scale, DA cycling frequency as well the length of the accumulation window for the lightning data, (ii) aggregate statistics from real time CONUS-scale experiments over the Spring 2020 and (iii) preliminary results employing ensemble of 3DEnVARs with hybrid (static + flow dependent) background error covariances. <br>Aggregate statistical results from all deterministic 3DVAR exercises in (i) and (ii) revealed that the assimilation of either radar (radial wind and reflectivity factor) or total lightning (GLM) resulted in overall notably more skillful, shorter term (0-3 h) forecast of composite reflectivity fields, accumulated rainfall, as well as individual storm tracks – with optimal skill obtained when both radar and lightning data were assimilated. In (iii) forecast impacts related to the following will be summarized: (1) the respective weights assigned to the flow-dependent component and static components of the background error covariances, (2) the inclusion of three time-level sampling for each member during each cycle and (3) the usage of Gaussian noise coupled with a fixed 3 to 12 h spin-up period prior to the beginning of the cycled 3DVAR.</p>


2020 ◽  
Vol 148 (3) ◽  
pp. 1005-1028 ◽  
Author(s):  
Junjun Hu ◽  
Alexandre O. Fierro ◽  
Yunheng Wang ◽  
Jidong Gao ◽  
Edward R. Mansell

Abstract The recent successful deployment of the Geostationary Lightning Mapper (GLM) on board the Geostationary Operational Environmental Satellite R series (GOES-16/17) provides nearly uniform spatiotemporal measurements of total lightning (intracloud plus cloud to ground) over the Americas and adjacent vast oceanic regions. This study evaluates the potential value of assimilating GLM-derived water vapor mixing ratio on short-term (≤6 h), cloud-scale (dx = 1.5 km) forecasts of five severe weather events over the Great Plains of the United States using a three-dimensional variational (3DVAR) data assimilation (DA) system. Toward a more systematic assimilation of real GLM data, this study conducted sensitivity tests aimed at evaluating the impact of the horizontal decorrelation length scale, DA cycling frequency, and the time window size for accumulating GLM lightning observations prior to the DA. Forecast statistics aggregated over all five cases suggested that an optimal forecast performance is obtained when lightning measurements are accumulated over a 10-min interval and GLM-derived water vapor mixing ratio values are assimilated every 15 min with a horizontal decorrelation length scale of 3 km. This suggested configuration for the GLM DA together with companion experiments (i) not assimilating any data, (ii) assimilating radar data only, and (iii) assimilating both GLM and radar data were evaluated for the same five cases. Overall, GLM data have shown potential to help improve the short-term (<3 h) forecast skill of composite reflectivity fields and individual storm tracks. While this result also held for accumulated rainfall, longer-term (≥3 h) forecasts were generally characterized by noteworthy wet biases.


2014 ◽  
Vol 142 (8) ◽  
pp. 2687-2708 ◽  
Author(s):  
Tammy M. Weckwerth ◽  
Lindsay J. Bennett ◽  
L. Jay Miller ◽  
Joël Van Baelen ◽  
Paolo Di Girolamo ◽  
...  

Abstract A case study of orographic convection initiation (CI) that occurred along the eastern slopes of the Vosges Mountains in France on 6 August 2007 during the Convective and Orographically-Induced Precipitation Study (COPS) is presented. Global positioning system (GPS) receivers and two Doppler on Wheels (DOW) mobile radars sampled the preconvective and storm environments and were respectively used to retrieve three-dimensional tomographic water vapor and wind fields. These retrieved data were supplemented with temperature, moisture, and winds from radiosondes from a site in the eastern Rhine Valley. High-resolution numerical simulations with the Weather Research and Forecasting (WRF) Model were used to further investigate the physical processes leading to convective precipitation. This unique, time-varying combination of derived water vapor and winds from observations illustrated an increase in low-level moisture and convergence between upslope easterlies and downslope westerlies along the eastern slope of the Vosges Mountains. Uplift associated with these shallow, colliding boundary layer flows eventually led to the initiation of moist convection. WRF reproduced many features of the observed complicated flow, such as cyclonic (anticyclonic) flow around the southern (northern) end of the Vosges Mountains and the east-side convergent flow below the ridgeline. The WRF simulations also illustrated spatial and temporal variability in buoyancy and the removal of the lids prior to convective development. The timing and location of CI from the WRF simulations was surprisingly close to that observed.


2019 ◽  
Author(s):  
Shizhang Wang ◽  
Zhiquan Liu

Abstract. A reflectivity forward operator and its associated tangent linear and adjoint operators (together named RadZIceVar) were developed for variational data assimilation (DA). RadZIceVar can analyze both rainwater and ice-phase species (snow and graupel) by directly assimilating radar reflectivity observations. The results of three-dimensional variational (3DVAR) DA experiments with a 3 km grid mesh setting of the Weather Research and Forecasting (WRF) model showed that RadZIceVar was effective at producing an analysis of reflectivity pattern and intensity similar to the observed data. Two to three outer loops with 50–100 iterations in each loop were needed to obtain a converged 3D analysis of rainwater, snow, and graupel, including the melting layers with mixed-phase hydrometeors. The deficiencies in the analysis using this operator could be caused by the poor quality of the background fields and the use of the static background error covariance, and these issues can be partially resolved by using radar-retrieved hydrometeors in a preprocessing step and tuning the spatial correlation length scales of the background errors. The direct radar reflectivity assimilation using RadZIceVar also improved the short-term (2 h–5 h) precipitation forecasts compared to those of the experiment without DA.


2018 ◽  
Vol 33 (1) ◽  
pp. 71-88 ◽  
Author(s):  
Shibo Gao ◽  
Juanzhen Sun ◽  
Jinzhong Min ◽  
Ying Zhang ◽  
Zhuming Ying

Abstract Radar reflectivity observations contain valuable information on precipitation and have been assimilated into numerical weather prediction models for improved microphysics initialization. However, low-reflectivity (or so-called no rain) echoes have often been ignored or not effectively used in radar data assimilation schemes. In this paper, a scheme to assimilate no-rain radar observations is described within the framework of the Weather Research and Forecasting Model’s three-dimensional variational data assimilation (3DVar) system, and its impact on precipitation forecasts is demonstrated. The key feature of the scheme is a neighborhood-based approach to adjusting water vapor when a grid point is deemed as no rain. The performance of the scheme is first examined using a severe convective case in the Front Range of the Colorado Rocky Mountains and then verified by running the 3DVar system in the same region, with and without the no-rain assimilation scheme for 68 days and 3-hourly rapid update cycles. It is shown that the no-rain data assimilation method reduces the bias and false alarm ratio of precipitation over its counterpart without that assimilation. The no-rain assimilation also improved humidity, temperature, and wind fields, with the largest error reduction in the water vapor field, both near the surface and at upper levels. It is also shown that the advantage of the scheme is in its ability to conserve total water content in cycled radar data assimilation, which cannot be achieved by assimilating only precipitation echoes.


2008 ◽  
Vol 136 (3) ◽  
pp. 945-963 ◽  
Author(s):  
Jidong Gao ◽  
Ming Xue

Abstract A new efficient dual-resolution (DR) data assimilation algorithm is developed based on the ensemble Kalman filter (EnKF) method and tested using simulated radar radial velocity data for a supercell storm. Radar observations are assimilated on both high-resolution and lower-resolution grids using the EnKF algorithm with flow-dependent background error covariances estimated from the lower-resolution ensemble. It is shown that the flow-dependent and dynamically evolved background error covariances thus estimated are effective in producing quality analyses on the high-resolution grid. The DR method has the advantage of being able to significantly reduce the computational cost of the EnKF analysis. In the system, the lower-resolution ensemble provides the flow-dependent background error covariance, while the single-high-resolution forecast and analysis provides the benefit of higher resolution, which is important for resolving the internal structures of thunderstorms. The relative smoothness of the covariance obtained from the lower 4-km-resolution ensemble does not appear to significantly degrade the quality of analysis. This is because the cross covariance among different variables is of first-order importance for “retrieving” unobserved variables from the radar radial velocity data. For the DR analysis, an ensemble size of 40 appears to be a reasonable choice with the use of a 4-km horizontal resolution in the ensemble and a 1-km resolution in the high-resolution analysis. Several sensitivity tests show that the DR EnKF system is quite robust to different observation errors. A 4-km thinned data resolution is a compromise that is acceptable under the constraint of real-time applications. A data density of 8 km leads to a significant degradation in the analysis.


2014 ◽  
Vol 142 (10) ◽  
pp. 3756-3780 ◽  
Author(s):  
Yujie Pan ◽  
Kefeng Zhu ◽  
Ming Xue ◽  
Xuguang Wang ◽  
Ming Hu ◽  
...  

Abstract A coupled ensemble square root filter–three-dimensional ensemble-variational hybrid (EnSRF–En3DVar) data assimilation (DA) system is developed for the operational Rapid Refresh (RAP) forecasting system. The En3DVar hybrid system employs the extended control variable method, and is built on the NCEP operational gridpoint statistical interpolation (GSI) three-dimensional variational data assimilation (3DVar) framework. It is coupled with an EnSRF system for RAP, which provides ensemble perturbations. Recursive filters (RF) are used to localize ensemble covariance in both horizontal and vertical within the En3DVar. The coupled En3DVar hybrid system is evaluated with 3-h cycles over a 9-day period with active convection. All conventional observations used by operational RAP are included. The En3DVar hybrid system is run at ⅓ of the operational RAP horizontal resolution or about 40-km grid spacing, and its performance is compared to parallel GSI 3DVar and EnSRF runs using the same datasets and resolution. Short-term forecasts initialized from the 3-hourly analyses are verified against sounding and surface observations. When using equally weighted static and ensemble background error covariances and 40 ensemble members, the En3DVar hybrid system outperforms the corresponding GSI 3DVar and EnSRF. When the recursive filter coefficients are tuned to achieve a similar height-dependent localization as in the EnSRF, the En3DVar results using pure ensemble covariance are close to EnSRF. Two-way coupling between EnSRF and En3DVar did not produce noticeable improvement over one-way coupling. Downscaled precipitation forecast skill on the 13-km RAP grid from the En3DVar hybrid is better than those from GSI 3DVar analyses.


2010 ◽  
Vol 3 (5) ◽  
pp. 4459-4495 ◽  
Author(s):  
C. López Carrillo ◽  
D. J. Raymond

Abstract. In this work, we describe an efficient approach for wind retrieval from dual Doppler radar data. The approach produces a gridded field that not only satisfies the observations, but also satisfies the anelastic mass continuity equation. The method is based on the so-called three-dimensional variational approach to the retrieval of wind fields from radar data. The novelty consists in separating the task into steps that reduce the amount of data processed by the global minimization algorithm, while keeping the most relevant information from the radar observations. The method is flexible enough to incorporate observations from several radars, accommodate complex sampling geometries, and readily include dropsonde or sounding observations in the analysis. We demonstrate the usefulness of our method by analyzing a real case with data collected during the TPARC/TCS-08 field campaign.


Sign in / Sign up

Export Citation Format

Share Document