FR 2953021 A1 20110527 FR 0958385 A 20091126 Domaine technique L'invention concerne le domaine des méthodes de mesure des émissions de gaz à effet de serre (GES). L'invention porte en particulier sur une méthode de mesure des émissions hebdomadaires et annuelles d'un gaz à effet de serre sur une aire géographique donnée. Elle concerne encore un système de mesure permettant la mise en oeuvre de la méthode de mesure. Art antérieur Les GES sont les gaz de l'atmosphère naturels et anthropiques, qui absorbent et émettent des rayonnements à des longueurs d'onde du spectre de rayonnement infrarouge thermique. Ce sont principalement le dioxyde de carbone CO2, le méthane CH4, le protoxyde d'azote N2O, Les oxydes d'azote NOx, les hydrofluorocarbures HFC, les hydrochlorofluorocarbures HCFC, les chlorofluorocarbures CFC, les perfluorocarbures PFC, l'hexafluorure de soufre SF6, l'ozone troposphérique 03, la vapeur d'eau H2O, le monoxyde de carbone CO et le dihydrogène H2. Le CO2 est généralement le gaz de référence. Lorsqu'ils absorbent le rayonnement infrarouge thermique émis par la surface de la terre, par l'atmosphère et les nuages, le rayonnement atmosphérique est émis de tous les côtés, jusqu'à la surface de la terre. Les GES diffèrent dans leur forçage radiatif sur le système climatique en raison de leurs différentes propriétés radiatives et de leur durée de vie dans l'atmosphère. Ils retiennent la chaleur au sein de la surface-troposphère, ce qui est communément appelé « effet de serre » et une augmentation de leur concentration peut engendrer un renforcement de cet effet avec un réchauffement. Les sources naturelles de CO2 sont beaucoup plus importantes que les sources anthropiques, mais sur de longues périodes de temps, les sources naturelles sont équilibrées par les puits 1 naturels. La concentration atmosphérique de CO2 est restée entre 260 et 280 parties par million (ppm) dans l'atmosphère depuis l'holocène, c'est-à-dire les derniers 10'000 ans interglaciaires, mais, depuis l'ère industrielle, l'activité humaine a augmenté sa concentration de l'ordre de 100ppm. La communauté scientifique a récemment reconnu que l'effet de serre issu des GES anthropiques a produit un forçage positif sur la température de la surface d'environ 1°C supérieur à la moyenne depuis le milieu du 20e siècle. Il est donc probable que le réchauffement anthropique dû à une élévation des niveaux de GES ait influencé les systèmes physiques et biologiques naturels. Les changements attendus sur les facteurs climatiques sont notamment d'impacter les ressources en eau douce, l'industrie, l'alimentation et la santé. La stabilisation des concentrations à un niveau qui empêcherait toute perturbation anthropique dangereuse du système climatique est donc devenue une priorité pour la communauté internationale. Le CO2 est la forme de carbone la plus répandue dans l'atmosphère, et la principale source de carbone dans les matières organiques. Il est issu des échanges entre l'atmosphère et la biosphère, l'atmosphère et les océans, les perturbations de la biosphère et la production anthropique. Le déséquilibre entre l'absorption et le dégagement conduit à une augmentation nette dans l'atmosphère. Dans la méthode d'évaluation appelée « CarbonTracker », la loi de conservation de masse est utilisée pour évaluer le flux atmosphérique FCO2(t) en assumant que la masse de carbone dans l'atmosphère est égale à l'effet net de tous les sources et les puits à un moment donné t. Ce flux signifie à la fois l'échange de masse par unité de surface et également 2 l'échange de masse sur une aire intégrée dans le contexte de superficies finies. On a donc : Fco2 (t) = Face (t) + Fbia (t) + Fff(t) + Ffeu(t) est le flux accumulé atmosphérique net de CO2 est le flux net des échanges atmosphère/océan est le flux net des échanges atmosphère/biosphère est le flux net des sources anthropiques est le flux net des sources liées aux incendies Le flux net d'accumulation atmosphérique de CO2 dans l'atmosphère est généralement exprimé en Petagrammes de Carbone par an (PgC/an) ou GTCO2 par an aux échelles régionales. A l'échelle d'une installation humaine, le flux est exprimé en TCO2/an ou en TCO2 équivalent/an en additionnant aux émissions de CO2, celles des autres GES en fonction de leur potentiel de réchauffement global comparé au CO2. En raison de sa longue espérance de vie, les concentrations atmosphériques de CO2 sont estimées assez uniformes et leur variation contribue à estimer les échanges de flux. Pendant l'holocène, les concentrations indiquées par les analyses de carottes de glace de Vostok et Taylor Dome étaient d'environ 275ppm, très inférieures à celles actuelles. En 2007 et 2008, la moyenne de concentration de CO2 dans l'atmosphère était respectivement d'environ 383.71ppm et 385.57ppm selon les observations du Mauna Loa. Avec une masse molaire de l'air d'environ 28.84g.mol-1 et une masse de l'atmosphère d'environ 5.137xlO18kg, lppm de CO2 représente environ 2.137PgC, ce qui permet de calculer un flux global annuel net 2008/2007 d'accumulation atmosphérique provenant des échanges naturels et anthropiques de: dCcoz(tzooa/zoos) Mc matin 3.97PgC ( 15 GTCO2 ) dt2008/zoo? Mair 106 FCOZZOOS/zoos 3 Le Protocole de Kyoto de la CCNUCC (Convention-Cadre des nations Unies sur les Changements Climatiques) est entré en vigueur le 16 février 2005 et cet accord international définit des objectifs astreignant 37 pays industrialisés et la Communauté Européenne à réduire les émissions de CO2, CH4, N2O, HFCs, PFCs et SF6 au moins de 5% au-dessous des niveaux de 1990 dans la période d'engagement 2008-2012. Les pays doivent répondre à leurs objectifs principalement par le biais de mesures nationales et ils se concentrent sur la diminution de la demande en biens et services intensifs en émissions, le développement de technologies en faible teneur de carbone, l'accroissement de leur efficacité énergétique et la réduction d'usage de combustibles fossiles. Pour vérifier leurs objectifs, les émissions sont surveillées et un reporting est effectué par les pays en soumettant des inventaires annuels d'émissions. Ces inventaires nationaux sont une liste détaillée des estimations des sources et des puits de GES nationaux par secteur et ils servent à la fois de base à la mise en place de mesures d'atténuation efficaces et à s'assurer que les tendances sont conformes aux engagements. Le protocole offre également des moyens supplémentaires aux pays pour atteindre ces objectifs tels que les marchés d'échanges de quotas d'émissions (Ex: European Union Emission Trading Scheme). Dans ces mécanismes, une autorité centrale définit un quota pour toutes les sources et les permis d'émission sont ensuite livrés aux installations, qui s'engagent à tenir un nombre équivalent de permis (ou crédits), qui représentent le droit d'émettre un montant précis. Le montant total des émissions ne peut dépasser le quota, limitant ainsi les émissions à ce niveau. Les installations sont autorisées à acheter et vendre leurs permis entre elles, ceci visant à stimuler l'investissement 4 écologique et à réduire leurs niveaux avec le meilleur rapport coût-efficacité. Afin que ces nouveaux marchés soient efficients, les participants doivent avoir confiance dans l'exactitude des données reportées pour établir les niveaux de référence d'émissions. Parallèlement, les autorités réglementaires ont également pour souci que les méthodes de surveillance, de reporting et de vérification (MRV, Monitoring, Reporting and Verification) des inventaires de GES par installation aient un haut degré de fiabilité afin de pouvoir démontrer leur conformité. Ces méthodes MRV sont essentiellement effectuées par le biais de calculs et d'observations « ascendantes » ou « bottom-up » des émissions de GES par installation. Malgré leurs améliorations continuelles, des incertitudes significatives restent, en particulier sur certaines catégories de sources ou les facteurs d'émissions peuvent être assez variables et les processus de mesure peuvent manquer d'homogénéité d'une installation à l'autre. Un objectif d'incertitude sur les résultats est d'environ 5% selon les standards internationaux et industriels et il est reconnu que les installations peuvent rencontrer certaines difficultés à atteindre ce niveau d'excellence. Les méthodes de mesure ascendantes des inventaires sont généralement effectuées via des méthodes de calcul et d'observation sur chaque installation. Les méthodes de calcul permettent de déterminer les sources d'émissions en utilisant des données d'activité et sont obtenues en combinant systèmes de mesure et paramètres venant d'analyses de laboratoire ou de coefficients types sous la forme suivante: TCO2 CO2 émissions ( ) = données d' activité • facteur d'émission • facteur d'oxydation an Pour les émissions de combustion, les données d'activité sont fondées sur la consommation de carburant. La quantité de 5 carburant utilisée est généralement exprimée en termes de contenu énergétique et de facteur d'émission. Lorsqu'un carburant est consommé, seule une partie du carbone est oxydée en CO2, ce qui est pris en compte dans le facteur d'oxydation. CO2 émissions (TCO2) ll an = flux de carburant (T ou Nm3) • valeur nette calorifique (TI ou Nm3) /TCO2 • facteur d'émission I TI ) • facteur d'oxydation Pour les émissions de processus, les données d'activité sont fondées sur la consommation de matières, le débit ou la production et le facteur d'émission. Le carbone contenu dans les matières premières, non converti en CO2 est pris en compte dans le facteur de conversion. TCO2 CO2 émissions ( ) an (TCO2 TCO2) = données d'activité (T ou Nm3) • facteur d'émission T ou Nm3 ) • facteur de conversion Des logiciels de calcul basés sur l'activité d'une installation se développent actuellement sur ce même principe. D'autres méthodes de mesure ascendantes des émissions par installation déterminent, quant à elles, les émissions à leur source par la mesure en continu, en un point représentatif, des concentrations de GES dans les conduits et la fumée des gaz. Le flux de gaz Qe est calculé au moyen d'un bilan massique tenant compte des matières entrantes, du débit d'air à l'entrée, de l'efficacité du processus et à la sortie de la production et des concentrations de GES. heures opérations par an TCO2 CO2 émissions ( ) = an CO2concentrationi • Qei 6 Ces calculs et observations ascendants doivent cependant être effectués sur chaque installation et malgré leur amélioration constante, les résultats démontrent une certaine hétérogénéité d'une installation à l'autre en fonction des paramètres et processus utilisés. Ils se focalisent sur des points précis de l'installation et peuvent également omettre des sources adjacentes. Un deuxième type de méthode de mesure des émissions, appelé descendante, se focalise sur la compréhension du cycle du carbone afin de déterminer les sources et puits de carbone à différentes échelles géographiques et de calculer, en agrégeant les flux, les inventaires locaux. La fraction molaire de CO2 (ppm), définie comme le nombre de moles de CO2 divisée par le nombre total de moles (excepté l'eau) d'une parcelle d'air donnée, est communément utilisée car c'est une quantité conservative, qui ne dépend pas de la pression, de la température, de la vapeur d'eau ou de la teneur en eau condensée, qui sont très variables. Moins variable, elle dépend uniquement des échanges entre sources et puits de CO2, presque tous causés par les processus de surface. Elle reflète la somme des échanges de CO2 et forme le résultat ultime des influences humaines et naturelles combinées. Dans cette approche, le CarbonTracker est une référence internationale utilisée pour mieux comprendre la variabilité du cycle naturel du carbone et estimer les contributions naturelles et humaines. Il estime les échanges atmosphériques de CO2 en combinant modélisation et observation et son principe est similaire à d'autres systèmes d'assimilation de données. Il commence par prévoir les fractions molaires atmosphériques de CO2 sur le globe par une combinaison de modèles d'échanges (module océanique, module biosphère, module incendie et module fossile) avec un modèle de 7 transport atmosphérique, piloté par des prévisions météorologiques. La distribution du CO2 en 3D est ensuite échantillonnée en temps et lieu ou les observations sont disponibles et la différence entre les observations et les prévisions du modèle sont minimisées avec un ensemble d'assimilation de données. Cette minimisation est obtenue par paramétrage d'une série de facteurs de réduction linéaire qui contrôlent l'ampleur des flux de surface pour obtenir des flux finaux optimisés de résolution l°xl° sur l'Amérique du Nord et l'Europe. Une mesure « descendante » ou « top-down », précise, in-situ et indépendante des sources et des puits naturels de GES évaluant de l'échelle planétaire jusque localement les inventaires de GES des installations peut compléter et corréler les méthodes ascendantes actuelles. Elle peut aider à confirmer que les actions actuelles d'atténuation, entreprises par les pays et installations, réduisent efficacement les niveaux et renforcer la confiance et la crédibilité dans les marchés d'émission et dans la valeur des droits à polluer lorsque dans le contexte actuel, le prix du carbone reste relativement volatil et les niveaux de concentrations sont historiquement élevés. Résumé de l'invention Un premier but de l'invention est de fournir une méthode de mesure des inventaires nets de GES par aire géographique et/ou installation, corrigés des interférences avec les aires avoisinantes ou lointaines. La présente invention propose une méthode de mesure améliorée par rapport au CarbonTracker pour évaluer avec précision les inventaires de GES par aire géographique, en particulier par aire géographique représentative d'une installation anthropique, notamment les inventaires de CO2, mais aussi CH4, N2O, NOx, HCFC, HFC, CFC, 8 PFC, SF6, 03, H2O, CO et H2. La méthode est initialement présentée sur le CO2 et le même processus est utilisé pour les autres GES. Un premier avantage de l'invention est de fournir une méthode de mesure des émissions de gaz à effet de serre plus précise que les méthodes de mesure actuelles, notamment plus précise que le Carbontracker. En particulier, une première caractéristique de l'invention est de résoudre de plus petites échelles spatiales pour obtenir des flux nets anthropiques, notamment de CO2, en Kg/m2/s mesurés depuis l'échelle mondiale jusqu'au niveau des installations émettrices avec une résolution de 0.l°x0.1° (~l00km2). Avec la mesure des flux des autres GES, un but de la méthode de mesure selon l'invention est donc de déterminer, sur des surfaces pouvant être comprises entre lkm2 et 10000km2, les émissions en (TCO2/an), (TGES/an) et (TCO2eq/an) avec une précision supérieure à 5%. Avec une telle précision, cette mesure descendante, uniforme et globale des inventaires de GES a les avantages de comparer les résultats sur l'ensemble des installations pour vérifier les incohérences, éviter les omissions de sources, réduire l'incertitude par installation et compléter les évaluations ascendantes actuelles. En résumé, un premier but de la méthode de mesure selon l'invention est donc de fournir une méthode de mesure permettant de mesurer, du haut vers le bas, de l'échelle planétaire jusqu'au niveau de l'installation, les inventaires de GES afin de fournir une mesure précise des émissions. Un deuxième but de l'invention est de fournir un système de mesure des émissions de GES qui puisse être combiné, notamment avec les systèmes de gestion de production des installations, pour permettre le pilotage d'installations afin de limiter les émissions de combustion et/ou de processus et d'automatiser leur réduction. Des moyens 9 matériels et logiciels spécifiques implémentant la méthode de mesure selon l'invention et assurant l'interface avec des systèmes de gestion de production sont installés au sein d'installations émettrices en fonction de leur activité (énergie, procédés industriel, usages de produit...), des processus mis en ouvres par celles-ci et des GES émis. Le système de mesure selon l'invention permet alors de calibrer et d'optimiser le processus de chaque installation en fonction des niveaux et types d'émissions mesurés (ex : pics de pollution). Ceci permet d'obtenir une réduction d'émission automatisée sur chaque installation, de contrôler au fur et à mesure du temps son efficacité et de rester en conformité avec les normes régulatoires et environnementales. Un deuxième but de l'invention est donc de fournir un système de mesure qui puisse être interfacé directement au sein d'une installation émettrice afin d'optimiser les processus de production tout en mettant en place divers processus ou techniques de réduction, ceci permettant alors de calibrer les installations en contrôlant et en optimisant les niveaux d'émissions en fonction des mesures effectuées. Selon l'invention, la méthode de mesure des émissions hebdomadaires et annuelles d'un gaz à effet de serre générées sur une aire géographique déterminée comprend les étapes suivantes . - effectuer des mesures journalières de concentration dudit gaz à effet de serre en une première pluralité de lieux répartis sur la totalité du globe terrestre et enregistrer lesdites mesures journalières de concentration dans un module d'observation, - effectuer des mesures journalières des flux dudit gaz à effet de serre en une deuxième pluralité de lieux répartis sur la totalité du globe, et enregistrer i0 lesdites mesures journalières des flux dans ledit module d'observation, - effectuer des mesures de paramètres satellites, de paramètres météorologiques, de paramètres marins et de paramètres d'écosystèmes en une troisième pluralité de lieux répartis sur le globe terrestre et enregistrer les dites mesures de paramètres dans ledit module d'observation, - extraire, au moyen d'un module d'extraction, des données de prévisions météorologiques d'au moins une source de données, - effectuer une modélisation de l'évolution des flux dudit gaz sur le globe au moyen d'un module d'échange modélisant les sources et les puits naturels et anthropiques, - effectuer une modélisation des émissions anthropiques hebdomadaires dudit gaz à effet de serre au moyen d'un module d'inventaires ascendants, ledit module intégrant des données brutes d'émissions pour une pluralité d'installations, - effectuer, en utilisant ladite modélisation de l'évolution des flux, ladite modélisation des émissions anthropiques hebdomadaires et les dites données de prévisions météorologiques, une modélisation du transport atmosphérique dudit gaz à effet de serre au moyen d'un module de transport, - calculer les flux finaux dudit gaz à effet de serre, au moyen d'un module d'inversion et d'assimilation de données utilisant ladite modélisation des flux effectuée par le module d'échange, ladite modélisation des émissions anthropiques hebdomadaires effectuée par le module d'inventaires ascendants, ladite modélisation du transport atmosphérique 11 effectuée par le module de transport et lesdites mesures enregistrées dans ledit module d'observation, - pondérer, au moyen d'un module de pondération, les dits flux finaux de façon à fournir des flux finaux pondérés, - calculer, en utilisant les dits flux finaux pondérés et ladite modélisation des émissions anthropiques hebdomadaires effectuée par le module d'inventaires ascendants, les émissions hebdomadaires dudit gaz à effet de serre de ladite aire géographique, au moyen d'un module de géocodage comprenant au moins un système d'information géographique, - extrapoler, à partir des dites émissions hebdomadaires, les émissions annuelles dudit gaz à effet de serre de ladite aire géographique. Tout d'abord et avant toute chose, il convient d'établir qu'au sens de la présente invention et dans toute la suite de la présente description, il convient d'interpréter le mot « module » au sens informatique du terme. En effet, tous les modules de la méthode de mesure selon la présente invention, et notamment le module d'observation, sont, de préférence, implémentés sous forme de software, hardware ou d'une combinaison des deux. Chaque module de la méthode peut avantageusement, en fonction de son rôle, être implémenté en utilisant des moyens informatiques, notamment des moyens de calculs (ordinateurs, serveurs dédiés, mainframes, etc), des systèmes de communication (WAN, LAN, INTERNET), mais également des logiciels, notamment des systèmes de gestion de base de données, des logiciels de modélisation, de calcul etc. La méthode de mesure peut également être implémentée sous forme d'un logiciel unique, éventuellement accessible online via Internet. 12 L'initialisation de la méthode commence donc par la prise de mesures journalières de concentration du gaz à effet de serre considéré en une première pluralité de lieux répartis sur la totalité du globe terrestre et par l'enregistrement de ces mesures dans un module d'observation. Parallèlement vient une seconde étape durant laquelle on effectue des mesures journalières des flux du gaz à effet de serre considéré en une deuxième pluralité de lieux répartis sur la totalité du globe, mesures qui sont également enregistrée dans le module d'observation. Dans une troisième étape, on effectue des mesures de paramètres satellites, de paramètres météorologiques, de paramètres marins et de paramètres d'écosystèmes en une troisième pluralité de lieux répartis sur le globe terrestre et on enregistre également ces mesures de paramètres dans le module d'observation. La quatrième étape de la méthode consiste ensuite à extraire, au moyen d'un module d'extraction, permettant un transfert automatisé de données, des prévisions météorologiques d'au moins une source de données. Les mesures effectuées lors des trois premières étapes de la méthode sont effectuées au moyen d'une pluralité de satellites, d'avions, de stations de mesures atmosphériques, de stations de mesures marines, de navires et/ou de stations de mesures d'écosystèmes qui permettent d'effectuer des mesures en des lieux distincts répartis sur la totalité du globe comme cela sera décrit plus en détail dans la suite. De plus, les moyens de mesures peuvent également comprendre des capteurs, capteurs marins, capteurs d'écosystèmes, etc. Les dites première, deuxième et troisième pluralités de lieux peuvent donc se recouvrir dans une large mesure en fonction de l'équipement de mesure local. 13 L'étape suivante de la méthode de mesure selon l'invention consiste ensuite en l'utilisation d'un module d'échange modélisant l'évolution des flux du gaz considéré sur le globe en modélisant les sources et les puits naturels et anthropiques. Ensuite, un module d'inventaires ascendants est utilisé afin de modéliser les inventaires d'émissions des pays et installations à l'échelle l°xl°, avec leur saisonnalité, en Kg/m2/semaine. L'étape qui suit consiste en une modélisation du transport atmosphérique du gaz à effet de serre considéré, cette modélisation étant effectué au moyen d'un module de transport, sur la base de la modélisation de l'évolution des flux effectuée au moyen du module d'échange, sur la base de la modélisation des dites émissions hebdomadaires effectuée au moyen du module d'inventaires ascendants et sur la base des données de prévisions météorologiques. On obtient alors une distribution atmosphérique actuelle des fractions molaires de CO2 et autres GES sur le globe, à comparer avec les données enregistrées dans le module d'observation comme cela sera décrit plus en détail ci-dessous. L'étape suivante de la méthode de mesure selon l'invention consiste en l'utilisation d'un module d'inversion et d'assimilation de données, initialisé avec les résultats étalonnés de l'holocène fournis par le module d'échange, les résultats fournis par le module d'inventaires ascendants et le module de transport. Les observations et les inventaires ascendants sont intégrés dans ce module et la distribution atmosphérique du CO2 est échantillonnée en temps et en lieu où les observations de fractions molaires atmosphériques sont disponibles. Les flux naturels et anthropiques modélisés sont redimensionnés avec des facteurs scalaires afin de corriger 14 les flux avec les observations réelles pour obtenir des flux finaux en Kg/m2/semaine sur des mailles de l°xl°. Lors de l'étape suivante, un module de pondération permet ensuite de déterminer des flux finaux pondérés, en utilisant une modélisation des activités de production et du marché des émissions, et de valider les résultats fournis par le module d'inversion et d'assimilation de données sur des échelles continentales, régionales et nationales. L'étape suivante consiste en l'utilisation d'un module de géocodage comprenant un système d'information géographique, permettant de corriger les inventaires ascendants sur la base des dits flux finaux pondérés afin d'obtenir les émissions hebdomadaires sur des mailles de 0.l°x0.1°en Kg/m2/semaine. Lors de l'étape finale, à partir des émissions hebdomadaires, on peut ensuite extrapoler les émissions annuelles sur des mailles de 0.l°x0.1°en TCO2/an, TGES/an et TCO2eq/an avec une précision supérieure à 5% des émissions en TCO2/an, TGES/an et TCO2eq/an. La méthode de mesure selon l'invention échantillonne donc les inventaires avec une approche descendante dans l'ordre suivant : planète, continents, régions continentales, états/pays, régions locales, jusqu'aux installations émettrices. Les différentes étapes permettent la vérification à chaque niveau géographique du total des inventaires anthropiques et de réduire les incertitudes venant d'omissions de sources et de puits, de facteurs d'émissions ou de flux latéraux. Selon l'invention, la surface de ladite aire géographique peut être comprise entre lkm2 et 10000km2, en particulier ladite aire géographique peut inclure au moins une source anthropique donnée. 15 Selon l'invention, le dit gaz à effet de serre peut être choisi parmi le dioxyde de carbone (CO2), le méthane (CH4), le protoxyde d'azote (N2O), les oxydes d'azote (NOx), les hydrofluorocarbures (HFC), les hydrochlorofluorocarbures (HCFC), les chlorofluorocarbures (CFC), les perfluorocarbures (PFC), l'hexafluorure de soufre (SF6), l'ozone (03), la vapeur d'eau (H2O), le monoxyde de carbone (CO) et le dihydrogène (H2). Selon l'invention, les dites mesures journalières de concentration dudit gaz à effet de serre sur le globe, les dites mesures journalières des flux dudit gaz à effet de serre sur le globe, les dites mesures de paramètres satellites, de paramètres météorologiques, de paramètres marins et de paramètres d'écosystèmes peuvent être effectuées au moyen d'une pluralité de satellites, d'avions, de stations de mesures atmosphériques, de stations de mesures marines, de navires et/ou de stations de mesures d'écosystèmes permettant d'effectuer des mesures sur la totalité du globe. Selon l'invention, ledit module d'échange peut effectuer ladite modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module solaire modélisant le rayonnement solaire en utilisant les paramètres orbitaux de la géométrie terrestre avec un calcul de l'excentricité de la terre déterminé proportionnellement à l'excentricité de Mars. Selon l'invention, ledit module d'échange peut effectuer ladite modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module énergétique modélisant le rayonnement ondes-courtes, en incluant réflectivité, absorptivité et transmissivité de l'atmosphère, l'absorption par les gaz à effet de serre et les nuages, les variations d'albédo planétaire et l'influence 16 du trou de la couche d'Ozone, ledit module énergétique modélisant également le rayonnement ondes-longues, en utilisant l'équation de Schwartzschild, la méthode des émissivités et en incluant l'absorption et l'émission par les gaz à effet de serre et les nuages de rayonnements ondes-longues, les flux de chaleur latente, les flux de chaleur sensible, les flux de conduction et la température de surface. Selon l'invention, ledit module d'échange peut effectuer ladite modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module océan modélisant l'effet net des échanges atmosphère-océan sur la base du modèle MOM3 combiné avec lesdites données de prévisions météorologiques et prenant en compte l'effet de tampon, l'absorption par altération chimique suivant le modèle CDIAC DBl0l2 et le dégagement par évaporation. Selon l'invention, ledit module d'échange peut effectuer la modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module biosphère modélisant l'effet net des échanges atmosphère-biosphère sur la base du modèle JSBACH et incluant les types de plantes de la biosphère, l'indice foliaire, la lumière, l'albédo, la photosynthèse C3 et C4, l'addition du taux de photosynthèse brute limité, la respiration autotrophique, la respiration hétérotrophique et/ou la modification anthropique de la couverture de végétation depuis au moins le dernier millénaire. Selon l'invention, ledit module biosphère peut utiliser un module incendie modélisant les perturbations dues aux incendies sur la base de données extraites de la base de données Global Fire Emission Database (GFEDv2) intégrées au modèle JSBACH. 17 Selon l'invention, ledit module d'échange peut effectuer ladite modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module fossile modélisant les émissions anthropiques fossiles à l'échelle globale sur la base des statistiques de production de pétrole et de charbon de l'Energy Information Administration (EIA) et les estimations de Etemad et al. Selon l'invention, ledit module d'inventaires ascendants peut extraire des inventaires d'émissions de la base de données EDGAR 4.0 en incluant un calcul de la variabilité temporelle des émissions. Selon l'invention, ledit module de transport atmosphérique peut utiliser le modèle de transport TM5 combiné avec les dites données de prévisions météorologiques pour calculer le transport atmosphérique des flux dudit gaz à effet de serre sur le globe. Selon l'invention, ledit module d'inversion et d'assimilation de données peut utiliser, pour calculer les dits flux finaux, une inversion de synthèse avec la fonction de Green pour les grandes régions et le filtre d'ensemble de Kalman. Selon l'invention, ledit module de pondération peut utiliser, pour pondérer les dits flux finaux, une analyse des activités de production des pays et des régions du monde accompagnée d'une modélisation des marchés des émissions basée sur le modèle de biens publics produits par le privé. Selon l'invention, ledit module de géocodage peut utiliser des coefficients correcteurs. Le système de mesure selon l'invention, effectuant la méthode de mesure telle que décrite ci-dessus, comprend : 18 - des moyens de mesure des concentrations et des flux des gaz à effet de serre, - des moyens de mesures de paramètres satellites, météorologiques, marins et d'écosystèmes, - au moins une base de données centralisée comprenant un module d'observation, - des moyens d'extraction et de transfert de données automatisés, - des moyens de calcul comprenant au moins un module d'échange, au moins un module d'inventaires ascendants, au moins un module de transport, au moins un module d'inversion et d'assimilation de données, au moins un module de pondération, - au moins un module de géocodage comprenant un système d'information géographique permettant de géocoder les résultats fournis par les dits moyens de calcul, - une plateforme internet centralisée permettant de visualiser et d'analyser les émissions de gaz à effet de serre d'une pluralité d'aires géographiques données. Selon l'invention, le système de mesure peut comprendre des moyens matériels et logiciels d'interface avec un système de gestion de production d'une installation. L'invention sera mieux comprise par l'homme du métier grâce à 25 la description détaillée de modes d'exécution en relation avec les dessins accompagnants, dans lesquels: la figure 1 représente un diagramme en bloc illustrant le processus et les composantes de la méthode, la figure 2 représente un diagramme en bloc illustrant 30 l'échantillonnage des flux anthropiques de CO2, la figure 3 représente un tableau illustrant l'espérance de vie et les potentiels de réchauffement globaux des GES, 19 la figure 4 représente un diagramme conceptuel d'échantillonnage du satellite GOSAT, la figure 5 représente les caractéristiques des observations satellites, la figure 6 représente les sites d'observations atmosphériques, la figure 7 représente le réseau océanique de mesure de pCO2 de surface, la figure 8 représente les sites d'observations d'écosystèmes, la figure 9 représente un diagramme en bloc illustrant le module d'échange, le module d'inventaires ascendants et le module de transport, la figure 10 représente un schéma illustrant l'orbite 15 terrestre autour du soleil, la figure 11 représente un schéma illustrant la radiation solaire sur le haut de l'atmosphère, la figure 12 représente une figure illustrant le module énergétique, 20 la figure 13 représente les inventaires EDGAR 4.0 sur des mailles 0.1°x0.1° en TCO2eq, la figure 14 représente un diagramme en bloc illustrant le module d'observations et le module d'inversion et d'assimilation de données, 25 la figure 15 représente un diagramme présentant trois cycles d'assimilation de données, la figure 16 représente un diagramme en bloc illustrant le module de pondération, la figure 17 représente un diagramme en bloc illustrant le 30 module de géocodage, la figure 18 représente un diagramme en bloc illustrant un système de mesure selon l'invention. Description détaillée de l'invention 20 La figure 1 représente, de manière générale, les différentes étapes de la méthode de mesure selon l'invention. L'invention est liée à une méthode de mesure et un système de mesure précis d'inventaires de GES incluant le CO2, CH4, N2O, NOx, HFC, HCFC, CFC, PFC, SF6, 03, H20, CO et H2 à partir de leurs sources et de leurs puits naturels et anthropiques sur une aire géographique déterminée, en particulier sur une aire dont la surface est comprise entre lkm2 et 10000km2. La méthode est initialement présentée sur le CO2 et le même processus est utilisé pour les autres GES. 1. CO2 Pour chacun des gaz à effet de serre considérés, et plus particulièrement pour le CO2, la méthode de mesure selon l'invention comprend des mesures in-situ de CO2 effectuées à partir d'une combinaison d'observations (Fig.l Bloc 100) associant des mesures satellites (Fig.l4 Bloc 101), des mesures aériennes (Bloc 102), des mesures atmosphériques (Bloc 103), des mesures marines (Bloc 104) et des mesures d'écosystèmes (Bloc 105). Elle comprend également la modélisation des flux de CO2 sur le globe, initialisé à partir de l'holocène, en utilisant un module d'échange (Fig. 1 Bloc 200). Ce module d'échange utilise un module solaire (Fig. 9 Bloc 201), un module énergétique (Bloc 202), un module océan (Bloc 203), un module biosphère (Bloc 204), un module incendie (Bloc 205) et un module fossile (Bloc 206). Un module d'inventaires ascendants permet d'obtenir une distribution spatiale précise des inventaires maillés d'émissions anthropiques (Bloc 301). Les fractions molaires sont modélisées sur l'atmosphère entière avec un module de transport (Bloc 400). La différence entre les observations et les prévisions du modèle est minimisée par le module d'inversion et d'assimilation de données (Bloc 21 0) en ajoutant une inversion de synthèse avec la fonction de Green (Bloc 501), suivie par le filtre d'ensemble de Kalman (Bloc 502) pour obtenir des flux finaux comme cela sera décrit plus en détail ci-dessous. La méthode de mesure selon l'invention ajoute une validation et une modulation des flux finaux par un module de pondération (Bloc 600) pour obtenir des flux finaux pondérés, puis utilise enfin un système de coefficients correcteurs permettant d'obtenir, sur la base des flux finaux pondérés (bloc 700), des émissions calculées à l'échelle des installations en TCO2/an, TGES/an et TCO2eq/an. Un aperçu des résultats des différents modules de la méthode est le suivant : Blocs Modules Résultat Maillage 100 Module d'observation Concentrations, flux et paramètres Local .......................................................................... .......................................................................... .......................................................................... ............. . 201 Module Solaire W/m2/jour 202 Module Energétique W/m2/jour 203 Module Océan PgC/mois 5°x4° .......................................................................... .......................................................................... .......................................................................... ............. . 204 Module Biosphère Kg/m2/s l°xl° 205 Module Incendie Kg/m2/mois l°xl° .......................................................................... .......................................................................... .......................................................................... ............. . 206 Module Fossile T/an Global .......................................................................... .......................................................................... .......................................................................... ............. . 300 Module d'inventaires ascendants Kg/m2/semaine l°xl° .......................................................................... .......................................................................... .......................................................................... ............. . 401 Module de Transport ppm/s l°xl° .......................................................................... .......................................................................... .......................................................................... ............. . 500 Module d'inversion et d'assimilation Kg/m2/semaine l°xl° .......................................................................... .......................................................................... .......................................................................... ............. . 600 Régions Module de pondération T/semaine GES en Kg/m2/semaine, TCO2/an, 700 et pays Module de géocodage TGES/an et TCO2eq/an 0.1°x0.1° La méthode est tout d'abord proposée sur le CO2, puis présente les sources et les puits des autres GES pour lesquels le même processus est appliqué. Elle inclue notamment pour le CH4, le modèle utilisé pour modéliser les émissions du Pergélisol et du fond des océans provenant des 22 hydrates de CH4 comme cela sera décrit plus en détail au point II. 1. Module d'observation Des observations in-situ en continu, de haute précision et à long terme sont nécessaires afin de comprendre les processus d'échanges du cycle du carbone et de réduire les incertitudes des estimations. La méthode améliore le Carbontracker en effectuant une combinaison de mesures satellites, aériennes, atmosphériques, d'écosystèmes et marines (Bloc 100, Fig. 14) afin d'obtenir une observation globale de la planète avec une couverture des différentes couches de l'atmosphère et des surfaces de la planète. Ceci permet d'obtenir une cartographie complète quasi temps-réel des sources et des puits de GES des échelles mondiales, continentales, étatiques, nationales, locales, jusqu'au niveau des installations pour refléter la réalité des niveaux d'émissions (Fig. 2). De plus, la réunion de ces cinq techniques complémentaires et essentielles d'observation permet l'obtention d'un traitement de concentrations et de flux homogènes sur l'ensemble du globe avec une visibilité locale détaillé pour obtenir des résultats d'inversion et d'assimilation précis à l'échelle 1°xl°. Toutes les mesures acquises sont enregistrées dans le module d'observation. • Observations satellites Les observations satellites (Bloc 101, Fig. 14) incluent préférablement la combinaison du satellite Japonais GOSAT (Maksyutov et al. 2008) et du satellite Européen ENVISAT (Bovensmann et al. 1999) pour obtenir une couverture de mesure globale des différentes couches de l'atmosphère, de l'échelle planétaire jusqu'au niveau des installations. Ces satellites mesurent le rayonnement solaire proche infrarouge réfléchi par la surface de la terre et l'atmosphère, ce qui 23 permet de détecter l'absorption atmosphérique des GES dans ces régions spectrales. Ceci requiert une forte sensibilité des mesures jusqu'à la surface où les signaux de sources et de puits sont les plus forts et le satellite GOSAT avec ses deux instruments TANSO- FTS et TANSO-CAI actuellement en orbite (Fig. 4) offre cette information avec une précision relative de mesure de l'ordre de 1% (4 ppm) et une empreinte de 10km de diamètre (Fig. 5). Les grandes sources fossiles locales telles que les installations émettrices augmentent les concentrations de CO2 dans l'atmosphère de 1 à 10 ppm sur la source et sont généralement dispersés sur quelques dizaines de kilomètres autour. Le satellite GOSAT démontre ses capacités de surveillance par installation et les échantillonnages fournissent des données indépendantes à comparer avec celles du module d'inventaires ascendants (Bloc 300). Avec de 100 à 500 grandes sources locales dans les pays de forte émission, il est possible d'obtenir un échantillon statistique de mesures des panaches de CO2 émis par ces grandes sources dans ces pays. La méthode utilise la combinaison des deux algorithmes Full Physics (FP) et Apparent Optical Path Difference (AOPD) (Boland et al. 2009) (Bloc 106) afin de récupérer les colonnes de CO2 (XCO2) des rayonnements mesurés par GOSAT. Avec une exactitude de 0.3% à 0.5% (1 à 2 ppm), une zone d'échantillonnage inférieure à 3km2, le satellite américain OCO (Crisp et al. 2004) a également été conçu pour mesurer l'accroissement de concentration au dessus des sources et des puits locaux. La méthode complètera les mesures de CO2 avec le successeur du satellite OCO lorsqu'il sera disponible, OCO ayant manqué son lancement début 2009. La méthode utilise également le spectromètre SCIAMACHY d'ENVISAT actuellement en orbite et la combinaison des deux algorithmes WFM-DOAS et BESD suivant Buchwitz et al. (2008) 24 (Bloc 106) afin de récupérer les colonnes de CO2 (XCO2) à partir des rayonnements mesurés. Le développement de ces algorithmes progresse et atteint actuellement 2-3% de précision selon Schneising et al. (2008) avec une résolution horizontale de 30kmx60km (Fig. 5). Les recherches ont l'objectif d'atteindre 1% de précision relative, ce qui est suffisant car un offset constant est pris en compte dans l'inversion et l'assimilation de données et un haut niveau de précision relative est requis pour valider les modèles. Les produits de données de GOSAT et de SCIAMACHY sont les colonnes moyennes de fractions molaires de CO2 en air sec (XCO2, ppm). Pour la mesure des autres GES et la validation des données de CO2, les satellites GOSAT et ENVISAT (SCIAMACHY) sont complétés par AIRS, IASI, TES et OMI (Fig. 5). La méthode de mesure selon l'invention utilise également l'Atmospheric Infrared Sounder (AIRS) (Aumann et al., 2003) qui est un sondeur infrarouge multi spectral de haute résolution sur le satellite AQUA conçu pour apporter des données précises sur l'atmosphère, la surface et les océans et offre des mesures sur la température atmosphérique, les profils d'humidité, de température de surface et de GES tels que 03, CO, CO2, CH4 et H2O. La méthode de mesure selon l'invention acquiert également des mesures grâce à l'Infrared Atmospheric Sounding Interferometer (IASI) (Crevoisier et al. 2009) qui est un spectromètre à transformée de fourrier sur le satellite METOP et apporte des mesures infrarouges des profils de température dans la troposphère et la stratosphère basse, les profils d'humidité dans la troposphère et de GES tels que CO2, CH4, N2O, CO, H2O et 03. 25 La méthode de mesure selon l'invention effectue par ailleurs des mesures grâce au Tropospheric Emission Spectrometer (TES) (Luo et al., 2007) qui est un spectromètre à transformée de Fourier à bord d'EOS AURA offrant une discrimination des espèces moléculaires radiativement active dans le bas de l'atmosphère. TES emploie à la fois les émissions naturelles thermiques de la surface et de l'atmosphère et la lumière du soleil réfléchie apportant une couverture jour-nuit sur le globe avec des mesures de CO2, CO, CH4, 03, H2O et NO2. En outre, la méthode de mesure selon l'invention effectue des mesures au moyen de l'Ozone Monitoring Instrument (OMI) (Levelt et al. 2000) qui est un spectromètre à bord d'EOS AURA mesurant le spectre de longueurs d'ondes ultraviolet/visible/proche infrarouge avec une haute résolution spectrale. OMI fournit notamment les colonnes totales de mesures troposphériques et stratosphériques d'O3, H2O et NO2 ainsi que les profils stratosphériques d'O3, l'albédo de la surface, les aérosols et les paramètres de nébulosité. La méthode utilise également de préférence les données de l'instrument MODIS des satellites Terra et Aqua qui apporte des données objectives du changement de couverture des terres (ALCC, Anthropogenic Land Cover Change). Les données satellitaires de GES sont préférablement validées par les réseaux de spectromètres à Transformée de Fourrier au sol, Network for the Detection of Atmospheric Composition Change (NDACC) (Kurylo, 1991) et Total Carbon Column Observing Network (TCCON) (Toon, 2009). Ces stations enregistrent les spectres solaires directs dans la région spectrale proche-infrarouge avec pour NDACC la mesure de 03, CO, CO2, N2O, CH4 et pour TCCON, celle de CO2, CH4, N2O, CO et H2O. 26 Les paramètres mesurés et leur fréquence lorsqu'ils sont disponibles incluent: Mesures Fréquence GES (dont CO2, CH4,CO3N20, NOx, H2O, 03) Continu .......................................................................... .......................................................................... .... . ALCC Continu .......................................................................... .......................................................................... ..... . Albédo Continu .......................................................................... .......................................................................... ................................................................... . • Observations aériennes La méthode de mesure selon l'invention complète les mesures satellites par des mesures effectuées grâce à des observations aériennes, mesures qui sont également enregistrées dans le module d'observation. Les observations satellites sont complétées par les observations aériennes disponibles effectuées avec le NOAA ESRL Carbon Cycle Greenhouse Gases group (CCGG) Air sampling, ainsi que par les mesures du programme In-service Aircraft for a Global Observing System - European Research Infrastructure (IAGOSERI). Le NOAA ESRL Air sampling permet d'effectuer des mesures de profils verticaux de CO2, CH4, N20, CO, H2 et SF6. IAGOS-ERI a pour origine le programme MOZAIC (Marenco et al. 1998) et inclue le programme CARIBIC (Schuck et al. 2009) et apporte des observations GES in-situ de haute-qualité dans la tropopause incluant CO2, CH4, CO, N20, H20, 03, CFC, HFC et HCFC. Mesures Fréquence GES (dont CO2, CH4, CO, N20, H20, 03, CFC, Continu et échantillonnage HFC, HCFC, SF6, H2) hebdomadaire • Observation atmosphériques La méthode de mesure selon l'invention effectue aussi des mesures de concentration atmosphérique et d'échantillons prélevés sur les sites NOAA ESRL Cooperative Global Air 27 Sampling Network et le CSIRO Air Sampling Network pour chaque année. Elle utilise également des séries chronologiques in situ quasi continues des tours et observatoires NOAA ESRL. Ces observations sont étalonnées sur le standard mondial (WMO-2005). La méthode complète ces observations atmosphériques (Bloc 103, Fig. 14) en incluant préférablement celles actuelles et en développement du Global Atmosphere Watch (GAW) incluant les stations WDCGG, de 1' International Global Atmospheric Chemistry Observations (IGACO), du GCOS Reference Upper-Air Network (GRUAN), du Network for the Detection of Atmospheric Composition Change (NDACC) dont les stations LIDAR, de l'Integrated Carbon Observation System (ICOS), du System for Observation of Halogenated Greenhouse Gases in Europe (SOGE) et du réseau ALE/GAGE/AGAGE (Fig. 6). Chaque station est un observatoire qui mesure continuellement la variabilité régionale et mondiale des concentrations de CO2 (ppm), de GES ainsi que les paramètres météorologiques. Elles sont utilisées pour détecter les changements long termes dans les tendances de concentrations et la variabilité interannuelle associée à des émissions anthropiques et anomalies climatiques. Certaines stations sont également équipées des instruments de mesures de flux. Chaque station est représentative d'une superficie généralement de plus de 100 km. Les mesures de concentration de CO2 sont effectuées idéalement avec une précision inférieure à 1 ppm et des échantillons d'air sont également collectés, de préférence sur une base hebdomadaire puis analysés. Les paramètres mesurés et leur fréquence lorsqu'ils sont disponibles incluent: Mesures Fréquence CO2 Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ...... . Paramètres météorologiques (pression, température, humidité relative, vent) Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ....... Hauteur de la couche limite Continu (30 min) 28 Flux CO2 Continu (30 min) .......................................................................... .......................................................................... .............................................. . GES (dont CH4, CO, N20, NOx, H20, 03, CFC, Continu et échantillonnage HFC, HCFC, SF6, H2, PFC) hebdomadaire • Observations marines La méthode de mesure selon l'invention effectue également des mesures via des observations marines (Bloc 104, Fig. 14) effectuées au moyen d'un réseau de navires instrumentés naviguant les océans et de stations fixes (Fig. 7). Les navires sont en général des navires commerciaux, ferrys, porte-conteneurs et bateaux citernes opérant sur des itinéraires répétés. Les stations fixes sont des sites sur l'océan sur lesquels des observations temporelles continues sont enregistrées par le biais d'amarres et de navires de recherche. La couverture doit être suffisante pour quadriller les flux océaniques air-mer des régions océaniques (Pacifique, Atlantique, Indien, Sud, Arctique). La méthode inclue les observations des programmes International Ocean Carbon Coordination Project (IOCCP), IOCCP underway liner, JCOMM VOS, IOCCP time series (Oceansites), IOCCP Hydrography (GO-SHIP), CarbonOcean-IP, SOLAS-IMBER Carbon group (SIC), Carbon Dioxide Information Analysis Center Ocean CO2 Center (CDIAC), National Oceanic & Atmospheric Administration (NOAA) VOS, Climate Variability and Predictability Research (CLIVAR) et Integrated Carbon Observing System (ICOS). Les navires et stations fixes sont équipés d'instruments automatisés mesurant la concentration atmosphérique et la pression partielle de CO2 de la surface, la température de surface et la salinité. Certains navires et stations marines sont équipés d'instruments pour la mesure de concentration atmosphérique de GES additionnels répétées à intervalles journalières et mensuelles et des échantillons d'air sont régulièrement recueillis puis analysées. Les flux air-mer 29 sont calculés à partir des mesures de pression partielle de CO2, tel qu'effectué dans le Carbontracker en utilisant le principe d'inversion de Jacobson et al. (2007) (Bloc 107). Les paramètres mesurés et leur fréquence lorsqu'ils sont disponibles incluent: Mesures Fréquence CO2 atmosphérique continu (30 min) .......................................................................... .......................................................................... ......................................................................... . pCO2 océan, pression totale atmosphérique continu (30 min) .......................................................................... .......................................................................... ......................................................................... . Salinité et température de surface des océans continu (30 min) .......................................................................... .......................................................................... ......................................................................... . Paramètres Météorologiques 4-heures .......................................................................... .......................................................................... ......................................................................... . GES (dont CO2, CH4, CO, N20, H2O, HFC, SF6, Continu et échantillonnage H2) mensuel • Observations d'écosystèmes La méthode effectue également des mesures via des observations d'écosystèmes (Bloc 105, Fig. 14) du programme Fluxnet/iLEAPS (Baldocchi et al. 2001) qui est un réseau de réseaux régionaux et inclue préférablement les stations d'écosystèmes actuelles et en développement de Carboeurope-IP, CarboAfrica, Asiaflux, Afriflux, Ozflux, Large-Scale Biosphere-Atmosphere (LBA), US-China Carbon Consortium (USCCC), Nordic Centre for Studies of Ecosystem Carbon Exchange and its Interactions with the Climate System (NECC), TCOS-Siberia, ChinaFlux, Ameriflux, Fluxnet-Canada, KoFlux ainsi que l'Integrated Carbon Observation System (ICOS) (Fig. 8). Chaque station mesure en continu les flux de CO2, les flux d'eau et d'énergie entre les écosystèmes terrestres et l'atmosphère ainsi que des variables d'écosystème telles que les variables météorologiques, budgets hydrologiques et de rayonnement et les réservoirs de carbone dans la végétation et les sols. Les stations transfèrent les données collectées des flux d'écosystèmes, de préférence quotidiennement.30 Certaines stations sont également équipées des instruments de mesures des concentrations des stations atmosphériques. Les données sont utilisées pour paramétrer et valider les modèles de carbone appliqués aux échelles continentales, détecter les changements à long terme dans les puits et sources et identifier l'impact des différences de gestion du budget de carbone. L'empreinte de chaque tour est en moyenne entre 200 et 1000 mètres. Les flux (Kg/m2/s) sont mesurés en utilisant la méthode de covariance des turbulences (Eddy covariance) à partir de mesures directes de la vitesse du vent vertical et des concentrations de CO2 pour déterminer les flux verticaux turbulents au sein des couches atmosphériques limites (Bloc 108). Des échantillons d'air pour la mesure de GES sont régulièrement recueillis puis analysés. Les paramètres mesurés et leur fréquence lorsqu'ils sont disponibles incluent: Mesures Fréquence Flux de chaleur sensible, CO2, H2O Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ...... . Profil vertical de CO2 Continu (30 min) Radiation globale nette réfléchie et diffuse Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ....... Profil de température de l'air et du sol Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ....... Profil de vitesse de vent Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ...... . Profil de contenu d'eau du sol Continu (30 min) chute et hauteur de neige au sol Continu (30 min) Flux de chaleur du sol Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ....... Contenu de Carbone du sol Echantillonnage sur 5 ans Biomasse Annuel . .......................................................................... .......................................................................... .......................................................................... ...... Gestion et perturbations Annuel . .......................................................................... .......................................................................... .......................................................................... ...... Flux de CH4 Continu (30 min) / Journalier .......................................................................... .......................................................................... .......................................................................... ....... . Flux de N20 Continu (30 min) / Journalier .......................................................................... .......................................................................... .......................................................................... ....... . Température de la canopée Continu (30 min) Réflectance Spectrale Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ...... . Radiation Photosynthétique Active sous la Continu (30 min) canopée Niveau d'eau du sol Continu (30 min) .......................................................................... .......................................................................... .......................................................................... ....... .......................................................................... .......................................................................... .......................................................................... ...... . Précipitations, 31 Flux de Sève Continu (30 min / 3 heures) .......................................................................... .......................................................................... .......................................................................... ...... . Respiration du Sol Continu (3 heures) .......................................................................... .......................................................................... .......................................................................... ...... . Caméra de phénologie Journalier .......................................................................... .......................................................................... .......................................................................... ....... Déposition de N Bihebdomadaire .......................................................................... .......................................................................... .......................................................................... ....... . Contenu des feuilles et du sol en N Bihebdomadaire Déposition de litière Mensuel .......................................................................... .......................................................................... .......................................................................... ...... . Import et export de C et N due à la gestion de la biosphère Annuel GES (dont CO2, CH4, CO, 03, N20, NOx, H2O, SF6, Continu et échantillonnage H2, HFC, HCFC, PFC) hebdomadaire 2. Module d'échange Parallèlement à la prise de mesures, la méthode selon l'invention effectue une modélisation de l'évolution des flux de GES, notamment CO2, à partir de l'holocène, en utilisant un module d'échange (Fig.9). Le module d'échange comprend un module solaire, un module énergétique, un module océan, un module biosphère, un module incendie et un module fossile. a. Module Solaire La méthode de mesure selon l'invention améliore le calcul du rayonnement solaire du Carbontracker en utilisant un module solaire (Bloc 201, Fig. 9) qui modélise le rayonnement solaire avec une influence plus précise sur les échanges entre l'atmosphère, les océans et la biosphère. L'insolation solaire est le montant du rayonnement solaire atteignant la terre par latitude et par saison et se réfère au rayonnement arrivant sur le haut de l'atmosphère (TOA, Top Of Atmosphere). D'après la théorie orbitale des paléoclimats, les variations de l'orbite terrestre dans le temps ont contribuées à changer le montant du rayonnement solaire reçu par la Terre pour chaque saison et ont dirigé les alternances des périodes glaciaires et interglaciaires. Selon les cycles de Milankovitch, trois paramètres de la géométrie terrestre 32 orbitale sont utilisés pour évaluer le forçage orbital: l'obliquité, qui est l'inclinaison de l'écliptique comparativement à l'équateur céleste avec un cycle d'environ 40 milliers d'année (Ka), l'excentricité de l'orbite terrestre autour du soleil avec un cycle d'environ 100Ka et la précession climatique, liée à la distance Terre/Soleil au solstice d'été avec un cycle d'environ 26 Ka. Suivant cette théorie, les périodes interglaciaires ont tendance à se produire pendant les périodes de plus intense rayonnement solaire estival dans l'hémisphère Nord et depuis environ 11'700 ans, la Terre est entrée dans un nouveau cycle interglaciaire appelé l'Holocène. Comme nouvelle approche, la méthode de mesure selon l'invention débute donc la modélisation des échanges de CO2 à t=0 à partir du début de l'holocène afin d'obtenir une base stable et étalonnée des échanges naturels pour déterminer avec plus de précision l'influence future des émissions anthropiques de l'ère industrielle. Pour le calcul de l'insolation solaire, les paramètres orbitaux de la géométrie terrestre sont obtenus en utilisant la théorie des Variations Séculaires des Orbites Planétaires de Bretagnon (1987). L'inventeur de la méthode ajoute également une modification du calcul de l'excentricité, car à sa connaissance, aucune influence précise de la lune sur la dissipation de marées solides et océaniques de la terre, n'a été prise en compte pour calculer la perturbation sur l'excentricité de la terre et l'insolation solaire incidente. Phobos, un des satellites de Mars est utilisé pour évaluer cette influence car c'est le meilleur cas connu d'évolution orbitale rapide d'un satellite dans le système solaire avec une période orbitale de seulement 7,65 heures, comparativement à 27,3 jours pour la lune. Son mouvement orbital a été intensivement étudié depuis 33 sa découverte en 1877 ou il a exécuté depuis, environ 145'500 orbites, équivalent à une période de 10'880 ans pour la lune. L'instrument Mars Orbiter Laser Altimeter (MOLA) sur le satellite Mars Global Surveyor a observé les transits de l'ombre de Phobos sur la surface de Mars, et a directement mesuré la distance avec Phobos afin de vérifier si les positions observées de Phobos et de son ombre sont en bon accord avec les modèles. Compte-tenu de la longue période et de la précision des observations, l'accélération séculaire de Phobos est utilisée pour déterminer le facteur de qualité (Q) de Mars, qui exprime le taux relatif de dissipation d'énergie et qui est associé au nombre de Love (k2), décrivant les propriétés élastiques des planètes. La mesure précise de ces paramètres permet de déterminer l'effet de la dissipation d'énergie sur l'excentricité de Mars. Avec sa forte proximité avec Mars, l'orbite de Phobos connaît une accélération séculaire orbitale qui est utilisée pour évaluer celle de la lune avec la terre. Sur la terre, la Lune exerce une attraction gravitationnelle causant des marées solides et océaniques. La terre induit une accélération séculaire, qui a un effet cumulatif sur la position de la lune lorsque extrapolée au fil des siècles. L'effet de cette accélération séculaire de la Lune est assez mal connu car les enregistrements de ses déviations remontent à environ un siècle. Une mesure précise de l'excentricité de Mars permet ainsi de déduire la perturbation de la lune sur l'excentricité de la Terre. La méthode évalue ce facteur comme une proportion de la distance Phobos-Mars, de la distance Lune-Terre et des excentricités respectives: eTerre eMars DLuneùTerre DPhobosùMars L'insolation instantanée est définie comme l'énergie reçue par unité de temps et de surface sur un plan horizontal à TOA 34 et la méthode suit l'approche de Liou (2002) pour son calcul. La trajectoire de la terre autour du Soleil est une ellipse (Fig. 10). Le point le plus proche de l'orbite terrestre au soleil est appelée le Perihelion, tandis que l'Aphelion est le plus éloigné. La forme en ellipse est caractérisée par son excentricité e = J(a2ùb2)/a. La distance (r) de la Terre au Soleil est calculée comme une fonction de y, la véritable anomalie de l'ellipse selon la première loi de Kepler. a(1ùe2) r= 1+ecosv. Le montant de la radiation solaire incidente par unité de 10 surface à TOA est une fonction de r et la distance moyenne soleil-terre (ro) est définie selon la seconde loi de Kepler: r 2 = a2J(1-e2) a2 En moyenne, le montant d'énergie solaire incidente en dehors de l'atmosphère terrestre est la multiplication de la constante solaire So par la surface qui intercepte les rayons 15 solaires. Sr est le montant du rayonnement solaire par unité de surface mesuré sur la surface extérieure de l'atmosphère dans un plan perpendiculaire aux rayons à une distance (r) du soleil et est fonction de Sa. A TOA, une surface à la distance moyenne Terre-Soleil perpendiculaire aux rayons 20 reçoit Sr = So rZ/r2 . Le montant d'énergie solaire reçu par unité de temps sur une unité de surface horizontale à TOA est fonction de 6o, l'angle solaire zénithal (Fig. 11). Sh est déduit en fonction de l'orientation des rayons solaires et de la normale à la surface terrestre selon : r2 Sh = Sr cos B0 = So 2 cos 00 r2 25 L'axe de rotation de la Terre n'est pas perpendiculaire à son plan orbital et est incliné relativement au plan céleste équatorial avec un angle s. L'équinoxe de printemps est utilisé comme une référence pour définir la véritable 35 longitude 2 avec w, la longitude du Perihelion mesurée de l'équinoxe d'automne (Fig. 11). De la trigonométrie sphérique, l'angle solaire zénithal dépend de la latitude 0 du point de la Terre, de la déclinaison solaire ô et de l'angle horaire h, selon: cos Bo = sin sin' + cos cos' cos h où h indique le temps depuis lequel, le soleil était au méridien local, mesuré à partir de l'observateur du méridien vers l'ouest. ô est défini comme l'angle entre une ligne du centre de la terre vers le soleil et l'équateur céleste. H représente une demi-journée et est défini par cos H = - tan ç tan S . Connaissant la véritable longitude et l'obliquité, ô varie au gré des saisons selon sin' = sins.sinil et l'énergie solaire reçue par unité de surface par jour est calculée selon: r 2 Sh,day ''s'' ( r ) (HsinOsin 8 +cos cos' sin H ) Au cours de l'holocène, cette énergie solaire est principalement influencée par la précession, puis par l'obliquité. De ce calcul, la précession était à son point culminant au début de l'Holocène contribuant majoritairement à ce stade à la plus haute insolation et a diminué jusqu'à son minimum vers environ 1300AD. Depuis lors, la précession augmente jusqu'à un maximum vers environ lOKa AD finissant son cycle. Cette insolation augmente donc depuis 1300AD et exerce une influence importante sur les flux de CO2 des océans et de la biosphère. Cette modélisation de la radiation solaire en W/m2/jour sur des mailles l°xl° est utilisée dans le module énergétique (Bloc 202). Le calcul de la radiation solaire TOA est initialisé au début de l'holocène avec une périodicité de 50 ans et un calcul global sur la planète. 36 2953021 b. Module énergétique La méthode de mesure selon l'invention améliore également le Carbontracker en utilisant un module énergétique (Bloc 202, Fig. 9) qui modélise les rayonnements ondes-courtes et ondes- 5 longues avec une influence plus précise sur les échanges entre l'atmosphère, les océans et la biosphère. Il inclut notamment un calcul plus précis de l'absorption du rayonnement solaire par les GES, de l'influence du trou de la couche d'Ozone et de l'effet de serre. 10 Le calcul du rayonnement ondes-courtes incident est effectué selon la méthode d'Huybers et al. (2007). Ce rayonnement, principalement d'origine solaire, est le résultat de multiples processus de dispersion et d'absorption impliquant essentiellement les molécules H20, 03, les aérosols, les 15 nuages, l'air et la surface sous-jacente (Fig. 12). Ceci est modelé en utilisant l'insolation solaire du module solaire (Bloc 201), la réflectivité (R), l'absorptivité (A), la transmissivité (T) de l'atmosphère et des nuages et l'albédo (a). À l'équilibre énergétique, A+R+T=l avec A, R et T, des 20 paramètres évoluant en fonction du climat. Le dénominateur prend en compte l'absorption et la réflexion du rayonnement par la surface de multiples fois selon: T(1 ùa) Ss 1 Sh,day a. R L'albédo est déterminé à partir du Carbontracker pour les océans et l'atmosphère, à partir du nouveau module biosphère 25 pour la surface terrestre et est validé par les observations satellites. Son calcul est important car il varie principalement en fonction de la nébulosité, de la neige, de la glace, de la surface foliaire et des changements de couverture terrestre. 37 La méthode de mesure selon l'invention améliore également le calcul du rayonnement ondes-courtes en ajoutant l'influence du trou de la couche d'ozone car l'ozone est un excellent absorbeur de rayons UV. En 1970, on a découvert que l'ozone est détruit par des radicaux dont l'hydrogène, l'azote, le chlore et le brome. Avec l'épuisement de la couche d'ozone, le filtre de protection prévu par l'atmosphère est progressivement réduit. Pendant la période du solstice d'hiver, le trou de la couche d'ozone est principalement situé au-dessus de l'Antarctique et l'augmentation de l'insolation sur ce pôle a fortement contribué à l'augmentation du rayonnement ondes-courtes aggravant le réchauffement. L'ozone est un indicateur clé intégré dans la méthode pour mesurer cette augmentation et est calculé avec sa transmissivité TO3 selon Tripathi et al. (2000) : T03 = e Où a est le coefficient d'absorption de l'ozone, µ est le ratio des longueurs des chemins actuels et verticaux à travers la couche d'ozone et 52 est la concentration en ozone. Pour atteindre l'équilibre climatique, l'énergie solaire incidente absorbée par le système terre/atmosphère est contrebalancé par un montant égal d'énergie IR thermique émis (Fig. 12). La terre a une atmosphère qui absorbe et émet un rayonnement ondes-longues et cet effet de serre sert à maintenir la chaleur près de la surface. L'absorption est dépendante de la longueur d'onde et est déterminé par la composition atmosphérique, les nuages, les aérosols et les concentrations de GES. Le traitement du rayonnement ondes- longues s'appuie sur l'équation de Schwartzschild selon Washington et al. (2005) et prend en compte l'absorption et 38 l'émission avec les lois de Lambert et de Kirchhoff avec le changement d'intensité de radiation exprimé selon: dl = -Ikpdz + B(T)kpdz où k est le coefficient d'absorption, p est la densité du médium et B(T) est la fonction de Planck. L'intégration sur tous les angles de l'hémisphère sur une surface horizontale transforme les intensités en flux ascendants et descendants. La simplification du modèle est effectuée en utilisant la méthode d'émissivité dans laquelle l'intégration sur de relativement larges bandes résulte dans le calcul des flux ascendants Fî et descendants F;1 avec les émissivités e' et e qui sont fonctions de la vapeur d'eau, de la pression et de la température sur le chemin que traverse le rayonnement. L'influence des GES tels que l'O3, CH4, NO2 et CO2 est également modélisée dans l'absorption (A). FS 1 (z) = nB(0) + fo s' (z, z')d(riB(z )) ; FS 1 (z) = nB(z0)E(z,w) + fZO «z, z »d(nB(z )) = , dB«z ,) = Bi z (z,z ) a,(z,z )dBz,) ~(z,z ) A,(z,z ) B(z,) i i Le sol, la surface et l'atmosphère échangent de la chaleur par contact direct entre la surface et l'air (chaleur sensible H9T), par l'évaporation et la transpiration (chaleur latente LVET) calculés selon la méthode utilisée par Xing et al. (2007) et par absorption dans le sol (conduction G9 ) selon la loi générale de Fick. Sans transfert de chaleur latente et sensible, la surface de la terre aurait une température beaucoup plus élevée. Lorsque l'évaporation prend place à la surface, la chaleur latente nécessaire à la phase de transition est prise hors de la surface résultant d'un refroidissement. Lors de la formation des nuages, la vapeur d'eau condense et la chaleur latente est libérée dans l'atmosphère. Ceci conduit à un transfert net de chaleur de la surface à l'atmosphère, l'un des principaux moteurs de la 39 circulation atmosphérique. Le ratio des chaleurs sensibles et latentes est appelé Bowen ratio (Bo=Hs/LvE) . Le flux de conduction advient sur des surfaces solides telles que le sol et la glace et pour l'océan, il est lié à la dynamique de mélange des couches océaniques. aT p e _ aT Hs T = Pacpch aZ ; LvE 1 = P Läc1 âaZ ; Gs 1ù x â où cp est la capacité thermique spécifique de l'air, pa la densité de l'air, ch le coefficient de transfert turbulent de chaleur, T la température, z la hauteur, Lä la chaleur latente de vaporisation, e la pression de la vapeur de l'air, s le ratio des poids moléculaires d'air humide et sec, et le coefficient de transfert de vapeur d'eau et x la conductivité thermique du médium. L'inclusion des composantes de flux ondes-courtes et ondes-15 longues avec C9 la capacité thermique de la couche de surface permet de calculer le budget radiatif à la surface du globe. aTs CS at SSI+FS1-FS1-H 1-LvE1-G,1 A l'équilibre thermodynamique CSaTs/at = 0, ce qui permet de calculer la température de surface T9 pour déterminer son influence sur les océans et la biosphère. Le budget et ses 20 composantes sont calculés en W/m2/jour sur des mailles l°xl° et sont utilisés dans les modules océan et biosphère. Les modélisations météorologiques du module de transport atmosphérique (nuages, GES) sont intégrées dans le module énergétique pour la mise à jour des données. La mesure des 25 paramètres in-situ de flux de chaleur, de rayonnement et de GES lors des observations (Bloc 100) permet de valider les données du module avec des mesures réelles. Le calcul du budget énergétique est initialisé au début de l'holocène avec une périodicité de 50 ans et un calcul global sur la planète. 40 2953021 c. Module Océan Afin d'obtenir une absorption océanique plus précise, la méthode de mesure selon l'invention améliore le module actuel du Carbontracker en utilisant un module océan incluant 5 l'addition du dégagement de CO2 par évaporation, l'absorption par altération chimique et l'effet de tampon (Bloc 203, Fig. 9). Les océans sont les plus grands puits de carbone long-terme de part leur forte capacité de stockage et de redistribution 10 au sein du système. La méthode de mesure selon l'invention modélise l'absorption principalement par la dissolution du CO2 atmosphérique entre l'air et les océans par la différence de pression partielle de CO2 (pCO2), la vitesse du vent et la température des eaux, ceci modifiant son équilibre carbonate vers un état plus acide. Elle modélise également le dégagement selon les températures locales, l'activité biologique, la vitesse du vent, et la circulation océanique. Les échanges de CO2 sont calculés à partir du transfert de masse : Foce (t) = Foce (t) 1 ùFoce (t) 1 Une augmentation de CO2 dans l'atmosphère provoque une augmentation de la pression partielle, ce qui accroît le taux auquel il se dissout dans l'eau. La pression partielle de CO2 obéit à la loi de Henry, où Ko est le coefficient de solubilité du CO2 dans l'eau. [CO2]eau de mer Ko(S,T)• CCO2 (t) Selon la loi de Fick, l'absorption océanique J est déterminée en calculant le flux de diffusion de CO2, généralement décrit comme le produit de la vélocité de transfert du gaz k,, avec le gradient de concentration de CO2 de l'eau et de l'air marin. Il est également fonction de la profondeur zm de 41 l'océan qui tient compte des eaux peu profondes particulièrement près des cotes. kW J = ù ([CO2] saturé ù [CO2] eau de mer) Zm La variabilité spatiale et temporelle de l'échange air-mer de CO2 dépend ainsi de la distribution de la vitesse du vent, de la température, de la concentration de CO2 dissous et de la solubilité K0. Cette absorption est calculée dans la méthode en utilisant la différence de pression partielle de CO2 entre l'air et l'océan combinée avec une vitesse de transfert du gaz. Les niveaux de pCO2 sont déterminés en utilisant différentes configurations du modèle Princeton/GFDL MOM3, puis en divisant par une vitesse de transfert des gaz calculée selon les prévisions du modèle météorologique ECMWF (European Center for Medium-Range Weather Forecast), ERA40, intégrées dans le module de transport (Bloc 400). La vitesse de transfert du gaz est paramétrée comme une fonction quadratique de la vitesse du vent, en utilisant la formulation des vents instantanés. Le transfert air-mer est entravé par la présence de glace et les flux sont redimensionnés dans chaque grille par la fraction quotidienne de glace fournie par les prévisions de données du ECMWF. Une augmentation de concentration de CO2 dans l'atmosphère entraîne une augmentation de la quantité de CO2 absorbée par les océans, qui diminue progressivement. La méthode ajoute cette réduction de la capacité tampon du système par la méthode de Wolf-Gladrow (1994). La somme des espèces carbonates dissoutes est définie comme le total de carbone inorganique dissous [DIC]. Le CO2 dans l'océan forme un acide carbonique faible, le H2CO3 qui se dissocie sous la forme dominante de carbone inorganique stocké, l'ion bicarbonate HCO3 puis en ion carbonate CO32 avec K1 = [HCO3-][H+]/[CO2(aq)] et K2 = [CO32-][H+]/[HCO3-], les constantes de dissociation selon : 42 CO2(aq) + H20(1) -> H2CO3(aq) K1 H+ + HCO3- K2 2H+ + CO32 La concentration de CO2 dans l'océan est tributaire de la solubilité et de la pression partielle du CO2 dans l'atmosphère et le total [DIC] en solution est déduit à partir des constantes de dissociation et de la concentration en ions hydrogènes. K K K2) [DIC] = 1 CO2 = [CO2] + [HCO3-] + [CO32-] pco2(t)•Ko(s,T)(1 + [H+] + [H+]z) La fraction du flux de CO2 de l'atmosphère vers la couche de mélange qui va réagir est fonction du facteur tampon ç, qui est la fraction de changement de CO2 atmosphérique divisé par la fraction de changement de [DIC] après établissement de l'équilibre. Ç dépend de la température, du [DIC], de la salinité et de 1' alcalinité (Alk [HCO3-] + 2[CO32-] + [OH-] + [B(OH)4-] - [H+]) . (a [CO2] /a[DIC])eau [CO2] )atm / l [DIC] de mer Suivant la méthode de Trenbeth (1992), le flux est corrigé selon l'équation suivante : [CO2]eau de mer kw T ù [DIC] Zm ([CO2]saturé ù [CO2]eaudemer) La méthode ajoute également une absorption induite par l'altération chimique du CO2 où l'altération des roches sur la surface continentale consomme du CO2 atmosphérique pour produire de l'alcalinité. L'alcalinité est ensuite transportée par les fleuves et les rivières et précipitée en carbonate dans les océans, qui sont déposés par sédimentation. Cette altération est intégrée suivant le modèle global CDIAC DB1012 de résolution 1°xl° de Suchet et al. (1995) qui contient des estimations du flux net surface/atmosphère de CO2 (moles/km2/an) ainsi que du transport de bicarbonate (HCO3-) des rivières à l'océan. Le 43 modèle est fondé sur un ensemble de relations empiriques entre le flux de CO2 (Fco2alt) et le ruissellement sur les principaux types de roches affleurant les continents. L'absorption océanique est modélisée selon: oce(t) l~ [CO21eau zW F ([CO2] [CO2] mersaturé ù eau de mer CO2alt• m La méthode complète également le module océan avec l'évaporation de CO2. Combiné avec le carbonate de calcium insoluble, la réaction de dissolution du CO2 produit une solution de bicarbonate de calcium Ca(HCO3)2, qui s'accumule dans les océans. H2CO3(acide faible) + CaCO3(calcaire) - Ca(HCO3)2(solunon) A la surface des océans, le vent et la marée provoquent des vagues et des embruns, accompagnés de fines pulvérisations et de mousse. Le rayonnement solaire incident crée directement un rayonnement IR par chaleur latente, ce qui provoque une hausse des températures océaniques et l'évaporation de ces pulvérisations et de la mousse. La vapeur d'eau est en outre décomposée en précipitations et dégagement de chaleur latente. Chaleur latente + Ca(HCO3)2(soiution) - CaCO3(caica,re) + CO2(g0.z) +H20(vapeur) H20(vapeur) 4 H20(liquide) + Chaleur latente Dans l'océan, la forme dominante de carbone inorganique stocké est l'ion bicarbonate. En solution de Ca(HCO3)2, deux molécules HCO3 stockent une molécule de CO2 et le montant Kc de CO2 stocké qui peut être potentiellement libéré par évaporation est Kc = [HCO3 i/2 • Mco2 . La concentration de CO2 stockée en Ca(HCO3)2 et libérable par évaporation est également plus élevée que la concentration de CO2 en solution. Il est considéré que le taux moyen d'évaporation sur l'eau est plus important que sur terre ayant moins d'eau exposée. Les observations IR spectrales par satellite 44 indiquent également que presque tous les rayonnements IR émis par les océans proviennent de la vapeur d'eau au dessus de la surface (chaleur latente). A l'équilibre radiatif, il est supposé que presque 100% du rayonnement solaire absorbé par les océans provoque l'évaporation d'eau. Avec une fraction couverte d'environ 70,8%, la fraction d'énergie solaire qui évapore l'eau des océans est considérée équivalente à Fp.~ 0.7. La chaleur latente est déterminée à partir de l'équation du module énergétique (Bloc 202) selon: L,E1=Fp(SS1+FS1ùFSTùHSîùGS1) Il est également supposé que presque toute l'évaporation de la surface de l'océan est due à celle d'embruns ou de mousse des vagues, soit F9 1. Afin d'évaluer la libération de CO2 avec le taux d'évaporation E à température d'eau constante, E est proportionnelle aux rayonnements et à LV, la chaleur latente de vaporisation. La méthode évalue le dégagement par évaporation selon la proportion suivante: Fp E. Kc. Fs (Ss 1 +Is 1 ùIs 1 ùIls 1 ù Gs 1) . Lp . Kc. Fs L'absorption des océans est ainsi principalement induite par l'augmentation de concentration atmosphérique de CO2, leur capacité de tampon, l'altération chimique et l'évaporation qui est influencée par l'insolation solaire. Face (t) -- (Ss 1 +Is 1 ùIs 1 ùHs T ù Gs 1) . Fp . Kc. Fs - [CO21eau e mer kw ([CO2] ù [CO21 eau saturé eau de mer) m + FCO2alt) Le résultat du module océan est une cartographie des échanges océaniques 5°x4° en PgC/mois. La mesure des paramètres océaniques in-situ des observations marines (Bloc 104) (pCO2 océan, pression atmosphérique, salinité, température) permet de valider les données du module avec des mesures réelles. Le 45 calcul du flux océanique est initialisé au début de l'holocène avec une périodicité de 50 ans et un calcul global sur la planète. d. Module Biosphère Les échanges air/biosphère sont traités par le modèle biosphère JSBACH de Raddatz et al. (2007) (Bloc 204, Fig. 9) à la place du modèle CASA utilisé par le Carbontracker. JSBACH a notamment l'avantage de traiter les changements anthropiques de couverture de surface terrestre. L'absorption de CO2 est régie par la photosynthèse et le dégagement par les respirations et les perturbations. Les échanges de CO2 sont modélisés à partir du transfert de masse. Les plantes absorbent le CO2 lors de la photosynthèse par diffusion au sein des stomates qui sont les pores des feuilles et des tiges par lesquels le CO2 est pris et converti sous l'influence de rayonnement actif visible en hydrates de carbone. lumière du soleil 6 CO2 + 12 H20 > C6H1206 + 6 02 + 6H20 Les facteurs biotiques affectant la photosynthèse incluent la forme de croissance, le type de feuille, la voie photosynthétique (C3, C4) et la longévité. C3 est la photosynthèse de la plupart des plantes alors que C4 est une adaptation à des conditions arides avec une meilleure utilisation de l'eau. Les types de végétation sont modélisés en utilisant différents types fonctionnels de plantes (PFT, Plant Functional Types) et une représentation des différents biomes (forêts, arbustes, tourbières, prairies C3 et C4, marais, toundra, terres cultivées, glaciers...). La photosynthèse est modelée par les équations décrivant les flux d'émissions de CO2 et de vapeur d'eau au niveau de la feuille et redimensionnée à la canopée. L'indice foliaire est 46 calculé comme le ratio de surface supérieure totale des feuilles de la végétation divisé par la superficie de la surface sur laquelle elle croît. surface folière totale A= surface du sol La surface foliaire est calculée interactivement avec le 5 climat et les saisons de croissance et de décroissance sont modélisées selon: dA A \ ù= k(1ù IA ù pA dt Amax/ Où k est le taux de croissance avec k=0 pour NPP -0 et pA est le taux de décroissance (perte des feuilles). La lumière influence la photosynthèse de la 10 canopée et varie en fonction de son architecture lors de son passage. La fraction de surface ombragée fombre est déduite de la relation Beer-Lambert avec une orientation aléatoire des feuilles. A (ombre = 1 ù e L'albédo est calculé en fonction de l'indice foliaire variant 15 avec les saisons et l'enneigement. acanopée ù fombreaveg + (1 ù fombre) asurf ace La distribution verticale de rayonnement photosynthétiquement actif (PAR) est calculé en utilisant la radiation ondes-courtes du module énergétique (Bloc 202): I=I1+I1+Ss1 e-xnz Où I', I.1,, sont les rayonnements diffus ascendants et 20 descendants dans la canopée, S9 est le rayonnement incident sur le haut de la canopée, K est le chemin optique du rayon direct, fonction de la distribution des feuilles et AZ est l'indice foliaire mesuré du haut de la canopée. Le taux d'assimilation (A) de CO2 est modélisé comme le minimum du 47 taux de carboxylation (Jc) et de transfert d'électrons (Je). La méthode complète JSBACH avec une limitation par le taux de photosynthèse brute (Js) limité par la capacité à transporter les produits photosynthétiques pour les plantes C3 et la capacité de CO2 limitée pour les plantes C4 selon : A = min fic, Je, Js} J( Ci-r*) avec J = al Jm ' Je 4 Ciù2r, z z z JJm +a I où ci est la pression partielle de CO2 dans le chloroplaste, F la pression partielle de compensation photorespiratoire de CO2, Vm le taux maximum photosynthétique d'activité Rubisco, Km la constante de Michaelis-Menten, J le taux potentiel de transport d'électrons, Jm le potentiel maximum de photosynthèse limité par saturation de lumière et fonction de l'azote de la feuille, I le rayonnement solaire pénétrant la feuille et p est la pression de l'air de surface. La Production Primaire Brute (GPP), soit la quantité de carbone fixé par photosynthèse est calculée selon: GPP = A • A Les plantes relâchent du carbone, sous forme de respiration autotrophique, soit une oxydation des composés organiques en CO2 et H2O. 6CH20 + 602 - 6CO2 + 6H20+ énergie Le montant de carbone absorbé par les plantes et incorporé dans de nouveaux tissus végétaux est appelé Production Primaire Nette (NPP) et comprend les augmentations de la biomasse dans les feuilles, tiges, branches, racines, et organes de reproduction. La partie restante est perdue par respiration autotrophique (Ra) incluant respiration de maintenance (Rm) et de croissance (Rg). Rm sert à garder les tissus en vie et est une fonction de l'indice foliaire et de la respiration sombre. Rg sert à synthétiser de nouveaux = C1 -r. Je Vm C{ùKm J 0.5Vm Js=12.104*Vm*~` P 48 matériaux et est corrélée avec la croissance totale des plantes. R. = f J feul 111eRrespiration sombre. A NPP + Rg CCcoûts de construction ù NPP NPP = GPP ù Ra = GPP ù Rm ù Rg La plupart du carbone fixé par NPP retourne à l'atmosphère par respiration hétérotrophique du sol (Rh) qui est la décomposition des matières organiques par des bactéries et des champignons. Elle est fortement dépendante de la température du sol lorsque l'humidité est disponible et est calculée en utilisant le coefficient thermique Ql0 du sol, en intégrant la taille du réservoir de carbone avec l'humidité a et le temps de renouvellement t à 10°c. dC _ QTsol/1o°c C ùa. dt io z Les données écologiques suggèrent que la respiration du sol et la NPP sont corrélées positivement l'une avec l'autre. Une forte productivité primaire nette produit plus de litière, ce qui encourage plus de décomposition et de respiration, favorise la minéralisation et produit plus d'azote pour la croissance des plantes. Ce carbone du sol est partitionné en réservoir vert, réservoir de bois et réservoir de réserve avec différentes teneur en carbone, compositions chimiques et compositions en bactéries et champignons selon: = NPPpc ù FL1tière (Réservoir vert) dt dCW CW dt = NPP w ù (Réservoir de bois) zW d CR CR dt = NPPpR ù z (Réservoir de réserve) R Ces réservoirs libèrent le CO2 dans un réservoir avec un court (1 an) et un long (100 ans) temps de renouvellement, d CG 49 ceci conduisant au calcul de la productivité nette des écosystèmes (NEP) qui représente la quantité de carbone annuellement stockée dans la biosphère terrestre. dCF = CR k Tsol/1o°c CF dt TR + FLitière ù RF Ou RF = a . Q10 TF (Réservoir rapide ) dCs Cw k Tso1/1ooc Cs = dt T1,1, + (1 ù fF>A)RF ù Rs où Rs = oc. Q10 TS (Réservoir lent) NEP = NPPù Rh =A•71ùRmùR9ùRFùRs La production nette de la biosphère (NBP) est formalisée en intégrant les perturbations telles que les incendies et les changements dus à l'utilisation des sols. Le module incendie 10 (Bloc 205) prend en compte les incendies et est intégré à JSBACH. Les perturbations liées à l'utilisation anthropique des sols (ALCC, Anthropogenic Land Cover Change) sont calculées selon la méthode de Pongratz (2009) depuis le dernier millénaire. Les manipulations de la surface des 15 terres sont essentiellement causées par l'expansion ou l'abandon de la superficie agricole, incluant terres cultivées et pâturages qui modifient la couverture des sols et altère l'absorption de la biosphère. L'urbanisation influence également le climat via une demande alimentaire et 20 en énergies de substitution (biocarburants) croissante de la population, augmentant la demande en zones agricoles. Les forêts, herbes naturelles et zones arborées sont aussi touchées par cette expansion agricole. Le montant de carbone de ces perturbations émis directement dans l'atmosphère par 25 les trois réservoirs de végétation est traité comme suit: F,,A = (c ancien in ù ouveau) (fG,,A• CG,i + fW,,A• CG,i + fRDA• cru) ica Où fGDA, fWtA, et fRoA, sont les fractions de carbone ALCC relâchées dans l'atmosphère en fonction des trois réservoirs de carbone (vert, bois et réserve), canclen ù nouveau est la variation quotidienne de changement de la fraction de 50 couverture pour chaque type fonctionnel de plante qui perd la surface (a-) due au changement anthropique de couverture des terres et CG,±, CW,i et CR,I sont les densités de carbone des trois réservoirs. Pour la réallocation du carbone dans les réservoirs rapides et lents, le carbone des réservoirs vert et de réserve est transféré dans le réservoir rapide, tandis que le réservoir de carbone du bois est transféré au réservoir lent selon : FoF = Eicaù(ancien ù cnouveau) ((1 ù fGpA). CG,i + (1 ù fR>A). 4,3 (Réservoir rapide) F,,S _ E ica ancien ù cnouveau)(1 ù fwpA). CW,i (Réservoir lent) ù ù (, Le carbone de la végétation est perdu d'un PFT en raison de la diminution de sa surface, alors que les densités de carbone ne sont pas affectées. Le carbone perdu est ensuite transféré aux réservoirs de carbone respectifs du sol des PFT 15 en expansion, distribuées proportionnellement à leurs nouvelles fractions de couverture et les densités de carbone de PFT sont ajustées en conséquence. La surface foliaire, la respiration autotrophique et l'albédo sont également modifiés en fonction de la modification des types de couverture des 20 terres et de la végétation. La NBP représente le montant net d'absorption sur de longues périodes de temps en incluant les perturbations. Fbho(t) = -NBP = -(GPP - Rh- Rh - perturbations) Le résultat du module JSBACH est une cartographie des flux de la biosphère en Kg/m2/s avec un maillage l°xl°. La mesure des 25 paramètres in-situ des observations d'écosystèmes (Bloc 104) permet de valider les données du module avec des mesures réelles. Les mesures d'ALCC par les observations satellites (Bloc 101) sont intégrées afin de valider le changement de couverture des surfaces. Le flux de la biosphère est remonté 30 pour l'initialisation depuis le début de l'holocène avec une périodicité de 50 ans et un calcul global sur la planète. 51 2953021 e. Module Incendie Avant que l'homme n'ait utilisé les incendies pour défricher les terres et fertiliser les sols, la plupart des écosystèmes ont été soumis à des incendies naturels qui permettent de 5 rajeunir les forêts anciennes et d'apporter des minéraux. Le module incendie est intégré à JSBACH pour estimer le flux des perturbations Ffeu(t) liées aux incendies (Bloc 205, Fig. 9). Les résultats sont calculés à partir de la base de données Global Fire Emission Database GFEDv2 de Randerson et al. 10 (2007) et les données se composent de mesures 1°xl° des superficies brûlées maillées, des charges de carburant, de l'efficacité de combustion, et des émissions de GES incluant le CO2 (Kg/m2/mois). Après avoir estimé globalement les zones brûlées, la végétation évoluant saisonnièrement, les stocks de biomasse des sols de JSBACH sont brûlés en fonction de l'estimation de la superficie brûlée et convertis en émissions atmosphériques de GES en utilisant les estimations d'efficacité de combustion, de complétude et de charges de carburant. Les données du GFEDv2 sont calculées depuis 1997. f. Module Fossile Un nouveau module fossile (Bloc 206, Fig. 9) est intégré à la place du module fossile du Carbontracker pour mesurer les émissions globales venant de la production de combustibles fossiles afin de s'assurer que l'ensemble de la production soit prise en compte. Au cours des deux derniers siècles, suivant la révolution industrielle et l'augmentation de la population mondiale, la combustion d'énergie fossile est devenue la plus grande source anthropique de CO2 pour notamment la production d'électricité, le transport, le chauffage et les procédés industriels. Le taux de rejet du CO2 fossile est calculé en additionnant les contributions du 52 charbon et du pétrole à partir des inventaires de production de l'Energy Information Administration (EIA) depuis 1990 : Fff(t) = Fcharbon(t) + Fpétrole (t) Le charbon est composé de près de 100% de carbone et en considérant que presque toute cette masse de carbone entre dans l'atmosphère, les émissions mondiales de charbon en CO2 sont obtenues en additionnant les productions des principaux producteurs: Production de Charbon (millions tonnes/année) x Fcharbon(t) charbon(t) Mc Le pétrole se compose de chaînes carbonées avec pour chaque atome de carbone environ 2 atomes d'hydrogène joints, soit environ 86% de carbone. La production annuelle mondiale de pétrole est obtenue en additionnant la production journalière en barils. Il est à noter que tout le carbone n'est pas rejeté dans l'atmosphère et qu'une partie est utilisée pour produire de l'asphalte et des résines. Toutefois, cette partie non-atmosphérique est très probablement compensée par la production de combustibles fossiles non déclarée, en particulier celle du charbon. Fpe -trole(t) Production de pétrôle (millions tonnes/année) x 86% x Mcoz Le calcul des données du module fossile est remonté depuis le début de l'ère industrielle jusqu'à 1800 selon les estimations de productions de Etemad et al. (1991) afin d'obtenir une mesure de flux anthropiques fossiles à l'échelle planétaire sur l'ère industrielle en T/an. 3. Module d'inventaires ascendants Suite à la modélisation de l'évolution des flux effectuée par le module d'échange, la méthode de mesure selon l'invention effectue ensuite une modélisation des émissions anthropiques Mc 53 hebdomadaires (Fig. 9, bloc 301), ou inventaires ascendants, au moyen d'un module d'inventaires ascendants. Ces inventaires ascendants sont utilisés comme des estimations à priori des flux anthropiques et proviennent du modèle EDGAR 4.0 http://edgar.jrc.ec.europa.eu/ de la Commission Européenne, Joint Research Centre (JRC) et du Netherlands Environmental Assessment Agency (PBL) (Bloc 300, Fig. 9). La répartition spatiale actuelle et en développement des émissions dans la série de données est effectuée en fonction des émissions annuelles déclarées pour la période 1970-2005 et incluent GES, gaz acidifiants et particules (CO2, CH4, N20, HFCs, PFCs, SF6, CO, NMVOC, S02, NOX, NH3, PM2.5, PM10, OC, BC, HCFC, CFC) (Fig. 13). Les inventaires ascendants sont déterminées en utilisant les données spécifiques de chaque pays qui sont organisés en régions du monde et les catégories principales de sources incluent: énergie, procédés industriels et usage de produit, agriculture, utilisation des terres, leur changement et la forêt, les déchets et autres sources anthropiques. Les inventaires ascendants d'émissions sont calculés pour chaque secteur et pays avec un facteur d'émission basé sur la technologie. Les données des paramètres de l'équation suivante sont inclues pour chaque combinaison pays/secteur considérée: Emissionx, (an) = 1{AD,S(an) * TECHAD cS(an) * EOPAD,c,s,TECH(an)} * EFAD,c,s,TECH «1 ù REDEOP) où (x) est le composé,(c) le pays, (s) le secteur, (an) l'année, (AD) la donnée d'activité en TJ/an (ex : charbon utilisé dans un pays pour la production de chaleur), (TECH) la technologie, (EOP) le pourcentage des technologies qui sont contrôlées par des mesures de réduction, (EF) le facteur 54 d'émission incontrôlé par secteur et technologie en KT/TJ et (RED) la réduction en pourcentage du facteur d'émission incontrôlé par la mesure de réduction. Ces données sont en partie fondées sur les rapports d'inventaires, les rapports industriels, les directives d'inventaires et la littérature scientifique. Les inventaires ascendants d'émissions par pays sont allouées sur une grille spatiale de haute résolution 0.l°x0.1° (100km2) qui peut être agrandie à de plus faibles résolutions de 0.5°x0.5° à l°xl° et sont intégrées aux différentes résolutions en utilisant les techniques de GIS (Geographic Information System) pour la conversion, le ré-échantillonnage et l'agrégation. Chaque grille spatiale est liée avec la grille du pays de référence construite sur la base de données des «populations maillées du monde» (GPWv3). La répartition des émissions par secteur pour chaque cellule d'un pays est effectuée selon: Fff(x, y, t) = Emission(xj, yj)AD,c,an = Emission(C)AD,an x EIndicateur(x, y)AD,c,an Le couple (xi, yi) représente le coin inférieur gauche de chaque maille 0.1, (AD) la donnée d'activité (ex: gaz naturel d'une centrale électrique), (C) le pays, (an) l'année et l'indicateur de grille spatiale. La méthode modélise en outre la variabilité temporelle de ces inventaires (Bloc 302) selon le modèle de Gurney (2009) pour les Etats-Unis qui est étendu sur le reste du monde avec l'équation de Gurney (2005) avant que les saisonnalités de Gurney (2009) ne soient disponibles à l'échelle mondiale: cos( 7t(m ù 1)) Fm=Fm+AkF,nsin0j l 12 où Fi7im, est le nouveau flux sur chaque cellule de la grille en utilisant, (i) l'index de longitude, (j) l'index de latitude et (m) l'index du mois. Fjm est le flux original et (Ak) est le facteur d'amplitude qui représente le pourcentage Indicateur (xi, yi)AD,c,n 55 d'accroissement des émissions fossiles d'origine. Ces saisonnalités permettent une décomposition temporelle des inventaires effectuée sur une base hebdomadaire. La méthode utilise ces inventaires dans le module de transport avec une résolution l°xl° en Kg/m2/semaine et ces résultats sont redimensionnés sur l'année en cours proportionnellement aux totaux du module fossile pour obtenir la même base étalonnée. 4. Module de transport Suite à la modélisation des émissions anthropiques, la méthode de mesure selon l'invention effectue ensuite une modélisation du transport atmosphérique (fig.9, bloc 401) au moyen d'un module de transport. La méthode utilise pour simuler les vents et le climat, le modèle de transport TM5 décrit par Krol et al. (2005) piloté par le modèle de prévisions météorologiques du European Center for Medium range Weather Forecast (ECMWF). Le transport du CO2 dans l'atmosphère (Bloc 400, Fig. 9) permet de lier les observations de CO2 des différentes couches de l'atmosphère (Bloc 100) aux flux de CO2 à la surface terrestre (Bloc 200, Bloc 300). Les tempêtes, complexes nuageux et conditions météorologiques sont à l'origine de vents qui transportent le CO2 et l'influence de dégagements, d'absorptions et d'événements locaux peuvent avoir des impacts à des endroits éloignés. Ce modèle complexe en 3D simule les concentrations de CO2 dans l'atmosphère à partir des flux en utilisant les champs de prévisions météorologiques de l'ECMWF. C'est un modèle de zoom atmosphérique à grilles imbriquées dont les régions pour lesquelles les simulations de haute résolution sont souhaitées peuvent être imbriquées dans une grille couvrant le domaine global. Il a l'avantage d'effectuer des simulations de transport avec un focus régional sans la 56 nécessité de mettre en place des conditions limites telles que dans d'autres modèles. Ceci permet aux mesures en dehors du domaine zoomé de contraindre les flux régionaux dans l'inversion et l'assimilation des données, et de s'assurer que les estimations régionales sont cohérentes avec les contraintes globales. TM5 fonctionne à une résolution horizontale de 6°x4° et le zoom descend jusqu'à une résolution l°xl° (Europe, Amérique du Nord, Amérique du Sud, Asie, Australie), zones qui sont imbriquées dans des régions 3°x2° pour assurer une transition en douceur entre les différents domaines. Le TM5 simule séparément l'advection, la convection et la diffusion verticale dans la couche limite planétaire et la troposphère. TM5 s'exécute sur un pas de temps de trois heures et les processus à plus fine échelle sont répétés toutes les 10 minutes (découpage et résolution raffinée dans les grilles imbriquées). La résolution verticale est de 25 niveaux hybrides de pression sigma, inégalement espacés avec plus de niveaux près de la surface. Le TM5 modélise les concentrations en ppm/s jusqu'à l'échelle l°xl°. La mesure in-situ des paramètres météorologiques par les observations atmosphériques, marines et d'écosystème (Bloc 300) permet de valider la modélisation du TM5 et d'affiner sa paramétrisation. 5. Module d'inversion et d'assimilation de données A la suite de ces modélisations, la méthode de mesure selon l'invention calcule ensuite des flux finaux au moyen d'un module d'inversion et d'assimilation de données (fig.14). Ce module (Bloc 500, Fig. 14) utilise les observations pour déduire la distribution spatiale 1°xl° des flux terrestres, océaniques et anthropiques. L'atmosphère est représentée par un vecteur d'état représentant le flux net surface-atmosphère 57 à déterminer ou 2\.r et off représentent un ensemble de facteurs scalaires linéaires appliqués aux flux et estimés chaque semaine. Fco, (x, y, t) = Àr. Fm, (x, y, t) + Àr . Foce (x, y, t) + Àf f . Fff (x, y, t) + Ffeu (x, y, t) Ce module combine une nouvelle inversion de synthèse (fonction de Green) (Bloc 501) simple et robuste pour l'analyse aux échelles régionales (planète, continents, régions continentales et pays), suivie par l'ensemble d'assimilation de données (Bloc 502) pour estimer les flux avec un niveau de détail plus précis sur les pays industrialisés ou les mesures 1°xl° sont souhaitées (ex : Europe, Amérique, Asie). Une répartition géographique des surfaces du globe basée sur le modèle de décomposition DB1016 de Li (1990) est ajoutée pour agréger les émissions des mailles 1°xl° par régions et pays. Les champs de concentration et de flux sont amorcés avec une moyenne d'estimation dans le vecteur initial d'état à partir des équations du module d'échange et leur lancement depuis l'holocène permet d'obtenir une base étalonnée au temps des premières observations. Dans l'inversion et l'assimilation, le terme « à priori » ou « p » est réservé aux flux créés initialement et sont fixés, « a » se réfère aux quantités analysées des étapes précédentes, « b » ou « d'arrière-plan » sont les flux résultants et contiennent des informations tirées au fur et à mesure des analyses de cycles antérieurs et les flux « à posteriori » sont les flux finaux. Les flux à priori sont donc redimensionnés par ce module dans lequel les observations sont utilisées pour déduire les flux à postériori. La méthode définit un critère d'optimalité avec une fonction de coût qui est minimisée afin que le processus de modélisation soit comparable à celui de l'atmosphère sur une échelle moyenne de temps et d'espace. 58 Les estimations sont affinées au fur et à mesure des observations et les valeurs uniques de Àr et de off résultent de la plus petite différence aux moindres carrés entre observations et modélisations. Pour la classification des surfaces, chaque facteur 2\.r est associé à une région (r) particulière du domaine global. L'océan est divisé en 30 grands bassins océaniques englobant les caractéristiques de circulation et la biosphère est divisée en fonction des types d'écosystème ainsi que de leur situation géographique. Chacune des 11 régions TransCom de la terre contient un maximum de 19 types d'écosystèmes et l'approche conduit à un nombre total de r=11xl9+30=239 facteurs d'échelle 2\.r optimisables chaque semaine sur le globe. Même avec un seul paramètre Àr disponible à l'échelle, chaque grille 1°xl° a un flux F(x,y,t) différent selon la moyenne des flux modelés par les modules océan et biosphère. Les flux liés aux incendies ne sont pas redimensionnés. La méthode ajoute le paramètre Àff au Carbontracker afin de redimensionner les résultats du module d'inventaires ascendants (Bloc 300). Chaque off est apposé sur les émissions à priori et est optimisable sur chaque grille 1°xl° en fonction de la corrélation entre les observations effectuées et les flux anthropiques d'origine. Ceci est motivé par le fait que la distribution géographique des émissions à priori est bien connue dans le module d'inventaires ascendants. La qualité du module d'échange et la densité des observations avec notamment la couverture, la résolution et la précision des observations satellites, permettent de contraindre spatialement les flux naturels et de distinguer la composante anthropique. De plus, l'analyse de la variabilité spatiale et temporelle entre les estimations individuelles des sources de flux permet de 59 déduire la composante anthropique qui est obtenue des différences modélisées entre les océans, la biosphère et les incendies. Les flux naturels ont une variabilité plus forte et corrélée, alors que les flux anthropiques ont des variations plus faibles et ne sont pas corrélés. De larges fluctuations interannuelles reflètent les échanges naturels des écosystèmes terrestres induits par les évolutions météorologiques et climatiques à grande échelle et ne sont pas expliquées par la variabilité des émissions fossiles. Ces corrélations permettent d'évaluer lorsque les émissions sont liées à des sources anthropiques ou naturelles et l'échantillonnage géographique affine cette distinction en réduisant les échelles progressivement. Dans le module, le principe de comparaison entre les observations, le module d'échange et le module d'inventaires ascendants est le même pour les observations aériennes, atmosphériques, d'écosystèmes et marines que dans le Carbontracker. Pour comparer les observations de concentrations aux flux, le TM5 prend une distribution initiale de concentrations de CO2 et la propage vers l'avant dans le temps en utilisant les prévisions météorologiques tout en altérant les concentrations à la surface par les flux à optimiser. Cette distribution est ensuite, comparée aux temps et lieux des observations. La comparaison des observations de flux aux flux modélisés est effectuée tel que pour Ameriflux dans le Carbontracker. Les modélisations de flux sont intégrées aux temps et lieux des observations de flux et l'inversion et l'assimilation permet de minimiser les différences entre les flux modélisés et ceux observés. La méthode ajoute toutefois une modification pour les données satellites XCO2 de GOSAT suivant la méthode de Feng et al. (2009). Les champs de concentration en 3D sont modelés par le 60 TM5 à partir du module d'échange et des inventaires ascendants puis intégrés aux temps et lieux des observations pour chaque mesure en utilisant leurs orbites. Des fonctions de probabilité de densité (PDFs) des nuages et des profondeurs optiques des aérosols (AODs) sont dérivées pour récupérer les données à ciel-clair et la comparaison des flux de surface avec les données XCO2 est effectuée ensuite en appliquant des noyaux moyennant (Column Averaging Kernels) pour prendre en compte la sensibilité verticale de chaque satellite et cartographier les profils de concentration l-D CO2 aux colonnes moyennes observées. Pour la comparaison avec les observations de XCO2 de SCIAMACHY, la méthode s'appuie sur Buchwitz et al. (2005a) qui permet sur le même principe, l'obtention des colonnes moyennes par l'utilisation de noyaux moyennant. • Inversion de synthèse (Fonction de Green) L'inversion de synthèse est effectuée suivant la méthode d'Enting (2002) pour optimiser les émissions des grandes régions. Seuls les totaux régionaux sont calculés par la somme des flux des mailles l°xl° et la fraction atmosphérique de CO2 est représentée comme une combinaison linéaire de lancements de modèle pour les émissions des différentes régions et des différentes semaines. La liste des symboles utilisés est la suivante : Symbole Nom Unité Dimension s Vecteur d'état kg/m2/s [s] .......................................................................... .......................................................................... ................................................ . z Vecteur d'estimations à priori de s kg/m2/s [s] .......................................................................... .......................................................................... .......................................................................... ..... . Vecteur contenant les fractions molaires [m] m modelées ppm .......................................................................... .......................................................................... .......................................................................... ..... . c Vecteur contenant les observations ppm [m] .......................................................................... .......................................................................... .......................................................................... ..... . Matrice de covariance inverse des R (ppm) 2 [m] x [m] données d'observation .......................................................................... .......................................................................... .......................................................................... ..... . 61 Matrice de covariance inverse des P (kg/m2/s)2 paramètres à priori .......................................................................... .......................................................................... .......................................................................... ..... . Opérateur d'observation (TM5) ppm ppm [m x m] .......................................................................... .......................................................................... .......................................................................... ..... . G Matrice de Green kg/m2/s ppm [m x s] La forme générale d'équation de transport du CO2 décrit le taux de changement avec le temps m(r,t), de la concentration atmosphérique modelée de CO2 à un point r et au temps t en fonction de la source locale s(r,t) à chaque point et de l'opérateur de transport J-C modélisant la contribution due au transport du gaz d'autres endroits et est sujette à des conditions limites spécifiques. a atm(r,t) = s(r,t) +3-C[m(r,t)] Elle définit une relation linéaire entre les concentrations m(r,t) et les sources s(r,t) afin de résoudre l'équation des concentrations observées (c). La solution de cette équation par la fonction de Green avec des conditions limites spécifiques (t, t') et (r, r') est exprimée selon: m(r,t) = mm(r,t) +J d3r' J G(r,t,r',t')s(r',t')dt' Où mo(r,t) décrit l'état initial, m(r,to), et la solution est calculée pour être équivalente à la solution d'intégration de am(r,t)/at. L'inversion de synthèse discrétise les sources en terme de processus µ comme des facteurs de dimensionnement sg inconnus, multipliés par la distribution de sources spécifiée 6µ (r,t), appelées « fonctions de base »: sl,(r,t) = sl,6l,(r,t) Ceci permet de relier les sources au module d'échange et au module d'inventaires ascendants pour chaque processus: 62 s(r, t) = sp(r, t) = sj,6j,(r, t) u u = Àr• Fbio (x, y, t) + i1 r . Face (x, y, t) + i1 ff . Fff (x, y, t) + Ffire (x, y, t) Puis à la fonction de Green tel que: m(r,t) _ E,sj,Gj,(r,t) avec Gj,(r,t) = f d3r' f dt'G(r,t,r',t')c (r',t') L'analyse formelle de la fonction de Green est exprimée de la relation générique discrétisée: cj = Gjj,sj, +s~ = mj+si u où cj est un élément de concentration observée, sµ est la source, mj est le modèle de prédiction de concentration pour cet élément, Ej est l'erreur dans cj et Gdµ est une discrétisation de G(r,t,r',t') reliant les concentrations aux sources. Pour chaque fonction de base 6(r,t), l'intégration du modèle de transport numérique produira une réponse G(r,t). Gdµ sont les réponses pour l'observation j d'une source définie par la distribution o (r,t) et les sources sont estimées en utilisant cette équation pour ajuster les coefficients sp soit 2\.r et Àff. Cette fonction qui utilise des composantes prédéfinies o (r,t) est appelée calcul de synthèse car l'estimation de la source est synthétisée à partir des composantes prédéfinies. De ces intégrations, les valeurs spatiales et temporelles spécifiques qui correspondent à chaque observation j peuvent êtres extraites pour produire la matrice Gdµ. Pour l'inversion, l'approche de Bayes est utilisée, incluant la connaissance des inventaires à priori avec la fonction de coût J. La solution optimale de vraisemblance de variable inconnues de flux de CO2 du vecteur d'état s est trouvée en minimisant : = [GTRG + P]ù1 [GTRc + Pz] 63 La matrice de covariance de s est [GTRG+P]-1 et les résultats de chaque fonction de base sont comparés avec les observations moyennées aux valeurs journalières. La solution du modèle final est calculée comme une superposition linéaire des modèles lancés sur les différentes régions et les différentes semaines. Une fois finalisée, les valeurs uniques de 2\.r et de off dans le vecteur d'état sont utilisées dans l'ensemble d'assimilation de Kalman. • Ensemble d'assimilation de Kalman (EnKF) L'assimilation est effectuée selon la méthode de Peters (2005) et progresse avec deux étapes distinctes par cycle, celle d'analyse et celle de prévisions. La première sert à trouver l'état du système qui est optimalement cohérent avec les observations et la deuxième décrit l'évolution dans le temps de cet état optimal lorsque de nouvelles observations sont disponibles. A partir de ce moment, l'état de prévisions sert de première approximation, ou « d'arrière plan » pour la prochaine étape d'analyse. Ces étapes sont ensuite combinées dans un cycle complet d'assimilation suivi par de nouvelles observations. La liste des symboles utilisés est la suivante : Symbole Nom Unité Dimension x Vecteur d'état kg/m2/s [s] .......................................................................... .......................................................................... .......................................................................... ..... . xi Déviations du vecteur d'état kg/m2/s [s] .......................................................................... .......................................................................... ...................... . p Matrice de covariance d'état (kg/m2/s)2 [s x s] X Matrice de déviation d'état kg/m2/s [s x N] .......................................................................... .......................................................................... .......................................................................... ..... . y° Vecteur d'observation ppm [m] Matrice de covariance R d'erreur d'observation ppm2 [m x m] .......................................................................... .......................................................................... .......................................................................... ..... . Opérateur d'observation (TM5) kg/m2/s ppm [s] [m] .......................................................................... .......................................................................... .......................................................................... ..... . Opérateur linéaire H d'observation (matrice) kg/m2/s ppm [s x m] .......................................................................... .......................................................................... ................................................................. . kg/m2/s Modèle dynamique kg/m2/s [s] [s] 64 Concentrations de CO2 grille TM5 x CO2i(x,y,z,t) d'arrière plan ppm N .......................................................................... .......................................................................... .......................................................................... ...... . Pour la première étape, la solution optimale de vraisemblance de variables inconnues de flux du vecteur d'état x est trouvée comme équilibre de la fonction de coût (J) suivante: Observations Paramètres arrièreùplan 1(x) = (yo û (x))TR-1(y° û (x)) + (x û xb)TP-1(x û xb5 L'opérateur J-C qui convertit l'état du modèle à l'espace d'observations échantillonne le vecteur d'état x et retourne un vecteur J-f(x) à comparer aux observations. Le vecteur d'observation y° contient les fractions molaires de CO2 observées moins celles d'arrière plan CO2(x,y,z,t) afin de comptabiliser les variations. Le vecteur d'état x qui minimise J est décrit par: xt = xb + K (yt - 3-C(e)) avec Pt = (1 - KH)Pb Où H est la forme linéaire de If, t est le temps et K, la matrice de gain Kalman, définie par: K = (PbHT)(HPbHT + R)-1 En créant un ensemble de N champs de flux de CO2 qui ont une moyenne z et qui couvre la structure de covariance P, les flux optimisés sont trouvés en utilisant une série d'observations de CO2 avec la covariance R en lançant le modèle de transport atmosphérique J-C, N fois vers l'avant et en l'échantillonnant en cohérence avec les observations. Pour créer les statistiques d'ensemble, la méthode étend le nombre de membres d'ensemble du Carbontracker à préférablement N=300. La structure de covariance P décrit l'ampleur de l'incertitude sur chaque paramètre, ainsi que leur corrélation dans l'espace et l'information dans P, d'arrière plan et analysée, est représentée dans les dimensions N d'un 65 ensemble de vecteurs d'états xi composé d'un état moyen () et de ses déviations (x';) tel que x;=+x'; où xi est fonction des paramètres 2\.r et off à optimiser. Les déviations x', sont créées telles que l'ensemble normalisé de déviations définis les colonnes de la matrice X qui est la racine carrée de la matrice de covariance P=XXT suivant: X N 1 (x ù , x2 ù X, ... , Xn ù )) L'ensemble des vecteurs d'état défini la Fonction de Densité de Probabilité Gaussienne (PDF) de x avec la covariance P. Le filtre de racine carré d'ensemble (EnSRF) du Carbontracker selon Whitaker et al. (2002) sert à calculer l'ensemble analysé et les paquets d'observation appartenant à une étape de temps du filtre sont traités un à un. K est calculé selon les approximations suivantes: HPHT N 1 1 (~C(x ' 1), 3-C(x' Z), ... , (x n)). (x(x ' 1), (x' 2), ... , 3-C(x n))T PHT N 1 1(x'1,x'2,...,x n).(3f(x'1),3-C(x' 2), (x n))T Chaque entrée N définie une colonne de vecteurs d'états d'ensemble ou d'ensemble de valeurs de CO2 modelées. K cartographie linéairement les quantités observées aux éléments de vecteurs d'état comme une moyenne sur l'ensemble des membres. Le vecteur moyen d'état et ses déviations sont mis à jour selon: xt = xb + K (yc û H(xb)) ; x= x' û k3-C(x' b) avec k = K. a = K (1 + ,\/R/(HPbHT + R)) La moyenne analysée et l'état d'ensemble d'une observation sert d'état d'arrière-plan de la prochaine. Elles iront également dans le calcul des prochaines observations de la matrice K. Le vecteur des concentrations échantillonnées est mis à jour de façon similaire au vecteur d'état en utilisant l'information moyennée d'ensemble de K. Chaque concentration 66 modelée d'une observation m a assimiler J-C(xt)m et ses déviations J-C(x'i)m sont mises à jour selon . (Xt )m = (Xt)m + HmK (Yt - (Xb) ; (X' )m = - Hmk3f(X; b) Après la mise à jour de l'ensemble de valeurs de CO2 modelées, cet algorithme se poursuit sur la prochaine observation jusqu'à ce que l'ensemble des observations soit traité pour aboutir à l'ensemble final analysé. Dans la deuxième étape, le modèle dynamique décrit l'évolution du vecteur d'état dans le temps. Il apporte une première estimation avant que de nouvelles observations soient introduites et est appliqué à la moyenne des valeurs de Àr et off selon: ~îbtt = ~Ct-2+3bt-1+1~p et /lbt = 1.ff,t_2+1bt-1+1p 3 Les valeurs de 2\.r et off pour un nouveau pas de temps sont choisies comme une combinaison entre les valeurs optimisées à partir des deux pas de temps précédents et une valeur fixe précédente. Ce lissage sur trois étapes réduit les variations de prévision de Àr et Àff sur le temps. L'inclusion de Àp et Àp à priori agit comme une régularisation afin que les paramètres reviennent à la valeur prédéterminée sans observations. Àp et Àp sont initialisés à 1. Pour les cycles d'assimilation, le vecteur d'état contient les estimations de flux de plusieurs étapes de temps, chacune correspondant à la moyenne d'une semaine. La figure 15 présente 3 cycles d'assimilation avec 5 semaines de flux composant le vecteur d'état indiquées par xi (0,..,4), ou (0,..,4), défini le nombre de fois ou une semaine particulière de flux a été estimée sur la base des observations des cycles précédents. i se réfère à chaque membre d'ensemble individuel et chaque case ombrée représente un ensemble [i=l,..,N] de flux de surface du globe. Les cases 67 claires montrent les flux d'arrière plan et les foncées, les flux postérieurs et le cycle s'exécute selon: (1) Le TM5 est lancé en avant à partir des champs de concentration initiaux d'arrière plan, de CO2(x,y,z,t) à CO2(x,y,z,t+5) forcés par les flux d'arrière plan xi (0,..,4). Il extrait les fractions molaires de CO2 aux temps et lieux des observations afin de construire un ensemble de concentrations modelées sur chaque site. (2) Les équations de mise à jour du vecteur moyen d'état et de ses déviations sont résolues pour donner un ensemble analysé de flux pour chaque élément du vecteur d'état et chaque semaine. (3) L'ensemble des flux finaux dans x'(5) ne sera plus estimé dans le prochain cycle. Il sera incorporé dans CO2i(x,y,z,t+l) en lançant le modèle TM5 une semaine en avant en débutant de CO2i(x,y,z,t) forcé avec l'ensemble final des flux x±(5) . (4) Chaque vecteur d'état analysé devient celui d'arrière-plan pour le prochain cycle (flèches verticales claires). Un nouveau flux moyen d'arrière-plan est créé pour aller dans z(0) par propagation avec le modèle M (flèches horizontales sombres). (5) Un nouvel ensemble de N déviations de flux xi(0)est obtenu de la structure de covariance d'arrière plan pour représenter le PDF Gaussien autour de la nouvelle moyenne de flux z(0) . (6) De nouvelles observations y° sont lues et un nouveau cycle commence. Une fois finalisé sur une année d'observation, les résultats finaux, autrement dit les flux finaux, sont les valeurs optimisées des paramètres Àr et Àff du vecteur d'état avec une 68 cartographie des flux naturels et anthropiques en Kg/semaine de résolution 1°xl°. 6. Module de pondération La méthode de mesure selon l'invention améliore ensuite le Carbontracker en ajoutant un nouveau module de pondération, qui apporte la validation que les résultats agrégés de flux anthropiques reproduisent fidèlement ceux de la planète, des continents, des régions continentales, des états et pays et fournit une vérification indépendante de la fiabilité des données scientifiques. Le module (Bloc 600, Fig. 16) est fondé sur la théorie des jeux dans les relations internationales (Luterbacher et al. 2001) et consiste en une modélisation macro-économique des activités de production des secteurs économiques (énergie, procédés industriels, usage de produits, agriculture, utilisation des terres, leurs changements et la forêt, déchets et autres sources) de chaque pays et de son utilisation en énergie fossile. Dans le principe, la majeure partie de l'énergie mondiale est produite par la combustion d'énergie fossile et lorsque la consommation augmente les émissions de CO2 suivent, et ce, même dans les pays produisant de l'électricité sans carbone ou les énergies fossiles jouent un rôle important dans les activités de production. Pour calculer les émissions, les fonctions de production représentent la valeur ajoutée de sortie de chaque secteur économique par zone et incluent énergie et mélanges de carburants utilisés (charbon, pétrole, gaz naturel, électricité). Le prix relatif de chaque carburant modifie le mélange de combustibles utilisé et suit le principe qu'un carburant plus cher est habituellement substitué par un carburant moins cher (Bloc 601) mais ce qui ne diminue pas nécessairement les émissions si la demande devient plus forte dans une énergie avec une proportion plus 69 importante de carbone. Les niveaux d'émissions sont ainsi définis à partir des prix relatifs des mélanges de combustibles utilisés (Bloc 602) et de la demande énergétique d'un processus de production donné (Bloc 603). Un facteur de saisonnalité module la consommation d'énergie afin d'obtenir des émissions calculées selon la saisonnalité du module d'inventaires ascendants (Bloc 604). Les corrections liées à l'efficience énergétique sont ajoutées car l'utilisation d'énergie diminue proportionnellement aux intrants d'un processus de production dans le temps (Bloc 605). La demande énergétique totale d'un processus de production avec les prix relatifs des carburants définit le mélange de combustibles, l'énergie utilisée et in fine les émissions totales du processus (Bloc 606). L'agrégation de l'ensemble des processus d'un pays et des régions donne les émissions totales de ces zones (Bloc 607). Une représentation du marché des émissions sur les zones réglementées telles que l'Europe permet de tenir compte des effets de changements technologiques et de la réduction des niveaux d'émission. Ces marchés sont différents des marchés d'échanges classiques car les actifs environnementaux tel que le contenu de carbone de l'atmosphère est le même pour tous. Ils ont des propriétés physiques indépendantes des institutions économiques car ce sont des biens publics qui ne sont pas rivaux en consommation mais produits de manière privée. La méthode effectue cette représentation suivant le modèle de marché des biens publics produits par le privé selon Chichilnisky et al. (2000). L'analyse est fondée sur les liens d'équité et d'efficience de ces marchés et deux facteurs importants sont pris en compte, les quotas d'émission de chaque pays et les prix du carbone. L'offre ou la demande en certificats de carbone est générée en fonction du niveau des quotas, ce qui établit les prix et induit les 70 changements technologiques réduisant le CO2 (Bloc 608). Les résultats du module de pondération, les flux finaux pondérés, sont les totaux d'émissions en TCO2/semaine. Les résultats du module d'inversion et d'assimilation de données sont ensuite corrigés en modulant les facteurs scalaires Àff pour obtenir une justesse économique à chaque sous-échelle géographique jusqu'aux niveaux nationaux, fournissant ainsi des flux finaux pondérés. 7. Module de géocodage La méthode de mesure selon l'invention améliore enfin la précision des mesures par rapport au Carbontracker en ajoutant également un nouveau module de géocodage (bloc 700, Fig. 17). Une fois validés, les flux de CO2 des mailles 1°xl° (kg/m2/semaine) et résultats de la méthode sont transmis vers le module de géocodage comprenant un système de coordonnées GIS (Geographic Information System)(bloc 702) permettant de géocoder les résultats et notamment, les données du module d'inventaires ascendants provenant d'EDGAR 4.0 (bloc 300) telles que l'emplacement géographique des installations énergétiques et manufacturières, des réseaux routiers, des routes commerciales, des densités de population humaines et animales et l'utilisation des terres agricoles. La répartition géographique des surfaces du globe par pays et régions est également effectuée selon le modèle de décomposition de Li (1990). Afin de modéliser les flux jusqu'à l'échelle de l'installation, les inventaires d'émissions anthropiques de CO2 déclarés dans le modèle EDGAR 4.0 du module d'inventaires ascendants sont distribués spatialement par maille 0.1°xO.1°. Puis, la méthode applique des coefficients correcteurs proportionnels (bloc 704) sur chacune des mailles 0.1°xO.1°. Ces coefficients correcteurs kid sont calculés 71 proportionnellement au total de chaque maille l°xl° obtenu à partir des flux finaux pondérés par le module de pondération. Les coefficients k±,j, sont déterminés linéairement sur chacune des mailles 0.l°x0.1° en fonction des inventaires de CO2 déclarés par installation et ce calcul de proportionnalité est effectué selon l'équation suivante : io io Indicateur (xi, yi)AD,c,an Fff(x,Y)1°X1° _ x Emission(C)AD,an x =i EIndicateur (x,)/ )AD,c,an En d'autres termes, la somme des inventaires corrigés par les coefficients de chaque maille 0.l°x0.1° est égale au total de la maille l°xl° obtenu à partir des flux finaux pondérés. Au final, la méthode de mesure selon l'invention permet de corriger les résultats par installation en fonction des observations scientifiques et d'obtenir les émissions en Kg/m2/semaine sur des mailles 0.l°x0.1° (100km2) avec une précision supérieure à 5% et une réduction des sources d'incertitudes liées aux biais venant de la consommation en énergie, des statistiques de production d'énergie, des facteurs d'émission, des ratios de consommation énergétique et des omissions de sources adjacentes. En addition, lorsque plusieurs installations sont présentes sur une maille 0.l°x0.1°, la méthode calcule les émissions de processus ou de combustion de ces installations en fonction de leurs activités respectives (énergie, procédés industriels...) à partir des données publiées (rapports d'activité, rapports annuels...) tel qu'effectué actuellement dans les méthodes ascendantes. Le calcul est effectué de sorte que le total des émissions de ces installations corresponde au montant total corrigé de la maille 0.l°x0.1° pour la période considérée, ce qui permet de déduire le montant pour chacune des installations. 72 2953021 II. Autres Gaz à Effet de Serre A l'instar des autres méthodes de mesure, dont le Carbontracker, La méthode complète les mesures de CO2 par celles de CH4, N2O, NOx, HCFC, HFC, CFC, PFC, SF6, 03 et H2O. 5 Le module d'échange des autres GES inclue une modélisation spécifique des sources et des puits pour chacun d'entre eux avec le calcul du flux suivant le bilan massique. Ces modules d'échanges sont différents du CO2 et sont fonction des sources et des puits modélisés mathématiquement pour chaque 10 GES. Les mesures de ces GES sont issues des observations satellites, aériennes, atmosphériques, d'écosystèmes et marines telles que présentées dans la méthode. La méthode de mesure selon l'invention s'applique sur ces GES aussi bien que sur le CO2 pour l'obtention des inventaires corrigés qui 15 sont calculés sur une base hebdomadaire en Kg/m2 sur des mailles 0.l°x0.1° extrapolées annuellement. Les modélisations de l'évolution des flux effectuées par les modules d'échange pour chacun des autres GES considérés sont présentées ci-dessous. 20 • Méthane CH4 Le méthane (CH4) est essentiellement produit par des processus anaérobiques de sources naturelles incluant les zones humides, les forêts, les termites, les océans et de sources anthropiques par la production et la combustion 25 d'énergies fossiles, la culture du riz, l'élevage, les décharges, la combustion de biomasse, le traitement des déchets et le fumier. La méthode prend en compte une source additionnelle naturelle considérable de CH4 car d'énormes volumes de CH4 sont stockés sous les océans et dans des couches profondes du pergélisol sous forme d'hydrates où le gaz est enfermé dans des cages 73 cristallines de glace, stables à des pressions élevées et de basses températures. Ces hydrates représentent une ressource énergétique potentielle importante et sont généralement situés dans une couche de roche souterraine ou de sédiments océaniques appelée zone de stabilité d'hydrates (HSZ, Hydrate Stability Zone). Sous la HSZ, le CH4 est trouvé en phase gazeuse mélangé avec de l'eau et des sédiments. Lorsque les températures atmosphériques augmentent, notamment par le réchauffement climatique, la HSZ se déplace vers le haut, laissant à sa place une couche de gaz libérée par déstabilisation des hydrates. La pression dans cette nouvelle couche s'accroît, forçant le gaz à passer à travers la HSZ vers la surface par le biais de veines et de fractures. Si le CH4 s'avère être libéré en quantités massives par ce biais, ce dernier accélèrera le réchauffement en piégeant le rayonnement thermique environ 25 fois plus efficacement que le CO2 (Fig. 3) et deviendra très certainement, bientôt, la préoccupation majeure à laquelle la société devra faire face. La méthode utilise le modèle essentiel de Jain et al. (2009) qui modélise, à l'échelle du grain, comment le CH4 souterrain du fond des océans s'échappe par des évents dans le sol marin à un rythme beaucoup plus rapide qu'attendu. Il utilise un modèle discret d'élement (DEM, Discrete Element Model) qui permet d'investiguer la migration ascendante du CH4 dans sa phase de gaz libre et identifie que les facteurs principaux contrôlant le mode de transport du gaz dans les sédiments sont la taille des grains et le stress effectif de confinement. Combiné avec des données sismiques et des échantillons, le modèle fournit une explication physique de la récente découverte par le NORA d'un panache de 1'400 mètres venant du plancher océanique de CH4 et d'hydrates au large de la marge continentale de la Californie du Nord (Gardner, 2009). Les conditions sédimentaires dans lesquelles 74 le mécanisme de migration du gaz CH4 domine, sont perméables dans la majeure partie de l'océan, ainsi que dans certaines régions de pergélisol et ce modèle est utilisé dans la méthode pour reproduire le dégagement de CH4 du fond des océans et du pergélisol. Le puits atmosphérique principal de CH4 est sa destruction troposphérique par les radicaux hydroxyles (OH) et en second lieu, il est principalement supprimé par absorption dans le sol avec oxydation par les bactéries et le transport vers la stratosphère où il réagit avec OH, Cl et O(1D). Son espérance de vie TCH4, est calculée par sa réaction avec les radicaux hydroxyles OH avec zoH = 1/(k14[OH]) ou kl est le coefficient de la réaction, [OH] la concentration et Taddltlonel défini son espérance de vie en fonction des autres puits additionnels mineurs selon: 1/TCH4 = 1/TOH + 1/Zaddtttonel Le flux est exprimé selon l'équation suivante FCH4(t) = FA(t) + Fblo(t) + Foce(t) + Fper(t) + Ffeu(t) ù CCH4(t)•L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère est le flux net des hydrates de CH4 des océans est le flux net des hydrates de CH4 du pergélisol est le flux net des sources liées aux incendies est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction 25 de l'espérance de vie 1/icH4• • Protoxyde d'azote (N20) Le protoxyde d'azote (N20) est essentiellement produit par des sources anthropiques (engrais azotés, procédés 75 industriels, transports, combustion de biomasse, combustion d'énergie fossile, aliments pour le bétail) et des mécanismes biologiques naturels dans les océans et les sols. Le principal puits de N2O est sa destruction par des réactions photochimiques dans la stratosphère impliquant la production d'oxydes d'azote et en second plan, la dénitrification par des bactéries du sol. FN20(t) = FA(t) + Fbto(t) + Foce(t) +Ffeu(t) - CN20(t)•L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère est le flux net des échanges atmosphère-océan est le flux net des sources liées aux incendies est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction 15 de l'espérance de vie 1/TN2O. • Oxydes d'azote (NOx=NO+NO2) Les oxydes d'azote(NOx=NO+NO2) sont essentiellement produits par la combustion d'énergie fossile, la combustion de biomasse, les émissions des sols, la foudre, l'oxydation de 20 l'ammoniaque et le trafic aérien. Le puits principal de NOx est son oxydation dans l'atmosphère et d'importants montants venant des sols sont utilisés par la canopée avant de s'échapper dans la troposphère. Les NOx sont aussi absorbés par déposition sèche sur les sols, dépôts pouvant ensuite 25 conduire à des émissions de N2O. Ils agissent comme GES indirects en produisant de l'O3 troposphérique par réactions photochimiques dans l'atmosphère. Ils ont aussi un effet sur l'abondance de radicaux OH car leur destruction donne lieu à une augmentation des OH, réduisant la durée de vie de 30 certains GES tels que le CH4. 76 FNOx(t) = FA(t) + Fbto(t) + Ffeu(t) - CNOx(t)• L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère est le flux net des sources liées aux incendies est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction de l'espérance de vie 1/TNOx. • Chlorofluorocarbures CFC Les chlorofluorocarbures CFC sont produits de manière anthropique (propulseurs d'aérosols, réfrigérants, nettoyants, climatiseurs, systèmes anti-incendie, procédés de fabrication). Ces molécules monte lentement dans la stratosphère et se dirigent vers les pôles ou ils sont décomposés par des processus photochimiques et ils détruisent l'ozone stratosphérique. FCFC (t) = FA(t) - CCFC (t) • L (t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction de l'espérance de vie 1/-t-CFC. • Les hydrofluorocarbures HFC et les hydrochlorofluorocarbures HCFC Les hydrofluorocarbures HFC, et hydrochlorofluorocarbures HCFC sont produits de manière anthropique (propulseurs d'aérosols, réfrigérants, nettoyants, climatiseurs, systèmes anti-incendie, procédés de fabrication, insolation, packaging) avec généralement une durée de vie de quelques 77 années et d'importants effets de serre. Ils réagissent avec OH dans la troposphère. FHFC (t) = FA(t) - CHFC (t)• L(t) FHCFC(t) = FA(t) ù CHCFC(t). L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques respectives 5 est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction de leur espérance de vie respective 1/THFC, 1/THCFC• • Les perfluorocarbures PFC Les perfluorocarbures PFC sont des GES presque entièrement 10 anthropiques (production d'aluminium, production trifluoroacétique ou TFA, fabrication de semi-conducteurs) et sont aussi issus de sources naturelles (fluorites). Un puits important est la destruction légère (photolyse) ou réactions ioniques dans la mésosphère. FPFC(t) = FA(t) + Fbio(t) ù CPFC(t)• L(t) 15 est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction 20 de l'espérance de vie 1/TPFC• • L'hexafluorure de soufre (SF6) L'hexafluorure de soufre (SF6) est un GES presque entièrement anthropiques (production de magnésium, fabrication de disjoncteurs et d'interrupteurs haute tension, de semi- 25 conducteurs, de solvants, utilisation dans les pneumatiques) et sont aussi issus de sources naturelles (fluorites). Le seul puits connu est la destruction légère (photolyse) ou 78 réactions ioniques dans la mésosphère. Le SF6 est un GES puissant et en raison de sa forte densité par rapport à l'air, il reste sur le bas de l'atmosphère, ceci limitant sa capacité de réchauffement. Un puits important est la destruction légère (photolyse) ou réactions ioniques dans la mésosphère. FSF6(t) ù FA(t) +Fbio(t) ù CSF6(t)•L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère 10 est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction de l'espérance de vie 1/CSF6• • L'ozone troposphérique (03) L'ozone troposphérique (03) est essentiellement issue de la 15 stratosphère et est également produite dans la troposphère par réactions photochimiques ou ses concentrations augmentent en relation avec des niveaux élevés de polluants atmosphériques provenant de sources anthropiques (combustion de biomasse, industrie, transport). Les puits photochimiques 20 dominants de l'ozone troposphérique sont le cycle catalytique de destruction incluant la réaction HO2+03 et la destruction photolytique impliquant la réaction de O(1D), un produit de la photodissociation de l'ozone. Un autre puits important est l'absorption par les plantes. Il agit aussi comme un GES 25 indirect car sa décomposition par la lumière du soleil produit des radicaux OH. F03(t) = FA(t) + Fb,o (t) ù CO3 (t) • L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère 79 est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction de l'espérance de vie 1/T03. • La vapeur d'eau (H20) La vapeur d'eau (H20) est fonction de la température, influencée par le climat. Dans la stratosphère, elle est essentiellement issue de l'oxydation du CH4, de l'augmentation du trafic aérien, des résidus de vapeur d'eau troposphérique et dans la troposphère, elle provient essentiellement de l'évaporation et de la transpiration de la végétation et des océans et est perdue par condensation et précipitation. Elle est également issue de sources anthropiques de l'industrie, des habitations et des transports. La vapeur d'eau, surtout stratosphérique, agit comme un GES puissant car une concentration plus élevée de vapeur d'eau absorbe davantage d'énergie IR thermique rayonnée par la terre et réchauffe l'atmosphère. La vapeur d'eau troposphérique est exprimée par : FH20(t) = FA(t) + Fbio(t) + Foce(t) - CH20(t). L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère est le flux net des échanges atmosphère-océan est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction 25 de l'espérance de vie 1/tin°. La vapeur d'eau stratosphérique est exprimée par : FH20(t) = FA(t) + FT(t) ù CH20(t). L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques 80 est le flux net issu des résidus de vapeur d'eau troposphérique est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction 5 de l'espérance de vie 1/ix2o• • Le monoxyde de carbone (CO) 10 Le monoxyde de CH4 et transports, combustion végétation de carbone (CO) provient de l'oxydation chimique d'autres hydrocarbures dans l'atmosphère, des de la combustion d'énergie fossile, de la de biomasse et de sources naturelles de la et des océans. Les puits de CO sont essentiellement sa réaction avec OH, ainsi que sa déposition sur le sol. Il a des effets de GES indirect importants en réagissant avec les radicaux OH dans l'atmosphère et conduit 15 également à la formation d'ozone troposphérique. Fco(t) = FA(t) + Fbio(t) + Foce(t) + Ffeu(t) - Cco(t)• L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère est le flux net des échanges atmosphère-océan 20 est le flux net des sources liées aux incendies est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction de l'espérance de vie 1/ico• • Le dihydrogène (H2) 25 Le dihydrogène (H2) est produit par l'oxydation de CH4 et principalement par la combustion d'énergie fossile. L'évolution future de génération d'électricité par hydrogène conduira potentiellement à une forte augmentation de ses émissions. Les puits de H2 sont essentiellement sa 81 suppression par la réaction avec OH et l'absorption par les microorganismes du sol. C'est aussi un GES indirect en réagissant avec les radicaux (OH). FH2(t) = FA(t) + Fb~o(t) + Ffeu(t) - CH2(t)•L(t) est le flux accumulé atmosphérique net est le flux net des sources anthropiques est le flux net des échanges atmosphère-biosphère est le flux net des sources liées aux incendies est la concentration atmosphérique est le taux moyen de perte dans l'atmosphère fonction de l'espérance de vie 1/TH2. III. Système de mesure Selon l'invention, la méthode de mesure est mise en oeuvre au moyen d'un système de traitement de données (Fig. 18, bloc 800) comprenant des moyens de mesure des concentrations et des flux des gaz à effet de serre tels que décrits ci-dessus (satellites, avions, stations de mesures atmosphériques, stations de mesures marines, navires et/ou stations de mesures d'écosystèmes, capteurs, capteurs d'écosystèmes, capteurs marins), au moins une base de données centralisée comprenant le module d'observation, des moyens d'extraction, comprenant des moyens de transfert de données automatisés, et assurant également l'interface nécessaire avec des réseaux de communication. Le système de mesure selon l'invention comprend également des moyens de calcul tels que une pluralité de serveurs d'information dédiés, des ordinateurs, mainframes, etc. Le système de mesure comprend en outre des moyens de reporting, une ou plusieurs interfaces graphiques et une ou plusieurs interface de pilotage d'installations. Comme il a été dit plus haut, chacun des modules de la méthode peut avantageusement être implémenté sous forme de software, hardware ou une combinaison des deux. De plus, 82 compte tenu de la relative complexité de la méthode de mesure selon l'invention, il est clair que le système de mesure qui l'implémente nécessite une forte puissance de calcul, des capacités de stockage de données importantes ainsi que des moyens de communication fiables et rapides. Comme il a été dit plus haut, l'invention a donc pour but de fournir des mesures affinées des émissions de GES par aire géographique donnée, et ce en effectuant la méthode de mesure selon l'invention. Ces mesures affinées peuvent ensuite, soit constituer le résultat final destiné à être pris en considération par des utilisateurs individuels ou institutionnels, soit être utilisés directement pour le pilotage technique d'installations industrielles. Dans le premier cas, une plateforme internet centralisée permet alors de visualiser et d'analyser les émissions de gaz à effet de serre d'une pluralité d'aires géographiques données couvrant la totalité du globe. Les mesures effectuées et les résultats sont transmis en continu, préférablement en temps réel, vers cette plateforme Internet. Les utilisateurs du système peuvent avantageusement mettre en place plusieurs axes d'analyse y compris, mais ne se limitant pas aux types de GES, coordonnées (latitude, longitude), valeurs de flux, temps, incertitude ainsi que les champs de résultats du module d'observation, du module d'échange, du module d'inventaires ascendants, du module de transport, du module d'inversion et d'assimilation, du module de pondération et du module de géocodage pour effectuer des analyses détaillées. La plateforme est accessible par Internet aux utilisateurs munis d'un ordinateur personnel ou d'un équipement similaire connecté, et ce, de préférence avec un accès sécurisé via une interface graphique. Cette interface permet aux utilisateurs de naviguer sur la carte à travers ces mailles en les 83 redimensionnant tel « GoogleEarth » et de visualiser l'évolution des flux en temps réel. Les droits d'accès aux données sont attribués en fonction des profils des utilisateurs et peuvent être limités géographiquement afin de préserver la confidentialité des inventaires (Bloc 802). Un reporting peut être effectué en fonction de la zone géographique souhaitée (monde, continents, états, pays, régions et installations), de la période de temps souhaitée (année, mois, semaine) et des types de GES. L'utilisateur sélectionne ensuite les sources ou groupes de sources anthropiques désirés qui sont géocodées sur la carte et le système agrège la somme des flux sur la zone et la période de temps considérée. Des rapports d'inventaires de GES, destiné aux exploitants des installations, peuvent être générés à tout instant. Ils incluent préférablement les inventaires des différents types de GES, le détail des mesures réalisées avec le type d'observation, la précision, la résolution et la continuité ainsi que des statistiques actualisées sur les niveaux historiques, actuels et les tendances. Les résultats sont en TCO2/an, puis TCO2eq/an après avoir appliqué la méthode sur les autres GES et obtenu les TGES/an (Bloc 803). Le système fournit des observations in situ avec une cartographie quasi temps-réel des sources et des puits de GES aux échelles mondiales, continentales, étatiques, nationales, locales, jusqu'au niveau des installations pour refléter la réalité des niveaux d'émissions. Une utilisation prévue est de fournir cet accès aux données aux exploitants des installations, qui souhaitent de leur plein gré, ou si elles sont réglementées, mesurer et gérer leurs inventaires de GES. La précision, la continuité et l'uniformité des mesures, leur permet de compléter les processus actuels de surveillance, de reporting et de 84 vérification (MRV) effectués par les vérificateurs privés notamment ceux accréditées par la Commission Européenne sur le marché EU-ETS des émissions. Les exploitants peuvent au moyen du système, accéder, visualiser et obtenir des rapports détaillés sur les sources et les puits locaux de GES relatifs à leurs installations afin de vérifier continuellement l'évolution des niveaux et vérifier l'efficacité des technologies d'atténuation mises en place. Le système permet également à celles réglementées sur les marchés des émissions, de planifier leur budget de GES en fonction de leurs inventaires actuels, de l'évolution des prix des permis de GES échangés sur les marchés et d'évaluer les crédits d'émissions dont elles auront besoin chaque année pour rester en conformité avec les autorités. Dans l'autre cas, le système de mesure selon l'invention peut également être interfacé directement au sein d'une installation émettrice, notamment avec un système de gestion de production, pour permettre le pilotage de l'installation afin de limiter les émissions de combustion et/ou de processus et d'automatiser leur réduction. Un logiciel spécifique est installé par installation en fonction de son activité (énergie, procédés industriel, usages de produit...), de ses processus et des GES émis. Le système de mesure selon l'invention permet alors de calibrer et d'optimiser directement le processus de chaque installation en fonction des niveaux et types d'émissions mesurés (ex : pics de pollution). Ceci permet d'obtenir une réduction d'émission automatisée sur chaque installation, de contrôler au fur et à mesure du temps son efficacité et de rester en conformité avec les normes régulatoires et environnementales. 85 A titre d'exemple, dans l'énergie, les installations sont à la recherche de moyens permettant de préserver la qualité de l'air en opérant des unités plus productives. Dans une centrale à charbon, plus la combustion et la production d'électricité augmente, plus les émissions de CO2, de NOx et de CO sont dégagées. En interfaçant alors le système de mesure selon l'invention avec des logiciels permettant d'automatiser les processus de production d'électricité, on peut alors optimiser le rapport entre réduction d'émissions et efficacité des brûleurs. Des logiciels transmettent des paramètres d'exploitation de différents points de la centrale dans un ordinateur de contrôle principal et d'autres, interfacés au système de mesure selon l'invention, optimisent la combustion en ajustant les flux d'air et de combustible dans les brûleurs, stabilisant ainsi les niveaux de GES. Interfacé avec le système de gestion de production d'une centrale à charbon, le système identifie et quantifie précisément les types de gaz émis (CO2, CH4, N2O, NOx, HFC, HCFC, CFC, PFC, SF6, 03, H20, CO, H2) puis contrôle et optimise au moyen de logiciels ad-hoc les processus de production d'électricité en influant, par exemple, sur le rapport entre réduction d'émissions et efficacité des brûleurs. L'homme du métier sait adapter les interfaces nécessaires 25 entre le système de mesure selon la présente invention et le/les systèmes de pilotage d'installations de production. 86 Revendications1. Méthode de mesure des émissions hebdomadaires et annuelles d'un gaz à effet de serre générées sur une aire géographique déterminée, caractérisée en ce qu'elle comprend les étapes suivantes : - effectuer des mesures journalières de concentration dudit gaz à effet de serre en une première pluralité de lieux répartis sur la totalité du globe terrestre et enregistrer lesdites mesures journalières de concentration dans un module d'observation, - effectuer des mesures journalières des flux dudit gaz à effet de serre en une deuxième pluralité de lieux répartis sur la totalité du globe, et enregistrer lesdites mesures journalières des flux dans ledit module d'observation, - effectuer des mesures de paramètres satellites, de paramètres météorologiques, de paramètres marins et de paramètres d'écosystèmes en une troisième pluralité de lieux répartis sur le globe terrestre et 20 enregistrer les dites mesures de paramètres dans ledit module d'observation, - extraire, au moyen d'un module d'extraction, des données de prévisions météorologiques d'au moins une source de données, 25 - effectuer une modélisation de l'évolution des flux dudit gaz sur le globe au moyen d'un module d'échange modélisant les sources et les puits naturels et anthropiques, - effectuer une modélisation des émissions 30 anthropiques hebdomadaires dudit gaz à effet de serre au moyen d'un module d'inventaires ascendants, ledit module intégrant des données brutes d'émissions pour une pluralité d'installations, 15 87- effectuer, en utilisant ladite modélisation de l'évolution des flux, ladite modélisation des émissions anthropiques hebdomadaires et les dites données de prévisions météorologiques, une modélisation du transport atmosphérique dudit gaz à effet de serre au moyen d'un module de transport, - calculer les flux finaux dudit gaz à effet de serre, au moyen d'un module d'inversion et d'assimilation de données utilisant ladite modélisation des flux effectuée par le module d'échange, ladite modélisation des émissions anthropiques hebdomadaires effectuée par le module d'inventaires ascendants, ladite modélisation du transport atmosphérique effectuée par le module de transport et lesdites mesures enregistrées dans ledit module d'observation, - pondérer, au moyen d'un module de pondération, les dits flux finaux de façon à fournir des flux finaux pondérés, - calculer, en utilisant les dits flux finaux pondérés et ladite modélisation des émissions anthropiques hebdomadaires effectuée par le module d'inventaires ascendants, les émissions hebdomadaires dudit gaz à effet de serre de ladite aire géographique, au moyen d'un module de géocodage comprenant au moins un système d'information géographique, - extrapoler, à partir des dites émissions hebdomadaires, les émissions annuelles dudit gaz à effet de serre de ladite aire géographique. 2. Méthode de mesure selon la revendication 1, caractérisée en ce que la surface de ladite aire géographique est comprise entre lkm2 et 10000km2, en particulier que 88ladite aire géographique inclue au moins une source anthropique donnée. 3. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que le dit gaz à effet de serre est choisi parmi le dioxyde de carbone (CO2), le méthane (CH4), le protoxyde d'azote (N2O), les oxydes d'azote (NOx), les hydrofluorocarbures (HFC), les hydrochlorofluorocarbures (HCFC), les chlorofluorocarbures (CFC), les perfluorocarbures (PFC), l'hexafluorure de soufre (SF6), l'ozone (03), la vapeur d'eau (H2O), le monoxyde de carbone (CO) et le dihydrogène (H2). 4. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que les dites mesures journalières de concentration dudit gaz à effet de serre sur le globe, les dites mesures journalières des flux dudit gaz à effet de serre sur le globe, les dites mesures de paramètres satellites, de paramètres météorologiques, de paramètres marins et de paramètres d'écosystèmes sont effectuées au moyen d'une pluralité de satellites, d'avions, de stations de mesures atmosphériques, de stations de mesures marines, de navires et/ou de stations de mesures d'écosystèmes permettant d'effectuer des mesures sur la totalité du globe. 5. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module d'échange effectue ladite modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module solaire modélisant le rayonnement solaire en utilisant les paramètres orbitaux de la 89géométrie terrestre avec un calcul de l'excentricité de la terre déterminé proportionnellement à l'excentricité de Mars. 6. Méthode de mesure selon l'une des précédentes, caractérisée en ce que d'échange effectue ladite modélisation de absorptivité et transmissivité de l'absorption par les gaz à effet de serre revendications ledit module l'évolution des de l'holocène, modélisant le réflectivité, l'atmosphère, et les nuages, flux dudit gaz à effet de serre, à partir en utilisant un module rayonnement ondes-courtes, énergétique en incluant les variations d'albédo planétaire et l'influence du trou de la couche d'ozone, ledit module énergétique modélisant également le rayonnement ondes-longues, en utilisant l'équation de Schwartzschild, la méthode des émissivités et en incluant l'absorption et l'émission par les gaz à effet de serre et les nuages de rayonnements ondes-longues, les flux de chaleur latente, les flux de chaleur sensible, les flux de conduction et la température de surface. 7. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module 25 d'échange effectue ladite modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module océan modélisant l'effet net des échanges atmosphère-océan sur la base du modèle MOM3 combiné avec lesdites données de prévisions 30 météorologiques et prenant en compte l'effet de tampon, l'absorption par altération chimique suivant le modèle CDIAC DBl0l2 et le dégagement par évaporation. 20 90. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module d'échange effectue la modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module biosphère modélisant l'effet net des échanges atmosphère-biosphère sur la base du modèle JSBACH et incluant les types de plantes de la biosphère, l'indice foliaire, la lumière, l'albédo, la photosynthèse C3 et C4, l'addition du taux de photosynthèse brute limité, la respiration autotrophique, la respiration hétérotrophique et/ou la modification anthropique de la couverture de végétation depuis au moins le dernier millénaire. 9. Méthode de mesure selon la revendication précédente, caractérisée en ce que ledit module biosphère utilise un module incendie modélisant les perturbations dues aux incendies sur la base de données extraites de la base de données Global Fire Emission Database (GFEDv2) intégrées au modèle JSBACH. 10. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module d'échange effectue ladite modélisation de l'évolution des flux dudit gaz à effet de serre, à partir de l'holocène, en utilisant un module fossile modélisant les émissions anthropiques fossiles à l'échelle globale sur la base des statistiques de production de pétrole et de charbon de l'Energy Information Administration (EIA) et les estimations de Etemad et al. 11. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module 91d'inventaires ascendants extrait des inventaires d'émissions de la base de données EDGAR 4.0 en incluant un calcul de la variabilité temporelle des émissions. 12. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module de transport atmosphérique utilise le modèle de transport TM5 combiné avec les dites données de prévisions météorologiques pour calculer le transport des flux dudit gaz à effet de serre sur le globe. 13. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module d'inversion et d'assimilation de données utilise, pour calculer les dits flux finaux, une inversion de synthèse avec la fonction de Green pour les grandes régions et le filtre d'ensemble de Kalman. 14. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module de pondération utilise, pour pondérer les dits flux finaux, une analyse des activités de production des pays et des régions du monde accompagnée d'une modélisation des marchés des émissions basée sur le modèle de biens publics produits par le privé. 15. Méthode de mesure selon l'une des revendications précédentes, caractérisée en ce que ledit module de géocodage utilise des coefficients correcteurs. 16. Système de mesure effectuant la méthode selon l'une des revendications 1 à 15 comprenant : des moyens de mesure (801) des concentrations et des flux des gaz à effet de serre, 92- des moyens de mesures (801) de paramètres satellites, météorologiques, marins et d'écosystèmes, - au moins une base de données centralisée (803) comprenant un module d'observation, - des moyens d'extraction (802) et de transfert de données automatisés, - des moyens de calcul (805) comprenant au moins un module d'échange, au moins un module d'inventaires ascendants, au moins un module de transport, au moins un module d'inversion et d'assimilation de données, au moins un module de pondération, - au moins un module de géocodage (804) comprenant un système d'information géographique permettant de géocoder les résultats fournis par les dits moyens de calcul, - une plateforme internet centralisée permettant de visualiser et d'analyser les émissions de gaz à effet de serre d'une pluralité d'aires géographiques données. 17. Système de mesure selon la revendication précédente, caractérisé en ce qu'il comprend des moyens d'interface (808) avec un système de gestion de production d'une installation. 93