Search for: All records

Abstract The geography of the Southeast Asian Islands (SEAI) has changed over the last 15 million years, as a result of tectonic processes contributing to both increased land area and high topography. The presence of the additional land area has been postulated to enhance convective rainfall, facilitating both increased silicate weathering and the development of the modern‐day Walker circulation. Using an Earth System Model in conjunction with a climate‐silicate weathering model, we argue instead for a significant role of SEAItopographyfor both effects. SEAI topography increases orographic rainfall over land, through intercepting moist Asian‐Australian monsoon winds and enhancing land‐sea breezes. Large‐scale atmospheric uplift over the SEAI region increases by ∼14% as a consequence of increased rainfall over the SEAI and enhancement through dynamical ocean‐atmosphere feedback. The atmospheric zonal overturning circulation over the Indo‐Pacific increases modestly arising from dynamical ocean‐atmosphere feedback, more strongly over the tropical Indian Ocean. On the other hand, the effect of the SEAI topography on global silicate weathering is substantial, resulting in a ∼109 ppm reduction in equilibriumpCO₂and decrease in global mean temperature by ∼1.7ºC. The chemical weathering increase comes from both enhanced physical erosion rates and increased rainfall due to the presence of SEAI topography. The lowering ofpCO₂by SEAI topography also enhances the Indo‐Pacific atmospheric zonal overturning circulation. Our results support a significant role for the progressive emergence of SEAI topography in global cooling over the last several million years. 

Abstract. The numerical model GEOCLIM, a coupled Earth system model for long-term biogeochemical cycle and climate, has been revised. This new version (v 7.0) allows a flexible discretization of the oceanic module, for any paleogeographic configuration, the coupling to any General Circulation Model (GCM), and the determination of all boundary conditions from the GCM coupled to GEOCLIM, notably, the oceanic water exchanges and the routing of land-to-ocean fluxes. These improvements make GEOCLIM7 a unique, powerful tool, devised as an extension of GCMs, to investigate the Earth system evolution at timescales, and with processes that could not be simulated otherwise. We present here a complete description of the model, whose current state gathers features that have been developed and published in several articles since its creation, and some that are original contributions of this article, like the seafloor sediment routing scheme, and the inclusion of orbital parameters. We also present a detailed description of the method to generate the boundary conditions of GEOCLIM, which is the main innovation of the present study. In a second step, we discuss the results of an experiment where GEOCLIM7 is applied to the Turonian paleogeography, with a 10 Myr orbital cycle forcings. This experiment focus on the effects of orbital parameters on oceanic O2 concentration, particularly in the proto-Atlantic and Arctic oceans, where the experiment revealed the largest O2 variations. 

Miocene strata of the Claremont, Orinda, and Moraga formations of the Berkeley Hills (California Coast Ranges, USA) record sedimentation and volcanism during the passage of the Mendocino triple junction and early evolution of the San Andreas fault system. Detrital zircon laser ablation–inductively coupled plasma–mass spectrometry (LA-ICP-MS) age spectra indicate a change in sedimentary provenance between the marine Claremont formation (Monterey Group) and the terrestrial Orinda and Moraga Formations associated with uplift of Franciscan Complex lithologies. A sandstone from the Claremont formation produced a detrital zircon chemical abrasion–isotope dilution–thermal ionization mass spectrometry (CA-ID-TIMS) maximum depositional age of 13.298 ± 0.046 Ma, indicating younger Claremont deposition than previously interpreted. A trachydacite tuff clast within the uppermost Orinda Formation yielded a CA-ID-TIMS U-Pb zircon date of 10.094 ± 0.018 Ma, and a dacitic tuff within the Moraga Formation produced a CA-ID-TIMS U-Pb zircon date of 9.974 ± 0.014 Ma. These results indicate rapid progression from subsidence in which deep-water siliceous sediments of the Claremont formation were deposited to uplift that was followed by subsidence during deposition of terrestrial sediments of the Orinda Formation and subsequent eruption of the Moraga Formation volcanics. We associate the Orinda tuff clast and Moraga volcanics with slab-gap volcanism that followed the passage of the Mendocino triple junction. Given the necessary time lag between triple junction passage and the removal of the slab that led to this volcanism, subsidence associated with ca. 13 Ma Claremont sedimentation and subsequent Orinda to Moraga deposition can be attributed to basin formation along the newly arrived transform boundary. 

The ~2,000-km-long Central Range of New Guinea is a hotspot of modern carbon sequestration due to the chemical weathering of igneous rocks with steep topography in the warm wet tropics. These high mountains formed in a collision between the Australian plate and ophiolite-bearing volcanic arc terranes, but poor resolution of the uplift and exhumation history has precluded assessments of the impact on global climate change. Here, we develop a palinspastic reconstruction of the Central Range orogen with existing surface geological constraints and seismic data to generate time–temperature paths and estimate volumes of eroded material. New (U-Th)/He thermochronology data reveal rapid uplift and regional denudation between 10 and 6 Mya. Erosion fluxes from the palinspastic reconstruction, calibrated for time with the thermochronological data, were used as input to a coupled global climate and weathering model. This model estimates 0.6 to 1.2 °C of cooling associated with the Late Miocene rise of New Guinea due to increased silicate weathering alone, and this CO₂sink continues to the present. Our data and modeling experiments support the hypothesis that tropical arc-continent collision and the rise of New Guinea contributed to Neogene cooling due to increased silicate weathering. 

Large igneous provinces (LIPs) can potentially cause cooling on tens- to thousand-year timescales via injection of sulfur aerosols to the tropo-sphere, and on million-year timescales due to the increase of global weatherability. The ca. 719-Ma Franklin LIP preceded onset of the Sturtian Snowball Earth glaciation by less than two million years, consistent with CO2 drawdown due to weathering of Ca- and Mg-rich LIP basalts, which may have contributed to cooling past a critical runaway ice-albedo threshold. A relatively cool background climate state and Franklin LIP emplacement near a continental margin in the warm wet tropics may have been critical factors for pushing the Earth’s climate past the threshold of runaway glaciation. 

Version 2.0 of GEOCLIM-DynSoil-steady-states with input and configuration files used for the experiments presented in The rise of New Guinea and the fall of Neogene global temperatures; (Martin et al., PNAS, 2023) 

Abstract. The warmer early Pliocene climate featured changes to global sea surface temperature (SST) patterns, namely a reduction in the Equator–pole gradient and the east–west SST gradient in the tropical Pacific, the so-called “permanent El Niño”. Here we investigate the consequences of the SST changes to silicate weathering and thus to atmospheric CO2 on geological timescales. Different SST patterns than today imply regional modifications of the hydrological cycle that directly affect continental silicate weathering in particular over tropical “hotspots” of weathering, such as the Maritime Continent, thus leading to a “weatherability pattern effect”. We explore the impact of Pliocene-like SST changes on weathering using climate model and silicate weathering model simulations, and we deduce CO2 and temperature at carbon cycle equilibrium between solid Earth degassing and silicate weathering. In general, we find large regional increases and decreases in weathering fluxes, and the net effect depends on the extent to which they cancel. Permanent El Niño conditions lead to a small amplification of warming relative to the present day by 0.4 ∘C, suggesting that the demise of the permanent El Niño could have had a small amplifying effect on cooling from the early Pliocene into the Pleistocene. For the reducedEquator–pole gradient, the weathering increases and decreases largely cancel, leading to no detectable difference in global temperature at carbon cycle equilibrium. A robust SST reconstruction of the Pliocene is needed for a quantitative evaluation of the weatherability pattern effect. 

This dataset stores the data of the article The effect of Pliocene regional climate changes on silicate weathering: a potential amplifier of Pliocene-Pleistocene cooling P. Maffre, J. Chiang & N. Swanson-Hysell, Climate of the Past). This study uses a climate model (GCM) to reproduce an estimate of Pliocene Sea Surface Temperature (SST). The main GCM outputs of this modeling (with a slab ocean model) are stored in "GCM_outputs_for_GEOCLIM/", as well as the climatologies from ERA5 reanalysis. The other GCM outputs that were used in intermediary steps (coupled ocean-atmosphere, and fixed SST simulations) are stored in "other_GCM_outputs/". The forcing files (Q-flux) and other boundary conditions to run the "main" GCM simulations can be found in "other_GCM_outputs/Q-flux_derivation/", as well as the scripts used to generate them. Secondly, the mentioned study uses the GCM outputs in "GCM_outputs_for_GEOCLIM/" as inputs for the silicate weathering model GEOCLIM-DynSoil-Steady-State (https://github.com/piermafrost/GEOCLIM-dynsoil-steady-state/tree/PEN), to investigate weathering and equilibrium CO2 changes due to Pliocene SST conditions. The results of these simulations are stored in "GEOCLIM-DynSoil-Steady-State_outputs/". The purpose of this dataset is to provide the raw outputs used to draw the conclusions of Maffre et al. (2023), and to allow the experiments to be reproduced, by providing the scripts to generate the boundary conditions.