EC&V Pty Ltd: TROPOMI CSF: Introduction

Introduction

Numerous studies have used TROPOMI satellite observations to analyse CH4 emissions. One such study is Sadavarte, et al., 2021 , which concludes that emissions from the Hail Creek open cut coal mine in central Queensland are exceptionally high.

The basic data for the Hail Creek mine reported in Table 1 of Sadavarte, et al., 2021 are reproduced below in Table 1.

Table 1. Extract of Hail Creek data from Sadavarte, et al., 2021 . Download
DetailsHail Creek
mine typesurface
mining typeDragline, truck and shovel
total raw coal production in million tonnes 2018-19: 7.7
2019-20: 5.8
longitude, latitude148.380oE, 21.490oS
number of clear-sky observations in TROPOMI32
annual emissions using the CSF method [μ ± 2σ] 230 ± 50 (Ggyear)

In the Supporting Information (Table S4), the authors report an Implied Emission Factor (IEF) of 34.12 g CH4 per kg of raw coal production, at least an order of magnitude higher than the second highest estimate reported by Kholod, et al., 2020 . For convenience, the data relevant to Hail Creek from the original Table S4 are reproduced below in Table 2.

Table 2. Emission factors (Replicating Table S4 from Sadavarte, et al., 2021 Supporting Information .) Download
SourceReferenceEmission yearIEFComment
Australia EDGAR v4.3.2 Australia NIR CRF 2017 20120.55 Split EDGAR underground and surface based on 2012 ratio NIR 2017. EDGARv4.3.2 CH4 emissions for 2012 and used national raw coal production from common reporting format table for 2012.
Australia Australia NIR CRF 2020 20120.45Table 1.B.1
Australia Australia NIR CRF 2020 20170.47Table 1.B.1
Australia Australia NIR CRF 2020 20180.53Table 1.B.1
QueenslandNIR (2020)20180.73 Using emissions and raw coal production details from common reporting format table for respective year.
Hail Creek Sadavarte, et al., 2021 20180.73Reconstructed bottom-up
Global Kholod, et al., 2020 2.03 - 3.38below 200m
Hail Creek Sadavarte, et al., 2021 34.12CSF

The IEF was calculated by dividing the total estimated methane emissions for 2018-2019 (460 (Gg)) by the combined raw coal production over the same period (7.7 + 5.8 (Mtyear) ), yielding a two year average IEF of 34.07 (gkg). The original paper reported IEF of 31.12 (gkg), the minor difference possibly due to use of the full precision data.

The default IPCC emission factors ( IPCC 2006 Guidelines Volume 2 Chapter 4 , p. 4.18) are expressed as volumetric emissions per unit mass (m3t), with a recommended conversion factor of 0.67 to (gkg).

Table 3. IPCC 2006 Recommended emission factors. IPCC 2006 Guidelines Volume 2 Chapter 4 . Download
(m3t)(gkg)Comment
0.30.20 Low Emission Factor. Value in Sadavarte et all marked as 0 - 200m
1.20.88 Average Emission Factor. Value in Sadavarte et all is slightly higher than conversion and marked as 200-400m
2.01.49 High Emissions Factor. Value in Sadavarte et all is slightly higher than conversion and marked as 400+ m

Between May 2022 and September 2023, aircraft measurements were conducted on four days Borchardt, et al., 2025 , yielding annualised emission estimates approximately half of those reported by Sadavarte, et al., 2021 .

The findings of Sadavarte, et al., 2021 and Borchardt, et al., 2025 have been widely cited in the media, including coverage by the Australian Financial Review (30, Nov, 2021) , Bloomberg (25, Mar, 2025) , and a range of advocacy, research, and policy organisations such as Australian Centre for Corporate Responsibility , Australian Conservation Foundation , Climate Council , Methane Zero , Move Beyond Coal , Queensland Conservation Council , Renew Economy , UNEP , UNSW .

Anticipating the availability of observations from the TROPOMI spectrometer and the MERLIN lidar mission, Jacob, et al., 2016 proposed several approaches for estimating emission rates based on probabilistic inversion frameworks. These include the analytical method, adjoint method, and Markov Chain Monte Carlo (MCMC) method. All approaches assume an initial estimate of the emission field, use a CTM to simulate atmospheric concentrations, and determine the optimal emission rates through error weighted statistical fitting between simulated and observed concentrations. An implementation of the adjoint approach, referred to as the Ensemble Kalman Filter Inverse Method, was applied to estimate emissions from the Hail Creek mine by Palmer, et al., 2021 . More recently, a similar inversion framework has been used by the Superpower Institute to estimate emissions for January-June 2023 as part of the Open Methane (beta) project.

Varon, et al., 2018 discuss several approaches to the plume inversion problem based on instantaneous plume observations, in which each satellite image is analysed independently to estimate emissions. These methods include the Gaussian plume inversion method, the source pixel method, the Cross Sectional Flux method (CSF) , and the Integrated Mass Enhancement method (IME) . All of these approaches require identification of the satellite pixels constituting the plume, as well as an estimate of the background concentration field. The CSF method was applied by Sadavarte, et al., 2021 .

The TROPOMICH4 retrieval algorithm has been extended through the introduction of the Weighting Function Modified Differential Optical Absorption Spectroscopy (WFMD) approach by Schneising, et al., 2023 . Version 1.8 of this retrieval additionally incorporates a de-striping procedure to mitigate instrumental artefacts. More recently, Schneising, et al., 2024 developed an adjacent method that enables the simultaneous analysis of CO and CO2, which has since been applied to quantify emissions from steel plants in Germany.

The TROPOMI WFMD data, was subsequently used by Borchardt, et al., 2025 as a benchmark for comparison with aircraft based measurements over Hail Creek.

In addition, preliminary (beta) emission estimates for January-June 2023 have been produced by the Superpower Institute as part of the Open Methane (beta) project.

The remainder of this Introduction provides an overview of the key concepts required for the subsequent discussion. It begins with a description of the observational datasets, starting with Automatic Weather Station (AWS) measurements from Moranbah Airport. This is followed by an overview of ERA5 reanalysis data obtained from the Copernicus Data Store, which are required inputs for the CSF algorithm, and a brief introduction to the TROPOMI dataset. The concept of the Stable Nocturnal Boundary Layer (SNBL) is then introduced, followed by a short overview of the HYSPLIT model. Finally, the basic equations underlying the CSF, TM and IME methods are derived to clarify the assumptions under which each approach is valid.


© 2026 Copyright EC&V Pty Ltd
ABN 23 114 756 688
Last Updated: