EnKF and 4D-Var Data Assimilation with a Chemistry Transport Model
Bibliographic record
Abstract
Abstract. We compare two optimized chemical data assimilation systems, one based on the ensemble Kalman filter (EnKF) and the other based on four-dimensional variational (4D-Var), using a comprehensive stratospheric chemistry transport model (CTM). The work is an extension of the Belgian Assimilation System for Chemical ObsErvations (BASCOE), initially designed to work with a 4D-Var data assimilation. A strict comparison of both methods in the case of chemical tracer transport was done in a previous study and indicated that both methods provide essentially similar results. In the present work, we assimilate observations of ozone, HCl, HNO3, H2O and N2O from EOS Aura-MLS data into the BASCOE CTM with a full description of stratospheric chemistry. Two new issues related to the use of full chemistry model with EnKF are taken into account. One issue concerns to a large number of error variance parameters that need to be optimized. We estimate an observation error parameter as function of pressure level for each observed species using the Desroziers' method. For comparison reasons, we apply the same estimate procedure in the 4D-Var data assimilation, where we keep both estimates: the background and observation error variances. However in EnKF, the background error covariance is modelled using the full chemistry model and a model error term. We found that it is adequate to have a single model error based on the chemical tracer formulation that is applied for all species. This is an indication that the main source of model error in chemical transport model is due to the transport. The second issue in EnKF with comprehensive atmospheric chemistry models is the sampling errors between species. When species are weakly chemically related, cross-species sampling noise errors occur at the same location. These errors need to be filtered out, in addition to a localization based on distance. The performance of two data assimilation methods was assessed through an eight-month long assimilation of limb sounding observations from EOS Aura-MLS. The paper discusses the differences in results and their relation to stratospheric chemical processes. Generally speaking, EnKF and 4D-Var provide results of comparable quality but differ substantially in presence of model error or observation biases. If the erroneous chemical modelling is associated with not too small chemical life-times, then EnKF performs better, while 4D-Var develops spurious increments in the chemically related species. If, on the other hand, the observation biases are significant, then 4D-Var is more robust and is able to reject erroneous observations, while EnKF does not.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.001 | 0.002 |
| Meta-epidemiology (narrow) | 0.001 | 0.000 |
| Meta-epidemiology (broad) | 0.001 | 0.001 |
| Bibliometrics | 0.000 | 0.001 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.001 | 0.001 |
| Open science | 0.001 | 0.001 |
| Research integrity | 0.001 | 0.001 |
| Insufficient payload (model declined to judge) | 0.001 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".