ALBERT

Description: A fully coupled atmosphere–ocean–ice model has been used to produce global weather forecasts at Environment and Climate Change Canada (ECCC) since November 2017. Currently, the system relies on four uncoupled data assimilation (DA) components for initializing the fully coupled global atmosphere–ocean–ice forecast model: atmosphere, ocean, sea ice and sea surface temperature (SST). The goal of the present study is to implement a weakly coupled data assimilation (WCDA) between the atmosphere and ocean components and evaluate its performance against uncoupled DA. The WCDA system uses coupled atmosphere–ocean–ice short-term forecasts as background states for the atmospheric and the ocean DA components that independently compute atmospheric and ocean analyses. This system leads to better agreement between the coupled atmosphere–ocean analyses and the coupled atmosphere–ocean–ice forecasts than between the uncoupled analyses and the coupled forecasts. The use of WCDA improves the atmospheric forecast score near the surface, but a slight increase in the atmospheric temperature bias is observed. A small positive impact from using the short-term SST forecast on the satellite radiance observation-minus-forecast statistics is noted. Ocean temperature and salinity forecasts are also improved near the surface. The next steps toward stronger DA coupling are highlighted.

Description: This paper presents a reanalysis of the atmospheric chemical composition from the upper troposphere to the lower mesosphere from August 2004 to December 2017. This reanalysis is produced by the Belgian Assimilation System for Chemical ObsErvations (BASCOE) constrained by the chemical observations from the Microwave Limb Sounder (MLS) on board the Aura satellite. BASCOE is based on the ensemble Kalman filter (EnKF) method and includes a chemical transport model driven by the winds and temperature from the ERA-Interim meteorological reanalysis. The model resolution is 3.75∘ in longitude, 2.5∘ in latitude and 37 vertical levels from the surface to 0.1 hPa with 25 levels above 100 hPa. The outputs are provided every 6 h. This reanalysis is called BRAM2 for BASCOE Reanalysis of Aura MLS, version 2. Vertical profiles of eight species from MLS version 4 are assimilated and are evaluated in this paper: ozone (O3), water vapour (H2O), nitrous oxide (N2O), nitric acid (HNO3), hydrogen chloride (HCl), chlorine oxide (ClO), methyl chloride (CH3Cl) and carbon monoxide (CO). They are evaluated using independent observations from the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS), the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) and N2O observations from a different MLS radiometer than the one used to deliver the standard product and ozonesondes. The evaluation is carried out in four regions of interest where only selected species are evaluated. These regions are (1) the lower-stratospheric polar vortex where O3, H2O, N2O, HNO3, HCl and ClO are evaluated; (2) the upper-stratospheric–lower-mesospheric polar vortex where H2O, N2O, HNO3 and CO are evaluated; (3) the upper troposphere–lower stratosphere (UTLS) where O3, H2O, CO and CH3Cl are evaluated; and (4) the middle stratosphere where O3, H2O, N2O, HNO3, HCl, ClO and CH3Cl are evaluated. In general BRAM2 reproduces MLS observations within their uncertainties and agrees well with independent observations, with several limitations discussed in this paper (see the summary in Sect. 5.5). In particular, ozone is not assimilated at altitudes above (i.e. pressures lower than) 4 hPa due to a model bias that cannot be corrected by the assimilation. MLS ozone profiles display unphysical oscillations in the tropical UTLS, which are corrected by the assimilation, allowing a good agreement with ozonesondes. Moreover, in the upper troposphere, comparison of BRAM2 with MLS and independent observations suggests a positive bias in MLS O3 and a negative bias in MLS H2O. The reanalysis also reveals a drift in MLS N2O against independent observations, which highlights the potential use of BRAM2 to estimate biases between instruments. BRAM2 is publicly available and will be extended to assimilate MLS observations after 2017.

Description: 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) data assimilation, using a comprehensive stratospheric chemistry transport model (CTM). This 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 the full chemistry model with EnKF are taken into account. One issue is a large number of error variance parameters that need to be optimized. We estimate an observation error variance parameter as a function of pressure level for each observed species using the Desroziers method. For comparison purposes, we apply the same estimate procedure in the 4D-Var data assimilation, where both scale factors of the background and observation error covariance matrices are estimated using the Desroziers method. However, in EnKF the background error covariance is modelled using the full chemistry model and a model error term which is tuned using an adjustable parameter. We found that it is adequate to have the same value of this parameter based on the chemical tracer formulation that is applied for all observed 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 noise in the cross-covariance between species that occurs when species are weakly chemically related 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 8-month long assimilation of limb sounding observations from EOS Aura MLS. This 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 the presence of model error or observation biases. If the erroneous chemical modelling is associated with moderately fast chemical processes, but whose lifetimes are longer than the model time step, then EnKF performs better, while 4D-Var develops spurious increments in the chemically related species. If, however, the observation biases are significant, then 4D-Var is more robust and is able to reject erroneous observations while EnKF does not.

Description: Ensemble Kalman filter methods are typically used in combination with one of two localization techniques. One technique is covariance localization, or direct forecast error localization, in which the ensemble-derived forecast error covariance matrix is Schur multiplied with a chosen correlation matrix. The second way of localization is by domain decomposition. Here, the assimilation is split into local domains in which the assimilation update is performed independently. Domain localization is frequently used in combination with filter algorithms that use the analysis error covariance matrix for the calculation of the gain like the ensemble transform Kalman filter (ETKF) and the singular evolutive interpolated Kalman filter (SEIK). However, since the local assimilations are performed independently, smoothness of the analysis fields across the subdomain boundaries becomes an issue of concern. To address the problem of smoothness, an algorithm is introduced that uses domain localization in combination with a Schur product localization of the forecast error covariance matrix for each local subdomain. On a simple example, using the Lorenz-40 system, it is demonstrated that this modification can produce results comparable to those obtained with direct forecast error localization. In addition, these results are compared to the method that uses domain localization in combination with weighting of observations. In the simple example, the method using weighting of observations is less accurate than the new method, particularly if the observation errors are small. Domain localization with weighting of observations is further examined in the case of assimilation of satellite data into the global finite-element ocean circulation model (FEOM) using the local SEIK filter. In this example, the use of observational weighting improves the accuracy of the analysis. In addition, depending on the correlation function used for weighting, the spectral properties of the solution can be improved.

Description: This paper presents a reanalysis of the atmospheric chemical composition from the upper troposphere to the lower mesosphere from August 2004 to December 2017. This reanalysis is produced by the Belgian Assimilation System for Chemical ObsErvations (BASCOE) constrained by the chemical observations from the Microwave Limb Sounder (MLS) onboard the Aura satellite. BASCOE is based on the Ensemble Kalman Filter (EnKF) method and includes a chemical transport model driven by the winds and temperature from the ERA-Interim meteorological reanalysis. The model resolution is 3.75° in longitude, 2.5° in latitude and 37 vertical levels from the surface to 0.1 hPa with 25 levels above 100 hPa. The outputs are provided every 6 hours. This reanalysis is called BRAM2 for BASCOE Reanalysis of Aura MLS, version 2. Vertical profiles of eight species from MLS version 4 are assimilated and are evaluated in this paper: ozone (O3), water vapour (H2O), nitrous oxide (N2O), nitric acid (HNO3), hydrogen chloride (HCl), chlorine oxide (ClO), methyl chloride (CH3Cl) and carbon monoxide (CO). They are evaluated using independent observations from the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACEFTS), the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES), N2O observations from another MLS radiometer than the one used to deliver the standard product and ozonesondes. The evaluation is done in four regions of interest where only selected species are evaluated. These regions are (1) the lower stratospheric polar vortex where O3, H2O, N2O, HNO3, HCl and ClO are evaluated, (2) the upper stratospheric lower mesospheric polar vortex where H2O, N2O, HNO3 and CO are evaluated, (3) the tropical tropopause layer (TTL) where O3, H2O, CO and CH3Cl are evaluated and (4) the middle stratosphere where O3, H2O, N2O, HNO3, HCl, ClO and CH3Cl are evaluated. In general BRAM2 reproduces MLS observations within their uncertainties and agrees well with independent observations, with several limitations discussed in this paper (see the summary in Sect. 5.5). In particular, ozone is not assimilated at altitudes above (i.e. pressures lower than) 4 hPa due to a model bias that cannot be corrected by the assimilation. MLS ozone profiles display unphysical oscillations in the TTL which are corrected by the assimilation, allowing a good agreement with ozonesondes. Moreover, in the upper troposphere, comparison of BRAM2 with MLS and independent observations suggests a positive bias in MLS O3 and a negative bias in MLS H2O. The reanalysis also reveals a drift in MLS N2O against independent observations which highlights the potential use of BRAM2 to estimate biases between instruments. BRAM2 is publicly available and will be extended to assimilate MLS observations post 2017.

Description: This paper discusses assimilation experiments of methane (CH4) and nitrous oxide (N2O) profiles retrieved from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Here we focus on data versions 6 and 7 provided by the ESA processor. These data sets have been assimilated by the Belgian Assimilation System for Chemical ObsErvations (BASCOE). The CH4 and N2O retrieved profiles can oscillate, especially in the tropical lower stratosphere. Using the averaging kernels of the observations and a background error covariance matrix, which has previously been calibrated, allows the system to partly remedy this issue and provide assimilated fields that are more regular vertically. In general, there is a good agreement between the BASCOE analyses and independent observations from ACE–FTS (CH4 and N2O) and MLS (N2O), demonstrating the general good quality of CH4 and N2O retrievals provided by MIPAS ESA. Nevertheless, this study also identifies two issues in these data sets. First, time series of the observations show unexpected discontinuities due to an abrupt change in the gain of MIPAS band B, generally occurring after the instrument decontamination. Since the calibration is performed weekly, the abrupt change in the gain affects the measurements until the subsequent calibration is performed. Second, the correlations between BASCOE analyses and independent observations are poor in the lower stratosphere, especially in the tropics, probably due to the presence of outliers in the assimilated data. In this region, we recommend using MIPAS CH4 and N2O retrievals with caution.

Description: A fully coupled atmosphere-ocean-ice model has been used to produce global weather forecasts at Environment and Climate Change Canada (ECCC) since November 2017. Currently, the system relies on four uncoupled data assimilation (DA) components for initializing the fully coupled global atmosphere-ocean-ice forecast model: atmosphere, ocean, sea ice and sea surface temperature (SST). The goal of the present study is to implement a weakly coupled data assimilation (WCDA) between the atmosphere and the ocean components and evaluate its performance against uncoupled DA. The WCDA system uses coupled atmosphere-ocean-ice short-term forecasts as background states for the atmospheric and the ocean DA components that independently compute atmospheric and ocean analyses. This system leads to better agreement between the coupled atmosphere-ocean analyses and the coupled atmosphere-ocean-ice forecasts than between the uncoupled analyses and the coupled forecasts. The use of WCDA improves the atmospheric forecast score near the surface, but a slight increase in the atmospheric temperature bias is observed. A small positive impact from using the short-term SST forecast on the satellite radiance Observation-minus-Forecast statistics is noted. Ocean temperature and salinity forecasts are also improved near the surface. Next steps toward stronger DA coupling are highlighted.

Description: This paper discusses assimilation experiments of methane (CH4) and nitrous oxide (N2O) profiles observed by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS). Here we focus on data versions 6 and 7 retrieved by the ESA processor. These datasets have been assimilated by the Belgian Assimilation System for Chemical ObsErvations (BASCOE). The CH4 and N2O profiles can be noisy especially in the tropical lower stratosphere. Using the averaging kernels of the observations and a background error covariance matrix – the B matrix, which has been previously calibrated, allows the system to partly remedy this issue and provide assimilated fields that are more regular vertically. In general, there is a good agreement between the BASCOE analyses and independent observations demonstrating the general good quality of these two retrievals provided by MIPAS ESA. Nevertheless, this study also identifies two issues in these datasets. First, time-series of the observations show unexpected discontinuities, due to the calibration method used for the level-1 data. Second, the correlations between BASCOE analyses and independent observations are poor in the lower stratosphere, especially in the tropics, probably due to the presence of outliers in the assimilated data. In this region, we recommend using MIPAS CH4 and N2O observations with caution.