Articles | Volume 21, issue 16
Research article
19 Aug 2021
Research article |  | 19 Aug 2021

Shipborne measurements of methane and carbon dioxide in the Middle East and Mediterranean areas and the contribution from oil and gas emissions

Jean-Daniel Paris, Aurélie Riandet, Efstratios Bourtsoukidis, Marc Delmotte, Antoine Berchet, Jonathan Williams, Lisa Ernle, Ivan Tadic, Hartwig Harder, and Jos Lelieveld

The increase of atmospheric methane (CH4) and carbon dioxide (CO2), two of the main anthropogenic greenhouse gases, is largely driven by fossil sources. Sources and sinks remain insufficiently characterized in the Mediterranean and Middle East areas, where very few in situ measurements are available. We measured the atmospheric mixing ratios of CH4 and CO2 by ship in the region in July and August 2017. High mixing ratios were observed over the Suez Canal, Red Sea and Arabian Gulf, while generally lower mixing ratios were observed over the Gulf of Aden and Gulf of Oman. We probe the origin of the CO2 and CH4 excess mixing ratio by using correlations with light alkanes and through the use of a Lagrangian model coupled to two different emission inventories of anthropogenic sources. We find that the CO2 and especially the CH4 enhancements are mainly linked to nearby oil and gas (OG) activities over the Arabian Gulf and a mixture of other sources over the Red Sea. The isomeric ratio of pentane is shown to be a useful indicator of the OG component of atmospheric CH4 at the regional level. Upstream emissions linked to oil in the northern Arabian Gulf seem to be underestimated, while gas-related emissions in the southern Gulf are overestimated in our simulations. Our results highlight the need for improvement of inventories in the area to better characterize the changes in magnitude and the complex distribution of the OG sources in the Middle East.

1 Introduction

Methane (CH4) and carbon dioxide (CO2) are potent anthropogenic greenhouse gases (GHGs). The CH4 atmospheric mole fraction has increased by 150 % since the preindustrial era (Saunois et al., 2020). Over a 100-year horizon, CH4 has a global warming potential 28 times larger than CO2. Roughly half of CH4 sources are of natural origin (mainly from wetlands, with contributions from fires and geologic sources). The remainder is anthropogenic, mainly linked to fossil fuels, agriculture (including enteric fermentation in ruminants, manure management and rice paddies) and waste management. The increasing anthropogenic emissions are driven equally by fossil fuel sources and agricultural sources (Jackson et al., 2020).

CO2 concentration in the atmosphere has increased by 47 % since the preindustrial era and reached 407.38 ± 0.10 ppm in 2018 (Dlugokencky and Tans, 2018; Le Quéré et al., 2018; Friedlingstein et al., 2019). Its increase is caused primarily by the use of fossil fuel, cement production and land use change and is partly mitigated by uptake in terrestrial ecosystems and the ocean. Over the last decade (2009–2018) CO2 sources have been dominated by fossil fuels (9.5 ± 0.5 Gt C yr−1), with a significant source from land use change (1.5 ± 0.7 Gt C yr−1). CO2 is removed from the atmosphere by sinks in the ocean (2.5 ± 0.6 Gt C yr−1) and land ecosystems (3.2 ± 0.6 Gt C yr−1; Friedlingstein et al., 2020). Atmospheric CO2 is expected to continue its growth despite widespread adoption of climate policies as oil and gas consumption trends suggest a continued increase in fossil CO2 emission (Jackson et al., 2019).

The 2015 Paris Agreement has set the objective to limit global temperatures below 1.5 C. To reach this objective, a corridor of compliant emission pathways has been designed that require strong reductions not only in CO2 emissions but also in anthropogenic CH4 emissions (Jones et al., 2018; Nisbet et al., 2020). Regarding CH4, a wide array of measures is necessary (Nisbet et al., 2020), with a potentially high mitigation impact in the energy, agricultural and waste management sectors (Rogelj et al., 2018). Rapid and efficient mitigation options are readily available for methane emissions. A total of 80 % of economically feasible abatements reside in the oil and gas (OG) sector (Ocko et al., 2021). IEA (2020) reviewed a number of efficient and cost-effective abatement options for the OG sector and found that 40 % of OG methane emissions could be avoided at a zero net cost.

Fossil fuel production and use is responsible for the release of 112 Mt CH4 yr−1 to the atmosphere, representing 33 % of the total anthropogenic emission of CH4 (Saunois et al., 2020). Considering fossil fuel emissions alone, 68 % of the emissions are linked to OG, while the rest are associated with coal mining. Emission of CH4 arises at each step from the production site to the consumption site. For the oil industry, methane emissions occur essentially as indirect emissions, i.e., as venting or incomplete combustion during flaring and during transport and refining (IEA, 2020). For the gas industry, emission of CH4 occurs at all stages of the life cycle as fugitive emissions (leaks from valves, connectors, and compressors and intentional venting) but also as incomplete combustion during flaring (GIE-MARCOGAZ, 2019). Alvarez et al. (2018) found that in the United States a large fraction of the net emission is associated with production, transport and processing.

Large uncertainties remain that are associated with the magnitude and spatial and temporal distribution of CH4 sources (Saunois et al., 2020). After a pause between 2000 and 2007, CH4 in the atmosphere has resumed its increase. Constraints from current available observations on the respective contribution to its sources and sinks do not allow a definitive explanation of this pattern (Saunois et al., 2017, 2020; Turner et al., 2019; Nisbet et al., 2020). Schwietzke et al. (2016) highlighted the strong underestimation of OG emissions in current inventories (see also Saunois et al., 2020, and references therein). Based on ice core 13CH4 measurements, Hmiel et al. (2020) found that fossil fuel methane sources could be underestimated by as much as 25 %–40 %. The largest uncertainty that remains associated with CH4 emissions is found during the extraction of OG, with global estimates ranging from 46 to 98 Mt CH4 yr−1 (Höglund-Isaksson, 2012). This uncertainty is critically related to the diversity of country-specific and site-specific emission factors used. In the US, emissions associated with OG production were underestimated by a factor 2 (Alvarez et al., 2018), which is mainly linked to underestimation in upstream (production) emissions. Furthermore, the spatial and temporal distribution of sources is poorly known. In the Barnett Shale area (USA), Zavala-Araiza et al. (2015) found that 10 % of the OG facilities accounted for 90 % of the emissions. Assessment of regional emission distributions can be provided by mobile measurements targeting facilities.

CO2 emissions linked to OG production and use represent 53 % of total CO2 emissions over the period 2008–2017 (Le Quéré et al., 2018; not including land use change). Unlike CH4, CO2 emissions are to a large extent induced by fuel usage for energy consumption rather than production. In the near future, CO2 emissions are expected to continue increasing, driven by a strong demand in Asia (Jackson et al., 2019).

A large fraction of global OG extraction occurs in the eastern Mediterranean and Middle East region (EMME). The Middle East is the main crude oil production region (32 % of world total), with 24.16 million barrels per day, and 17,5 % of global natural gas production, with 701.12 billion standard cubic meters produced in 2019 (OPEC, 2020). The main OG fields are located in the Arabian Gulf and neighboring countries, as well as in the Gulf of Suez and the Nile Delta. Gas fields have recently been discovered in the Levantine Sea. CH4 emissions reported by the EDGAR inventory for countries neighboring the Mediterranean Sea and Middle East amount to 16 Mt CH4 yr−1 (6.0 % of global emissions) and 21 Mt CH4 yr−1 (8.3 % of global emissions) (Janssens-Maenhout et al., 2017), respectively. Across the different inventories reported in Saunois et al. (2020) for the Middle East, the spread represents 18 % of the mean emission intensity, reflecting a significant uncertainty on country-level emissions. Increasing CH4 emissions in the Middle East and Africa have been proposed as a contribution to the post-2007 CH4 increase (McNorton et al., 2018).

According to the inventory data compiled by Friedlingstein et al. (2019), Mediterranean countries and Middle East CO2 emissions represent respectively 6.8 % and 5.6 % of global emissions. CO2 emissions for the EMME have increased from 662 to 831 Mt C yr−1 over the period 2009–2018, essentially driven by fossil fuel use in Middle Eastern countries. This represent an annual growth rate of 2.5 %, which is significantly higher than the 1.3 % yr−1 increase rate of global emissions. CO2 emissions from the Middle East mostly arise from domestic consumption, while the emission transfers linked to international trade are negligible (2 ± 7 Mt C yr−1; Peters et al., 2012). However, inventories may underestimate emissions from the main urban centers of the Middle East by a factor 2 (Yang et al., 2020).

The EMME region is a transitional zone between midlatitude climates and subtropical areas, located in the high-pressure subtropical ridge (Lelieveld et al., 2012). The northern part of the EMME is under a westerly regime, with the eastern Mediterranean located in the outflow from European air masses in the lower troposphere (Lelieveld et al., 2002), while the southern part experiences trade winds. The EMME region includes a large amount of both desert and densely populated areas. Despite the region's important contribution to anthropogenic GHG emission, very few atmospheric measurements of the distribution of GHG are available in the area (Ricaud et al., 2018), limiting the possibility to reduce uncertainties in regional emissions rates. The Middle East remains undersampled, especially by the in situ surface networks (Ciais et al., 2010). It also offers fewer cloudy days than other midlatitude locations, thus enhancing the potential for passive satellite measurements (Yang et al., 2020).

What are the typical CO2 and CH4 mixing ratios around the Arabian Peninsula? What are the main drivers for the observed variability, and can we link this variability to sources? To what extent can we confirm or inform the inventories based on the measurements and more specifically the OG component of these inventories? This paper aims to obtain a better understanding the drivers of variability for these species in the area and to relate this variability to regional sources.

To address these questions, CO2 and CH4 mixing ratio measurements were performed during a ship cruise across the Mediterranean Sea, Red Sea and Arabian Gulf in the summer of 2017. Large-scale, shipborne measurements of CH4 and CO2 atmospheric concentrations have been previously performed for a wide variety of purposes, including measuring emissions from OG platforms (see, e.g., Yacovitch et al., 2020), assessing nation-wide emissions in the UK and Ireland (Helfter et al., 2019), and disentangling regional sources in the Arctic (Berchet et al., 2020).

Here we build on the work of Bourtsoukidis et al. (2019), who characterized light alkanes in relation to the various hydrocarbon sources in the region. We compare these results to a simulation of the anthropogenic component of CO2 and CH4 excess during the cruise using a Lagrangian particle dispersion model.

Section 2 details the AQABA (Air Quality and climate change in the Arabian Basin) campaign, as well as the measurement and modeling methodology. Section 3 presents measurement data and discusses how these data compare with simulations of CH4 and CO2 enhancements linked to anthropogenic activities.

2 Methods

2.1 Campaign and platform

The AQABA campaign took place in the summer of 2017 from 24 June to 3 September. The scientific purpose of the campaign was to investigate the atmospheric composition and chemical processes over the Arabian Basin (Pfannerstill et al., 2019; Wang et al., 2020; Bourtsoukidis et al., 2019, 2020; Tadic et al., 2020; Celik et al., 2020). The ship both departed from and returned to La Seyne-sur-Mer, near Toulon, France. Figure 1 shows the ship's route and calling ports. The ship traveled through the Suez Canal and the Red Sea, then around the Arabian Peninsula and through the Arabian Gulf (also called the Persian Gulf) to Kuwait, where it anchored for 3 d (from 31 July to 3 August) at the port. The ship eventually returned and reached La Seyne-sur-Mer by approximately the same route on 2 September. Several calls occurred at various ports on the way. The ship used was the R/V Kommandor Iona, a UK-based, 76 m long dynamic positioning research and survey vessel. The mean speed of the vessel during the campaign was 3.4 ± 1.8 m s−1.

Figure 1Ship track during the campaign. Regional limits are indicated. Background map: © OpenStreetMap contributors 2021. Distributed under the Open Data Commons Open Database License (ODbL) v1.0.

2.2 Measurements

To minimize contamination of measurements by the ship's own emissions, the measurement inlet for CO2 and CH4 was located at the front of the ship, 15 m above sea level. Air was drawn through 15 m of Synflex tubing using a 2 L min−1 flush pump to prevent degradation of the system response time. The CO2 and CH4 mole fractions were measured in situ using a cavity ring-down spectroscopy analyzer (Picarro G2401; Santa Clara, USA) located in an air-conditioned container on the front deck. CO measurements were also collected. The air was not dried prior to analysis, and water vapor effects on CO2 and CH4 were determined specifically for this instrument and corrected (Hazan et al., 2016). A filter and a peristaltic pump ensured that no dust or liquid water entered the analyzer. Four different calibration cylinders of compressed air bracketing typical ambient concentrations were injected into the analyzer before departure, every 15 d during the trip and upon return to Toulon for calibration and quality assurance. Calibration gases were prepared in Luxfer aluminum cylinders. Calibration gas concentrations ranged between 377.63 ± 0.01 ppm and 474.84 ± 0.01 ppm for CO2, between 1791.8 ± 0.1 and 2791.0 ±​​​​​​​ 0.2 ppb for CH4, and between 127 ± 1 and 346 ± 1 ppb for CO. The injection sequence consisted of four 20 min injections of each of the four gases. An additional target gas of ambient concentration was injected twice daily to assess measurement accuracy. Precision of the G2401 (expressed as continuous measurement repeatability) is typically better than 0.03 ppm CO2, 0.3 ppb CH4 and 8 ppb CO (Yver-Kwok et al., 2015). Calibrations were quite stable and consistently exhibited a within-run standard deviation of 0.03 ppm CO2, 0.3 ppb CH4 and 10 ppb CO with little change during the trip. The data have been processed and quality-controlled following ICOS (Integrated Carbon Observing System) standard procedure (Hazan et al., 2016), including the propagation of the calibration- and threshold-based filters. The concentrations measured from calibration cylinders between two calibration sequences show a mean drift of 0.05 ppm and 0.5 ppb CH4, significantly below the drifts typically observed at fixed observatories (Hazan et al., 2016). Target injections showed a small residual bias (after calibration) below 0.05 ppm CO2 and 0.3 ppb CH4. The measured target values for CO vary within ±10 ppb SD. The processed data are reported as 1 min averages.

Alkanes were measured using a gas chromatography–flame ionization detector (GC-FID; Bourtsoukidis et al., 2019). GC-FID sampling intervals were 10–30 min, with uncertainties ranging between 5 % and 10 %. A detailed description of the GC-FID measurements can be found in Bourtsoukidis et al. (2019).

Meteorological parameters such as wind speed and direction, pressure, temperature, and GPS position and course were acquired from a meteorological station (Shipborne European Common Automatic Weather Station, EUCAWS) at the starboard side of the front deck of the ship.

2.3 Stack contamination data filter

As the sampling inlet was situated at the front of the boat, measurements can be occasionally influenced by the ship's stack emissions. We assume that this influence mainly depends on the relative wind direction, where a tailwind is likely to bring smoke from the boat's chimney back onto the instrument and thus contaminate the measurements. We flagged 1 min measured concentrations using a binary index that indicated a possible presence or absence of stack contamination as in Tadic et al. (2020). Different values of angle sector around the stern were tested to filter data using relative wind direction and speed. A total of 16 % of the data were flagged as potentially contaminated by the ship and not considered in the following analysis unless specified otherwise.

2.4 Lagrangian modeling

Atmospheric transport was investigated using the FLEXible PARTicle dispersion model (FLEXPART) v9 Lagrangian particle dispersion model (Pisso et al., 2019) in receptor-oriented mode (Seibert and Frank, 2004). FLEXPART calculates the trajectories of a large number of tracer particles using the mean winds interpolated from the analysis fields plus random motions representing turbulence (Stohl and Thomson, 1999). Results presented here use European Centre for Medium-Range Weather Forecasts (ECMWF) analysis fields at 1 resolution. Here, the backward method is used to analyze transport pathways from source regions to the receptor position. Each simulation consists of 10 000 particles released every hour. Released particles were followed up to 14 d backward in time. Potential emission sensitivity is considered when particles reside below the boundary layer height as retrieved from the ECMWF analysis. Potential emission sensitivities are then convolved with gridded surface fluxes from emission inventories to simulate mixing ratios at the ship position during the cruise.

2.5 Emission inventories

For simplicity, our simulation uses surface emissions of CO2 and CH4 that include only anthropogenic emissions. Data are taken from the Emission Database for Global Atmospheric Research (EDGAR) v4.3.2 for the year 2012 (Janssens-Maenhout et al., 2017, 2019). EDGAR estimates country total emissions for a variety of species including CO2 and CH4 based on international emission factors and activity data. Maps are eventually generated at a resolution of 0.1× 0.1 according to spatial proxies (see Fig. A1 for CH4 and Fig. A2 for CO2 in Appendix A). For example, for fossil fuels exploitation, spatial proxies combine observed flaring and cartography of infrastructures. Our FLEXPART simulations calculate contributions detailed by sector and by country following the typology of the EDGAR database.

The Evaluating the Climate and Air Quality Impacts of Short-Lived Pollutants (ECLIPSE) baseline version CLE 5a (Höglund-Isaksson, 2012) is also used for comparative simulations of CH4 using FLEXPART (CO2 is not provided as gridded data). ECLIPSE is derived from the GAINS model used for scenario analysis. ECLIPSE maps emissions at a resolution of 0.5× 0.5 according to spatial proxies that are different from EDGAR.

3 Results and discussion

3.1 Spatial distribution of CH4 and CO2 along the cruise

Figure 2 shows the concentrations observed during the two legs of the cruise. GPS data acquisition was missing in part of the western Mediterranean, and therefore parts of the data are not shown on the map. Figure 3a shows the measured CH4 mixing ratios statistics for each region during the cruise, excluding data influenced by the ship's own exhaust and data acquired when the ship is at rest. Median values in the western and eastern Mediterranean were 1903 ppb (interquartile range, IQR, 1896–1907 ppb) and 1906 ppb (IQR 1879–1929 ppb), respectively, but significant differences were observed for each region between the two legs. For the eastern Mediterranean, Leg 1 median CH4 concentration is 1878 ppb under a general northern wind influence, whereas Leg 2 median is 1929 ppb with a generally southern wind. Closer to the Suez Canal, wind comes more generally from the south. The Gulf of Aden, the Arabian Sea and the Gulf of Oman were markedly lower in CH4, suggesting a different air mass origin. This is confirmed by a back-trajectory analysis using the FLEXPART model, showing a markedly different regional pattern in sensitivity in particle residence time in the boundary layer over Ethiopia, Somalia and the surrounding ocean (Fig. A3a). This is also consistent with the very low total OH reactivity observed in the area during the campaign (Pfannerstill et al., 2019). Wind direction in this area is mostly parallel to the ship track.

Figure 2Measurements of CH4 (a, c) and CO2 (b, d) collected during Leg 1 from Toulon to Kuwait between 24 June and 31 July 2017 (a, b) and Leg 2 from Kuwait to Toulon between 3 August and 2 September 2017 (c, d). Background map: © OpenStreetMap contributors 2021. Distributed under the Open Data Commons Open Database License (ODbL) v1.0.

Figure 3b shows the CO2 mixing ratio statistics for each region and each leg observed during the campaign. For each area, CO2 concentrations were higher during Leg 1 compared to Leg 2, likely reflecting hemispheric-scale drawdown in CO2 associated with uptake by Northern Hemisphere vegetation (see, e.g., Ramonet et al., 2010). The highest median concentrations are observed in the Arabian Gulf (410.27 ppm, IQR 407.16–415.36 ppm), whereas the lowest concentrations were encountered in August over the eastern Mediterranean, with 403.30 ppm (IQR 399.12–407.64 ppm), which was the closest value to the Mace Head mean mixing ratio in July and August. Ramonet et al. (2010) proposed that the minor changes in concentrations differences observed over eastern Europe compared to the marine baseline are due to biospheric uptake in Europe being offset by anthropogenic (fossil fuel) emissions along air mass pathways. Here, air mass origin is over eastern Europe's boundary layer over the previous 10 d according to FLEXPART back-trajectories (Fig. A3b).

Figure 3Methane (a) and carbon dioxide (b) mixing ratios observed during the cruise. See Fig. 1 for region definition. Data from when the ship was at rest or contaminated by the ship's own exhaust are excluded. The box plot shows the interquartile range (25 %–75 %) and median values for each grouping of data. The bars extend to the lowest and highest values, (ignoring outliers beyond 1.5 times the interquartile range). The statistics of the two legs are represented separately. The horizontal black line is the mean methane mixing ratio at Mace Head, Ireland (NOAA), during July and August 2017.


The median mixing ratio of CO2 over the Aden and Arabian Sea are comparable to Red Sea and Mediterranean values, contrary to CH4 where mixing ratios were distinctly lower. CO2 mixing ratios do, however, exhibit a much lower variability (interquartile range) compared to other regions.

In the remainder of this paper, enhancements of CH4 and CO2 are defined as excess concentrations over a background that is defined for each region and each leg, according to Eq. (1).

(1) Δ X = [ X ] - [ X ] background

The background is defined for the whole campaign as the 10th percentile of measured mole fractions.

3.2 Significance of fugitive sources based on light alkane measurements

In this section we investigate the origin of CH4 enhancements in the regions where the highest median mixing ratio were observed. These regions, the Arabian Gulf, the Suez Canal and the Red Sea, are all areas where both offshore and onshore OG exploration are important (Fig. 4). Although the ship's track did not target specifically offshore platforms, the cruise collected a significant amount of strong CH4 enhancements coinciding with crossing of offshore platforms. Bourtsoukidis et al. (2019) showed that emissions from OG were significant contributions to non-methane hydrocarbons measured during the cruise over several areas around the Arabian Peninsula. According to the EDGAR emission inventory (Janssen-Maenhout et al., 2017), the OG industry is responsible for up to 80 % of CH4 emissions in the Middle East. Local enhancements of mixing ratio in air masses are typically linked to local or regional sources. The most likely origin of the high CH4 and CO2 concentrations over the Arabian Gulf and the Suez Canal is anticipated to be anthropogenic emissions. In this section we investigate the origin of CH4 enhancements by assessing their correlation with alkanes measured onboard the ship. In this analysis, all CH4 and CO2 measurements are averaged specifically over the sampling interval of the GC-FID.

Figure 4Ship track superimposed with flaring sources detected by NASA/NOAA Visible Infrared Imaging Radiometer Suite active in 2017. Dot size is proportional to the highest observed radiative heat (in MW), which in turn is related to the volume of gas flared (Elvidge et al., 2016). Flaring occurs mostly at the upstream production sites and is therefore a proxy for the presence of extraction and production sites. Flaring spots retrieved from (last access: 16 April 2019). Background map: © OpenStreetMap contributors 2021. Distributed under the Open Data Commons Open Database License (ODbL) v1.0.

Over the whole campaign, we observed a significant but small correlation between CH4 and C2H6 (r=0.386, p< 10−15). The highest C2H6 concentrations are observed over the Arabian Gulf (Fig. 5). In this region, a significant correlation between C2H6 and CH4 (r=0.332, p=0.004) was derived, with a region-wide regression slope of 0.062 ppb ppb−1. Similarly, CH4 was correlated with propane (r= 0.356, p= 0.002) and the sum of i-butane and n-butane (r= 0.360, p= 0.002). CH4 was also correlated with CO2 (r= 0.53, p= 2 × 10−12) and CO (r= 0.76, p< 10−15) in the Arabian gulf, suggesting that CH4 variability is dominated by emissions from OG activities colocated with combustion sources. Flaring may typically contribute to such an emission profile. However, CO2 showed no significant correlation with any of the three alkanes mentioned above. This suggests that CO2 variability in the area is dominated by other factors than OG sources.

Figure 5Scatterplot of ethane and methane for the Gulf of Oman, Arabian Gulf, Red Sea and Suez Canal. The regression line and associated 95 % confidence interval are plotted for the Arabian Gulf, Red Sea and Suez Canal regions.


Over the Red Sea, CH4 enhancements of about 200 ppb above background have been observed. A clear correlation between C2H6 and CH4 (r= 0.71, p< 10−15) is associated with a regression slope of 0.047 ppb ppb−1. No significant correlation is found for CH4 against CO2 (r=0.08, p= 0.1368) or CO (r= 0.23, p= 2 × 10−15).

For the Suez Canal area, the regression slope between C2H6 and CH4 is 0.006 ppb ppb−1(r= 0.75, p< 10−15). Strong correlations are found between CH4 and CO2 (r= 0.83, p< 10−15) and CH4 and CO (r= 0.87, p< 10−15). This indicates that combustion sources are strongly contributing to CH4 variability in the Suez Canal. The air mass originates from the eastern Mediterranean and Egypt according to FLEXPART back-trajectory simulations (Fig. A3c). Differences in C2H6 to CH4 regression slopes between the Suez Canal and the Red Sea may therefore reflect either an aggregation of different OG emission sources or different fraction of OG in observed CH4 enhancements.

Ethane is a valuable proxy for fossil fuel emissions at the global scale and has provided evidence of the reduction of fossil fuel emission in the 2000s (Simpson et al., 2012). As a proxy at the site scale, it contributes to the separation of thermogenic and biogenic emissions (e.g., Rella et al., 2015; Assan et al., 2017; Lowry et al., 2020, Defratyka et al., 2021). However, the range of C2H6 : CH4 ratio can be highly variable at regional scale, from emission site to site, and along the production chain. Measuring downwind of > 100 offshore platforms in the Gulf of Mexico, Yacovitch et al. (2020) found C2H6 : CH4 ratio ranging from 0.16 % to 17 %. The gas associated with the Deepwater Horizon oil spill sampled in the hydrocarbon plume rising from the seafloor showed a C2H6 : CH4 ratio between 8.1 % and 8.3 % (Reddy et al., 2012). This ratio depends on the temperature of formation of natural gas (Whitticar, 1994). The ratio can be further modified if a significant fraction of CH4 is emitted through flaring, with different depletion during combustion (Yacovitch et al., 2020). Hence, due to its inherent variability from field to field, the C2H6 : CH4 ratio alone is insufficient without additional constraints to quantify the relative contribution of OG to CH4 emissions at the regional scale.

Figure 6 shows the i-pentane to n-pentane ratio for selected regions. This ratio is another tracer of the contribution of natural gas, vehicle exhaust or fuel evaporation emissions to the atmospheric burden of alkanes and associated species (Gilman et al., 2013; Bourtsoukidis et al., 2019). For natural gas the ratio is expected to be 0.86 (indicated as a horizontal line in Fig. 6). Since the ratio is conserved against oxidation by OH in the atmosphere, the tracer is conserved during atmospheric transport. Isomerization of oil after distillation is widely used to obtain gasoline with higher octane indices. This process results in a decreased nC5 concentration in gasoline and hence a higher iC5 / nC5 ratio in urban environments. In natural gas processing, pentane is part of the condensate separated from the commercial natural gas. In the Arabian Gulf, where the mean CH4 enhancement was 67 ppb and the maximum was 255 ppb, the linear regression slope was found to be 0.939 ± 0.023 ppb ppb−1. The mean ratio of the isomers was 1.00 ± 0.12, suggesting a clear dominance of sources linked to OG extraction and production.

Figure 6Scatterplot of i-pentane /n-pentane and methane (a) and i-pentane/n-pentane and CO2 (b) for the Oman Gulf, Arabian Gulf, Red Sea and Suez Canal. The value of the i-pentane to n-pentane ratio associated with oil and gas emission (0.86) is drawn as a horizontal line in the two panels.


3.3 Photochemical age and remoteness of sources

Alkanes have different atmospheric lifetimes, and (if species are co-emitted at the same time and with a known emission ratio) measured changes in their ratios can help trace the photochemical age of an air mass (Purvis et al., 2003), and hence they can provide an indication of a transport time from the source to the receptor. The ratio of species is less sensitive to dilution and air mass mixing than the concentration of individual species. Transport time and photochemical age correlate particularly well for lighter alkanes (Parrish et al., 2007). For the AQABA campaign, Pfannerstill et al. (2019) used decreases in the [C3H6O] / [C3H8] ratio and simultaneous increases in the toluene / benzene ratio in the Arabian Gulf as an indicator of nearby sources. Wang et al. (2020) analyzed the mismatch between the toluene / benzene measured ratio and emission ratio. Here we compare, for several regions, (1) the decorrelation that occurs between CH4 and light alkanes (propane and pentane) during early photochemical aging due to differences in kinetic rates against OH and (2) the photochemical age of air masses exhibiting significant correlations between propane and pentane.

Propane (C3H8) has a lifetime of ca. 5.5 d, while pentane (C5H12) has a lifetime of 0.5–1 d with respect to OH radical oxidation (Bourtsoukidis et al., 2019). Figure 7 shows scatterplots of CH4 to propane and pentane (the latter as the sum of n-pentane and i-pentane). Correlation between CH4 and propane is relatively homogeneous across regions, albeit with very different slopes. For the Arabian Gulf the correlation coefficient is r= 0.37, for the Gulf of Oman it is r= 0.46, and for the Suez Canal we obtain r= 0.50. The correlation of CH4 to pentane is more contrasted. For the Arabian Gulf and Suez Canal the correlation coefficient remains significantly high (r= 0.40 and r= 0.67, respectively), which is comparable to the propane to CH4 correlations. The significant correlations suggest that the air mass is affected by emissions with constant ratios, and hence there is a regional signature of local sources. This corroborates the results of Wang et al. (2020), who found limited departure between the measured ratio and emission ratio for these two regions. In contrast, in the Gulf of Oman there is no correlation between CH4 and propane (r= 0.003). This suggests that high CH4 concentrations observed in the Arabian Gulf and Suez Canal are mainly determined by relatively nearby sources co-emitting CH4, pentane and propane. This approach suggests that excess CH4 observed over the Gulf of Oman is emitted from remote sources or without co-emission of these species.

Figure 7Scatterplot of (a) pentane and (b) propane against methane for the Gulf of Oman, Arabian Gulf, Red Sea and Suez Canal. The regression line and associated uncertainty are plotted for each of these four regions.


As a next step we evaluate the average photochemical age of air masses for the Suez Canal and Arabian Gulf regions. For that purpose, we examine, at the regional scale, the change in measured propane to (n+i)-pentane ratio relative to the initial emission ratio. We follow the definition of photochemical age against OH defined in Parrish et al. (2007), given here in Eq. (2) below.

(2) t a [ OH ] = t = t E t R [ OH ] d t = - 1 k A - k B ln A B - ln A 0 B 0 ,

where ta is the average photochemical age against OH, 〈[OH]〉 is the mean concentration of OH along the pathway, tE refers to time of emission and tR time of reception, A and B are the two species propane and pentane, and A0 and B0 are the initial ratio at the time of emission. We evaluate the mean ratio of concentration between species at the regional scale as the regression slope between the two species. The regression slope of pentane against propane is 0.149 mol mol−1 for the Suez Canal and 0.378 mol mol−1 for the Arabian Gulf. The ratio of species at the time of origin is approximated by the emission ratio from the EDGAR v4.3.2 inventory for volatile organic compounds (Huang et al., 2017). The inventory-based emission ratio is 0.460 mol mol−1 for the total emissions of Middle East countries and 0.435 mol mol−1 for Egypt. Overall this value varies within a limited span (ranging from 0.347 mol mol−1 for Kuwait to 0.682 mol mol−1 for Iran) across countries with significant emissions of these two species. We use kinetic rate coefficients as reported by Pfannerstill et al. (2019, their Table S1) and assume a mean OH concentration of 7.5 × 106 mol cm−3 following Wang et al. (2020). This approach yields a photochemical age for our air masses of 1.5 d for the Suez Canal and 0.28 d for the Arabian Gulf. At a median true wind speed of 3.9 m s−1 (as observed over the Arabian Gulf and over the Suez Canal), this suggests that sources are distant on average by 206 and 38 km, respectively, for the Suez Canal and the Arabian Gulf.

However, this crude approach to the assessment of remoteness of sources of CH4 excesses in the Arabian Gulf and the Red Sea is not sufficient on its own and should be compared to tagged transport simulations.

3.4 Tagged tracer simulation of anthropogenic methane

We simulated the anthropogenic excess mixing ratio of CH4 and CO2 using FLEXPART and the EDGAR emission inventory. Figure 8 shows a comparison of simulated excess CH4 time series at the receptor position and ΔCH4 measured from the ship. The model reproduces the region-wise variability of the signal reasonably well. In the distinctly low CH4 concentration area over the Gulf of Aden and Arabian Sea, the model appropriately simulates enhancements close to zero. In the high-concentration area of the Arabian Gulf (from 28 July to 6 August) the model captures the presence of strong enhancements, despite being poorly positioned along the ship track. Over the Red Sea, the model captures the variability correctly. The regional baseline of the signal is, however, challenging to reproduce with the model since measurements are expressed as excess CH4 with a background calculated by leg and across all regions. This sharply highlights the difference between the low-concentration areas (Gulf of Aden and Arabian Sea) versus the other areas more influenced by anthropogenic emissions.

Figure 8Comparison of simulated excess CH4 time series at the receptor position and hourly average ΔCH4 measured from the ship (black line). All data at harbors are included. The simulation combines FLEXPART potential emission sensitivity and EDGAR inventory (stacked areas). The main source sectors are indicated using a color code. The potential emission sensitivity is also combined with the ECLIPSE inventory for comparison (blue line). The stops at each harbor are indicated as shaded gray areas, with the names of the harbors indicated at the top.


The OG emissions contribute to 35 % of the simulated excess CH4 overall. Over the Arabian Gulf and Gulf of Oman, it accounts, respectively, for 81 % and 63 % of CH4 excess on average, consistent with the regional emission pattern and alkane measurements presented above. Over the Red Sea, OG emission of CH4 explains 37 % of the simulated variability on average. Waste management represents on average 27 % of the simulated CH4 excess for the whole campaign. In the Arabian Gulf, waste management represents only 12 % of CH4 variability.

The observations at harbor locations (gray shading in Fig. 8) exhibit a strong diurnal pattern that is well reproduced by the model, probably due to the proximity of the land mass and the diurnal cycle of the atmospheric boundary layer development. However, the amplitude of the diurnal cycle is rather underestimated by the model, notably for the stop-over in Jeddah between 10 and 13 July. At harbors, the OG signal remains dominant, especially in Kuwait. The simulated share of waste management sources to the total CH4 enhancement is twice as much as that found over open sea. However, due to the proximity of sources within the inventory spatial resolution, the model may not fully resolve the ratio between these contributions. Closer to Kuwait, back-trajectories show a stronger sensitivity to potential sources in Iran.

The highest CH4 concentrations over the Arabian Gulf are mostly observed in its northern part, while the simulation predicts higher concentrations in the southern part, especially in the vicinity of the South Pars (North Dome) gas field. However, the CH4 enhancement observed in the immediate vicinity of this gas field is much lower than the simulated one. The model points to an “OG exploitation” source type in the northwestern half of the Arabian Gulf for the strongest CH4 enhancements observed. In a wide US survey of 8000 facilities, Lyon et al. (2016) found that super-emitter sites were 3 times more frequently associated to oil than the gas assets and that 90 % of sources were from tank vents and hatches used in oil storage. Our study would suggest an underestimation of oil upstream emissions by the inventories and an overestimation of leak rates in the upstream natural gas industry.

Figure 9Comparison of simulated excess CO2 time series at the receptor position in the Gulf of Oman and Arabian Gulf and hourly average ΔCO2 measured from the ship (black line). Flagged, hourly average data are shown in light red. The simulation combines FLEXPART potential emission sensitivity and the EDGAR inventory (stacked areas). The main source sectors are indicated using a color code. The stops at each harbor are indicated as shaded gray areas, with the names of the harbors indicated at the top.


Figure 9 shows the measured excess CO2 while navigating in the Gulf of Oman and the Arabian Gulf. The simulated excess is dominated by the power industry sector, accounting for most of the variability, with contributions from the oil refineries and oil transformation sectors and local maxima off the coast of United Arab Emirates. This area hosts the Jebel Ali oil refinery, which is dedicated to the production of liquefied petroleum gas, naphtha and a variety of fuel types, as well as the Takreer Abu Dhabi oil refinery, which is located further south (retrieved from, last access: 30 June 2021).

Figure 10Simulation of ΔCH4 by country of origin over the Arabian Gulf. The solid line shows the measured DCH4. The simulated excess CH4 is tagged and color coded by country or area of origin. The regions are indicated by the following abbreviations: Iraq_N is northern Iraq, Iraq_S is southern Iraq, PG_N is the northern half of the Persian (Arabian) Gulf, PG_S is the southern half of the Persian (Arabian) Gulf, RoW is the rest of the world, SA is Saudi Arabia, and UAE is the United Arab Emirates. The inset shows the trajectory of the ship corresponding to the simulation, and the labels A, B, C and D highlight points of the ship's route discussed in the main text. Red points correspond to flaring detected from space (see Fig. 4) at the time of the campaign. Inset background map: © OpenStreetMap contributors 2021. Distributed under the Open Data Commons Open Database License (ODbL) v1.0.

Figure 10 shows the contribution of countries or areas of origin to simulated CH4 in the Arabian Gulf during Leg 1, where emissions are dominated by fossil fuel exploitation (see discussion on alkanes above and simulation in Fig. 8). By tagging the simulated CH4 enhancements by source region, it is possible to link the contributions of source regions to the ability of the model to reproduce the observed signal.

Off the coast of Dubai and the United Arab Emirates (transect A–B in Fig. 10), the simulation is in good agreement with the observations and points to a strong contribution from the nearby state (United Arab Emirates) and the local offshore activities in the Arabian Gulf. The highest CH4 concentration of the section is measured (2109 ppb) within a thin plume lasting 15 min (labeled B in Fig. 10, 18:05–18:20 UTC). The model successfully captures this enhancement with an average of 150 ppb excess CH4, although in this case it is spread over 2 h. It is simulated as coming mostly (80 %) from local offshore emissions (Arabian Gulf South in Fig. 10) while the ship is moving through the Fateh oil field off Dubai. Approaching the North Dome field off the coast of Qatar (point C, and during 29 July), however, CH4 is strongly overestimated by the model, with simulated enhancements of more than 400 ppb having no equivalent pattern in the measurement. In the model this strong enhancement is due to local offshore emissions (southern part of the Arabian Gulf) associated with the South Pars (North Dome) gas field. South Pars (North Dome) field is the largest gas field in the world and is shared between Iran and Qatar (Conti et al., 2016). In the northern part of the Arabian Gulf (transect from C to D), significant CH4 enhancements are measured, and the model fails to reproduce these enhancements, although it suggests local emissions (northern Arabian Gulf and Kuwait) are dominant, consistent with the photochemical age of air masses (Sect. 3.3).

3.5 Model–data comparison

Can we confirm or verify the inventories based on the measurements? The entire dataset (excluding data flagged as contaminated and stationary measurements at harbors) has a correlation of r= 0.13. The agreement improves during nighttime with r= 0.22. Overall, the simulation represents the variability of CH4 well at the synoptic scale. The first-order discrepancies arise from the difference in background between regions, since here the background used to offset the CH4 measurements is calculated for the whole campaign. While the measured ΔCH4 is an excess over a pre-defined background, the simulated excess CH4 integrates anthropogenic surface sources during 14 d prior to reaching the receptor. Therefore, in order to compare these two quantities, the mismatch between the boundary conditions of the model (i.e., CH4 of the air mass prior to the 14 d before measurement) and the background defined for the measurement data introduces a varying offset. However, the simulated time series can be considered an indication of the geographical or sectoral origin of the observed excess mixing ratio.

Figure 11 provides a comparison of the fractional contribution of OG exploitation sources to simulated CH4 excess with the i-pentane to n-pentane ratio (see Sect. 3.2). We find that an exponential decrease model can fit the data with a significant correlation r=−0.38 (p< 10−16). This exponential model yields an estimated i-pentane to n-pentane ratio of 0.99 when the simulation fraction of OG is 100 %, whereas an absence of OG exploitation would yield a ratio of 1.84. This confirms that the fraction of emissions linked to OG exploitation in the model is consistent with our expectations according to the isomeric ratio of pentane (Bourtsoukidis et al., 2019).

Figure 11Scatterplot of the fraction of simulated methane due to OG sources against the i-pentane to n-pentane ratio. An exponential model is fitted to the data (red line).


The inventory used for the simulation also integrates a large part of uncertainty in terms of intensity and location of the sources. To test the sensitivity to the emission inventory itself, we simulated CH4 excess using the ECLIPSE inventory for CH4 (light blue line in Fig. 8). The two inventories agree to a large extent for all regions apart from the Arabian Gulf, where ECLIPSE yields a mean excess that is 65 % higher than the EDGAR simulation (mean 255 and 156 ppb, respectively). For the Arabian Gulf, none of the inventories yield significant correlation, suggesting a poor representation of spatial distribution of emissions. ECLIPSE also simulates much higher CH4 excess in the vicinity of the Fateh oil field, with a local maximum of 566 ppb compared to the EDGAR simulation that yields 150 ppb. In the northern part of the Arabian Gulf, ECLIPSE also simulates higher excess concentrations than EDGAR, but the ECLIPSE simulation shows a much better agreement with the observed CH4 variability, especially at the Kuwait coast on 30 and 31 July. The enhanced divergence among simulations using different inventories is most likely due to the increasing differences between inventories at finer spatial scales, as has been highlighted by Ciais et al. (2010) for CO2 in Europe. Furthermore, the high skewness in the regional distribution of emission for each facility may not be properly represented in the inventories (Zavala-Araiza et al., 2015), leading to strong discrepancies between simulation and measurements when measuring directly near such facilities. Overall, in the southern part of the Arabian Gulf, where gas fields represent a larger part of the OG extraction, EDGAR tends to accurately represent specific plumes (see point B in Fig. 10), which is not the case for ECLIPSE. However, on average both EDGAR and ECLIPSE tend to overestimate the signal (Fig. 8). In the northern part of the Arabian Gulf, where oil extraction is comparatively more represented, the EDGAR inventory tends to underestimate the measurements, while ECLIPSE underestimates the measurements for Leg 1 and overestimates it for Leg 2.

Accurately simulating plumes measured in the vicinity of point sources is highly dependent on the injection position of emissions in the inventory. Here we used inventory maps in the model with a resolution of 0.1× 0.1 for EDGAR; hence, a point source is spread within 10 km. This may result in discrepancies between observations and simulated fields. This is particularly true if an important source point is within the distance corresponding to the spatial resolution of the wind fields driving the model (here, 1× 1, i.e., approximately 100 km). Capturing local sea breeze patterns might be challenging for the model's driving wind field. This may cause the strong overestimations of CH4 enhancements over the southeastern part of the Arabian Gulf.

Unaccounted for time variations may also play a role in model–observation mismatch. Since the inventory is static, any daily or weekly pattern in emissions would not be reproduced in the model, which could be important when the ship is in harbor. Moreover, venting, incomplete flaring combustion and maintenance activities are typically leading to intermittent CH4 emission. In a study of black carbon emissions from flaring in Siberia, Petäjä et al. (2020, their Sect. 3.8.2) found that accounting for actual times of flaring instead of using annual means can significantly improve the simulation of downwind atmospheric measurements.

Since our study is focusing on anthropogenic emissions and given the low natural emissions in the eastern Mediterranean and Middle East area (Saunois et al., 2020), we used only anthropogenic inventories in our study. However, some level of sensitivity to natural fluxes can be expected. The CO2 sink would affect concentrations downwind of large forested areas, for example around the eastern Mediterranean basin. However, biogenic sinks are not likely to play a role in the desert-dominated regions of the Arabian Peninsula.

OG extraction at the country level varies from year to year, and the life cycle of individual productive fields evolves over several years. As a result, changes over time of the emission spatial pattern and intensity as reported in inventories for specific years may affect the accuracy of the simulation. Here, our simulation uses the EDGAR inventory for 2012. Oil production in the Middle East increased by 10.8 % between 2012 and 2017, while gas production increased by 19.2 % (Dudley, 2019). According to the most recent available update of the inventory (EDGAR v5.0), OG emissions have increased by 8.61 % for countries bordering the Arabian Gulf between 2012 and 2015, the latest available year. In order to investigate the potential evolution of the oil extraction activity, the evolution of flaring activity between 2012 and 2017 has been investigated. Flaring data are obtained from the VIIRS Skytruth nightlight product (Elvidge et al., 2016) for the Middle East. Flaring has shown a 2-fold increase in intensity (in terms of number of hotspots detected from space) over the period, which suggests that the 8.61 % increase in inventory OG emissions from the neighboring countries is not sufficient to match the increase in extraction activity. The spatial distribution of the flaring has not varied significantly over the period. In addition to flaws in the inventory for its base year, accounting for increased activity over the period 2012-2017 would lead to further overestimations between model and observations linked to OG activity in the Arabian Gulf. This therefore does not contribute to explaining the model overestimation, especially in the southern part of the Arabian Gulf.

4 Conclusion

The AQABA campaign provided the first overview of the regional distribution of CH4 and CO2. Three distinct elements are identified in the distribution of GHG during the campaign. The Mediterranean Sea is dominated by European emissions. The Gulf of Aden and the Arabian Sea are in an air mass relatively poor in GHG, with air masses originating from eastern Africa. The Red Sea, the Gulf of Oman and Arabian Gulf showed high CO2 and CH4 concentrations.

The C2H6 : CH4 and i-pentane to n-pentane ratios suggest that over the Red Sea, Gulf of Oman and Arabian Gulf, CH4 enhancements are originating from OG emissions, especially in the latter area. Repeated CH4 enhancements over the Suez Canal, Red Sea and Arabian Gulf are also unambiguously identified as emitted from local OG extraction and exploitation. This is clearly supported by a Lagrangian simulation based on the EDGAR inventory, showing that 81 % of excess CH4 over the Arabian Gulf is due to the exploitation of OG. The CO2 variability is dominated in this area by anthropogenic emissions, with a dominant contribution from the power industry and the oil-refining and transformation sectors.

While the model predicts a dominant contribution from OG exploitation to CH4 enhancements, in agreement with the light alkanes measurements, only weak quantitative agreement has been found between modeled and measured CH4. This is most likely explained by a combination of factors, including error in the inventory, poor dilution in the model in the vicinity of the sources, poor distribution of the point source intensities and a lack of representation of temporal variations in emission patterns. A similar simulation using a different inventory (ECLIPSE) tended to overestimate the measured CH4 close to sources in the southern Arabian Gulf. Despite these compounded uncertainties, our study provides strong indication that inventories overestimate part of the regional upstream OG emissions in Middle Eastern countries neighboring the Arabian Gulf, especially linked to gas extraction. In contrast, in the northern Arabian Gulf there are comparatively more oil fields, and in that region the measured methane is generally underestimated by the simulations. Increases in OG emissions in the Middle East compared to the reference year of the inventories would further enhance overestimations but could partly explain the underestimation for the northern Arabian Gulf area. Our study also shows that the inventories must be improved for the spatial distribution of emissions in this area, where emissions are dominated by the sector of OG extraction, transport and transformation. Separating onshore from offshore emissions would require a dedicated modeling study investigating local atmospheric circulation combined with atmospheric composition measurements.

More targeted measurements specifically investigating offshore and onshore extraction sites will enable a better understanding of the distribution of emission in the area. Assessing emissions from individual wells and processing facilities with dedicated measurements and combining these estimates at the regional level would be necessary to further improve our knowledge of actual OG emissions in the Middle East. Combining CH4 and alkane measurements offers the potential to (1) separate the OG component of measured CH4 from other sources and (2) to investigate the fraction of emissions linked to each phase of OG value chain at the regional and local level, from extraction to end usage through storage, transport and processing.

Appendix A

Figure A1CH4 emissions predicted by EDGAR v4.3.2 for year 2012 in the Mediterranean and Middle East, expressed in log(kg CH4 m−2 s−1). The ship trajectory is shown in blue.

Figure A2CO2 emissions predicted by EDGAR v4.3.2 for year 2012 in the Mediterranean and Middle East, expressed in log 10(kg CO2 m−2 s−1). The ship trajectory is shown in blue.

Figure A3Footprints for ship positions on (a) 8 August 2017 at 06:00 UTC, (b) 27 August 2017 at 08:00 UTC and (c) 4 July 2017 at 07:00 UTC. See the main text for details. Footprints are expressed as log 10 of residence time (s) of particles below 2000 m. The ship's position is indicated by the gray star.

Data availability

The data used in this work are available upon reasonable request directly from the authors.

Author contributions

JDP conceived this study and wrote the manuscript. AR, AB and JDP performed and analyzed the FLEXPART simulations. EB, JW and LE did the alkane measurements. JDP and MD did the greenhouse gas measurements. IT provided complementary data. HH was the chief of mission. JL was the PI of the campaign. All authors contributed to the manuscript.

Competing interests

The authors declare that they have no conflict of interest.


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


The Max Planck Institute for Chemistry, Mainz, organized the AQABA campaign and chartered the ship. Thanks to Marcel Dorf, Dieter Scharffe, Charlotte Beall, Claus Koeppel, and the many other participants of the cruise for their crucial contributions. Jean Sciare enabled the presence of our measurements during the campaign. Olivier Laurent, Céline Lett and Laurence Vialettes (ICOS teams at LSCE) were key contributors in designing and implementing the measurement setup. Thanks to Hays Ships and the R/V Kommandor Iona crew for their successful and friendly implementation of the cruise.

Financial support

This research has been supported by the EMME-CARE project from the European Union's Horizon 2020 Research and Innovation Programme (grant agreement no. 856612), and the Government of the Republic of Cyprus.

Review statement

This paper was edited by Eduardo Landulfo and reviewed by Carlos Oliveira and two anonymous referees.


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. 

Assan, S., Baudic, A., Guemri, A., Ciais, P., Gros, V., and Vogel, F. R.: Characterization of interferences to in situ observations of δ13CH4 and C2H6 when using a cavity ring-down spectrometer at industrial sites, Atmos. Meas. Tech., 10, 2077–2091,, 2017. 

Berchet, A., Pison, I., Crill, P. M., Thornton, B., Bousquet, P., Thonat, T., Hocking, T., Thanwerdas, J., Paris, J.-D., and Saunois, M.: Using ship-borne observations of methane isotopic ratio in the Arctic Ocean to understand methane sources in the Arctic, Atmos. Chem. Phys., 20, 3987–3998,, 2020. 

Bourtsoukidis, E., Ernle, L., Crowley, J. N., Lelieveld, J., Paris, J.-D., Pozzer, A., Walter, D., and Williams, J.: Non-methane hydrocarbon (C2–C8) sources and sinks around the Arabian Peninsula, Atmos. Chem. Phys., 19, 7209–7232,, 2019. 

Bourtsoukidis, E., Pozzer, A., Sattler, T., Matthaios, V. N., Ernle, L., Edtbauer, A., Fischer, H., Könemann, T., Osipov, S., Paris, J.-D., Pfannerstill, E. Y., Stönner, C., Tadic, I., Walter, D., Wang, N., Lelieveld, J., and Williams, J.: The Red Sea Deep Water is a potent source of atmospheric ethane and propane, Nat. Commun., 11, 447,, 2020. 

Celik, S., Drewnick, F., Fachinger, F., Brooks, J., Darbyshire, E., Coe, H., Paris, J.-D., Eger, P. G., Schuladen, J., Tadic, I., Friedrich, N., Dienhart, D., Hottmann, B., Fischer, H., Crowley, J. N., Harder, H., and Borrmann, S.: Influence of vessel characteristics and atmospheric processes on the gas and particle phase of ship emission plumes: in situ measurements in the Mediterranean Sea and around the Arabian Peninsula, Atmos. Chem. Phys., 20, 4713–4734,, 2020. 

Ciais, P., Paris, J. D., Marland, G., Peylin, P., Piao, S. L., Levin, I., Pregger, T., Scholz, Y., Friedrich, R., Rivier, L., Houwelling, S., Schulze, E. D., and members of the CARBOEUROPE synthesis team: The European carbon balance, Part 1: fossil fuel emissions, Global Change Biol., 16, 1395–1408, 2010. 

Conti, J., Holtberg, P., Diefenderfer, J., LaRose, A., Turnure, J. T., and Westfall, L.: International Energy Outlook 2016 With Projections to 2040, USDOE Energy Information Administration, United States: N. p.,, 2016. 

Defratyka, S. M., Paris, J.-D., Yver-Kwok, C., Loeb, D., France, J., Helmore, J., Yarrow, N., Gros, V., and Bousquet, P.: Ethane measurement by Picarro CRDS G2201-i in laboratory and field conditions: potential and limitations, Atmos. Meas. Tech., 14, 5049–5069,, 2021. 

Dlugokencky, E. and Tans, P.: Trends in atmospheric carbon dioxide, National Oceanic & Atmospheric Administration, Earth System Research Laboratory (NOAA/ESRL), available at:, last access: 26 July 2018. 

Dudley, B.: Statistical Review of World Energy 68th Edition, Technical Report, BP p.l.c., London, UK, 2019. 

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. 

Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Hauck, J., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Bakker, D. C. E., Canadell, J. G., Ciais, P., Jackson, R. B., Anthoni, P., Barbero, L., Bastos, A., Bastrikov, V., Becker, M., Bopp, L., Buitenhuis, E., Chandra, N., Chevallier, F., Chini, L. P., Currie, K. I., Feely, R. A., Gehlen, M., Gilfillan, D., Gkritzalis, T., Goll, D. S., Gruber, N., Gutekunst, S., Harris, I., Haverd, V., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Joetzjer, E., Kaplan, J. O., Kato, E., Klein Goldewijk, K., Korsbakken, J. I., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Marland, G., McGuire, P. C., Melton, J. R., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Neill, C., Omar, A. M., Ono, T., Peregon, A., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Séférian, R., Schwinger, J., Smith, N., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Werf, G. R., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2019, Earth Syst. Sci. Data, 11, 1783–1838,, 2019. 

Friedlingstein, P., O'Sullivan, M., Jones, M. W., Andrew, R. M., Hauck, J., Olsen, A., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Canadell, J. G., Ciais, P., Jackson, R. B., Alin, S., Aragão, L. E. O. C., Arneth, A., Arora, V., Bates, N. R., Becker, M., Benoit-Cattin, A., Bittig, H. C., Bopp, L., Bultan, S., Chandra, N., Chevallier, F., Chini, L. P., Evans, W., Florentie, L., Forster, P. M., Gasser, T., Gehlen, M., Gilfillan, D., Gkritzalis, T., Gregor, L., Gruber, N., Harris, I., Hartung, K., Haverd, V., Houghton, R. A., Ilyina, T., Jain, A. K., Joetzjer, E., Kadono, K., Kato, E., Kitidis, V., Korsbakken, J. I., Landschützer, P., Lefèvre, N., Lenton, A., Lienert, S., Liu, Z., Lombardozzi, D., Marland, G., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Niwa, Y., O'Brien, K., Ono, T., Palmer, P. I., Pierrot, D., Poulter, B., Resplandy, L., Robertson, E., Rödenbeck, C., Schwinger, J., Séférian, R., Skjelvan, I., Smith, A. J. P., Sutton, A. J., Tanhua, T., Tans, P. P., Tian, H., Tilbrook, B., van der Werf, G., Vuichard, N., Walker, A. P., Wanninkhof, R., Watson, A. J., Willis, D., Wiltshire, A. J., Yuan, W., Yue, X., and Zaehle, S.: Global Carbon Budget 2020, Earth Syst. Sci. Data, 12, 3269–3340,, 2020. 

GIE-MARCOGAZ: Potential ways the gas industry can contribute to the reduction of methane emissions, Report for the Madrid Forum, 5–6 June 2019, available at: (last access: 14 August 2021), 2019. 

Gilman, J. B., Lerner, B. M., Kuster, W. C., and de Gouw, J. A.: Source Signature of Volatile Organic Compounds from Oil and Natural Gas Operations in Northeastern Colorado, Environ. Sci. Technol., 47, 1297–1305,, 2013. 

Hazan, L., Tarniewicz, J., Ramonet, M., Laurent, O., and Abbaris, A.: Automatic processing of atmospheric CO2 and CH4 mole fractions at the ICOS Atmosphere Thematic Centre, Atmos. Meas. Tech., 9, 4719–4736,, 2016. 

Helfter, C., Mullinger, N., Vieno, M., O'Doherty, S., Ramonet, M., Palmer, P. I., and Nemitz, E.: Country-scale greenhouse gas budgets using shipborne measurements: a case study for the UK and Ireland, Atmos. Chem. Phys., 19, 3043–3063,, 2019. 

Hmiel, B., Petrenko, V. V., Dyonisius, M. N., Buizert, C., Smith, A. M., Place, P. F., Harth, C., Beaudette, R., Hua, Q., Yang, B., Vimont, I., Michel, S. E., Severinghaus, J. P., Etheridge, D., Bromley, T., Schmitt, J., Faïn, X., Weiss, R. F., and Dlugokencky, E.: Preindustrial 14 CH4 indicates greater anthropogenic fossil CH4 emissions, Nature, 578, 409–412,, 2020. 

Höglund-Isaksson, L.: Global anthropogenic methane emissions 2005–2030: technical mitigation potentials and costs, Atmos. Chem. Phys., 12, 9079–9096,, 2012. 

Huang, G., Brook, R., Crippa, M., Janssens-Maenhout, G., Schieberle, C., Dore, C., Guizzardi, D., Muntean, M., Schaaf, E., and Friedrich, R.: Speciation of anthropogenic emissions of non-methane volatile organic compounds: a global gridded data set for 1970–2012, Atmos. Chem. Phys., 17, 7683–7701,, 2017. 

IEA: Methane Tracker 2020, IEA, Paris, France, available at: (last access: 14 August 2021), 2020. 

Jackson, R. B., Friedlingstein, P., Andrew, R. M., Canadell, J. G., Le Quéré, C., and Peters, G. P.: Persistent fossil fuel growth threatens the Paris Agreement and planetary health, Environ. Res. Lett., 14, 121001,​​​​​​​, 2019. 

Jackson, R. B., Saunois, M., Bousquet, P., Canadell, J. G., Poulter, B., Stavert, A. R., Bergamaschi, P., Niwa, Y., Segers, A., and Tsuruta, A.: Increasing anthropogenic methane emissions arise equally from agricultural and fossil fuel sources, Environ. Res. Lett., 15, 071002,, 2020. 

Janssens-Maenhout, G., Crippa, M., Guizzardi, D., Muntean, M., Schaaf, E., Olivier, J. G. J., Peters, J. A. H. W., and Schure, K. M.: Fossil CO2 and GHG emissions of all world countries, EUR 28766 EN, Publications Office of the European Union, Luxembourg, ISBN 978-92-79-73207-2,, 2017. 

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. 

Jones, A., Haywood, J. M., and Jones, C. D.: Can reducing black carbon and methane below RCP2.6 levels keep global warming below 1.5 C?, Atmos. Sci. Lett., 19, e821,, 2018. 

Lelieveld, J., Berresheim, H., Borrmann, S., Crutzen, P. J., Dentener, F. J., Fischer, H., Feichter, J., Flatau, P. J., Heland, J., Holzinger, R., Korrmann, R., Lawrence, M. G., Levin, Z., Markowicz, K. M., Mihalopoulos, N., Minikin, A., Ramanathan, V., de Reus, M., Roelofs, G. J., Scheeren, H. A., Sciare, J., Schlager, H., Schultz, M., Siegmund, P., Steil, B., Stephanou, E. G., Stier, P., Traub, M., Warneke, C., Williams, J., and Ziereis, H.: Global air pollution crossroads over the Mediterranean, Science, 298, 794–799,, 2002. 

Lelieveld, J., Hadjinocolaou, P., Kostopoulou, E., Chenoweth, J., El Maayar, M., Giannakopoulos, C., Hannides, C., Langa, M. A., Tanarhte, M., Tyrlis, E., and Xoplaki, E.: Climate change and impacts in the eastern Mediterranean and Middle East, Clim. Change, 114, 667–687, 2012. 

Le Quéré, C., Andrew, R. M., Friedlingstein, P., Sitch, S., Hauck, J., Pongratz, J., Pickers, P. A., Korsbakken, J. I., Peters, G. P., Canadell, J. G., Arneth, A., Arora, V. K., Barbero, L., Bastos, A., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Doney, S. C., Gkritzalis, T., Goll, D. S., Harris, I., Haverd, V., Hoffman, F. M., Hoppema, M., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Johannessen, T., Jones, C. D., Kato, E., Keeling, R. F., Goldewijk, K. K., Landschützer, P., Lefèvre, N., Lienert, S., Liu, Z., Lombardozzi, D., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S., Neill, C., Olsen, A., Ono, T., Patra, P., Peregon, A., Peters, W., Peylin, P., Pfeil, B., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rocher, M., Rödenbeck, C., Schuster, U., Schwinger, J., Séférian, R., Skjelvan, I., Steinhoff, T., Sutton, A., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Laan-Luijkx, I. T., van der Werf, G. R., Viovy, N., Walker, A. P., Wiltshire, A. J., Wright, R., Zaehle, S., and Zheng, B.: Global Carbon Budget 2018, Earth Syst. Sci. Data, 10, 2141–2194,, 2018. 

Lowry, D., Fisher, R. E., France, J. L., Coleman, M., Lanoisellé, M., Zazzeri, G., Nisbet, E. G., Shaw, J. T., Allen, G., Pitt, J., and Ward, R. S.: Environmental baseline monitoring for shale gas development in the UK: Identification and geochemical characterization of local source emissions of methane to atmosphere, Sci. Total Environ., 708, 134600,, 2020. 

Lyon, D. R., Alvarez, R. A., Zavala-Araiza, D., Brandt, A. R., Jackson, R. B., and Hamburg, S. P.: Aerial Surveys of Elevated Hydrocarbon Emissions from Oil and Gas Production Sites, Environ. Sci. Technol., 50, 4877–4886,, 2016. 

McNorton, J., Wilson, C., Gloor, M., Parker, R. J., Boesch, H., Feng, W., Hossaini, R., and Chipperfield, M. P.: Attribution of recent increases in atmospheric methane through 3-D inverse modelling, Atmos. Chem. Phys., 18, 18149–18168,, 2018. 

Nisbet, E. G., Fisher, R. E., Lowry, D., France, J. L., Allen, G., Bakkaloglu, S., Broderick, T. J., Cain, M., Coleman, M., Fernandez, J., Forster, G., Griffiths, P. T., Iverach, C. P., Kelly, B. F. J., Manning, M. R., Nisbet-Jones, P. B. R., Pyle, J. A., Townsend-Small, A., al-Shalaan, A., Warwick, N., and Zazzeri, G.: Methane mitigation: methods to reduce emissions, on the path to the Paris agreement, Rev. Geophys., 58, e2019RG000675,, 2020 

Ocko, I. B., Sun, T., Shindell, D., Oppenheimer, M., Hristov, A. N., Pacala, S. W., Mauzerall, D. L., Xu, Y., and Hamburg, S. P.: Acting Rapidly to Deploy Readily Available Methane Mitigation Measures by Sector Can Immediately Slow Global Warming, Environ. Res. Lett., 16, 054042,, 2021. 

OPEC: Annual statistical bulletin, 2020, available at:, last access: 9 December 2020. 

Parrish, D. D., Stohl, A., Forster, C., Atlas, E. L., Blake, D. R., Goldan, P. D., Kuster, W. C., and de Gouw, J. A.: Effects of mixing on evolution of hydrocarbon ratios in the troposphere, J. Geophys. Res., 112, D10S34,, 2007. 

Petäjä, T., Duplissy, E.-M., Tabakova, K., Schmale, J., Altstädter, B., Ancellet, G., Arshinov, M., Balin, Y., Baltensperger, U., Bange, J., Beamish, A., Belan, B., Berchet, A., Bossi, R., Cairns, W. R. L., Ebinghaus, R., El Haddad, I., Ferreira-Araujo, B., Franck, A., Huang, L., Hyvärinen, A., Humbert, A., Kalogridis, A.-C., Konstantinov, P., Lampert, A., MacLeod, M., Magand, O., Mahura, A., Marelle, L., Masloboev, V., Moisseev, D., Moschos, V., Neckel, N., Onishi, T., Osterwalder, S., Ovaska, A., Paasonen, P., Panchenko, M., Pankratov, F., Pernov, J. B., Platis, A., Popovicheva, O., Raut, J.-C., Riandet, A., Sachs, T., Salvatori, R., Salzano, R., Schröder, L., Schön, M., Shevchenko, V., Skov, H., Sonke, J. E., Spolaor, A., Stathopoulos, V. K., Strahlendorff, M., Thomas, J. L., Vitale, V., Vratolis, S., Barbante, C., Chabrillat, S., Dommergue, A., Eleftheriadis, K., Heilimo, J., Law, K. S., Massling, A., Noe, S. M., Paris, J.-D., Prévôt, A. S. H., Riipinen, I., Wehner, B., Xie, Z., and Lappalainen, H. K.: Overview: Integrative and Comprehensive Understanding on Polar Environments (iCUPE) – concept and initial results, Atmos. Chem. Phys., 20, 8551–8592,, 2020. 

Peters, G. P., Davis, S. J., and Andrew, R.: A synthesis of carbon in international trade, Biogeosciences, 9, 3247–3276,, 2012. 

Pfannerstill, E. Y., Wang, N., Edtbauer, A., Bourtsoukidis, E., Crowley, J. N., Dienhart, D., Eger, P. G., Ernle, L., Fischer, H., Hottmann, B., Paris, J.-D., Stönner, C., Tadic, I., Walter, D., Lelieveld, J., and Williams, J.: Shipborne measurements of total OH reactivity around the Arabian Peninsula and its role in ozone chemistry, Atmos. Chem. Phys., 19, 11501–11523,, 2019. 

Pisso, I., Sollum, E., Grythe, H., Kristiansen, N. I., Cassiani, M., Eckhardt, S., Arnold, D., Morton, D., Thompson, R. L., Groot Zwaaftink, C. D., Evangeliou, N., Sodemann, H., Haimberger, L., Henne, S., Brunner, D., Burkhart, J. F., Fouilloux, A., Brioude, J., Philipp, A., Seibert, P., and Stohl, A.: The Lagrangian particle dispersion model FLEXPART version 10.4, Geosci. Model Dev., 12, 4955–4997,, 2019. 

Purvis, R. M., Lewis, A. C., Carney, R. A., McQuaid, J. B., Arnold, S. R., Methven, J., Barjat, H., Dewey, K., Kent, J., Monks, P. S., Carpenter, L. J., Brough, N., Penkett, S. A., and Reeves, C. E.: Rapid uplift of nonmethane hydrocarbons in a cold front over central Europe, J. Geophys. Res.-Atmos., 108, 4224,, 2003. 

Ramonet, M., Ciais, P., Aalto, T., Aulagnier, C., Chevallier, F., Cipriano, D., Conway, T., Haszpra, L., Kazan, V., Meinhardt, F., Paris, J.-D., Schmid, M., Simmonds, P., Xueref-Rémy, I., and Necki, J.: A recent build-up of atmospheric CO2 over Europe. Part 1: observed signals and possible explanations, Tellus B, 62, 1–13,, 2010. 

Reddy, C. M., Arey, J. S., Seewald, J. S., Sylva, S. P., Lemkau, K. L., Nelson, R. K., Carmichael, C. A., McIntyre, C. P., Fenwick, J., Ventura, G. T., van Mooy, B. A. S., and Camilli, R. Composition and fate of gas and oil released to the water column during the Deepwater Horizon oil spill, P. Natl. Acad. Sci. USA, 109, 20,229–20,234,, 2012. 

Rella, C. W., Hoffnagle, J., He, Y., and Tajima, S.: Local- and regional-scale measurements of CH4, δ13CH4, and C2H6 in the Uintah Basin using a mobile stable isotope analyzer, Atmos. Meas. Tech., 8, 4539–4559,, 2015. 

Ricaud, P., Zbinden, R., Catoire, V., Brocchi, V., Dulac, F., Hamonou, E., Canonici, J., El Amraoui, L., Massart, S., Piguet, B., Dayan, U., Nabat, P., Sciare, J., Ramonet, M., Delmotte, M., di Sarra, A., Sferlazzo, D., Di Iorio, T., Piacentino, S., Cristofanelli, P., Mihalopoulos, N., Kouvarakis, G., Pikridas, M., Savvides, C., Mamouri, R.-E., Nisantzi, A., Hadjimitsis, D., Attié, J.-L., Ferré, H., Kangah, Y., Jaidan, N., Guth, J., Jacquet, P., Jaidan, J., Guth, P., Robert, C., Bourdon, A., Bourdinot, J.-F., Etienne, J.-C., Krysztofiak, G., Theron, P., and Chevrier, S.: The GLAM airborne campaign across the Mediterranean Basin, Bull. Am. Met. Soc., 99.2, 361–380,, 2018. 

Rogelj, J., Shindell, D., Jiang, K., Fifita, S., Forster, P., Ginzburg, V., Handa, C., Kheshgi, H., Kobayashi, S., Kriegler, E., Mundaca, L., Seferian, R., and Vilarino, M. V.: Mitigation Pathways Compatible with 1.5 C in the Context of Sustainable Development, in: Global Warming of 1.5 C. An IPCC Special Report on the impacts of global warming of 1.5 C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty, edited by: Masson-Delmotte, V., Zhai, P., Pörtner, H.-O., Roberts, D., Skea, J., Shukla, P. R., Pirani, A., Moufouma-Okia, W., Péan, C., Pidcock, R., Connors, S., Matthews, J. B. R., Chen, Y., Zhou, X., Gomis, M. I., Lonnoy, E., Maycock, T., Tignor, M., and Waterfield, T., available at: (last access: 14 August 2021), 2018. 

Saunois, M., Bousquet, P., Poulter, B., Peregon, A., Ciais, P., Canadell, J. G., Dlugokencky, E. J., Etiope, G., Bastviken, D., Houweling, S., Janssens-Maenhout, G., Tubiello, F. N., Castaldi, S., Jackson, R. B., Alexe, M., Arora, V. K., Beerling, D. J., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Crevoisier, C., Crill, P., Covey, K., Frankenberg, C., Gedney, N., Höglund-Isaksson, L., Ishizawa, M., Ito, A., Joos, F., Kim, H.-S., Kleinen, T., Krummel, P., Lamarque, J.-F., Langenfelds, R., Locatelli, R., Machida, T., Maksyutov, S., Melton, J. R., Morino, I., Naik, V., O'Doherty, S., Parmentier, F.-J. W., Patra, P. K., Peng, C., Peng, S., Peters, G. P., Pison, I., Prinn, R., Ramonet, M., Riley, W. J., Saito, M., Santini, M., Schroeder, R., Simpson, I. J., Spahni, R., Takizawa, A., Thornton, B. F., Tian, H., Tohjima, Y., Viovy, N., Voulgarakis, A., Weiss, R., Wilton, D. J., Wiltshire, A., Worthy, D., Wunch, D., Xu, X., Yoshida, Y., Zhang, B., Zhang, Z., and Zhu, Q.: Variability and quasi-decadal changes in the methane budget over the period 2000–2012, Atmos. Chem. Phys., 17, 11135–11161,, 2017. 

Saunois, M., Stavert, A. R., Poulter, B., Bousquet, P., Canadell, J. G., Jackson, R. B., Raymond, P. A., Dlugokencky, E. J., Houweling, S., Patra, P. K., Ciais, P., Arora, V. K., Bastviken, D., Bergamaschi, P., Blake, D. R., Brailsford, G., Bruhwiler, L., Carlson, K. M., Carrol, M., Castaldi, S., Chandra, N., Crevoisier, C., Crill, P. M., Covey, K., Curry, C. L., Etiope, G., Frankenberg, C., Gedney, N., Hegglin, M. I., Höglund-Isaksson, L., Hugelius, G., Ishizawa, M., Ito, A., Janssens-Maenhout, G., Jensen, K. M., Joos, F., Kleinen, T., Krummel, P. B., Langenfelds, R. L., Laruelle, G. G., Liu, L., Machida, T., Maksyutov, S., McDonald, K. C., McNorton, J., Miller, P. A., Melton, J. R., Morino, I., Müller, J., Murguia-Flores, F., Naik, V., Niwa, Y., Noce, S., O'Doherty, S., Parker, R. J., Peng, C., Peng, S., Peters, G. P., Prigent, C., Prinn, R., Ramonet, M., Regnier, P., Riley, W. J., Rosentreter, J. A., Segers, A., Simpson, I. J., Shi, H., Smith, S. J., Steele, L. P., Thornton, B. F., Tian, H., Tohjima, Y., Tubiello, F. N., Tsuruta, A., Viovy, N., Voulgarakis, A., Weber, T. S., van Weele, M., van der Werf, G. R., Weiss, R. F., Worthy, D., Wunch, D., Yin, Y., Yoshida, Y., Zhang, W., Zhang, Z., Zhao, Y., Zheng, B., Zhu, Q., Zhu, Q., and Zhuang, Q.: The Global Methane Budget 2000–2017, Earth Syst. Sci. Data, 12, 1561–1623,, 2020. 

Schwietzke, S., Sherwood, O. A., Bruhwiler, L. M. P., Miller, J. B., Etiope, G., Dlugokencky, E. J., Michel, S. E., Arling, V. A., Vaughn, B. H., White, J. W. C., and Tans, P. P.: Upward revision of global fossil fuel methane emissions based on isotope database, Nature, 538, 88–91,, 2016. 

Seibert, P. and Frank, A.: Source-receptor matrix calculation with a Lagrangian particle dispersion model in backward mode, Atmos. Chem. Phys., 4, 51–63,, 2004. 

Simpson, I. J., Sulbaek Andersen, M. P., Meinardi, S., Bruhwiler, L., Blake, N. J., Helmig, D., Rowland, F. S., and Blake, D. R.: Long-term decline of global atmospheric ethane concentrations and implications for methane, Nature, 488, 490–494,, 2012. 

Stohl, A. and Thomson, D. J.: A density correction for Lagrangian particle dispersion models, Bound.-Lay. Meteorol., 90, 155–167, 1999. 

Tadic, I., Crowley, J. N., Dienhart, D., Eger, P., Harder, H., Hottmann, B., Martinez, M., Parchatka, U., Paris, J.-D., Pozzer, A., Rohloff, R., Schuladen, J., Shenolikar, J., Tauer, S., Lelieveld, J., and Fischer, H.: Net ozone production and its relationship to nitrogen oxides and volatile organic compounds in the marine boundary layer around the Arabian Peninsula, Atmos. Chem. Phys., 20, 6769–6787,, 2020. 

Turner, A. J., Frankenberg, C., and Kort, E. A.: Interpreting contemporary trends in atmospheric methane, P. Natl. Acad. Sci. USA, 116, 2805–2813,, 2019. 

Wang, N., Edtbauer, A., Stönner, C., Pozzer, A., Bourtsoukidis, E., Ernle, L., Dienhart, D., Hottmann, B., Fischer, H., Schuladen, J., Crowley, J. N., Paris, J.-D., Lelieveld, J., and Williams, J.: Measurements of carbonyl compounds around the Arabian Peninsula: overview and model comparison, Atmos. Chem. Phys., 20, 10807–10829,, 2020. 

Whiticar, M.: Correlation of Natural Gases with Their Sources, in: The Petroleum System – From Source to Trap, edited by: Magoon, L. B. and Dow, W. G., AAPG, Tulsa, USA, 261–283, 1994. 

Yacovitch, T. I., Daube, C., and Herndon, S. C., Methane emissions from offshore oil and gas platforms in the Gulf of Mexico, Environ. Sci. Technol., 54, 3530–3538,, 2020. 

Yang, E. G., Kort, E. A., Wu, D., Lin, J. C., Oda, T., Ye, X., and Lauvaux, T.: Using space-based observations and Lagrangian modeling to evaluate urban carbon dioxide emissions in the Middle East, J. Geophys. Res., 125, e2019JD031922,, 2020. 

Yver Kwok, C., Laurent, O., Guemri, A., Philippon, C., Wastine, B., Rella, C. W., Vuillemin, C., Truong, F., Delmotte, M., Kazan, V., Darding, M., Lebègue, B., Kaiser, C., Xueref-Rémy, I., and Ramonet, M.: Comprehensive laboratory and field testing of cavity ring-down spectroscopy analyzers measuring H2O, CO2, CH4 and CO, Atmos. Meas. Tech., 8, 3867–3892,, 2015.  

Zavala-Araiza, D., Lyon, D. R., Alvarez, R. A., Davis, K. J., Harriss, R., Herndon, S. C., Karion, A., Kort, E. A., Lamb, B. K., Lan, X., Marchese, A. J., Pacala, S. W., Robinson, A. L., Shepson, P. B., Sweeney, C., Talbot, R., Townsend-Small, A., Yacovitch, T. I., Zimmerle, D. J., and Hamburg, S. P.: Reconciling divergent estimates of oil and gas methane emissions, P. Natl. Acad. Sci. USA, 112, 15597–15602,, 2015. 

Short summary
We measured atmospheric methane and CO2 by ship in the Middle East. We probe the origin of methane with a combination of light alkane measurements and modeling. We find strong influence from nearby oil and gas production over the Arabian Gulf. Comparing our data to inventories indicates that inventories overestimate sources from the upstream gas industry but underestimate emissions from oil extraction and processing. The Red Sea was under a complex mixture of sources due to human activity.
Final-revised paper