Au cours du 20ème} siècle, les enregistrements marégraphiques suggèrent une hausse du niveau de la mer de 1.8 mm/an. Plus récemment, les observations spatiales indiquent une hausse de 3.3 mm/an sur la période 1993-2009. Cette augmentation au cours du temps est attribuée au réchauffement global de la planète enregistré depuis plusieurs années maintenant. Durant cette thèse, nous analysons les observations et les causes de la hausse moyenne globale du niveau de la mer. Nous estimons les variations stériques du niveau de la mer grâce aux données du projet international Argo et, les variations de masse des océans liés aux apports d'eau des continents à l'aide des mesures de la mission GRACE. Une autre étude se concentre sur la variabilité du stock d'eaux continentales des plus grands bassins hydrologiques de la planète, à l'aide des données GRACE, et l'impact de cette composante à la hausse du niveau de la mer. Puis, nous analysons l'impact de la variabilité interannuelle du stock d'eaux continentales aux variations du niveau de la mer sur diverses périodes. Enfin, nous étudions le bilan des contributions climatiques à la hausse observée du niveau de la mer sur la période altimétrique totale et pour les années récentes. Dans une deuxième partie, nous étudions les variations régionales du niveau de la mer. Nous établissons dans un premier temps, les causes de la variabilité régionale du niveau de la mer et nous interprétons le signal résiduel issu de la différence entre le niveau de la mer observé et l'expansion thermique des océans. Nous étudions ensuite les variations passées des structures régionales du niveau de la mer sur les dernières décennies (1950-2003). Le but est de reconstruire les variations passées du niveau de la mer en 2-D et ainsi d'avoir une connaissance plus approfondie de l'évolution régionale du niveau moyen des mers. Ce travail a pour objectif ultime de contraindre les modèles climatiques couplés utilisés par l'IPCC (Intergovernmental Panel on Climate Change) pour prédire l'évolution future du niveau de la mer au cours du 21ème siècle. Cette analyse nous permet de détecter un signal basse fréquence dans la variabilité régionale du niveau de la mer que nous observons non seulement dans les données in situ mais aussi dans les modèles climatiques couplés.

Au cours du XXème siècle, les mesures marégraphiques ont permis d'estimer la hausse du niveau de la mer global à 1.7 mm.a-1. Depuis deux décennies, les observations faites par les satellites altimétriques indiquent une hausse du niveau de la mer plus rapide, de 3.2 mm. a-1 sur la période 1993-2011. Grâce à leur couverture quasi-globale, les observations spatiales ont aussi révélé une forte variabilité régionale dans la hausse du niveau de la mer qui dépasse de beaucoup la hausse moyenne globale dans de nombreuses régions du globe. Cette composante régionale qui s'ajoute à la hausse globale pour donner le niveau de la mer total local, est essentielle dans l'étude des impacts de la hausse du niveau de la mer sur les régions côtières et les îles basses. Dans cette thèse, nous analysons les observations de la variabilité régionale de la hausse du niveau de la mer, nous proposons une reconstruction de cette variabilité régionale depuis 1950 (i.e. avant l'avènement de l'altimétrie spatiale) et nous étudions ses causes et ses origines. Tout d'abord, nous proposons une reconstruction de la variabilité régionale du niveau de la mer dans le passé (avant la période altimétrique) en combinant des données marégraphiques avec les structures spatiales propres de l'océan déduites des modèles d'océan. Cette méthode permet de reconstruire le niveau de la mer en 2 dimensions depuis 1950, sur la majeure partie du globe, avec une résolution proche de celle de l'altimétrie spatiale. Ensuite, nous appliquons la méthode de reconstruction pour estimer la variabilité régionale de la hausse du niveau de la mer passée dans trois régions sensibles au réchauffement climatique : le Pacifique tropical, la mer Méditerranée et l'océan Arctique. Nous en déduisons pour ces régions la hausse totale ( régionale plus moyenne globale) du niveau de la mer local au cours des dernières décennies. Pour les sites où l'on dispose de mesures du mouvement de la croûte terrestre, nous évaluons la hausse local du niveau de la mer relatif (i.e. hausse du niveau de la mer totale plus mouvement de la croûte local) depuis 1950. Le but est de permettre les études de l'impact local de la hausse du niveau de la mer aux échelles climatiques. Enfin, nous analysons l'origine de la variabilité régionale de la hausse du niveau de la mer pour déterminer si elle est due à l'activité anthropique ou si elle résulte de la variabilité naturelle du système climatique. Nous nous focalisons sur le Pacifque tropical qui est marqué par une très forte variabilité régionale de la hausse du niveau de la mer depuis 1993. Grâce a la reconstruction du niveau de la mer depuis 1950, nous montrons que cette variabilité régionale récente (17 dernières années) n'est pas stationnaire dans le temps mais qu'elle fluctue en lien avec une basse fréquence du mode de variabilité ENSO. Avec les modèles de climat du projet CMIP3, nous montrons de plus que cette variabilité régionale est essentiellement d'origine naturelle (variabilité interne du système climatique) et que l'impact anthropique y est trop faible pour l'instant pour y être détecté.

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

La région du Pamir - Karakoram - Himalaya (PKH) constitue la plus grande réserve de glace terrestre après les régions polaires. Cependant, l'évolution récente de ces glaciers, indicateurs privilégiés du changement climatique en haute altitude, reste encore mal connue, du fait notamment de difficultés d'accès et de conditions climatiques qui rendent délicate l'acquisition de mesures in situ. L'objectif de cette thèse est de contribuer à l'amélioration des connaissances sur l'évolution globale des glaces du PKH au cours de la dernière décennie, en s'appuyant sur des images satellite et des modèles numériques de terrain (MNTs). Une premièreméthodologie a été développée pour assurer le suivi automatique de la distribution spatiale et de l'évolution temporelle des lacs glaciaires à partir d'images Landsat entre 1990 et 2009 sur sept zones d'études réparties le long du PKH. Ainsi, une certaine disparité des types, tailles et évolutions des lacs entre la partie orientale et occidentale du PKH a été mise en évidence. Sur la période de temps considérée, la superficie des lacs a légèrement diminué à l'ouest (Karakoram et Hindu Kush), a été en très nette augmentation à l'est (Népal et Bouthan) et relativement stable sur la partie centrale (Inde du nord-ouest). Le bilan de masse des glaciers a ensuite été calculé, à partir des variations d'épaisseurs mesurées en comparant deuxMNTs, acquis à deux dates différentes, et issus de lamission SRTM et du satellite SPOT5. Cette méthode implique un certain nombre de corrections et d'ajustements au préalable, afin de garantir des mesures les moins biaisées possible. Ainsi, la différence de résolution spatiale initiale des MNTs peut être à l'origine d'un biais fonction de l'altitude, de même que la pénétration des ondes radar de la mission SRTM dans la neige et la glace est à prendre en compte le cas échéant, pour ne pas sous-estimer les altitudes sur les glaciers. Là encore, on observe des disparités entre les différents bilans de masse régionaux sur la période 1999-2011, avec des pertes de masse modérées sur l'Himalaya central et oriental(-0.30±0.08 m a-1 w.e.), plus accentuées sur l'Himalaya occidental (-0.43±0.09 m a-1 w.e.) et des gains de masse plus à l'ouest, pour les glaciers des massifs du Pamir (+0.14±0.11 m a-1 w.e.) et du Karakoram (+0.10±0.20 m a-1 w.e.). Ces résultats confirment donc l'anomalie des glaciers du Karakoram et suggèrent des comportements similaires au Pamir. Le bilan de masse global des glaciers du PKH est estimé à -0.13±0.06 m a-1 w.e.

La région du Pamir - Karakoram - Himalaya (PKH) constitue la plus grande réserve de glace terrestre après les régions polaires. Cependant, l'évolution récente de ces glaciers, indicateurs privilégiés du changement climatique en haute altitude, reste encore mal connue, du fait notamment de difficultés d'accès et de conditions climatiques qui rendent délicate l'acquisition de mesures in situ. L'objectif de cette thèse est de contribuer à l'amélioration des connaissances sur l'évolution globale des glaces du PKH au cours de la dernière décennie, en s'appuyant sur des images satellite et des modèles numériques de terrain (MNTs). Une premièreméthodologie a été développée pour assurer le suivi automatique de la distribution spatiale et de l'évolution temporelle des lacs glaciaires à partir d'images Landsat entre 1990 et 2009 sur sept zones d'études réparties le long du PKH. Ainsi, une certaine disparité des types, tailles et évolutions des lacs entre la partie orientale et occidentale du PKH a été mise en évidence. Sur la période de temps considérée, la superficie des lacs a légèrement diminué à l'ouest (Karakoram et Hindu Kush), a été en très nette augmentation à l'est (Népal et Bouthan) et relativement stable sur la partie centrale (Inde du nord-ouest). Le bilan de masse des glaciers a ensuite été calculé, à partir des variations d'épaisseurs mesurées en comparant deuxMNTs, acquis à deux dates différentes, et issus de lamission SRTM et du satellite SPOT5. Cette méthode implique un certain nombre de corrections et d'ajustements au préalable, afin de garantir des mesures les moins biaisées possible. Ainsi, la différence de résolution spatiale initiale des MNTs peut être à l'origine d'un biais fonction de l'altitude, de même que la pénétration des ondes radar de la mission SRTM dans la neige et la glace est à prendre en compte le cas échéant, pour ne pas sous-estimer les altitudes sur les glaciers. Là encore, on observe des disparités entre les différents bilans de masse régionaux sur la période 1999-2011, avec des pertes de masse modérées sur l'Himalaya central et oriental(-0.30±0.08 m a-1 w.e.), plus accentuées sur l'Himalaya occidental (-0.43±0.09 m a-1 w.e.) et des gains de masse plus à l'ouest, pour les glaciers des massifs du Pamir (+0.14±0.11 m a-1 w.e.) et du Karakoram (+0.10±0.20 m a-1 w.e.). Ces résultats confirment donc l'anomalie des glaciers du Karakoram et suggèrent des comportements similaires au Pamir. Le bilan de masse global des glaciers du PKH est estimé à -0.13±0.06 m a-1 w.e. / The Pamir - Karakoram - Himalaya (PKH) mountain range is considered to be the largest terrestrial ice reservoir outside polar regions. However, the recent evolution of these glaciers, recognized as valuable high-altitude climatic indicators, remains poorly known, mainly because of accessibility issues and harsh meteorologic conditions that hamper field work and in situ observations. The aim of this thesis is therefore to improve the knowledge of glacier changes in PKH and study their evolution over the past decade, based on satellite images and digital elevation models (DEMs). We first developed automatic classification algorithms to monitor the spatial distribution and temporal evolution of glacial lakes between 1990 and 2009 along the PKH, based on Landsat images. We thereby highlight different types, sizes and evolutions of glacial lakes between eastern and western PKH. During the study period the total glacial lake area slightly decreased in the west (Hindu Kush and Karakoram), greatly increased in the east (Nepal, Bhutan) and remained stable in the central part (north-west India). We then computed the mass balance of PKH glaciers from elevation changes measured by comparing two DEMs, acquired ∼10 years appart, by the SRTM mission and the SPOT5 satellite. This method relies on a precise relative adjustment (horizontal and vertical) of the DEMs to remove possible systematic biases within glacier elevation changes. The difference in the original spatial resolution of the DEMs can result in an elevation-dependent bias, as well as the radar penetration into snow and ice can seriously underestimate glacier elevation in the case of a DEM derived from radar data such as SRTM. The spatial pattern of regional glacier mass balances between 1999 and 2011 turns out to be contrasted, with moderatemass losses in eastern and central Himalaya (-0.30±0.08m yr-1 w.e.), stronger in western Himalaya (-0.43±0.09 m yr-1 w.e.) and mass gains further west, for Pamir(+0.14±0.11 m yr-1 w.e.) and Karakoram glaciers (+0.10±0.20 m yr-1 w.e.). The global mass balance of PKH glaciers is estimated at -0.13±0.06 m yr-1 w.e.

Dans quelle mesure les variations du niveau marin actuelles agissent-elles sur la mobilité du trait de côte? Cette question est difficile en raison du caractère lacunaire des jeux de données côtières disponibles. Cette thèse montre tout d'abord qu'il est possible d'évaluer l'élévation relative du niveau de la mer en combinant l'interférométrie radar satellitaire l'analyse de données géodésiques ponctuelles telles que le GPS. Elle examine ensuite le cas de littoraux ayant fait l'expérience d'une élévation du niveau de la mer sensiblement différente de la moyenne globale lors des 50 dernières années. Dans le cas d'atolls de Polynésie Française, les données disponibles montrent le rôle majeur des vagues saisonnières et cycloniques pour contrôler l'évolution du trait de côte. Dans certains secteurs, les effets des vagues sont suspectés se combiner avec l'élévation du niveau de la mer pour favoriser l'érosion de secteurs abrités. Dans le cas de la base de données côtière européenne Eurosion, environ 17 000 observations côtières sont disponibles et peuvent être analysées par une méthode d'exploration de données basée sur les réseaux Bayésiens. Une partie de la variabilité spatiale de l'évolution du trait de côte s'explique par le fait que les littoraux de Fennoscandie sont majoritairement en accrétion, alors qu'ils sont affectés par une baisse du niveau de la mer liée au rebond post-glaciaire. D'une manière générale, ces résultats suggèrent qu'il est encore trop tôt pour observer des effets érosifs évidents de l'élévation du niveau de la mer d'origine climatique. Une analyse probabiliste de l'équation du bilan sédimentaire côtier montre que si cette équation est vérifiée, les effets de l'élévation du niveau de la mer deviendront perceptibles à partir de la seconde moitié du XXIe siècle si les objectifs de réduction des gaz à effet de serre ne sont pas atteints. Finalement, cette thèse ouvre un champ de recherches dans le domaine des mathématiques appliquées pour l'étude des conséquences du changement climatique et de l'élévation du niveau de la mer dans les zones côtières. / To which extent does present day sea-level rise affect shoreline changes? This question remains largely open due to the lack of knowledge and data regarding coastal hydrosedimentary processes and relative sea-level changes (i.e., including vertical ground motions). This PhD Thesis first addresses the question of measuring relative sea-level changes using a combination of geodetic data, including satellite-radar interferometry and permanent GPS stations. Then, it examines the case of coasts that experienced sea-level changes significantly different from the global average over the last 50 years. In the case of atoll islands of French Polynesia, the available data illustrate the major role of southern, trade and cyclonic waves in controlling shoreline changes. For some sheltered coasts, waves are suspected to combine with sea-level rise to favour shoreline retreat. In the case of the European coastal database Eurosion, about 17 000 coastal observations are available and have been analysed using a Bayesian network. Here, part of the spatial variability of shoreline changes can be explained by the fact that Fennoscandian coasts are not only accreting, but also uplifting due to the post-glacial rebound. Overall, these results suggest that it is still too early to observe obvious effects of climate-induced sea-level rise. A probabilistic analysis of the coastal sedimentary equation shows that under common assumptions on hydrosedimentary processes, sea-level rise impacts for beaches should become observable during the second half of the 21st century, if the efforts to maintain climate warming below the 2 C threshold fail. Finally, this PhD thesis offers opportunities for stimulating future research in the field of mathematics applied to the question of climate change and sea-level rise impacts to coastal zones.

L’objectif de la présente thèse est de proposer un cadre conceptuel d’approche desmouvements de population liés aux risques environnementaux. Depuis une vingtained’années, un ensemble de termes a émergé afin de caractériser ces populations: réfugiésenvironnementaux, éco-réfugiés, réfugiés climatiques, éco- migrant pour ne citer que les plusconnus. La caractérisation et le statut accordé à ces populations sont dépendants de la nature,la force et la rapidité de l’évènement environnemental contribuant à leur déplacement. Il estdonc indispensable de définir des critères précis pour éviter tout risque de confusion etprocéder à de véritables évaluations des flux liés à ces catastrophes. Notre analyse sefocalisera sur les risques liés à la hausse du niveau de la mer: risques de submersion etd’inondation sur le littoral et dans les espaces deltaïques. Face à ces risques croissants, il estimpérieux de mettre en place une «stratégie proactive», basée sur le triptyque protection -prévention et anticipation. Cette stratégie exige au-delà de la délimitation des zonesgéographiques concernées, une quantification et qualification des populations potentiellementtouchées. C’est la conjonction de deux formes de vulnérabilité qui est au coeur de notreproblématique: la vulnérabilité face aux risques d’inondations et la vulnérabilité face à lapression démographique croissante. C’est dans cet esprit que nous étudierons deux zonessoumises au risque d’inondation, l’une en France (Estuaire de la Gironde) et l’autre en Grèce(Deltas d’Axios-Loudias-Aliakmonas en Macédoine Centrale) afin d’estimer l’importanceéventuelle des mouvements de population à l’horizon 2025 et 2050. / The present thesis seeks to define a conceptual framework in order to examine the populationmovements related to environmental risks. For about twenty years, several terms haveemerged to describe thedisplaced people such as: environmental refugees, eco-refugees,climate refugees or eco-migrants. The terms and the status attributed to these populations arelargely dependent on the nature, strength and speed of the environmental event contributing totheir displacement. It is therefore necessary to specify the criteria in order to avoid anyconfusion and to be able to assess the physical flows generated by the environmentaldisasters. Our analysis is focused on the risks of submersion and inundation along the coastand in the delta areas associated with rise in the sea levels. In view of these increasing risks, it is absolutely necessary to set up a "proactive strategy", based on the triptych protection -prevention - anticipation. For each type of risk considered, this strategy requires, beyond thedelimitation of the geographical areas concerned, quantification and qualification of thepotentially affected populations. Consequently, the main aspect of our problematic combinestwo types of vulnerability:the firstdue to environmental disasters’ risks such as floods and thesecondto increasing population pressure. Finally, we chose to examine in our study two areashighly subjected to flood risks such as the Gironde estuaryin France and the Axios -Loudias-Aliakmonas Deltas in Central Macedoniain Greece trying to estimate the potential size ofpopulation movements in 2025 and 2050.

La flèche littorale du Cap Ferret s'est édifiée au cours des derniers 3 000 ans, à la faveur du ralentissement de l'élévation du niveau global des océans. Depuis 250 ans, les cartes marines révèlent une importante instabilité de l'extrémité libre de la flèche, à laquelle les travaux présentés dans ce mémoire tentent d'apporter une explication physique. Dans un premier temps le croisement des données cartographiques avec des données environnementales révèle la synchronisation des déplacements de la flèche avec d'une part l'Oscillation Nord Atlantique (NAO), et, d'autre part, les structuations décennales à multidécennales de la vitesse de l'élévation du niveau de l'Atlantique Nord-est. Les processus sous-jacents à cette synchronisation sont alors détaillés avec un modèle numérique de transport sédimentaire. Puis, l'érosion chronique de la flèche entre 1986 et 2016, ainsi que les mécanismes d'évolution de la plage océanique de la Pointe du Cap Ferret sont décrits et quantifiés à l'aide de données topographiques. A l'échelle locale, ces résultats apportent des preuves tangibles du rôle du climat dans l'érosion actuelle de la Pointe du Cap Ferret. D'un point de vue plus général, ces travaux contribuent à la meilleure description des processus contrôlant l'évolution des systèmes de barrières holocènes, dans le contexte actuel d'accélération de l'élévation du niveau de la mer. / Low rates of global sea level rise have allowed the Cap Ferret coastal barrier-spit to develop since about 3 000 years ago. Over the last 250 years, navigation charts attest of a significant instability of the distal end of the spit, to which the present study aims to provide a physical explanation. First, the crossing of cartographic data with meteoceanic data reveals a strong relationship of the spit-end behaviour with, on the one hand, the North Atlantic Oscillation (NAO), and, on the other hand, decadal to multidecadal variations in the rate of sea level rise in the Northeast Atlantic. Underlying processes were then investigated using a numerical model for sediment transport. Finally, the chronic erosion of the spit between 1986 and 2016, as well as the evolution mechanisms of the oceanic beach of the spit-end are described and quantified by means of topographic data. Locally, these results support the role of climate in the ongoing erosion of the spit-end. From a wider perspective, this study brings new insigths on the processes controling the evolution of Holocene barrier systems, in the present context of an accelerating sea level rise.

Ces travaux de thèse s'intéressent à la variabilité océanique intrinsèque générée spontanément sous forçage atmosphérique saisonnier dans des simulations océaniques à haute-résolution (1/4 de degré et 1/12 de degré ), en présence de tourbillons de méso-échelle (environ 100 km). L'empreinte de la variabilité intrinsèque sur le niveau de la mer (SLA) et sur la température de surface (SST) est caractérisée dans plusieurs gammes d'échelles spatio-temporelles par des méthodes de filtrage. En particulier, il est montré que la variabilité petite-échelle inter-annuelle à décennale de SLA observée par les satellites altimétriques, a un caractère stochastique et est majoritairement intrinsèque. Aux échelles de temps multi-décennales, la variabilité intrinsèque de SLA dans les régions océaniques turbulentes est d'amplitude comparable à la variabilité interne simulée dans les modèles climatiques couplés (dénués de turbulence océanique), et pourrait constituer une source additionnelle d'incertitudes dans la détection de l'élévation régionale du niveau de la mer d'origine anthropique. Une analyse spectrale montre enfin que l'advection non-linéaire de vorticité relative transfère spontanément de l'énergie cinétique d'ondes frontales hautes-fréquences vers des tourbillons de méso-échelle plus lents, générant in fine de la variabilité intrinsèque basse-fréquence. Les diagnostiques ont été réalisés à l'aide d'outils développés pour traiter de façon optimale les données océaniques haute-résolution. Ces travaux témoignent plus généralement du caractère chaotique de l'océan turbulent, dont l'interaction avec l'atmosphère est encore mal comprise. / This work focuses on the intrinsic oceanic variability spontaneously generated under seasonal atmospheric forcing in high-resolution oceanic simulations (1/4 ? and 1/12 ? ) including mesoscale eddies (~ 100 km). The imprint of intrinsic variability on sea-level (SLA) and sea surface temperature (SST) is characterized in various spatio-temporal ranges using filtering methods. In particular, the small-scale interannual-to-decadal SLA variability observed by satellite altimetry, is stochastic and mostly intrinsic. At multi-decadal timescales, SLA intrinsic variability is comparable to the internal climate variability simulated in climate coupled models (devoid of oceanic turbulence), and may yield additional uncertainties in the detection of human-induced regional sea-level change. A spectral analysis eventually shows that non-linear advection of relative vorticity spontaneously transfers kinetic energy from high-frequency frontal waves to slower mesoscale eddies, ultimately generating low-frequency intrinsic variability. The analyses have been performed using tools developed for optimal processing of high-resolution oceanic dataset. This work generally provides evidence of the chaotic behavior of the turbulent ocean, whose interaction with the atmosphere is still poorly-known.

Le principal objectif de cette étude est de mieux comprendre et de quantifier les paramètres et les processus de contrôle de la morphologie et de l’évolution stratigraphique du nord du bassin de Santos pour les derniers 49,5 my. Il s’agit de déconvoluer le signal climatique et tectonique dans l’enregistrement sédimentaire à travers de l’analyse sismique à diverses résolutions (~7-8 s et ~400 ms de penetration), de l’analyse lithologique et le calage chrono stratigraphique des limites de séquence majeures, à partir de 19 forages industriels. La modélisation stratigraphique par Dionisos combinée à nos observations nous a permis de quantifier les paramètres majeurs qui vont contrôler l’architecture sédimentaire à différentes échelles de temps. 9 séquences stratigraphiques sont ainsi identifiées et regroupées en 3 mégaséquences (Mégaséquences A, B et C): (i) La Mégaséquence A est caractérisée par des systèmes deltaïques qui viennent remplir l’espace d’accommodation créé par le fluage du sel vers le bassin entre la baie de Guanabara et Cabo Frio. Le flux sédimentaire qui la nourrit est mis en relation avec le paléo fleuve du Paraíba do Sul qui se deplace vers le bassin de Campos à la fin du development de la Mégaséquence et resulte dans une drastique reduction des flux du bassin de Santos; (ii) La Mégasequence B est caractérisée par la rétrogradation de tous les systèmes sédimentaires comme une réponse directe à la diminution des flux sédimentaires. Cette megasequence est aussi caracterisée par l’implantation des depots contouritiques entre 30 et 21 Ma, problablement liés à des changements dans la circulation oceanique dû a l’ouverture de la passage de Drake et de la Tasmanie (entre 37 et 28 Ma); (iii) La Mégasequence C est caractérisée par une progradation et le déplacement des offlap breaks de 70 km vers l’océan, depuis les derniers 11.5 my comme réponse de l’augmentation de flux sédimentaire qui caractérise le passage à des cycles glacio-eustatiques marqués dans le Quaternaire. / The main objective of this study is to better understand and to quantify experimentally the main parameters that control the geomorphological and stratigraphic evolution of the Cenozoic sedimentary succession of the northeast Santos basin, Brazil, in the last 49.5 my. For this purpose, a variety of combining methologies were used, such as: seismic analysis of a dataset composed of different resolution levels (~7-8 s and ~400 ms penetration), chronostratigraphic and lithological data from 19 exploratory wells and stratigraphic modeling. Seismic analysis and stratigraphic modelling led to the identification of 9 main stratigraphic sequences that could be grouped into 3 megasequences (Megasequences A, B and C): (i) Megasequence A is characterized by deltaic features developed into a ramp feature located between Guanabara Bay and Cabo Frio. The formation of this delta resulted in a sedimentary overcharge that induced the basinward salt migration. The sedimentary input that feeds this delta was interpreted as related to an ancient Paraíba do Sul River that diversion to the Campos basin at the end of the megasequence development, resulting in a dramatically reduce the sedimentary influxes into Santos basin; (ii) Megasequence B is characterized by the retrogradational trend of sedimentary systems that offset offlap break of surfaces developed at 30 and 21 Ma about 78 km landward. This megasequence is also characterized by the beginning of conturites deposition related to the changes of oceanic circulation due to the opening of Drake and Tasmanian passages; (iii) Megasequence C is characterized by progradational systems that resulted in the coastline shift about 70 km oceanwards, in the last 11.5 my. The dominant progradational trend of this sequence is interpreted as a response of the increase in the sedimentary influx rate intensified by the spectral characteristics and amplitude of the glacial-induced eustatic signal that characterizes the Quaternary period.