scholarly journals Comparison between Local Ensemble Transform Kalman Filter and PSAS in the NASA finite volume GCM – perfect model experiments

2008 ◽  
Vol 15 (4) ◽  
pp. 645-659 ◽  
Author(s):  
J. Liu ◽  
E. J. Fertig ◽  
H. Li ◽  
E. Kalnay ◽  
B. R. Hunt ◽  
...  

Abstract. This paper compares the performance of the Local Ensemble Transform Kalman Filter (LETKF) with the Physical-Space Statistical Analysis System (PSAS) under a perfect model scenario. PSAS is a 3D-Var assimilation system used operationally in the Goddard Earth Observing System Data Assimilation System (GEOS-4 DAS). The comparison is carried out using simulated winds and geopotential height observations and the finite volume Global Circulation Model with 72 grid points zonally, 46 grid points meridionally and 55 vertical levels. With forty ensemble members, the LETKF obtains analyses and forecasts with significantly lower RMS errors than those from PSAS, especially over the Southern Hemisphere and oceans. This observed advantage of the LETKF over PSAS is due to the ability of the 40-member ensemble LETKF to capture flow-dependent errors and thus create a good estimate of the evolving background uncertainty. An initial decrease of the forecast errors in the Northern Hemisphere observed in the PSAS but not in the LETKF suggests that the LETKF analysis is more balanced.

2019 ◽  
Vol 12 (7) ◽  
pp. 2899-2914
Author(s):  
Yun Liu ◽  
Eugenia Kalnay ◽  
Ning Zeng ◽  
Ghassem Asrar ◽  
Zhaohui Chen ◽  
...  

Abstract. We developed a carbon data assimilation system to estimate surface carbon fluxes using the local ensemble transform Kalman filter (LETKF) and atmospheric transport model GEOS-Chem driven by the MERRA-1 reanalysis of the meteorological field based on the Goddard Earth Observing System model, version 5 (GEOS-5). This assimilation system is inspired by the method of Kang et al. (2011, 2012), who estimated the surface carbon fluxes in an observing system simulation experiment (OSSE) as evolving parameters in the assimilation of the atmospheric CO2, using a short assimilation window of 6 h. They included the assimilation of the standard meteorological variables, so that the ensemble provided a measure of the uncertainty in the CO2 transport. After introducing new techniques such as “variable localization”, and increased observation weights near the surface, they obtained accurate surface carbon fluxes at grid-point resolution. We developed a new version of the local ensemble transform Kalman filter related to the “running-in-place” (RIP) method used to accelerate the spin-up of ensemble Kalman filter (EnKF) data assimilation (Kalnay and Yang, 2010; Wang et al., 2013; Yang et al., 2012). Like RIP, the new assimilation system uses the “no cost smoothing” algorithm for the LETKF (Kalnay et al., 2007b), which allows shifting the Kalman filter solution forward or backward within an assimilation window at no cost. In the new scheme a long “observation window” (e.g., 7 d or longer) is used to create a LETKF ensemble at 7 d. Then, the RIP smoother is used to obtain an accurate final analysis at 1 d. This new approach has the advantage of being based on a short assimilation window, which makes it more accurate, and of having been exposed to the future 7 d observations, which improves the analysis and accelerates the spin-up. The assimilation and observation windows are then shifted forward by 1 d, and the process is repeated. This reduces significantly the analysis error, suggesting that the newly developed assimilation method can be used with other Earth system models, especially in order to make greater use of observations in conjunction with models.


2017 ◽  
Author(s):  
Yun Liu ◽  
Eugenia Kalnay ◽  
Ning Zeng ◽  
Ghassem Asrar ◽  
Zhaohui Chen ◽  
...  

Abstract. We developed a Carbon data assimilation system to estimate the surface carbon fluxes using the Local Ensemble Transform Kalman Filter and atmospheric transfer model of GEOS-Chem driven by the MERRA-1 reanalysis of the meteorological fields based on the Goddard Earth Observing System Model, Version 5 (GEOS-5). This assimilation system is inspired by the method of Kang et al. [2011, 2012], who estimated the surface carbon fluxes in an Observing System Simulation Experiment (OSSE) mode, as evolving parameters in the assimilation of the atmospheric CO2, using a short assimilation window of 6 hours. They included the assimilation of the standard meteorological variables, so that the ensemble provided a measure of the uncertainty in the CO2 transport. After introducing new techniques such as variable localization, and increased observation weights near the surface, they obtained accurate carbon fluxes at grid point resolution. We developed a new version of the LETKF related to the Running-in-Place (RIP) method used to accelerate the spin-up of EnKF data assimilation [Kalnay and Yang, 2010; Wang et al., 2013, Yang et al., 2014]. Like RIP, the new assimilation system uses the no-cost smoothing algorithm for the LETKF [Kalnay et al., 2007b], which allows shifting at no cost the Kalman Filter solution forward or backward within an assimilation window. In the new scheme a long observation window (e.g., 7-days or longer) is used to create an LETKF ensemble at 7-days. Then, the RIP smoother is used to obtain an accurate final analysis at 1-day. This analysis has the advantage of being based on a short assimilation window, which makes it more accurate, and of having been exposed to the future 7-days observations, which accelerates the spin up. The assimilation and observation windows are then shifted forward by one day, and the process is repeated. This reduces significantly the analysis error, suggesting that this method could be used in other data assimilation problems.


2019 ◽  
Vol 9 ◽  
pp. A30 ◽  
Author(s):  
Sean Elvidge ◽  
Matthew J. Angling

The Advanced Ensemble electron density (Ne) Assimilation System (AENeAS) is a new data assimilation model of the ionosphere/thermosphere. The background model is provided by the Thermosphere Ionosphere Electrodynamics General Circulation Model (TIE-GCM) and the assimilation uses the local ensemble transform Kalman filter (LETKF). An outline derivation of the LETKF is provided and the equations are presented in a form analogous to the classic Kalman filter. An enhancement to the efficient LETKF implementation to reduce computational cost is also described. In a 3 day test in June 2017, AENeAS exhibits a total electron content (TEC) RMS error of 2.1 TECU compared with 5.5 TECU for NeQuick and 6.8 for TIE-GCM (with an NeQuick topside).


2008 ◽  
Vol 60 (1) ◽  
pp. 113-130 ◽  
Author(s):  
Istvan Szunyogh ◽  
Eric J. Kostelich ◽  
Gyorgyi Gyarmati ◽  
Eugenia Kalnay ◽  
Brian R. Hunt ◽  
...  

2008 ◽  
Vol 25 (9) ◽  
pp. 1638-1656 ◽  
Author(s):  
Ross N. Hoffman ◽  
Rui M. Ponte ◽  
Eric J. Kostelich ◽  
Alan Blumberg ◽  
Istvan Szunyogh ◽  
...  

Abstract Data assimilation approaches that use ensembles to approximate a Kalman filter have many potential advantages for oceanographic applications. To explore the extent to which this holds, the Estuarine and Coastal Ocean Model (ECOM) is coupled with a modern data assimilation method based on the local ensemble transform Kalman filter (LETKF), and a series of simulation experiments is conducted. In these experiments, a long ECOM “nature” run is taken to be the “truth.” Observations are generated at analysis times by perturbing the nature run at randomly chosen model grid points with errors of known statistics. A diverse collection of model states is used for the initial ensemble. All experiments use the same lateral boundary conditions and external forcing fields as in the nature run. In the data assimilation, the analysis step combines the observations and the ECOM forecasts using the Kalman filter equations. As a control, a free-running forecast (FRF) is made from the initial ensemble mean to check the relative importance of external forcing versus data assimilation on the analysis skill. Results of the assimilation cycle and the FRF are compared to truth to quantify the skill of each. The LETKF performs well for the cases studied here. After just a few assimilation cycles, the analysis errors are smaller than the observation errors and are much smaller than the errors in the FRF. The assimilation quickly eliminates the domain-averaged bias of the initial ensemble. The filter accurately tracks the truth at all data densities examined, from observations at 50% of the model grid points down to 2% of the model grid points. As the data density increases, the ensemble spread, bias, and error standard deviation decrease. As the ensemble size increases, the ensemble spread increases and the error standard deviation decreases. Increases in the size of the observation error lead to a larger ensemble spread but have a small impact on the analysis accuracy.


Sign in / Sign up

Export Citation Format

Share Document