scholarly journals Short-Term Forecasting of a Midlatitude Convective Storm by the Assimilation of Single–Doppler Radar Observations

2009 ◽  
Vol 137 (12) ◽  
pp. 4115-4135 ◽  
Author(s):  
Kao-Shen Chung ◽  
Isztar Zawadzki ◽  
M. K. Yau ◽  
Luc Fillion

Abstract The McGill University radar data assimilation system is used to initialize a convective storm at high resolution (1 km) from single–Doppler radar observations. In this study, the background term in the assimilation system is improved. Specifically, by assuming the correlation of the errors of the control variables to be isotropic and homogeneous, the background error covariance matrix is modeled by a recursive filter. In addition, a 3-h-prior high-resolution model forecast is used as the background field. The analysis fields from the assimilation system successfully trigger the convective storms in the radar-observed regions from a single assimilation window. Without data assimilation, the modeled storms did not occur at the right time and place. To account for the rapid evolution of the convective storms and to correct the forecast errors with time, a cycling process is applied for a very short-term forecast. It is found that the first assimilation window can maintain the prediction of the storms for less than 1 h. The cycling process helps to maintain the intensity of the storm cells for a longer period of time. However, a comparison of radar observations with the 90-min simulation indicates an error in the position of the convective cells. The error of the radial component of the wind field between the observation and the simulation is larger at the upper levels. A wavelet analysis between the observation and simulated reflectivities indicates that the forecast is able to adequately predict the convective scale (∼10–20 km) during the first 20 min, whereas the simulation has more predictability at the longer scale (>30 km) beyond 20 min.

2016 ◽  
Vol 144 (9) ◽  
pp. 3159-3180 ◽  
Author(s):  
Rebecca M. Cintineo ◽  
Jason A. Otkin ◽  
Thomas A. Jones ◽  
Steven Koch ◽  
David J. Stensrud

This study uses an observing system simulation experiment to explore the impact of assimilating GOES-R Advanced Baseline Imager (ABI) 6.95-μm brightness temperatures and Weather Surveillance Radar-1988 Doppler (WSR-88D) reflectivity and radial velocity observations in an ensemble data assimilation system. A high-resolution truth simulation was used to create synthetic radar and satellite observations of a severe weather event that occurred across the U.S. central plains on 4–5 June 2005. The experiment employs the Weather Research and Forecasting Model at 4-km horizontal grid spacing and the ensemble adjustment Kalman filter algorithm in the Data Assimilation Research Testbed system. The ability of GOES-R ABI brightness temperatures to improve the analysis and forecast accuracy when assimilated separately or simultaneously with Doppler radar reflectivity and radial velocity observations was assessed, along with the use of bias correction and different covariance localization radii for the brightness temperatures. Results show that the radar observations accurately capture the structure of a portion of the storm complex by the end of the assimilation period, but that more of the storms and atmospheric features are reproduced and the accuracy of the ensuing forecast improved when the brightness temperatures are also assimilated.


2014 ◽  
Vol 142 (6) ◽  
pp. 2118-2138 ◽  
Author(s):  
Weiguang Chang ◽  
Kao-Shen Chung ◽  
Luc Fillion ◽  
Seung-Jong Baek

Abstract An 80-member high-resolution ensemble Kalman filter (HREnKF) is implemented for assimilating radar observations with the Canadian Meteorological Center’s (CMC’s) Global Environmental Multiscale Limited-Area Model (GEM-LAM). This system covers the Montréal, Canada, region and assimilates radar data from the McGill Radar Observatory with 4-km data thinning. The GEM-LAM operates in fully nonhydrostatic mode with 58 hybrid vertical levels and 1-km horizontal grid spacing. As a first step toward full radar data assimilation, only radial velocities are directly assimilated in this study. The HREnKF is applied on three 2011 summer cases having different precipitation structures: squall-line structure, isolated small-scale structures, and widespread stratiform precipitation. The short-term (<2 h) accuracy of the HREnKF analyses and forecasts is examined. In HREnKF, the ensemble spread is sufficient to cover the estimated error from innovations and lead to filter convergence. It results in part from a realistic initiation of HREnKF data assimilation cycle by using a Canadian regional EnKF system (itself coupled to a global EnKF) working at meso- and synoptic scales. The filter convergence is confirmed by the HREnKF background fields gradually approaching to radar observations as the assimilation cycling proceeds. At each analysis step, it is clearly shown that unobserved variables are significantly modified through HREnKF cross correlation of errors from the ensemble. Radar reflectivity observations are used to verify the improvements in analyses and short-term forecasts achievable by assimilating only radial velocities. Further developments of the analysis system are discussed.


Atmosphere ◽  
2020 ◽  
Vol 11 (9) ◽  
pp. 1007
Author(s):  
Shixuan Zhang ◽  
Zhaoxia Pu

The feasibility of a hurricane initialization framework based on the Gridpoint Statistical Interpolation (GSI)-based four-dimensional ensemble-variational (GSI-4DEnVar) hybrid data assimilation system for the Hurricane Weather Research and Forecasting model (HWRF) model is evaluated in this study. The system considers the temporal evolution of error covariances via the use of four-dimensional ensemble perturbations that are provided by high-resolution, self-consistent HWRF ensemble forecasts. It is different from the configuration of the GSI-based three-dimensional ensemble-variational (GSI-3DEnVar) hybrid data assimilation system, similar to that used in the operational HWRF, which employs background error covariances provided by coarser-resolution global ensembles from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) ensemble Kalman filtering data assimilation system. In addition, our proposed initialization framework discards the empirical intensity correction in the vortex initialization package that is employed by the GSI-3DEnVar initialization framework in operational HWRF. Data assimilation and numerical simulation experiments for Hurricanes Joaquin (2015), Patricia (2015), and Matthew (2016) are conducted during their intensity changes. The impacts of two initialization frameworks on the HWRF analyses and forecasts are compared. It is found that GSI-4DEnVar leads to a reduction in track, minimum sea level pressure (MSLP), and maximum surface wind (MSW) forecast errors in all of the HWRF simulations, compared with the GSI-3DEnVar initialization framework. With assimilating high-resolution observations within the hurricane inner-core region, GSI-4DEnVar can produce the initial hurricane intensity reasonably well without the empirical vortex intensity correction. Further diagnoses with Hurricane Joaquin indicate that GSI-4DEnVar can significantly alleviate the imbalances in the initial conditions and enhance the performance of the data assimilation and subsequent hurricane intensity and precipitation forecasts.


2013 ◽  
Vol 6 (2) ◽  
pp. 3581-3610
Author(s):  
S. Federico

Abstract. This paper presents the current status of development of a three-dimensional variational data assimilation system. The system can be used with different numerical weather prediction models, but it is mainly designed to be coupled with the Regional Atmospheric Modelling System (RAMS). Analyses are given for the following parameters: zonal and meridional wind components, temperature, relative humidity, and geopotential height. Important features of the data assimilation system are the use of incremental formulation of the cost-function, and the use of an analysis space represented by recursive filters and eigenmodes of the vertical background error matrix. This matrix and the length-scale of the recursive filters are estimated by the National Meteorological Center (NMC) method. The data assimilation and forecasting system is applied to the real context of atmospheric profiling data assimilation, and in particular to the short-term wind prediction. The analyses are produced at 20 km horizontal resolution over central Europe and extend over the whole troposphere. Assimilated data are vertical soundings of wind, temperature, and relative humidity from radiosondes, and wind measurements of the European wind profiler network. Results show the validity of the analysis solutions because they are closer to the observations (lower RMSE) compared to the background (higher RMSE), and the differences of the RMSEs are consistent with the data assimilation settings. To quantify the impact of improved initial conditions on the short-term forecast, the analyses are used as initial conditions of a three-hours forecast of the RAMS model. In particular two sets of forecasts are produced: (a) the first uses the ECMWF analysis/forecast cycle as initial and boundary conditions; (b) the second uses the analyses produced by the 3-D-Var scheme as initial conditions, then is driven by the ECMWF forecast. The improvement is quantified by considering the horizontal components of the wind, which are measured at a-synoptic times by the European wind profiler network. The results show that the RMSE is effectively reduced at the short range (1–2 h). The results are in agreement with the set-up of the numerical experiment.


2017 ◽  
Author(s):  
Wei He ◽  
Ivar R. van der Velde ◽  
Arlyn E. Andrews ◽  
Colm Sweeney ◽  
John Miller ◽  
...  

Abstract. We have implemented a regional carbon dioxide data assimilation system based on the CarbonTracker Data Assimilation Shell (CTDAS) and a high-resolution Lagrangian transport model, the Stochastic Time-Inverted Lagrangian Transport model driven by the Weather Forecast and Research meteorological fields (WRF-STILT). With this system, named as CTDAS‑Lagrange, we simultaneously optimize terrestrial biosphere fluxes and four parameters that adjust the lateral boundary conditions (BCs) against CO2 observations from the NOAA ESRL North America tall tower and aircraft Programmable Flask Packages (PFPs) sampling program. Least-squares optimization is performed with a time-stepping ensemble Kalman smoother, over a time window of 10 days and assimilating sequentially a time series of observations. Because the WRF-STILT footprints are pre-computed, it is computationally efficient to run the CTDAS-Lagrange system. To estimate the uncertainties of the optimized fluxes from the system, we performed sensitivity tests with various a priori biosphere fluxes (SiBCASA, SiB3, CT2013B) and BCs (optimized mole fraction fields from CT2013B and CTE2014, and an empirical data set derived from aircraft observations), as well as with a variety of choices on the ways that fluxes are adjusted (additive or multiplicative), covariance length scales, biosphere flux covariances, BC parameter uncertainties, and model-data mismatches. In pseudo-data experiments, we show that in our implementation the additive flux adjustment method is more flexible in optimizing NEE than the multiplicative flux adjustment method, and that the CTDAS-Lagrange system has the ability to correct for the potential biases in the lateral boundary conditions and to resolve large biases in the prior biosphere fluxes. Using real observations, we have derived a range of estimates for the optimized carbon fluxes from a series of sensitivity tests, which places the North American carbon sink for the year 2010 in a range from −0.92 to −1.26 PgC/yr. This is comparable to the TM5-based estimates of CarbonTracker (version CT2016, −0.91 ± 1.10 PgC/yr) and CarbonTracker Europe (version CTE2016, −0.91 ± 0.31 PgC/yr). We conclude that CTDAS-Lagrange can offer a versatile and computationally attractive alternative to these global systems for regional estimates of carbon fluxes, which can take advantage of high-resolution Lagrangian footprints that are increasingly easy to obtain.


2020 ◽  
Vol 21 (9) ◽  
pp. 2023-2039
Author(s):  
Dikra Khedhaouiria ◽  
Stéphane Bélair ◽  
Vincent Fortin ◽  
Guy Roy ◽  
Franck Lespinas

AbstractConsistent and continuous fields provided by precipitation analyses are valuable for hydrometeorological applications and land data assimilation modeling, among others. Providing uncertainty estimates is a logical step in the analysis development, and a consistent approach to reach this objective is the production of an ensemble analysis. In the present study, a 6-h High-Resolution Ensemble Precipitation Analysis (HREPA) was developed for the domain covering Canada and the northern part of the contiguous United States. The data assimilation system is the same as the Canadian Precipitation Analysis (CaPA) and is based on optimal interpolation (OI). Precipitation from the Canadian national 2.5-km atmospheric prediction system constitutes the background field of the analysis, while at-site records and radar quantitative precipitation estimates (QPE) compose the observation datasets. By using stochastic perturbations, multiple observations and background field random realizations were generated to subsequently feed the data assimilation system and provide 24 HREPA members plus one control run. Based on one summer and one winter experiment, HREPA capabilities in terms of bias and skill were verified against at-site observations for different climatic regions. The results indicated HREPA’s reliability and skill for almost all types of precipitation events in winter, and for precipitation of medium intensity in summer. For both seasons, HREPA displayed resolution and sharpness. The overall good performance of HREPA and the lack of ensemble precipitation analysis (PA) at such spatiotemporal resolution in the literature motivate further investigations on transitional seasons and more advanced perturbation approaches.


2017 ◽  
Vol 32 (3) ◽  
pp. 1185-1208 ◽  
Author(s):  
Phillipa Cookson-Hills ◽  
Daniel J. Kirshbaum ◽  
Madalina Surcel ◽  
Jonathan G. Doyle ◽  
Luc Fillion ◽  
...  

Abstract Environment and Climate Change Canada (ECCC) has recently developed an experimental high-resolution EnKF (HREnKF) regional ensemble prediction system, which it tested over the Pacific Northwest of North America for the first half of February 2011. The HREnKF has 2.5-km horizontal grid spacing and assimilates surface and upper-air observations every hour. To determine the benefits of the HREnKF over less expensive alternatives, its 24-h quantitative precipitation forecasts are compared with those from a lower-resolution (15 km) regional ensemble Kalman filter (REnKF) system and to ensembles directly downscaled from the REnKF using the same grid as the HREnKF but with no additional data assimilation (DS). The forecasts are verified against rain gauge observations and gridded precipitation analyses, the latter of which are characterized by uncertainties of comparable magnitude to the model forecast errors. Nonetheless, both deterministic and probabilistic verification indicates robust improvements in forecast skill owing to the finer grids of the HREnKF and DS. The HREnKF exhibits a further improvement in performance over the DS in the first few forecast hours, suggesting a modest positive impact of data assimilation. However, this improvement is not statistically significant and may be attributable to other factors.


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.


Sign in / Sign up

Export Citation Format

Share Document