L'augmentation importante de la population mondiale, et par conséquent de ses besoins, exerce une pression de plus en plus importante sur les surfaces forestières. L'outil le mieux adapté au suivi des forêts, à l'échelle du globe, est la télédétection. C'est dans ce contexte que se situe ce travail de thèse, qui vise à améliorer l'estimation des paramètres biophysiques des arbres à partir de données de télédétection. L'originalité de ce travail a été d'étudier cette estimation des paramètres biophysiques en menant plusieurs études de sensibilité avec une démarche expérimentale sur des données expérimentales et sur des données simulées. Tout d'abord, l'étude s'est portée sur des séries temporelles de mesures de diffusiométrie radar obtenues sur deux sites : l'un constitué d'un cèdre en zone tempérée et l'autre d'une parcelle de forêt tropicale. Puis, cette étude de sensibilité a été poursuivie en imageant, avec une résolution élevée, plusieurs parcelles aux configurations différentes à l'intérieur d'une forêt de pin. Enfin, des données optiques et radars simulées ont été fusionnés afin d'évaluer l'apport de la fusion de données optique et radar dans l'inversion des paramètres biophysiques.

Titre de l'écran-titre (visionné le 25 mars 2024) / Au cours des dernières décennies, la télédétection a grandement facilité la documentation des tendances associées au raccourcissement de la période d'englacement des lacs en réponse au changement climatique. Bien que ces investigations soient diversifiées, certaines lacunes persistent, notamment l'exclusion des petits lacs et des réservoirs des études globales, et ce malgré leur abondance dans l'hémisphère nord. Ainsi, l'objectif de ce projet est de caractériser l'englacement des lacs et réservoirs à travers l'hémisphère nord grâce à l'imagerie satellitaire Sentinel-2A. Il vise également à examiner les tendances spatiales des dates de phénologie et à évaluer l'influence des caractéristiques géographiques, morphométriques et climatiques associées aux plans d'eau. Ces travaux reposent sur le développement d'une vaste base de données regroupant les dates de phénologie de la glace de 3702 lacs et de 1028 réservoirs à travers l'hémisphère nord de 2019 à 2023. En outre, cette base de données a permis la calibration d'un modèle prédictif de l'englacement des plans d'eau. Ce modèle, de type forêt aléatoire, a démontré une bonne performance avec une erreur absolue moyenne de 11 et 13 jours pour la prise et la fonte de la glace, respectivement. Cette marge d'erreur s'explique principalement par les causes d'incertitudes associées à la détection de glace par imagerie optique, notamment la gestion de la couverture nuageuse, mais aussi la détection de glace noire en climat aride. Parmi les variables d'entrée du modèle, la latitude et l'accumulation maximum des degrés-jours de gel constituent les paramètres déterminants dans la prédiction des dates de gel et de fonte de la glace. Bien que l'utilisation des capteurs optiques présente des défis spécifiques, le potentiel de la télédétection reste significatif dans le cadre d'études globales sur la glace, offrant également la possibilité d'améliorer les modèles climatiques actuels. L'étude de cas du réservoir de la Romaine-2 illustre un exemple concret de la synergie entre la télédétection et la modélisation. / In recent decades, remote sensing has considerably facilitated the documentation of trends related to the shortening of lake ice duration in response to climate change. Although these studies are diverse, certain gaps remain, in particular the exclusion of small lakes and reservoirs from global studies, despite their abundance in the Northern Hemisphere. Using Sentinel-2 satellite imagery, this project aims to characterize lake and reservoir freeze-up period across the Northern Hemisphere, and to assess the influence of their geographic, morphometric and climatic characteristics. This work is based on the development of an extensive database of ice phenology dates for 3702 lakes and 1028 reservoirs across the Northern Hemisphere from 2019 to 2023. This database was also used to calibrate a predictive model of lake and reservoir ice conditions. This random forest model performed well, with mean absolute errors of 11 and 13 days for freeze-up and break-up, respectively. This error is mainly due to uncertainties associated with ice detection from optical imagery, particularly cloud cover management, but also black ice detection in arid climates. Among the model input variables, latitude and maximum accumulation of freezing degree-days were the determining parameters in predicting freeze-up and break-up dates. Although the use of optical sensors presents specific challenges, the potential of remote sensing in the context of global ice studies remains significant and offers the possibility of improving current climate models. The case study of the Romaine-2 reservoir provides a concrete example of the synergy between remote sensing and modeling.

Titre de l'écran-titre (visionné le 18 mars 2024) / La fréquence et l'ampleur des dommages causés par les inondations au Canada augmentent chaque année, nécessitant des solutions préventives. Les étages inférieurs des bâtiments sont particulièrement vulnérables, justifiant la nécessité de détecter avec précision les ouvertures comme les portes et les fenêtres. Les nuages de points obtenus par télémétrie mobile sont bien adaptés pour identifier la présence de ces ouvertures sur les façades de bâtiments. Toutefois, la présence d'occlusions dans les nuages de points causées par des objets obstruant le faisceau LiDAR rend cette tâche complexe. Le présent projet de maîtrise vise à répondre à cette problématique en proposant une approche pour détecter et caractériser les occlusions dans les scènes acquises dans les milieux résidentiels et ruraux. La conception de la méthode est articulée autour de deux volets : la segmentation sémantique pour étiqueter le nuage de points et la détection d'occlusions par lancer de rayons. La méthode a été validée à partir de données simulées car elles permettent un contrôle total sur l'environnement et facilitent la visualisation des résultats. Cette analyse a ensuite été transposée sur les résultats obtenus à partir du nuage de points réel. Les résultats obtenus par la méthode proposée ont démontré une capacité à détecter la présence d'occlusions dans les milieux résidentiels et ruraux. La combinaison de la segmentation sémantique pour l'étiquetage du nuage de points avec la détection d'occlusions par lancers de rayons a permis une identification robuste des occlusions. Ces résultats soulignent l'efficacité de la solution dans des contextes réalistes, renforçant ainsi sa pertinence pour la détection précise des ouvertures sur les façades, une étape cruciale pour évaluer les risques liés aux inondations. La transposition réussie de l'analyse des données simulées aux résultats du nuage de points réel valide la robustesse de l'approche et son applicabilité dans des scénarios du monde réel. / The frequency and extent of flood damage in Canada increases every year, necessitating preventive solutions. The lower floors of buildings are particularly vulnerable, justifying the need for accurate detection of openings such as doors and windows. Point clouds obtained by mobile telemetry are suitable for identifying the presence of such openings on building facades. However, the presence of occlusions in the point clouds caused by objects obstructing the LiDAR beam makes this a complex task. The present Master's project aims to address this issue by proposing an approach for detecting and characterizing occlusions in scenes acquired in residential and rural environments. The method was designed around two components: semantic segmentation to label the point cloud, and occlusion detection using ray tracing. The method was validated using simulated data, as these allow total control over the environment and facilitate visualization of the results. This analysis was then transposed to results obtained from the real point cloud. The results obtained by the proposed method demonstrated an ability to detect the presence of occlusions in residential and rural environments. Combining semantic segmentation for point cloud labeling with occlusion detection by ray-tracing enabled robust occlusion identification. These results underline the effectiveness of the solution in realistic contexts, reinforcing its relevance for the accurate detection of facade openings, a crucial step in assessing flood risks. The successful transposition of simulated data analysis to real point cloud results validates the robustness of the approach and its applicability in real-world scenarios.

Le nerprun bourdaine (Rhamnus frangula L.) est une espèce exotique qui envahit plusieurs régions du sud du Québec, et plus particulièrement la région administrative de l'Estrie. Actuellement, on connaît encore peu l'écologie de l'espèce dans le contexte québécois et il n’existe pas de portrait d’ensemble de sa distribution dans les forêts tempérées de cette région. Dans ce contexte, le premier objectif du projet était de cartographier par télédétection la distribution du nerprun bourdaine dans deux secteurs de l'Estrie. Un second objectif était d'évaluer les variables environnementales déterminantes pour expliquer le recouvrement de nerprun bourdaine. La phénologie du nerprun bourdaine diffère de celle de la plupart des espèces indigènes arborescentes puisque ses feuilles tombent plus tard en automne. Cette caractéristique a permis de cartographier, par démixage spectral, la probabilité d'occurrence du nerprun bourdaine grâce à une série temporelle d'images du capteur OLI de Landsat 8. Le recouvrement du nerprun bourdaine a été calculé dans 119 placettes sur le terrain. La cartographie résultante a montré un accord de 69% avec les données terrain. Une image SPOT-7, dont la résolution spatiale est plus fine, a ensuite été utilisée, mais n’a pas permis d'améliorer la cartographie, puisque la date d’acquisition de l’image n’était pas optimale dû à un manque de disponibilité. Concernant le second objectif de la recherche, la variable la plus significative pour expliquer la présence de nerprun bourdaine était la densité du peuplement, ce qui suggère que l’ouverture de la couverture forestière pourrait favoriser l’envahissement. Néanmoins, les résultats tendent à démontrer que le nerprun bourdaine est une espèce «généraliste» qui s’adapte bien à plusieurs conditions environnementales. / Glossy buckthorn (Rhamnus frangula L.) is an exotic species invading many areas in southern of Quebec, particularly in the Eastern Townships. Currently, we do not know very much about the species ecology and no thorough study of its distribution in temperate forest has been performed. Therefore, the first objective of the project was to map the spatial distribution of glossy buckthorn in two areas of the Eastern Townships, using remote sensing techniques. The second objective was to evaluate the environmental variables, or predictors, best explaining the presence of glossy buckthorn. The phenology of glossy buckthorn differs from most of the indigenous tree species found in this area because its leaves fall later in autumn. This characteristic allowed to map, using spectral unmixing, the probability of occurrence of glossy buckthorn, with temporal Landsat 8 (OLI) imagery data series. Glossy buckthorn coverage was calculated on 119 plots on the field. The resulting maps showed an agreement of 69% with field data. A SPOT-7 image, which has a finer resolution than Landsat 8 (OLI), was then used but it did not improve the quality of the map, since its acquisition date was not optimal, due to a lack of availability. Concerning the second objective of the research, the best variable explaining the presence of glossy buckthorn was stand density, which leads to believe that forest cover openings could ease the establishment of buckthorn. However, the results tend to show that glossy buckthorn is a generalist species, easily adapting to various environmental conditions.

L'augmentation importante de la population mondiale, et par conséquent de ses besoins, exerce une pression de plus en plus importante sur les surfaces forestières. L'outil le mieux adapté au suivi des forêts, à l'échelle du globe, est la télédétection. C'est dans ce contexte que se situe ce travail de thèse, qui vise à améliorer l'estimation des paramètres biophysiques des arbres à partir de données de télédétection. L'originalité de ce travail a été d'étudier cette estimation des paramètres biophysiques en menant plusieurs études de sensibilité avec une démarche expérimentale sur des données expérimentales et sur des données simulées. Tout d'abord, l'étude s'est portée sur des séries temporelles de mesures de diffusiométrie radar obtenues sur deux sites : l'un constitué d'un cèdre en zone tempérée et l'autre d'une parcelle de forêt tropicale. Puis, cette étude de sensibilité a été poursuivie en imageant, avec une résolution élevée, plusieurs parcelles aux configurations différentes à l'intérieur d'une forêt de pin. Enfin, des données optiques et radars simulées ont été fusionnés afin d'évaluer l'apport de la fusion de données optique et radar dans l'inversion des paramètres biophysiques. / The significant increase of the world population, and therefore its needs, pushes increasingly high in forest areas. The best tool for monitoring forest across the globe is remote sensing. It is in this context that this thesis, which aims to improve the retrieval of biophysical parameters of trees from remote sensing data, takes place. The originality of this work was to study the estimation of biophysical parameters across multiple sensitivity studies on experimental data and simulated data. First, the study focused on the time series of radar scatterometry measurements obtained on two sites: one characterized by a cedar in the temperate zone and the other by a forest plot of rainforest. Then, the sensitivity analysis was continued by imaging with high resolution, several forest plots with different configurations within a pine forest. Finally, simulated radar and optical data were combined to evaluate the contribution of optical and radar data fusion in the inversion of biophysical parameters.

Avec son échantillonnage spatial élevé, IASI (Interféromètre de sondage atmosphérique dans l'infrarouge) embarqué sur la plateforme MetOp, fournit 1.300.000 spectres par jour et une couverture globale bi-journalière. Cette spécificité du sondeur, qui offre la possibilité de suivre l'évolution des concentrations des espèces présentes en faible quantité dans l'atmosphère, est déterminante pour notre étude de la chimie dans les panaches de feux. Elle s'appuie<p>sur les premières observations des feux par IASI, qui ont révélé la possibilité de mesurer, en plus du CO habituellement utilisé comme traceur, plusieurs espèces chimiquement actives dans les panaches, dont NH3 et divers composés organiques. Via une analyse approfondie des spectres, nous montrons dans notre travail qu'au total 24 espèces gazeuses sont détectées dans les spectres IASI, ce qui représente une avancée majeure dans la détection par les sondeurs infrarouges au nadir.<p>Ces observations inattendues ouvrent la voie à un examen approfondi de la réactivité dans les panaches, au cours de leur transport et en fonction des conditions météorologiques locales. Nous nous sommes attachés à cette étude en ciblant principalement deux événements intenses et isolés. Le premier événement est celui des feux d'Australie de Février 2009, qui était particulier de par la hauteur d'injection du panache (haute troposphère voire basse stratosphère) mais aussi suite à la forte quantité d'aérosols émise. Le second événement est celui des feux de Russie centrale de l'été 2010, qui se caractérisent par la durée de l'épisode (d'environ un mois) mais aussi par sa très large étendue.<p>Nous avons déterminé, pour le premier événement, les masses totales d'une série d'espèces émises dans les panaches et nous les avons confrontées à celles issues du modèle de chimie-transport CHIMERE et à d'autres valeurs de la littérature. Nous montrons une évolution temporelle fortement différente pour les espèces plus réactives et supposons pour certains composés, dont HCOOH,<p>une vraisemblable source secondaire au sein du panache.<p>Pour les feux en Russie, nous avons déterminé les rapports d'enrichissement ∆X/∆CO (où X est l'espèce réactive), particulièrement utiles pour suivre la<p>chimie. Nous les comparons aux valeurs publiées dans la littérature en fonction du type de végétation brûlée. Certains aspects de la réactivité dans le panache sont mis en évidence, notamment la formation secondaire de HCOOH. Nous<p>avons également estimé les masses totales et les flux d'émission de chacune des 3 espèces émises (CO, NH3, et HCOOH), et les comparons à d'autres<p>estimations publiées dans la littérature scientifique.<p>Sur base de ces résultats, pour des événements locaux, nous menons une brève première analyse de la variabilité des émissions de gaz réactifs par les feux à l'échelle globale, en utilisant plus particulièrement les distributions de concentration de CO, NH3, et HCOOH. Nos résultats pour HCOOH/CO suggèrent une grande variabilité des rapports tant spatialement que temporellement. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Les traits fonctionnels foliaires sont des paramètres clés des processus écologiques dans les forêts. Malgré les progrès réalisés dans la cartographie de traits foliaires à partir d’imagerie hyperspectrale, il subsiste un besoin de construire des modèles spectres-traits spécifiques à un environnement donné, entraînés à partir de mesures à l’échelle de l’arbre, afin d'améliorer la précision des cartes locales de traits. Nous avons cartographié 12 traits fonctionnels foliaires dans une forêt tempérée mixte à partir d’imagerie hyperspectrale aéroportée. Des échantillons foliaires de cimes d'arbres (n = 166), couvrant un total de 16 espèces, ont été recueillis à l'aide d'une plateforme de drone pour mesurer des traits foliaires d’arbres individuels, à partir desquels la réflectance spectrale de la couronne a également été mesurée. Des modèles de régression des moindres carrés partiels (PLSR) ont été utilisés pour prédire les traits foliaires à partir des spectres de réflectance à l’échelle de l’arbre (400-2400 nm). Ces modèles ont prédit la masse foliaire par unité de surface (LMA), la surface foliaire spécifique (SLA) et l'épaisseur en eau (EWT) avec une haute précision (R2 > 0.8, %RMSE < 15). Les modèles de concentration de pigments, de l'azote et de la cellulose ont montré une performance modérée (R2 = 0.53–0.68, %RMSE = 17.24–21.31). Les performances les plus faibles ont été observées pour la lignine, le carbone, le contenu en matière sèche (LDMC) et l'hémicellulose (R2 = 0.24–0.44, %RMSE = 20.67–26.13). Des cartes à haute résolution spatiale (1.25 m pixel-1) de traits foliaires ont été produites pour l'ensemble de la zone d'étude de 16 km2. Notre étude s'ajoute aux recherches approfondies visant à utiliser la télédétection pour surveiller la biodiversité des traits fonctionnels à plus grande échelle et fournit des modèles qui saisissent la variation intraspécifique de nombreuses espèces d'arbres d'une forêt tempérée mixte de l'est du Canada. / Foliar functional traits are key drivers of ecological processes in forests. Despite progress in forest foliar trait mapping from imaging spectroscopy, there is a need to build environment-specific, spectra-trait models trained from tree-level measurements to improve the accuracy of local trait maps. We mapped 12 foliar functional traits in a mixed temperate forest using airborne imaging spectroscopy. Top-of-canopy foliar samples from tree crowns (n = 166), representing a total of 16 species, were collected using a drone platform to measure foliar traits for individual trees, from which tree-level crown spectra were also determined. Partial least squares regression (PLSR) models were used to predict foliar traits from tree-level reflectance spectra (400-2400 nm). These models predicted leaf mass per area (LMA), specific leaf area (SLA) and equivalent water thickness (EWT) with high accuracy (R2 > 0.8, %RMSE < 15). Models for pigment, nitrogen and cellulose concentrations showed a moderate performance (R2 = 0.53–0.68, %RMSE = 17.24–21.31). Poorest performance was observed for lignin, carbon, leaf dry mass content (LDMC) and hemicellulose (R2 = 0.24–0.44, %RMSE = 20.67–26.13). High-resolution (1.25 m pixel-1) foliar trait maps were produced for the entire 16-km2 study area. Our study adds to the extensive research aiming to use remote sensing to monitor forest functional trait biodiversity at larger scales and provides models that capture intraspecific variation across many tree species from a mixed temperate forest in eastern Canada.

Avec l'augmentation des menaces asymétriques partout sur la planète depuis le début des années 1990s impliquant des armes produisant de forts dommages par rapport à leurs coûts, plusieurs types d'armes biologiques sont devenus une menace d'importance croissante dans le théâtre militaire ainsi que dans le monde civil. Les armes biologiques peuvent êtres facilement dissimulées sur une personne. Elles sont difficilement détectables et leurs déploiements peuvent causer des dommages psychologiques incommensurables sur les populations affectées, assmant un ravage souvent loin au delà de leur emplacement physique. Pour essayer de diminuer l'efficacité de ces armes, des méthodes pour détecter et distinguer en retrait les agents biologiques parmi les autres aérosols naturels présents dans l'atmosphère sont actuellement en développement. Recherche et Développement pour la Défense Canada (RDDC) et le gouvernement du Canada explorent actuellement les LIDARs comme la technologie principale visant la détection à distance d'organismes biologiques. Toutes informations pouvant êtres récoltées à propos des armes biologiques en utilisant cette technologie sont utiles à l'avancement de la cause de la biodétection à distance. Ce rapport présente les résultats d'un projet ayant pom but de décrire les effets de l'âge des bactéries sur leurs propriétés spectrales de fluorescence en utilisant la technologie LIDAR. Les expériences ont été exécutées dans une chambre hermétique LIDAR de laboratoire sur des simulants d'agents biologiques connues, trois desquelles (2 souches de Bacillus subtilis var. niger identifiées comme BGnew et BGold et une souche de Bacillus thuringiensis identifié comme BT) étant des simulants de la bactérie Bacillus anthracis (l'anthrax), sans doute l'agent biologique le plus probable à être utilisé, et l'autre (Erwinia herbicola identifié comme EH) étant un simulant de Yersinia pestis (la peste). Les bactéries étaient mélangées avec leur milieu de culture, soit le ± Tryptic soy broth ¿ (TSB). Bien que plus de tests soient nécessaires afin de pouvoir finaliser des conclusions plus robustes, certaines tendances ont été observées. Le mélange de Erwinia herbicola a montré une fluorescence induite ayant une tendance vers les courtes longueurs d'ondes lorsque la culture était plus âgée. Le mélange de Bacillus thuringiensis (BT) a démontré une petite tendance vers les longues longueurs d'ondes avec l'âge. Les deux différentes souches de Bacillus globigii ont montrées des fluorescences induites ayant des tendances vers les courtes longueurs d'ondes après un mois et vers les longues longueurs d'ondes après trois mois. En bref, la signature de fluorescence semble dépendre de l'âge du mélange bactérie-bouillon de culture bien que cet effet varie en fonction du type mélange. Certains mélanges étaient plus efficaces à produire de la fluorescence que d'autres. Le mélange de BGnew-ljour a démontré la plus haute efficacité de fluorescence tandis que BG_old-3mois, BGold-ljour et EH-ljour ont démontrés les plus faibles. Nos résultats démontrent que la signature spectrale d'une bactérie est influencée par le temps d'incubation dans son milieu de culture. Cette information est très importante lors de l'élaboration de banques de signatures standardisées servant à la classification bactérienne à distance dans le cadre d'intervention sur le terrain. Dans des contextes militaires et de sécurité publique, cette information serait aussi un outil important pour retrouver l'origine des agents biologiques (ex: s'ils parviennent d'un labo ou d'un lieu d'entreposage).

En passant de 50 000 à plus d'un million d'individus en moins de 40 ans, la Grande Oie des neiges exerce une pression si forte sur son habitat que l'équilibre écologique de ce dernier est menacé.En effet, lors de ses haltes migratoires automnales et printanières, elle dévore tout sur son passage, en particulier les rhizomes de scirpe américain. Cette situation est fort préoccupante, particulièrement pour la préservation de l'intégrité écologique des marais. Il s'avère donc essentiel de déterminer avec le plus de précision possible à quel point l'évolution de la végétation des marais à scirpe est importante. Pour ce faire, une analyse spatio-temporelle a été effectuée pour les sites traditionnels d'halte migratoire, c'est-à-dire Cap Tourmente, Montmagny, l'Isle-aux-Grues et Cap-Saint-Ignace. Aux images satellitaires IKONOS datant de 2002 s'ajoutent des photographies aériennes prises au cours des quatre dernières décennies. Des classifications servant à identifier les groupements végétaux et les taux de recouvrement du scirpe sont réalisées à partir du logiciel eCognition ( Definiens ) afin de tirer profit de toutes les caractéristiques spectrales, texturales et contextuelles des classes identifiées. De nombreuses données prises sur le terrain servent à la fois à l'entraînement et à la validation des classifications. Les résultats tirés des analyses des secteurs de Cap-St-Ignace, Cap Tourmente et Montmagny révèlent une diminution des proportions de scirpe particulièrement importante à partir des années 1980 à l'intérieur du bas marais. Cette transition s'effectue généralement au profit de la zizanie aquatique dont les populations se sont sensiblement accrues au cours de la même période. De plus, un phénomène important d'érosion a été observé à la limite entre le haut et le bas marais pour pratiquement l'ensemble des secteurs étudiés. Le taux d'érosion peut atteindre 1,5 m/année à plusieurs endroits, ce qui met en péril l'équilibre écologique de ces milieux.

Proliferation of cyanobacteria is a growing problem in lacustrine environment that results in rapid degradation of water quality. Moreover, certain cyanobacteria species produce harmful toxins. Phycocyanin (PC) is a photosynthetic pigment typical of cyanobacteria and affects the water color: it is therefore possible to study them using remote sensing. At least three algorithms to estimate PC concentration ([PC]) have been published, but their relative errors are important, especially for lower concentration. In this study, we are presenting the results of a new algorithm that uses the second order variability (anomalies) of water's reflectance spectrum to estimate [PC]. This method has never been used in lacustrine environment. The dataset used to develop and validate the algorithm was obtained between 2001 and 2005 in 57 different lakes and reservoirs of the Netherlands and Spain. The performance of the second order algorithm is equivalent or better than the three previously published algorithms. For the subset were [PC] > 32 mg m[superscript -3], the contribution of the second order term (R[superscript 2] =0.68 and RMSE=0.25) seems to improve considerably the first order algorithm (R[superscript 2] =0.50 and RMSE=0.35). The accuracy of the second order algorithm for [PC] > 32 mg m[superscript -3] is superior to the one calculated for the whole dataset (R[superscript 2] =0.69 and RMSE=0.44). The algorithm can also be adapted to the bands of satellite sensor MERIS for the study of cyanobacteria. The application of this algorithm to a MERIS image acquired the 29 August 2010 taken over the Missisquoi Bay (Quebec, Canada) demonstrates the potential of this new algorithm for a future cyanobacteria' monitoring system. Note that all the statistical results presented above are for the logarithm of [PC] and the units of the RMSE are log(mg/m[superscript 3]).