ALBERT

Description: During the Archaean, the Sun's luminosity was 18 to 25% lower than the present day. One-dimensional radiative convective models (RCM) generally infer that high concentrations of greenhouse gases (CO2, CH4) are required to prevent the early Earth's surface temperature from dropping below the freezing point of liquid water and satisfying the faint young Sun paradox (FYSP, an Earth temperature at least as warm as today). Using a one-dimensional (1-D) model, it was proposed in 2010 that the association of a reduced albedo and less reflective clouds may have been responsible for the maintenance of a warm climate during the Archaean without requiring high concentrations of atmospheric CO2 (pCO2). More recently, 3-D climate simulations have been performed using atmospheric general circulation models (AGCM) and Earth system models of intermediate complexity (EMIC). These studies were able to solve the FYSP through a large range of carbon dioxide concentrations, from 0.6 bar with an EMIC to several millibars with AGCMs. To better understand this wide range in pCO2, we investigated the early Earth climate using an atmospheric GCM coupled to a slab ocean. Our simulations include the ice-albedo feedback and specific Archaean climatic factors such as a faster Earth rotation rate, high atmospheric concentrations of CO2 and/or CH4, a reduced continental surface, a saltier ocean, and different cloudiness. We estimated full glaciation thresholds for the early Archaean and quantified positive radiative forcing required to solve the FYSP. We also demonstrated why RCM and EMIC tend to overestimate greenhouse gas concentrations required to avoid full glaciations or solve the FYSP. Carbon cycle–climate interplays and conditions for sustaining pCO2 will be discussed in a companion paper.

Description: Paleoelevation reconstructions of mountain belts have become a focus of modern science since surface elevation provides crucial information for understanding both geodynamic mechanisms of Earth’s interior and influence of mountains growth on climate. Stable oxygen isotopes paleoaltimetry is one of the most popular techniques nowadays, and relies on the difference between δ18O of paleo-precipitation reconstructed using the natural archives, and modern measured values for the point of interest. Our goal is to understand where and how complex climatic changes linked with the growth of mountains affect δ18O in precipitation. For this purpose, we develop a theoretical expression for the precipitation composition and we use the isotope-equipped atmospheric general circulation model LMDZ-iso. Experiments with reduced height over the Tibetan Plateau and the Himalayas have been designed. Our results show that the isotopic composition of precipitation is very sensitive to climate changes related with the growth of the Himalayas and Tibetan Plateau, notably changes in relative humidity and precipitation amount. The relative contribution of controlling factors and their magnitude differ depending on the uplift stage and the region considered. Thus future paleoaltimetry studies should take into account constraints on climatic factors to avoid misestimating ancient altitudes.

Description: The 58–51 Ma interval was characterized by a long-term increase of global temperatures (+4 to +6 °C) up to the Early Eocene Climate Optimum (EECO, 52.9–50.7 Ma), the warmest interval of the Cenozoic. It was recently suggested that sustained high atmospheric pCO2, controlling warm early Cenozoic climate, may have been released during Neo-Tethys closure through the subduction of large amounts of pelagic carbonates and their recycling as CO2 at arc volcanoes. To analyze the impact of Neo-Tethys closure on early Cenozoic warming, we have modeled the volume of subducted sediments and the amount of CO2 emitted along the northern Tethys margin. The impact of calculated CO2 fluxes on global temperature during the early Cenozoic have then been tested using a climate carbon cycle model (GEOCLIM). We show that CO2 production may have reached up to 1.55 × 1018 mol Ma−1 specifically during the EECO, ~ 4 to 37 % higher that the modern global volcanic CO2 output, owing to a dramatic India-Asia plate convergence increase. The subduction of thick Greater Indian continental margin carbonate sediments at ~ 55–50 Ma may also have led to additional CO2 production of 3.35 × 1018 mol Ma−1 during the EECO, making a total of 85 % of the global volcanic CO2 outgassed. However, climate modeling demonstrates that timing of maximum CO2 release only partially fits with the EECO, and that corresponding maximum pCO2 values (750 ppm) and surface warming (+2 °C) do not reach values inferred from geochemical proxies, a result consistent with conclusions arising from modeling based on other published CO2 fluxes. These results demonstrate that CO2 derived from decarbonation of Neo-Tethyan lithosphere may have possibly contributed to, but certainly cannot account alone for early Cenozoic warming. Other commonly cited sources of excess CO2 such as enhanced igneous province volcanism also appear to be up to 1 order of magnitude below fluxes required by the model to fit with proxy data of pCO2 and temperature at that time. An alternate explanation may be that CO2 consumption, a key parameter of the long-term atmospheric pCO2 balance, may have been lower than suggested by modeling. These results call for a better calibration of early Cenozoic weathering rates.

Description: The timing of the onset of the Antarctic Circumpolar Current (ACC) is a crucial event of the Cenozoic because of its cooling and isolating effect over Antarctica. It is intimately related to the glaciations occurring throughout the Cenozoic from the Eocene–Oligocene (EO) transition (≈34 Ma) to the middle Miocene glaciations (≈13.9 Ma). However, the exact timing of the onset remains debated with evidence for a late Eocene set up contradicting others data pointing to an occurrence closer to the Oligocene–Miocene (OM) boundary. In this study, we show the potential impact of the Antarctic ice sheet on the initiation of a proto-ACC at the EO boundary. Our results reveal that the regional cooling effect of the ice sheet increases the sea ice formation, which disrupts the meridional density gradient in the Southern Ocean and leads to the onset of a circumpolar current and its progressive strengthening. We also suggest that subsequent variations in atmospheric CO2, ice sheet volumes and tectonic reorganizations may have affected the ACC intensity after the Eocene–Oligocene transition, which in turn may provide an explanation for the second initiation of the ACC at the Oligocene–Miocene boundary and may reconcile evidence supporting both early Oligocene and early Miocene onset of the ACC.

Description: For a long time, evaporitic sequences have been interpreted as indicative of an arid climate. Such systematic interpretations led to the suggestion that the central segment of the South Atlantic (20°–0°) was characterized by an arid climate during the upper Aptian. Indeed, synchronous to this period that corresponds to the rifting and to the opening of this part of the South Atlantic, a large evaporitic sequence spreads out from the equator to 20° S. Using the fully ocean atmosphere coupled model FOAM, we test the potential for the Aptian geography to produce an arid area over the central segment. Sensitivity to the altitude of the rift shoulders separating the Africa and the South America cratons, to the water depth of the central segment and to the drainage pattern have been performed. Using seawater salinity as a diagnostic, our simulations show that the southern part of the central segment is characterized by very high salinity in the case of catchment areas draining the water out of the central segment. Conversely, whatever the boundary conditions used, the northern part of the central segment remains humid and salinities are very low. Hence, we conclude that the evaporites deposited in the southern part of the central segment may have been controlled by the climate favouring aridity and high saline waters. In contrast, the evaporites of the northern part can hardly be reconciled with the climatic conditions occurring there and may be due to hydrothermal sources. Our interpretations are in agreement with the gradient found in the mineralogical compositions of the evaporites from the north to the south, i.e. the northern evaporites are at least 4 times more concentrated than the southern one.

Description: The Ordovician Period (485–443 Ma) is characterized by abundant evidence for continental-sized ice sheets. Modeling studies published so far require a sharp CO2 drawdown to initiate this glaciation. They mostly used non-dynamic slab mixed-layer ocean models. Here, we use a general circulation model with coupled components for ocean, atmosphere, and sea ice to examine the response of Ordovician climate to changes in CO2 and paleogeography. We conduct experiments for a wide range of CO2 (from 16 to 2 times the preindustrial atmospheric CO2 level (PAL)) and for two continental configurations (at 470 and at 450 Ma) mimicking the Middle and the Late Ordovician conditions. We find that the temperature-CO2 relationship is highly non-linear when ocean dynamics are taken into account. Two climatic modes are simulated as radiative forcing decreases. For high CO2 concentrations (≥ 12 PAL at 470 Ma and ≥ 8 PAL at 450 Ma), a relative hot climate with no sea ice characterizes the warm mode. When CO2 is decreased to 8 PAL and 6 PAL at 470 and 450 Ma, a tipping point is crossed and climate abruptly enters a runaway icehouse leading to a cold mode marked by the extension of the sea ice cover down to the mid-latitudes. At 450 Ma, the transition from the warm to the cold mode is reached for a decrease in atmospheric CO2 from 8 to 6 PAL and induces a ~9 °C global cooling. We show that the tipping point is due to the existence of a 95% oceanic Northern Hemisphere, which in turn induces a minimum in oceanic heat transport located around 40° N. The latter allows sea ice to stabilize at these latitudes, explaining the potential existence of the warm and of the cold climatic modes. This major climatic instability potentially brings a new explanation to the sudden Late Ordovician Hirnantian glacial pulse that does not require any large CO2 drawdown.

Description: The timing of the onset of the Antarctic Circumpolar Current (ACC) is a crucial event of the Cenozoic because of its cooling and isolating effect over Antarctica. It is intimately related to the glaciations occurring throughout the Cenozoic from the Eocene–Oligocene (EO) transition (≈ 34 Ma) to the middle Miocene glaciations (≈ 13.9 Ma). However, the exact timing of the onset remains debated, with evidence for a late Eocene setup contradicting other data pointing to an occurrence closer to the Oligocene–Miocene (OM) boundary. In this study, we show the potential impact of the Antarctic ice sheet on the initiation of a strong proto-ACC at the EO boundary. Our results reveal that the regional cooling effect of the ice sheet increases sea ice formation, which disrupts the meridional density gradient in the Southern Ocean and leads to the onset of a circumpolar current and its progressive strengthening. We also suggest that subsequent variations in atmospheric CO2, ice sheet volumes and tectonic reorganizations may have affected the ACC intensity after the Eocene–Oligocene transition. This allows us to build a hypothesis for the Cenozoic evolution of the Antarctic Circumpolar Current that may provide an explanation for the second initiation of the ACC at the Oligocene–Miocene boundary while reconciling evidence supporting both early Oligocene and early Miocene onset of the ACC.

Description: In this contribution, we continue our exploration of the factors defining the Mesozoic climatic history. We improve the Earth system model GEOCLIM designed for long term climate and geochemical reconstructions by adding the explicit calculation of the biome dynamics using the LPJ model. The coupled GEOCLIM-LPJ model thus allows the simultaneous calculation of the climate with a 2-D spatial resolution, the coeval atmospheric CO2, and the continental biome distribution. We found that accounting for the climatic role of the continental vegetation dynamics (albedo change, water cycle and surface roughness modulations) strongly affects the reconstructed geological climate. Indeed the calculated partial pressure of atmospheric CO2 over the Mesozoic is twice the value calculated when assuming a uniform constant vegetation. This increase in CO2 is triggered by a global cooling of the continents, itself triggered by a general increase in continental albedo owing to the development of desertic surfaces. This cooling reduces the CO2 consumption through silicate weathering, and hence results in a compensating increase in the atmospheric CO2 pressure. This study demonstrates that the impact of land plants on climate and hence on atmospheric CO2 is as important as their geochemical effect through the enhancement of chemical weathering of the continental surface. Our GEOCLIM-LPJ simulations also define a climatic baseline for the Mesozoic, around which exceptionally cool and warm events can be identified.

Description: In this contribution, we continue our exploration of the factors defining the Mesozoic climatic history. We improve the Earth system model GEOCLIM designed for long term climate and geochemical reconstructions by adding the explicit calculation of the biome dynamics using the LPJ model. The coupled GEOCLIM-LPJ model thus allows the simultaneous calculation of the climate with a 2-D spatial resolution, the coeval atmospheric CO2, and the continental biome distribution. We found that accounting for the climatic role of the continental vegetation dynamics (albedo change, water cycle and surface roughness modulations) strongly affects the reconstructed geological climate. Indeed the calculated partial pressure of atmospheric CO2 over the Mesozoic is twice the value calculated when assuming a uniform constant vegetation. This increase in CO2 is triggered by a global cooling of the continents, itself triggered by a general increase in continental albedo owing to the development of desertic surfaces. This cooling reduces the CO2 consumption through silicate weathering, and hence results in a compensating increase in the atmospheric CO2 pressure. This study demonstrates that the impact of land plants on climate and hence on atmospheric CO2 is as important as their geochemical effect through the enhancement of chemical weathering of the continental surface. Our GEOCLIM-LPJ simulations also define a climatic baseline for the Mesozoic, around which exceptionally cool and warm events can be identified.

Description: The 58–51 Ma interval was characterized by a long-term increase of global temperatures (+4 to +6 °C) up to the Early Eocene Climate Optimum (EECO, 52.9–50.7 Ma), the warmest interval of the Cenozoic. It was recently suggested that sustained high atmospheric pCO2, controlling warm early Cenozoic climate, may have been released during Neo-Tethys closure through the subduction of large amounts of pelagic carbonates and their recycling as CO2 at arc volcanoes ("carbonate subduction factory"). To analyze the impact of Neo-Tethys closure on early Cenozoic warming, we have modeled the volume of subducted sediments and the amount of CO2 emitted at active arc volcanoes along the northern Tethys margin. The impact of calculated CO2 fluxes on global temperature during the early Cenozoic have then been tested using a climate carbon cycle model (GEOCLIM). We first show that CO2 production may have reached up to 1.55 × 1018 mol Ma−1 specifically during the EECO, ~ 4 to 37 % higher that the modern global volcanic CO2 output, owing to a dramatic India–Asia plate convergence increase. In addition to the background CO2 degassing, the subduction of thick Greater Indian continental margin carbonate sediments at ~ 55–50 Ma may also have led to additional CO2 production of 3.35 × 1018 mol Ma−1 during the EECO, making a total of 85 % of the global volcanic CO2 outgassed. However, climate modelling demonstrates that timing of maximum CO2 release only partially fit with the EECO, and that corresponding maximum pCO2 values (750 ppm) and surface warming (+2 °C) do not reach values inferred from geochemical proxies, a result consistent with conclusions arise from modelling based on other published CO2 fluxes. These results demonstrate that CO2 derived from decarbonation of Neo-Tethyan lithosphere may have possibly contributed to, but certainly cannot account alone for early Cenozoic warming, including the EECO. Other commonly cited sources of excess CO2 such as enhanced igneous province volcanism also appear to be up to one order of magnitude below fluxes required by the model to fit with proxy data of pCO2 and temperature at that time.