1. Introduction
Resources invested in enhancing wastewater treatment plant (WWTP) performance have long been focused on improving the effluent quality [
1,
2,
3,
4]. In recent years, new challenges have been addressed to improve the sustainability of WWTP operation. Energy consumption and greenhouse gas (GHG) emissions are the two critical factors that an increasing number of water utilities are beginning to evaluate when planning for major capital upgrades and operational modifications [
5,
6]. Recent studies have identified wastewater treatment plants as important emission sources of anthropogenic GHG, contributing to climate change and air pollution [
7,
8]. Studies indicate that methane (CH
4) emissions from municipal and industrial wastewater management account for approximately 5% of global methane production [
9]. In addition, WWTPs generate and release carbon dioxide (CO
2) and nitrous oxide (N
2O) from the treatment processes. CO
2 is also emitted during the combustion of the fossil fuels necessary for plant operation [
10,
11]. All of these contribute to the fugitive GHG emissions from WWTPs.
Methane is a commonly reported GHG from WWT, with a global warming potential (GWP) 25 times higher than CO
2 [
12]. At a typical WWTP, Daelman et al. [
13] have reported that around one percent of the chemical oxygen demand (COD) that was brought into the WWTP was released as methane. CH
4 emissions have been identified from various process units in WWTPs, even under toxic conditions [
14,
15]. Numerous studies have reported that the primary sources of methane are associated with the sludge process, including sludge holding, thickening, dewatering, anaerobic digestion, and the exhaust line from cogeneration [
16,
17]. The sludge process accounts for approximately 72% of methane emissions from wastewater treatment, while the remaining emissions are associated with the wastewater processes at locations where anaerobic conditions are developed and relatively higher organic loads are present [
18].
Identifying and quantifying CH
4 emissions is essential to the establishment of sustainable wastewater treatment practices for mitigating and reducing GHG emissions. Traditional emission factor-based estimation has been widely criticized for its lack of accuracy [
19]. Recently, short-term, sporadic measurements for CH
4 have been attempted in order to quantify the total emission from the entire WWTP [
20]. However, at a typical WWTP, factors like weather, influent flows, pollutants loads, and operational conditions are always changing, causing the methane concentrations to fluctuate. As a result, sporadic measurements that only catch isolated “snapshots” tend to underestimate the temporal fluctuations and spatial heterogeneity of methane emissions from the WWTPs [
21]. Consequently, extrapolating snapshot measurements for estimating longer-term, e.g., one year, emissions are expected to result in misleading results, which eventually will create challenges for evaluating the effectiveness of mitigation strategies.
Therefore, developing an approach for long-term monitoring provides the advantage of capturing the temporal and spatial variations of methane emissions at WWTPs. Satellite imagery has been an important tool for identifying and characterizing methane emitting facilities, such as oil and gas well pads and surface coal mines [
22,
23,
24,
25,
26].
Satellite technology allows for the mapping of emission sources within a site, offering the potential to identify localized hotspots that may be disproportionately responsible for the overall emissions. Satellite data also allow the analysis of methane variations on different temporal scales, providing a more comprehensive understanding of the emission dynamics, which is critical for a detailed quantification, identification, and trending of the emissions, upon the basis of which mitigation strategies are then developed. In addition, utilizing high-frequency monitoring data means that timely detection of, and response to, unexpected emission events can be achieved [
27]. Without the need for human assistance, automated detection systems powered by deep learning algorithms can identify methane emissions by analyzing multispectral satellite imagery, greatly increasing the efficiency and accuracy of the analysis [
28]. Due to the various benefits of satellite imagery, this approach has been applied to the estimation of methane emissions from a range of sources, including agricultural practices such as paddy fields, landfills, and wetlands [
29,
30,
31,
32]. Although the GWP of N
2O is higher than that of CH
4, the use of satellite imagery to monitor nitrous oxide emissions remains a developing area of study. Nitrous oxide is a notable greenhouse gas, but its remote monitoring is complicated by its low atmospheric concentration and spectral characteristics that make it difficult to distinguish from other gases in satellite observations.
Various methods for N
2O detection are under investigation, primarily using satellites equipped with hyperspectral or thermal infrared sensors capable of capturing the specific wavelengths absorbed by N
2O. Satellites like Sentinel-5P and missions such as GOSAT have the spectral capability to detect and quantify trace gases, including N
2O; however, they are limited by their spatial and temporal resolution [
33,
34]. The scarcity of data for N
2O is due to its lower atmospheric concentration and weaker absorption in critical spectral bands. Additionally, N
2O can interfere with other gases, complicating the analysis of satellite imagery. These limitations reduce the effectiveness of satellite imagery in investigating smaller emission sources, such as wastewater treatment plants.
The tropospheric monitoring instrument (TROPOMI) aboard the Sentinel-5 precursor satellite provides daily global methane measurements with up to 5.5 × 7 km
2 pixel resolution [
35]. Normally, satellite images of this resolution are considered too coarse to resolve smaller sources that are spatially clustered and produce a plume extending shorter than 1 km in distance [
36]. The GHGSat microsatellite instruments are specifically designed to detect methane sources using fine pixel resolution (25–50 m) over limited domains (12 × 12 km
2) and with relatively higher precision (∼1–15%) [
37]. Hyperspectral imaging spectrometers designed to observe land surfaces at 1–10 nanometers spectral resolution with 30 m pixel resolution can detect large methane plumes, as was recently demonstrated with the Italian Space Agency’s PRISMA instrument [
38,
39]. In addition, the frequency of revisiting a specific location with targeting instruments is constrained by multiple factors, including the satellite’s spatial coverage, operational limitations regarding tasking, and the finite number of satellites available for observation. Each satellite is limited to observing a specific area during its overhead passage and, due to competing tasking priorities, imaging cannot be exclusively assigned to a single location. Furthermore, the limited number of satellites in orbit complicates the achievement of regular monitoring. A constellation of satellites is essential to increased revisit frequency, because of the improved coverage and more continuous observations across various locations.
To date, the use of satellite imagery for analyzing methane emissions from wastewater treatment has been extremely rare, if not nonexistent, for which the likely reason is that the low pixel resolution of earlier satellite images is not able to capture the complex, heterogeneous distribution of the GHG within a relatively smaller footprint. However, as previously described, recent advances in satellite technology have made significant progress meaning that high-resolution tracking capability has become available for methane monitoring. This allows for the development of image processing algorithms to capture the methane plumes from facilities of smaller footprints such as WWTPs [
40].
In this study, an approach has been developed, based on Sentinel-2 images, for characterizing 5-year fugitive CH4 emissions from a wastewater treatment plant located in Canada. The analysis provides, for the first time, a comprehensive and overall mapping of fugitive methane emissions at the WWTP. The results reveal both the time-varying nature and the spatial heterogeneity of fugitive methane emissions from the wastewater treatment, indicating the importance of continuous measurements for quantifying the annual emissions for compliance reporting.
2. Literature Review
The detection and quantification of methane emissions using satellite imagery have significantly improved in recent years, driven by the need for more accurate and higher-resolution data to inform climate policy and guide mitigation efforts. Various satellite technologies and methodologies have been investigated in different studies.
Satellite-based methane detection typically relies on identifying the absorption in the short-wave infrared (SWIR) region of backscattered sunlight, specifically in the spectral regions where methane is known to absorb light. Hyperspectral satellites, with their high spectral resolution in the SWIR bands, enable a more precise determination of methane column concentrations [
24,
41,
42]. These satellite images have either smaller spatial coverage but high spatial resolution (e.g., target-mode satellites, such as PRISMA) or low spatial resolution but broader coverage (e.g., Sentinel-5). The use of multiple Sentinel satellites in tandem for detecting and quantifying methane leaks has been covered in a recent study by Pandey et al., in which Sentinel-3 was found to enhance Sentinel-5p’s coverage to detect methane leaks of up to 8–20 tonnes per hour [
30]. This study demonstrates the value of using different satellite instruments synergistically to achieve a more comprehensive characterization of methane emissions.
The study by Wecht et al. (2014) utilized a satellite-based approach for quantifying methane emissions across large regions. This research estimated that the total methane emission from U.S. livestock is approximately 40% greater than was previously estimated in the emission inventories. However, a similar discrepancy for U.S. oil and gas emissions was not observed [
40]. Varon et al. (2021) have explored how Sentinel-2 data could be used to locate and quantify significant methane point sources. Using the reflectance differences of Sentinel-2 bands 11 and 12 with three retrieval methods, this research monitored methane point sources in the Hassi Messaoud Oil Field in Algeria and the Korpezhe Oil and Gas Field in Turkmenistan over several years [
27].
Most recently, deep learning techniques have been integrated in to the analysis of satellite-based methane detection. Rouet-Leduc and Hulbert (2024), introduced a deep learning-based tool that overcomes the limitation of multispectral satellite data, thus enabling high-resolution methane detection. This model, once validated against airborne measurements, was found to be capable of detecting a methane point source as small as 0.01 km
2, emitting approximately 200–300 kg of CH
4 per hour. This approach allows for the automatic detection of methane emissions at a higher resolution, and thus improving overall monitoring capability [
28]. Radman et al.’s (2023) deep learning approach involved the creation of the EfficientNet-V2L model, which was optimized to estimate the methane emission rates from Sentinel-2 data. With a Pearson R of 0.957, this approach achieved a much higher accuracy, outperforming conventional physical-based techniques such as the integrated mass enhancement (IME) method [
23].
To address the limitation of hyperspecial satellites in observing small scale sources, there has been increasing interest in developing methane detection techniques using data from multispectral satellites like ESA’s Sentinel-2 constellation. Sentinel-2 provides high spatial resolution and global coverage, scanning the entire Earth every 2 to 5 days. However, multispectral satellites offer high spatial and temporal resolution at the expense of much less spectral information, with only a dozen spectral bands compared with the hundreds available on hyperspectral satellites [
26,
27].
Despite this, a number of studies have been recently performed using satellite imagery for methane emission quantification from industrial facilities and super emitters. Currently, there is no published research that was specifically focused on using satellite imagery to investigate methane emissions sources for smaller footprint facilities like WWTPs. Methane emission-related studies have been focused mostly on the oil and gas sector, landfills, and farming activities. Because of their higher emission rates and much larger footprint, the plumes are easier to detect using satellite imaging. In regard of WWTPs, the fugitive nature of the emissions generally creates lower intensity plumes, hence, detection and quantification can be challenging using the existing satellite data. In addition, most of the existing satellite sensors used for methane detection are incapable of achieving a spatial and temporal resolution that is high enough to capture and localize the highly heterogeneous methane emissions from wastewater treatment. Satellite remote sensing of methane provides valuable data over extensive areas by detecting specific infrared wavelengths that are absorbed by atmospheric methane. However, these measurements typically capture column-integrated concentrations from the surface to the upper atmosphere, which may not accurately represent ground-level values. Discrepancies arise from factors such as vertical methane distribution, atmospheric scattering, cloud cover, instrument calibration, and the reflectance characteristics of various surfaces. Additionally, satellites collect data at predetermined intervals and resolutions may not adequately capture small, localized sources. To mitigate these uncertainties, we focus on low-cloud days and emphasize the importance of validating satellite observations with in situ ground data to ensure precise correlation with ground-level concentrations.
This study is the first evaluation of methane emissions from WWTPs utilizing the latest satellite technologies. With the recent development and launch of advanced satellites, higher-resolution sensors have been deployed on the Sentinel-2 constellation and other specialized missions like GHGSat. The Sentinel-2 satellites were selected in this study. Sentinel-2 was originally designed for environmental risk management, land cover classification, land change detection, and terrestrial mapping. The high pixel resolution of Sentinel-2 images allows for detailed investigation of fugitive methane emissions within complex environments with smaller footprints, such as those found at WWTPs. Further- more, Sentinel-2’s frequent revisit interval provides a larger number of images, offering greater insight into the temporal variations of CH4 emissions than other satellite products.
3. Materials and Methods
3.1. Location
In this project, a WWTP in Canada was selected as the study site. It is one of the largest wastewater treatment plants in Canada.
The WWTP shown in
Figure 1 operates with an average capacity of 396 million liters per day (ML/d). It employs biological nutrient removal (BNR) technology, making it not only one of the largest BNR plants in the country but also the largest cold-weather BNR plant globally. The plant must function under challenging environmental conditions, including winter air temperatures that can drop as low as −40 °C and influent temperatures that reach a minimum of 10 °C. These conditions provide a unique setting for analyzing methane emission under cold-weather operational conditions.
3.2. Satellite Image Data and Characteristics
The data for this research project were obtained from Sentinel-2 satellites. Sentinel-2 is a part of the Copernicus Programme and is an Earth observation mission that captures high-resolution optical imagery of land and coastal waters. It is operated by the European Space Agency, with the satellites built by a consortium in Germany [
43,
44]. The Sentinel2 constellation is composed of two satellites (S2A and S2B), positioned 180° apart in a shared sun-synchronous orbit, with a local solar time of 10:30 AM at the equator during the descending pass. Sentinel-2A was launched in June 2015, followed by Sentinel-2B in March 2017. Both satellites are equipped with a multispectral instrument (MSI) that continuously captures the Earth’s surface in 13 spectral bands, covering wavelengths from the visible to the shortwave infrared (SWIR) regions. The pixel resolution ranges from 10 to 60 m, spanning a 290 km cross-track swath [
43]. An overview of the available bands of Sentinel-2 is available in
Table 1.
The twin-satellite configuration of Sentinel-2 enables complete global coverage every 5 days, with revisit intervals of 2–3 days at the mid-latitudes. This study demonstrates that Sentinel-2’s SWIR bands 11 (~1560–1660 nm) and 12 (~2090–2290 nm), with a 20 m pixel resolution, are effective for detecting methane plumes from large point sources and quantifying their emission rates. These bands capture radiance within methane’s absorption features at 1650 nm and 2300 nm. Although Band 12 overlaps with the stronger and broader absorption feature at 2300 nm and is more sensitive to methane than band 11, both bands have comparable signal-to-noise ratios (SNRs) at their respective reference radiances, despite band 12 receiving lower solar irradiance [
45]. Band 11 can thus serve as a proxy for the continuum, given its spectral proximity to band 12 and generally similar surface reflectance. While the multispectral instrument has a spectral resolution of approximately 100–200 nm, which is too coarse for standard hyperspectral retrieval of methane column concentrations in the SWIR, methane columns can still be estimated by analyzing reflectance differences between the spectral bands and across multiple satellite passes. To use the data sets downloaded for 5 years, preprocessing was undertaken using the semi-automatic classification (SCP) plug-in, with the default values, in QGIS software (3.28.3-Firenze), and the built-in python (version 3.9.5) console in the software [
45]. The next stage of processing in this project involves applying principles from optical spectroscopy.
3.3. Methane Column Concentration Retrievals
Beer’s law is a fundamental principle in spectroscopic quantitative analysis, establishing a direct relationship between the concentration of a sample and the intensity of radiation measured by a spectrometer. This principle is essential for quantifying the concentration of a substance based on its interaction with light.
The intensity of the incident radiation was denoted as I
0, referring to band 11 reflections, while the intensity of the transmitted radiation was denoted as I, which in this study was assumed to be band 12 reflections. Absorbance (A) was calculated as the logarithm of the ratio of these intensities.
To relate absorbance to the concentration of the sample, Beer’s Law—also known as the Beer–Lambert Law or Bouguer–Lambert–Beer Law, was applied. In this context,
The molar absorption coefficient (ϵ), also known as the molar absorptivity, measures how strongly a chemical species absorbs light at a specific wavelength. This coefficient is an intrinsic property of the substance. The path length (b) represents the distance the light travels through the sample and is typically measured in meters. The molar concentration (C) of the absorbing species in the sample is determined based on the absorbance and the known values of ϵ and b.
Additionally, the absorption cross-section (σ), quantifies a molecule’s ability to absorb photons. It is derived from the molar absorption coefficient and Avogadro’s number.
For accurate concentration determination, the absorption cross-section must be incorporated into the calculations. To achieve this, gas adsorption constants under isothermal conditions were evaluated using the Langmuir equation. The HITRAN application programming interface (HAPI) was employed to compute the integrated cross-section area of methane within band 12 wavelengths. HAPI, a python-based set of routines, provides remote access to HITRAN online data and functionality, enabling the calculation of various spectral functions, including absorption, transmittance, and radiance, based on the Hartmann–Tran (HT) profile [
46]. The absorption cross-section is pivotal for these calculations, as it underlies the simulation of spectral interactions. With Equations (1)–(3), the enhancement in methane column concentration ∆
Xch4 (referred to as “column concentration” hereafter) can be calculated using Equation (4).
The radiative transfer model incorporates factors such as variable surface height, solar zenith angle (SZA), and instrument viewing zenith angle (VZA). For each scene, the SZA and VZA are extracted from the metadata accompanying the Sentinel-2 image tiles. This metadata provides the essential geometric information required to accurately calculate the column concentration of each pixel. Our Sentinel-2 retrievals utilize a general approach of determining methane column enhancements by examining fractional changes in the top-of-atmosphere (TOA) reflectance for band 12, when compared with a reference image where there is no, or significantly less, shortwave infrared absorption by a methane plume. By analyzing these fractional changes, methane column concentrations can be derived from any radiometric measurement proportional to TOA radiance. For this study, TOA reflectance from the Sentinel-2 Level-1C (L1C) products with cloud percentages of below 30% was used. The reason 30% percent cloudiness was selected as the threshold is to use higher quality images with the least amount of noise that derives from the presence of water vapor and clouds interacting with the light on the location. One additional criterion for processing the images was that they must all have captured the entire WWTP.
Gorroño et al. had previously validated the methane retrieval process by applying it to simulated emissions across a variety of locations, with a particular emphasis on point sources with notably higher emissions [
25]. Furthermore, this retrieval method has been implemented in numerous investigations to investigate methane emissions from a variety of sources [
47,
48,
49]. The concentration patterns from fugitive and multi-source emissions with lower emission levels at the WWTP were analyzed using this retrieval approach in our study. The use of satellite imagery for exhaustive GHG emission reporting and legislative applications could be advanced by the future validation of smaller, sporadic emission sources.
In this study, 181 images of the WWTP site over five years were gathered from 2019 to 2023. For each year, depending on the quality of the images and number of revisits, a minimum of 33 and a maximum of 44 images were processed, converting the pixel values as described to the mol per meter square column concentration. Phase one of the processing was performed on the images to determine the column concentration of the methane (mol/pixel) for each pixel using Equation (4). In the second phase of data processing, the seasonal, temporal, and annual characteristics of the dataset were investigated using image processing tools.
3.4. Seasonal and Temporal Analysis
The methane column concentration data obtained from the first phase of data processing were in the form of a tensor V with dimensions h × w × c, where h and w represent the height and width of the area of interest, respectively, and c represents different time captures. Because the intervals between captures are inconsistent, VI was generated by interpolating V along the c axis with dimensions h × w × t, where t denotes the number of days since the first capture, covering the entire period from the first to the last day of capture.
The mean value of each pixel in
VI is calculated over the
t axis, yielding
Mi,j ∈ R
h×w. This is expressed in Equation (5).
To derive the trends of column concentration for the entire WWTP using the whole data set,
TW,i ∈ R
t, is calculated by averaging
VI across the
h and
w axes, yielding the mean emission for all pixels over time. This is expressed in Equation (6).
To calculate the seasonal average column concentration similar to the average emission of all data, the seasonal average
Si was calculated by averaging the column concentration
Si,j,k in each pixel (
i,j) over the entire plant,
Si,j,k ∈ Rh×w, where
i is the season index (1–4, corresponding to spring through winter). The seasonal average map was generated by averaging each pixel’s value over the time period within the respective season. This is formally expressed in Equation (7).
Di represents the days within season
i.
A flowchart model of the entire data processing procedure is shown in
Figure 2. There are two unique parts to the workflow: the first phase includes preprocessing and cleaning of the raw data, and the second phase is the analysis and interpretation of the processed data. This figure illustrates each step as well as the relationships between the various phases of data handling.
4. Results and Discussion
This study demonstrates the feasibility of obtaining long-term, high resolution methane variation within a WWTP through top-of-atmosphere observation. While the results may not precisely reflect ground-level concentrations, they reveal the heterogenous nature of source distribution. When combined with ground measurement, these findings can help identify major sources by pinpointing emissions hotspots within the WWTPs. The results illustrate the seasonal patterns in emissions and identify the locations with relatively higher methane emissions by characterizing the temporal variation, spatial distribution, and intensity of methane column concentrations from 2019 to 2023, using the data from the Sentinel-2 satellite. The findings of this study provide valuable information for identifying potential methane emission sources in order to develop efficient mitigation approaches.
Figure 3 shows the year-to-year location and strength of methane hotspots within the WWTP from 2019 to 2023. To better visualize the emission locations, the color scale was set to be 6–12, meaning that the pixels with values below 6 and above 12 are assigned the same color as for 6 and 12, respectively. The color map in
Figure 4 also follows the same color scale approach. The values of column concentrations were calculated for each pixel, based on bands 11 and 12 of Sentinel-2 with a spatial resolution of 20 × 20 m (i.e., 400 m2 per pixel).
Figure 3 shows that the spatial distribution of methane hotpots within the plant was similar across multiple years. Methane hotspots were primarily located in the older sludge process area, as well as the southeast and southwest areas where the activated sludge tanks and clarifiers are.
To perform a thorough evaluation of the overall dynamics of methane column concentration at the WWTP site, both seasonally weighted and non-seasonally weighted averages were calculated based on the mean column concentration of each pixel over the five-year period. Due to variations in the quantity of available data sets for each season, we implemented weighting factors to account for seasonal discrepancies. This method guaranteed that each season contributed proportionately, facilitating a more balanced and precise seasonal comparison. This analysis enabled a more detailed comparative mapping of methane emissions between seasons. For example, while the seasonally weighted and non-seasonally weighted averages were found to indicate that the locations of hotspots in the dataset were generally consistent across methods, a key difference emerged in the magnitude of the column concentrations. The weighted averages provided a more accurate representation of the strength of these hotspots.
These calculations were performed using the formulae outlined in the methodology section. It is important to note that the amount of available satellite data varied by season, being higher in spring and summer and lower in fall and winter. This discrepancy is likely due to the increased frequency of cloudy days during the cooler months, which obscure the satellite imagery and reduce data acquisition opportunities. Consequently, this seasonal variation in data availability could have biased the observed trends, particularly for methane concentration in the fall and winter seasons.
As
Figure 4 indicates, the methane hotspots over the ground surface were primarily confined to the older sludge process area in the northwest and to the activated sludge tanks and clarifiers located in the southeast and southwest corners of the plant, respectively. To gain deeper insight, a thorough investigation into the processes or operations of each location could provide valuable information on how to reduce methane emissions.
To further investigate the seasonal changes in the strengths of the hotspots and emissions based on column concentrations, the seasons were categorized by month: spring (April to June), summer (July to September), fall (October to December), and winter (January to March). Please note the definition of the seasons has one month shift from that defined by air temperature, (i.e., December to February as winter). The reason for this shift is based on the consideration of the delay in temperature change in the sewage following the temperature change of the weather. In addition, the impact of temperature change on the biological process, including both the microbial communities and the metabolic activities, is also expected to take some time to be revealed. Therefore, re-defining the season allows for a realistic observation of seasonal variations in methane emission hotspots throughout the study period.
The results are presented in
Figure 5, which illustrates the monthly average and 3-month moving average of column concentrations in hotspot areas in each month as well as 3-month moving average wastewater temperature. A hotspot area is defined as an area where pixels indicate that the column concentration is two standard deviations above the average value for all of the pixels over the five-years period. Identifying hotspots that vary significantly across months can indicate periods of higher-than-normal emission activities, which may correlate with operational cycles or environmental factors such as temperature and humidity. This refined method for estimating emissions provides a crucial framework for optimizing monitoring efforts.
Figure 5 illustrates the temporal variation of methane column concentrations from 2019 to 2023. Higher column concentrations tend to peak in the fall, reaching as high as 25 mol/pixel.
Figure 5 reveals that methane emissions peaked approximately 2–3 months after the wastewater temperature peaked. Because methane generation is a result of microbial activities within the wastewater, the observed “lagging” of methane emissions in response to the change of the wastewater temperature is expected. A thorough interpretation of the results requires a detailed evaluation of the plant influent, effluent, and operational data, as well as the physical–chemical and biological processes, and thus is beyond the scope of this paper. The overall observation is that the temporal variations of CH
4 column concentrations are positively correlated with the wastewater temperature, which is in agreement with the wastewater biological processes. During the warmer seasons when the temperatures are relatively higher, the biological activities are at a higher level, thus producing more CH
4.
In addition to temporal variation, the spatial distribution of CH
4 concentration over the plant site has been observed to be extremely heterogeneous and time varying. To evaluate how the CH
4 emission hotspots change over time, the images were grouped by season and the CH
4 strength of the images taken in the same season was averaged, as demonstrated in
Figure 6. The same color scale approach as
Figure 3 was used to better visualize the emission locations.
Figure 6 shows that the primary hotspots in spring are located near the southeast secondary clarifiers and the older anaerobic digester and primary sludge fermenter tanks. In summer, the hotspots extend to cover more process units including the southwest activated sludge tanks and clarifiers and the northeast secondary clarifiers. The fall emission locations are similar to summer, except that additional emission sources are shown around the primary clarifiers in the northwest corner. The winter emission pattern shows a more scattered distribution of hotpots. These findings unveil the highly dynamic nature of CH
4 emissions, both temporally and spatially.
These findings emphasize the necessity of conducting continuous multi-point monitoring of fugitive methane from wastewater treatment plants for capturing the temporal–spatial variation in methane emissions. In addition, the observed seasonal shifting in hotspot locations and emission intensities suggests that, when planning for mitigation measures, the emission dynamics should be considered to enhance the effectiveness.
5. Discussion and Conclusions
The results demonstrate the value of satellite remote sensing images in characterizing the spatial and temporal variation of methane emission from relatively smaller-scale industrial GHG sources like WWTPs. For the first time, this study provides a more complete, measurement-based methane emission characterization for WWTPs, revealing both the spatial and temporal variations of methane emissions.
Potential errors are introduced by factors such as atmospheric conditions, vertical methane distribution, and instrument resolution when satellite data are employed to infer surface-level methane concentrations. Satellite observations are still valuable in environmental monitoring due to their broad spatial and temporal coverage, which is not possible with ground-based measurements alone, despite the fact that these limitations influence predictive accuracy. The method can be justified in regulatory contexts by explicitly communicating these uncertainties and validating findings with ground data. Although the use of satellite data may not provide precise predictions, it effectively identifies large-scale trends and emissions anomalies, making it a powerful complementary instrument for regulatory and monitoring efforts, rather than a standalone solution.
The high spatial resolution of Sentinel-2 allows for the detection and quantification of methane column concentrations from emission sources within the plant site, enabling a detailed assessment of their environmental footprint. It is particularly useful in places where ground monitoring would either be cost prohibitive or logistically impossible to perform. Being able to monitor and reduce the emissions of small sources, like WWTPs, is of importance when seeking to meet global carbon neutrality objectives because of the large number of such types of sources [
50].
By closely examining the operational dynamics at specific hotspots, potential operational inefficiencies or leaks can be identified. The findings of this study also provide regulatory implications in multiple ways. For example, instead of setting annual emission standards, the limitations can be set for different seasons, or only for the high emission months, to alleviate the costs associated with both monitoring and implementation of mitigation strategies.
Analyzing the five-year satellite imagery of methane emissions at the WWTP site provides valuable insights into the temporal and spatial dynamics of methane generation within the plant. This indicates the importance of performing continuous measurements to estimate the annual CH4 emissions.
Finally, although Sentinel-2 satellite imagery is valuable in assisting in localizing methane emission sources and estimating plumes based on the column concentrations, whether or not this approach, by itself, is capable of quantitatively evaluating the emissions still needs to be investigated. Seasonal data variability, affected by cloud cover and other atmospheric interferences, particularly impacts data availability in cooler months. These seasonal gaps may influence observed trends, potentially underrepresenting methane levels in fall and winter. Additionally, while Sentinel-2 provides high-resolution spatial data, quantitative methane estimations carry uncertainties without ground-based validation. Future research incorporating meteorological data and ground-level measurements would strengthen the reliability of satellite-based methane quantification, enabling more precise assessment and monitoring of WWTP emissions.