Next Article in Journal
CARNet: Context-Aware Residual Learning for JPEG-LS Compressed Remote Sensing Image Restoration
Next Article in Special Issue
Automatic Mapping of Burned Areas Using Landsat 8 Time-Series Images in Google Earth Engine: A Case Study from Iran
Previous Article in Journal
Slip Models of the 2016 and 2022 Menyuan, China, Earthquakes, Illustrating Regional Tectonic Structures
Previous Article in Special Issue
Impacts of the Urmia Lake Drought on Soil Salinity and Degradation Risk: An Integrated Geoinformatics Analysis and Monitoring Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Tailored Approach for the Global Gas Flaring Investigation by Means of Daytime Satellite Imagery

1
Institute of Methodologies for Environmental Analysis, National Research Council, 85050 Tito Scalo, Italy
2
School of Engineering, University of Basilicata, 85100 Potenza, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(24), 6319; https://doi.org/10.3390/rs14246319
Submission received: 27 October 2022 / Revised: 5 December 2022 / Accepted: 8 December 2022 / Published: 13 December 2022
(This article belongs to the Special Issue Natural Hazard Mapping with Google Earth Engine)

Abstract

:
The Daytime Approach for gas Flaring Investigation (DAFI), running in Google Earth Engine (GEE) environment, exploits a Normalized Hotspot Index (NHI), analyzing near-infrared and short-wave infrared radiances, to detect worldwide high-temperature gas flaring sites (GFs). Daytime Landsat 8—Operational Land Imager (OLI) observations, of 2013–2021, represents the employed dataset. A temporal persistence criterion is applied to a gas flaring customized NHI product to select the GFs. It assures the 99% detection accuracy of more intense and stable GFs, with a very low false positive rate. As a result, the first daytime database and map of GF sites, operating during the last 9 years at global scale, has been generated. For each site, geographical metadata, frequency of occurrence and time persistence levels, at both monthly and annual scale, may be examined, through the specific developed GEE App. The present database will complement/integrate existing gas flaring maps. The joint use of global scale daytime and nighttime GFs inventories, in fact, will allow for tracking gas flaring dynamics in a timely manner. Moreover, it enables a better evaluation of GF emissions into the atmosphere. Finally, the next DAFI implementation on Landsat 9 and Sentinel 2 data will further improve our capabilities in identifying, mapping, monitoring and characterizing the GFs.

Graphical Abstract

1. Introduction

Gas flaring sites (GFs) have a strong impact on Earth’s system [1,2,3], since they represent an important source of greenhouse gas (GHG) emission [4,5,6].
Controlling anthropogenic GHGs is important to mitigate global warming for human sustainable development [7]. Therefore, an accurate investigation of GFs is critical to evaluate the effects of emission reduction and control measures, as established by the World Bank’s “Zero Routine Flaring by 2030” initiative [8]. For this reason, it is essential to identify GFs at global scale and to monitor their operating conditions. However, “on site” GF measurements are generally scarce, and require extensive costs [9]. In addition, direct measurements and reports are often handled by the oil companies, and they are generally not public and thus hard to find and access to [10].
Satellite sensors, providing multispectral data at different spatial and temporal resolution, are largely used to investigate and monitor gas flaring sites, which are energy-intensive, heavy pollution industrial sectors, and large heat releasers. However, if on the one hand the spatial invariance and temporal persistence of GFs favor their discrimination from other high-temperature features (e.g., forest fires), on the other hand, these sources, operating under different environmental and climatic conditions, may not always be easy to identify by satellite.
To date, the major inventory of GFs has been realized by the Earth Observation Group (EOG), which developed the Visible Infrared Imaging Radiometer Suite (VIIRS) Nightfire algorithm (VNF) [11,12].
The algorithm exploits multi-spectral nighttime VIIRS data to identify and quantify gas flaring sources, with temperatures between 1600–2200 K, at a global scale. Parameters such as energy, flaring temperature and flared gas volume (https://eogdata.mines.edu/products/vnf/global_gas_flare.html, accessed on 26 October 2022) are generally provided together with the GF’s location.
However, the coarse spatial resolution of VIIRS data (hundred meters) does not enable an accurate characterization (in terms of location and spatial extent) of these high-temperature sources, which are highly radiant in the SWIR (short-wave infrared) region [13,14,15]. Mid–high (i.e., tens of meters) spatial resolution satellite data, from sensors such as OLI (Operational Land Imager) and MSI (Multispectral Instrument) may complement VIIRS observations [2]. Indeed, these sensors are largely used to detect and monitor high-temperature targets, associated with biomass burning and thermal volcanic activity [2,3,16,17,18,19,20].
More recently, daytime Landsat-8 (L8) OLI and Sentinel-2 (S2) MSI data have also been used to investigate gas flaring sources from space (e.g., [2,3,21]). Hot targets are easier to identify, in the SWIR band, in nighttime (as performed by VNF) than daylight conditions, because of the absence of the solar-reflected radiance from the background (e.g., [22]). Nevertheless, the joint usage of SWIR and NIR (near infrared, ~0.8 μm) bands may enable an efficient identification of gas flaring sources even in daylight conditions (e.g., [2,21]). Indeed, a high temperature source significantly increases the SWIR signal measured by satellite, whereas spectral radiance emitted by very hot targets is mostly insignificant in the NIR region (e.g., [3,23]).
By taking the aforementioned into account, a fixed threshold applied to a tri-spectral thermal anomaly index (TAI), exploiting NIR and SWIR S2-MSI reflectance images, was used to identify offshore gas flaring sites on a global scale, showing a detection reliability rate (visually determined) of about 88.0% [2]. The TAI also showed a great potential onshore (overall accuracy of >95% on Texas area), combining MSI and OLI observations [24]. A set of empirical thresholds, applied to SWIR to NIR ratios (both as reflectance and radiance) for L8-OLI and S2-MSI scenes, was proposed by [3], and automatically processed within Convolutional Neural Network and decision tree-based classification methods, showing an overall accuracy of 97.9% in detecting worldwide hot sources.
The Normalized Hotspot Indices (NHI) [25], guarantees a reliable identification and monitoring of thermal anomalies associated with onshore and offshore gas flaring activity, by using NIR and SWIR radiances, as demonstrated in [21], where six different test sites were investigated over time (i.e., 2015–2020).
In this work, a more in-depth analysis of these indices is performed. Moreover, a GF-customized configuration of NHI, named Daytime Approach for gas Flaring Identification (DAFI), and the first global GF sites inventory, obtained by analyzing daytime L8-OLI data collected during 2013–2021, are shown and discussed. Additionally, a Google Earth Engine (GEE) App, designed and implemented for the interactive analysis of the GF sites, is described.
The paper is organized as follows:
  • Section 2 describes the used datasets;
  • Section 3 explains how the DAFI works;
  • Section 4 and Section 5 detail the results achieved with the proposed method, along with the website hosting the Earth Engine App;
  • Section 6 concludes the work, discussing current challenges and future research directions.

2. Datasets

2.1. OLI-Landsat 8 Images

In this work, we used the L8 Collection 1 (C1) provided by the USGS (U.S. Geological Survey) and made available in GEE, which is a cloud-computation platform for processing large amounts of Earth Observation data [26]. Satellite scenes, placed in all C1 Tiers, were considered for the time-series processing analysis [27].
L8-OLI scenes, acquired in the NIR (Band 5: 0.85–0.88 μm), SWIR1 (Band 6: 1.57–1.65 μm) and SWIR2 (2.11–2.29 μm) bands, with the best Quality Assessment (QA) accuracy (QA = 9), and a sun elevation greater than zero degrees, were analyzed. A total of ~2.1 million scenes (for each band), covering the five continents in the temporal window 1 April 2013–31 December 2021 represents the investigated dataset.
The radiance values, obtained for each pixel, by applying the “ee.Algorithms.Landsat.calibratedRadiance” method within Earth Engine [28], were used to compute the DAFI-based metrics.

2.2. Other Datasets

The following additional datasets were used:
(i)
High-spatial resolution optical satellite data (up to a few meters, depending on data sources [2]), from Google Earth and Google Maps, providing spatially detailed structures and distinguishable marks of different industrial categories (as support and validation of industrial hotspot classification);
(ii)
Global industrial heat sources inventory, from VNF (~750 m) time series to compare two independent satellite products, operating in opposite solar lighting conditions (i.e., DAFI by day, and VNF by night). The VNF detections from 2012 to 2021 were downloaded as KML files from https://eogdata.mines.edu/products/vnf/global_gas_flare.html (accessed on 6 July 2022) and clustered both spatially and temporally, to enhance their visualization.

3. The DAFI Method

The NHI algorithm was developed by [25] to detect and map volcanic thermal anomalies, through the following normalized indices:
  NHI SWIR = L SWIR 2   L SWIR 1 / L SWIR 2 + L SWIR 1
NHI SWNIR = L SWIR 1   L NIR / L SWIR 1 + L NIR
where LSWIR2, LSWIR1, LNIR are the TOA (Top of Atmosphere) radiances (W m−2 sr−1 μm−1) measured, for each pixel of the scene, at 2.2 μm, 1.6 μm and 0.8 μm wavelengths, in the relative OLI/MSI spectral channels. Pixels showing positive values of one or both indices are flagged as “hot” by NHI [25].
While positive values of the index in Equation (1) enable the identification of hot targets emitting strongly in the SWIR2 band, positive values of the index in Equation (2) guarantee a more accurate identification of intense thermal features (e.g., GFs more radiant in the SWIR1 band), which frequently saturate the OLI/MSI SWIR2 channel.
A further step is used to recover the very intense thermal features, named “Extreme Pixel”, yielding the SWIR1 band saturation (see [29]).
As previously demonstrated in [25], and following studies (e.g., [29]), the NHI algorithm enables an effective discrimination of hot targets from clouds, the latter leading to negative values of both indices. Therefore, we did not filter out cloudy images. This is an important feature of the NHI algorithm, which should guarantee a low false positive rate at global scale (e.g., [19]), even when GF sites are analyzed. Indeed, cloud detection is an unresolved issue in the thermal anomaly identification in daylight. If on the one side a cloud-masking is generally required to filter out clouds in the SWIR bands, on the other side, gas flaring sites may be undetected if clouds are overestimated; moreover, hot targets may be mislabeled as clouds in clear-sky conditions (e.g., [2,30]).
The “Test for Hot Spot Pixels” and the “Test for Extreme Pixels” described in [29], implemented operationally in GEE, enable the identification of a Hot Pixels (HP) according to the following criterion:
NHI SWIR > 0   OR   NHI SWNIR > 0   OR   Extreme   Pixel > 0 = 1   H P otherwise = 0   no   HP
The algorithm (which also includes an initial test on SWIR2 radiance [29]) assigns the flag “1” (i.e., TRUE) to a pixel when at least one of the three indicators (NHISWNIR, NHISWIR, Extreme Pixel) is greater than zero and the flag “0” (i.e., FALSE) is otherwise. Therefore, for each scene of the satellite collection, the final output is a binary mask, where “1” means a hot pixel.
This detection scheme was first applied to detect GFs in [21]; in this work, the OLI-based NHI indices are deeply investigated to develop a gas flaring tailored approach.
From a satellite perspective, gas flaring sources may exhibit a similar spectral behavior when compared to fires or active lava flows, but differently from them, they are stable/fixed in space and (mostly) continuously operating in time (except for maintenance/interruption periods) [31]. Therefore, the binary masks computed for each pixel, at 30 m spatial resolution, in the globe have been summed up for the given time interval (i.e., 2013–2021), and an HP temporally cumulative mask obtained. As volcanoes may also experience long-term activities, a volcanic mask was defined to exclude possible volcanic thermal anomalies from the performed analyses. A buffer area, with a radius ranging from 3 to 30 km, depending on the volcanic cone extension, was then applied to the volcanoes listed by the Global Volcanism Program (https://volcano.si.edu/volcanolist_holocene.cfm, accessed on 1 June 2022). This exclusion map is expected to not affect the sensitivity in GF identification at all, as gas flares are certainly not installed within (or even close to) active volcanic areas.
The derived values in the cumulative mask were converted in an occurrence frequency (OF), to weigh the NHI contribution with respect to the number of OLI available images, depending on the revisit cycle of the satellite [32]. For each pixel, the OF, expressed as a percentage, is computed as:
OF = 2013 2021 HP N       [ % ]
where N counts all (both cloudy and clear) OLI images acquired in the investigated time period.
To select only stable (in space and time) hot targets, we collected the OF values retrieved in correspondence with a GF testing dataset (recognizable in Google Earth). A minimum threshold value of 10% was found and set to identify GFs only.
The final OF raster map was vectorized as a centroid geometry. As a result, 2084 candidate GF sites were identified (Figure 1). Each site may correspond to a single OLI pixel (30 × 30 m) or to a cluster of adjacent and contiguous pixels, which satisfy the imposed condition (i.e., OF ≥ 10%). To better describe the GF thermal behavior, a 100 m buffer (3-pixel), around the detected centroids was used to recover the OF maximum value within the cluster. In this way, a more precise localization of each gas flare site is guaranteed.
To assess the efficiency of the initial assumption, we performed a careful visual inspection of the detected GF candidates, by using the high-resolution images available in Google Earth. The detected centroids were classified into five categories (see Figure 1): gas flares (37%), industrial sites (44%) (including hotspots related to incineration plants, landfill, gas power plants, metal processing factories, liquefied natural gas terminals, chemical factories, wood production factories, wastewater treatment plants, cement factories), quarries/mines (7%), nude soils (2%, mostly in Australia) and rooftops (10%, mostly distributed in United States and Japan). While thermal anomalies could also be associated with quarries/mines (due to sewage/wastewater treatment plants), the last two categories (i.e., nude soils and rooftops) clearly indicate false positives. Therefore, this classification, retrieved from the standard NHI configuration did not enable the discrimination of gas flares from other heat sources (i.e., industrial sites) with a suitable accuracy. Moreover, a few false positives (e.g., in correspondence with rooftops) were generated, requiring the use of an ancillary dataset (e.g., land cover or land use map). Based on those results, we carried out a more in-depth analysis of the NHI indices. In particular, the HP value (the one computed by summing the flags “1” over the years 2013–2021), measured for each of the 2084 centroids, was broken down to weigh each single indicator. Table 1 provides an overview of all possible indicator combinations, from I to VI. The used notation refers to the HP’s sum over the analyzed time period, while the subscript describes the indicator used to identify the hotspot (e.g., “EP” stands for Extreme Pixel).
A hot pixel is then flagged when all the indicators are positive (condition I), two of them are positive (conditions II–III) or only a single indicator is positive (conditions IV–VI).
By taking into account these conditions, the 2084 centroids were grouped (see second column of Table 2). Additionally, the single centroid of each condition was attributed to the corresponding category (e.g., gas flares, industrial sites) after overlapping the centroids on the very high-resolution images of Google Earth.
Table 2 shows that each condition discriminates GFs from other hot sources (e.g., industrial sites) and highly reflective (e.g., rooftops) targets in a different way (see the third column); the condition VI is not statistically significant, as it refers to 4 sites only (see the second column). Therefore, if a high confidence level is aimed at, the condition II (highlighted in grey in Table 2) should be selected. This condition, which takes into account the peak of thermal emissions at 1.6 μm and the possible saturation of SWIR1 channel in presence of gas flaring sources (e.g., [6,11,13]), guarantees an overall detection reliability around 99%. The latter was estimated by the spatial correspondence with the centroids with the gas flaring sites, revealed by satellite images at very high-spatial resolution. Moreover, condition II seems to automatically filter out highly reflective features (e.g., roofing materials [2,3]). The latter may lead to high radiance values in the SWIR2 band, which may be misinterpreted as hot pixels, generating commission errors (e.g., condition V in Table 2). However, it should be pointed out that gas flares operating at lower temperature (less than 1400 K) may be discarded using the condition II. Indeed, this condition makes the database biased towards high-temperature gas flaring sources. In addition, a threshold based on the temporal persistence of the signals may further reduce the number of detectable GFs.
Consequently, the NHI algorithm was rearranged by developing the “Daytime Approach for gas Flaring Investigation” (DAFI) method. This method consists of a three-module cascade shown in Figure 2, where “HPdaily” refers to hot pixels identified using the condition II of Table 2. Equation (4) is used to sum up the HPdaily values and compute the OF for each pixel.
Hence, DAFI works on a single OLI scene searching for hot pixels through the NHISWNIR index and the test on saturated pixels in the SWIR1 band (see module (a)). This operation was performed for all the scenes of the time series, preserving the “detection history” of each location. The OF was calculated from the HPdaily masks and the number of analyzed images, summed in the investigated time window. Pixels showing OF values ≥ 10% were classified as gas flares (module (b)). Finally, the persistence level of each GF was computed through a quartiles analysis, which allowed us to define four GF persistence levels (module (c)), shown in Table 3. The lower the level, the less persistent the activity.

4. Results

The DAFI inventory consists of 1711 high-temperature gas flaring sites, identified over the period 2013–2021 (see Figure 3).
The distribution of the global GF sites is spatially concentrated: 1237 points are located onshore (72%), while the remaining 474 (28%) were flagged offshore. About 70% of the onshore GFs belongs to Asia (~50%, located in Russian Federation, Iran and Iraq), and ~15% to Africa (Algeria, South Africa) and Americas (United States of America, Mexico, Venezuela). The distribution of the global offshore GFs is concentrated into the North Sea (~9%), Gulf of Guinea (~22%), Persian Gulf (~15%) and South China Sea (~20%); they represent about ~67% of the total [30]. As far as both onshore and offshore GFs are considered together, the Asia continent is the most populated (Figure 3), followed by the American and African ones, where GF sites show a quite similar distribution in both offshore and onshore conditions. Europe and Oceania have a marginal role in terms of oil and gas facilities (less than ~15% and ~5%, respectively).
The countries more affected by the gas flaring sites, selected as the ones accounting for at least the 50% of overall continental detections, are listed in Figure 4.
The Russian Federation, Iran and Iraq occupy the podium, accounting for about 40% of the worldwide detected GFs, followed by Algeria, Mexico and USA. These results are in good agreement with the last report of the Global Gas Flaring Reduction Partnership [8], stating that Russia, Iraq, Iran, the United States, Algeria, Venezuela and Nigeria are among the top 10 flaring countries. The latter have held this position consistently for the last 10 years, accounting for roughly two-thirds (65%) of global gas flaring.
Figure 5 shows the distribution of the persistence levels for both onshore and offshore GFs categories, emphasizing their different behavior in terms of operability.
In particular, the figure reveals that while the onshore sites have a mid–low persistence, the marine platforms are the most operational, with about 50% of them belonging to the high-persistence class.
Flaring may occur in the oil and gas industry for many reasons, ranging from initial start-up testing of a facility to unplanned equipment malfunctions. The combustion in gas flares can be almost constant over the time, but can also be very variable, hence the flaring can be reduced or turned off. Moreover, a gas flare may be located on a sea platform or on the ground flares with open pits and multiple nozzles, and elevated individual stacks [33]. All these features affect the time duration of the process, which needs to be considered for detection purposes. Thus, thermal anomalies identified on each site by DAFI have been temporally disaggregated, on an annual and monthly basis, to acquire some information about the source type (i.e., continuous or non-continuous flaring). Some examples are shown in Section 5.
In sum, the DAFI results, retrieved from daytime L8-OLI data, provide an accurate information about the global distribution, localization and mapping of the most intense gas flaring sites.

5. The DAFI Website

The DAFI information can be accessed online, at the website https://sites.google.com/view/flaringsitesinventory (accessed on 26 October 2022), without any authentication. Data can be visualized and navigated at two levels: static and dynamic.
The first level (named “The global gas flaring DAFI map”) allows users to visualize the spatial distribution of the gas flaring sites, navigating a custom map created using Google My Maps (Figure 6). By clicking on the single site of interest, the following information appears on the left-side of the map (Figure 6):
-
Continent, region, country;
-
GF site longitude and latitude (in decimal degree format);
-
Occurrence Frequency, OF (%);
-
Long-term Persistence Level (from low to high, depending on the OF value);
-
Shore (gas flare location, classified in onshore and offshore).
An identification code is assigned to each gas flaring site. It consists of the country acronym (i.e., Alpha-3 country code) followed by a 4-digit progressive number (e.g., DZA_0357, see example in Figure 6).
The dynamic characterization of the GF sites runs automatically in GEE. A schematic picture of the developed DAFI App is shown in Figure 7.
By approaching this level, named “The DAFI App for gas flaring thermal investigation”, two steps are required to investigate the gas flaring source:
(1)
Choose the gas flare to inspect.
(2)
Run the DAFI to get the GF behavior. As previously mentioned, by exploring the map, users can select the gas flare of interest. By directly obtaining the ID number (in the drop-down menu) or copying/tasting it from the map the time series of SWIR1 and NIR radiances, measured over the selected GF site, is displayed. Thermal anomalies that are aggregated both in terms of years and months (i.e., OF value during a given calendar year/month), can be visualized through line charts, highlighting the periods in which the thermal activity was dominant, increased or declined at site level. Each graph can be downloaded in CSV, SVG or PNG format.
Some examples of the tool outputs are shown below and refer to both onshore and offshore sites (Figure 6). The criterion adopted to select them is the following:
  • GF sites belonging to the same persistence category but located in different countries;
  • GF sites located in the same country but showing different persistence category.
This analysis may provide interesting information about differences/similarities in the GF dynamics and thermal behavior. For each site, a picture displaying the GF-affected area (white circle), with the thermally anomalous pixels, depicted in different colors according to their OF value, is shown, along with the annual and monthly temporal profile chart of the GF occurrence.
Figure 8 displays the GF sites belonging to the first aforementioned category (i.e., mid–high): DZA_0313 and RUS_1310.
Figure 8 shows that some differences characterize the GFs sites within the same long-term persistence class. The Algerian gas flare (DZA_0313), with a 12-pixel cluster (Figure 8a), seems to operate with an increasing trend (with percentages of occurrence lower than 22%) up to 2019, where the peak is recorded (~61%), followed by a reduced continuity in the next two years (see Figure 8b, purple line). The Russian site (RUS_1310), composed by three pixels (Figure 8c), shows a quite stable flaring activity between 2013–2018, apart from a slight peak recorded in 2015 (Figure 8b, purple dashed line), while it halves in the following years. Additionally, at the monthly scale, the OFs distribution of the two sites follows an opposite trend, which is a convex trend for the Algerian flare and quite concave for the Russian one (see Figure 8d).
By looking at the selected offshore gas flaring sites, both are classified as mid–high persistence objects (Figure 9), located in the Caspian Sea (TKM_1518, Figure 9a) and Persian Gulf (ARE_0041, Figure 9c), 30 km and 70 km from the Turkmen and United Arab Emirates shores. Both the marine platforms, part of a flaring sites archipelago, are clearly visible from space.
The oil and gas platform in the Caspian Sea has not been operational since 2017 (purple line in Figure 9b), with an unstable operation in the previous years. The gas flaring has been routinely practiced at the ARE_0041 (purple dashed line in Figure 9b) from 2014 to 2020, with two peaks in 2016 and 2018. As the previous one, it stopped in 2021. At a monthly scale (Figure 9d), the two sites work over the whole 12-month period (except in March for TKM_1518), with a higher persistence level for the Arab platform in most of the months (higher than 10%).
Regarding the second category detailed above, the Iraqi sites, IRQ_0795 and IRQ_0766, flagged, respectively, as high and low long-term persistence sources (the cluster size likely reflects this), are represented in Figure 10a,c.
As highlighted by the line chart in Figure 10b, the level of annual persistence is significantly different between the two sites over the analyzed years (with a mean bias of 50%). Moreover, observing the IRQ_0795 “high” source, its persistence is gradually increasing during the first four years (2013–2016), with a growth rate of over 70% (the peak is reached in 2016, around 83%) (Figure 10b, purple line). Since 2016, the level declines for touching the starting one in 2018, when a new rise is recorded. IRQ_0766 appears quite stable, with OF percentages ranging from 5% to 20% over the years 2013–2019, with a peak of 30% in 2020 (Figure 10b, purple dashed line). The monthly profiles of the Iraqi sites’ occurrences (Figure 10d) indicate that the gas flaring activity mainly occurs between September and February/March. As expected, the most intense levels (more than 50% on average) can be observed in correspondence with the “high” source.
Figure 10d puts in evidence the same concave trend also characterizing the Algerian and Russian sites (see Figure 8d). Further investigations are required to better explain those observations.

6. Conclusions

DAFI uses a temporal persistence criterion, applied to a customized NHI-based indicator, to identify the GF sites which were active over the last nine years (i.e., 2013–2021). This way, DAFI is capable of detecting very intense and persistent gas flaring sources, allowing for a deep characterization, in terms of temporal dynamics (at both annual and monthly scale) and spatial spreading of hot pixels at the local scale.
As a final result, the first global map of high-temperature gas flaring sites based on daytime OLI-L8 observations was produced. A GF detection reliability of 99% was achieved by setting the 10% threshold persistence.
It should be pointed out that DAFI information is not exhaustive, owing to the imposed restrictions (see Section 3). Indeed, while commission errors, associated only with other industrial plants, are negligible and do not represent an issue, the GF database mainly includes the high-temperature gas flares operating in time. In particular, missed detections may occur in the presence of mid–low intensity flares (e.g., GF operating at lower temperatures, and emitting more strongly in the SWIR2 than SWIR1 band) and/or of extremely recent (or presently dismissed) flares, which do not meet the selected criteria (i.e., OF slightly less than 10%). Both factors may affect the completeness of the GF dataset. Moreover, these cases may be further emphasized in case of cloud cover, which hinders or masks the gas flare signals (especially the less intense ones). Moreover, the temporal investigation of the GF occurrence frequency may be misunderstood in geographic areas where clouds persist in some months of the year. In such a case, a measured null (or very small) value of OF maynot be directly related to the gas flare operational features but to the cloud coverage, preventing an effective gas flare site characterization in time domain.
Despite those limitations, the DAFI global daytime map may provide additional information about GFs when integrated with the nighttime VIIRS observations. Figure 11 displays the spatial distribution of the GFs from the nighttime VNF method, independently provided by the EOG, and the daytime DAFI one.
A general satisfying spatial agreement between daytime and nighttime GF datasets is evident, with the flaring regions identified by DAFI which are also captured by VNF. However, the different spatial resolution, revisit interval, acquisition time and spectral bands used by the two methods make a point-by-point comparison between them impracticable (and even out of the scope of this work). On the other hand, their integration may represent an added value in terms of detection, monitoring and characterization capabilities of GF from space. Indeed, daytime and nighttime observations, if integrated, by increasing the temporal sampling, may allow for a more effective identification of the GFs. Moreover, the different GF typology and characteristics may be better revealed when different spectral/spatial and temporal domains are investigated. A visual investigation of DAFI and VNF maps has revealed relative performance and limitations, as put in evidence by some examples reported in Figure 12 showing:
o
The higher sensitivity of OLI to small-scale hot targets, undetected by VIIRS (e.g., Figure 12a);
o
The greater spatial accuracy offered by OLI than VIIRS in localizing and mapping the GF areas, useful information for accurately characterizing the site in terms of flared volumes and GHG emissions (Figure 12b);
o
The lower sensitivity of DAFI than VNF in detecting some GFs because of the selected persistence threshold (Figure 12c) (shale oil/gas included).
Figure 12 clearly shows the advantage of a gas flaring mapping from a complementary satellite-based inventory.
The current DAFI map will be updated annually to follow the GF spatiotemporal dynamics. While the aforementioned limitations probably led to an underestimation of gas flares in this study, they will be likely (at least partially) overcome by exporting DAFI on additional satellite data at mid–high spatial resolution. In particular, the DAFI porting on Landsat 9 and Sentinel 2A/B (at 20m spatial resolution) data is under development, by taking into account the preliminary outcomes shown in [21]. The combination of Sentinel-2A/B MSI and L8/9 OLI/OLI-2 data would guarantee an increased frequency of observation, up to a revisit time of 1–3 days [34,35,36,37], with a positive impact on the detection sensitivity. The OF percentage levels lower than 10% will also be explored using L9 and S2 collections. Moreover, the different behavior of the onshore and offshore sites, observed by applying the long-term persistence class, will be better analyzed, evaluating the DAFI performance in recognizing hot pixels over different areas. Estimates of the radiative power, for each gas flare, will be also performed and made available on the DAFI website. This metric will be crucial to carefully retrieve both the volume of gas flared and the CO2 equivalent emissions by each identified GF over one calendar year/month.
Given the growing international focus on global methane emissions, and acknowledging the important role the oil and gas industry could play towards achieving climate goals [8], a global overview of the GF sites, as exhaustive as possible, becomes crucial for monitoring the countries’ progress in decarbonization strategies [8].
Considering that bottom-up emission inventories for the oil and gas industry are often fragmentary, due to challenges in quantifying the emissions from gas flaring or to incomplete reporting by individual companies, satellite-based observations may then be offering a valuable alternative to independently acquire more complete, accurate and timely emission estimates [7].
A global GF emissions scenario that is as exhaustive as possible may help governments and operators in tracking commitments and efforts in oil-producing countries to abate flaring, as well as in conducting air quality monitoring for better understanding of exposures, and reducing potential human health impacts [31,38].

Author Contributions

Conceptualization, M.F., N.G. and N.P.; methodology, M.F. and F.M.; formal analysis, M.F. and N.G.; investigation, M.F.; writing—original draft preparation, M.F.; writing—review and editing, F.M. and N.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ismail, O.S.; Umukoro, G.E. Global impact of gas flaring. Energy Power Eng. 2012, 4, 290–302. [Google Scholar] [CrossRef] [Green Version]
  2. Liu, Y.; Zhi, W.; Xu, B.; Xu, W.; Wu, W. Detecting high-temperature anomalies from Sentinel-2 MSI images. ISPRS J. Photogramm. Remote Sens. 2021, 177, 174–193. [Google Scholar] [CrossRef]
  3. Kato, S.; Miyamoto, H.; Amici, S.; Oda, A.; Matsushita, H.; Nakamura, R. Automated classification of heat sources detected using SWIR remote sensing. Int. J. Appl. Earth Obs. Geoinf. 2021, 103, 102491. [Google Scholar] [CrossRef]
  4. Fawole, O.G.; Cai, X.-M.; MacKenzie, A. Gas flaring and resultant air pollution: A review focusing on black carbon. Environ. Pollut. 2016, 216, 182–197. [Google Scholar] [CrossRef]
  5. Caseiro, A.; Gehrke, B.; Rücker, G.; Leimbach, D.; Kaiser, J.W. Gas flaring activity and black carbon emissions in 2017 derived from the Sentinel-3A Sea and Land Surface Temperature Radiometer. Earth Syst. Sci. Data 2020, 12, 2137–2155. [Google Scholar] [CrossRef]
  6. Faruolo, M.; Caseiro, A.; Lacava, T.; Kaiser, J.W. Gas Flaring: A Review Focused on Its Analysis from Space. IEEE Geosci. Remote Sens. Mag. 2021, 9, 258–281. [Google Scholar] [CrossRef]
  7. Ialongo, I.; Stepanova, N.; Hakkarainen, J.; Virta, H.; Gritsenko, D. Satellite-based estimates of nitrogen oxide and methane emissions from gas flaring and oil production activities in Sakha Republic, Russia. Atmos. Environ. X 2021, 11, 100114. [Google Scholar] [CrossRef]
  8. Global Gas Flaring Reduction Partnership (GGFR)—The World Bank. 2022 Global Gas Flaring Tracker Report. Available online: https://thedocs.worldbank.org/en/doc/1692f2ba2bd6408db82db9eb3894a789-0400072022/original/2022-Global-Gas-Flaring-Tracker-Report.pdf (accessed on 3 October 2022).
  9. Giwa, S.O.; Nwaokocha, C.N.; Kuye, S.I.; Adama, K.O. Gas flaring attendant impacts of criteria and particulate pollutants: A case of Niger Delta region of Nigeria. J. King Saud Univ. Eng. Sci. 2019, 31, 209–217. [Google Scholar] [CrossRef]
  10. Zhang, Z.; Sherwin, E.D.; Brandt, A.R. Estimating global oilfield-specific flaring with uncertainty using a detailed geographic database of oil and gas fields. Environ. Res. Lett. 2021, 16, 124039. [Google Scholar] [CrossRef]
  11. Elvidge, C.D.; Zhizhin, M.; Hsu, F.-C.; Baugh, K.E. VIIRS Nightfire: Satellite pyrometry at night. Remote Sens. 2013, 5, 4423–4449. [Google Scholar] [CrossRef]
  12. Elvidge, C.D.; Zhizhin, M.; Baugh, K.E.; Hsu, F.-C.; Ghosh, T. Methods for global survey of natural gas flaring from visible infrared imaging radiometer suite data. Energies 2016, 9, 14. [Google Scholar] [CrossRef]
  13. Fisher, D.; Wooster, M.J. Shortwave IR adaption of the mid-infrared radiance method of fire radiative power (FRP) retrieval for assessing industrial gas flaring output. Remote Sens. 2018, 10, 305. [Google Scholar] [CrossRef] [Green Version]
  14. Elvidge, C.D.; Zhizhin, M.; Baugh, K.; Hsu, F.; Ghosh, T. Extending nighttime combustion source detection limits with short wavelength VIIRS data. Remote Sens. 2019, 11, 395. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, P.; Yuan, C.; Sun, Q.; Liu, A.; You, S.; Li, X.; Zhang, Y.; Jiao, X.; Sun, D.; Sun, M.; et al. Satellite-Based Detection and Characterization of Industrial Heat Sources in China. Env. Sci Technol 2019, 53, 11031–11042. [Google Scholar] [CrossRef]
  16. Hu, X.; Ban, Y.; Nascetti, A. Sentinel-2 MSI data for active fire detection in major fire-prone biomes: A multi-criteria approach. Int. J. Appl. Earth Obs. Geoinf. 2021, 101, 102347. [Google Scholar] [CrossRef]
  17. Ngadze, F.; Mpakairi, K.S.; Kavhu, B.; Ndaimani, H.; Maremba, M.S. Exploring the utility of Sentinel-2 MSI and Landsat 8 OLI in burned area mapping for a heterogenous savannah landscape. PLoS ONE 2020, 15, e0232962. [Google Scholar] [CrossRef] [PubMed]
  18. Marchese, F.; Genzano, N.; Nolde, M.; Falconieri, A.; Pergola, N.; Plank, S. Mapping and characterizing the Kīlauea (Hawaii) lava lake through Sentinel-2 MSI and Landsat-8 OLI observations of December 2020–February 2021. Environ. Model. Softw. 2022, 148, 105273. [Google Scholar] [CrossRef]
  19. Marchese, F.; Genzano, N. Global volcano monitoring through the NHI (Normalized Hotspot Indices) system. J. Geol. Soc. 2022, 180, jgs2022-014. [Google Scholar] [CrossRef]
  20. Mazzeo, G.; Ramsey, M.S.; Marchese, F.; Genzano, N.; Pergola, N. Implementation of the NHI (Normalized Hot Spot Indices) Algorithm on Infrared ASTER Data: Results and Future Perspectives. Sensors 2021, 21, 1538. [Google Scholar] [CrossRef] [PubMed]
  21. Faruolo, M.; Falconieri, A.; Genzano, N.; Lacava, T.; Marchese, F.; Pergola, N. A Daytime Multisensor Satellite System for Global Gas Flaring Monitoring. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–17. [Google Scholar] [CrossRef]
  22. Schroeder, W.; Prins, E.; Giglio, L.; Csiszar, I.; Schmidt, C.; Morisette, J.; Morton, D. Validation of GOES and MODIS active fire detection products using ASTER and ETM + data. Remote Sens. Environ. 2008, 112, 2711–2726. [Google Scholar] [CrossRef]
  23. Giglio, L.; Csiszar, I.; Restás, Á.; Morisette, J.T.; Schroeder, W.; Morton, D.; Justice, C.O. Active fire detection and characterization with the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER). Remote Sens. Environ. 2008, 112, 3055–3063. [Google Scholar] [CrossRef]
  24. Wu, W.; Liu, Y.; Rogers, B.M.; Xu, W.; Dong, Y.; Lu, W. Monitoring gas flaring in Texas using time-series sentinel-2 MSI and landsat-8 OLI images. Int. J. Appl. Earth Obs. Geoinf. 2022, 114, 103075. [Google Scholar] [CrossRef]
  25. Marchese, F.; Genzano, N.; Neri, M.; Falconieri, A.; Mazzeo, G.; Pergola, N. A multi-channel algorithm for mapping volcanic thermal anomalies by means of sentinel-2 MSI and Landsat-8 OLI data. Remote Sens. 2019, 11, 2876. [Google Scholar] [CrossRef] [Green Version]
  26. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  27. Google Earth Engine Landsat Algorithms. Available online: https://developers.google.com/earth-engine/guides/landsat#landsatcollection-structure (accessed on 3 October 2022).
  28. Chander, G.; Markham, B.L.; Helder, D.L. Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors. Remote Sens. Environ. 2009, 113, 893–903. [Google Scholar] [CrossRef]
  29. Genzano, N.; Pergola, N.; Marchese, F. A Google Earth engine tool to investigate, map and monitor volcanic thermal anomalies at global scale by means of mid-high spatial resolution satellite data. Remote Sens. 2020, 12, 3232. [Google Scholar] [CrossRef]
  30. Lu, W.; Liu, Y.; Wang, J.; Xu, W.; Wu, W.; Liu, Y.; Zhao, B.; Li, H.; Li, P. Global proliferation of offshore gas flaring areas. J. Maps 2020, 16, 396–404. [Google Scholar] [CrossRef]
  31. Franklin, M.; Chau, K.; Cushing, L.J.; Johnston, J.E. Characterizing Flaring from Unconventional Oil and Gas Operations in South Texas Using Satellite Observations. Environ. Sci Technol. 2019, 53, 2220–2228. [Google Scholar] [CrossRef]
  32. Wulder, M.A.; White, J.C.; Loveland, T.R.; Woodcock, C.E.; Belward, A.S.; Cohen, W.B.; Fosnight, E.A.; Shaw, J.; Masek, J.G.; Roy, D.P. The global Landsat archive: Status, consolidation, and direction. Remote Sens. Environ. 2016, 185, 271–283. [Google Scholar] [CrossRef]
  33. Labeyrie, H.; Rocher, A. Reducing Flaring and Improving Energy Efficiency: An Operator’s View. In Proceedings of the SPE International Conference on Health, Safety and Environment in Oil and Gas Exploration and Production, Rio de Janeiro, Brazil, 12 April 2010. [Google Scholar]
  34. Li, J.; Chen, B. Global Revisit Interval Analysis of Landsat-8 -9 and Sentinel-2A -2B Data for Terrestrial Monitoring. Sensors 2020, 20, 6631. [Google Scholar] [CrossRef] [PubMed]
  35. Li, J.; Roy, D.P. A Global Analysis of Sentinel-2A, Sentinel-2B and Landsat-8 Data Revisit Intervals and Implications for Terrestrial Monitoring. Remote Sens. 2017, 9, 902. [Google Scholar] [CrossRef] [Green Version]
  36. Jiang, R.; Sanchez-Azofeifa, A.; Laakso, K.; Xu, Y.; Zhou, Z.; Luo, X.; Huang, J.; Chen, X.; Zang, Y. Cloud Cover throughout All the Paddy Rice Fields in Guangdong, China: Impacts on Sentinel 2 MSI and Landsat 8 OLI Optical Observations. Remote Sens. 2021, 13, 2961. [Google Scholar] [CrossRef]
  37. Ustin, S.L.; Middleton, E.M. Current and near-term advances in Earth observation for ecological applications. Ecol. Process. 2021, 10, 1–57. [Google Scholar] [CrossRef]
  38. Zhizhin, M.; Matveev, A.; Ghosh, T.; Hsu, F.-C.; Howells, M.; Elvidge, C. Measuring Gas Flaring in Russia with Multispectral VIIRS Nightfire. Remote Sens. 2021, 13, 3078. [Google Scholar] [CrossRef]
Figure 1. Candidate GF sites identified by setting OF ≥ 10% (see text).
Figure 1. Candidate GF sites identified by setting OF ≥ 10% (see text).
Remotesensing 14 06319 g001
Figure 2. DAFI conditions for (a) detecting a potential sub-pixel hotspot over the single daytime OLI image (first module); (b) flagging the pixel as a gas flaring site (second module); (c) assigning a persistence level using a long-term temporal scale to the GF (third module).
Figure 2. DAFI conditions for (a) detecting a potential sub-pixel hotspot over the single daytime OLI image (first module); (b) flagging the pixel as a gas flaring site (second module); (c) assigning a persistence level using a long-term temporal scale to the GF (third module).
Remotesensing 14 06319 g002
Figure 3. Distribution of GF sites detected by DAFI (purple dots) among the five continents. The onshore and offshore percentages computed for each continent is also stated.
Figure 3. Distribution of GF sites detected by DAFI (purple dots) among the five continents. The onshore and offshore percentages computed for each continent is also stated.
Remotesensing 14 06319 g003
Figure 4. The worldwide countries including the greatest number of detected gas flares. Each percentage indicates the ratio of the national detections over the total detected gas flares (i.e., 1711).
Figure 4. The worldwide countries including the greatest number of detected gas flares. Each percentage indicates the ratio of the national detections over the total detected gas flares (i.e., 1711).
Remotesensing 14 06319 g004
Figure 5. OFs distribution among the four persistence classes for the (a) onshore and (b) offshore sites.
Figure 5. OFs distribution among the four persistence classes for the (a) onshore and (b) offshore sites.
Remotesensing 14 06319 g005aRemotesensing 14 06319 g005b
Figure 6. The global gas flaring DAFI map: onshore (orange flames) and offshore (blue flames) gas flaring sites. On the left side, a panel summarizing the GF features (e.g., DZA_0357).
Figure 6. The global gas flaring DAFI map: onshore (orange flames) and offshore (blue flames) gas flaring sites. On the left side, a panel summarizing the GF features (e.g., DZA_0357).
Remotesensing 14 06319 g006
Figure 7. DAFI GEE App picture.
Figure 7. DAFI GEE App picture.
Remotesensing 14 06319 g007
Figure 8. Features of the selected onshore GF sites, labeled as (a) DZA_0313 and (c) RUS_1310, retrieved from the DAFI GEE tool; the insets summarize the GF geographic information, along with its long-term persistence category. The (b) annual and (d) monthly temporal persistence for DZA_0313 (purple line) and RUS_1310 (purple dashed line).
Figure 8. Features of the selected onshore GF sites, labeled as (a) DZA_0313 and (c) RUS_1310, retrieved from the DAFI GEE tool; the insets summarize the GF geographic information, along with its long-term persistence category. The (b) annual and (d) monthly temporal persistence for DZA_0313 (purple line) and RUS_1310 (purple dashed line).
Remotesensing 14 06319 g008
Figure 9. Features of the selected onshore GF sites, labeled as (a) TKM_1518 and (c) ARE_0041, retrieved from the DAFI GEE tool; the insets summarize the GF geographic information, along with its long-term persistence category. The (b) annual and (d) monthly temporal persistence for TKM_1518 (purple line) and ARE_0041 (purple dashed line).
Figure 9. Features of the selected onshore GF sites, labeled as (a) TKM_1518 and (c) ARE_0041, retrieved from the DAFI GEE tool; the insets summarize the GF geographic information, along with its long-term persistence category. The (b) annual and (d) monthly temporal persistence for TKM_1518 (purple line) and ARE_0041 (purple dashed line).
Remotesensing 14 06319 g009
Figure 10. Features of the selected onshore GF sites, labeled as (a) IRQ_0795 and (c) IRQ_0766, retrieved from the DAFI GEE tool; the insets summarize the GF geographic information, along with its long-term persistence category. The (b) annual and (d) monthly temporal persistence for IRQ_0795 (purple line) and IRQ_0766 (purple dashed line).
Figure 10. Features of the selected onshore GF sites, labeled as (a) IRQ_0795 and (c) IRQ_0766, retrieved from the DAFI GEE tool; the insets summarize the GF geographic information, along with its long-term persistence category. The (b) annual and (d) monthly temporal persistence for IRQ_0795 (purple line) and IRQ_0766 (purple dashed line).
Remotesensing 14 06319 g010aRemotesensing 14 06319 g010b
Figure 11. DAFI (purple dots) and VNF (pink dots) gas flaring sites detections.
Figure 11. DAFI (purple dots) and VNF (pink dots) gas flaring sites detections.
Remotesensing 14 06319 g011
Figure 12. Some examples of gas flares: (a) identified by DAFI, not by VNF; (b) identified by both; (c) identified by VNF, not by DAFI (because of OF < 10%).
Figure 12. Some examples of gas flares: (a) identified by DAFI, not by VNF; (b) identified by both; (c) identified by VNF, not by DAFI (because of OF < 10%).
Remotesensing 14 06319 g012aRemotesensing 14 06319 g012b
Table 1. Indicator combination contributing to the HP value.
Table 1. Indicator combination contributing to the HP value.
ConditionIndicator Combination
I 2013 2021 ( HP SWNIR > 0   AND   HP SWIR > 0   AND   HP EP > 0 )
II 2013 2021 ( HP SWNIR > 0   AND   HP SWIR = 0   AND   HP EP > 0 )
III 2013 2021 ( HP SWNIR > 0   AND   HP SWIR > 0   AND   HP EP = 0 )
IV 2013 2021 ( HP SWNIR > 0   AND   HP SWIR = 0   AND   HP EP = 0 )  
V 2013 2021 ( HP SWNIR = 0   AND   HP SWIR > 0   AND   HP EP = 0 )
VI 2013 2021 ( HP SWNIR = 0   AND   HP SWIR = 0   AND   HP EP > 0 )
Table 2. HP conditions for each identified category (see Figure 1).
Table 2. HP conditions for each identified category (see Figure 1).
ConditionCentroids belonging to Each ConditionGas FlaresIndustrial SitesQuarries/
Mines
Nude SoilsRooftops
I32986%9%5%0%0%
II21699%1%0%0%0%
III98519%67%11%0%2%
IV13552%3%14%30%1%
V4151%56%2%0%41%
VI4100%0%0%0%0%
Table 3. Criteria of the persistence classes, predicted from the 25th (Q1), 50th (Q2) and 75th (Q3) quartiles of OFs.
Table 3. Criteria of the persistence classes, predicted from the 25th (Q1), 50th (Q2) and 75th (Q3) quartiles of OFs.
Persistence Level
LowMid–LowMid–HighHigh
OF < 15%15% ≤ OF < 20%20% ≤ OF < 30%OF ≥ 30%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Faruolo, M.; Genzano, N.; Marchese, F.; Pergola, N. A Tailored Approach for the Global Gas Flaring Investigation by Means of Daytime Satellite Imagery. Remote Sens. 2022, 14, 6319. https://doi.org/10.3390/rs14246319

AMA Style

Faruolo M, Genzano N, Marchese F, Pergola N. A Tailored Approach for the Global Gas Flaring Investigation by Means of Daytime Satellite Imagery. Remote Sensing. 2022; 14(24):6319. https://doi.org/10.3390/rs14246319

Chicago/Turabian Style

Faruolo, Mariapia, Nicola Genzano, Francesco Marchese, and Nicola Pergola. 2022. "A Tailored Approach for the Global Gas Flaring Investigation by Means of Daytime Satellite Imagery" Remote Sensing 14, no. 24: 6319. https://doi.org/10.3390/rs14246319

APA Style

Faruolo, M., Genzano, N., Marchese, F., & Pergola, N. (2022). A Tailored Approach for the Global Gas Flaring Investigation by Means of Daytime Satellite Imagery. Remote Sensing, 14(24), 6319. https://doi.org/10.3390/rs14246319

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop