scholarly journals Examination of Analysis and Forecast Errors of High-Resolution Assimilation, Bias Removal, and Digital Filter Initialization with an Ensemble Kalman Filter

2012 ◽  
Vol 140 (12) ◽  
pp. 3992-4004 ◽  
Author(s):  
Brian C. Ancell
Keyword(s):  
Kalman Filter ◽  
High Resolution ◽  
Digital Filter ◽  
Ensemble Kalman Filter ◽  
Weather Prediction ◽  
Surface Wind ◽  
Forecast Errors ◽  
Temperature Bias ◽  
High Resolution Analysis ◽  
Bias Removal

Abstract Mesoscale atmospheric data assimilation is becoming an integral part of numerical weather prediction. Modern computational resources now allow assimilation and subsequent forecasting experiments ranging from resolutions of tens of kilometers over regional domains to smaller grids that employ storm-scale assimilation. To assess the value of the high-resolution capabilities involved with assimilation and forecasting at different scales, analyses and forecasts must be carefully evaluated to understand 1) whether analysis benefits gained at finer scales persist into the forecast relative to downscaled runs begun from lower-resolution analyses, 2) how the positive analysis effects of bias removal evolve into the forecast, and 3) how digital filter initialization affects analyses and forecasts. This study applies a 36- and 4-km ensemble Kalman filter over 112 assimilation cycles to address these important issues, which could all be relevant to a variety of short-term, high-resolution, real-time forecasting applications. It is found that with regard to surface wind and temperature, analysis improvements gained at higher resolution persist throughout the 12-h forecast window relative to downscaled, high-resolution forecasts begun from analyses on the coarser grid. Aloft, however, no forecast improvements were found with the high-resolution analysis/forecast runs. Surface wind and temperature bias removal, while clearly improving surface analyses, degraded surface forecasts and showed little forecast influence aloft. Digital filter initialization degraded temperature analyses with or without bias removal, degraded wind analyses when bias removal was used, but improved wind analyses when bias removal was absent. No forecast improvements were found with digital filter initialization. The consequences of these results with regard to operational assimilation/forecasting systems on nested grids are discussed.

scholarly journals Verification of 24-h Quantitative Precipitation Forecasts over the Pacific Northwest from a High-Resolution Ensemble Kalman Filter System

2017 ◽  
Vol 32 (3) ◽  
pp. 1185-1208 ◽  
Author(s):  
Phillipa Cookson-Hills ◽  
Daniel J. Kirshbaum ◽  
Madalina Surcel ◽  
Jonathan G. Doyle ◽  
Luc Fillion ◽  
...  
Keyword(s):  
Kalman Filter ◽  
High Resolution ◽  
Data Assimilation ◽  
Pacific Northwest ◽  
Ensemble Kalman Filter ◽  
Ensemble Prediction ◽  
Forecast Errors ◽  
Ensemble Prediction System ◽  
The Pacific ◽  
Quantitative Precipitation Forecasts

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.

scholarly journals Analysis of a Tornadic Mesoscale Convective Vortex Based on Ensemble Kalman Filter Assimilation of CASA X-Band and WSR-88D Radar Data

2011 ◽  
Vol 139 (11) ◽  
pp. 3446-3468 ◽  
Author(s):  
Nathan Snook ◽  
Ming Xue ◽  
Youngsun Jung
Keyword(s):  
Kalman Filter ◽  
High Resolution ◽  
Ensemble Kalman Filter ◽  
Weather Prediction ◽  
Radar Data ◽  
Convective System ◽  
Foundation Engineering ◽  
X Band ◽  
Regional Prediction ◽  
Mesoscale Convective

Abstract One of the goals of the National Science Foundation Engineering Research Center (ERC) for Collaborative Adaptive Sensing of the Atmosphere (CASA) is to improve storm-scale numerical weather prediction (NWP) by collecting data with a dense X-band radar network that provides high-resolution low-level coverage, and by assimilating such data into NWP models. During the first spring storm season after the deployment of four radars in the CASA Integrated Project-1 (IP-1) network in southwest Oklahoma, a tornadic mesoscale convective system (MCS) was captured by CASA and surrounding Weather Surveillance Radars-1988 Doppler (WSR-88Ds) on 8–9 May 2007. The MCS moved across northwest Texas and western and central Oklahoma; two tornadoes rated as category 1 on the enhanced Fujita scale (EF-1) and one tornado of EF-0 intensity were reported during the event, just to the north of the IP-1 network. This was the first tornadic convective system observed by CASA. To quantify the impacts of CASA radar data in storm-scale NWP, a set of data assimilation experiments were performed using the Advanced Regional Prediction System (ARPS) ensemble Kalman filter (EnKF) system configured with full model physics and high-resolution terrain. Data from four CASA IP-1 radars and five WSR-88Ds were assimilated in some of the experiments. The ensemble contained 40 members, and radar data were assimilated every 5 min for 1 h. While the assimilation of WSR-88D data alone was able to produce a reasonably accurate analysis of the convective system, assimilating CASA data in addition to WSR-88D data is found to improve the representation of storm-scale circulations, particularly in the lowest few kilometers of the atmosphere, as evidenced by analyses of gust front position and comparison of simulated Vr with observations. Assimilating CASA data decreased RMS innovation of the resulting ensemble mean analyses of Z, particularly in early assimilation cycles, suggesting that the addition of CASA data allowed the EnKF system to more quickly achieve a good result. Use of multiple microphysics schemes in the forecast ensemble was found to alleviate underdispersion by increasing the ensemble spread. This work is the first assimilating real CASA data into an NWP model using EnKF.

Assimilation of Simulated Doppler Radar Observations with an Ensemble Kalman Filter*

2003 ◽  
Vol 131 (8) ◽  
pp. 1663-1677 ◽  
Author(s):  
Chris Snyder ◽  
Fuqing Zhang
Keyword(s):  
Kalman Filter ◽  
Radial Velocity ◽  
Ensemble Kalman Filter ◽  
Doppler Radar ◽  
Weather Prediction ◽  
Forecast Model ◽  
Forecast Errors ◽  
Radar Observations ◽  
Cloud Models ◽  
Reference Simulation

Abstract Assimilation of Doppler radar data into cloud models is an important obstacle to routine numerical weather prediction for convective-scale motions; the difficulty lies in initializing fields of wind, temperature, moisture, and condensate given only observations of radial velocity and reflectivity from the radar. This paper investigates the potential of the ensemble Kalman filter (EnKF), which estimates the covariances between observed variables and the state through an ensemble of forecasts, to assimilate radar observations at convective scales. In the basic experiment, simulated observations are extracted from a reference simulation of a splitting supercell and assimilated using the EnKF and the same numerical model that produced the reference simulation. The EnKF produces accurate analyses, including the unobserved variables, after roughly 30 min (or six scans) of radial velocity observations. Additional experiments, in which forecasts are made from the ensemble-mean analysis, reveal that forecast errors grow significantly in this simple system, so that the ability of the EnKF to track the reference solution is not simply because of stable system dynamics. It is also found that the covariances between radial velocity and temperature, moisture, and condensate are important to the quality of the analyses, as is the initialization chosen for the ensemble members prior to assimilating the first observations. These results are promising, especially given the ease of implementing the EnKF. A number of important issues remain, however, including the initialization of the ensemble prior to the first observation, the treatment of uncertainty in the environmental sounding, the role of error in the forecast model (particularly the microphysical parameterizations), and the treatment of lateral boundary conditions.

scholarly journals Assimilation of Tropical Cyclone Track and Structure Based on the Ensemble Kalman Filter (EnKF)

2010 ◽  
Vol 67 (12) ◽  
pp. 3806-3822 ◽  
Author(s):  
Chun-Chieh Wu ◽  
Guo-Yuan Lien ◽  
Jan-Huey Chen ◽  
Fuqing Zhang
Keyword(s):  
Kalman Filter ◽  
High Resolution ◽  
Tropical Cyclone ◽  
Ensemble Kalman Filter ◽  
Wind Profile ◽  
Surface Wind ◽  
Wrf Model ◽  
Wind Structure ◽  
Two Parameters ◽  
Satellite Analysis

Abstract A new tropical cyclone vortex initialization method based on the ensemble Kalman filter (EnKF) is proposed in this study. Three observed parameters that are related to the tropical cyclone (TC) track and structure—center position, velocity of storm motion, and surface axisymmetric wind structure—are assimilated into the high-resolution Weather Research and Forecasting (WRF) model during a 24-h initialization period to develop a dynamically balanced TC vortex without employing any extra bogus schemes. The first two parameters are available from the TC track data of operational centers, which are mainly based on satellite analysis. The radial wind profile is constructed by fitting the combined information from both the best-track and the dropwindsonde data available from aircraft surveillance observations, such as the Dropwindsonde Observations for Typhoon Surveillance near the Taiwan Region (DOTSTAR). The initialized vortex structure is consistent with the observations of a typical vertical TC structure, even though only the surface wind profile is assimilated. In addition, the subsequent numerical integration shows minor adjustments during early periods, indicating that the analysis fields obtained from this method are dynamically balanced. Such a feature is important for TC numerical integrations. The results here suggest that this new method promises an improved TC initialization and could possibly contribute to some high-resolution numerical experiments to better understand the dynamics of TC structure and to improve operational TC model forecasts. Further applications of this method with sophisticated data from The Observing System Research and Predictability Experiment (THORPEX) Pacific Asian Regional Campaign (T-PARC) will be shown in a follow-up paper.

Assessing Predictability with a Local Ensemble Kalman Filter

2007 ◽  
Vol 64 (4) ◽  
pp. 1116-1140 ◽  
Author(s):  
David Kuhl ◽  
Istvan Szunyogh ◽  
Eric J. Kostelich ◽  
Gyorgyi Gyarmati ◽  
D. J. Patil ◽  
...  
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Ensemble Kalman Filter ◽  
Initial Conditions ◽  
Forecast Error ◽  
Weather Prediction ◽  
Forecast Errors ◽  
Data Assimilation Scheme ◽  
Ensemble Mean ◽  
Assimilation Scheme

Abstract In this paper, the spatiotemporally changing nature of predictability is studied in a reduced-resolution version of the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS), a state-of-the-art numerical weather prediction model. Atmospheric predictability is assessed in the perfect model scenario for which forecast uncertainties are entirely due to uncertainties in the estimates of the initial states. Uncertain initial conditions (analyses) are obtained by assimilating simulated noisy vertical soundings of the “true” atmospheric states with the local ensemble Kalman filter (LEKF) data assimilation scheme. This data assimilation scheme provides an ensemble of initial conditions. The ensemble mean defines the initial condition of 5-day deterministic model forecasts, while the time-evolved members of the ensemble provide an estimate of the evolving forecast uncertainties. The observations are randomly distributed in space to ensure that the geographical distribution of the analysis and forecast errors reflect predictability limits due to the model dynamics and are not affected by inhomogeneities of the observational coverage. Analysis and forecast error statistics are calculated for the deterministic forecasts. It is found that short-term forecast errors tend to grow exponentially in the extratropics and linearly in the Tropics. The behavior of the ensemble is explained by using the ensemble dimension (E dimension), a spatiotemporally evolving measure of the evenness of the distribution of the variance between the principal components of the ensemble-based forecast error covariance matrix. It is shown that in the extratropics the largest forecast errors occur for the smallest E dimensions. Since a low value of the E dimension guarantees that the ensemble can capture a large portion of the forecast error, the larger the forecast error the more certain that the ensemble can fully capture the forecast error. In particular, in regions of low E dimension, ensemble averaging is an efficient error filter and the ensemble spread provides an accurate prediction of the upper bound of the error in the ensemble-mean forecast.

scholarly journals Data Assimilation of the High-Resolution Sea Surface Temperature Obtained from the Aqua-Terra Satellites (MODIS-SST) Using an Ensemble Kalman Filter

2013 ◽  
Vol 5 (6) ◽  
pp. 3123-3139 ◽  
Author(s):  
Yasumasa Miyazawa ◽  
Hiroshi Murakami ◽  
Toru Miyama ◽  
Sergey Varlamov ◽  
Xinyu Guo ◽  
...  
Keyword(s):  
Kalman Filter ◽  
High Resolution ◽  
Data Assimilation ◽  
Surface Temperature ◽  
Sea Surface Temperature ◽  
Ensemble Kalman Filter ◽  
Sea Surface

Evaluation of Surface Analyses and Forecasts with a Multiscale Ensemble Kalman Filter in Regions of Complex Terrain

2011 ◽  
Vol 139 (6) ◽  
pp. 2008-2024 ◽  
Author(s):  
Brian C. Ancell ◽  
Clifford F. Mass ◽  
Gregory J. Hakim
Keyword(s):  
Kalman Filter ◽  
High Resolution ◽  
Data Assimilation ◽  
Ensemble Kalman Filter ◽  
Complex Terrain ◽  
Error Variance ◽  
Small Scale ◽  
Observation Error ◽  
Grid Spacing ◽  
Background Error

Abstract Previous research suggests that an ensemble Kalman filter (EnKF) data assimilation and modeling system can produce accurate atmospheric analyses and forecasts at 30–50-km grid spacing. This study examines the ability of a mesoscale EnKF system using multiscale (36/12 km) Weather Research and Forecasting (WRF) model simulations to produce high-resolution, accurate, regional surface analyses, and 6-h forecasts. This study takes place over the complex terrain of the Pacific Northwest, where the small-scale features of the near-surface flow field make the region particularly attractive for testing an EnKF and its flow-dependent background error covariances. A variety of EnKF experiments are performed over a 5-week period to test the impact of decreasing the grid spacing from 36 to 12 km and to evaluate new approaches for dealing with representativeness error, lack of surface background variance, and low-level bias. All verification in this study is performed with independent, unassimilated observations. Significant surface analysis and 6-h forecast improvements are found when EnKF grid spacing is reduced from 36 to 12 km. Forecast improvements appear to be a consequence of increased resolution during model integration, whereas analysis improvements also benefit from high-resolution ensemble covariances during data assimilation. On the 12-km domain, additional analysis improvements are found by reducing observation error variance in order to address representativeness error. Removing model surface biases prior to assimilation significantly enhances the analysis. Inflating surface wind and temperature background error variance has large impacts on analyses, but only produces small improvements in analysis RMS errors. Both surface and upper-air 6-h forecasts are nearly unchanged in the 12-km experiments. Last, 12-km WRF EnKF surface analyses and 6-h forecasts are shown to generally outperform those of the Global Forecast System (GFS), North American Model (NAM), and the Rapid Update Cycle (RUC) by about 10%–30%, although these improvements do not extend above the surface. Based on these results, future improvements in multiscale EnKF are suggested.

scholarly journals A Global Carbon Assimilation System using a modified EnKF assimilation method

2014 ◽  
Vol 7 (5) ◽  
pp. 6519-6547
Author(s):  
S. Zhang ◽  
X. Zheng ◽  
Z. Chen ◽  
B. Dan ◽  
J. M. Chen ◽  
...  
Keyword(s):  
Kalman Filter ◽  
Ensemble Kalman Filter ◽  
Carbon Assimilation ◽  
Terrestrial Ecosystem ◽  
Carbon Fluxes ◽  
Ecosystem Model ◽  
Atmospheric Co2 ◽  
Forecast Errors ◽  
Global Carbon ◽  
Assimilation System

Abstract. A Global Carbon Assimilation System based on Ensemble Kalman filter (GCAS-EK) is developed for assimilating atmospheric CO2 abundance data into an ecosystem model to simultaneously estimate the surface carbon fluxes and atmospheric CO2 distribution. This assimilation approach is based on the ensemble Kalman filter (EnKF), but with several new developments, including using analysis states to iteratively estimate ensemble forecast errors, and a maximum likelihood estimation of the inflation factors of the forecast and observation errors. The proposed assimilation approach is tested in observing system simulation experiments and then used to estimate the terrestrial ecosystem carbon fluxes and atmospheric CO2 distributions from 2002 to 2008. The results showed that this assimilation approach can effectively reduce the biases and uncertainties of the carbon fluxes simulated by the ecosystem model.

scholarly journals Time Series High-Resolution Land Surface Albedo Estimation Based on the Ensemble Kalman Filter Algorithm

2019 ◽  
Vol 11 (7) ◽  
pp. 753 ◽  
Author(s):  
Guodong Zhang ◽  
Hongmin Zhou ◽  
Changjing Wang ◽  
Huazhu Xue ◽  
Jindi Wang ◽  
...  
Keyword(s):  
Time Series ◽  
Kalman Filter ◽  
High Resolution ◽  
Spatial Resolution ◽  
Ensemble Kalman Filter ◽  
Land Surface ◽  
Surface Albedo ◽  
Data Sets ◽  
Broadleaf Forest ◽  
Kalman Filter Algorithm

Continuous, long-term sequence, land surface albedo data have crucial significance for climate simulations and land surface process research. Sensors such as the Moderate-Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer (VIIRS) provide global albedo product data sets with a spatial resolution of 500 m over long time periods. There is demand for new high-resolution albedo data for regional applications. High-resolution observations are often unavailable due to cloud contamination, which makes it difficult to obtain time series albedo estimations. This paper proposes an “amalgamation albedo“ approach to generate daily land surface shortwave albedo with 30 m spatial resolution using Landsat data and the MODIS Bidirectional Reflectance Distribution Functions (BRDF)/Albedo product MCD43A3 (V006). Historical MODIS land surface albedo products were averaged to obtain an albedo estimation background, which was used to construct the albedo dynamic model . The Thematic Mapper (TM) albedo derived via direct estimation approach was then introduced to generate high spatial-temporal resolution albedo data based on the Ensemble Kalman Filter algorithm (EnKF). Estimation results were compared to field observations for cropland, deciduous broadleaf forest, evergreen needleleaf forest, grassland, and evergreen broadleaf forest domains. The results indicated that for all land cover types, the estimated albedos coincided with ground measurements at a root mean squared error (RMSE) of 0.0085–0.0152. The proposed algorithm was then applied to regional time series albedo estimation; the results indicated that it captured spatial and temporal variation patterns for each site. Taken together, our results suggest that the amalgamation albedo approach is a feasible solution to generate albedo data sets with high spatio-temporal resolution.

Sign in / Sign up

Export Citation Format

Share Document