Persistent Hot Spot Detection and Characterisation Using SLSTR

Gas flaring is a disposal process widely used in the oil extraction and processing industry. It consists in the burning of unwanted gas at the tip of a stack and due to its thermal characteristic and the thermal emission it is possible to observe and to quantify it from space. Spaceborne observations allows us to collect information across regions and hence to provide a base for estimation of emissions on global scale. We have successfully adapted the Visible Infrared Imaging Radiometer Suite (VIIRS) Nightfire algorithm for the detection and characterisation of persistent hot spots, including gas flares, to the Sea and Land Surface Temperature Radiometer (SLSTR) observations on-board the Sentinel-3 satellites. A hot event at temperatures typical of a gas flare will produce a local maximum in the night-time readings of the shortwave and mid-infrared (SWIR and MIR) channels of SLSTR. The SWIR band centered at 1.61 μm is closest to the expected spectral radiance maximum and serves as the primary detection band. The hot source is characterised in terms of temperature and area by fitting the sum of two Planck curves, one for the hot source and another for the background, to the radiances from all the available SWIR, MIR and thermal infra-red channels of SLSTR. The flaring radiative power is calculated from the gas flare temperature and area. Our algorithm differs from the original VIIRS Nightfire algorithm in three key aspects: (1) It uses a granule-based contextual thresholding to detect hot pixels, being independent of the number of hot sources present and their intensity. (2) It analyses entire clusters of hot source detections instead of individual pixels. This is arguably a more comprehensive use of the available information. (3) The co-registration errors between hot source clusters in the different spectral bands are calculated and corrected. This also contributes to the SLSTR instrument validation. Cross-comparisons of the new gas flare characterisation with temporally close observations by the higher resolution German FireBIRD TET-1 small satellite and with the Nightfire product based on VIIRS on-board the Suomi-NPP satellite show general agreement for an individual flaring site in Siberia and for several flaring regions around the world. Small systematic differences to VIIRS Nightfire are nevertheless apparent. Based on the hot spot characterisation, gas flares can be identified and flared gas volumes and pollutant emissions can be calculated with previously published methods.


Introduction
Gas flaring (GF) is part of the upstream oil and gas industry processes as a means of disposing of unwanted natural gas through high temperature oxidation at the tip of a stack.GF is a problem of local and global concern.GF impacts the local environment [1] through: noise [2,3], visual pollution [4,5], heat stress [4,6] and the emission of air pollutants like black carbon, polycyclic aromatic hydrocarbons, volatile organic compounds and acid rain precursors [7][8][9][10].Flaring produces greenhouse gases (GHG) and black carbon as the main by-products of the combustion.In terms of the global GHG budget gas flaring produced an estimated yearly average emission of 304 Tg CO 2 between 2003 and 2012, representing 0.6% of the global carbon dioxide equivalent anthropogenic emissions [11].Regarding the short-lived climate forcer black carbon, GF may be of regional prime importance.GF has been identified as the main input of black carbon in boreal regions [12,13], with implications for the albedo of snow-covered surfaces, the earth's radiative balance [14,15] and the Arctic amplification phenomenon [16], being therefore of global relevance [17].Again in the Arctic region, GF's contribution to the NO 2 concentrations is important, and has been increasing in the past decade [18].
The information on flared volumes and emissions is sparse and methodologically inconsistent due to technical difficulties and the lenient reporting requirements and guidelines of some jurisdictions.
Since GF creates a persistent or intermittent thermal signal, remote sensing offers the possibility of a globally consistent and independent monitoring of flaring.Important applications of GF monitoring are the estimation of the burnt gas volume in billion cubic metres (BCM) and the pollutant emissions to the atmosphere.Such estimations have already been performed with the help of a conversion factor that scales observed radiative energy release to BCM and emission factors that convert BCM into the amount of different chemical smoke species [19].The emission factors have been measured in-situ and reported in the literature [20][21][22].The conversion factor has been calculated as fraction of reported BCM and the observed radiative energy in entire countries [19], at the regional scale or for individual flares [23,24].
The identification of flaring in the night-time observations in the visual and near infra-red (Vis-NIR) spectral range [25,26] allowed for the first semi-automatic monitoring of flares from space [27][28][29].With the advent of vegetation fire detection products based on mid and thermal infrared bands (MIR and TIR), flaring was highlighted as a main source of false alarms [30,31].This feature was exploited to study flaring [32,33] and later the additional evaluation of short-wave infrared (SWIR) bands allowed for a more accurate detection [19,[34][35][36][37], which is the current state-of-the-art.
Algorithms exploiting the IR part of the spectrum for flare detection and characterisation have been developed for a number of sensors.Casadio et al. [35] considered the radiances at four wavelengths (in the SWIR, MIR and TIR bands of (A)ATSR, (Advanced) Along Track Scanning Radiometer) to be a linear combination of black body radiances from two areas with different temperatures within the satellite pixel footprint (actively flaming and background).Besides using the 1.61 µm channel for the detection, the authors also discriminated between persistent (at least 4 detections a year at a given location) and non-persisten signals, and attributed persistent detections to gas flaring.Elvidge et al. [36] developed the Nightfire algorithm, in which a Planck curve is fitted to the night-time Vis-IR measurements of the Visible Infrared Imaging Radiometer Suite instrument (VIIRS, on-board the Suomi-National Polar Partnership satellite) to retrieve the hot spot temperature.They considered an emission scaling factor to estimate the flare size.The methodology was further developed into a dual Planck curve fitting (for the background and the flare) to the observations in five bands (NIR, SWIR and MIR) of VIIRS [19].It was used to retrieve information on the global distribution and characteristics of gas flaring and the results are made public by NOAA's National Centre for Environmental Information as daily global fields (https://ngdc.noaa.gov/eog/viirs/download_viirs_fire.html).It was later demonstrated that the SWIR radiance could be used by itself to estimate the flaring radiative power [38].Gas flaring in Africa was also monitored using Landsat [37] and MODIS [23] imagery.The former algorithm used a thresholds series for the NIR, SWIR and TIR bands, and the latter used a combination of a fixed threshold and spatial filtering.The BIRD algorithm was developed to apply the bi-specral method [39] to MIR and TIR data from the Hot Spot Recognition System (HSRS) instrument on board the bi-spectral Infrared Detection (BIRD) Experimental Small Satellite (2001--04) [40][41][42].It was designed to retrieve effective temperature, effective area and effective radiative power of sub-pixel hot sources, namely fires [43,44].The methodology was later ported to the successor FireBIRD mission, with a similar sensor in IR and a modified Vis-NIR payload on board of the TET-1 (Technologieerprobungsträger-1) spacecraft [45].
The Sea and Land Surface Temperature Radiometer (SLSTR) instrument on board ESA's Sentinel-3 features night-time observations in two SWIR (S5: 1.61 µm and S6: 2.25 µm), one MIR (S7: 3.74 µm) and two TIR bands (S8: 10.85 µm and S9: 12.0 µm).The instrument also measures in two fire-dedicated bands (F1: 3.74 µm and F2:10.85 µm) with the same central wavelength and band width as S7 and S8, but extended dynamical ranges to prevent saturation over active fires.[46] The distribution of the SLSTR spectral channels with SWIR, MIR and 2 channels in TIR should allow for the detection and characterisation of GF and other hot spots via the SWIR detection and dual Planck curve fitting methodology.Compared to VIIRS, SLSTR is observing in one additional SWIR channel at night-time and it has better signal-to-noise requirements [47,48].SLSTR can therefore be expected to detect even smaller and cooler hot targets than VIIRS.However, SLSTR is missing the observations in the visible spectral range, i.e.VIIRS' DNB channel, which may make the temperature retrieval less accurate.
Furthermore, the long-term commitment of the EU Copernicus programme, which funds the Sentinel satellites, would warrant data availability well into the 2030s.
In this paper, we present a new algorithm for the detection and characterisation of persistent hot spots, including gas flares, detection and characterization from SLSTR observations and apply it to actual SLSTR data.Section 2 presents the data used.Section 3 describes the developed methodology in detail.In Section 4 we present the results of the application of the newly developed algorithm at the regional level in four regions of interest (West Africa, The North and Caspian Seas and the Persian Gulf) and evaluate our methodology against the VIIRS Nightfire product.The performance of our algorithm is further evaluated against VIIRS Nightfire and retrievals based on HSRS on-board the German small satellite TET-1 at the level of a single gas flaring site on the Yamal peninsula, Northern Siberia.Finally, in Section 5 we present our conclusions.

Data
SLSTR Level 1b version 2 products obtained from the ESA's Sentinel Expert User's Hub in January 2017 were used for this work.The downloaded products were selected using a geographic criterion (four regions of interest: West Africa, the Caspian and the North Seas and the Persian Gulf) and a time of day criteria (only night-time acquisitions).The products were sampled in the second half of 2016.For the SWIR bands (S5 and S6), the data available in the product are top of the atmosphere (TOA) radiances, while for the MIR (S7 and F1) and TIR (S8, S9 and F2) bands, the available data are brightness temperatures.The latter were converted back to TOA radiances using lookup tables provided by the European Space Agency.
FireBird TET-1 level-2 co-registered data of TET-1 night time acquisition mode (only MIR and TIR bands at 170 m spatial resolution) were obtained for near coincident or temporally close SLSTR and TET-1 observations during polar night conditions in 2016/17 over the Yamal peninsula and other areas in Northern Siberia.Due to solar light contamination, only data North of about 70 • latitude could be used, which excluded the other regions.Probably due to the extremely cold background, the standard level-2 fire processor of FireBird [45] failed to identify valid background pixels, and thus did not detect any hot clusters.The co-registered data were therefore reprocessed with the BIRD night-time algorithm [42] which was adapted to the spatial resolution of TET-1 and the cold background.The algorithm for fire detection and characterisation output includes an estimate of the fire area together with their uncertainties [42] which were used for comparison with the temporally close SLSTR retrievals.The BIRD algorithm artificially sets the upper temperature bounds to 1500 K to avoid unrealistically high bi-spectral retrievals in the case of overestimations of the TIR background radiance.This limit was adjusted to 2000 K to account for the higher temperature of gas flares within this study.VIIRS Nightfire data [36] were downloaded from the NOAA website .No further processing was applied to this dataset.
Table 1 summarizes the specifications of the used sensors, while Table 2 summarizes the data used.

Methodology
A general flowchart of the developed methodology is shown in Figure 1.The individual processing steps are described in the following sections.
The ability of SLSTR to detect hot spots in the SWIR and MIR bands was evaluated according to the instrument's requirements and is shown in Figure 3 and compared to VIIRS.In the 1.6 µm channel, SLSTR is expected to perform similarly as VIIRS.VIIRS does not record in the 2.25 µm channel at night.
In the 3.7 µm channels, VIIRS shows more sensitivity to cooler hot spots but SLSTR is, thanks to the extended dynamic range of F1 which complements S7, theoretically able to quantify larger and hotter hot spots.
Remote detection of gas flares with space-borne SWIR observations has been using a priori fixed thresholds on one or several bands [34,35,37] or determining the threshold based on the pixel's surroundings, i.e. contextual thresholding [19,23,36].In the present work, we use contextual thresholding.Due to the representation of the observations as digital numbers (DN) in the satellite as the lowest multiple of the step size that is not reported for any pixel: where } is the ordered set of the 1000 largest radiances in the granule (typically SWIR bands: 2400 × 3000 pixels per granule, MIR and TIR bands: 1200 × 1500 pixels per granule), and s f λ is the scale factor for the product and that band, which is calculated from the smallest interval between the recorded values.Radiance values below the threshold are interpreted as smoothly varying background plus "typical" noise, while radiance values above the threshold are interpreted as potential hot spot detections and are processed further.
We have tested several other approaches on the ordered set of the radiances of a product: • The cut-off value was set as the lowest radiance which produced two statistically distinct groups (using the Mann-Whitney-Wilcoxon U-test).This approach failed to produce thresholds for the majority of products.
• The cut-off value was the lowest radiance which was not included within the confidence interval of the linear regression of the previous points.This methodology lead to the inclusion of irrealistically low radiance pixels as hot pixels and the subsequent production of extremely large clusters.
• Also based on the linear regression, we examined where the slope changed significantly and whether that point could be used as a threshold.This was highly dependent on the assumptions on the linear regression and considered as non robust.
• Otsu's method [52] produced large thresholds and therefore left out real hot spots.and noise-equivalent differential temperature levels (S7 and F1) [47] and were 1.5×10 −2 , 8.4×10 −3 , 2.6×10 −4 and 2.1×10 −1 W m −2 µm −1 sr −1 for S5, S6, S7 and F1 respectively.A reflectance of 30% as given in the original works [47] was considered, together with irradiance data from SCIAMACHY http://www.iup.uni-bremen.de/UVSAT/Datasets/solar-reference-data[51], to compute the noise levels.A background temperature of 285 K is assumed.Noise and saturation for VIIRS bands were computed using the signal-to-noise ratios given by Cao et al. [48].Brightness temperature to radiance conversions were all computed using Planck's law.Temperatures between 600 K and 1500 K represent, roughly, the lower and upper bounds for biomass burning, volcanoes and industrial sources.The range 1300-2000 K is typical for gas flaring temperatures.[  The threshold is set contextually as the lowest radiance value whose difference to the closest inferior radiance value is larger than the scale factor.
• An iterative clustering method which minimized differences between radiances of data points within two clusters also lead to large thresholds.
• Based on the fact that the night-time background is rather homogenous, the variance of the radiances of the product should be increased strongly by hot spots, while the background pixels should have a small variance.This approach was very sensitive on the a priori explained variance expected and was therefore abandoned.
The selected approach was the most robust, identifying hot pixels in different contexts with low variation between different products.This approach captures hot pixels when there is a gradual or an abrupt increase (upper and lower plots of Figure 4, respectively).While not fixed, and therefore adaptable to product-specific conditions, this contextual methodology is not statistical and therefore not directly sensitive to the number of hot events sampled and their intensity.
At the end of this step, a collection of i hot pixels are registered, for each of the four SWIR and MIR bands, with the following information: • hot pixel location (x i , y i ) as pixel index pair (the the pixel length on the x and y axis is computed as average of the ground distance between the pixel and its direct neighbours in those axes).Peer-reviewed version available at Remote Sens. 2018, 10, 1118; doi:10.3390/rs10071118

Clustering
The clustering of contiguous hot pixels is necessary because the signal of a single hot spot may influence more than one pixel [33,53].It may also be the case that a single hot spot is detected in two or more adjacent pixels.[23,37] Indeed, in the case of flaring a facility may comprise arrays of many flares with sizes of up to several hundred metres, while individual stations are several kilometres apart.
The misregistration between SLSTR channels is not an issue here because the IR bands are treated independently up to the step right before the characterisation (see section 3.6).
A cluster is defined as a set of hot pixels which are adjacent spacially.A cluster may comprise only one pixel if none of its adjacent pixels are hot.
A background area is defined as extending two pixels beyond the cluster's limits in any direction, including diagonally and excluding the cluster of hot pixel(s) itself .For each single pixel i of the n cluster's pixels, the radiance B λ,i and area A λ,i are registered.Likewise, for each single pixel j of the m background pixels, B λ,j and A λ,j are registered.
For each cluster, the following quantities are then computed: the average radiance and its standard deviation in [Wm −2 µm −1 sr −1 ], the average background radiance and its standard deviation in [Wm −2 µm −1 sr −1 ], the cluster area in [m 2 ], its x-position as across-track index and y-position as along-track index, and the number of cloud-flagged pixels (using the standard SLSTR cloud product) for both the cluster and the background.

Misregistration characterisation
A hot event at temperatures typical of a gas flare, a steel mill or an iron smelter, will produce a local maximum in the SWIR and MIR channels (see Figure 2), which will be translated into a cluster in those bands.The MIR fire channels S7/F1 of SLSTR have been designed with a different footprint size and position than the corresponding SWIR channels.Because of these differences the clusters in the different bands cannot be spatially superimposed.
We determine the misregistration of channels S6, S7 and F1 relative to channel S5 by assuming that the majority of hot spots are isolated point sources: For each hot spot cluster in S6, S7 and F1, taken from a random set of 500 SLSTR scenes, the closest S5 cluster is found.For each pair, the distances of the cluster centres in along-track and across-track directions are calculated in coordinates of product image grid indeces.These distances are plotted over the across-track pixel position in Figure 5.
A systematic dependency on the across-track position is obvious in several of the plots.We, therefore, fit parabolas to the data.They represent our best estimate of the misregistration between S5 and S6 resp.S7 resp.F1.The fitted parameters are given in Figure 5.In order to determine confidence intervals for the misregistration estimates, the fitted parabolas are shifted vertically such that 10% of the data points are below resp.above.

Misregistration correction: building multi-band clusters
Multi-band clusters C hotspot are subsequently constructed.They consist of the single-band clusters C λ with λ in S5, S6, S7 and F1 that observe the same hot source.The SWIR cluster C S5 is used as reference.Then those clusters from the remaining bands (C λ with λ in S6, S7 and F1) that are closest to after correction of the misregistration are added.In doing so, only clusters in the confidence intervals are considered: where C λ is any of the λ-band clusters (λ in S6, S7 and F1), is the λ-band cluster within the multi-band cluster C hotspot (λ in S6, S7 and F1), For the TIR bands S8, S9 and F2, the flaring high-temperature event is not expected to impact the radiance (Figure 2).The average radiance for each TIR band 2 pixels around the reference S5 cluster position is then associated to the multi-band cluster.

Radiance corrections
The SWIR radiances in S5 and S6 exhibit a systematic overestimation of 11 and 20%, respectively.
[54] Following recommendations by ESA, all values in the SWIR bands were corrected accordingly.
This correction needs to be verified and possibly updated for future versions of SLSTR products.
After building the multi-band clusters, the MIR S7 channel is checked for saturation.If the S7 radiance is above the upper limit of the linearity range (0.56 Wm −2 µm −1 sr −1 , corresponding to a brightness temperature of 306K) in any of the hot pixels within the cluster, then it is discarded.If the S7 cluster is discarded or not present, the F1 cluster is used.The F1 band is a fire dedicated channel of SLSTR which measures at the same wavelength as the S7 channel, but with a larger dynamic range and lower sensitivity.Only F1 clusters where all of the pixels have a brightness temperature within the range 300-480 K are considered, as this is the range where the Level 1 quantification is accurate.

Super Cluster definition
The areas of the individual clusters in the SWIR and MIR bands within a multi-band cluster may vary.To overcome this, hypothetical super clusters, with perfect coregistration and identical footprint areas in all bands, are built.The largest of the individual band clusters' areas A λ is chosen as area A cluster of the super cluster (Figure 6).
For each band's cluster C hotspot λ within a multi-band cluster C hotspot , the observed cluster radiance B obs λ is calculated as area-weighted average of the observed radiances of the cluster and the background: It is assumed that the background radiance B bg λ is constant in the vicinity of the cluster.For the TIR bands (S8, S9 and F2) the weighting is not necessary as the hot spot signal is not expected to be distinguishable from the background in that spectral region and the average radiance of background and cluster, is used.

Planck curve fitting
In order to determine the temperature and the area of the flaring event, the sum of two Planck curves is fitted to the radiance data of the multiband cluster, as established by Elvidge et al. [19].The two Planck curves represent the two contributors for the IR radiance measured by the sensor at night, i.e. the flaring event and the background, each weighted by its respective relative area:

Radiative power
The radiative power (RP) of the hot source may then be computed using the Stefan-Boltzmann equation and assuming a perfect emissivity (in W): where σ SB is the Stefan-Boltzmann constant (5.670373 × 10 −8 Wm −2 K −4 ).The respective uncertainty is computed by propagation of errors.
Elvidge et al. [19] used the VIIRS M10 band as a primary band for the detection of gas flares and found that "there is a strongly coherent linear relationship between top of the atmosphere and

Regional study: 4 flaring regions
The described algorithm was applied to observations in the four regions of West Africa, the North sea, the Caspian sea and the Persian gulf, cf.Table 2.

Detection thresholds
Figure 8 shows the relationship between the retrieved temperature and area of hot spots detected by SLSTR and compares it with the noise-equivalent radiance levels for the detection bands (S5, S6, S7 and F1) and with the limit of detection derived for VIIRS Nightfire [36].The cloud of points generally follows the M10 detection limit derived by Elvidge et al. [36], but also expands slightly to lower threshold values in the regions of typical industrial temperatures (around 1000K) and gas flaring (1500-2000K).
Hot spots detected in the S5 channel only (grey crosses in Figure 8) represent 5% of all the hot spots within the study regions.They exhibit temperatures in the lowest end (600-1000K) of the range of retrieved temperatures and areas in the upper end (10 2 -10 5 ) of the range of retrieved areas.These values are highly uncertain because only one constraint exists in the SWIR and MIR range so that the TIR radiances are expected to be overfitted.The retrievals are thus treated as "low-accuracy" (see Section 4.1.4).Clusters detected in the MIR (green and blue crosses in Figure 8) represent 28% of all the hot spots within the study regions (3% S5 and MIR and 25% in S5, S6 and MIR).The lower limit of the cloud follows an isoradiance line that represents a larger radiance value (approximately 5 W m −2 sr −1 µm −1 ) than the noise-equivalent for S7 or F1 from the instruments specifications.We attribute this to the larger variability of the background in the MIR.Hot spots detected in both SWIR channels (S5 and S6) but not in the MIR channels (pink points in Figure 8) make up the most of the detections (67%).Hot spots in the region below the M10 limit of detection (blue dashed line) as derived by Elvidge et al. [36] were detected in both SWIR channels.They are mostly in the expected range for gas flares (1300-2000K [19]).This suggest some capability, thanks to the availability of both SWIR channels, to detect smaller gas flares than those detected by VIIRS.SLSTR can also characterise the smaller gas flares more effectively using the Planck curve fitting method.

Temperature and area retrievals
Figure 9 shows the hot spot detections of SLSTR with temperatures in the range [500 K, 5000 K], which we consider to relate to actual hot spots on the ground .The methodology identifies hot spots in regions where flaring is known to occur (e.g. the North Sea oil fields, the mouths of the Congo and Niger rivers, the Persian Gulf and the Tigris and Euphrates rivers valleys).Observing fire temperature on a global scale, even only at night-time, is complementary to the widely used burnt area and fire radiative power observations and has great potential for reducing errors in the current vegetation fire emission datasets.
Hot spot characteristics evaluated from space could be dependent on the viewing angle.This was shown to be relevant for VIIRS [19].However, SLSTR has a narrower range of viewing angles and the effect is not expected to be as pronounced.Since we fit an unconstrained Planck curve, this effect would be noticed in both the retrieved temperature and area.We therefore present the retrieved hot spot radiative power variation dependending on the across track index in Figure 11.There is no apparent effect.

Persistence
Hot or bright events that can be observed from space in the IR part of the spectrum at night include gas flares, wildfires, auroras, industry (e.g.steel mills) and volcanoes.[34,36] Noise and an increased exposition of the sensor to radiation due to the South Atlantic Anomaly may also generate spurious hot spots.[19,34] It is possible to discriminate between these events based on retrieved characteristics, such as the temperature.[19,55] In the present work, we opt for an analysis of the persistence of the signal at the location of a given hot spot as it has been used previously to filter out noise and ephemeral phenomena [19,34].
We base the definition of our criterion on the work of Casadio et al. [34], who used a threshold of at least 4 detections a year for (A)ATSR, the heritage instrument of SLSTR.Since the swath of (A)ATSR is roughly one third of SLSTR's, the criterion would correspond to 12 times a year for SLSTR.The sampling period is roughly 2 months, corresponding to 2 observations within that period.This was raised to 3 in order to filter out noise.
Although not suitable for a thorough analysis due to the short sampling period, the persistent hot spots (at least 3 detections within a spatial accuracy of ±0.02 degrees in longitude and latitude over the sampling period) mainly correspond to locations within known flaring regions (see Figure 12).In future works where this methodology will be used on data from longer periods, this threshold might need to be raised or otherwise adapted.For known interferences, persistent or semi-persistent hot sources such as volcanoes, the use of a mask will also be useful.
The persistent locations thus derived were compared with high resolution imagery (Google Earth, for the Caspian Sea test region) in order to check for the existence of a hot spot: at 96 (65%) of the 148 persistent hot spot locations a gas flare was visible, 34 (23%) locations were offshore or onshore but lacking sufficient resolution, 12 (8%) were onshore locations without an industrial area nor a visible gas flare nearby and 6 (4%) were at an industrial site without a visible flare.

Selection of persistent hot spots and radiative power computations
In order to filter low-quality determinations, cloud cover and overfitting were considered.
Clouds may interfere with the amount of radiation measured by the SLSTR instrument, in turn interfering with the temperature, area and radiative power retrievals.However, simply discarding cloudy pixels could incur a large omission error since pixels containing gas flares are frequently marked as clouds, namely as isolated cloudy patches (e.g. for the VIIRS instrument [19] ).For this reason, we   analyse the cloudiness of the background, as defined in section 3.2, and discard only observations with less than 3 cloud-free background pixels.
In order to avoid too much dependence on the TIR channels, which are not affected by small hot sources with the temperature of gas flares, we discard persistent hot spots for which no clusters were detected in S6, S7 and F1.In the work of Elvidge et al. [19], signals detected in the 1.6 µm channel only are not quantitatively evaluated and information from nearby hot spots is used.Although we compute the hot spot temperature and fraction area and the background temperature with the dual Planck curve fitting (section 3.7) using data from the super cluster (S5 and TIR channel, section 3.6), we opt to discard those results (clustered in a low temperature/large area region, Figure 8) and to not assume any characteristics for that particular detection.Since these detections are only 5% of the total hot spot persistent detections, we assume that they do not significantly influence hot spots characterization within the study regions.
Figure 13 shows the distributions of the retrieved temperature, area and radiative power for detections that have been filtered with the persistence, cloud cover and S6/S7/F7 availability criteria described above.We will label these "high-accuracy persistent" below.A range of temperatures around 1800 K can be considered as characteristic for flaring [19,36], and the distributions of the high-accuracy persistent hot spots fall mostly within the expected range for gas flares, but observations in the 500-1500 K temperature range are also important.The retrieved temperatures show a clearly unimodal distribution approximately centered at 1600 K for the Persian Gulf and the Caspian Sea test regions.For the North Sea test region, the distribution is bi-modal, with modes around 1000 K and 1600 K.While the latter is probably associated with gas flaring, the former is more likely associated with industry (see Figure 12, the colder persistent detections are onshore.).[19] There were very few high-accuracy persistent hot spots for the West Africa test region due to unfavourable cloud conditions.

Comparison with VIIRS Nightfire
The SLSTR hot spot at persistent locations were compared to the VIIRS Nightfire data for a same time period and a roughly similar area.The VIIRS Nightfire results were subject to the same spatial persistency analysis as the SLSTR data (3 times within a spatial accuracy of 0.02 degrees in longitude and latitude).VIIRS detects more gas flaring locations than SLSTR, which can be traced back to its wider swath (3040 km against 1420 km), and thus a shorter revisit time.Another reason for the larger number of detections by VIIRS is that the Nightfire algorithm processes single pixels, selected as the local maxima of pixels above the threshold, the method here aggregates contiguous hot pixels into clusters.Despite these differences , all but one persistent locations detected by SLSTR were also detected by VIIRS (Table 3).The VIIRS Nightfire algorithm produces a slightly larger fraction of detections in the 1.6 µm channel only (7% for Nightfire against 5% for the present work).But when considering Nightfire detections in the 1.6 µm and visible (0.5-0.9 µm, Day-Night Band, DNB) channels, the fraction is very important: 55%, as is the fraction of SLSTR hot spots detected in both SWIR channels only (67%).The MIR channels, both of VIIRS and SLSTR, produce less detections than the SWIR channels.We interprete this as a larger variability of the background in the MIR region, as well as a lower sensitivity to high temperature sources (Figure 2).The clustering of contiguous hot pixels used in the present method also explains why the VIIRS Nightfire temperature retrievals tend to peak at higher values than the SLSTR ones (Figure 14).
Selecting only pixels which are local maxima, as the VIIRS Nightfire algorithm does, will produce detections with higher radiances and mostly higher temperatures will be retrieved.

Single site study: Bovanenkovo, Yamal peninsula
As another test case, a flaring site in the Bovanenkovo field was observed in more detail.This analysis is performed in order to test the capabilities of the presented algorithm, not only in statistical terms (as conducted when comparing with VIIRS Nightfire), but also at the single site level.The Bovanenkovo field, located on the western shore of the Yamal peninsula in northern Siberia, Russia, produces natural gas and condensates.The flaring location was repeatedly observed by the sensors SLSTR, HSRS and VIIRS over a sampling period of 19 days (between 15/12/2016 and 2/12/2017), with very high temporal co-location between SLSTR and HSRS.The detections are clustered around 4 distinct locations (see Figure 15) and are summarized in Table 4 Due to clustering in the SLSTR and HSRS data analysis, these sensors may detect one cluster in a situation, in which the VIIRS product lists two or more flares in very close proximity.For SLSTR, the observing conditions were very close to twilight (the night-only mask was not considered in order to have the opportunity to observe the same flaring site by the three sensors), and thus not ideal.
Figure 16 shows the temperature, area and radiative power determinations at Locations 2 and 3 for the 3 instruments.
At Location 2, the single SLSTR observation agrees well with the determinations based on VIIRS and HSRS for the temperature.At Location 3, the temperature determinations based on SLSTR are lower than those based on VIIRS but in the same range as those based on HSRS.The gas flare temperature depends on the flared gas composition and the completeness of combustion, which depends on flare design and meteorological conditions [56][57][58][59].Thus the flaring characteristics may vary in time.The methodological differences may also explain lower temperatures in the SLSTR product as explained above.
While the retrieved areas from TET-1 are larger than those from SLSTR, area determinations based on VIIRS and SLSTR are close for both locations except the highly uncertain SLSTR observation on 2/1/2017.The shape of the flame, and thus its area as observed by a satellite, is dependent on the flow rate and the meteorological conditions [57,60] and therefore also variable in time as its temperature.
Despite the variability seen between sensors due to inherent variability of flares and methodological differences, the radiative power determinations, the variable to which the emissions computations are linked, compare well between sensors.

Conclusions
We have adapted the VIIRS Nightfire algorithm [19] for the detection and characterisation of persistent hot spots using the SLSTR instrument on-board the Copernicus Sentinel-3 satellite [46].The

Figure 1 .Figure 2 .
Figure 1.General algorithm flowchart for the detection of persistent hot spots and their characterisation: temperature (T), area (A) and radiative power (RP).The dotted line represents the one-time misregistration determination.The parameterization was then used when building the multi-band cluster.n cloud− f ree bg> 3 represents at least three cloud-free pixels in the background.> 4λ means there is radiance data at least 4 of the 5 wavelengths considered by the algorithm.We do not believe spurious signals are hot sources, but noise such as e.g.lunar reflection or the Southern Atlantic anomaly.Transient hot spots are e.g.vegetaton fires.Persistent hot spots are e.g.gas flares, volcanoes, industry.We believe low accuracy persistent hot spots are hot sources, but the fitted parameters (T and A), and therefore RP, are relatively uncertain.

PreprintsFigure 3 .
Figure 3. Isoradiance lines for the noise (SWIR bands S5 and S6 and MIR bands S7 and F1) and saturation levels (MIR bands S7 and F1) for SLSTR and noise and saturation for the SWIR and MIR bands (M10, M12 and M13) of VIIRS used in the Nightfire algorithm.The SLSTR noise levels were computed using the End-of-Life noise-equivalent differential reflectance (S5 and S6) and noise-equivalent differential temperature levels (S7 and F1)[47] and were 1.5×10 −2 , 8.4×10 −3 , 2.6×10 −4 and 2.1×10 −1 W m −2 µm −1 sr −1 for S5, S6, S7 and F1 respectively.A reflectance of 30% as given in the original works[47] was considered, together with irradiance data from SCIAMACHY http://www.iup.uni-bremen.de/UVSAT/Datasets/solar-reference-data[51],to compute the noise levels.A background temperature of 285 K is assumed.Noise and saturation for VIIRS bands were computed using the signal-to-noise ratios given by Cao et al.[48].Brightness temperature to radiance conversions were all computed using Planck's law.Temperatures between 600 K and 1500 K represent, roughly, the lower and upper bounds for biomass burning, volcanoes and industrial sources.The range 1300-2000 K is typical for gas flaring temperatures.[19]

Figure 4 .
Figure 4. Four examples of thresholding for the S5 band.The horizontal line marks the threshold value.The threshold is set contextually as the lowest radiance value whose difference to the closest inferior radiance value is larger than the scale factor.

1 2Figure 5 .
Figure 5. Parameterizations of the distances between a cluster in S5 and the closest cluster in S6, S7 and F1 as a function of the scan index.The solid line shows the best fit (second order polynomial).The dashed lines are parallel to the best fit so that 10% of the points are above resp.below.They represent the confidence interval.

Figure 7 .
Figure 7.An example of the temperature and area retrieval by dual Planck curve fitting.
atmospherically corrected RH [radiant heat]".The authors thus conducted their work with the uncorrected TOA radiances.Since the S5 band of SLSTR is within the same clear atmospheric window, the same choice was made for the present work.Preprints (www.preprints.org)| NOT PEER-REVIEWED | Posted: 24 June 2018 doi:10.20944/preprints201805.0020.v2Peer-reviewed version available at Remote Sens. 2018, 10, 1118; doi:10.3390/rs100711184. Results Figure 10 shows the temperature distributions in nine specific areas: four areas with know gas flaring, one area with biomass burning and four areas without any known hot objects.Hot spots detected offshore where no oil extraction activities are present show a clearly lower temperatures than for the other study regions.Hot spots were also detected in the plateau of central Angola and the Southern D. R. Congo.The intermediate temperatures retrieved in those regions indicate the presence of biomass burning.

Figure 8 .Figure 9 .Figure 10 .
Figure 8. Relationship between the temperatures and areas of all retrieved hot spots.The solid lines represent the noise-equivalent radiance levels for SLSTR SWIR and MIR channels computed from the instruments specifications (see Figure3).The dashed line represents the limit of detection for the VIIRS M10 channel derived by Elvidge et al.[36].

PreprintsFigure 11 .Figure 12 .
Figure 11.Retrieved radiative power for hot spots as a function of the across track index (ATI).

Figure 13 .
Figure13.Violin plots for the temperature, area and the radiative power retrieval distributions of the high-accuracy persistent hot spots.The number of high-accuracy persistent hot spots is also indicated.

Figure 14 .
Figure 14.Comparison between temperature, area and radiative power retrieved by SLSTR and VIIRS at persistent locations.

Table 1 .
Comparison between the used sensors.

Table 2 .
Data used in this study.The number of VIIRS products was derived from the VIIRS Nightfire dataset from NOAA, corresponding to the products which exhibited detections co-located in space and time.

Table 3 .
Summary of the detections based on SLSTR and spatial comparison with VIIRS Nightfire.

Table 4 .
Summary of the detections at the Bovanenkovo gas and condensates field in the Yamal Peninsula, Siberia.