Articles | Volume 23, issue 13
Research article
 | Highlight paper
11 Jul 2023
Research article | Highlight paper |  | 11 Jul 2023

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, Zhen Qu, Zachary R. Barkley, Natasha L. Miles, Scott J. Richardson, Kenneth J. Davis, Sudhanshu Pandey, Xiao Lu, Alba Lorente, Tobias Borsdorff, Joannes D. Maasakkers, and Ilse Aben

We quantify weekly methane emissions at 0.25× 0.3125 (≈25× 25 km2) resolution from the Permian Basin, the largest oil production basin in the US, by inverse analysis of satellite observations from the TROPOspheric Monitoring Instrument (TROPOMI) from May 2018 to October 2020. The mean oil and gas emission from the region (± standard deviation of weekly estimates) was 3.7 ± 0.9 Tg a−1, higher than previous TROPOMI inversion estimates that may have used biased prior emissions or background assumptions. We find strong week-to-week variability in emissions superimposed on longer-term trends, and these are consistent with independent inferences of temporal emission variability from tower, aircraft, and multispectral satellite data. New well development and natural gas spot price were significant drivers of variability in emissions over our study period but the concurrent 50 % increase in oil and gas production was not. The methane intensity (methane emitted per unit of methane gas produced) averaged 4.6 % ± 1.3 % and steadily decreased from 5 %–6 % in 2018 to 3 %–4 % in 2020. While the decreasing trend suggests improvement in operator practices during the study period, methane emissions from the Permian Basin remained high, with methane intensity an order of magnitude above the industry target of <0.2 %. Our success in using TROPOMI satellite observations for weekly estimates of emissions from a major oil production basin shows promise for application to near-real-time monitoring in support of climate change mitigation efforts.

1 Introduction

Sharp reductions in atmospheric methane emissions are needed to slow near-term climate change (IPCC, 2022). Emissions from the oil and gas industry have been identified as a high priority because implementing new practices can mitigate a significant quantity of emissions at no net cost (IEA, 2022). But the magnitude of oil and gas methane emissions is heavily disputed. National inventories use bottom-up methods to estimate emissions from infrastructure data, limited device-level measurements, and engineering models. This is difficult because emissions from oil and gas production typically involve thousands of individual point sources, some of which may be individually large (Irakulis-Loixate et al., 2021) and others individually small but accumulating to large regional totals (Omara et al., 2022). Complicating matters is the strong and poorly understood temporal variability of emissions, both in terms of the intermittency of individual point sources (Cusworth et al., 2021a) and regional fluctuations across oil and gas production basins (Lin et al., 2021; Cusworth et al., 2022). Observations of atmospheric methane have been used extensively for top-down estimates of methane emissions through inversion of atmospheric transport models relating emissions to concentrations (Houweling et al., 2017). Here we show how satellite observations can be used to quantify weekly temporal variability in oil and gas methane emissions from a major production basin (the Permian) over a  30-month period, and we show that this temporal variability can be explained by specific activity drivers.

Satellite instruments can sense atmospheric methane from backscattered sunlight in the shortwave infrared (SWIR) to monitor methane emissions (Jacob et al., 2016, 2022). The TROPOspheric Monitoring Instrument (TROPOMI) launched in 2017 on the Sentinel-5 Precursor (S5P) satellite has been transformative in enabling high-resolution mapping of methane emissions at regional scales (de Gouw et al., 2020). TROPOMI provides daily global coverage at 5.5 × 7 km2 nadir pixel resolution (Hu et al., 2018; Schneising et al., 2019; Lorente et al., 2021a). A number of studies have used TROPOMI observations for inverse analyses of emissions at the scale of individual oil and gas basins (Zhang et al., 2020; Schneising et al., 2020; Liu et al., 2021; Shen et al., 2021, 2022; McNorton et al., 2022) but have generally focused on seasonal or annual estimates. Emissions can vary on shorter timescales due to economic drivers (Lyon et al., 2021), accidents (Pandey et al., 2019; Cusworth et al., 2021b), and the lifecycle of facilities (Cardoso-Saldaña and Allen, 2020, 2021; Allen et al., 2022). This short-term variability is important to understand for interpreting field studies and projecting future trends.

The Permian Basin in Texas and New Mexico is the largest oil-producing basin in the US, accounting for more than 40 % of US oil production (FRBD, 2022). Bottom-up and top-down estimates of methane emissions from the Permian vary widely, from 0.6 Tg a−1 as reported in the gridded version of the 2012 US Environmental Protection Agency (EPA) greenhouse gas inventory (GHGI; Maasakkers et al., 2016) extrapolated to 2018 (Shen et al., 2022) to more than 3.7 Tg a−1 as inferred by Shen et al. (2022) from 2018–2020 TROPOMI satellite observations. Chen et al. (2022) reported combined emissions of 194 t h−1 (1.7 Tg a−1) from point sources surveyed aerially across the New Mexico Permian Basin from October 2018 through January 2020, corresponding to 9.4 % of the region's gross gas production, a methane intensity more than twice that inferred from basin-wide satellite observations (Zhang et al., 2020; Schneising et al., 2020; Liu et al., 2021; Shen et al., 2022). These discrepancies may be due in part to spatiotemporal emission variability as has been reported for different sub-regions of the Permian during the 2020 COVID-19 shutdowns (Lyon et al., 2021) and between individual 2019–2021 aircraft campaigns (Cusworth et al., 2022). High-frequency TROPOMI satellite observations offer a means to probe this variability at the basin scale.

Here we conduct an inverse analysis of TROPOMI satellite observations to characterize the magnitude and spatiotemporal variability of methane emissions from the Permian Basin on a weekly basis, starting from the best available bottom-up prior estimates of emissions and using Bayesian inverse modeling to obtain optimized posterior estimates assimilating the information from TROPOMI. We use a suboptimal Kalman filter (Todling and Cohn, 1994) to quantify weekly basin-wide emissions at 0.25× 0.3125 (≈25× 25 km2) resolution for 127 weeks from 1 May 2018 to 5 October 2020. The combination of fine spatiotemporal resolution and long-term basin-wide coverage allows us to directly compare our results with previous space-based and aerial emission estimates, which cover different time periods and sub-regions within the basin. To evaluate the Kalman filter's ability to resolve weekly variability, we also compare our results with independent tower and aircraft measurements, inverse analyses of those measurements, and Sentinel-2 satellite detections of a strong methane point source. Finally, we interpret our weekly estimates by comparing them with oil and gas activity data and economic variables to gain insight into the factors controlling methane emissions from the Permian Basin, as a step towards using satellite observations to guide improvements in bottom-up inventories and identify opportunities to decrease emissions.

2 Materials and methods

2.1 TROPOMI observations, GEOS-Chem forward model, and prior emission estimates

TROPOMI was launched into sun-synchronous orbit on the S5P satellite in October 2017. It operates in a nadir-viewing push-broom configuration with a 2600 km wide swath and  13:30 local overpass time (Veefkind et al., 2012). The TROPOMI operational data record begins in May 2018 and provides daily global observations of dry-air column-average methane mixing ratios (XCH4) at 5.5 × 7 km2 nadir pixel resolution (7 × 7 km2 before August 2019). Here we use 2.5 years of version 02.02.00 observations from the scientific data product (Lorente et al., 2021a), from 1 May 2018 to 5 October 2020, and select only high-quality retrievals with quality assurance values greater than 0.5 (Hasekamp et al., 2022). On average (5th, 95th percentile) we obtain 19 346 (5674, 39 502) observations per week over our full inversion domain (25–38 N, 96–110 W), including 3062 (414, 6707) per week within the Permian itself (Fig. 1).

We use the nested version of the GEOS-Chem chemical transport model version 12.7.1 as forward model for the inversion. The simulation is driven by NASA GEOS Fast Processing (GEOS-FP) meteorological fields at 0.25× 0.3125 resolution over 47 vertical layers from the surface to the mesopause. We use dynamic 3 h boundary conditions from a global 4× 5 simulation corrected with spatially and temporally smoothed TROPOMI data as described by Shen et al. (2021). A 1-month spin-up simulation starting from these boundary conditions is used for initialization. TROPOMI observation of the GEOS-Chem atmosphere is simulated by applying the TROPOMI instrument averaging kernels and prior methane vertical profiles for individual retrievals to the GEOS-Chem vertical profiles of methane mixing ratios, as described by Varon et al. (2022).

Figure 1TROPOMI methane data, prior and posterior emission estimates, and state vector emission elements for our Permian inversion. The thick black contour defines the region of interest for the inversion as the geological extent of the Permian Basin. (a) TROPOMI dry-air column-average methane mixing ratios (XCH4) for a typical week (1–8 May 2018), with 23 059 observations plotted on a 0.1× 0.1 grid, including 4937 within the Permian. (b) Prior emission estimates on the 0.25× 0.3125 GEOS-Chem forward model grid, including the Environmental Defense Fund (EDF) inventory (Zhang et al., 2020) for the Permian. (c) State vector emission elements to be optimized in the inversion, including 235 elements at 0.25× 0.3125 resolution within the Permian and eight coarse buffer elements to correct boundary conditions. (d) Mean posterior emissions from the weekly inversions for the May 2018 to October 2020 study period, on the same grid as (b).

Our prior emission estimates for the Permian are from the 2018 Environmental Defense Fund (EDF) bottom-up inventory of Zhang et al. (2020) (2.7 Tg a−1; also referred to as the EIME inventory), which provides a much better representation of atmospheric methane observations than the gridded version of the US EPA GHGI (Maasakkers et al., 2016). The EDF inventory features two emission maxima in the Delaware and Midland sub-basins also seen in the TROPOMI observations (Fig. 1). It attributes 94 % of Permian emissions to oil and gas activity. Prior anthropogenic emissions for the rest of the inversion domain in the US (state vector buffer elements outside the Permian) are from the gridded EPA GHGI. Anthropogenic emissions for northern Mexico are from the Global Fuel Exploitation Inventory (GFEI; Scarpelli et al., 2020) for oil, gas, and coal emissions and from the EDGAR v4.3.2 inventory (Janssens-Maenhout et al., 2019) for non-fossil sources. Monthly wetland emissions are from the WetCHARTS v1.2.1 inventory ensemble mean (Bloom et al., 2017), and other natural emissions and sinks are as described by Lu et al. (2022).

2.2 Kalman filter inversion methodology

We quantify weekly average methane emissions from the Permian Basin by Kalman filter inversion of TROPOMI XCH4 observations for 127 weeks from 1 May 2018 to 5 October 2020. In this framework, the posterior emission estimates for each week determine the prior estimates for the next week, enabling the inversion to quantify temporally variable emissions by providing regular updates as new observations are introduced.

For a given week, the GEOS-Chem forward model F relates the weekly average methane emissions (state vector x) to TROPOMI XCH4 observations (observation vector y) as y=Fx+ϵO, with observational error ϵO stemming from uncertainties in both the forward model and observations. Here x (length n=243) contains spatially gridded emission elements for the Permian and its surroundings, including 235 elements at 0.25× 0.3125 resolution within the basin and eight additional coarse-resolution buffer elements adjacent to the basin (Fig. 1c) to mitigate boundary-condition errors (Shen et al., 2021; Varon et al., 2022). The observation vector y (length m) assembles the TROPOMI observations for the week, to be ingested sequentially in the inversion. The inversion optimizes x to match the observations y subject to constraints from prior information, using optimal estimation by Bayesian synthesis with Gaussian error statistics to retrieve posterior emission estimates x+ of the form (Brasseur and Jacob, 2017)

(1) x + = argmin x x - x - T S - 1 x - x - + γ y - K x T R - 1 y - K x .

Here x contains the prior emission estimates, S (n×n) is the prior error covariance matrix describing uncertainty in the prior, R (m×m) is the observational error covariance matrix describing combined uncertainty in the forward model and observations, γ is a regularization parameter to prevent overfitting, and the Jacobian matrix K=F(x)/x defines the sensitivities of the observations to the emissions as described by the forward model y=F(x). F(x) is linear so that K fully defines the forward model for the purpose of the inversion. We compute K column by column via n forward model simulations, perturbing one state vector element at a time. The error covariance matrices S and R are assumed diagonal with uniform error standard deviations of 50 % and 15 ppb, respectively. An observational error standard deviation of 15 ppb is typical for TROPOMI as determined by the residual error method (Heald et al., 2004; Qu et al., 2021). Ignoring error correlations, which are difficult to quantify, exaggerates the information content of the observations and underestimates the prior error on emissions when aggregated on the regional scale (Worden et al., 2022; Chen et al., 2023). We effectively correct for this through the regularization factor γ (Lu et al., 2021). The standard Kalman filter would allow S to vary from week to week with x; we discuss the rationale for holding it fixed below. We set γ=0.25 following Shen et al. (2021) and investigate the sensitivity to this choice through additional inversions with γ=0.10 and γ=0.30.

Our Kalman filter to quantify weekly methane emissions involves a prediction step and an update step for each week. In the prediction step, prior estimates xk- of the emissions during week k are made from the most recent posterior estimates xk-1+. We assume persistence of emissions to advance the Kalman filter from week to week. Normally this would be done via direct substitution (xk-=xk-1+), but we find that this can cause some emission elements to become locked at very low values because the prior error (set as a relative 50 %) then becomes very low as well. To avoid this problem, in the prediction step we retain information from the original prior by computing xk- as

(2) x k - = λ ( α x k - 1 + + 1 - α x 1 - ) .

Here x1- contains the prior estimates for the first week (Fig. 1b), α=0.9 nudges towards the original prior estimates x1- with 90 % weight to the latest posterior estimate xk-1+, λ is a scale factor to ensure that xk- has the same mean as xk-1+, and in the first week (k=1) we have x0+=x1-. We only use Eq. (2) to compute prior estimates for emissions within the basin; emissions in the external buffer elements (Fig. 1) are assigned the same prior estimates (from x1-) each week. We examine the sensitivity of our results to the choice of α through additional inversions with α=0.8 and α=0.95.

The standard Kalman filter prediction step also estimates the prior error covariance matrix Sk- for the state at time k, by applying a state evolution model to the most recent posterior Sk-1+. We instead adopt a suboptimal Kalman filter with fixed (diagonal) error covariance matrix Sk-=S for each week, effectively adding error to the state evolution model described by Eq. (2) such that the prior error variance does not decrease with time. This prevents the prior error estimates from gradually becoming smaller and weakening further corrections to emissions.

In the update step, the Kalman filter evaluates Eq. (1) for week k to generate posterior estimates of the emissions (Rodgers, 2000)

(3) x k + = x k - + ( γ K k T R k - 1 K k + S - 1 ) - 1 γ K k T R k - 1 ( y k - K k x k - )

and associated error covariance matrix

(4) S k + = ( γ K k T R k - 1 K k + S - 1 ) - 1 .

The information content of the week's inversion is then computed from the averaging kernel matrix

(5) A k = I n - S k + S - 1 ,

where In is the n-dimensional identity matrix. The diagonal elements of Ak are called averaging kernel sensitivities and take values between 0 and 1. They describe the extent to which each emission element in the state vector is constrained by the observations of week k as opposed to the prior estimate. An emission element with unit sensitivity is independent of the prior and entirely constrained by the observations; one with zero sensitivity is not at all constrained by the observations and assumes the prior value. The sum of the averaging kernel sensitivities (trace of Ak) defines the degrees of freedom for signal (DOFS). The DOFS value measures the information content of the inversion as the number of independent pieces of information on emissions that can be constrained by the observations.

Figure 2Flowchart of our Kalman filter inversion methodology for estimating weekly emissions. Starting from a set of initial conditions and prior emission estimates for the inversion domain (Fig. 1), a 1-week inversion produces posterior emission estimates for week k. The posterior emission estimates are then nudged towards the original EDF inventory emissions before substitution to serve as prior estimates for week k+1 and are used in a 1-week simulation for week k to generate initial conditions for week k+1.


Figure 2 illustrates the procedure for advancing the Kalman filter in time. Starting from a set of initial conditions and prior emission estimates xk- for week k, a 1-week inversion of TROPOMI observations generates posterior emission estimates xk+ (Eq. 3). Initial conditions for the first week (1–8 May 2018) are obtained from the 1-month spin-up simulation. The posterior estimates are then nudged towards the EDF emission inventory (Fig. 1b) following Eq. (2), and the nudged emissions are used as prior estimates for week k+1. Finally, a 1-week GEOS-Chem simulation is performed for week k, driven by the posterior emission estimates xk+, to obtain new initial conditions for the following week. This procedure is repeated for the 127 weeks of the May 2018 to October 2020 study period.

3 Results and discussion

3.1 Permian methane emissions

Figure 3 shows our weekly emission estimates for the Permian region of interest (Fig. 1) from 1 May 2018 to 5 October 2020 (127 weeks). Here and elsewhere, we assume that 94 % of total emissions in the Permian are from the oil and gas sector, as given by the prior estimate, and this is not a significant source of error because emissions are mainly from the oil and gas sector in any case. Mean weekly emissions for the period expressed as Tg a−1 are (mean ± standard deviation) 3.9 ± 1.0 Tg a−1, including 3.7 ± 0.9 Tg a−1 from oil and gas. The long-term trend over the 127-week record shows a weak but statistically significant decrease (p<0.01) of −3.5 Gg a−1 per week, with a gradual increase of about 20 Gg a−1 per week from May 2018 to mid-June 2019 followed by a gradual decrease of about 14 Gg a−1 per week through September 2020. Superimposed on these slow trends are strong fluctuations at multi-week scales, often by a factor of 1.5–2 or more, which we examine further in Sect. 3.6. This includes a rapid reduction in emissions in the months following the March 2020 COVID-19 shutdown and several earlier periods of strong variability. Running GEOS-Chem with our weekly posterior emission estimates improves the mean model–TROPOMI bias over the Permian from −9.0 ppb to −2.0 ppb.

Figure 3Weekly methane emissions from the Permian Basin from 1 May 2018 to 5 October 2020. The Permian is defined by the thick black contour in Fig. 1. Emissions are inferred by inversion of weekly TROPOMI satellite observations using a suboptimal Kalman filter. The mean emission (± standard deviation) is 3.9 ± 1.0, with 3.7 ± 0.9 Tg a−1 from oil and gas. The dashed lines are 4-week moving averages of the weekly emission estimates from four sensitivity inversions varying inversion parameters γ and α and a fifth one (with low initial values) decreasing the prior emission inventory.


The dashed lines in Fig. 3 are from four sensitivity inversions investigating the influence of the parameters γ (0.10–0.30) and α (0.80–0.95) and a fifth that scales down the EDF inventory to use a prior estimate for regional total emissions of 1.2 Tg a−1 as done by Zhang et al. (2020) instead of 2.7 Tg a−1. The posterior mean (3.9–4.1 Tg a−1) and standard deviation (0.8–1.0 Tg a−1) vary little across these sensitivity inversions. The error standard deviation of the weekly emission estimates, expressed as the standard deviation of residuals between the base inversion and the five sensitivity inversions, is 0.45 Tg a−1 ( 10 % of the mean emission). Increasing γ tends to increase the range of weekly emissions by assigning more weight to the observations and thus allowing the optimization more freedom to depart from the prior estimate. Increasing α produces a similar effect by assigning less weight to the EDF bottom-up inventory during the Kalman filter nudging operation (Eq. 2). The inversion with lower prior inventory converges to the other inversions after 7–8 weeks of burn-in, demonstrating the ability of the observations to quantify mean emissions independently of the prior estimate. We discuss the ability to quantify weekly variability in what follows.

3.2 Observation density and information content

Our ability to quantify the weekly variability of basin-wide emissions depends on the information content of the Kalman filter as determined by the number of TROPOMI observations and the resulting DOFS. Fig. 4 shows how the information content of our weekly inversions varies with the number of TROPOMI observations. Weekly observations for the Permian region of interest (Fig. 1) range from 5 to 13 288, with a mean of 3062, for our 127 weeks. The resulting weekly DOFS within the Permian range from 0.02 to 14.09, with a mean of 4.67. The number of DOFS increases with observation count by roughly 1.4 DOFS per 1000 TROPOMI observations. Shen et al. (2022) identified DOFS >0.5 as a practical minimum to estimate total basin methane emissions with 2σ error ≤30 % from inverse analysis of TROPOMI observations. Inversions with low DOFS are mainly constrained by the prior emission estimate and may introduce smoothing error in the inference of weekly temporal trends. Here 124 of our 127 weekly inversions meet the DOFS >0.5 criterion. The three remaining inversions are 45 and 7 weeks apart so that we do not observe low DOFS in consecutive weeks. TROPOMI thus provides sufficient observational information to support weekly estimation of total Permian emissions for the vast majority of our study period.

Figure 4Weekly information content of TROPOMI observations over the Permian Basin. The figure shows the dependence of the weekly degrees of freedom for signal (DOFS) for the inversion on the weekly number of TROPOMI observations in the Permian region of interest (Fig. 1), for 127 weeks of inversions. The DOFS represent the number of independent pieces of information on the emission distribution that can be quantified by the inversion independently from the prior estimate. The top histogram shows the distribution of weekly observation counts in the Permian. The right histogram shows the distribution of weekly DOFS in the Permian. Also shown is the reduced-major-axis regression line, indicating 1.4 DOFS per 1000 TROPOMI observations.


3.3 Verification with independent tower and aircraft observations

In situ methane measurements from the Permian Methane Analysis Project (PermianMAP; Lyon et al., 2021) provide an opportunity for independent evaluation of our results. PermianMAP is an initiative started in November 2019 that gathers methane information from ground-based, airborne, and satellite measurement platforms to monitor emissions from the Permian, with particular attention to a highly active 100 km × 100 km region of the Delaware sub-basin (31.4–32.4 N, 103.2–104.2 W). We begin here by comparing our weekly inversions with PermianMAP measurements from the Permian Basin tall tower network (Monteiro et al., 2022). The network comprises five towers, and the three with the most continuous hourly measurements of methane mole fraction since March 2020 are located near Fort Stockton (30.867 N, 102.815 W), Carlsbad (32.178 N, 104.441 W), and Notrees (31.966 N, 102.770 W), Texas. We compare measurements from these three towers with GEOS-Chem methane simulations driven by (1) constant prior emissions from the EDF bottom-up inventory of Zhang et al. (2020) and (2) our weekly posterior emission estimates.

Figure 5Independent evaluation of our weekly methane emission estimates with tall tower measurements of methane concentrations in the Permian. The figure compares weekly mean afternoon (20:00–24:00 UTC) methane concentrations simulated by GEOS-Chem using prior (constant) and posterior (weekly updated) emissions to tall tower measurements by Monteiro et al. (2022) from March through September 2020 at Fort Stockton (30.867 N, 102.815 W), Carlsbad (32.178 N, 104.441 W), and Notrees (31.966 N, 102.770 W), Texas. The insets show the locations of the three towers (stars) within the Permian (thick black contour). Also shown are the mean bias (MB) and root-mean-square error (RMSE). The 1:1 line is shown in black.

Figure 5 shows the locations of the three towers in the basin and compares their weekly mean methane concentrations with the GEOS-Chem simulations in the nearest model grid cell for afternoon hours (20:00–24:00 UTC) when the mixed layer is fully developed (Barkley et al., 2023). GEOS-Chem underestimates methane concentrations at the three towers, but in all cases our weekly posterior emission estimates improve both the mean bias and RMSE. The improvement in the RMSE is limited by GEOS-Chem representation and transport error, which Lu et al. (2021) previously estimated at 36–68 ppb for surface and tower data.

We further compare in Fig. 6 our posterior weekly emissions for the PermianMAP region in January–October 2020 with independent time-resolved estimates from Scientific Aviation aircraft flights on different days and an inverse analysis of measurements from all five towers in the Permian Basin tower network, both reported by Barkley et al. (2023). The satellite and tower inversions generally agree within errors and show a similar decreasing emission trend over the March–October 2020 period. The mean satellite-inferred emission (0.72 Tg a−1) is 20 % lower than the mean tower (0.88 Tg a−1) and Scientific Aviation (0.89 Tg a−1) estimates during the period of overlap (March through September 2020). However, linear interpolation of our weekly inversions yields very close agreement with the individual Scientific Aviation estimates (R=0.6, RMSE = 0.18 Tg a−1, mean deviation = 0.004 Tg a−1).

Figure 6Weekly methane emissions from the Lyon et al. (2021) PermianMAP study region (red box in inset: 31.4–32.4 N, 103.2–104.2 W). Results from our weekly TROPOMI inversions are compared to monthly tower inversions and aircraft model-based scaling-factor estimates reported by Barkley et al. (2023). The grey shaded area and orange bars represent sensitivity ranges.

3.4 Comparison with Sentinel-2 point source observations

Satellite observations of methane point sources in the Permian provide an additional opportunity to evaluate our inversion results. Single detections of Permian Basin point sources have been reported by Irakulis-Loitxate et al. (2021) and Sánchez-García et al. (2022) using the PRISMA, GF-5, ZY1, and WV-3 satellite instruments. More recently, Ehret et al. (2022) tracked intermittent emissions exceeding  5 t h−1 from a super-emitting compressor station in the Midland sub-basin (31.731 N, 102.042 W) over 3 months (July–September 2020) using the Sentinel-2 and Landsat 8 satellites. Here we reconstruct the emissions from that specific point source over the full TROPOMI record by applying the multi-band–multi-pass retrieval method (Varon et al., 2021) to Sentinel-2 observations made over that location every 2–3 d and quantifying source rates using the integrated mass enhancement (IME) method (Varon et al., 2018, 2021) with local GEOS-FP wind speed data. We obtain a detection threshold of about 2 t h−1 as our lowest detected emission.

Figure 7 compares the resulting time series of point source emission rates from the compressor station to the results of our TROPOMI inversion for the corresponding 0.25× 0.3125 grid cell. The point source was detected on 24 out of 269 non-cloudy passes (9 % persistence), with four periods of intermittent activity: May–August 2018, March–May 2019, August–September 2019, and July–September 2020. Our posterior estimates also show emission spikes during these periods, although we do not capture the high values in July–August 2020. Direct comparison between the two time series is challenging because of differences in spatiotemporal scale; there are hundreds of facilities within the 0.25× 0.3125 grid cell, and the Sentinel-2 observations are instantaneous snapshots, whereas our TROPOMI analysis averages emissions over the week. Nevertheless, the mean TROPOMI emissions for that grid cell (9100 kg h-1=3.1× 10−9 kg m−2 s−1) are among the highest in the Permian (Fig. 1d).

Figure 7Temporal variability of methane emissions from the Permian (Midland) Parks compressor station (31.731 N, 102.042 W; red star in inset map) over the May 2018 to October 2020 period. Point source detections from the Sentinel-2 satellite instrument are compared to weekly TROPOMI estimates for the corresponding 0.25× 0.3125 GEOS-Chem grid cell. The red and grey vertical markers indicate detections and non-detections on non-cloudy passes, respectively.

3.5 Comparison with previous TROPOMI estimates of Permian emissions

Table 1 compares our inversion results to previous basin-scale top-down estimates of Permian emissions using TROPOMI observations. The definition of the Permian varies between studies, which may introduce some discrepancies. The studies are also for different periods, but we account for that with our resolved weekly emissions. Zhang et al. (2020) were the first to use TROPOMI to estimate total methane emissions from the basin, and they reported 2.9 ± 0.5 Tg a−1 for an 11-month inversion of 2018–2019 data. Schneising et al. (2020) applied a Gaussian integral technique to TROPOMI methane data for 2018–2019 (all 24 months) and estimated total basin emissions of 3.2 ± 1.1 Tg a−1. Liu et al. (2021) developed a divergence method for inferring methane emissions from TROPOMI observations and reported 3.1 (2.8, 3.8) Tg a−1 for the Permian in 2019. McNorton et al. (2022) obtained total emissions of 2.3 ± 0.5 Tg a−1 from a 4D-Variational (4D-Var) inversion of TROPOMI data for 15 months during 2019–2020. Shen et al. (2022) reported 3.7 ± 0.5 Tg a−1 from oil and gas for a 22-month analytical inversion of TROPOMI observations from May 2018 to February 2020.

Table 1Comparison of our work with previous TROPOMI estimates of Permian methane emissions.a

a All authors reported combined emissions for all sectors in the basin, except Shen et al. (2022), who reported oil and gas emissions only. Oil/gas accounts for 94 % of Permian emissions in our work. b Our weekly emission estimates are matched to these different periods, except for Schneising et al. (2020; see footnote). c Values represent mean ± standard deviation of weekly estimates for the period of overlap with the previous study. d The definition of the reported uncertainty varies between studies. e Comments refer to the previous studies. f This includes all 24 months in 2018–2019; we show our period-average result from May 2018 through 2019 for comparison. g This is split across two discrete periods. h Shen et al. (2022) performed a 22-month analytical inversion with 16 state vector buffer elements.

Download Print Version | Download XLSX

Our mean emission for the Permian (3.9 Tg a−1; 3.7 Tg a−1 for oil and gas only) is 20 %–35 % higher than these previous TROPOMI-based estimates, except for the 3.7 Tg a−1 oil and gas emission reported by Shen et al. (2022). In contrast, our emission for the PermianMAP 100 km × 100 km sub-region in 2020 is on average 20 % lower than estimated by inverse analysis of tower and aircraft observations (Barkley et al., 2023; Sect. 3.3). The TROPOMI estimates of Table 1 are for the most part consistent within errors, but the discrepancy between our results and previous (lower) ones appears to reflect biases in background specification and the spatial distribution of prior emission estimates. Zhang et al. (2020) used a prior estimate of 1.2 Tg a−1 with spatial distribution from the 2012 GHGI in their base inversion and included a spatially uniform model bias term in their state vector to remove the effects of boundary-condition (background) errors on the regional inversion. Permian emissions in the 2012 GHGI are known to be too low and spatially biased (Zhang et al., 2020), and this would impede the inversion's ability to obtain the correct solution (Yu et al., 2021). In an alternative inversion with the higher and more spatially accurate EDF inventory (2.7 Tg a−1) as prior estimate, Zhang et al. (2020) reported 3.2 Tg a−1. Here we use the EDF prior but with TROPOMI-derived boundary conditions and eight state vector buffer elements to guard against boundary-condition error. Our inversion is largely insensitive to downscaling of the EDF inventory to 1.2 Tg a−1 (Fig. 3) thanks to preservation of the inventory's accurate spatial distribution through the Kalman filter nudging operation (Fig. 2). Shen et al. (2022), with whom our results agree closely, used a similar prior estimate (2.2 Tg a−1 with accurate spatial distribution) and the same buffer approach but with 16 elements, and when they instead used a 0.6 Tg a−1 prior estimate, they reported a moderately lower emission of 2.9 Tg a−1, on par with the Zhang et al. (2020) result despite the much lower prior estimate. Thus, the range of estimates from these three inversion studies can be explained by differences in prior emission estimates (mainly spatial distribution) and background specification. McNorton et al. (2022) used a prior estimate of 2.0 Tg a−1 but obtained only 2.3 Tg a−1 from a 4D-Var inversion optimizing both emissions and model transport, which could reflect smoothing error if the observations have insufficient information to constrain both.

The Gaussian integral method of Schneising et al. (2020) and divergence method of Liu et al. (2021) do not rely on prior estimates but are subject to large errors from uncertainty in the local background subtraction. Schneising et al. (2020) defined the background from a 2× 4 upwind box and identified numerous days with negative or very low (<1 Tg a−1) emissions despite stringent filtering criteria for the background observations. Liu et al. (2021) defined the methane background from the local neighborhood of individual model grid cells, which aliased Permian enhancements into the background field.

3.6 Attribution of methane emission trends

A variety of factors can influence methane emissions from oil and gas activity on weekly to monthly scales, including economic drivers (Lyon et al., 2021), regulatory changes (Cardoso-Saldaña and Allen, 2021), accidents (Pandey et al., 2019; Cusworth et al., 2021b), infrastructure development (Lyon et al., 2021), and the lifecycle of facilities (Cardoso-Saldaña and Allen, 2020, 2021; Allen et al., 2022). Figure 8 compares our weekly TROPOMI-derived emission estimates (Fig. 8a) with a range of economic and activity variables (Fig. 8b–f) to investigate possible causes of Permian emission variability.

Figure 8b shows monthly data for Permian oil and gas production from Enverus Drillinginfo (Enverus, 2021). Rising production could lead to higher emissions from increasing activity but also to lower emissions if more by-produced gas is brought to market rather than being vented or inefficiently flared. Permian oil and gas production increased steadily by about 50 % from May 2018 through March 2020, but this was not associated with an increase in emissions, which may reflect increased gas marketing. Production declined sharply in April–May 2020 following the COVID-19 shutdown and gradually recovered from June through September 2020. Lyon et al. (2021) reported concurrent declines in methane emissions and production levels for the PermianMAP 100 km × 100 km study region during the COVID-19 shutdown period, but like Barkley et al. (2023) we find that basin-wide emissions continued to decline after production recovered.

Figure 8c describes monthly trends in new wells entering first production (Enverus, 2021). New well development can drive variability in emissions because new wells tend to produce high methane emissions in the first few months from well completion and tank flashing that then decline rapidly over time (Cardoso-Saldaña and Allen, 2020; Allen et al., 2022). The number of new wells entering production in the Permian was relatively steady at  600 new wells per month leading up to fall 2019 but spiked to  1200 wells per month in October 2019 and then again in March 2020 before plunging to  300 wells per month through September 2020. Both spikes precede a local maximum in emissions by  1 month, and the sharp reduction in well development in April 2020 matches the drop in emissions beginning  1 month later. This confirms the importance of new wells for driving emissions as was reported by Lyon et al. (2021) for the PermianMAP study region during the 2020 COVID-19 shutdown. If the increases in monthly mean emissions during October 2019 and March 2020 are attributed entirely to the spikes in new wells, then we estimate 5500–6300 kg d−1 per new well during the first month of production. This is 1–2 orders of magnitude higher than bottom-up estimates for new wells in the Barnett Shale and Eagle Ford basins using the Methane Emission Estimation Tool (MEET; Allen et al., 2022) and may reflect routine flaring and venting of by-produced gas at new well sites not yet connected to the Permian gas pipeline network (Rystad Energy, 2021). Point sources >1000 kg h−1 have frequently been detected in previous airborne and satellite surveys of the Permian (Cusworth et al., 2021a, 2022; Irakulis-Loitxate et al., 2021; Chen et al., 2022).

Figure 8d shows monthly data for Permian natural gas pipeline takeaway capacity and outflows to external markets (S&P Global Platts, 2022). Outflows reflect natural gas production volumes minus gas consumed within the basin. Increasing pipeline capacity and outflows could reduce emissions by increasing the marketed proportion of by-produced gas. Pipeline capacity and outflows increased over the study period by 34 % and 57 %, respectively, in pace with production, and a step change in capacity occurred in September 2019 when the 2 Bcf d−1 Gulf Coast Express (GCX) pipeline entered commercial service (S&P Global Platts, 2019). Paradoxically, in the 4 months preceding this step change, total Permian outflows appear to have exceeded total capacity. This is likely the result of preparatory (commissioning and line packing) flows in the months leading up to the commercial deployment of the GCX pipeline. Partial capacity up to 2 Bcf d−1 of GCX capacity may have become practically available as early as May 2019 (dashed line in Fig. 8d), in which case the steep emission reduction we observe in June 2019, from  5.5 Tg a−1 to  3 Tg a−1, may reflect increased gas marketing facilitated by the new pipeline.

Figure 8e shows the number of monthly flaring targets across the Permian detected by the Visible Infrared Imager Radiometer Suite (VIIRS) satellite instrument using the Elvidge et al. (2016) flare-monitoring methodology (Lyon et al., 2021). Flaring activity produces methane emissions due to incomplete combustion of flared gas (Plant et al., 2022), and unlit flares vent methane directly to the atmosphere. Lyon et al. (2021) reported combustion issues and unlit flares at, respectively, 6 % and 5 % of flare sites observed during three Permian surveys performed in 2020. The number of VIIRS-detected flaring targets in the Permian followed an upward trend from May 2018 through August 2019 and then a downward trend through the end of our study period. This is consistent with the gradual increase in emissions seen from May 2018 to mid-June 2019 and subsequent decrease through September 2020 and may reflect increased gas marketing post-GCX. The  25 % spike in emissions observed in November 2018 coincides with a short-lived  40 % increase in flaring targets. The sharply higher flaring activity during this period may suggest increased methane emissions from flaring.

Figure 8f shows daily spot prices for the West Texas Intermediate Cushing (WTI-Cushing) oil benchmark (EIA, 2022a), the Henry Hub natural gas benchmark (EIA, 2022b), and the local Permian Waha Hub natural gas benchmark (S&P Global Platts, 2022). Trends in oil price could influence methane emissions by motivating more or less activity (e.g., drilling), and trends in natural gas price could change the incentive to market by-produced gas. In contrast with Lyon et al. (2021), we find that trends in WTI oil price do not track basin-wide emission trends, even during the early 2020 COVID-19 shutdown period. Henry Hub natural gas spot prices correlate positively with increased methane emissions in November 2018, but there is no obvious explanatory mechanism. The steep decline of the Waha Hub natural gas benchmark into negative price territory in March 2019 could explain the concurrent increase in emissions if operators were more motivated to vent or flare by-produced gas, and Fig. 8e shows a moderate uptick in flaring activity that month. This dynamic (low Waha Hub price, high flaring, high emissions) is also present in November 2018, and the sharp May 2019 decline in Waha Hub price matches an increase in emissions as well.

We apply linear regression analysis to the variables in Fig. 8 to quantify their influence on emissions. Regressions targeting our individual weekly estimates (black points in Fig. 8a) yield relatively low R2 scores <0.3. Targeting the 4-week moving-average emissions (red line in Fig. 8a) produces better predictive power, in part because many of the predictors are monthly quantities but also due to smoothing of inversion noise and high-frequency variability in emissions from large intermittent point sources such as the Parks compressor station of Fig. 7. We present three multiple linear regressions: one using all the variables as predictors to estimate the moving-average emissions and two using only the statistically significant predictors (p<0.05) from the previous regressions. In all cases, the predictors are first resampled to weekly values and standardized to equal variance (mean = 0, standard deviation = 1). The regressions use Tikhonov (L2) regularization with leave-one-out cross-validation of the regularization parameter to mitigate collinearity between predictors (Hastie et al., 2009).

Table 2 summarizes the models and results. The first regression (full model) uses one predictor for combined oil and gas production (barrels of oil equivalent), two predictors for new well development (number of new wells for the target week and number with a 4-week lag), and pipeline takeaway capacity assuming that the GCX pipeline began operating in May 2019 (dashed line in Fig. 8) for a total of nine predictors. Only four predictors are statistically significant: new wells (current and lagged) and the Henry Hub and Waha Hub natural gas prices. New wells and the Henry Hub price correlate positively with emissions, and the Waha Hub price correlates negatively. A simplified model uses these four predictors to estimate weekly emissions and shows similar predictive power to the full model (R2=0.43, from 0.46) but with a much lower Akaike information criterion (AIC = 2.78, from 12.87), indicating a better model. Lagged new wells remain statistically significant, confirming their importance as drivers of emissions, but p increases to 0.065 for the current new wells, suggesting weaker predictive power. A second simplified model, with R2=0.33 and AIC =−1.55, drops this predictor along with the Henry Hub gas price, since we lack a mechanism to explain its positive correlation with emissions during the study period. Lagged new wells and the local Waha Hub gas price alone can therefore explain roughly one-third of Permian emission variability.

To sum up, Fig. 8 shows that trends in Permian methane emissions reflect multiple factors of variability. New well development and natural gas price are likely the most important factors during our study period. Stepping through Fig. 8a in reverse chronology, the April 2020 emission peak and subsequent COVID-19 drop were likely caused by trends in new wells, as were the October 2019 rise and December 2019 decline in emissions. The June 2019 decline may have been the result of a step change in pipeline takeaway capacity from commissioning and line packing of the 2 Bcf d−1 GCX pipeline, and the earlier March 2019 rise in emissions may have been caused by plunging Waha Hub natural gas prices. The spike in emissions observed in late 2018 may reflect low Waha Hub prices and sharply higher flaring activity. Long-term trends beneath this short-term variability are much weaker, with relatively flat emissions over the study period despite surging production, possibly due to increased gas marketing.

Table 2Multiple linear (ridge) regression of weekly Permian methane emissions onto the predictors from Fig. 8.

a The predictors are sampled as weekly values and standardized (mean = 0, standard deviation = 1). The model coefficients are thus responses to 1σ increases from the predictor means. Predictors in bold are statistically significant (p<0.05) in the full model regression. b The full model uses nine predictors, representing all the variables of Fig. 8, plus an intercept term, to predict weekly emissions. c The simplified model A uses only the statistically significant predictors from the full model, plus an intercept term. d The simplified model B uses only lagged new wells, the Waha Hub gas price, and an intercept term as predictors. e Combined monthly production in barrels of oil equivalent. f Monthly new wells for target week. g Monthly new wells with 4-week lag.

Download Print Version | Download XLSX

Figure 8Time series of weekly Permian emission estimates and different oil and gas activity variables. (a) Weekly methane emissions, reproduced from Fig. 3. We have excluded the low prior sensitivity inversion from the ensemble, as it introduces bias early in the period as shown in Fig. 3. (b) Monthly oil and gas production rates. (c) Monthly number of new wells entering production. (d) Monthly total pipeline outflows and takeaway capacity. The dashed line assumes that the Gulf Coast Express (GCX) pipeline became available for commissioning and line packing flows in May 2019 (see text). (e) Monthly number of VIIRS flaring targets. (f) Daily WTI-Cushing spot prices for oil and Henry Hub and Waha Hub spot prices for natural gas.


3.7 Temporal trend of Permian methane intensity

Figure 9 shows weekly estimates of methane intensity for the Permian, defined as total methane emitted from the oil and gas sector per unit of methane gas produced. The methane intensity is a metric for the amount of methane emitted that could have instead been taken to market and thus measures the potential for emission reductions. Values in Fig. 9 are derived from our posterior emissions (Fig. 3) and the monthly production levels of Fig. 8b, assuming 80 % methane content for Permian natural gas (Alvarez et al., 2018). We find a methane intensity of (mean ± standard deviation) 4.6 % ± 1.3 %, many times higher than the <0.2 % intensity target of the industry-based Oil & Gas Climate Initiative based on upstream processes (OGCI, 2020). There is an overall decreasing trend of about 0.02 % per week or 1 % a−1, which may reflect increased gas marketing as discussed in Sect. 3.6. The strong short-term variability we observe in intensity mirrors that in emissions and is consistent with expectations for rapidly evolving oil and gas basins (Cardoso-Saldaña and Allen, 2021).

Figure 9Weekly estimates of Permian methane intensity (oil and gas methane emissions divided by methane gas produced), assuming 80 % methane content for Permian natural gas. Monthly production is resampled to weekly values for comparison with our emission estimates. The dashed lines are moving averages of the weekly methane intensity estimates from the four sensitivity inversions of Fig. 3 varying inversion parameters γ and α. The low prior sensitivity inversion with low bias early in the study period is excluded.


Mean methane intensity in the New Mexico Permian (5.9 % ± 2.7 %) was higher than in the Texas Permian (4.2 % ± 1.2 %) during the study period, possibly due to lower production per well in New Mexico (Omara et al., 2022). These intensities are higher than the 3 %–4 % values reported by Lu et al. (2023) for the Permian in 2018–2019 based on inverse analysis of surface and GOSAT satellite data because like some of the previous TROPOMI studies discussed in Sect. 3.5, those authors used the EPA inventory with biased spatial distribution and very low emissions <1 Tg a−1 as prior estimate. An aerial survey of the New Mexico section of the Permian Basin from October 2018 through January 2020 by Chen et al. (2022) yielded total regional emissions of 194 t h−1 (1.7 Tg a−1; 1.1–2.3 Tg a−1 confidence interval), representing a very high 9.4 % (6.1 %–12.9 %) methane intensity. For the New Mexico Permian during their study period, we find an emission of (mean ± standard deviation) 1.1 ± 0.5 Tg a−1 from oil and gas and a corresponding methane intensity of 6.1 % ± 2.8 %, 35 % lower than their result but consistent within errors.

4 Conclusions

We quantified weekly methane emissions from the Permian Basin for the May 2018 to October 2020 period with a suboptimal Kalman filter inversion of TROPOMI satellite data at 0.25× 0.3125 resolution. We estimated total emissions of (mean ± standard deviation) 3.9 ± 1.0 Tg a−1 over the study period, with 3.7 ± 0.9 Tg a−1 from oil and gas, representing a mean methane intensity of 4.6 % ± 1.3 % (defined as total oil and gas methane emissions per unit of methane gas produced). These levels are nearly 4 times higher than reported in the gridded version of the US EPA greenhouse gas inventory (GHGI) for 2012, more than 6 times higher than the extrapolated GHGI for 2018, and an order of magnitude higher than methane intensity targets recently announced by the Oil & Gas Climate Initiative (OGCI). Our mean emission estimate is  20 %–35 % higher than some previous TROPOMI-based estimates (Schneising et al., 2020; Zhang et al., 2020; Liu et al., 2021), likely due to differences in prior emission estimates and methane background specification, but still within reported uncertainties of those results. For the New Mexico section of the Permian during the Chen et al. (2022) study period, we find a mean emission of 1.1 ± 0.5 Tg a−1 (6.1 % ± 2.8 % methane intensity),  35 % lower than their aerial survey reports but consistent within errors and supporting their finding of elevated intensity in that region.

We showed how TROPOMI can quantify both high-frequency and long-term variability in Permian methane emissions. Observations from the Permian Basin tower network, Scientific Aviation aircraft, and the Sentinel-2 satellites provided independent evaluation of our results. Our weekly inversions improve the fit between GEOS-Chem and the tower observations, and our estimated emissions are 20 % lower than aircraft- and tower-based emission estimates by Barkley et al. (2023) for the PermianMAP 100 km × 100 km study region of the Delaware Basin (Lyon et al., 2021). Comparing intermittent Sentinel-2 detections of a compressor station point source with our posterior emission estimates for the corresponding 0.25× 0.3125 grid cell shows roughly consistent variability.

We identified several periods of strong emission variability and investigated their causes by comparing our weekly estimates with oil and gas activity data and economic variables. Periods of rapid growth in emissions were explained by enhanced development of new wells and sharp reductions in local natural gas spot prices. Periods of rapid decline corresponded with reduced new well development (including during the 2020 COVID-19 shutdown) and rapid expansion of Permian pipeline capacity through the Gulf Coast Express project. Linear regression of our weekly emission estimates onto these and other explanatory variables confirmed that new well development and natural gas price were important (statistically significant) predictors of emission variability during the study period. The sharp emission increases inferred from TROPOMI data in October 2019 and March 2020 following spikes in new well development indicate much higher first-month emissions (5500–6300 kg d−1) per new well than previously estimated for the Barnett Shale and Eagle Ford basins using bottom-up methods (Allen et al., 2022).

We find that methane intensity in the Permian decreased over our study period, from 5 %–6 % in 2018 to 3 %–4 % in 2020. This suggests an improvement in operating practices, most apparent in the absence of a large emission increase as production increased by 50 % during the first part of the study period. The methane intensity remained an order of magnitude larger than the 0.2 % goal set by industry (OGCI) for oil and gas production by 2025.

The recently passed US Inflation Reduction Act includes a methane emissions fee for oil and gas companies and requires the US EPA to update the Greenhouse Gas Reporting Program based on measurement data. Our work demonstrates the potential role of satellite observations in validating reported emissions under such policies. This is a step towards near-real-time satellite monitoring of methane emissions as a strategy to evaluate bottom-up inventories, monitor progress on climate goals, and identify sudden increases in emissions. Such monitoring can help improve our understanding of the factors controlling methane emissions and how emissions respond to economic forcing and policy changes.

Data availability

The TROPOMI methane data are available at (last access: 11 October 2022) (DOI:, Lorente et al., 2021b). The Permian Basin tower network data are available at (Monteiro et al., 2021). The Sentinel-2 satellite data are available through the Copernicus Open Access Hub at (ESA, 2022). Our weekly emission estimates are available at (Varon, 2022).

Author contributions

DJV and DJJ conceptualized the study. DJV, DJJ, MS, LS, DP, HN, and ZQ contributed to the method development. SP performed the Sentinel-2 point source analysis. DJV led the data analysis and interpretation of results with input from all authors. DJV wrote the original draft, and all authors reviewed and edited the paper.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


This work was supported by the Environmental Defense Fund and by NASA's Carbon Monitoring System. We thank Jack Winters and James Robinson of S&P Global Platts for sharing natural gas throughput and pricing data. Alba Lorente acknowledges funding from the TROPOMI national program through NSO.

Financial support

This research has been supported by the National Aeronautics and Space Administration (Carbon Monitoring System – NASA – grant no. 80NSSC21K1057) and the Environmental Defense Fund (grant no. 1063-000000-10800).

Review statement

This paper was edited by Gabriele Stiller and reviewed by two anonymous referees.


Allen, D. T., Cardoso-Saldaña, F. J., Kimura, Y., Chen, Q., Xiang, Z., Zimmerle, D., Bell, C., Lute, C., Duggan, J., and Harrison, M.: A Methane Emission Estimation Tool (MEET) for predictions of emissions from upstream oil and gas well sites with fine scale temporal and spatial resolution: Model structure and applications, Sci. Total Environ., 829, 154277,, 2022. 

Alvarez, R. A., Zavala-Araiza, D., Lyon, D. R., Allen, D. T., Barkley, Z. R., Brandt, A. R., Davis, K. J., Herndon, S. C., Jacob, D. J., Karion, A., Kort, E. A., Lamb, B. K., Lauvaux, T., Maasakkers, J. D., Marchese, A. J., Omara, M., Pacala, S. W., Peischl, J., Robinson, A. L., Shepson, P. B., Sweeney, C., Townsend-Small, A., Wofsy, S. C., and Hamburg, S. P.: Assessment of methane emissions from the U.S. oil and gas supply chain, Science, 361, 186–188,, 2018. 

Barkley, Z., Davis, K., Miles, N., Richardson, S., Deng, A., Hmiel, B., Lyon, D., and Lauvaux, T.: Quantification of oil and gas methane emissions in the Delaware and Marcellus basins using a network of continuous tower-based measurements, Atmos. Chem. Phys., 23, 6127–6144,, 2023. 

Bloom, A. A., Bowman, K. W., Lee, M., Turner, A. J., Schroeder, R., Worden, J. R., Weidner, R., McDonald, K. C., and Jacob, D. J.: A global wetland methane emissions and uncertainty dataset for atmospheric chemical transport models (WetCHARTs version 1.0), Geosci. Model Dev., 10, 2141–2156,, 2017. 

Brasseur, G. P. and Jacob, D. J.: Modeling of Atmospheric Chemistry, Cambridge University Press, Cambridge, UK, ISBN 9781316544754, 2017. 

Cardoso-Saldaña, F. J. and Allen, D. T.: Projecting the Temporal Evolution of Methane Emissions from Oil and Gas Production Sites, Environ. Sci. Technol., 54, 14172–14181,, 2020. 

Cardoso-Saldaña, F. J. and Allen, D. T.: Projecting the Temporal Evolution of Methane Emissions from Oil and Gas Production Basins, Environ. Sci. Technol., 55, 2811–2819,, 2021. 

Chen, Y., Sherwin, E. D., Berman, E. S. F., Jones, B. B., Gordon, M. P., Wetherley, E. B., Kort, E. A., and Brandt, A. R.: Quantifying Regional Methane Emissions in the New Mexico Permian Basin with a Comprehensive Aerial Survey, Environ. Sci. Technol., 56, 4317–4323,, 2022. 

Chen, Z., Jacob, D. J., Gautam, R., Omara, M., Stavins, R. N., Stowe, R. C., Nesser, H. O., Sulprizio, M. P., Lorente, A., Varon, D. J., Lu, X., Shen, L., Qu, Z., Pendergrass, D. C., and Hancock, S.: Satellite quantification of methane emissions and oil/gas methane intensities from individual countries in the Middle East and North Africa: implications for climate action, EGUsphere [preprint],, 2023. 

Cusworth, D. H., Duren, R. M., Thorpe, A. K., Olson-Duvall, W., Heckler, J., Chapman, J. W., Eastwood, M. L., Helmlinger, M. C., Green, R. O., Asner, G. P., Dennison, P. E., and Miller, C. E.: Intermittency of large methane emitters in the Permian Basin, Environ. Sci. Tech. Let., 8, 567–573,, 2021a. 

Cusworth, D. H., Duren, R., Thorpe, A. K., Pandey, S., Aben, I., Jervis, D., Varon, D., Jacob, D. J., Randles, C. A., Gautam, R., Omara, M., Schade, G., Dennison, P. E., Frankenberg, C., Gordon, D., Lopinto, E., and Miller, C. E.: Multi-satellite imaging of a gas well blowout provides new insights for methane monitoring, Geophys. Res. Lett., 48, e2020GL090864,, 2021b. 

Cusworth, D. H., Thorpe, A. K., Ayasse, A., Stepp, D., Heckler, J., Asner, G. P., Miller, C. E., Yadav, V., Chapman, J. W., Eastwood, M. L., Green, R. O., Hmiel, B., Lyon, D. R., and Duren, R. M.: Strong methane point sources contribute a disproportionate fraction of total emissions across multiple basins in the United States, P. Natl. Acad. Sci. USA, 119, e2202338119,, 2022. 

de Gouw, J. A., Veefkind, J. P., Roosenbrand, E., Dix, B., Lin, J. C., Landgraf, J., and Levelt, P. F.: Daily Satellite Observations of Methane from Oil and Gas Production Regions in the United States, Sci. Rep., 10, 1379,, 2020. 

Ehret, T., De Truchis, A., Mazzolini, M., Morel, J.-M., d'Aspremont, A., Lauvaux, T., Duren, R., Cusworth, D., and Facciolo, G.: Global Tracking and Quantification of Oil and Gas Methane Emissions from Recurrent Sentinel-2 Imagery, Environ. Sci. Tech. Let., 56, 10517–10529,, 2022. 

Elvidge, C. D., Zhizhin, M., Baugh, K., Hsu, F.-C., and Ghosh, T.: Methods for Global Survey of Natural Gas Flaring from Visible Infrared Imaging Radiometer Suite Data, Energies, 9, 14,, 2016. 

Enverus: Drillinginfo,, last access: 15 December 2021. 

European Space Agency (ESA): Copernicus Open Access Hub,, last access: 11 October 2022. 

Federal Reserve Bank of Dallas (FRBD): Energy in the Eleventh District: Permian Basin: Oil production,, last access: 12 October 2022. 

Hasekamp, O., Lorente, A., Hu, H., Butz, A., aan de Brugh, J., and Landgraf, J.: Algorithm Theoretical Baseline Document for Sentinel-5 Precursor Methane Retrieval, SRON-S5P-LEV2-RP-001, (Issue 2.4.0), European Space Agency technical document, European Space Agency, 1–64, (last access: 21 June 2023), 2022. 

Hastie T., Tibshirani R., and Friedman J.: The elements of statistical learning: data mining, inference, and prediction, 2nd edn., Springer Science & Business Media, ISBN 9781282126749, 2009. 

Heald, C. L., Jacob, D. J., Jones, D. B. A., Palmer, P. I., Logan, J. A., Streets, D. G., Sachse, G. W., Gille, J. C., Hoffman, R. N., and Nehrkorn, T.: Comparative inverse analysis of satellite (MOPITT) and aircraft (TRACE-P) observations to estimate Asian sources of carbon monoxide, J. Geophys. Res., 109, D23306,, 2004. 

Houweling, S., Bergamaschi, P., Chevallier, F., Heimann, M., Kaminski, T., Krol, M., Michalak, A. M., and Patra, P.: Global inverse modeling of CH4 sources and sinks: an overview of methods, Atmos. Chem. Phys., 17, 235–256,, 2017. 

Hu, H., Landgraf, J., Detmers, R., Borsdorff, T., Aan de Brugh, J., Aben, I., Butz, A., and Hasekamp, O.: Toward global mapping of methane with TROPOMI: First results and intersatellite comparison to GOSAT, Geophys. Res. Lett., 45, 3682–3689,, 2018. 

International Energy Agency: Global Methane Tracker 2022, Paris, IEA,, last access: 11 October 2022. 

IPCC: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, in press,, 2021. 

Irakulis-Loitxate, I., Guanter, L., Liu, Y.-N., Varon, D. J., Maasakkers, J. D., Zhang, Y., Thorpe, A. K., Duren, R. M., Frankenberg, C., Lyon, D., Hmiel, B., Cusworth, D. H., Zhang, Y., Segl, K., Gorron, J., Sanchez-García, E., Sulprizio, M. P., Cao, K., Zhu, H., Liang, J., Li, X., Aben, I., and Jacob, D. J.: Satellite-based survey of extreme methane emissions in the Permian Basin , Sci. Adv., 7, eabf4507,, 2021. 

Jacob, D. J., Turner, A. J., Maasakkers, J. D., Sheng, J., Sun, K., Liu, X., Chance, K., Aben, I., McKeever, J., and Frankenberg, C.: Satellite observations of atmospheric methane and their value for quantifying methane emissions, Atmos. Chem. Phys., 16, 14371–14396,, 2016. 

Jacob, D. J., Varon, D. J., Cusworth, D. H., Dennison, P. E., Frankenberg, C., Gautam, R., Guanter, L., Kelley, J., McKeever, J., Ott, L. E., Poulter, B., Qu, Z., Thorpe, A. K., Worden, J. R., and Duren, R. M.: Quantifying methane emissions from the global scale down to point sources using satellite observations of atmospheric methane, Atmos. Chem. Phys., 22, 9617–9646,, 2022. 

Janssens-Maenhout, 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,, 2019. 

Lin, J. C., Bares, R., Fasoli, B., Garcia, M., Crosman, E., and Lyman, S.: Declining methane emissions and steady, high leakage rates observed over multiple years in a western US oil/gas production basin, Scientific Reports, 11, 22291,, 2021. 

Liu, M., van der A, R., van Weele, M., Eskes, H., Lu, X., Veefkind, P., de Laat, J., Kong, H., Wang, J., Sun, J., Ding, J., Zhao, Y., and Weng, H.: A New Divergence Method to Quantify Methane Emissions Using Observations of Sentinel-5P TROPOMI, Geophys. Res. Lett., 48, e2021GL094151,, 2021. 

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,, 2021a. 

Lorente, A., Borsdorff, T., aan de Brugh, J., Landgraf, J., and Hasekamp, O.: SRON S5P – RemoTeC scientific TROPOMI XCH4 dataset, Zenodo [data set],, 2021b. 

Lu, X., Jacob, D. J., Zhang, Y., Maasakkers, J. D., Sulprizio, M. P., Shen, L., Qu, Z., Scarpelli, T. R., Nesser, H., Yantosca, R. M., Sheng, J., Andrews, A., Parker, R. J., Boesch, H., Bloom, A. A., and Ma, S.: Global methane budget and trend, 2010–2017: complementarity of inverse analyses using in situ (GLOBALVIEWplus CH4 ObsPack) and satellite (GOSAT) observations, Atmos. Chem. Phys., 21, 4637–4657,, 2021. 

Lu, X., Jacob, D. J., Wang, H., Maasakkers, J. D., Zhang, Y., Scarpelli, T. R., Shen, L., Qu, Z., Sulprizio, M. P., Nesser, H., Bloom, A. A., Ma, S., Worden, J. R., Fan, S., Parker, R. J., Boesch, H., Gautam, R., Gordon, D., Moran, M. D., Reuland, F., Villasana, C. A. O., and Andrews, A.: Methane emissions in the United States, Canada, and Mexico: evaluation of national methane emission inventories and 2010–2017 sectoral trends by inverse analysis of in situ (GLOBALVIEWplus CH4 ObsPack) and satellite (GOSAT) atmospheric observations, Atmos. Chem. Phys., 22, 395–418,, 2022. 

Lu, X., Jacob, D. J., Zhang, Y., Shen, L., Sulprizio, M. P., Maasakkers, J. D., Varon, D. J., Qu, Z., Chen, Z., Hmiel, B., Parker, R. J., Boesch, H., Wang, H., He, C., and Fan, S.: Observation-derived 2010–2019 trends in methane emissions and intensities from US oil and gas fields tied to activity metrics, P. Natl. Acad. Sci. USA, 120, e2217900120,, 2023. 

Lyon, D. R., Hmiel, B., Gautam, R., Omara, M., Roberts, K. A., Barkley, Z. R., Davis, K. J., Miles, N. L., Monteiro, V. C., Richardson, S. J., Conley, S., Smith, M. L., Jacob, D. J., Shen, L., Varon, D. J., Deng, A., Rudelis, X., Sharma, N., Story, K. T., Brandt, A. R., Kang, M., Kort, E. A., Marchese, A. J., and Hamburg, S. P.: Concurrent variation in oil and gas methane emissions and oil price during the COVID-19 pandemic, Atmos. Chem. Phys., 21, 6605–6626,, 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,, 2016. 

McNorton, J., Bousserez, N., Agustí-Panareda, A., Balsamo, G., Cantarello, L., Engelen, R., Huijnen, V., Inness, A., Kipling, Z., Parrington, M., and Ribas, R.: Quantification of methane emissions from hotspots and during COVID-19 using a global atmospheric inversion, Atmos. Chem. Phys., 22, 5961–5981,, 2022. 

Monteiro, V., Miles, N. L., Richardson, S. J., Barkley, Z. R., Haupt, B. J., and Davis, K. J.: Permian Basin: in-situ tower greenhouse gas data, Penn State [data set],, 2021. 

Monteiro, V. C., Miles, N. L., Richardson, S. J., Barkley, Z., Haupt, B. J., Lyon, D., Hmiel, B., and Davis, K. J.: Methane, carbon dioxide, hydrogen sulfide, and isotopic ratios of methane observations from the Permian Basin tower network, Earth Syst. Sci. Data, 14, 2401–2417,, 2022. 

OGCI: Methane Intensity Target,, (last access: 7 September 2022), 2020. 

Omara M., Zavala-Araiza D., Lyon D. R., Hmiel B., Roberts K. A., and Hamburg S. P.: Methane emissions from us low production oil and natural gas well sites, Nat. Commun., 13, 2085,, 2022. 

Pandey, S., Gautam, R., Houweling, S., Denier van der Gon, H., Sadavarte, P., Borsdorff, T., Hasekamp, O., Landgraf, J., Tol, P., van Kempen, T., Hoogeveen, R., van Hees, R., Hamburg, S. P., Maasakkers, J. D., and Aben, I.: Satellite observations reveal extreme methane leakage from a natural gas well blowout, P. Natl. Acad. Sci. USA, 116, 23676–23681,, 2019. 

Plant, G., Kort, E. A., Brandt, A. R., Chen, Y., Fordice, G., Gorchov Negron, A. M., Schwietzke, S., Smith, M., and Zavala-Araiza, D.: Inefficient and unlit natural gas flares both emit large quantities of methane, Science, 377, 1566,, 2022. 

Qu, Z., Jacob, D. J., Shen, L., Lu, X., Zhang, Y., Scarpelli, T. R., Nesser, H., Sulprizio, M. P., Maasakkers, J. D., Bloom, A. A., Worden, J. R., Parker, R. J., and Delgado, A. L.: Global distribution of methane emissions: a comparative inverse analysis of observations from the TROPOMI and GOSAT satellite instruments, Atmos. Chem. Phys., 21, 14159–14175,, 2021. 

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding: Theory and Practice, World Scientific, River Edge, USA, ISBN 9789810227401, 2000. 

Rystad Energy: Permian Basin Flaring Outlook: Condensed Report, (last access: 26 October 2022), 2021. 

S&P Global Platts: Homepage, (last access: 10 October, 2022), 2019. 

S&P Global Platts: Platts North American Natural Gas Analytics,, last access: 27 January 2022. 

Sánchez-García, E., Gorroño, J., Irakulis-Loitxate, I., Varon, D. J., and Guanter, L.: Mapping methane plumes at very high spatial resolution with the WorldView-3 satellite, Atmos. Meas. Tech., 15, 1657–1674,, 2022. 

Scarpelli, T. R., Jacob, D. J., Octaviano Villasana, C. A., Ramirez Hernandez, I. F., Cardenas Moreno, P. R., Cortes Alfaro, E. A., Garcia Garcia, M. A., and Zavala-Araiza, D.: A gridded inventory of anthropogenic methane emissions from Mexico based on Mexico's National Inventory of Greenhouse Gases and Compounds, Environ. Res. Lett., 15, 105015,, 2020. 

Schneising, O., Buchwitz, M., Reuter, M., Bovensmann, H., Burrows, J. P., Borsdorff, T., Deutscher, N. M., Feist, D. G., Griffith, D. W. T., Hase, F., Hermans, C., Iraci, L. T., Kivi, R., Landgraf, J., Morino, I., Notholt, J., Petri, C., Pollard, D. F., Roche, S., Shiomi, K., Strong, K., Sussmann, R., Velazco, V. A., Warneke, T., and Wunch, D.: A scientific algorithm to simultaneously retrieve carbon monoxide and methane from TROPOMI onboard Sentinel-5 Precursor, Atmos. Meas. Tech., 12, 6771–6802,, 2019. 

Schneising, O., Buchwitz, M., Reuter, M., Vanselow, S., Bovensmann, H., and Burrows, J. P.: Remote sensing of methane leakage from natural gas and petroleum systems revisited, Atmos. Chem. Phys., 20, 9169–9182,, 2020. 

Shen, L., Zavala-Araiza, D., Gautam, R., Omara, M., Scarpelli, T., Sheng, J., Sulprizio, M. P., Zhuang, J., Zhang, Y., Lu, X., Hamburg, S. P., and Jacob, D. J.: Unravelling a large methane emission discrepancy in Mexico using satellite observations, Remote Sens. Environ., 260, 112461,, 2021. 

Shen, L., Gautam, R., Omara, M., Zavala-Araiza, D., Maasakkers, J. D., Scarpelli, T. R., Lorente, A., Lyon, D., Sheng, J., Varon, D. J., Nesser, H., Qu, Z., Lu, X., Sulprizio, M. P., Hamburg, S. P., and Jacob, D. J.: Satellite quantification of oil and natural gas methane emissions in the US and Canada including contributions from individual basins, Atmos. Chem. Phys., 22, 11203–11215,, 2022. 

Todling, R. and Cohn, S. E.: Suboptimal schemes for atmospheric data assimilation based on the Kalman filter, Mon. Weather Rev., 122, 2530–2557,<2530:SSFADA>2.0.CO;2, 1994. 

U. S. Energy Information Administration (EIA): Data for Petroleum and Other Liquids, (last access: 15 December 2021), 2022a. 

U. S. Energy Information Administration (EIA): Data for Natural Gas, (last access: 15 December 2021), 2022b. 

Varon, D.: Replication Data for Continuous weekly monitoring of methane emissions from the Permian Basin by inversion of TROPOMI satellite observations, V1, Harvard Dataverse [data set],, 2022. 

Varon, D. J., Jacob, D. J., McKeever, J., Jervis, D., Durak, B. O. A., Xia, Y., and Huang, Y.: Quantifying methane point sources from fine-scale satellite observations of atmospheric methane plumes, Atmos. Meas. Tech., 11, 5673–5686,, 2018. 

Varon, D. J., Jervis, D., McKeever, J., Spence, I., Gains, D., and Jacob, D. J.: High-frequency monitoring of anomalous methane point sources with multispectral Sentinel-2 satellite observations, Atmos. Meas. Tech., 14, 2771–2785,, 2021. 

Varon, D. J., Jacob, D. J., Sulprizio, M., Estrada, L. A., Downs, W. B., Shen, L., Hancock, S. E., Nesser, H., Qu, Z., Penn, E., Chen, Z., Lu, X., Lorente, A., Tewari, A., and Randles, C. A.: Integrated Methane Inversion (IMI 1.0): a user-friendly, cloud-based facility for inferring high-resolution methane emissions from TROPOMI satellite observations, Geosci. Model Dev., 15, 5787–5805,, 2022. 

Veefkind, J. P., Aben, I., McMullan, K., Förster, H., De Vries, J., Otter, G., Claas, J., Eskes, H. J., De Haan, J. F., Kleipool, Q., and Van Weele, M.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83,, 2012.  

Worden, J. R., Cusworth, D. H., Qu, Z., Yin, Y., Zhang, Y., Bloom, A. A., Ma, S., Byrne, B. K., Scarpelli, T., Maasakkers, J. D., Crisp, D., Duren, R., and Jacob, D. J.: The 2019 methane budget and uncertainties at 1° resolution and each country through Bayesian integration Of GOSAT total column methane data and a priori inventory estimates, Atmos. Chem. Phys., 22, 6811–6841,, 2022. 

Yu, X., Millet, D. B., and Henze, D. K.: How well can inverse analyses of high-resolution satellite data resolve heterogeneous methane fluxes? Observing system simulation experiments with the GEOS-Chem adjoint model (v35), Geosci. Model Dev., 14, 7775–7793,, 2021. 

Zhang, Y., Gautam, R., Pandey, S., Omara, M., Maasakkers, J. D., Sadavarte, P., Lyon, D., Nesser, H., Sulprizio, M. P., Varon, D. J., Zhang, R., Houweling, S., Zavala-Araiza, D., Alvarez, R. A., Lorente, A., Hamburg, S. P., Aben, I., and Jacob, D. J.: Quantifying methane emissions from the largest oil-producing basin in the United States from space, Sci. Adv., 6, 17,, 2020. 

Executive editor
Methane is a potent greenhouse gas with a global warming potential much greater than carbon dioxide. In order to understand and limit its global warming effect, it is important to have have observation systems to estimate its anthropogenic emissions. This paper analyses in a novel way TROPOMI satellite observations to track methane emissions from the largest oil production basin in the United States over a 2-year period. The analysis shows that emission variability and trends are driven by multiple factors, among which new well development and natural gas spot price are the most significant ones. The work is an excellent demonstration of the potential of satellite observations for near-real-time monitoring of methane emissions. It opens a broad range of applications, helping science and policy to understand and mitigate climate change.
Short summary
We use TROPOMI satellite observations to quantify weekly methane emissions from the US Permian oil and gas basin from May 2018 to October 2020. We find that Permian emissions are highly variable, with diverse economic and activity drivers. The most important drivers during our study period were new well development and natural gas price. Permian methane intensity averaged 4.6 % and decreased by 1 % per year.
Final-revised paper