the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Efficient use of a Lagrangian Particle Dispersion Model for atmospheric inversions using satellite observations of column mixing ratios
Abstract. Satellite instruments for measuring atmospheric column mixing ratios have improved significantly over the past couple of decades with increases in pixel resolution and accuracy. As a result, satellite observations are being increasingly used in atmospheric inversions to improve estimates of emissions of greenhouse gases (GHGs), particularly CO2 and CH4, and to constrain regional and national emission budgets. However, in order to make use of the increasing resolution in inversions, the atmospheric transport models used need to be able to represent the observations at these finer resolutions. Here, we present a new and computationally efficient methodology to model satellite column average mixing ratios with a Lagrangian Particle Dispersion Model (LPDM) and calculate the Jacobian matrices describing the relationship between surface fluxes of GHGs and atmospheric column average mixing ratios, as needed in inversions. We present a case study using this methodology in the LMPD, FLEXPART, and the inversion framework, FLEXINVERT, to estimate CH4 fluxes over Siberia using column average mixing ratios of CH4 (XCH4) from the TROPOMI instrument onboard the Sentinel-5P satellite. The results of the inversion using TROPOMI XCH4 are evaluated against results using ground-based observations.
- Preprint
(3432 KB) - Metadata XML
-
Supplement
(4489 KB) - BibTeX
- EndNote
Status: open (until 08 Apr 2025)
-
RC1: 'Comment on egusphere-2025-147', Anonymous Referee #1, 20 Mar 2025
reply
General comments
This study present an innovative method for estimating total column methane (XCH4) using a Lagrangian Particle Dispersion Model (LPDM), FLEXPART, and way to assimilate the data in an atmospheric inverse model, FLEXINVERT. The case study is carried out for Siberia for 2022. The comparison against the “traditional” ground-based inversion showed broad agreement with the inversion using TROPOMI data, and consequently reliability and a good potential in the presented method for estimation of regional CH4 fluxes. The method sounds applicable for other LPDMs, and could contribute significantly to the atmospheric inverse modelling community with new ways to infer greenhouse gas flux information from satellite data. The fact that LPDMs can be run in much higher resolution than Eulerian transport models will be an advantage for incorporating information from future satellites with much higher spatial resolution. Therefore, this paper is worth of prompt publication after considering a few points below.
- The authors found large differences in modelled XCH4 depending on background initial mixing ratios. The boundary conditions were optimised to somewhat discriminate the “errors”, but how can it be sure that the signals within the domain is not over constrained by the background? In L280, it is said that the correction of modelled XCH4 was “largely due to the improvement to the background estimate”, and the posterior fluxes from the TROPOMI inversion did not change much from the prior. I suppose that with less uncertainty in the background, the fluxes would change more. How do you know what is a good balance?
- High northern latitudes are challenging regions such that satellite retrievals are associated with various biases, especially those related to seasonal variations may plan an significant role. I wonder how much of the differences between the ground-based and TROPOMI inversions were due to these biases. Please discuss.
- The inverse model results are associated with various uncertainties. You have discussed and did sensitivity tests on background mixing ratios, but other optimization setups, such as choice of retrieval products, pre-processing methods of the satellite data, prior fluxes and prior uncertainties for observations and fluxes are also important. As this paper do not present variety of sensitivity tests, the conclusion of the paper about CH4 flux estimates should be presented carefully that the results may change significantly depending on the setups.
Specific comments
L160: Why do you use area-weighed averages? I understand it is somewhat reasonable for mixing ratios, but for averaging kernels and presser weighting, I do not fully understand how the area would be affected. Could you also explain how did you take into account differences in number of observations within the aggregated cells?
Section 3: Did you include temporal correlation of the state vectors? Please add information somewhere.
L195: Could you add a figure on spatial resolution? Where were lowest and highest resolutions? Were the resolutions same for the TROPOMI and ground-based inversions, despite the fact that they would have differences in “how strongly the fluxes influence the observations” due to differences in locations and quantity (surface vs total column) of the observations? Please clarify.
L209-213: In later sections, it is said that the background mixing ratios were also optimised. What were the uncertainties in the boundary conditions?
L218: Why did you chose the grid cell sizes of 0.25° and 0.5°? FLEXPART is run at 0.5° and smallest optimisation spatial resolution is also 0.5°, so why did you chose to have observations at higher spatial resolutions? Can FLEXPART resolve differences well (or what is done) if there are more than one observations within a 0.5° x 0.5° grid cell?
L223: I suppose number of observation vary a lot within the study period. Please add information about number of observations also perhaps in 14-days temporal resolution, which is your flux optimisation resolution. I would also like to see for both TROPOMI and ground-based observations.
L225-232:
- The source information for JR-STATIONS are available in Data availability section, but how about other ground-based data?
- Did you process/filter these data at all?
- For the TROPOMI data, you mentioned that the observation uncertainties were 14-20 ppb. How about for these ground-based observations?
Section 3.1.3:
- Did all the prior fluxes had estimates for 2020? If not, what did you do? What were the original resolution of the prior fluxes? Did you do any interpolation when original resolution was lower than 0.5°?
- Do I understand it correctly that you include ocean emissions, but do not optimise them? How large were the contribution of ocean fluxes to the total fluxes of this domain?
Section 3.2: I understand that you only optimise total fluxes, but as you find some spatial differences in the flux increments between the TROPOMI and ground-based inversions, can you speculate whether emissions from oil and gas sources have different seasonal patters in the two inversions?
L286, L306: Could you add uncertainty estimates as well?
L295-299: Are the number of TROPOMI observations less than those from the ground-based stations? Do you argue that number of the observations was persistent for all months? I suppose flux uncertainties are larger in summer (as a whole domain) due to contribution of wetlands (although it is not so clear from Figure 6)? How would the retrieval biases possibly play a role that were discussed in e.g. Lindqvist et al. (2024)?
L318-321: Related to questions above, why do you think that the flux were not as well constrained in the TROPOMI inversions? Satellite data are suppose to have good spatial coverage compared to ground-based data, and with much large number of data, it should, in principal, constrain the fluxes better than the ground-based data. But it is not the case here. Is this a general feature or something specific to high northern latitudes?
Technical comments
Introduction: Please add information about focus/simulation years
L67: Please add references to FLEXPART.
Equations: Please use bold fonts for vectors and matrices.
L210: ...ERA5 at 0.5° x 0.5° and… ?
Table 1. What are the altitudes here? Elevation of the site or height from which FLEXPART trajectories were calculated?
Figure 2: Are these of super-observations? Please clarify.
Figure 4:
- You could perhaps consider adding number of observations here?
- Could you also consider adding ranges?
Please consider combining Figures 5 and 8, and Figures 7 and 9. It would be easy to compare between TROPOMI-based and ground-based inversions that way.
References
Lindqvist, H., et al.: Evaluation of Sentinel-5P TROPOMI Methane Observations at Northern High Latitudes, Remote Sensing, 16, 2979, https://doi.org/10.3390/rs16162979, 2024.
Citation: https://doi.org/10.5194/egusphere-2025-147-RC1
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
114 | 28 | 6 | 148 | 10 | 4 | 4 |
- HTML: 114
- PDF: 28
- XML: 6
- Total: 148
- Supplement: 10
- BibTeX: 4
- EndNote: 4
Viewed (geographical distribution)
Country | # | Views | % |
---|---|---|---|
United States of America | 1 | 50 | 29 |
Sweden | 2 | 17 | 10 |
China | 3 | 12 | 7 |
Germany | 4 | 11 | 6 |
Finland | 5 | 9 | 5 |
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
- 50