the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Continuous weekly monitoring of methane emissions from the Permian Basin by inversion of TROPOMI satellite observations
Daniel J. Varon
Daniel J. Jacob
Benjamin Hmiel
Ritesh Gautam
David R. Lyon
Mark Omara
Melissa Sulprizio
Lu Shen
Drew Pendergrass
Hannah Nesser
Zachary R. Barkley
Natasha L. Miles
Scott J. Richardson
Kenneth J. Davis
Sudhanshu Pandey
Alba Lorente
Tobias Borsdorff
Joannes D. Maasakkers
Ilse Aben
Download
- Final revised paper (published on 11 Jul 2023)
- Preprint (discussion started on 10 Nov 2022)
Interactive discussion
Status: closed
-
RC1: 'Comment on acp-2022-749', Anonymous Referee #1, 24 Jan 2023
The authors have used TROPOMI methane observations to obtain weekly estimates of methane emissions in the Permian Basin at a resolution of 0.25° x 0.3125°. For the period from 1 May 2018 to 5 October 2020, they estimated mean weekly methane emissions of 3.9 ± 1.0 Tg a-1. They also estimated a statistically significant trend in the emissions of -3.5 Gg a-1 per week, and suggested that the trend and variability in the emissions are driven by a combination of several different factors, with the development of new wells and local natural gas prices being the statistically significant contributors. The manuscript is well written and is a great demonstration of the potential of satellite observations from instrument such as TROPOMI for quantifying and monitoring methane emissions. The ability to quantify the variability in methane emissions on a weekly timescale opens a broad range of possible policy applications. I therefore recommend the manuscript for publication after revisions to address my comments below.
Main Comments
- The authors describe the inversion approach as using a Kalman filter, but it is unclear to me in what way is this a Kalman filter. On page 6, lines 156-157, they mention that instead of using a “full Kalman filter” they used a “suboptimal Kalman filter with fixed (diagonal) error covariance matrix.” If the covariance is fixed, and the model dynamics is not used to update the state nor the errors, how is the scheme a Kalman filter? It is unclear to me what is the difference between this scheme and what is referred to as a Bayesian synthesis inversion? I would describe this as a time-dependent Bayesian synthesis inversion.
- Page 3, lines 70-72: The text states that “…starting from best available bottom-up prior estimates of emissions and using Bayesian synthesis to obtain optimized posterior estimates assimilating the information from TROPOMI. We use a Kalman filter to quantify weekly basin-wide emissions.” Based on this description, it sounds as though both the Bayesian synthesis inversion and the Kalman filter schemes are being used. It would be helpful if the authors could better explain this. Perhaps the distinction regarding how the two schemes are used can be added to the schematic in Figure 2?
- Page 5, line 126: The text states that the observation vector assembles the observations for the week. Are the observations ingested sequentially during the week or does the inversion ingest weekly mean observations? This is somewhat unclear.
- Table 2: The discussion in Section 3.5 is focused on comparing the estimated Permian emissions with previous TROPOMI-based estimates (which are given in Table 2). Are the differences in the emission estimates that are compared in Table 2 really meaningful? It seems to me that the errors for these emission estimates all overlap, with the exception of the McNorton et al. (2022) results. Thus, I am not sure how to interpret the discussion in Section 3.5.
- Page 13, lines 309-311: The text here states that the range of reported estimates between this study and the Zhang et al. and Shen et al. studies can be explained by differences in the prior emissions and the background specification. Since this study has produced averaging kernels, the authors could consider substituting the EDF prior emissions with those used in the Zhang et al. and Shen et al. studies to see how much they impact the posterior emissions. The posterior emissions are given by x+ = Ax + (I – A)x-, with the contribution from the prior given by (I – A)x-, where x- is the prior and A is the averaging kernel matrix. How does this contribution change when the other priors are used?
Technical comments
- Figure 4 caption. Change “Figure” to “figure. Also, DOFS was already defined on Page 6, line 171.
- Page 16, line 398: I think the paragraph should start with “Table 2 summarizes…”
Citation: https://doi.org/10.5194/acp-2022-749-RC1 -
RC2: 'Comment on acp-2022-749', Anonymous Referee #2, 24 Feb 2023
The MS[ACP-2022-749] by Varon et al., used satellite-based xCH4 observations and inversion modeling to constrain CH4 emissions from oil and gas production hotspot in Permian Basin, the largest oil production basin in the United States at a weekly scale. It’s very important to quantify CH4 leakage in these fossil fuel production areas. In general, this MS is easy to follow and well written, considering a few similar studies have been conducted in the same region, some clarification is needed to be highlighted, especially for the most important improvement in this study as displayed below. And it can be accepted after addressing the following comments.
Main comments:
In general, there are a few published papers in the same region that used similar approach (i.e. Zhang et al., 2020), although this submission is conducted at a weekly scale instead of a monthly or annual scale, the authors still need to address the improvement of their study from previous ones, i.e. approach in the inversion framework? or found the weekly relationship between CH4 emissions and other activity indexes? or used more robust observations? or prior emissions and background?
Line 19-21 “The mean oil and gas emission from the region (± standard deviation of weekly estimates) was 3.7 ± 0.9 Tg higher than previous TROPOMI inversion estimates that may have used too-low prior emissions or biased background assumptions”. It seems the inversion results are sensitive to prior emissions, have you tested or quantified this potential bias of using different prior emissions to your results?
Line 88-89, “19346 ± 13073 observations per week over our full inversion domain (96°–110°W, 25°–38°N), including 3062 ± 2314”, The standard deviation of available data numbers in each week is so high and is comparable with averages, which indicates there are not enough data in some weeks, please address what the potential bias for emission inversion in these weeks with lower available observation numbers. And can you display the time series of available CH4 observation numbers in each week?
For some rainy or cloudy weeks, the available observation data can be sparse, leading to a large missing data gap in the study domain, and how will this situation affect inversion results for this large region? The reason to mention this comment is that your following analysis of the relationship between CH4 emissions and activity indexes ignored the influence of available data.
Line 92-94, “We use dynamic 3-hour boundary conditions from a global 4°×5° simulation corrected with spatially and temporally smoothed TROPOMI data as described by Shen et al. (2021). A one-month spin-up simulation starting from these boundary conditions is used for initialization”. As we know that CH4 background uncertainty (bias) will be carried on to calculated CH4 enhancement, which is directly related to posterior CH4 emissions, what the bias of CH4 background in this study and potential uncertainty in deriving CH4 emissions?
Line 110, “It attributes 94% of Permian emissions to oil and gas activity, and we assume the same fraction for our posterior emission estimates.”, Can you clarify what is the potential uncertainty of using the constant fraction of 94% to oil and gas activity in inversion results? As I know most inversion studies have the ability to constrain posterior emissions from different categories.
Section 2.2, As displayed in the reference list, there are some other inversion studies in the same domain, (i.e. Zhang et al., 2020), it's better to illustrate what the main improvement of your study when compared with these previous studies, because it's very hard for audiences to remember and distinguish the method difference between all related studies.
Line 125, “mitigate boundary-condition errors (Shen et al., 2021; Varon et al., 2022).”, have you assessed the improvement of CH4 background with observations?
Line 136-137, “The error covariance matrices and are assumed diagonal with uniform error standard deviations of 50% and 15 ppb, respectively”, As I understand, the 50% uncertainty for prior inventory may represent regional averages not all grid cells in study domain, which can be much larger than 50%, the same as 15 ppb for observations and GEOS-Chem simulations, so whether the inversion results are sensitive to different values of 50% and 15 ppb, if you assign a slightly larger extent (i.e. 80%, and 20 ppb), how much will the results change?
In figure 3, the weekly emission changes can vary by 100%, indicating the potential bias or uncertainty of CH4 emissions at weekly scale can be much larger than 50%.
Line 190, It seems the use of proportion 94% will largely influence your results of CH4 emissions from oil and gas. I am curious why the inversion model cannot constrain CH4 emissions from each category?
Data displayed in Figure 5 for model simulated xCH4 and observation.
Overall, why tower based CH4 concentrations seem higher than simulations with both prior and posterior emissions (scatter plot is below 1:1 line)? Whether it indicates the posterior CH4 emissions are still underestimated? How about plotting time series of concentration?
Whether one of the reasons for the large difference between model simulation and tower observations is the vertical gradient in the lowest GEOS-Chem model? And what is the height of the lowest model level? or aggregation error for spatial resolution between the point scale and regional scale(25km)?
Line 256-257, “The mean satellite-inferred emission (0.72 Tg is 20% lower than the mean tower (0.88 Tg and Scientific Aviation (0.89 Tg estimates during the period of overlap”, From the above concentration comparisons, I also agree that the satellite-inferred emissions are obviously underestimated.
Line 398-410, It's better to display the comparisons between atmospheric inversions and multiple linear regression with figures instead of only using tables and numbers.
Line 438, “assuming 80% methane content for Permian”, whether the assumption of using 80% is reasonable, and what the general extent of this value in the study region?
Technical comments
Line 398, “summarizes the models and results” I just guess the first author forget to delete this sentence of comment from other co-authors.
- AC1: 'Comment on acp-2022-749: Responses to reviews', Daniel Varon, 23 Mar 2023