4D Models Generated with UAV Photogrammetry for Landfill Monitoring Thermal Control of Municipal Solid Waste (MSW) Landfills

: The management of the increasing volume of municipal solid waste is an essential activity for the health of the environment and of the population. The organic matter of waste deposited in landfills is subject to aerobic decomposition processes, bacterial aerobic decomposition, and chemical reactions that release large amounts of heat, biogas, and leachates at high temperatures. The control of these by-products enables their recovery, utilization, and treatment for energy use, avoiding emissions to the environment. UAVs with low-cost thermal sensors are a tool that enables the representation of temperature distributions for the thermal control of landfills. This study focuses on the development of a methodology for the generation of 3D thermal models through the projection of TIR image information onto a 3D model generated from RGB images and the identification of thermal anomalies by means of photointerpretation and GIS analysis. The novel methodological approach was implemented at the Meruelo landfill for validation. At the facility, a 4D model (X,Y,Z-temperature) and a 13.8 cm/px GSD thermal orthoimage were generated with a thermal accuracy of 1.63 ◦ C, which enabled the identification of at least five areas of high temperatures associated with possible biogas emissions, decomposing organic matter, or underground fires, which were verified by on-site measurements and photointerpretation of the RGB model, in order to take and assess specific corrective measures.


Introduction
Urbanization and industrial development lead to an increase in municipal solid waste (MSW) [1], and managing the growing volume of waste is a worldwide problem.Waste disposal is the last of the alternatives in the hierarchy of MSW management alternatives [2].However, despite the promotion of recycling and incineration initiatives to reduce the volume of waste, the most commonly used alternative is landfilling [3].Considering the growth in the volume of waste [4], as well as the potential risk to the environment and human health [5], the control and monitoring of landfills is an essential activity to avoid the materialization of such risks.
Landfills can be defined as bioreactors [6] where the stored organic matter reacts as a consequence of storage conditions, exposure to aerobic decomposition processes, aerobic bacterial decomposition, and chemical reactions [7].Exothermic reactions release large amounts of heat, which is the main by-product of the spontaneous decomposition of organic matter in landfills, along with leachate and landfill gas [8].The percentage and volume of each by-product generated depends on several factors: the characteristics of the waste, such as its age and composition, the chemical products present, and environmental factors related to temperature, humidity, and oxygen content [9].The gas generated by the decomposition of organic matter is called biogas or landfill gas (LFG).It is composed of 45-60% methane (CH 4 ), 40-60% carbon dioxide (CO 2 ), 2-5% nitrogen (N2), 0.1-1% oxygen (O 2 ), and other gas fractions in lower percentages, including ammonia (NH 3 ), carbon monoxide (CO), hydrogen (H), sulfides, and non-methane organic compounds (NMOC) such as benzene, vinyl chloride, and trichloroethylene [10].
The set of decomposition processes and exothermic reactions reach temperatures of 30-45 • C under normal landfill operating conditions [11].However, it is possible for the temperature to rise to values of 55-70 • C, and these temperature peaks can be maintained for several years [12].It is estimated that 40% of the LFG leaks into the atmosphere [13], releasing large amounts of methane, a greenhouse gas with a global warming potential 28 times higher than carbon dioxide [14].The hot biogas generated in the deep areas of the landfill migrates to the surface, passing through the waste and increasing the temperature in the landfill piles [15].As a consequence of the thermal processes typical of landfills that occur at depths between 20 and 60 m [8], their upper zone reaches temperatures up to 4-10 • C higher than those of neighboring areas [16].
The accumulated heat from exothermic reactions and the volatility of the gases generated make landfills prone to surface fires and the slow-burning of MSW below [10].High temperatures affect the proper functioning of landfills, influencing gas generation, leachate quality, slope stability, and landfill integrity, and may also affect human health and the environment due to gas emissions into the atmosphere and the release of incomplete combustion products as a result of fires [17,18].This requires monitoring and control of the landfill for at least 20 years after landfill closure [19].
The set of reactions and processes that take place in landfills induce variations in their surface temperature.These thermal anomalies of the surface zone associated with temperature increases generate hot spots that can be directly related to methane emissions to the atmosphere or to superficial or underground fires [20].The identification of anomalies is an essential activity during landfill monitoring to avoid the risk of fire and maximize gas recovery, avoiding environmental effects and guaranteeing the landfill's energetic use [21].Conventional techniques for the detection of thermal anomalies and landfill gas emissions are time consuming and have limited accuracy [22].In light of this problem, unmanned aerial vehicles (UAV) have emerged as an alternative capable of overcoming the difficulties of current techniques [23].The advantages of size, weight, flexibility, speed, and cost [24] offered by this type of equipment have enabled its implementation in the management of municipal solid waste landfills and dumpsites [25].UAVs have been used for landfill growth monitoring [26], topographic and volumetric control of landfills [27], safety and integrity inspections of facilities [28], monitoring of LFG emissions [29], and identification of possible fires [10], and assessing the evolution of emitted environmental pollution [30].This study proposes the control and identification of hot spots, an essential activity for MSW landfill managers, which enables the detection of LFG emissions and the channeling of them through collection wells for utilization or incineration.On the other hand, other techniques can be used to reduce LFG production, such as the aeration or gasification of organic matter prior to landfill disposal.
The miniaturization of sensors has enabled UAVs to carry thermographic sensors [31] to capture thermal images (TIR) that determine the temperature of the imaged surface.The combination of UAV images and photogrammetry enables the representation of temperature distributions in thermal heat maps, thermal orthomosaics [32], and thermal 3D models [33] from the reconstruction of the scene based on two-dimensional images [34].The combination of UAV thermal images and photogrammetry has been used for the identification of hot spots associated with different risk factors such as LGF emissions [1,22] and landfill fires [10,35] based on 2D orthoimages.However, 3D thermal models and representations of the surface are a tool with great potential since, in addition to providing temperature information, they enable representation of the position in space and the real geometry of the object, benefiting from the three-dimensional component of space, which favors the complete interpretation of the area and the approach of the proposed solution, since the studies are carried out using 4D (X,Y,Z-temperature) representations, in which the X,Y,Z dimensions represent the position of the point in space and the fourth dimension (T) represents the surface temperature of the point.
Similar issues relating to the thermal monitoring of landfills have been previously addressed by other authors.Firstly, Lewis et al. [13] analyzed the different existing methodologies for leakage detection and control, paying special attention to infrared thermography techniques.The authors evaluated the capability of infrared thermography for LFG leak detection in an investigation based on the analysis of images taken with hand-held sensors, concluding that the technique at that time could only be used as an approximate alternative and not as an accurate tool for detecting LFG leaks.Subsequently, Yan et al. [3] proposed a methodology for MSW disposal site monitoring using satellite imagery.The proposed technique, which is based on multi-temporal Landsat imagery, enables the monitoring of the land surface temperature of MSW landfills, showing how the temperature of an analyzed landfill is more than 4 • C higher than the immediately surrounding areas.A similar alternative was later proposed by Nazari et al. [35], who developed a methodology for the identification of surface fires and their migration in landfills using satellite thermal images of moderate spatial resolution.Their proposed technique is based on the processing of Landsat images using an in-house MATLAB code to generate multitemporal thermal maps of the study area to identify and flag hot spots and thermal anomalies.Along the lines proposed in the aforementioned research, the use of UAVs in this field has also been evaluated by other researchers.Lega and Napoli [19] introduced an integrated system based on a sensor array for a UAV, combining an infrared camera and air quality sensors to detect thermal emissions and physical and chemical parameters with an orientation toward landfill monitoring.The proposed system enables the monitoring of concentrations, spills, and leachates for the management of waste dumps and illegal discharges.The use of UAVs equipped with methane detectors was proposed by Emran et al. [36] to discover and monitor methane leaks over large areas with remote access in a rapid and economical manner.The proposed technique was based on the use of a multirotor UAV with an offthe-shelf laser-based methane detector for the generation of landfill methane concentration maps quickly and efficiently.Fjelted et al. [37] investigated and analyzed the relationship between ground surface temperature and LFG emissions in landfills under different atmospheric conditions.Their analysis was based on contrasting measurement campaigns under different conditions for capturing UAV images; complementing measurements of methane and carbon dioxide emissions were taken at each measuring point with a static flux chamber to assess the evolution of the results with GIS.Monster et al. [38] compiled different methodologies for the measurement of CH 4 emissions in a review, considering vertical soil gas concentration profiles, surface flux chambers, radial plume mapping, and mass balances using aerial measurements or differential absorption LiDAR, among other methods.The review collects alternatives to infrared sensors and CH 4 concentration instruments for UAVs, considering the present alternatives to not be very accurate in creating 2D maps of concentrations or carrying out field infrared surveys by capturing TIR images with UAVs.Tanda et al. [1] proposed a novel approach for UAV TIR image processing to generate heat maps for the detection of thermal anomalies in surface temperatures and their relationship with local biogas escaping in landfills that leave a thermal footprint.Their research focuses on the determination of methane flows in hot spots identified through visual analysis and photointerpretation of two-dimensional heat maps.More recently, Sliusar et al. [25] developed a comprehensive review of drone technology in MSW landfill management and control.Their paper mentions how most of the publications related to this technology focus on monitoring the emissions of landfill gas or its individual components, mainly methane.In relation to the subject matter, it highlights the capacity of aerial thermal imaging to generate heat maps, orthoimages, and 3D thermal models of the surface, but it does not detail the processing methodology or workflow to obtain them.Complementarily, the review highlights the scarce experience of the current technique in analyzing UAV photogrammetry-derived products using GIS tools in landfill areas despite the multiple advantages of this type of imagery [39].
The present study focuses on the development of a methodological proposal for the generation of 4D models of a landfill surface using UAV thermal images for the identification of thermal anomalies through the processing of images with Geographic Information Systems (GIS).The methodology is developed through its application to a specific case, with the modeling of an active MSW landfill for its control and periodic monitoring, enabling a follow-up of the proper functioning of the LFG collection systems.In case hot spots are identified, it allows their geolocation for the application of corrective measures by the landfill regulator.
Conventionally, the activities at a landfill consist of a set of operations for the disposal, compaction, covering, and control of sediments.Within this set of tasks is thermal control, which includes biogas control.In order to reduce LFG emissions, different strategies can be employed for MSW management in landfills.There are composting and gasification strategies for the organic matter fraction of MSW prior to landfill, which enable both the extraction of biogas and reductions in the volume of waste in the landfill.However, despite the various filters and treatments applied, it is impossible to eliminate all of the organic matter, which, once deposited, decomposes and generates LFG.At this point, biogas control also requires minimizing subsurface migration and atmospheric emissions, protecting groundwater, reducing the risk of fire, protecting structures, and gas collection for energy use.To this end, during waste compaction, pathways are created to allow the flow of biogas and leachate to drainage areas or catchment pits, minimizing risks.This system can be passive or active, depending on whether the migration is done naturally through pressure differences or is forced [21].
This set of processes induces increases in landfill temperatures.There are three main strategies for temperature management: heat extraction, regulation, and supplementation.The different strategies enable both the use of heat and biogas to cover energy needs or to produce energy as well as the regulation of the landfill temperature around 35-45 • C to maximize biogas production and its subsequent use [40].Under normal operating conditions, the landfill temperature is kept stable by the extraction and collection of LFG and the removal of excess heat.From all of the above, it can be deduced that thermal monitoring of MSW is fundamental for the correct operation of landfills, and this monitoring can be carried out by means of thermal orthoimages generated with images obtained with low-cost sensors installed in UAVs.

Methods
The development of UAV uses is limited by the gap between technology and the methodological approaches for application [25].Achieving a simultaneous development requires parallel research in both areas.This study proposes a methodology for the generation of 3D thermal models of the MSW landfill surface using RGB and TIR images captured with low-cost dual sensors for their subsequent interpretation using GIS, oriented to the quasi-automatic identification of thermal anomalies with potential fire risks to estimate their possible behavior.

Data and Thermal Imaging Capture with Low-Cost UAVs
Images are the fundamental basis of photogrammetry for the generation of threedimensional models or orthoimages [41].In the field of UAVs, images are generally captured on the basis of preprogrammed automatic flights that are designed according to the needs of the project objective [42].At present, despite the losses in sensitivity and accuracy associated with low-cost UAVs, they are the most widely used alternative in several investigations [43].This study proposes the use of dual sensors integrating an infrared thermal sensor and a non-metric visual sensor, enabling the capture of RGB and TIR images simultaneously from the same point of view [44].Working with this type of low-cost equipment requires careful work during the planning and execution of the activities in order to overcome their limitations and obtain accurate results [45].
To capture the images, the flight plan must be designed to adhere to the requirements of the model's objective and the specifications of the platform and low-cost dual sensor to be used.The whole set of flight parameters to be defined must take into account the resolution, sensitivity, and contrast limitations of the most limiting equipment, generally the thermal sensor [46].The flight design is usually performed by specific software automatically by inserting the initial flight parameters [47] described below.For the capture of the photogrammetric block of large areas such as landfills, it is advisable to use nadiral flights with flight heights high enough to achieve a ground sample distance (GSD) suitable for the application [48], in this case the identification of thermal anomalies, where the decimeter accuracies are sufficient to set the detection of a possible hot spot.On the other hand, the frontlap and sidelap used for TIR capture set high values with respect to UAV photogrammetry for visual images, proposing values higher than 70% [49] or even more in both cases [50].The last parameter to consider is the flight speed, which determines the total flight time, a factor to be taken into account due to UAV autonomy limitations.It is set below 3 m/s to avoid resolution, sharpness, and contrast issues associated with TIR images [51].
In order to obtain an accurate thermal model, which also enables the correct quasiautomatic identification of thermal anomalies, the planning also has to take into account the surface emissivity, the sensor stabilization, and the atmospheric conditions of the day of flight.The emissivity of the surface materials will depend on the nature of the MSW, so a suitable value needs to be adapted for the characterization of the entire landfill.Due to the instability of low-cost infrared sensors, it is recommended to set a sensor stabilization time of about 30 min before taking flight [52,53].Finally, atmospheric and climatological conditions determine relevant parameters such as the ambient temperature, wind, solar radiation, and cloudiness.When scheduling the flight, in order to minimize the effect of solar radiation and maximize the temperature difference between the hot spot and the rest of the landfill, the summer season should be avoided [1,4], and flying in the early hours of the day before sunrise or after sunset is preferable [13].
Finally, in order to correctly georeference the RGB and TIR images, ground control points (GCPs), which are one of the most influential parameters in the final accuracy of the model [54,55], are used in a standardized manner.For this purpose, metallic plates can be used, as they are visible in both types of images.These plates must be distributed homogeneously along the landfill to be represented [56,57].In addition, the use of thermal control points (TCPs) is proposed, which enables the contrasting of the temperature measured with the UAV sensor and the surface temperature in order to thermally calibrate the final model, adjusting it more closely to reality [58].For this purpose, at least three TCPs distributed along the surface should be used, from which the temperature is measured with a hand-held sensor or other more accurate device [52].

Photogrammetric Processing of RGB and TIR Images
The application of photogrammetry for the generation of 3D thermal and orthomosaic models with thermal images presents a large gap with respect to the quality of the representations obtained with RGB images, mainly due to the lack of resolution, contrast, and homogeneity in the TIR images [59].The standardized methods for photogrammetry with RGB images do not obtain the same results when working with TIR images.The representations with the temperature variable can be obtained by implementing different strategies [60], from the single application of thermal images to generate the model [61], the texturing of 3D models with TIR images [62], or the implementation of four-band images through the fusion of RGB and TIR images [63].
This study proposes a methodology resulting from the combination of different strategies based on the simultaneous capture of RGB and TIR images with low-cost dual sensors.Based on previous works, the proposal is based on solving the problems in the alignment of the thermal images of the structure from motion (SfM) algorithms used in modern photogrammetry [64] by implementing the optimized external orientation parameters (X, Y, Z coordinates and pitch, yaw, and roll angles) of the RGB images to position the TIR images [46].
The proposed method (Figure 1) is based on the processing of RGB images following the conventional workflow standardized for SfM photogrammetry software.First, the images are aligned by identifying homologous points through scale invariant feature transform (SIFT) algorithms, enabling the generation of sparse point clouds.The optimization of the orientation of the images and georeferencing of the model is achieved through GCP registration, which subsequently obtains a dense point cloud using multi-view stereo (MVS) algorithms [65].Finally, once the point cloud is obtained, the SfM photogrammetry software enables the generation of 3D models, meshes, orthomosaics, and textures.Then, using the thermal images, the model generated with the RGB images is taken as a basis.The TIR images are oriented on it using the optimized orientation parameters of the homologous RGB images to project the radiometric information on the model.Thus, the thermal 3D model is generated, which includes the temperature for each represented point and enables acquisition of a thermal orthoimage that shows the distribution of the surface temperatures of the landfill.
Appl.Sci.2023, 13,13164 6 of 22 This study proposes a methodology resulting from the combination of different strategies based on the simultaneous capture of RGB and TIR images with low-cost dual sensors.Based on previous works, the proposal is based on solving the problems in the alignment of the thermal images of the structure from motion (SfM) algorithms used in modern photogrammetry [64] by implementing the optimized external orientation parameters (X, Y, Z coordinates and pitch, yaw, and roll angles) of the RGB images to position the TIR images [46].
The proposed method (Figure 1) is based on the processing of RGB images following the conventional workflow standardized for SfM photogrammetry software.First, the images are aligned by identifying homologous points through scale invariant feature transform (SIFT) algorithms, enabling the generation of sparse point clouds.The optimization of the orientation of the images and georeferencing of the model is achieved through GCP registration, which subsequently obtains a dense point cloud using multi-view stereo (MVS) algorithms [65].Finally, once the point cloud is obtained, the SfM photogrammetry software enables the generation of 3D models, meshes, orthomosaics, and textures.Then, using the thermal images, the model generated with the RGB images is taken as a basis.The TIR images are oriented on it using the optimized orientation parameters of the homologous RGB images to project the radiometric information on the model.Thus, the thermal 3D model is generated, which includes the temperature for each represented point and enables acquisition of a thermal orthoimage that shows the distribution of the surface temperatures of the landfill.

D Thermal Model Calibration
The 3D thermal model obtained from the proposed workflow enables representation of the temperature for each point of the generated cloud.However, the temperature determined by the UAV sensor may differ from the real temperature of the modeled surface due to the large distance between the sensor and the captured point, the atmospheric conditions, the viewing angle, or the stability of low-cost equipment.To correct the error and obtain an accurate representation of the surface temperature distribution of the MSW piles of the landfill, the implementation of TCPs is proposed, which provides a higher accuracy than that obtained with the UAV sensor [58].
The temperature correction of the model is performed by determining the temperature deviation of the UAV model points with respect to the complementary sensor measurements, establishing an error that is applied through a correction in the model transformation.It can also be applied with a raster transformation in the GIS tool for the calibration of the thermal orthomosaic during its analysis.

Photointerpretation Analysis of Thermal Orthomosaics with GIS
The thermal orthoimage obtained is a two-dimensional raster representation showing the surface temperature gradient of the landfill in false color using a predetermined palette.The analysis and evaluation of the orthoimage for the identification of thermal anomalies associated with potential LFG emissions, fire hazards, or decomposing organic matter is performed through photointerpretation via visual analysis of the image, using GIS tools for correct visualization and editing.
Starting from the thermal orthomosaic and the RGB, a first photointerpretation of the image is proposed to identify hot spots according to the temperature gradient represented, which are associated with possible thermal anomalies.The identification of hot spots is made on the basis of a quantitative analysis via photointerpretation of the color code of the orthomosaic in false color.Areas whose color is associated with high temperatures are identified following the palette of the representation.This interpretation may be conditioned by the experience of the person in charge of this analysis and, therefore, the results may vary.
The area identified in GIS can be georeferenced and marked by alternating layers in the tool for its interpretation and the detection of possible causes in the RGB orthomosaic.The evaluation of the thermal anomaly in both images will enable justifying the determination of the cause of the temperature difference recorded, either through the existence of uncovered geomembranes, plastics, metallic elements, or pipes in the LFG collection system, as these elements are able to reach higher temperatures than the rest of the environment as a result of mere solar exposure without gasification; decomposition; or fire processes.In this way, the aim is to discard the areas with controlled behavior among thermal anomalies including potential risks such as fires, LFG emissions, etc., so that the appropriate corrective measures can be taken.

Method Validation
The methodology developed for the thermal control of landfills using 3D thermal models can be validated by implementing it in a real case study.Therefore, the application of the workflow to a MSW landfill for the identification of thermal anomalies was proposed.Once the thermal representation was generated and the hot spots with potential risk were identified, they were contrasted with the RGB representation to ensure the absence of a foreign element or surface that could be causing a false temperature measurements due to emissivity conditions different from the rest of the surface.Finally, a subsequent field visit verified the temperature measured by a complementary method and its subsequent in situ evolution analysis in order to corroborate the existence of the risk in order to take the necessary corrective measures, thus validating the proposed methodology.

Study Area
For the validation of the proposed methodology, the workflow was implemented in a landfill that met the appropriate conditions for its application.The Meruelo landfill is an integrated waste treatment center located in the autonomous community of Cantabria (Spain).It is a landfill whose waste is mainly organic (from MSW) or similar to urban waste, with a total capacity of 60,000,000 tons.The objective of the facility is to minimize the volume of MSW in landfill, for which it first implements waste pre-selection with a capacity of 65 t/h, separating recyclable materials such as cardboard, paper, and plastics, by means of magnetic and induction separators, from the organic matter.The organic matter from the recovery is subjected to waste energy recovery processes through degassing operations for the production of biogas and the generation of electrical and thermal energy with a cogeneration plant [66].The facilities have an area of about 60 ha, and are divided into the MSW treatment and valorization plant, the leachate treatment plant, and the non-hazardous waste landfill (Figure 2).The latter is the area to be modeled from the UAV images.
For the validation of the proposed methodology, the workflow was implemented in a landfill that met the appropriate conditions for its application.The Meruelo landfill is an integrated waste treatment center located in the autonomous community of Cantabria (Spain).It is a landfill whose waste is mainly organic (from MSW) or similar to urban waste, with a total capacity of 60,000,000 tons.The objective of the facility is to minimize the volume of MSW in landfill, for which it first implements waste pre-selection with a capacity of 65 t/h, separating recyclable materials such as cardboard, paper, and plastics, by means of magnetic and induction separators, from the organic matter.The organic matter from the recovery is subjected to waste energy recovery processes through degassing operations for the production of biogas and the generation of electrical and thermal energy with a cogeneration plant [66].The facilities have an area of about 60 ha, and are divided into the MSW treatment and valorization plant, the leachate treatment plant, and the non-hazardous waste landfill (Figure 2).The latter is the area to be modeled from the UAV images.

UAV Photogrammetric Flight
The photogrammetric flight for the capture of the RGB and TIR images was designed following the established recommendations for image capture with low-cost UAV equipment.The flight began at 9:00 a.m. with the image capture, after a previous wait of 30 min for the stabilization of the sensor.During the flight, the atmospheric conditions were

UAV Photogrammetric Flight
The photogrammetric flight for the capture of the RGB and TIR images was designed following the established recommendations for image capture with low-cost UAV equipment.The flight began at 9:00 a.m. with the image capture, after a previous wait of 30 min for the stabilization of the sensor.During the flight, the atmospheric conditions were favorable for data collection; the ambient temperature was 21 degrees and the skies were clear with no wind.
The platform used was a quadcopter multirotor, specifically a DJI Mavic 3T, which has an integrated dual M3T sensor (Figure 3).The main features and technical information of the UAV platform and the dual sensor cameras are shown in Tables 1-3.
The calibration of the RGB and thermal sensors was guaranteed by the company contracted for the photogrammetric flight.In addition, the self-calibration tool implemented by Agisoft Metashape was applied during the photogrammetric processing of the images for the generation of the 3D model.This step eliminates possible remaining lens distortions to ensure greater accuracy of the representation.
Appl.Sci.2023, 13, 13164 9 of 22 favorable for data collection; the ambient temperature was 21 degrees and the skies were clear with no wind.
The platform used was a quadcopter multirotor, specifically a DJI Mavic 3T, which has an integrated dual M3T sensor (Figure 3).The main features and technical information of the UAV platform and the dual sensor cameras are shown in Tables 1-3.The calibration of the RGB and thermal sensors was guaranteed by the company contracted for the photogrammetric flight.In addition, the self-calibration tool implemented by Agisoft Metashape was applied during the photogrammetric processing of the images for the generation of the 3D model.This step eliminates possible remaining lens distortions to ensure greater accuracy of the representation.
The automatic flight for image acquisition was designed using the proprietary software of the UAV to guarantee maximum compatibility.The sensor was configured with an inclination of 30° with respect to nadir to adequately capture the steep slopes of the  The automatic flight for image acquisition was designed using the proprietary software of the UAV to guarantee maximum compatibility.The sensor was configured with an inclination of 30 • with respect to nadir to adequately capture the steep slopes of the landfill.The automatic flight designed and its main specifications are shown in Figure 4 and Table 4.
Complementarily, in order to guarantee the correct measurement of the surface temperature of the landfill, the thermal sensor was configured with an emissivity of 0.94 [67], which can be considered the average for the entire surface of the landfill due to the homogeneity of the MSW arrangement.
In addition to capturing the images, a total of 17 GCPs were distributed along the landfill to orient the images, scale, and georeference the model.The GCPs were evenly distributed throughout the study area at different heights.The coordinates of the GCPs were taken with GPS using a Leica GS-15 GPS device (Leica Geosystems, San Galo, Switzerland) in an official reference system (ETRS89-UTM Zone 30).Additionally, during the flight, the temperature was measured with a hand-held sensor at 10 geolocated points distributed throughout the study area as ground truth for the calibration and thermal adjustment of the 4D model.landfill.The automatic flight designed and its main specifications are shown in Figure 4 and Table 4. Complementarily, in order to guarantee the correct measurement of the surface temperature of the landfill, the thermal sensor was configured with an emissivity of 0.94 [67], which can be considered the average for the entire surface of the landfill due to the homogeneity of the MSW arrangement.
In addition to capturing the images, a total of 17 GCPs were distributed along the landfill to orient the images, scale, and georeference the model.The GCPs were evenly distributed throughout the study area at different heights.The coordinates of the GCPs were taken with GPS using a Leica GS-15 GPS device (Leica Geosystems, San Galo, Switzerland) in an official reference system (ETRS89-UTM Zone 30).Additionally, during the flight, the temperature was measured with a hand-held sensor at 10 geolocated points distributed throughout the study area as ground truth for the calibration and thermal adjustment of the 4D model.

Results
The results of the application of the proposed methodology for the generation of 4D models for the thermal control of landfills are described below.

Results
The results of the application of the proposed methodology for the generation of 4D models for the thermal control of landfills are described below.

4D Thermal Model Generation
As a result of the UAV flight, a total of 1289 RGB and TIR images were obtained.Both sets of images were managed through the SfM photometric process, following the workflow proposed in Figure 1 and using Agisoft Metashape software (version 1.2.7).First, the set of RGB images in JPEG format was used to generate the 3D RGB model and the 3.46 cm/px GSD ortoimage shown in Figure 5, which were georeferenced using the GCPs for better accuracy.
From the previous 3D RGB model, the 3D thermal model was generated following the workflow shown in Figure 1.The orientation parameters of the RGB images were exported to align the thermal images.Subsequently, the thermal information from the oriented TIR images was projected onto the 3D RGB model, giving false color to each of the points of the representation, enabling acquisition of the thermal model (Figure 6) and the 13.8 cm/px GSD thermal orthoimage shown in Figure 7.The set of thermal images in R-JPEG format enabled acquisition of qualitative and quantitative data that represents the temperature gradient in false color.

4D Thermal Model Generation
As a result of the UAV flight, a total of 1289 RGB and TIR images were obtained.Both sets of images were managed through the SfM photometric process, following the workflow proposed in Figure 1 and using Agisoft Metashape software (version 1.2.7).First, the set of RGB images in JPEG format was used to generate the 3D RGB model and the 3.46 cm/px GSD ortoimage shown in Figure 5, which were georeferenced using the GCPs for better accuracy.From the previous 3D RGB model, the 3D thermal model was generated following the workflow shown in Figure 1.The orientation parameters of the RGB images were exported to align the thermal images.Subsequently, the thermal information from the oriented TIR images was projected onto the 3D RGB model, giving false color to each of the points of the representation, enabling acquisition of the thermal model (Figure 6) and the 13.8 cm/px GSD thermal orthoimage shown in Figure 7.The set of thermal images in R-JPEG format enabled acquisition of qualitative and quantitative data that represents the temperature gradient in false color.From the previous 3D RGB model, the 3D thermal model was generated following the workflow shown in Figure 1.The orientation parameters of the RGB images were exported to align the thermal images.Subsequently, the thermal information from the oriented TIR images was projected onto the 3D RGB model, giving false color to each of the points of the representation, enabling acquisition of the thermal model (Figure 6) and the 13.8 cm/px GSD thermal orthoimage shown in Figure 7.The set of thermal images in R-JPEG format enabled acquisition of qualitative and quantitative data that represents the temperature gradient in false color.The 4D model, unlike single images or 2D orthoimages, enables the interpretation of the morphology of the landfill, providing the analysis with altimetric information that favors the understanding of its berms, slopes, and elevations, given its strong steepness.In addition, relevant thermal information is provided.

GIS Analysis for Hot Spot Identification
The orthoimages obtained from the 3D models are analyzed with GIS in order to identify possible thermal anomalies and their causes.Through a first photointerpretation, different hot spots to be evaluated were identified, as shown in Figure 8.The identified points were chosen because these hot spots were considered the most relevant and there was no known prior justification.For example, it is possible to find hot spots in the lower left corner due to ash deposition from the incineration plant during image capture; these were not analyzed.The 4D model, unlike single images or 2D orthoimages, enables the interpretation of the morphology of the landfill, providing the analysis with altimetric information that favors the understanding of its berms, slopes, and elevations, given its strong steepness.In addition, relevant thermal information is provided.

GIS Analysis for Hot Spot Identification
The orthoimages obtained from the 3D models are analyzed with GIS in order to identify possible thermal anomalies and their causes.Through a first photointerpretation, different hot spots to be evaluated were identified, as shown in Figure 8.The identified points were chosen because these hot spots were considered the most relevant and there was no known prior justification.For example, it is possible to find hot spots in the lower left corner due to ash deposition from the incineration plant during image capture; these were not analyzed.
Analysis using GIS functions such as reclassifications, segmentations, or the application of enhancement filters for digital identification of thermal anomalies is possible.The application of enhancing or softening filters allows certain conclusions to be reached.Since the temperature values of the data of interest may be modified, the application of this type of analysis may screen out or eliminate useful information.
Photointerpretation of the analog RGB orthoimage provides the most extensive thermal anomalies identified, such as zones A and D, with temperatures of 41.6 °C and 31.2 °C, respectively.Most of the hot spots identified correspond to areas of the landfill where plastic liners and geotextiles are exposed, awaiting the deposition of new MSW to cover them, as can be seen in Figure 9.   Analysis using GIS functions such as reclassifications, segmentations, or the application of enhancement filters for digital identification of thermal anomalies is possible.The application of enhancing or softening filters allows certain conclusions to be reached.Since the temperature values of the data of interest may be modified, the application of this type of analysis may screen out or eliminate useful information.
Photointerpretation of the analog RGB orthoimage provides the most extensive thermal anomalies identified, such as zones A and D, with temperatures of 41.6 • C and 31.2 • C, respectively.Most of the hot spots identified correspond to areas of the landfill where plastic liners and geotextiles are exposed, awaiting the deposition of new MSW to cover them, as can be seen in Figure 9.In the lower right corner of Figure 8 is the hot spot G, an area of high temperatures (33.8 • C) that is shown in detail in Figure 10.The anomaly in this case may be due to two causes.First of all, the interpretation of the RGB orthoimage enables the deduction that the lack of vegetation is the main cause.Secondly, by analyzing the leachate collection system, it is possible to note that there is a set of subway pipes located in this specific area.The fluids from the decomposition of organic matter in the landfill may have high temperatures, and the buried collection pipe may heat up the soil as a result, or there may be leaks in the piping.
The other hot spots that were identified lack an a priori justified cause derived from the photointerpretation or analysis of the installations.In this case (Figure 11), zones B (33.3 • C) and C (27.3 • C) show homogeneous parts of the surface terrain where there are temperature differences that may indicate a heat source underneath.These anomalies may be due to possible LFG emission zones, decomposition of buried organic matter, or underground fires.
In the oldest zone of the landfill (right side of Figure 8), two other areas of thermal anomalies in need of justification can be identified.Figure 12 shows a stockpile of material that is at a higher temperature (26.4 • C) than the surrounding area, and it is even higher than that of the concrete access next to it, indicating an unusual thermal behavior that needs to be evaluated in depth.
In the lower right corner of Figure 8 is the hot spot G, an area of high temperatures (33.8 °C) that is shown in detail in Figure 10.The anomaly in this case may be due to two causes.First of all, the interpretation of the RGB orthoimage enables the deduction that the lack of vegetation is the main cause.Secondly, by analyzing the leachate collection system, it is possible to note that there is a set of subway pipes located in this specific area.The fluids from the decomposition of organic matter in the landfill may have high temperatures, and the buried collection pipe may heat up the soil as a result, or there may be leaks in the piping.The other hot spots that were identified lack an a priori justified cause derived from the photointerpretation or analysis of the installations.In this case (Figure 11), zones B (33.3 °C) and C (27.3 °C) show homogeneous parts of the surface terrain where there are temperature differences that may indicate a heat source underneath.These anomalies may be due to possible LFG emission zones, decomposition of buried organic matter, or underground fires.In the oldest zone of the landfill (right side of Figure 8), two other areas of thermal anomalies in need of justification can be identified.Figure 12 shows a stockpile of material that is at a higher temperature (26.4 °C) than the surrounding area, and it is even higher than that of the concrete access next to it, indicating an unusual thermal behavior that needs to be evaluated in depth.peratures, and the buried collection pipe may heat up the soil as a result, or there may be leaks in the piping.The other hot spots that were identified lack an a priori justified cause derived from the photointerpretation or analysis of the installations.In this case (Figure 11), zones B (33.3 °C) and C (27.3 °C) show homogeneous parts of the surface terrain where there are temperature differences that may indicate a heat source underneath.These anomalies may be due to possible LFG emission zones, decomposition of buried organic matter, or underground fires.In the oldest zone of the landfill (right side of Figure 8), two other areas of thermal anomalies in need of justification can be identified.Figure 12 shows a stockpile of material that is at a higher temperature (26.4 °C) than the surrounding area, and it is even higher than that of the concrete access next to it, indicating an unusual thermal behavior that needs to be evaluated in depth.Finally, Figure 13 shows a homogeneous zone (F) of high temperature (23.8 °C), similar to hot spots B and C, where the thermal difference cannot be justified by differences in the emissivity or heat capacity of the surface represented.Finally, Figure 13 shows a homogeneous zone (F) of high temperature (23.8 °C), similar to hot spots B and C, where the thermal difference cannot be justified by differences in the emissivity or heat capacity of the surface represented.

Validation of Results
The non-justified hot spots were evaluated in the R-JPEG formatted thermal images with the DJI Thermal Analysis Tool 3 software (version 3.2.0),which enables image analysis and temperature measurements of the photographed surface.The temperature measurement of the two areas of the current landfill (B and C) shows that the hot spots are at 33.3 or 27.3 °C compared to the 23.7 or 19.5 °C temperatures of the surroundings, thus verifying the existence of a hot spot which has to be further evaluated in order to take preventive measures.
In addition, the areas of both hot spots are close to some of the 79 collection wells located along the entire length of the landfill, which are shown in green in Figure 14.

Validation of Results
The non-justified hot spots were evaluated in the R-JPEG formatted thermal images with the DJI Thermal Analysis Tool 3 software (version 3.2.0),which enables image analysis and temperature measurements of the photographed surface.The temperature measurement of the two areas of the current landfill (B and C) shows that the hot spots are at 33.3 or 27.3 • C compared to the 23.7 or 19.5 • C temperatures of the surroundings, thus verifying the existence of a hot spot which has to be further evaluated in order to take preventive measures.
In addition, the areas of both hot spots are close to some of the 79 collection wells located along the entire length of the landfill, which are shown in green in Figure 14.The collection wells capture the LFG for energy recovery, avoiding emissions and other associated safety issues, while eliminating the LFG generated and temperature increases.Despite being close to some of the wells, the temperature in these two areas is high, which may indicate possible problems in the recovery of nearby wells or the need for new wells.
On the other hand, in the two hot spots on the right of the image (E and F), which are The collection wells capture the LFG for energy recovery, avoiding emissions and other associated safety issues, while eliminating the LFG generated and temperature increases.Despite being close to some of the wells, the temperature in these two areas is high, which may indicate possible problems in the recovery of nearby wells or the need for new wells.
On the other hand, in the two hot spots on the right of the image (E and F), which are around the oldest landfill, at a lower elevation, the absence of LFG collection wells can be observed.This condition may imply the emission of biogas to the atmosphere and a consequent increase in temperature as the gas migrates through the deposited MSW layers.In the identified area related to a possible material stockpile, the temperature is 26.4 • C, which is higher than the 24.9 • C of the nearby road and the 18.2 • C of the surrounding vegetation.This anomaly may therefore be due to the decomposition of organic matter, ashes from energy recovery deposited in that area of the landfill, etc.Finally, the hot spot of the homogeneous zone without vegetation (G) shows temperatures of 23.8 • C compared to 19.7 • C in the surrounding cold zones.In this case, in the absence of any other reason for or interpretation of the RGB image, a possible LFG emission can be intuited due to the absence of collection wells in this location.
In short, it can be assured that the 3D thermal model enables the monitoring of the landfill from a thermal point of view, providing managers with a tool that enables the detection of possible hot spots in MSW for evaluation and the application of the most appropriate corrective measures according to the incidence shown.

Discussion
Once the proposed methodology for thermal control by means of low-cost sensors installed in UAVs has been applied to the case considered, and in view of the results obtained from the analysis and evaluation of the 3D thermal model, the following interpretations of the research carried out can be taken into account.
The selection of the sensor conditions the results obtained.The storage format of the thermal images and the type of information that they contain determines the final properties of the model.JPEG or R-JPEG image formats provide false color images that enable qualitative 3D models to be obtained, from which only evaluations based on photointerpretations by the user can be developed, in addition to a quantitative analysis of the isolated images.However, TIFF image formats make it possible to generate quantitative thermal 3D models in which the temperature can be directly determined for each point of the model.The R-JPEG images obtained enable an individual quantitative analysis of each of them, but not the generation of a model with temperature data, resulting in a qualitative thermal 3D model.The selection of the platform conditions the capacity of the resulting model, enabling, in the case of quantitative models, a quasi-automatic analysis by means of GIS tools for the identification of thermal anomalies.The potential of this methodology has to be further evaluated in investigations with greater scope in order to analyze the results when using other more complex equipment; this is beyond the scope of this study, which is focused on a first iteration oriented to the generation and validation of 3D thermal models for the identification of anomalies.
Flight planning based on the atmospheric situation may condition the results.Prioritizing times or days of the year with low temperature and low solar radiation will guarantee a greater temperature difference between the anomalies and the environment, and therefore a greater gradient in the image, which facilitates anomaly identification.Similarly, cloudy days reduce the influence of the sun on the geomembranes, coatings, and pipes, keeping them at the temperature of the environment and avoiding their identification as possible pathologies.Despite the satisfactory results, the proposed method for MSW monitoring is also conditioned by the state of the facility during the UAV flight.In addition to the atmospheric and meteorological variables, the constant dumping of debris and its coverage with geotextile materials can modify the conditions, hiding and masking possible anomalies identified before they can be verified and corrected in the field.Therefore, the methodology has to be implemented periodically and frequently in order to act quickly and provide useful and effective solutions.
The generation of the thermal 3D models following the proposed workflow, which is based on the combination of blocks of images, is fast and simple, solving the conventional problems of TIR image processing.Despite obtaining a 3D thermal model with the highest quality, detail, precision, and density possible, it should be noted that the resulting thermal orthoimage on which to perform the photointerpretation has a lower resolution compared to the RGB orthoimage.This is due to the lower spatial resolution of the thermal image.However, the GSD of the thermal orthoimage is 13.8 cm/px, which is considered sufficient for this type of application, where very high accuracies are not required.In spite of this, the advantages of having the 3D thermal model enable a 3D global conception of the landfill to be obtained, which also enables evaluation of the altimetry, morphology, and metric magnitudes, such as the volumes of the facility.On the other hand, the 3D thermal model generated by following the proposed methodology offers an accuracy of 1.63 • C in comparison with ground-truth measurements carried out with a complementary hand-held sensor in the field.
The analysis of the orthoimage based on photointerpretation is a relatively complicated, time-consuming process that can be subject to experience.Visual analysis is the only alternative when dealing with a qualitative 3D thermal model, which is conditioned by the type of sensor used.In the case of quantitative models, semi-automatic analysis based on reclassifications with GIS tools can speed up the work, reducing time, favoring the identification of anomalies, and avoiding reliance on the interpretation of the naked human eye.
GIS tools are one of the main alternatives used for the management of UAV photogrammetry-derived products.In this case, the software has been used for the correct visualization of the orthoimages and the georeferencing of the anomalies identified, as well as of the LFG collection wells, favoring the model analysis.However, the potential of this tool can be enhanced when using quantitative models, where the potential of the functions that it offers can be exploited in order to automatically identify thermal anomalies based on temperature.
The validation of the results was carried out through a combination of photointerpretation of the thermal and RGB representations of the landfill, field visits, and information about the state of the facility.Although the characterization of the anomalies, their justification, and their correction are not the main objective of this study, more detailed and concise results can be obtained by using complementary sensors, such as methane concentration meters, boreholes for core recovery, and subsurface temperature probes.However, the use of these alternatives is much more expensive, contrasting with the low-cost profile of the proposed UAVs.The proposed method can be used as a decision-making tool in landfill management, providing complementary information when investing in tests involving more accurate and expensive profiling equipment.
The results obtained enable the implementation of 4D models in thermal control applications as a more complete tool than the heat maps proposed in previous publications.The proposed alternative enables the establishment of a methodology for the identification of thermal anomalies using only one type of sensor.In this way, the necessary resources and information capture can be optimized compared to alternatives that propose the use of different sensors [19].The proposed UAV terminal imaging application facilitates progress in the field of drone technology, providing a more efficient solution for the thermal monitoring of landfills compared to techniques based on hand-held sensors [13] or satellite imagery [3,35], where spatial resolution is compromised.Unlike other investigations, the obtained representations also enable the determination of possible relations between the anomaly and LFG emissions based on more complete 3D models than those alternatives based on the interpretation of isolated UAV images [37] or 2D representations, such as heat maps [1] or ortho-thermal images [1], enabling the interpretation of altimetric information for better understanding and decision making.Finally, unlike the rest of the reviewed research [25], the methodology of photogrammetric processing used to generate a 4D model is detailed, facilitating its application in other case studies.
The methodology is based on the interpretation of a single model based on one flight.The implementation of several periodic flights over time would allow the generation of dynamic models, showing the evolution of the landfill temperature over time, as well as a control of the anomalies identified for a better monitoring of the temperature.
The results obtained from the analysis of the model and the current situation of the LFG collection system have made it possible to identify areas where the thermal anomaly lacks a simple a priori explanation.This requires specific measures to be taken for an in-depth evaluation of each anomaly that enable its justification and the corresponding corrective measures to be taken, such as the analysis of discharges, in situ gas measurements, etc.

Conclusions
Landfills are essential facilities for MSW management.They not only enable reductions in the risk of environmental and human health impacts, but also the acquisition of a certain energy benefit through the valorization of waste or LFG.The thermal control of these facilities is an activity that guarantees their correct operation, and must be carried out during the whole life cycle of the facility in order to reduce environmental impacts and maximize energy production.The implementation of the methodological proposal developed in this study enables the application of 3D thermal models generated with images obtained with low-cost UAV sensors to the identification of thermal anomalies, favoring decision making for the management of MSW landfills, the sizing of corrective measures, and the prediction of fire evolution.
This methodology enables the acquisition of 3D thermal models and orthomosaics which represent the position, geometry, and temperature of the landfill surface quickly, through photogrammetric processing of homologous RGB and thermal images captured simultaneously with dual sensors.The 3D thermal model is generated by projecting the radiometric information of the TIR images, based on the orientation of the homologous RGB images, onto a previously generated RGB model.The analysis of the orthoimage and thermal model enables identification through the photointerpretation of hot spots where the temperature, higher than that of the environment or surrounding areas, may imply fires, LFG emissions to the atmosphere, or uncontrolled decomposition of organic matter, all of them compromising the stability and safety of the landfill.
The analysis, based on photointerpretation, can be complemented by georeferencing the anomalies in GIS tools, which also enables the evaluation of the possible causes through parallel photointerpretation of RGB orthoimages and the superposition of other landfill management elements, such as the pipes or wells of the LFG collection system.In this way, risk areas or points where collection systems are not operating correctly can be identified.In turn, the 3D model can be used to provide geometric information (X, Y, Z), which enables the analysis of the evolution of anomalies through dynamic models, as well as the modeling of fire propagation and the planning of corrective measures.The analysis based on 4D models enables a global conception of MSW as opposed to traditional 2D alternatives, such as heat maps or isolated TIR images, thus providing more complete information in four dimensions.
In short, the proposed methodology is a novel alternative for the implementation of UAVs and low-cost thermal sensors, combined with GIS tools, for thermal monitoring of landfills.The combination of these technologies enables a fast, cost-effective, accessible, safe, and simple assessment compared to conventional 2D techniques for the identification of thermal anomalies in MSW landfills.

Figure 1 .
Figure 1.Diagram of the proposed workflow for the 4D model generation and identification of landfill thermal anomalies.

Figure 1 .
Figure 1.Diagram of the proposed workflow for the 4D model generation and identification of landfill thermal anomalies.

Figure 2 .
Figure 2. General overview of the Meruelo MSW landfill and distribution of the different facilities: treatment and recovery plant, landfill, and leachate plant.

Figure 2 .
Figure 2. General overview of the Meruelo MSW landfill and distribution of the different facilities: treatment and recovery plant, landfill, and leachate plant.

Figure 3 .
Figure 3. Image of the DJI Mavic 3T UAV which was used for the flight over the landfill.

Figure 3 .
Figure 3. Image of the DJI Mavic 3T UAV which was used for the flight over the landfill.

Figure 4 .
Figure 4. Double grip UAV flight plan designed for the Meruelo MSW landfill.

Table 4 .
Definition of the parameters of the UAV flight plan over the landfill./px) 4.77 (RGB) 10.71 (IR) Number of photographs 1289 Flight time (min) 87

Figure 6 .
Figure 6.3D thermal model (X,Y,Z-temperature) of the MSW landfill.The 4D point cloud provides the position in space and the temperature for each point.

Figure 6 .
Figure 6.3D thermal model (X,Y,Z-temperature) of the MSW landfill.The 4D point cloud provides the position in space and the temperature for each point.

Figure 6 .
Figure 6.3D thermal model (X,Y,Z-temperature) of the MSW landfill.The 4D point cloud provides the position in space and the temperature for each point.

Figure 8 .
Figure 8. Hot spots identified through photointerpretation of the surface temperature distribution in the thermal orthoimage of the Meruelo MSW landfill.

Figure 8 .
Figure 8. Hot spots identified through photointerpretation of the surface temperature distribution in the thermal orthoimage of the Meruelo MSW landfill.

Figure 8 .
Figure 8. Hot spots identified through photointerpretation of the surface temperature distribution in the thermal orthoimage of the Meruelo MSW landfill.

Figure 9 .
Figure 9. Detail of the thermal and RGB photointerpretation for the hot spots A and D. These thermal anomalies are due to the exposure of geotextiles in the Meruelo MSW landfill.

Figure 9 .
Figure 9. Detail of the thermal and RGB photointerpretation for the hot spots A and D. These thermal anomalies are due to the exposure of geotextiles in the Meruelo MSW landfill.

Figure 10 .
Figure 10.Detail of the thermal and RGB photointerpretation of the hot spot G, which is associated with the leachate collection system.

Figure 11 .
Figure 11.Detail of the thermal and RGB photointerpretation of homogeneous terrain areas where the non-justified hot spots B and C are identified.

Figure 10 .
Figure 10.Detail of the thermal and RGB photointerpretation of the hot spot G, which is associated with the leachate collection system.

Figure 10 .
Figure 10.Detail of the thermal and RGB photointerpretation of the hot spot G, which is associated with the leachate collection system.

Figure 11 .
Figure 11.Detail of the thermal and RGB photointerpretation of homogeneous terrain areas where the non-justified hot spots B and C are identified.

Figure 11 .
Figure 11.Detail of the thermal and RGB photointerpretation of homogeneous terrain areas where the non-justified hot spots B and C are identified.Appl.Sci.2023, 13, 13164 15 of 22

Figure 12 .
Figure 12.Detail of the thermal and RGB photointerpretation of the non-justified hot spot due to stockpiling E.

Figure 12 .
Figure 12.Detail of the thermal and RGB photointerpretation of the non-justified hot spot due to stockpiling E.

Finally, Figure 13
shows a homogeneous zone (F) of high temperature (23.8 • C), similar to hot spots B and C, where the thermal difference cannot be justified by differences in the emissivity or heat capacity of the surface represented.

Figure 12 .
Figure 12.Detail of the thermal and RGB photointerpretation of the non-justified hot spot due to stockpiling E.

Figure 13 .
Figure 13.Detail of the thermal and RGB photointerpretation of non-justified hot spot F in the oldest area of the landfill.

Figure 13 .
Figure 13.Detail of the thermal and RGB photointerpretation of non-justified hot spot F in the oldest area of the landfill.

Figure 14 .
Figure 14.The overlapping of the identified hotspots and the LFG recovery well locations enables evaluation of the performance of the emissions control system.

Figure 14 .
Figure 14.The overlapping of the identified hotspots and the LFG recovery well locations enables evaluation of the performance of the emissions control system.