the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Estimation of CH_{4} emission based on an advanced 4DLETKF assimilation system
Jagat S. H. Bisht
Prabir K. Patra
Masayuki Takigawa
Takashi Sekiya
Yugo Kanaya
Naoko Saitoh
Kazuyuki Miyazaki
Methane (CH_{4}) is the second major greenhouse gas after carbon dioxide (CO_{2}) which has substantially increased during recent decades in the atmosphere, raising serious sustainability and climate change issues. Here, we develop a data assimilation system for in situ and columnaveraged concentrations using a local ensemble transform Kalman filter (LETKF) to estimate surface emissions of CH_{4}. The data assimilation performance is tested and optimized based on idealized settings using observation system simulation experiments (OSSEs), where a known surface emission distribution (the truth) is retrieved from synthetic observations. We tested three covariance inflation methods to avoid covariance underestimation in the emission estimates, namely fixed multiplicative (FM), relaxationtoprior spread (RTPS), and adaptive multiplicative. First, we assimilate the synthetic observations at every grid point at the surface level. In such a case of dense observational data, the normalized root mean square error (RMSE) in the analyses over global land regions is smaller by 10 %–15 % in the case of RTPS covariance inflation method compared to FM. We have shown that integrated estimated flux seasonal cycles over 15 regions using RTPS inflation are in reasonable agreement between true and estimated flux, with 0.04 global normalized annual mean bias. We then assimilated the columnaveraged CH_{4} concentration by sampling the model simulations at Greenhouse Gases Observing Satellite (GOSAT) observation locations and time for another OSSE. Similar to the case of dense observational data, the RTPS covariance inflation method performs better than FM for GOSAT synthetic observation in terms of normalized RMSE (2 %–3 %) and integrated flux estimation comparison with the true flux. The annual mean averaged normalized RMSE (normalized mean bias) in LETKF CH_{4} flux estimation in the case of RTPS and FM covariance inflation is found to be 0.59 (0.18) and 0.61 (0.23), respectively. The χ^{2} test performed for GOSAT synthetic observations assimilation suggests high underestimation of background error covariance in both RTPS and FM covariance inflation methods; however, the underestimation is much higher (>100 % always) for FM compared to RTPS covariance inflation method.
 Article
(5222 KB)  Fulltext XML

Supplement
(2856 KB)  BibTeX
 EndNote
Methane (CH_{4}) is the second major greenhouse gas, after carbon dioxide (CO_{2}), that has anthropogenic sources. According to the contemporary record of the global CH_{4} budget, the total of all CH_{4} sources ranged 538–593 Tg yr^{−1} during 2008–2017 (Saunois et al., 2020). The primary natural sources are from wetlands (∼ 40 %). The main anthropogenic CH_{4} emissions are from microbial emissions associated with ruminants (livestock and waste), rice cultivation, fugitive emissions (oil and gas production and use), and incomplete combustion of biofuels and fossil fuels. The major fraction of atmospheric CH_{4} sinks (range: 474–532 Tg yr^{−1}) occurs in the troposphere by oxidation via reaction with hydroxyl (OH) radicals (Patra, et al., 2011; Saunois et al., 2020); other loss processes include oxidation by soil and reactions with O^{1}D and Cl. The lifetime of CH_{4} in the atmosphere is estimated to be 9.1 ± 0.9 years (Szopa et al., 2021).
Regional CH_{4} emissions can be estimated from CH_{4} concentration fields and chemistry transport models using Bayesian synthesis approaches based on inverse modeling techniques (e.g., Enting, 2002). In such approaches, emissions are optimized on a coarse resolution (e.g., for a limited number of predefined regions) mostly using surfacebased observations. CH_{4} concentrations are provided by the NOAA cooperative air sampling network sites (Lan et al., 2022) and other networks by the World Data Centre for Greenhouse Gases (WDCGG) website, hosted by the Japan Meteorological Agency. In recent years, satellite measurements have been made by the Greenhouse Gases Observing Satellite (GOSAT) or the TROPOspheric Monitoring Instrument (TROPOMI) (Lorente et al., 2021), covering the globe with fine spatiotemporal scales. GOSAT has provided an extensive global observations of column CH_{4} concentrations since 2009 (Yoshida et al., 2013). Some of the inverse modeling studies utilize the satellite observations for CH_{4} flux estimation (Zhang et al., 2021; Maasakkers et al., 2016), but this requires enormous computational resources as a result of dealing with more flux regions and more observations.
Gridbased CH_{4} flux optimization is also performed using adjoint technique (4D Var data assimilation) and an ensemble Kalman filter (EnKF) but was limited to small sets of observations (Houweling et al., 1999; Meirink et al., 2008; Bruhwiler et al., 2014). Bruhwiler et al. (2014) followed the EnKF method of Peters et al. (2005) to estimate the CH_{4} surface fluxes that utilizes an offline atmospheric chemistry tracer model (ACTM) framework. Techniques such as 4D Var and EnKF are important to estimate CH_{4} fluxes since they can assimilate a large number of observations and manage highresolution fluxes. In the EnKF system, a flowdependent forecast error covariance structure is provided by ensemble model forecasts, while it does not need an adjoint model, which makes it a simple but powerful tool for flux estimation. One of the limitations of the EnKF method is the dependence of the resolution of state vector on ensemble size, which can give spurious results if the number of ensemble members is much smaller than the rank of the error covariance matrix (Houtekamer and Zhang, 2016).
A local ensemble transform Kalman filter (LETKF) is a type of squareroot EnKF that performs analysis locally in space without perturbing the observations (Ott et al., 2002, 2004; Hunt et al., 2007). LETKFs are computationally efficient since the observations are assimilated simultaneously and not serially; it is simple to account for observation error correlation. Miyazaki et al. (2011) and Kang et al. (2012) demonstrated the implementation of LETKF data assimilation system by coupling an ACTM for carboncycle research using atmospheric CO_{2} observations. It is also extensively applied for the emission estimation of shortlived species using satellite data (Skachko et al., 2016; Miyazaki et al., 2019; Sekiya et al., 2021). In this work, we will estimate the CH_{4} fluxes using a LETKF data assimilation system. Assimilation windows ranging from 6 h (Kang et al., 2012) to several months (Bruhwiler et al., 2014) have been used, depending on the desired time resolution of the estimated emissions, which is often limited by the observational data density. The time frame over which the system behaves linearly and in what time frame the observations respond to the control variables, such as atmospheric transport, as well as observation abundance, must also be taken into consideration. Within an assimilation window, where and when the fluxes would be constrained by specific observations is to be ascertained by the correlation between ensemble prior fluxes and the ensemble CH_{4} concentration simulation from a forward model (Liu et al., 2016).
The main objective of this work is to develop an advanced 4D data assimilation system based on a LETKF that simultaneously estimates atmospheric distributions and surface fluxes of CH_{4}. Observation system simulation experiments (OSSEs) are conducted to assess the performance of the LETKF since it is important to test the system against the known emissions or the truth. The OSSE LETKF setup of topdown CH_{4} flux estimation using an online ACTM is an essential step before implementation in real in situ and satellite observation.
We briefly describe the LETKF in the application of CH_{4} flux estimation, while detailed derivation of equations and code implementation are given elsewhere (Hunt et al., 2007; Miyazaki et al., 2011; Miyoshi et al., 2010). The notation used here for LETKF formulation is adopted from Kotsuki et al. (2017). In the LETKF, the background ensemble (columns of matrix x^{b}) in a local region evolved from a set of perturbed initial conditions. The background ensemble mean, ${\stackrel{\mathrm{\u203e}}{x}}^{\mathrm{b}}$, and its perturbation, X^{b}, are estimated from the ensemble forecast as
where m indicates the ensemble size. The background error covariance matrix P^{b} in the mdimensional ensemble is defined as
The analysis ensemble mean ${\stackrel{\mathrm{\u203e}}{x}}^{\mathrm{a}}$ is derived using background ensemble mean ${\stackrel{\mathrm{\u203e}}{x}}^{\mathrm{b}}$ and ensemble perturbations X^{b} as
where H, Y, R, and ${\stackrel{\mathrm{\u0303}}{P}}^{\mathrm{a}}$ denote the linear observation operator, ensemble perturbation matrix in the observation space (Y≡Hx), observation error covariance matrix, and analysis error covariance matrix in the ensemble space, respectively. The superscripts “o”, “b”, and “a” denote the observations, background (prior), and analysis (posterior), respectively. w^{a} defines the analysis increment (or analysis weight) in observation space and is derived using the information about observational increment ${y}^{\mathrm{o}}H{\stackrel{\mathrm{\u203e}}{x}}^{\mathrm{b}}$. The analysis error covariance matrix (${\stackrel{\mathrm{\u0303}}{P}}^{\mathrm{a}}$) in the mdimensional ensemble space is spanned by ensemble perturbation (Hunt et al., 2007) and defined as
Finally, the analysis ensemble perturbations X^{a} at the central grid point are derived such as
where $\mathit{\left\{}\right(m\mathrm{1}){\stackrel{\mathrm{\u0303}}{P}}^{\mathrm{a}}{\mathit{\}}}^{\mathrm{1}/\mathrm{2}}$ is a multiple of the symmetric square root of the local analysis error covariance matrix in ensemble space and could be computed by a singular vector decomposition method. The LETKF solves the analysis update (Eqs. 3 and 5) at every model grid point independently by assimilating local observations within the localization cutoff radius.
We have applied a gross error check as a quality control to exclude observations that are far from the first guess; the appropriate degrees of the gross error check are also examined. Figure 1 shows the schematic diagram of our LETKF setup with two ensemble members for three consecutive assimilation cycles with an 8 d assimilation window. The analysis is obtained at the midpoint time of the assimilation window (Fig. 1). The analyzed (updated) surface flux is used for the next data assimilation cycle starting from the midpoint time of the previous data assimilation window. The state vector augmentation approach is used to estimate the atmospheric CH_{4} surface flux (Kang et al., 2012; Miyazaki et al., 2011).
Assimilation window size and ensemble members are chosen based on computational efficiency and estimation accuracy. A larger assimilation window means fluxes are constrained by more observations; however, it requires handling of large matrix optimization which is difficult in cases of dense observation and introduces sampling errors related to transport errors. In this study, a few sensitivity experiments were performed to demonstrate the choice of assimilation window length and ensemble size when GOSAT synthetic observations are assimilated in Sect. 4.2.
2.1 Covariance inflation
The LETKF data assimilation needs variance inflation to mitigate the underdispersed ensemble. We tested three methods: fixed multiplicative (FM), relaxationtoprior spread (RTPS), and adaptive multiplicative covariance inflation.
The fixed multiplicative (FM) inflation method (Anderson and Anderson, 1999) inflates the prior ensemble by inflating the background error covariance matrix P^{b} defined in Eq. (2) as
where ${\mathbf{P}}_{\mathrm{tmp}}^{\mathrm{b}}$ represents the temporary background error covariance matrix, which is inflated by a factor γ.
The other inflation methods used to prevent the reduction of ensemble spread are relaxationtoprior perturbation (RTPP) (Zhang et al., 2004) and relaxationtoprior spread (RTPS) (Whitaker and Hamill, 2012). The RTPP method relaxes the reduction of the ensemble spread after updating the ensemble perturbations, which blends the background and analysis ensemble perturbations as
where α_{RTPP} denotes the relaxation parameter of the RTPP.
The RTPS inflation method relaxes the reduction of the ensemble spread by relaxing the analysis spread to prior spread as
where σ and α_{RTPS} denote the ensemble spread and relaxation parameter of the RTPS, respectively. The range of the α_{RTPS} parameter is bounded by [0, 1]. This study focuses mainly on the FM and RTPS covariance inflation methods.
In addition, Miyoshi (2011) applied adaptive inflation by determining the multiplicative inflation factors at every grid point at every analysis step using the observationspace statistics derived by Daley (1992) and Desroziers et al. (2005).
where the operator “<⋅>” denotes the statistical expectation and $d={y}^{\mathrm{o}}H{\stackrel{\mathrm{\u203e}}{x}}^{\mathrm{b}}$ (observation minus first guess), and R is the error observation covariance matrix.
The impact of using the adaptive multiplication inflation method is discussed in the GOSAT synthetic observation assimilation experiments in Sect. 4.2.
2.2 MIROC4ACTM
Model for Interdisciplinary Research on Climate, version 4.0 (MIROC4)based ACTM (hereafter referred to as MIROC4ACTM) (Patra et al., 2018; Bisht et al., 2021) is used here for CH_{4} concentration simulations. The model simulations have been performed at a horizontal grid resolution of approximately 2.8 × 2.8^{∘} latitude–longitude (T42 spectral truncations) and at hybrid vertical coordinates of 67 levels (Earth's surface to 0.0128 hPa; Watanabe et al., 2008). Bisht et al. (2021) performed multitracer analysis and demonstrated the importance of very well resolved stratosphere in the MIROC4ACTM that illustrates better extratropical stratospheric variabilities and simulated tropospheric dynamical fields. The meteorological fields in MIROC4ACTM are nudged to the JMA reanalysis (JRA55) data (Kobayashi et al., 2015).
3.1 Construction of known surface emissions (truth)
Present OSSEs intend to develop basic tuning strategies before the actual data to be assimilated, which is useful to accelerate the operational use of real observations. The OSSE has been discussed here by exploiting the known “truth”. The synthetic observations to be assimilated in the OSSE are generated from nature runs which use bottomup surface emission (true) data to simulate global 3D CH_{4} concentrations. The true surface CH_{4} emissions are prepared on the monthly scale using anthropogenic and natural sectors, minus the surface sinks due to bacterial consumption in the soil (Chandra et al., 2021). The anthropogenic emissions were obtained from the Emission Database for Global Atmospheric Research, version 4.3.2 inventory (EDGARv4.3.2) (JanssensMaenhout et al., 2019), which includes the emissions from the major sectors, such as fugitive sources, enteric fermentation and manure management, and solid waste and wastewater handling. The biomass burning emissions are taken from the Global Fire Database (GFEDv4s) (van der Werf et al., 2017) and Goddard Institute for Space Studies emissions (Fung et al., 1991). The wetland and rice emissions are taken from the processbased model of the terrestrial biogeochemical cycle, Vegetation Integrated Simulator of Trace gases (VISIT) (Ito, 2019), which is based on Cao et al. (1996). Other natural emissions, such as those from the ocean, termites, and mud volcanoes are taken from the TransComCH_{4} intercomparison experiment (Patra et al., 2011). The total emissions are taken as the truth for the OSSEs, and the concentration simulated by MIROC4ACTM will be referred to as synthetic observations.
3.2 Prior flux preparation and LETKF setting
Based on our understanding of CH_{4} inverse modeling, the uncertainty in regional flux estimation is found to be 30 % or lower (Chandra et al., 2021). Therefore, we attempted to reproduce the true flux by starting with a prior flux that is lower than the true flux by 30 % (prior flux has the same seasonal cycles as true flux). The MIROC4ACTM is initialized with a spinup of 3 years (2007–2009) with prior flux distribution. The initial CH_{4} distribution on 1 January 2007 was taken from an earlier simulation of 27 years. An initial perturbation with standard deviation of approximately 6 %–8 % spread is applied to the a priori flux as the initial ensemble spread, whereas no ensemble perturbation was applied to the initial CH_{4} concentration. The sensitivity of the initial ensemble spread to CH_{4} flux estimation is discussed in Sect. 4.2. The uncertainty to perturb prior fluxes is generated based on random positive values with normal distribution. The monthly scale prior emission is linearly interpolated at 6hourly intervals to be used in the MIROC4ACTM simulation for data assimilation. This study performs two LETKF data assimilation experiments. In these experiments, we provided initial perturbation on a regional basis over land (53 different land regions; Chandra et al., 2021), and at every grid over ocean, no spatial error correlation between grid points is considered among ensemble members. However, in Sect. 4.2.5, we also discussed the sensitivity of CH_{4} data assimilation by providing the initial ensemble spread at every grid by considering the horizontal spatial error correlation between grid points among ensemble members, with a global mean correlation of 20 %.
3.3 Experiment 1: synthetic dense observation formulation
The OSSE setting with very accurate and dense observation surface data is an attempt to demonstrate that the data assimilation system works reasonably in the estimation of the true surface flux. Errors in the estimated flux could arise due to the insufficient ensemble size and also the implemented inflation methods to overcome the undersampling, along with a simplified forecast process of the emissions. In real data assimilation, there are additional sources of potential errors, such as atmospheric transport and inappropriate prior or observation uncertainties. In our OSSEs, CH_{4} fluxes as mentioned in Sect. 3.2 are used as “true” fluxes in generating synthetic observations (CH_{4} concentrations). In Experiment 1, the simulated surface layer CH_{4} concentrations at each grid for the entire globe were used as synthetic observations. We added a constant measurement uncertainty of 5 ppb, which is typically achieved by the presentday measurement systems (e.g., Lan et al., 2022).
In this study, the CH_{4} observations are assimilated by applying the observation error covariance localization (Kotsuki et al., 2020) to reduce the spurious spatial correlation due to a smaller ensemble size than the degrees of freedom of the system ($\mathbf{R}\leftarrow \mathbf{R}\times \mathrm{exp}(\frac{\mathrm{1}}{\mathrm{2}}\mathit{\{}{\left({d}_{\mathrm{h}}/{\mathit{\sigma}}_{\mathrm{h}}\right)}^{\mathrm{2}}+{\left({d}_{\mathrm{v}}/{\mathit{\sigma}}_{\mathrm{v}}\right)}^{\mathrm{2}}\mathit{\left\}}\right)$), where d_{h} and d_{v} denote the horizontal distance (km) and vertical difference (log[Pa]) between the analysis model grid point and observation location. The tunable parameters σ_{h} and σ_{v} are the horizontal localization scale (km) and vertical localization scale (log[Pa]), respectively. Using the spatial localization technique, we have estimated the CH_{4} flux for each grid by choosing the CH_{4} observations that influence the grid point using an optimal cutoff radius ($\simeq \mathrm{3.65}{\mathit{\sigma}}_{\mathrm{h},\phantom{\rule{0.125em}{0ex}}\mathrm{v}}$; Miyoshi et al., 2007) with a horizontal covariance localization (σ_{h}) of 2200 km and a vertical covariance localization (σ_{v}) of 0.3 in the natural logarithmic pressure (log[Pa]) coordinate. The localization is performed to improve the signaltonoise ratio of ensemblebased covariance. Numerous sensitivity experiments have been performed by varying the horizontal and vertical localization length in order to obtain the optimized CH_{4} flux that best compares with the truth. The LETKF assimilates the observations within the specified radius to solve the analysis state at each grid point independently (Liu et al., 2016; Kotsuki et al., 2020). The state vector of the analysis includes the atmospheric CH_{4} concentration, which is the prognostic variable of forecast model, and the state vector is further augmented by the surface CH_{4} flux, which is not a model prognostic variable. This augmentation enables the LETKF to directly estimate the parameter through the background error covariance with observed variables (Baek et al., 2006). The state vector augmentation is implemented similar to that used by Miyazaki et al. (2011). This approach analyzes CH_{4} flux during the analysis step. The purpose of the simultaneous CH_{4} emission and concentration optimization is to reduce the uncertainty of the initial CH_{4} concentrations on the CH_{4} evolution during the assimilation window and to maximize the observations potential (Tian et al., 2014).
The atmospheric CH_{4} concentration is changed during both the analysis and forecast steps. A challenge of this scheme is that the analysis increment is added to the model state at each analysis step, without considering the global total CH_{4} mass conservation in the model but consistent with the observed local CH_{4} abundance.
In this case, the surface flux at every model grid point is analyzed with an 8 d assimilation window during the year 2010 with 100 ensemble members. The ensemble size and assimilation window are chosen based on the CH_{4} flux estimation accuracy calculated by performing a sensitivity experiment for the ensemble size (60, 80, and 100) and assimilation window (3 and 8 d), respectively (not shown).
3.4 Experiment 2: synthetic satellite observation formulation
One way to address the realworld CH_{4} flux estimation problem is to first make the OSSE dataset like real observations. In this OSSE, we have assimilated synthetic columnaveraged CH_{4} concentrations with a coverage mimicking GOSAT satellite observations. We prepared a modelsimulated columnaveraged CH_{4} concentration (XCH_{4}) dataset that is spatiotemporally sampled with GOSAT observations as follows:
where XCH_{4} is the columnaveraged modelsimulated CH_{4} concentration. XCH_{4(a priori)} is a priori columnaveraged concentration. CH_{4(ACTM)} and CH_{4(a priori)} are the CH_{4} profile from ACTM and a priori, respectively. h_{j} is the pressure weighting function (j is the vertical layer index), and a_{j} represents the averaging kernel matrix for the column retrieval, which is the sensitivity of the retrieved total column at the various (“j”) atmospheric levels. In the next step, we added the same retrieval (XCH_{4}) error as GOSAT to the XCH_{4} (ACTMsimulated) to make the OSSE more realistic and then attempt to estimate the true fluxes.
In this case, the CH_{4} flux has been estimated for each grid by choosing the CH_{4} observation with a cutoff radius ($\simeq \mathrm{3.65}\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{h},\mathrm{v}}$), with a horizontal covariance localization (σ_{h}) of 5000 km and a vertical covariance localization (σ_{v}) of 0.35 in the natural logarithmic pressure (log[Pa]) coordinate. The optimal horizontal and vertical covariance localization values are chosen based on a trialanderror method (those with the best fits to estimate the CH_{4} flux when compared with truth). A long cutoff radius has been chosen due to sparse observational coverage of GOSAT. Covariance localization is necessary to remove longrange erroneous correlations and for mitigating sampling errors in the ensemblebased error covariance with a limited ensemble size (Miyoshi et al., 2007; Greybush et al., 2011; Kotsuki et al., 2020). The surface flux is analyzed at every model grid point with an 8 d assimilation window and 100 ensemble members; they are chosen based on the sensitivity experiments discussed in Sect. 4.2.
4.1 Experiment with dense OSSEs
The time series of normalized root mean square error (RMSE; $\sqrt{{\sum}_{i=\mathrm{1}}^{n}({x}_{i}^{\mathrm{a}}{x}_{i}^{\mathrm{t}}{)}^{\mathrm{2}}/n}/{\stackrel{\mathrm{\u0303}}{x}}^{\mathrm{t}}$, where ${x}_{i}^{\mathrm{a}}$ and ${x}_{i}^{\mathrm{t}}$ are the analysis and true state at ith model grid point, n is the total number of grid points, and ${\stackrel{\mathrm{\u0303}}{x}}^{\mathrm{t}}$ represents the mean of true flux) in the analyses over the global landmass region is shown in Fig. 2. The normalized global RMSE is calculated using FM and RTPS inflation methods (Fig. 2) after assimilating synthetic observation at every grid (Sect. 3.4). Noteworthy is that the experiment with the FM inflation method shows 10 %–15 % larger error in estimating the atmospheric surface CH_{4} flux compared to the RTPS inflation method. One of the reasons of the better RMSE using the RTPS inflation method is the higher number of degrees of freedom provided by relaxation (α_{RTPS}) in the ensemble spread (Eq. 8) that could nudge the ensemble of CH_{4} concentrations towards observations. The initial flux analysis spread using RTPS and FM is shown in the Supplement (Fig. S1) and shows larger initial analysis flux spread over Brazil, tropical America, and Asia in RTPS inflation compared to the FM inflation method. We performed numerous sensitivity tests with the RTPS inflation method and found that uniform relaxation is not substantial for some of the regions. Figure 2 shows the RMSE for FM, fixed RTPS (α_{RTPS}=0.4; applied globally, the optimized value is obtained by manual finetuning), and conditional RTPS (α_{RTPS}=0.3–0.7 applied different α_{RTPS} values regionally by manual finetuning). In the case of conditional RTPS, the optimal values of α_{RTPS}, i.e., 0.6, 0.3, and 0.7 for the regions south of 20^{∘} S, 20^{∘} S–20^{∘} N, and north of 20^{∘} N, respectively, were obtained from data assimilation sensitivity calculations with varying α_{RTPS} values for the three regions separately to best match the true states. We find that the conditional RTPS method improves the accuracy by ∼ 5 % compared to fixed RTPS and 10 %–15 % compared to FM. In the following, we discuss the results obtained using the conditional RTPS and FM inflation methods.
We have also shown the RMSE (not normalized) of the surface flux in the Supplement (Fig. S2). The flux RMSE has been estimated globally for both the inflation methods and also for the region south of 20^{∘} N (by considering only those land grids which fall in the region south of 20^{∘} N; Fig. S2) for comparative purposes. It was noticed that (Fig. S2 in the Supplement), above north of 20^{∘} N, the flux estimation error is higher, specifically during spring–summer when CH_{4} emissions peak over most of the northern hemispheric regions (Fig. 3). The high uncertainty during spring–summer (Fig. S2) in the flux estimation over these regions could appear due to the attenuation of surface observations as a result of active vertical mixing. The RMSE during autumn (Fig. S2) is comparable in the case of the global region and the region south of 20^{∘} N, which indicates that the RMSE is arising from southern hemispheric regions, likely over Brazil, as it peaks during autumn (Fig. 3).
Figure 3 shows a regional total flux seasonal cycle comparison of the estimated fluxes for 15 terrestrial regions with the cycles of the prior and true fluxes. The estimated flux retrieved using RTPS inflation method over different regions agrees well with that of the true flux. We intend to show the capability of LETKFestimated fluxes over these regions using surface observations to mimic the true fluxes in our understanding of the terrestrial biosphere CH_{4} cycle. These results are consistent with Fig. 2, with an annual global normalized mean bias (${\sum}_{i=\mathrm{1}}^{n}({x}_{i}^{\mathrm{a}}{x}_{i}^{\mathrm{t}})/{\sum}_{i=\mathrm{1}}^{n}\left({x}_{i}^{\mathrm{t}}\right)$) of −0.04. It can also be noticed from Fig. 3 that estimated fluxes converge to true fluxes over most of the regions after about 2–3 months.
To see the degree of similarity in the flux distribution between the estimated and true fluxes, we show monthly mean spatial flux distribution for June and November in Figs. 4 and 5, respectively, along with the bias in the prior and estimated flux. As shown in Figs. 4 and 5, the general spatial patterns of the true flux are estimated well. These results suggest that our LETKF system is capable of reproducing continental spatial flux patterns by using such idealized dense surface observational data. However, some clear differences in flux estimation could be noticed from the FM and RTPS inflation method (Figs. 4 and 5); e.g., over the Eurasian and American continent, analysis with RTPS shows clear improvement compared to the FM covariance inflation method. We calculated the global mean normalized bias with the RTPS and FM covariance inflation method, which is found to be −0.04 and −0.11, respectively, over land regions, and this showed that RTPS significantly improved the flux estimation compared to the FM covariance inflation method.
4.2 Experiment by mimicking the real satellite observational dataset
In this section we discuss the LETKF flux estimation by assimilation of GOSAT synthetic CH_{4} concentration observations. Figure 6 shows the modelsimulated mean XCH_{4} concentration sampled spatiotemporally with GOSAT observations during January and July for the year 2010 (sampling method discussed in Sect. 3.4). In this case we have shown different LETKF sensitivity experiments, such as LETKF sensitivity to (1) FM, RTPS, and adaptive multiplicative inflation; (2) the assimilation window; (3) the ensemble size; (4) the χ^{2} test; and (5) the prior ensemble spread. In the LETKF sensitivity experiments from 1–4, the initial ensemble spread employed a similar method to Experiment 1, and conditional RTPS inflation method is used. A conditional RTPS method is also used in Sect. 4.2.6 for CH_{4} flux estimation.
4.2.1 LETKF sensitivity to FM, RTPS, and adaptive multiplicative inflation
This study mainly emphasizes FM and RTPS inflation methods used in CH_{4} LETKF data assimilation. The annual average normalized RMSE (absolute bias) with RTPS and FM covariance inflation is found to be 0.59 (0.18) and 0.64 (0.22), respectively. The RTPS inflation method performs better than the FM inflation method overall. In addition to RTPS inflation, a sensitivity test is also performed using an adaptive multiplicative inflation method.
In the adaptive inflation, we need to provide an initial multiplicative inflation factor at the beginning of data assimilation cycle (Cycle 1 in Fig. 1). Following the method of Miyoshi (2011), the multiplication inflation factor information calculated in the previous cycle (i.e., Cycle 1 in Fig. 1) is used for the next data assimilation cycle at every grid point (Cycle 2 in Fig. 1). We perform two sensitivity experiments. In the first (second) case, we provided 50 % (40 %) initial inflation in the beginning of Cycle 1 (Fig. 1). The normalized RMSE in the both the adaptive inflation sensitivity experiments is comparable (0.65, Supplement Fig. S3a) till July, but from the beginning of August, the RMSE increases exponentially in the first experiment. However, in terms of the χ^{2} distribution, CH_{4} flux estimation with the first sensitivity adaptive multiplicative inflation experiment (50 % initial inflation case) is better than with the second sensitivity experiment (Supplement Fig. S3b; χ^{2} test described in Sect. 4.2.4). To identify the regions of high estimated CH_{4} flux error, we have shown the background error spread in CH_{4} flux estimation over 15 regions (Supplement Fig. S3c) and found that the spread over west and southeast Asia rises exponentially postJuly, which indicates the rise of estimated CH_{4} flux error over these regions in the first sensitivity adaptive multiplicative inflation experiment. Our analysis suggests that CH_{4} flux estimation depends on the initial inflation factor provided in the beginning of the data assimilation cycle (Cycle 1, Fig. 1) in the adaptive multiplication method. Also, we need to be very careful to monitor the background error spread evolution with time to estimate the CH_{4} flux with adaptive inflation; the χ^{2} distribution analysis is not sufficient.
In the case of RTPP inflation, we found the parameter α_{RTPP} is very difficult to finetune due to its very high sensitivity to estimating the CH_{4} flux. We fail to obtain an optimized α_{RTPP} value to estimate the CH_{4} flux. Whitaker and Hamill (2012) also demonstrated the better accuracy of the LETKF meteorological data assimilation with RTPS compared to the RTPP covariance inflation method. They found the RTPP method produces very large errors if the inflation parameter exceeds the optimal value.
4.2.2 Assimilation window
The LETKF data assimilation window length determines the time span of the observations assimilated in each assimilation cycle. We have shown the sensitivity of two assimilation window size configurations, 3 and 8 d, in the Supplement Fig. S4. Our sensitivity experiments with window size configurations show that the 8 d long assimilation window estimates the CH_{4} flux with better accuracy (∼ 10 %) compared to the 3 d assimilation window because more observational information is incorporated into the system with the 8 d long assimilation window. This study uses an 8 d assimilation window for CH_{4} LETKF data assimilation.
4.2.3 Ensemble size
Figure 7a shows the RMSE using different ensemble members. The RMSE stabilizes gradually as the ensemble size increases from 60 to 80 to 100 ensemble members. The ensemble size dependency of flux estimation suggests the further scope of the improvement in flux estimation by increasing the ensemble members. In this study we stick to 100 ensemble members due to high computational cost while solving large covariance matrices. The larger error in flux estimation in the case of columnaveraged synthetic GOSAT CH_{4} observations assimilation compared to dense observations (Fig. 2) is likely due to the weaker constraint on surface fluxes provided by satellite observations and sparse observations.
4.2.4 χ^{2} test
We have carried out a χ^{2} test for the evaluation of background error covariance matrix (Miyazaki et al., 2012). For the χ^{2} test, the innovation statistics are diagnosed from the observation minus forecast (y^{o}−Hx^{b}), the estimated error covariance in the observation space (HP^{b}H^{T}+R), and the number of observations k as
Using this statistic, the χ^{2} is defined as follows:
The performance of the background error covariance matrix is determined based on the high and lower value of χ^{2}. The χ^{2} value should converge to 1; a value higher (lower) than 1 indicates underestimation (overestimation) of the background error covariance matrices. Our results suggest that the background error covariance matrix is highly underestimated in both RTPS and FM covariance inflation methods (Fig. 7b). However, the χ^{2} values' convergence towards 1 is better in the case of RTPS compared to the FM covariance inflation method, which indicates the improved representation of background errors and then more appropriate data assimilation corrections in the case of the RTPS inflation method. The χ^{2} distribution starts saturating after the month of March. PostMarch analysis shows the background error covariance matrix underestimation is much higher (>100 %) in the case of FM compared to the RTPS covariance inflation method.
4.2.5 CH_{4} LETKF sensitivity to the initial ensemble spread
A test case for CH_{4} LETKF data assimilation has been performed, where the initial spread is provided by considering the initial perturbation on each model grid with spatial error correlation between grid points among ensemble members, with a global mean correlation of 20 %. In this case, we found that the analysis fluxes are extremely sensitive to the initial ensemble spread if prior fluxes are perturbed with more than 5 % prior uncertainty. Therefore, we used initial ensemble perturbation with only 2 % prior uncertainty. Reducing the initial ensemble spread reduces the CH_{4} flux estimation sensitivity (>60 %). However, it also poses a challenge to mitigate the underdispersed background error covariance matrix. We performed LETKF data assimilations in this case with the RTPS covariance inflation method (α_{RTPS}=0.9 optimized value is used here uniformly) with an 8 d long assimilation window and 100 ensemble members and calculated the normalized RMSE between the analysis and true fluxes (Supplement Fig. S5). It is noteworthy that the estimated error between the analysis and true fluxes (Fig. S5) with this setting (gridwise initial ensemble spread) is still larger (25 %) than the case when the regionwise initial ensemble spread is used (Fig. 7a; 100 ensemble size). It suggests that initial ensemble spreads among ensemble members need to be meticulously chosen so that they best represent CH_{4} variability among ensembles to estimate the CH_{4} flux.
Note that the OSSEs used in this study did not consider the effects of model errors other than CH_{4} fluxes, such as model transport errors. In real situations, model errors can have a substantial impact on flux estimates (Locatelli et al., 2013), which needs to be taken into account in background covariances. Therefore, the optimal data assimilation setting can differ between the OSSEs presented in this study and real observation cases. Further efforts, e.g., by conducting a more comprehensive OSSE that accounts for various model errors and by performing various sensitivity calculations in real cases, would provide an improved understanding of the optimal inflation settings to improve CH_{4} flux estimates in following study.
4.2.6 Estimated CH_{4} flux analysis
Figure 8 shows the regional flux seasonal cycle comparison for the estimated fluxes over 15 terrestrial regions with the cycles of the prior and true fluxes. We have also shown assimilation results in the case of the FM inflation method in the Supplement (Fig. S6), which shows the flux estimation disagreement over more regions compared to the RTPS inflation method, e.g., for tropical and North America, the whole African continent, and Australia–New Zealand.
We have shown the GOSAT observations in Figs. 6 and S7. We found very marginal flux estimation improvement over Central Africa after May (Fig. 8), which could be associated with the lower GOSAT coverage over this region (Fig. 6). On the other hand, over northern Africa, no improvement in flux estimation is found. In the case of dense OSSEs too (Fig. 3), we did not find satisfactory flux estimation over northern Africa, which is most probably related to the insufficient initial spread among ensemble members over this region (we used the same initial ensemble spread in both OSSE cases). Over Europe, GOSAT observations are remarkably fewer, specifically for the first few months (January–April; Supplement Fig. S7). Therefore, the flux update over Europe would be influenced by the observations from neighboring regions falling under the chosen cutoff radius that are mainly in northern Africa, where the flux estimation itself not satisfactory. It could also be noticed that the retrieval error added in this OSSE case is high over Europe (September–October; Supplement Fig. S7) and its adjacent sea (Mediterranean Sea; June–August), which could also affect the surface CH_{4} flux estimation.
Figures 9 and 10 show spatial patterns of the true and estimated fluxes by assimilating the columnaveraged CH_{4} concentrations during June and November (Fig. 6). It may be noticed that the RTPS covariance inflation method is more able to estimate the true flux pattern compared to the FM covariance inflation method. The spatial pattern shown using the RTPS inflation method emphasizes the positive and negative bias in the estimated flux (Figs. 9 and 10) but generally agrees with the flux seasonal cycle plots shown in Fig. 8.
Our LETKF CH_{4} data assimilation experiment by assimilating GOSAT synthetic observation with the implementation of the advanced RTPS covariance inflation method better estimates the timeevolving surface CH_{4} fluxes compared to the FM covariance inflation method. The difficulty to estimate the surface CH_{4} flux over a few regions may be overcome by applying additional methodologies, such as the assimilation of surface observations simultaneously and the use of information about the CH_{4} flux climatology. A correction factor derived based on empirical formulation that could use CH_{4} flux climatology information is needed to apply to maintain the CH_{4} mass conservation. This could be implemented by checking the simulated CH_{4} burden gain between years in comparison with the observed CH_{4} growth rates.
In this study, we have introduced a 4DLETKF data assimilation system that utilizes MIROC4ACTM as a forward model for CH_{4} flux estimation. This study has extensively tested both FM and RTPS inflation methods for the LETKF CH_{4} flux estimation. We have conducted two experiments to demonstrate the ability of LETKF system to estimate the CH_{4} surface flux globally. In Experiment 1, we have assimilated the synthetic dense surface CH_{4} observations, while in Experiment 2, synthetic GOSAT CH_{4} observations are assimilated. Based on the results of the sensitivity tests using FM and RTPS inflation methods in Experiment 1, we have found that RTPS inflation produces significantly less normalized RMSE (10 %–15 %) compared to the FM inflation method. In Experiment 2, we discussed LETKF parameters, such as different inflation techniques, ensemble size, assimilation window, initial ensemble spread sensitivity, and χ^{2} test. The ensemble size (this study uses maximum 100 ensemble members) sensitivity test suggests that more ensemble members could help to accurately represent the covariance matrix with a higher number of degrees of freedom. The assimilation window sensitivity test shows that an 8 d assimilation window reduces the normalized flux RMSE by about 10 % compared to a 3 d assimilation window in the case of GOSAT synthetic observations assimilation.
Our approach of assimilation with RTPS inflation could provide a higher number of degrees of freedom to fit the ensemble of CH_{4} concentrations to the observed ones, resulting in the improved analyzed fluxes. The RTPS inflation method is capable of obtaining reasonable flux estimates with a normalized annual mean bias of 0.04 and 0.61 in the case of dense surface synthetic observations and GOSAT synthetic observations, respectively. We demonstrated in our sensitivity OSSE with synthetic GOSAT observations that, over American and African continents and also over Australia–New Zealand, the LETKF data assimilation with the FM inflation method does not show much improvement in the true flux estimation, but the RTPS inflation method reasonably estimates the true flux over most of these regions. One of the reasons for better flux estimates with the RTPS inflation method is the drastic prevention of analysis spread. In the CH_{4} LETKF flux estimation, the surface CH_{4} flux is not a prognostic state vector in the ACTM, which results in the continuous decay of spread in analysis steps. The RTPS inflation method could mitigate such an underdispersed spread problem. This study finds that spatially homogeneous relaxation is not sufficient. It needs to be finetuned and applied conditionally.
The sensitivity of LETKF CH_{4} flux estimation to the initial ensemble spread needs to be carefully dealt with when applied to real data assimilation system. A future OSSE with an additive covariance inflation technique could be interesting while applied with the RTPS inflation method for CH_{4} LETKF data assimilation since in additive covariance inflation, initial estimated flux error cannot propagate. The state vector augmentation technique used here updates the flux after each data assimilation cycle, but it does not conserve the total atmospheric CH_{4} amount, which is one of the limitations of this work. A correction factor needs to be implemented to conserve the total atmospheric CH_{4} amount after completion of a few data assimilation cycles. We have not accounted for the transport error due to meteorological fields in this work (Patra et al., 2011); in the case of real observation data assimilation, a weeklong window may introduce transport errors in CH_{4} analysis because of the nonlinear growth of ensemble perturbations.
The LETKF source codes can be accessed from https://doi.org/10.5281/zenodo.7127658 (Bisht et al., 2022a). All the scripts for running the LETKF data assimilation software and the input and output result data files are available at https://doi.org/10.5281/zenodo.7098323 (Bisht et al., 2022b). The CH_{4} ACTM simulation module coupled with MIROC4AGCM can be accessed from https://doi.org/10.5281/zenodo.7118365 (Bisht et al., 2022c). The source code of MIROC4AGCM is archived at https://doi.org/10.5281/zenodo.7274240 (Patra et al., 2022) with restriction because of the copyright policy of the MIROC developer community. This work did not contribute to the MIROC4 source code development.
The supplement related to this article is available online at: https://doi.org/10.5194/gmd1618232023supplement.
The LETKF data assimilation experiments were designed by JSHB. PKP, MT, and TS helped to set up the LETKF code on MIROC4ACTM for CH_{4} data assimilation. The manuscript was prepared by JSHB, and analysis interpretation input and feedback were provided by PKP, TS, and KM. All coauthors, KM, TS, PKP, NS, MT, and YK, contributed to the writing and revision of the paper.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We acknowledge Ryu Saito for the initial setup of the LETKF code on MIROC4ACTM for CH_{4}. We also thank Takemasa Miyoshi for his LETKF scheme, which served as the basis for the development of our data assimilation system.
This research has been supported by the Environment Research and Technology Development Fund (grant no. JPMEERF20182002) of the Environmental Restoration and Conservation Agency of Japan and the GOSATGW project fund.
This paper was edited by ShuChih Yang and reviewed by two anonymous referees.
Anderson, J. L. and Anderson, S. L.: A Monte Carlo Implementation of the Nonlinear Filtering Problem to Produce Ensemble Assimilations and Forecasts, Mon. Weather Rev., 127, 2741–2758, https://doi.org/10.1175/15200493(1999)127<2741:AMCIOT>2.0.CO;2, 1999.
Baek, S.J., Hunt, B. R., Kalnay, E., Ott, E., and Szunyogh, I.: Local ensemble Kalman filtering in the presence of model bias, Tellus A, 58, 293–306, https://doi.org/10.1111/j.16000870.2006.00178.x, 2006.
Bisht, J. S. H., Machida, T., Chandra, N., Tsuboi, K., Patra, P. K., Umezawa, T., Niwa, Y., Sawa, Y., Morimoto, S., Nakazawa, T., Saitoh, N., and Takigawa, M.: Seasonal Variations of SF_{6} , CO_{2} , CH_{4}, and N_{2}O in the UT/LS Region due to Emissions, Transport, and Chemistry, J. Geophys. Res.Atmos., 126, e2020JD033541, https://doi.org/10.1029/2020JD033541, 2021.
Bisht, J. S. H., Patra, P. K., Sekiya, T., and Miyazaki, K.: LETKF: CH4 data assimilation code, including the wrapper script to run the assimilation system, Zenodo [code], https://doi.org/10.5281/zenodo.7127658, 2022a.
Bisht, J. S. H., Patra, P. K., Takigawa, M., Sekiya, T., Kanaya, Y., Saitoh, N., and Miyazaki, K.: MIROC4ACTM: Model setup, input and output data for CH4 LETKF (Bisht et al., GMDD, 2022), Zenodo [data set], https://doi.org/10.5281/zenodo.7098323, 2022b.
Bisht, J. S. H., Patra, P. K., and Takigawa, M.: MIROC4ACTM code for CH4 simulation, Zenodo [code], https://doi.org/10.5281/zenodo.7118365, 2022c.
Bruhwiler, L., Dlugokencky, E., Masarie, K., Ishizawa, M., Andrews, A., Miller, J., Sweeney, C., Tans, P., and Worthy, D.: CarbonTrackerCH_{4}: an assimilation system for estimating emissions of atmospheric methane, Atmos. Chem. Phys., 14, 8269–8293, https://doi.org/10.5194/acp1482692014, 2014.
Cao, M., Marshall, S., and Gregson, K.: Global carbon exchange and methane emissions from natural wetlands: Application of a processbased model, J. Geophys. Res.Atmos., 101, 14399–14414, https://doi.org/10.1029/96JD00219, 1996.
Chandra, N., Patra, P. K., Bisht, J. S. H., Ito, A., Umezawa, T., Saigusa, N., Morimoto, S., Aoki, S., JanssensMaenhout, G., Fujita, R., Takigawa, M., Watanabe, S., Saitoh, N., and Canadell, J. G.: Emissions from the Oil and Gas Sectors, Coal Mining and Ruminant Farming Drive Methane Growth over the Past Three Decades, J. Meteorol. Soc. Japan. Ser. II, 99, 309–337, https://doi.org/10.2151/jmsj.2021015, 2021.
Daley, R.: The laggedinnovation covariance: A performance diagnostic for data assimilation, Mon. Weather Rev., 120, 178–196, 1992.
Desroziers, G., Berre, L., Chapnik, B., and Poli, P.: Diagnosis of observation, background and analysiserror statistics in observation space, Q. J. Roy. Meteor. Soc., 131, 3385–3396, https://doi.org/10.1256/qj.05.108, 2005.
Enting, I. G.: Inverse Problems in Atmospheric Constituent Transport, Cambridge University Press, Online ISBN 9780511535741, https://doi.org/10.1017/CBO9780511535741, 2002.
Fung, I., John, J., Lerner, J., Matthews, E., Prather, M., Steele, L. P., and Fraser, P. J.: Threedimensional model synthesis of the global methane cycle, J. Geophys. Res., 96, 13033, https://doi.org/10.1029/91JD01247, 1991.
Greybush, S. J., Kalnay, E., Miyoshi, T., Ide, K., and Hunt, B. R.: Balance and Ensemble Kalman Filter Localization Techniques, Mon. Weather Rev., 139, 511–522, https://doi.org/10.1175/2010MWR3328.1, 2011.
Houtekamer, P. L. and Zhang, F.: Review of the Ensemble Kalman Filter for Atmospheric Data Assimilation, Mon. Weather Rev., 144, 4489–4532, https://doi.org/10.1175/MWRD150440.1, 2016.
Houweling, S., Kaminski, T., Dentener, F., Lelieveld, J., and Heimann, M.: Inverse modeling of methane sources and sinks using the adjoint of a global transport model, J. Geophys. Res.Atmos., 104, 26137–26160, https://doi.org/10.1029/1999JD900428, 1999.
Hunt, B. R., Kostelich, E. J., and Szunyogh, I.: Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter, Phys. D Nonlinear Phenom., 230, 112–126, https://doi.org/10.1016/j.physd.2006.11.008, 2007.
Ito, A.: Methane emission from panArctic natural wetlands estimated using a processbased model, 1901–2016, Polar Sci., 21, 26–36, https://doi.org/10.1016/j.polar.2018.12.001, 2019.
JanssensMaenhout, G., Crippa, M., Guizzardi, D., Muntean, M., Schaaf, E., Dentener, F., Bergamaschi, P., Pagliari, V., Olivier, J. G. J., Peters, J. A. H. W., van Aardenne, J. A., Monni, S., Doering, U., Petrescu, A. M. R., Solazzo, E., and Oreggioni, G. D.: EDGAR v4.3.2 Global Atlas of the three major greenhouse gas emissions for the period 1970–2012, Earth Syst. Sci. Data, 11, 959–1002, https://doi.org/10.5194/essd119592019, 2019.
Kalnay, E. and Yang, S.C.: Accelerating the spinup of Ensemble Kalman Filtering, Q. J. Roy. Meteor. Soc., 136, 1644–1651, https://doi.org/10.1002/qj.652, 2010.
Kang, J.S., Kalnay, E., Miyoshi, T., Liu, J., and Fung, I.: Estimation of surface carbon fluxes with an advanced data assimilation methodology, J. Geophys. Res.Atmos., 117, D24101, https://doi.org/10.1029/2012JD018259, 2012.
Kobayashi, S., Ota, Y., Harada, Y., Ebita, A., Moriya, M., Onoda, H., Onogi, K., Kamahori, H., Kobayashi, C., Endo, H., Miyaoka, K., and Takahashi, K.: The JRA55 Reanalysis: General Specifications and Basic Characteristics, J. Meteorol. Soc. Japan. Ser. II, 93, 5–48, https://doi.org/10.2151/jmsj.2015001, 2015.
Kotsuki, S., Ota, Y., and Miyoshi, T.: Adaptive covariance relaxation methods for ensemble data assimilation: experiments in the real atmosphere, Q. J. Roy. Meteor. Soc., 143, 2001–2015, https://doi.org/10.1002/qj.3060, 2017.
Kotsuki, S., Pensoneault, A., Okazaki, A., and Miyoshi, T.: Weight structure of the Local Ensemble Transform Kalman Filter: A case with an intermediate atmospheric general circulation model, Q. J. Roy. Meteor. Soc., 146, 3399–3415, https://doi.org/10.1002/qj.3852, 2020.
Lan, X., Dlugokencky, E. J., Mund, J. W., Crotwell, A. M., Crotwell, M. J., Moglia, E., Madronich, M., Neff, D., and Thoning, K. W.: Atmospheric Methane Dry Air Mole Fractions from the NOAA GML Carbon Cycle Cooperative Global Air Sampling Network, 1983–2021, Version: 20221121, NOAA [data set], https://doi.org/10.15138/VNCZM766, 2022.
Liu, J., Bowman, K. W., and Lee, M.: Comparison between the Local Ensemble Transform Kalman Filter (LETKF) and 4DVar in atmospheric CO 2 flux inversion with the Goddard Earth Observing SystemChem model and the observation impact diagnostics from the LETKF, J. Geophys. Res.Atmos., 121, 13066–13087, https://doi.org/10.1002/2016JD025100, 2016.
Locatelli, R., Bousquet, P., Chevallier, F., FortemsCheney, A., Szopa, S., Saunois, M., AgustiPanareda, A., Bergmann, D., Bian, H., CameronSmith, P., Chipperfield, M. P., Gloor, E., Houweling, S., Kawa, S. R., Krol, M., Patra, P. K., Prinn, R. G., Rigby, M., Saito, R., and Wilson, C.: Impact of transport model errors on the global and regional methane emissions estimated by inverse modelling, Atmos. Chem. Phys., 13, 9917–9937, https://doi.org/10.5194/acp1399172013, 2013.
Lorente, A., Borsdorff, T., Butz, A., Hasekamp, O., aan de Brugh, J., Schneider, A., Wu, L., Hase, F., Kivi, R., Wunch, D., Pollard, D. F., Shiomi, K., Deutscher, N. M., Velazco, V. A., Roehl, C. M., Wennberg, P. O., Warneke, T., and Landgraf, J.: Methane retrieved from TROPOMI: improvement of the data product and validation of the first 2 years of measurements, Atmos. Meas. Tech., 14, 665–684, https://doi.org/10.5194/amt146652021, 2021.
Maasakkers, J. D., Jacob, D. J., Sulprizio, M. P., Turner, A. J., Weitz, M., Wirth, T., Hight, C., DeFigueiredo, M., Desai, M., Schmeltz, R., Hockstad, L., Bloom, A. A., Bowman, K. W., Jeong, S., and Fischer, M. L.: Gridded National Inventory of U.S. Methane Emissions, Environ. Sci. Technol., 50, 13123–13133, https://doi.org/10.1021/acs.est.6b02878, 2016.
Meirink, J. F., Bergamaschi, P., and Krol, M. C.: Fourdimensional variational data assimilation for inverse modelling of atmospheric methane emissions: method and comparison with synthesis inversion, Atmos. Chem. Phys., 8, 6341–6353, https://doi.org/10.5194/acp863412008, 2008.
Miyazaki, K., Maki, T., Patra, P., and Nakazawa, T.: Assessing the impact of satellite, aircraft, and surface observations on CO_{2} flux estimation using an ensemblebased 4D data assimilation system, J. Geophys. Res., 116, D16306, https://doi.org/10.1029/2010JD015366, 2011.
Miyazaki, K., Eskes, H. J., and Sudo, K.: Global NO_{x} emission estimates derived from an assimilation of OMI tropospheric NO_{2} columns, Atmos. Chem. Phys., 12, 2263–2288, https://doi.org/10.5194/acp1222632012, 2012.
Miyazaki, K., Sekiya, T., Fu, D., Bowman, K. W., Kulawik, S. S., Sudo, K., Walker, T., Kanaya, Y., Takigawa, M., Ogochi, K., Eskes, H., Boersma, K. F., Thompson, A. M., Gaubert, B., Barre, J., and Emmons, L. K.: Balance of Emission and Dynamical Controls on Ozone During the KoreaUnited States Air Quality Campaign From Multiconstituent Satellite Data Assimilation, J. Geophys. Res.Atmos., 124, 387–413, https://doi.org/10.1029/2018JD028912, 2019.
Miyoshi, T.: The Gaussian Approach to Adaptive Covariance Inflation and Its Implementation with the Local Ensemble Transform Kalman Filter, Mon. Weather Rev., 139, 1519–1535, https://doi.org/10.1175/2010MWR3570.1, 2011.
Miyoshi, T., Yamane, S., and Enomoto, T.: Localizing the Error Covariance by Physical Distances within a Local Ensemble Transform Kalman Filter (LETKF), SOLA, 3, 89–92, https://doi.org/10.2151/sola.2007023, 2007.
Miyoshi, T., Sato, Y., and Kadowaki, T.: Ensemble Kalman Filter and 4DVar Intercomparison with the Japanese Operational Global Analysis and Prediction System, Mon. Weather Rev., 138, 2846–2866, https://doi.org/10.1175/2010MWR3209.1, 2010.
Ott, E., Hunt, B. R., Szunyogh, I., Zimin, A. V., Kostelich, E. J., Corazza, M., Kalnay, E., Patil, D. J., and Yorke, J. A.: A Local Ensemble Kalman Filter for Atmospheric Data Assimilation, arXiv [preprint], https://doi.org/10.48550/arXiv.physics/0203058, 19 March 2002.
Ott, E., Hunt, B. R., Szunyogh, I., Zimin, A. V., Kostelich, E. J., Corazza, M., Kalnay, E., Patil, D. J., and Yorke, J. A.: A local ensemble Kalman filter for atmospheric data assimilation, Tellus A, 56, 415–428, https://doi.org/10.3402/tellusa.v56i5.14462, 2004.
Patra, P. K., Houweling, S., Krol, M., Bousquet, P., Belikov, D., Bergmann, D., Bian, H., CameronSmith, P., Chipperfield, M. P., Corbin, K., FortemsCheiney, A., Fraser, A., Gloor, E., Hess, P., Ito, A., Kawa, S. R., Law, R. M., Loh, Z., Maksyutov, S., Meng, L., Palmer, P. I., Prinn, R. G., Rigby, M., Saito, R., and Wilson, C.: TransCom model simulations of CH_{4} and related species: linking transport, surface flux and chemical loss with CH_{4} variability in the troposphere and lower stratosphere, Atmos. Chem. Phys., 11, 12813–12837, https://doi.org/10.5194/acp11128132011, 2011.
Patra, P. K., Takigawa, M., Watanabe, S., Chandra, N., Ishijima, K., and Yamashita, Y.: Improved Chemical Tracer Simulation by MIROC4.0based Atmospheric ChemistryTransport Model (MIROC4ACTM), SOLA, 14, 91–96, https://doi.org/10.2151/sola.2018016, 2018.
Patra, P. K., Takigawa, M., Watanabe, S., Chandra, N., Ishijima, K., and Yamashita, Y.: Improved Chemical Tracer Simulation by MIROC4.0based Atmospheric ChemistryTransport Model (MIROC4ACTM) (Version 1), Zenodo [code], https://doi.org/10.5281/zenodo.7274240, 2022.
Peters, W., Miller, J. B., Whitaker, J., Denning, A. S., Hirsch, A., Krol, M. C., Zupanski, D., Bruhwiler, L., and Tans, P. P.: An ensemble data assimilation system to estimate CO_{2} surface fluxes from atmospheric trace gas observations, J. Geophys. Res., 110, D24304, https://doi.org/10.1029/2005JD006157, 2005.
Saunois, M., Stavert, A. R., Poulter, B., Bousquet, P., Canadell, J. G., Jackson, R. B., Raymond, P. A., Dlugokencky, E. J., Houweling, S., Patra, P. K., Ciais, P., Arora, V. K., Bastviken, D., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Carlson, K. M., Carrol, M., Castaldi, S., Chandra, N., Crevoisier, C., Crill, P. M., Covey, K., Curry, C. L., Etiope, G., Frankenberg, C., Gedney, N., Hegglin, M. I., HöglundIsaksson, L., Hugelius, G., Ishizawa, M., Ito, A., JanssensMaenhout, G., Jensen, K. M., Joos, F., Kleinen, T., Krummel, P. B., Langenfelds, R. L., Laruelle, G. G., Liu, L., Machida, T., Maksyutov, S., McDonald, K. C., McNorton, J., Miller, P. A., Melton, J. R., Morino, I., Müller, J., MurguiaFlores, F., Naik, V., Niwa, Y., Noce, S., O'Doherty, S., Parker, R. J., Peng, C., Peng, S., Peters, G. P., Prigent, C., Prinn, R., Ramonet, M., Regnier, P., Riley, W. J., Rosentreter, J. A., Segers, A., Simpson, I. J., Shi, H., Smith, S. J., Steele, L. P., Thornton, B. F., Tian, H., Tohjima, Y., Tubiello, F. N., Tsuruta, A., Viovy, N., Voulgarakis, A., Weber, T. S., van Weele, M., van der Werf, G. R., Weiss, R. F., Worthy, D., Wunch, D., Yin, Y., Yoshida, Y., Zhang, W., Zhang, Z., Zhao, Y., Zheng, B., Zhu, Q., Zhu, Q., and Zhuang, Q.: The Global Methane Budget 2000–2017, Earth Syst. Sci. Data, 12, 1561–1623, https://doi.org/10.5194/essd1215612020, 2020.
Sekiya, T., Miyazaki, K., Ogochi, K., Sudo, K., Takigawa, M., Eskes, H., and Boersma, K. F.: Impacts of Horizontal Resolution on Global Data Assimilation of Satellite Measurements for Tropospheric Chemistry Analysis, J. Adv. Model. Earth Syst., 13, e2020MS002180, https://doi.org/10.1029/2020MS002180, 2021.
Skachko, S., Ménard, R., Errera, Q., Christophe, Y., and Chabrillat, S.: EnKF and 4DVar data assimilation with chemical transport model BASCOE (version 05.06), Geosci. Model Dev., 9, 2893–2908, https://doi.org/10.5194/gmd928932016, 2016.
Szopa, S., Naik, V., Adhikary, B., Artaxo, P., Berntsen, T., Collins, W. D., Fuzzi, S., Gallardo, L., KiendlerScharr, A., Klimont, Z., Liao, H., Unger, N., and Zanis, P.: ShortLived Climate Forcers, in: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 817–922, https://doi.org/10.1017/9781009157896.008, 2021.
Tian, X., Xie, Z., Liu, Y., Cai, Z., Fu, Y., Zhang, H., and Feng, L.: A joint data assimilation system (TanTracker) to simultaneously estimate surface CO_{2} fluxes and 3D atmospheric CO_{2} concentrations from observations, Atmos. Chem. Phys., 14, 13281–13293, https://doi.org/10.5194/acp14132812014, 2014.
Watanabe, S., Miura, H., Sekiguchi, M., Nagashima, T., Sudo, K., Emori, S., and Kawamiya, M.: Development of an atmospheric general circulation model for integrated Earth system modeling on the Earth Simulator, J. Earth Simulator, 9, 27–35, 2008.
van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen, Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–2016, Earth Syst. Sci. Data, 9, 697–720, https://doi.org/10.5194/essd96972017, 2017.
Whitaker, J. S. and Hamill, T. M.: Evaluating Methods to Account for System Errors in Ensemble Data Assimilation, Mon. Weather Rev., 140, 3078–3089, https://doi.org/10.1175/MWRD1100276.1, 2012.
Yoshida, Y., Kikuchi, N., Morino, I., Uchino, O., Oshchepkov, S., Bril, A., Saeki, T., Schutgens, N., Toon, G. C., Wunch, D., Roehl, C. M., Wennberg, P. O., Griffith, D. W. T., Deutscher, N. M., Warneke, T., Notholt, J., Robinson, J., Sherlock, V., Connor, B., Rettinger, M., Sussmann, R., Ahonen, P., Heikkinen, P., Kyrö, E., Mendonca, J., Strong, K., Hase, F., Dohe, S., and Yokota, T.: Improvement of the retrieval algorithm for GOSAT SWIR XCO_{2} and XCH_{4} and their validation using TCCON data, Atmos. Meas. Tech., 6, 1533–1547, https://doi.org/10.5194/amt615332013, 2013.
Zhang, F., Snyder, C., and Sun, J.: Impacts of Initial Estimate and Observation Availability on ConvectiveScale Data Assimilation with an Ensemble Kalman Filter, Mon. Weather Rev., 132, 1238–1253, https://doi.org/10.1175/15200493(2004)132<1238:IOIEAO>2.0.CO;2, 2004.
Zhang, Y., Jacob, D. J., Lu, X., Maasakkers, J. D., Scarpelli, T. R., Sheng, J.X., Shen, L., Qu, Z., Sulprizio, M. P., Chang, J., Bloom, A. A., Ma, S., Worden, J., Parker, R. J., and Boesch, H.: Attribution of the accelerating increase in atmospheric methane during 2010–2018 by inverse analysis of GOSAT observations, Atmos. Chem. Phys., 21, 3643–3666, https://doi.org/10.5194/acp2136432021, 2021.