Der regionalen Schwerefeldmodellierung in polaren Gebieten kommt in vielerlei Hinsicht eine besondere und wachsende Bedeutung zu. Einerseits sind hochauflösende und präzise Geoidmodelle eine wichtige Eingangsgröße bei der Untersuchung und Quantifizierung geophysikalischer, ozeanographischer bzw. glaziologischer Phänomene, z. B. bei der Bestimmung der mittleren dynamischen Ozeantopographie oder der Anwendung des Schwimmgleichgewichts im Bereich von Schelfeisen, Meereis oder subglazialen Seen. Zudem werden sie allgemein zur Referenzierung von Höhenmodellen benötigt. Andererseits sind, aufgrund der unvermeidbaren polaren Datenlücken von Satellitenbeobachtungen jenseits der Grenzbreite (sog. „Polloch“), terrestrische Schweredaten auch für die globale Schwerefeldmodellierung unerlässlich. Jedoch sind die verfügbaren terrestrischen (bodennahen) Schwerebeobachtungen insbesondere im Gebiet der Antarktis äußerst lückenhaft und heterogen. So entspricht das tatsächliche Auflösungsvermögen selbst aktueller kombinierter Schwerefeldmodelle wie EGM2008 oder EIGEN-6C über dem antarktischen Kontinent lediglich dem der reinen Satellitenmodelle aus GRACE bzw. GOCE. Des Weiteren sind Standardverfahren der regionalen Geoidmodellierung hier nicht ohne Weiteres anwendbar. Neben der Heterogenität der Daten als praktischer Herausforderung muss aus theoretischer Sicht dem zusätzlichen Dichtekontrast durch das Eis Rechnung getragen werden. Die vorliegende kumulative Dissertation greift diese Problematik auf. Während die einzelnen Publikationen die Ergebnisse ausgewählter regionaler Fallstudien präsentieren, soll die folgende zusammenfassende Abhandlung einen doppelten Bogen spannen, indem die geophysikalischen Phänomene gleichzeitig als zu untersuchende Anwendungsgebiete und als Einflussfaktoren im Kontext der regionalen Geoidmodellierung beschrieben werden. So wird am Beispiel der Weddellsee gezeigt, wie die Meereisbedeckung die Qualität und Zuverlässigkeit der mithilfe der Satellitenaltimetrie abgeleiteten Schwerefeldmodelle beeinträchtigt. Diese Modelle bilden derzeit die alleinige Datengrundlage für die hochauflösenden globalen Modelle im Gebiet des Antarktischen Ozeans. Zugleich wird anhand des verfeinerten regionalen Modells und daraus abgeleiteter geostrophischer Geschwindigkeiten demonstriert, dass selbst lückenhafte und heterogene terrestrische Daten hier einen wesentlichen Beitrag zur simultanen Kalibrierung und Vereinheitlichung des Datenbestandes leisten können. Im Ergebnis konnten in den küstennahen Gewässern Differenzen von mehreren Dezimetern gegenüber Geoidhöhen aus EGM2008 festgestellt werden, welche teils auf systematische Abweichungen und teils auf Rauschen im globalen Modell zurückzuführen sind. Über dem Festland erreicht dessen Vernachlässigungsfehler im quadratischen Mittel sogar 0,75 m und Maxima von über 3 m. Ein weiteres verfeinertes und, dank geeigneter Eingangsdaten, sehr genaues und hochauflösendes Geoidmodell wird für die Region um den Vostoksee in der Ostantarktis abgeleitet. In Kombination mit Eisoberflächenhöhen und Eisdicken gelingt es, das Schwimmgleichgewicht des subglazialen Sees nachzuweisen. Das gegenüber GOCE zusätzlich gewonnene Geoidsignal ist hier mit 0,56 m Standardabweichung zwar etwas kleiner, jedoch wird im Vergleich mit der residualen Auslenkung des Seespiegels (0,26 m Standardabweichung) auch für diese Anwendung der signifikante und gegenüber dem Auflösungsvermögen von GOCE auch notwendige Beitrag eines regionalen Geoidmodells deutlich. Für das hydrostatische Gleichgewicht eines subglazialen Sees ist streng genommen das tatsächliche Schwerepotential in Höhe des Seespiegels maßgeblich. Dessen Berechnung erfordert eine Fortsetzung des Störpotentials nach unten innerhalb der Topographie, welche konzeptionell in engem Zusammenhang mit dem bekannten Geoid-Quasigeoid-Separationsterm steht. Dessen oft angenommene Approximation mithilfe der Bougueranomalie kann, angesichts der heutigen Anforderungen an ein modernes zentimetergenaues Geoid, gerade in rauem Gelände zu ungenau sein. In Anlehnung an aktuelle Arbeiten auf diesem Gebiet wird ein verallgemeinerter und zugleich verfeinerter Ansatz zur praktischen Berechnung des Terms erarbeitet. Am Beispiel des Himalaya werden die einzelnen Anteile im Rahmen einer Simulationsstudie quantifiziert und insbesondere ihre Sensitivität gegenüber dem Integrationsradius der Topographie untersucht. Besonderes Augenmerk liegt ebenso auf dem indirekten Effekt der Topographie in Bezug auf das Potential, welcher, im Gegensatz zur Anwendung eines planaren Modells, in sphärischer Betrachtungsweise nicht verschwindet. / In many respects, regional gravity field modeling in polar areas is of special, and growing, interest. On the one hand, high-resolution and precise geoid models are an important input parameter to investigate and quantify manifold geophysical, oceanographical and glaciological phenomena, e.g., the determination of the mean dynamic ocean topography, or the application of the hydrostatic equilibrium condition in the areas of ice shelves, sea ice, or subglacial lakes. Moreover, geoid models are in general needed as a reference for height models. On the other hand, because of the unavoidable polar data gaps in satellite measurements due to the inclination (the so-called “polar gap”), terrestrial gravity data are indispensable also for global gravity field modeling. However, the available terrestrial (ground-based) gravity datasets, in particular of Antarctica, are very sparse and heterogeneous. For example, over the Antarctic continent the true resolution of even the most recent combined global geopotential models such as EGM2008 or EIGEN-6C only corresponds to that of the satellite-only models derived from GRACE and GOCE, respectively. Furthermore, standard techniques of regional geoid modeling cannot be readily used in this area. Apart from the heterogeneity of the data as a practical challenge the additional density contrast implied by the covering sheet needs to be accounted for from the theoretical point of view. This complex situation is the starting point for the present cumulative dissertation. Whereas the individual publications present the results of selected regional case studies, the intention of the following summary is to draw an integrated picture aiming at explaining the geophysical phenomena as both applications and influencing factors in the context of regional geoid modeling. Using the example of the Weddell Sea it is shown how sea-ice coverage affects the quality and reliability of marine gravity field models derived from radar satellite altimetry. At present, these models are the only input data to the high-resolution global geopotential models. At the same time, the refined regional model and geostrophic velocities derived thereof are employed to demonstrate how even sparse and heterogeneous terrestrial gravity data may provide a contribution to simultaneously calibrate and unify the available datasets. As a result, near the coast differences at the order of some decimeters could be observed in comparison with EGM2008, originating partly from systematic effects and noise in the global model. In the continental areas, its omission error even yields a standard deviation of 0.75 m and attains a maximum of more than 3 m. Another refined and, owing to appropriate input data, very precise and highly resolving geoid model is derived for the region around subglacial Lake Vostok. In combination with ice-surface heights and ice thickness data it is used to provide observational evidence that the lake is in a state of hydrostatic equilibrium. There, the additional geoid signal w.r.t. GOCE is a bit smaller (0.56 m standard deviation). However, considering the residual deviations of the apparent lake level (0.26 m standard deviation) the significant and necessary, as compared to the resolution of GOCE, contribution of a regional geoid model to this application is shown. In a strict sense, the relevant quantity to evaluate the hydrostatic equilibrium condition of a subglacial lake is the actual geopotential at the anticipated lake level. Its computation requires a downward continuation of the disturbing potential inside the topography, which is closely related to the concept of the well-known geoid-quasigeoid separation term. In the past, this term was frequently described as an approximation by means of the Bouguer anomaly. However, considering the modern requirements of the “one-centimeter geoid” this approximation may be too coarse over rough terrain. Following recent works in this field, a generalized yet refined approach for practical implementation of the term is developed. The individual constituents of the term are quantified. In particular, their sensitivity against the radius up to which topography is taken into account is investigated. For this simulation study, the Himalaya mountain region served as test area. Furthermore, special focus is given to the indirect of topography on the potential which, contrary to applying a planar model, does not vanish in the spherical approach.

Les changements du niveau de la mer sont un des indicateurs des changements climatiques qui affectent le monde depuis les dernière décennies. Dans les zones côtières, ce phénomène est à l'origine de l'érosion côtière, l'intrusion d'eau salée dans la nappe phréatique et les inondations. Les côtes du Canada ne sont pas épargnées. En effet, le niveau de la mer dans ces régions augmente actuellement et suit les tendances mondiales. Parmi ces zones côtières, celles du golfe du Saint-Laurent nous ont intéressées. Elles sont caractérisées par des propriétés physiques spécifiques comme l'isolement relatif par rapport l'océan Atlantique, la couverture de glace en hiver et la coulée de l'eau douce au printemps qui leurs confèrent un comportement particulier en termes de changements du niveau de la mer. Les changements du niveau de la mer peuvent être déterminés en intégrant les données des marégraphes et de l'altimétrie satellitaire que le résultat de cette procédure conduit à une détermination relative du niveau moyen de la mer. Cependant, parce que la structure de mesure de ces données est différente, certaines corrections doivent être appliquées. L'une des corrections parmi les plus essentielles est la longueur du Mouvement de Terrain Vertical (MTV) qui doit être prise en compte dans les données des marégraphes, car ces derniers sont connectés directement à la croûte terrestre. Différents phénomènes peuvent causer le MTV. À une échelle locale, la tectonique, l'affaissement et la sédimentation sont les facteurs qui créent le mouvement vertical et à une méso ou macro-échelle, la fonte des inlandsis qui entraîne une modification de la charge de masse de la surface de la terre, provoquent le mouvement des terres. Il existe divers modèles et méthodes de correction du MTV globaux et locaux, tels que l'Ajustement Isostatique Glaciaire (AIG), le Système Mondial de Navigation par Satellite (SMNS), l'Interférométrie par Radar Synthèse d'Ouverture (InRSO), Gravity Recovery and Climate Experiment (GRACE) et le calcul de la différence entre données marégraphiques et altimétriques par satellite. Parmi ces méthodes, l'utilisation de la vitesse verticale des données GPS est la plus précise. Cependant, en raison du faible nombre de stations GPS dans ou proximité des stations marégraphiques, l'utilisation d'une approche plus globale est inévitable. En utilisant cette correction, le niveau moyen absolu de la mer sera obtenu. Dans cette étude, l'évaluation du MTV se fait par un algorithme d'InRSO, appelé Diffuseurs Persistants InRSO (DPInRSO en anglais PSInSAR). Les images du RSO du satellite Sentinel-1 sont soumises à l'analyse DPInRSO. Afin d'évaluer la précision de cette approche, les résultats de DPInRSO (obtenus du logiciel Stanford Method for Persistent Scatterers (StaMPS)) pour deux régions de l'est du Canada (Halifax et St. John's) sont comparés à la GPS vitesse verticale en analysant les séries temporelles de données. Également, en utilisant la même méthode, le mouvement vertical des terres dans 24 stations marégraphiques du golfe du Saint-Laurent sera corrigé. L'étude du niveau moyen de la mer est mise en oeuvre en utilisant l'analyse spectrale (analyse harmonique et analyse spectrale des moindres carrés (ASMC)) afin de supprimer les effets à long terme des constituants de la marée (constituants annuels et semi-annuels). Aussi, pour évaluer le niveau moyen absolu de la mer (NMAM) dans la région considérée, le niveau moyen de la mer calculé est corrigé sur la base des résultats du MTV. En conclusion, la comparaison des résultats de la technique DPInRSO avec les données GPS indique le comportement similaire des tendances au cours de la période d'étude. En d'autres termes, la méthode DPInRSO mesure le mouvement vertical des terres en précision millimétrique dans la région d'étude (à l'exception d'une station) et elle peut être appliquée pour les régions qu'elles ne se connectent pas sur le positionnement GPS continu des marégraphes. Dans le cas d'un changement du niveau de la mer, l'incompatibilité entre les données marégraphiques et altimétriques en termes de temps rend la conclusion plus diffcile. Cependant, sur la base d'observations continues dans des séries temporelles de données altimétriques, nous avons obtenu une tendance plus homogène du niveau de la mer (niveau moyen de la mer) dans toutes les stations, mais pour évaluer le changement à long terme du niveau de la mer et en raison de l'emplacement exact de la station marégraphique, ces dernières données doivent être prises en compte (à la fois dans le calcul du niveau moyen de la mer et des constituants de la marée). Enfin, en raison du faible changement de MTV, les valeurs moyennes absolues du niveau de la mer sont très proches de celles relatives, ce qui indique que nous pouvons évaluer les changements du niveau de la mer dans le golfe du Saint-Laurent, en utilisant uniquement les données marégraphiques et altimétriques. En plus de faire des corrections, des cartes de la tendance du niveau de la mer et des différences de phase entre les stations marégraphiques seront représentées. La qualité des données des marégraphes et de l'altimétrie satellitaire sera discutée et rapportée sur les résultats tirés de l'étude. / Sea level changes are one of the indicators of climate change that has affected the world in recent decades. In coastal areas, this phenomenon has caused coastal erosion, saltwater intrusion into the groundwater and floods. These environmental changes can also be seen in Canada's coasts. Indeed, the sea level in these regions is currently increasing and follows the global trends. Among these coastal areas, those of the Gulf of St. Lawrence interested us. This area is characterized by specific physical properties such as relative isolation from the Atlantic Ocean, winter ice cover and spring freshwater flow that confer particular behavior in terms of sea level changes. Sea level changes can be determined by integrating the tide gauge and satellite altimetry data that the result of this procedure leads to a relative determination of mean sea level. However, because the measurement structures of these data are different, some corrections should be applied. One of the most essential corrections is the length of the Vertical Land Motion (VLM) which must be taken into account in tide gauge data, as these are connected directly to the Earth crust. Different phenomena can cause VLM. In the local scale, tectonics, subsidence and sedimentation are the factors that create vertical movement and in global scale, melting of ice sheets, which leads to change the mass loading of the earth's surface, cause the land motion. There are various global and local VLM correction models and methods such as, Glacial Isostatic Adjustment (GIA), Global Navigation Satellite Systems (GNSS), Interferometry Synthetic Aperture Radar (InSAR), Gravity Recovery and Climate Experiment (GRACE) and calculation of difference between satellite altimetry and tide gauge data. Between these methods, using the vertical velocity of GPS data is the most accurate one. However, due to the low number of GPS stations in or near the tide gauge stations, using a more comprehensive approach is inevitable. By applying this correction, the absolute mean sea level will be obtained. In this study, VLM is examined by one algorithm of InSAR which is called Persistent Scatterer InSAR (PSInSAR). The Sentinel-1 single-look SAR images are used for PSInSAR analysis. In order to accuracy assessment of this approach, the results of PSInSAR (obtained from Stanford Method for Persistent Scatterers software (StaMPS)) for two regions in east of Canada (Halifax and St. John's) are compared with GPS vertical velocity by analyzing their time series data. Equally, by using the same method, vertical land movement in 24 tide gauge station in Gulf of St. Lawrence will be corrected. Investigation of mean sea level is implemented by using the spectral analysis (Harmonic analysis and Least Square Spectral Analysis) in order to remove the long-term effects of tidal constituent (annual and semiannual constituents). Also, to assess the Absolute Mean Sea Level (AMSL) in considered region, calculated mean sea level is corrected based on the results of VLM. In conclusion, comparing of the results of PSInSAR technique with GPS data indicate the similar behavior of trends during the period of study. In other words, PSInSAR method measures vertical land motion in millimeter precision in study region (except of one station) and it can be applied for regions that they do not connect on continuous GPS positioning of tide gauges. In the case of sea level change, incompatibility between tide gauge and altimetry data in terms of time makes the conclusion harder. However, based on continues observations in time series of altimetry data, we obtained more homogeneous sea level trend (mean sea level) in all stations, but to evaluate the long-term sea level change and because of exact location of tide gauge station, this latter data must be considered (both in calculation of mean sea level and tidal constituents). Finally, because of small values of VLM, absolute mean sea level values are very close to relative one which indicate that, we can evaluate the sea level changes in Gulf of Saint Lawrence by using only the tide gauge and altimetry data. In addition to making corrections, the maps of sea level trend and phase differences between tide gauge stations will be represented. The quality of the tide gauge and satellite altimetry records will be discussed and reported over the drawn results of the study.

This thesis, structured in two parts, addresses a series of problems of relevance in the field of Spatial Geodesy. The first part delves into the application of satellite gravity data to enhance our understanding of water transport dynamics. Here, we present two significant contributions. Both are based on satellite gravity data but stem from different mission concepts with distinct objectives: time-variable gravity monitoring and high-resolution, accurate static geoid modelling. First, the fundamental notions about gravity are introduced and a brief summary is made of the different gravity satellite missions throughout history, with emphasis on the GRACE/GRACE-FO and GOCE missions, whose data are the basis of this work. The first application focuses on estimating water transport and geostrophic circulation in the Southern Ocean by leveraging a GOCE geoid and altimetry data. The Volume Transport across the Antartic Circumpolar Current is analyzed and the resulsts are validated validated using the in-situ data collected during the multiple campaigns in the DP. The second application uses time-variable gravity data from the GRACE and GRACE-FO missions to estimate the water cycle in the Mediterranean and Black Sea system, a critical region for regional climate and global ocean circulation. The analysis delves into the analysis of the different components of the hydrological cycle within this region, including the water flow across the Gibraltar Strait, examining their seasonal variations, climatic patterns, and their connection with the North Atlantic Oscillation Index. The second part of the thesis is more focused on data analysis, with the objective of developing mathematical methods to estimate the cross correlation function between two time series that are both unevenly spaced spaced (the sampling is not uniform over time) and observed at unequal time scales (the set of time points for the first series is not identical to the set of time points of the second series). Such time series are frequently encountered in geodetic surveys, especially when combining data from different sources. The estimation of the the cross correlation function for these time series presents unique challenges and requires the adaptation of traditional analysis methods designed for evenly spaced and synchronized time series. The two main contributions in this context are: (i) the study of the asymptotic properties of the Guassian Kernel estimator, that is the recommended estimator for the cross correlation function when the two time series are observed at unequal time scales; (ii) an extension of the stationary bootstrap that allows to construct bootstrap-based confidence intervals for the cross correlation function for unevenly spaced time series not sampled on identical time points.

Depuis le début des années 1990 on suit l'évolution globale du niveau de la mer grâce aux satellites altimétriques. Ils observent une hausse du niveau moyen global de la mer (GMSL) de 3.4 ± 0.4 mm/an sur la période 1993-2016 (ce qui représente le double de ce qui a été observé au cours du 20ème siècle par les marégraphes, hausse à 1.7 mm/an entre 1900 et 1990). Le GMSL présente aussi des fluctuations interannuelles qui peuvent atteindre quelques millimètres, surtout pendant les épisodes ENSO. Cette hausse n'est pas régionalement uniforme : elle a été 3 fois plus rapide que la hausse moyenne globale dans certaines zones entre 1993 et 2016. Au cours du 21ème siècle, on s'attend à une hausse accrue du GMSL pouvant aller jusqu'à 1 m à l'horizon 2100, avec une forte variabilité régionale. Il est donc important de comprendre l'évolution actuelle du niveau des océans qui constitue une menace sérieuse pour de nombreuses régions côtières basses souvent très peuplées. Cette thèse s'inscrit dans le contexte du projet niveau de la mer CCI (Climate Change Initiative) de l'Agence Spatiale Européenne (ESA) ayant pour objectif de fournir de meilleurs produits du niveau de la mer combinant les missions Topex/Poseidon, Jason-1/2, ERS-1/2 et Envisat. L'objectif premier de cette thèse est de valider ces produits SL_CCI du niveau de la mer en utilisant différentes approches, en particulier par l'étude du bilan (comparaison du GMSL observé avec la somme des différentes contributions : composante stérique, fonte des glaces continentales et transferts d'eau depuis les terres émergées). Un autre objectif est d'estimer les composantes du niveau de la mer mal connues, et tout particulièrement le contenu thermique de l'océan profond non mesurable par le système Argo, et la contribution du stock d'eau sur les continents. Ces travaux ont montré que la contribution de l'océan profond en dessous de 2000m est faible sur la période 2005-2013 et contenue dans la barre d'incertitudes des données (erreurs qui proviennent essentiellement, (1) des produits niveau de la mer altimétriques et des lacunes de la couverture géographique des données Argo dans la région Indonésienne pour la tendance et (2) des produits GRACE et Argo pour la variabilité interannuelle). Nos résultats et la méthode utilisée montrent que le niveau de la mer et ses composantes sont encore entachés d'erreurs importantes. Dans la deuxième partie, nous avons analysé l'influence du phénomène ENSO (El Niño et La Niña) sur les variations interannuelles du GMSL. Nous montrons que lors des évènements La Niña comme celui de 2010-2011, le déficit de précipitations sur l'océan (et l'excès sur les continents) conduit à une baisse temporaire de la masse de l'océan global et donc du niveau de la mer. C'est essentiellement la variation de masse de l'océan qui explique la variabilité interannuelle du niveau de la mer lors des évènements ENSO, et le déficit (La Niña) ou excès (El Niño) de masse se trouve confiné dans l'océan Pacifique tropical Nord. Pour finir, nous analysons l'évolution de la température moyenne de l'air et de l'océan en surface sur la période du "hiatus" (2003-2013). Nous montrons que ce hiatus, c'est à dire le ralentissement récent de la hausse de la température moyenne globale de la Terre est un phénomène quasi global, même si le Pacifique tropical Est s'est fortement refroidi. Cette "supposée" pause récente s'explique par la variabilité naturelle interne du climat. La Terre est toujours en état de déséquilibre énergétique dû à l'accumulation de gaz à effet de serre. Nous mettons en évidence le rôle de la variabilité naturelle à court terme sur les changements à plus long terme associés au réchauffement climatique anthropique. / Since the early 1990s sea level is routinely measured using high-precision altimeter satellites. These observe a rise in global mean sea level (GMSL) of 3.4 ± 0.4 mm/yr over the 1993-2016 period (which is twice what has been observed during the 20th century by the tide gauges, with a rise of 1.7 +/- 0.3 mm/yr). The interannual variability in the GMSL can reach several millimeters, especially during ENSO events. The rate of sea level rise is not regionally uniform. During the altimetry era, it was three times faster than the global mean in some areas. During the 21st century, we expect a greater rise of the GMSL than today, up to 1 m in 2100, with strong regional variability. It is therefore important to understand the current evolution of the sea level, since it represents a serious threat to many low coastal areas, often densely populated of the planet. My thesis research deals with the Sea Level CCI (Climate Change Initiative) project of the European Space Agency (ESA) which objective is to provide improved sea level products combining several altimetry missions, including Topex/Poseidon, Jason-1/2, ERS-1/2 and Envisat. The primary objective of my thesis was to validate the CCI sea level products using different approaches, in particular the sea level budget approach. It consists of comparing the observed GMSL with the sum of different contributions : the steric component, melting of continental ice and transfers of water between the land surface and oceans. Another objective was to estimate the poorly known components to sea level rise, in particular the heat content of the deep ocean not measurable by Argo, and the contribution of water storage on the land. My work has shown that the contribution of the deep ocean below 2000m to the rising sea level is small over the 2005-2013 periods and not significant compared to the data uncertainties. The main uncertainties come from: (1) -in terms of trend- the altimetry sea level products and gaps in the geographical coverage of Argo data in the Indonesian region, and (2) -in terms of interannual variability- the GRACE and Argo products. My results and the method used show that the sea level and its components are still affected by important errors. In the second part, I analyzed the influence of ENSO (El Niño and La Niña) on the interannual variations of the GMSL. I showed that during La Niña events, like that of 2010-2011, the rainfall deficit over the ocean (and excess over the continents) leads to a temporary decrease in the global ocean mass and therefore in the GMSL. This is essentially the ocean mass variation that explains the interannual variability of the GMSL during ENSO events. Furthermore, the deficit (La Niña) or excess (El Niño) ocean mass is confined in the north tropical Pacific Ocean. Finally, I analyzed the evolution of the average temperature of air and ocean surface over the period of the "hiatus" (2003-2013). I showed that this hiatus, i.e. the recent slowdown in the rise of the global mean Earth's temperature is an almost global phenomenon, though cooling of the tropical eastern Pacific has slightly contributed. This recent pause is attributable to natural internal climate variability. The Earth is indeed still in a state of energetic imbalance due to the accumulation of greenhouse gases. I highlighted the role of the natural variability that is superimposed to the anthropogenic global warming.