Low atmospheric CO2 levels before the rise of forested ecosystems

Nature Communications volume 13, Article number: 7616 (2022) Cite this article


The emergence of forests on Earth (~385 million years ago, Ma)1 has been linked to an order-of-magnitude decline in atmospheric CO2 levels and global climatic cooling by altering continental weathering processes, but observational constraints on atmospheric CO2 before the rise of forests carry large, often unbound, uncertainties. Here, we calibrate a mechanistic model for gas exchange in modern lycophytes and constrain atmospheric CO2 levels 410–380 Ma from related fossilized plants with bound uncertainties of approximately ±100 ppm (1 sd). We find that the atmosphere contained ~525–715 ppm CO2 before continents were afforested, and that Earth was partially glaciated according to a palaeoclimate model. A process-driven biogeochemical model (COPSE) shows the appearance of trees with deep roots did not dramatically enhance atmospheric CO2 removal. Rather, shallow-rooted vascular ecosystems could have simultaneously caused abrupt atmospheric oxygenation and climatic cooling long before the rise of forests, although earlier CO2 levels are still unknown.


Atmospheric CO2 is a greenhouse gas that has affected Earth’s climate throughout geological history2,3. Its variation in the past informs us about the natural long-term sources and sinks. In the absence of anthropogenic fossil fuel combustion, the dominant atmospheric CO2 source is volcanic outgassing, and this source is balanced mainly by the removal that occurs when CO2-bearing fluids chemically react and weather silicate rocks followed by deposition of carbonate in the oceans4. The dissolution of silicate minerals in the weathering zone occurs via interactions between the terrestrial ecosystem and geological processes that make fresh rock available at the surface for reaction. Yet, the role of biology and the CO2-sensitivity of the feedbacks governing global CO2 removal is debated5,6,7,8,9. Enhanced continental weathering is suggested to have caused a decline in atmospheric CO2 pressure (pCO2) from a level ~10 times higher than today’s concentration3,10,11 to near modern levels linked to the Devonian-Carboniferous transition from greenhouse to icehouse conditions in response to the afforestation of the continents. This process is accompanied by burial and preservation of organic matter that also influence atmospheric CO2 levels and, in turn, acts as the main long-term source of atmospheric O2. However, recent geochemical evidence and Earth system models12 suggests atmospheric oxygenation occurred well before trees evolved on the continents ~393–383 Ma1,3. Further, the temporal correlation between plant colonization and the Permo-Carboniferous glaciation has been disputed8. There is compelling evidence that Earth also transitioned into a glaciated state in the Ordovician-Silurian13. Yet, the link between glaciation and atmospheric CO2 is complicated, and palaeoclimate models shows that glaciations could persist even at 12–14 times pre-industrial atmospheric levels (PIAL, 280 ppmv)14. Therefore, a precise reconstruction of atmospheric pCO2 in relation to plant evolution is key to assess the impact of the terrestrial biota on Earth’s climate. Here, we show that atmospheric pCO2 was markedly lower than previously thought when trees and forests appeared on our planet.

In the canonical view, atmospheric CO2 concentrations were one order of magnitude above pre-industrial levels in the early Palaeozoic3,10,15,16,17, although more recent studies suggest levels much closer to today18. Traditionally, palaeo-CO2 estimates from proxy data come with large and sometimes unbounded uncertainty19,20. For example, in the Late Ordovician (~445 Ma), there is evidence for high pCO2 levels of 17 ± 4 PIAL (1 sd, standard deviation) from CO2 hosted in pedogenic goethite from the Neda Formation in Wisconsin, USA19. This paleosol shows a coupling between CO2 content in goethite and its isotope composition interpreted to reflect variable mixing of atmospheric CO2 and soil respired carbon substituted into the goethite mineral lattice. The range above reflects all analytical errors propagated through the calculation (see details in the supplementary information, SI), and we note that the error could still be larger if the CO2 surface adsorption properties21 on natural goethite deviate from that of phosphated goethite grown in the laboratory22. This has never been verified in modern soils, and the proxy has also never been applied at any other time in Earth history.

Further, a systematic decline of atmospheric CO2 levels through the Devonian from 5 to 0.7 PIAL has been inferred from the carbon isotope compositions of pedogenic carbonate2,16,17. Although the declining trend may be real, the absolute atmospheric CO2 level reported from pedogenic carbonates from this time interval have been adjusted further down23 and are likely systematically overestimated because of lower productivity in early Palaeozoic soils relative to modern soils17. Propagating the uncertainty associated with assumed model parameters that cannot be independently constrained from the rock record, such as the proportion of soil-respired CO2 in the soil and its isotope signature, shows that the absolute palaeo-pCO2 estimate obtained this way carries one order of magnitude uncertainty (see supplementary Fig. 15 and the supplementary information for details).

In addition, the first reported evidence for high Palaeozoic CO2 levels (~16 PIAL)20 comes from low stomatal density (mm−2) in some fossil plants (i.e. AglaophytonSawdonia). These anatomical features were interpreted as evidence that these plants had adapted to minimize water loss in a high CO2 atmosphere, but it is essential to compare to plants with similar gas exchange anatomy and behavioural control over water loss rate and CO2 uptake. This becomes very problematic with this group of enigmatic early vascular plants that lack living descendants. In comparison fossils of lycophytes that co-occur with these extinct plant groups that do have modern relatives with similar physiology display similar stomatal density as their modern descendants (Supplementary table 5) suggesting that the high CO2 predictions inferred from Aglaophyton and Sawdonia specimens could be erroneous and that early Palaeozoic atmospheric CO2 was much closer to the modern level24.

Recently, atmospheric CO2 levels were found to be only modestly elevated in the late Ordovician (~400–700 ppm) and Mid-Devonian (~700–1400 ppm) based on carbon isotope data of marine phytoplankton recorded via a diagenetic product of chlorophyll, phytane18. No phytane-CO2 data is reported between ~432 Ma and 390 Ma (supplementary Figs. 1921). The phytane proxy shows the expected response at elevated CO2 concentrations and predicts declining CO2 levels as a function of distance to a modern CO2 seep25 (albeit with larger errors than reported from Palaeozoic phytane). When applied to modern phytoplankton, the phytane proxy predicts a wide range of atmospheric CO2 levels even today (~300–1200 ppm)26. The accuracy of this proxy depends on variables that are not easily detectable from the geological record, including the ratio of atmospheric CO2 to dissolved CO2 of the seawater in which the phytoplankton grew, the nature of the phytoplankton species, other sources of phytane, the growth rate of phytoplankton and, thus, the isotopic response to ambient dissolved CO2 concentration of the phytoplankton from which phytane was derived26. Given the many uncertainties, the errors associated with the phytane CO2-proxy are likely large.

By the same fundamental principle plus an opportunity to constrain all involved parameters, carbon isotope fractionation in plant tissue from terrestrial plants (Δleaf) is sensitive to the internal CO2 concentration in the substomatal cavities that, in turn, depends on water availability and ambient pCO2 in the environment27. Recently, the combination of isotope data and stomatal parameters has led to the development of a mechanistic proxy for prediction of palaeo-pCO2 based around leaf-gas exchange27. Applying this approach to the fossil record shows that post-Devonian atmospheric CO2 levels were <1000 ppm most of the time, but CO2 estimates from the Lower and Middle Devonian using this model show considerable variation with estimates ranging from 530–2853 ppm with one outlier reaching 7320 ppm27.

Importantly, this mechanistic proxy has also allowed for both the accuracy and the precision of the prediction to be tested allowing for the development of realistic constraints around the predicted palaeo-pCO2 to be assessed. Thus, the framework that underpins the proxy allows for the re-evaluation of CO2 predictions from both long-term carbon cycle models which have large uncertainty bracketing their predictions11,12,15 and predictions from first generation palaeo-pCO2 proxies that generated estimates in pCO2, which were largely unconstrained.

The mechanistic model by Franks et al.27 uses both carbon isotope data and stomatal features (density and size) obtained from the fossils to predict palaeo-pCO2. Both characteristics are an expression of how the plant, when alive, was adapted to its local environment. Consequently, there is a strong relationship between these different “recorders” of plant climate interactions. For example, the carbon isotope signature of leaf tissue has been demonstrated to be strongly influenced by water availability28 and the utility of mechanistic proxies underpinned solely by isotope data29,30,31 to be questioned32,33,34,35,36. Due to the co-dependency of isotope and stomatal parameters it is recommended that these data are ideally extracted from the same sample whenever possible or from contemporaneous sedimentary deposits. In taking this approach, the accuracy of the prediction should be improved because it simultaneously accounts for variable CO2 assimilation rate and water conductance in response to changing water availability27. Here, we calibrate a leaf gas-exchange model for lycophytes and estimate atmospheric CO2 levels with uncertainties by applying it to some of the oldest representatives found in the fossil record.

Results and discussion

Calibration of a robust palaeo-CO2 barometer for lycophytes

To anchor the early atmospheric CO2 records, we calibrated a mechanistic leaf-gas exchange model27 in living representatives of the most ancient vascular plant taxa (Lycopodiaceae) constrained by their stomatal density, stomata pore size, and carbon isotope compositions, and applied it to fossil remains of Lower and Middle Devonian lycophytes with the same stomatal anatomy37. When trying to use the plant fossil record to estimate lower Palaeozoic CO2 levels at the time when vascular plants underwent their initial radiation, it becomes challenging to use paired sets of stomatal and isotope data. This situation arises due to the mode of fossil preservation, the need for destructive analysis for carbon isotopes, and the intrinsic scientific value of curated specimens. Therefore, we have further developed the mechanistic model27 to incorporate a solver routine that allows for the propagation of a wide range for each of the observable parameters required to deliver palaeo-pCO2 estimates. This approach has been developed as it allows for all variation between local environments and/or a temporal mismatch between data sources to be folded into the model predictions. Again, to ground truth that our solver approach delivers reasonable estimates of palaeo-pCO2, we have benchmarked the approach by applying the technique to our extant lycophyte dataset.

Initially, we calibrated the gas exchange model27 using two lycophyte species (Huperzia phlegmaria and H. squarrosa) presumed physiologically similar to the Devonian lycophytes and grown at known ambient CO2 levels for ~8 years in natural light within a soil substrate and under optimal Relative Humidity (RH) of ~80% in a greenhouse of the Botanical Garden in Copenhagen. Measurements of the leaf carbon isotope fractionation (Δleaf), stomatal density (SD), and stomatal pore length (p) (Supplementary table 2; Supplementary Data 1) suggest that their CO2 assimilation rate at reference CO2 (A0) was indistinguishable from that of other modern lycophytes (A0 = 3.7 ± 1.6 μmol m−2 s−1, 6 species). Using this average A0 value as representative for lycophytes, and other parameters fixed, such as the operational stomatal conductance efficiency (ζ) with small effect on pCO2 estimates (see supplementary text for details), we were able to reproduce measured CO2 levels of 448 ± 51 ppm (daily variation, 1 sd) in the greenhouse with a CO2 prediction of 407+24−22407−22+24 ppm and 545+69−56545−56+69 ppm (median ± quartiles) for the two distinct species of lycophytes, respectively.

Following this initial screening we further explored the mechanistic model to derive estimates of CO2 by developing a dataset based on non-paired data; i.e., where stomata and isotope data comes from H. phlegmaria and H. squarrosa from different growth conditions to mimic fossil samples that come from geographically/environmentally distinct sedimentary deposits of similar age. Based on cultivated lycophytes grown under dry and humid conditions, we find that unpaired data could carry up to ~100 ppm additional uncertainty, if stomatal and isotopic data were derived from plants that lived under vastly different conditions; i.e., stomatal data from slow-growing or even semi-epiphytic plants under extremely dry conditions with isotope data from plants grown under optimal, humid conditions or vice versa (Supplementary Data 1). When explored using this data framework, the revised predictions of pCO2 shows good agreement with both the actual growth CO2 and values of CO2 predicted using the standard model parameterization (Figs. 12; Supplementary Data 1). Thus, the simultaneous analyses of fossil lycophyte Δleaf, SD, and p should deliver both accurate and precise predictions of palaeo-pCO2 irrespective of whether the data are generated from paired samples (stomatal and isotope data from the same locality) or from different localities. Consequently, by applying the mechanistic model27 to lycophyte data, we are now able for the first time to resolve atmospheric pCO2 accurately and precisely prior to the emergence of forests.

Fig. 1: Atmospheric pCO2 levels are derived from leaf carbon isotope fractionation (Δleaf), Stomata Density (SD), and stomata pore length (p).

a Δleaf is determined from the carbon isotope composition (δ13C) of the fossil flora (green), marine carbonates (blue) and calculated atmospheric CO2 (light blue) (Supplementary Table 4). b Plants grown under higher ambient CO2 levels yields a higher Δleaf and/or lower stomata density. The curves represent gas-exchange model calibrated for modern lycophytes with a reference CO2 assimilation rate at modern atmospheric levels (A0) of 3.5 µmol m−2 s−1, an operational stomatal conductance efficiency (ratio of operational to maximal stomatal conductance) ζ = 0.2, stomata pore shape β = 0.6 and p = 18 ± 2 µm comparable to ancient and modern lycophytes (see supplement for details). Here, a fossil lycophyte from the Devonian is used as an example. Devonian lycophytes typically display Δleaf values close to 20‰ and stomata densities between 15 and 25 mm−1 (Supplementary Data 1).

Fig. 2: Error propagation of the mechanistic pCO2 proxy applied to Givetian Drepanophycus sp. from the Hamilton Group, NY, USA.

Probability distributions for model input parameters (Δleaf, SD, p; upper row) are sampled (N = 10,000) to calculate model parameters (middle row: stomata area, amax, the maximal conductance at full daylight, gc,max, and ratio of CO2 concentration in the substomatal cavities to atmospheric CO2, ci/ca😉 and posterior probability distributions for the model output parameters (lower row: atmospheric CO2, ca, CO2 assimilation rate An, and total conductance, gc,tot). Numerical solutions were obtained using Matlab’s fsolve function (see details in supplementary information). The initial guess output parameters are marked with red triangles, and the results do not depend on the initial guess. The median values and errors represented by 1st and 3rd quartiles are reported above each panel for the calculated properties.

Atmospheric CO2 levels at 410.8–382.7 Ma (Pragian through Givetian) were reconstructed from 66 fossil lycophytes representing three distinct genera (Asteroxylon, Baragwanathia and Drepanophycus) from 13 geological deposits at nine geographically distinct localities. From this, we derive four pCO2 estimates from paired data and six predictions from unpaired data (Supplementary Data 1 for location and data details). However, we note that the stomata densities from all localities were remarkably similar (±1.5 mm−2, 1 sd) and there is a degree of similarity in Δleaf (18.3 ± 1.1‰, 1 sd). Using paired data pCO2 predictions are constrained to between 525+139−101525−101+139 ppm and 695+99−73695−73+99 ppm (median ± quartiles of the probability distribution). Unpaired predictions give a range of pCO2 estimates constrained to between 532+77−76532−76+77 ppm and 715+140−102715−102+140 ppm. Our data yields consistent estimates of CO2 regardless of the method used to compile the raw data supporting the predictions (Fig. 3; Supplementary Data 1). Thus, we conclude that these data substantiate that atmospheric CO2 levels were only 1.9–2.6 times above pre-industrial levels during the ~30-million-year time interval when plants evolved tree stature and forests appeared on Earth.

Fig. 3: Summary of plant evolution and evolving atmospheric composition versus age.

Atmospheric pCO2 constraints from 410–380 Ma lycophytes (Asteroxylon, Baragwanathia, and Drepanophycus) with errorbars smaller than the size of the circles. Paired stomata and isotope data from the same locality are highlighted with black circles. Atmospheric pO2 is constrained by charcoal evidence for wildfire since ~423 Ma5455 (red horizontal line) and fossil roots (red circles)56 that sets minimum pO2 levels according to controlled calibrations in the laboratory5684. Modelled evolution of atmospheric pCO2 (blue curve) and atmospheric pO2 (red curve) derived by adjusting the COPSE Reloaded model for the coupled biogeochemical cycles such that continental weathering efficiency of early vegetated ecosystems scales with the physical weathering record of fines in terrestrial deposits46 and outgassing from the Earth scales with continental arc volcanism. The ‘COPSE reloaded’ model are shown for comparison (thin, dashed curves)12. For further details, see supplementary information. The emergence of land plants, vascular plants, arborescent plants with deep root systems and seed plants are shown with thin lines representing their origin by molecular clock estimates and thicker arrows representing fossil occurrences45. Llan Llandovery, W Wenlock, L Ludlow, P Pridoli. The uncertainty of absolute age assignments are defined biostratigraphically, approximately ±1.5 Ma.

Two potentially confounding factors on the magnitude of leaf carbon isotope fractionation are humidity and atmospheric O2 levels. First, plants down-regulate stomatal conductance and increase their water utilization efficiency in drier habitats38. Lycophytes have passive stomata control and can adapt in drier habitats by minimizing stomata density and size and/or survive at lower CO2 assimilation rate, which would then be expressed in lower Δleaf27. We verified this adaptation experimentally using H. squarrosa that grew semi-epiphytically (Supplementary Fig. 6) under drier conditions (RH~60%) resulting in significantly lower Δleaf (13.3 ± 0.2‰) than the Δleaf (20.1 ± 0.9‰) of plants grown under optimal humidity (RH~80%; Supplementary table 2). This effect is also observed on the isotope composition of natural populations of C3 plants, but is significant only in areas where the mean annual rainfall is below ~1000 mm/yr (Supplementary Fig. 10)28. In the drier greenhouse, H. squarrosa plants had a lower, albeit more variable, stomata density (16.1 ± 4.7 vs. 20.1 ± 2.8 mm−2) and a similar pore length (22.8 ± 1.0 µm vs. 21.6 ± 2.7 µm). Although, we do not know that these lycophytes had fully adapted to the drier greenhouse conditions, the pCO2 estimate derived from such plant material yields an under prediction of ambient glasshouse pCO2 level by -128+79−67128−67+79 ppm, which could be explained by a slightly lower operational to maximal stomatal conductance ratio (ζ) than plants grown under ideal (natural habitat) conditions (e.g. ζ of 0.14 vs. 0.20; see supplementary information for details).

To assess this humidity effect further, we used an atmosphere-ocean palaeoclimate model of intermediate complexity (CLIMBER-3α) to evaluate the Devonian climate with 500 ppm of atmospheric CO2, at the lower end of the range of our new CO2 estimate. The model was run to better constrain conditions where the early lycophytes lived (see supplementary information for details) and shows a temperate planet with mean tropical surface air temperatures of 24.1–24.6 °C. The lycophyte floras in Australia, Germany, and China grew in the monsoonal belt with high mean annual rainfall (MAR) exceeding ~1000 mm/yr. But, Asteroxylon (Rhynie Chert, Scotland) and Drepanophycus floras from Maine, USA; New Brunswick and Gaspé, Canada are found in the slightly drier subtropical zones. Due to its coarse spatial resolution, our model cannot accurately capture local variations in rainfall and humidity at the fossil sites, but the model outputs do suggest growth in either humid or relatively humid environments. Importantly, we see no coupling between Δleaf or predicted pCO2 from each fossil localities and the modelled MAR, RH%, or palaeolatitude. Also, the Δleaf data is not suggestive of arid growth conditions. Given these data, there is no indication that the Devonian lycophytes lived under semi-arid/ arid conditions or the palaeo-pCO2 proxy reported here should deliver a biased prediction as a function of a down regulation in stomatal conductance due to growth in an arid environment.

Secondly, atmospheric O2 can also affect the leaf carbon isotope fractionation because photosynthetic CO2 fixation competes with photosynthetic O2 fixation on the Rubisco enzyme39. This effect is more pronounced in modern plants grown at sub-ambient CO2 and super-ambient O2 levels40. The Devonian atmospheric pO2 levels were likely ~15 to 20 atm%12 and not much lower than today, so any O2 effect on our palaeo-pCO2 estimates will likely be small. Controlled growth experiments with vascular plants, including one species of lycophytes (Selaginella kraussiana), show a small positive shift in Δleaf (0.5 ± 0.4‰) when grown at sub-ambient O2 levels (16 atm%)40. This effect has been ascribed to O2-sensitivity of the photorespiratory compensation point41 and has a negligible effect on past atmospheric CO2 estimates compared to the O2-insensitive parameterisation of the mechanistic proxy (Supplementary table 4).

Importantly, our new and more precise CO2 predictions are radically lower than previous reported mean values, but they are within uncertainty envelope of other pCO2 proxies, including the pedogenic carbonate record when considering that Devonian soils were arguably less productive with lower soil respiration rates17 and lower CO2 concentration than modern soils (supplementary Fig. 14).

Reconfiguring the leaf-gas exchange model to solve for CO2 assimilation rate (A0) allows us to explore ecophysiological performance of early land plants that have no living relatives. This approach enables us to contextualize the very low stomatal density of these plants. Model results suggest that these plants had substantially lower CO2 assimilation rates than contemporaneous lycophytes (A0 < 3.7 μmol m−2 s−1); Aglaophyton (0.40 ± 0.10 μmol m−2 s−1), Rhynia (0.62 ± 0.13 μmol m−2 s−1), Horneophyton (0.84 ± 0.21 μmol m−2 s−1), Sawdonia (0.89 ± 0.04 μmol m−2 s−1) and Nothia (1.27 ± 0.29 μmol m−2 s−1) (Supplementary Data 1). Thus, it is plausible to suggest that the low rates of assimilation in this grouping of plants could be a factor in their eventual displacement. The calculation of low assimilation rates suggest that the sporophyte could be physiologically tied to the gametophyte, suggestive of a degree of matrotrophy, as suggested for some species of Cooksonia42.

Atmospheric CO2 and early afforestation

The atmospheric pCO2 curve (Fig. 3) illuminates how the early history of land plants and their colonization of the Earth’s land surface affected the Earth’s climate system. The origin of trees in the Mid-Devonian (393-383 Ma)1 represents a major evolutionary shift in which arborescence evolved independently in three distinct plant lineages (lycopsids, sphenopsids, ferns). Progymnosperm trees towering up to ~30 m diversified and became widespread from palaeo-equatorial to palaeo-boreal latitudes and from seasonally dry to wet habitats1,43. In contrast to the earlier lycophyte-dominated flora of typically ~0.1 m in height and with rhizoids penetrating only few cm into the ground, the evolving trees had both deeper roots for anchoring and a highly advanced root system akin to modern seed plants1. Still, atmospheric CO2 levels remained rather low and may have declined17 by, at most, a few hundred ppm during the Devonian emergence of forests.

The absolute magnitude of the early Palaeozoic pCO2 decline is still poorly constrained, and there is a conflict between some low pCO2 estimates obtained from Mid-Ordovician marine phytane records and the high pCO2 estimates from pedogenic goethite18. If we trust that the Late Ordovician atmosphere had a higher pCO2 than today, at least outside glacial maxima (supplementary Fig. 12), our new results points to a dramatic and abrupt decline from ~9 (unbound) PIAL to ~1.9 ± 0.3 PIAL that took place within a relatively short time interval ~445–410 Ma before forests appeared. Coincident with this, shallow shrub-like vascular plants spread on the continents (e.g. the Eophytic flora44) and a dramatic shift in the physical weathering regime driven by the evolving terrestrial ecosystems is recorded by the retention of fine grained sediment in continental deposits45,46,47. We propose that the earliest vascular vegetation promoted the exposure of more mineral surface area to weathering fluids and amplified global silicate weathering on the continents (far more than subsequent deep-rooted ecosystems could do) owing to a greater nutrient loss from less developed soils, and therefore, forcing a higher weathering demand.

To simulate the effect of the establishment of early vascular plant ecosystems, we used a dynamic model (COPSE Reloaded) for the coupled biogeochemical cycles to predict atmospheric pCO2 and pO2 trajectories (see Methods). Our model predicts a massive atmospheric CO2 decline from ~2500 ppm to ~500 ppm in only ~30 Myr in response to enhanced silicate weathering by early vascular plants (Fig. 3). To exemplify such a scenario, we updated several forcing functions in the latest COPSE model12,48,49 (section S5; Supplementary Figs. 1921). The weathering forcing (W) was adjusted to scale in proportion to the plant-induced effect on mudrock retention in continental deposits normalized to the Carboniferous average46. This is justified because mineral surface area is a key factor facilitating mineral dissolution during chemical weathering. Also, we scaled up plant evolution in concert with the radiation of vascular plants rather than non-vascular plants in previous COPSE models in an attempt to both capture a more extensive plant coverage50,51 and also accounting for selective P weathering to mimic a greater weathering demand of early vascular plants with primitive root-like systems49,50,51. Lastly, volcanic outgassing rates in the Lower and Middle Devonian were adjusted so that the CO2 flux emanating from Earth’s interior scales with subducted carbonate platforms rather than to global seafloor spreading rates as was assumed in previous models11,52. This was done by assuming a linear relationship between outgassing rate, continental arc volcanism and the relative abundance of young to older arc-derived grains in sedimentary deposits53. The Late Silurian to Middle Devonian (~430–370 Ma) outgassing rates is then ~1.0–1.2 times modern levels (in contrast to ~1.5 in previous models) characteristic of ages when supercontinents assembly.

The revised model simultaneously predicts that the evolving vascular ecosystems also caused a rise in atmospheric pO2 consistent with multiple lines of geochemical evidence for Earth’s oxygenation12 including wild fire evidence supporting atmospheric O2 levels above 15 atm% since ~420 Ma54,55 and fossil roots56 (Fig. 3). Previous models have suggested a two-stage transition with relatively high pCO2 and high pO2 in the Silurian12, but our new data and model offers a simpler solution where a synchronous atmospheric O2 rise and CO2 decline to near-modern conditions happened in the same event.

Ultimately, the composition of Earth’s atmosphere is governed by an interplay between biological and geological processes and how land plants and their root symbionts affect the physical and chemical weathering processes on land45,57. The difference in atmospheric CO2 by forestation post-vascular colonization was at most a few hundred ppm58. Thus, global climatic cooling is not strongly linked to afforestation59, but rather to how the terrestrial biosphere extracts and maintains nutrients from its planetary substrate.

The Mid-Palaeozoic climate revisited

The Devonian climate has previously been described mostly as a warm greenhouse that transited into a colder state with polar glaciations in the Late Devonian-Carboniferous10,15. At only 500 ppm CO2, however, our Early Devonian palaeoclimate model predicts a clement climate with global mean annual temperatures of 14.8–15.7 °C for different orbital configurations and a strong latitudinal temperature gradient not too dissimilar from today. Mean tropical surface air temperatures are 24.1–24.6 °C. The predicted sea surface temperatures in tropical Gondwana, South China, Western Laurentia, and Europe are consistent with Early and Middle Devonian palaeotemperature proxy records based on oxygen isotope compositions of phosphatic conodonts and well-preserved calcitic brachiopods (Supplementary Data 2)60,61. Our paleoclimate model with 500 ppm atmospheric CO2 predicts significantly lower temperatures in the subtropical and temperature zone where there are currently no precise palaeotemperature estimates (see supplementary information section S4.2 for further discussion)62. Polar sea ice and snowfall on Gondwana are predicted during winters (Fig. 4, Supplementary Figs. 1416). CLIMBER-3α is not coupled to a dynamic ice sheet model, so the extent and persistence of ice sheets is indirectly constrained by the extent of snow and sea ice. We find that the Earth system was climatically stable and that a cascade into a snowball-style glaciation due to the ice-albedo feedback is highly unlikely at this low CO2 level even under the weaker solar insolation, a finding that is also supported by model simulations bracketing the Early Devonian63,64,65. Therefore, our results suggest that Earth’s climate was in an icehouse state with partial glaciation on the south polar continent Gondwana above 60-80°S palaeolatitude (Fig. 4). We suggest from modern palaeoclimate models14,66,67 that the original conjecture10 that several thousands of ppm CO2 in the atmosphere were necessary to compensate for the ~3% weaker Palaeozoic Sun should be abandoned, in favour of much higher climate sensitivity of CO2 as supported by the new palaeo-CO2 proxy data3.

Fig. 4: Early Devonian (415 Ma) annual mean surface-air temperature for an atmospheric pCO2 of 500 ppmv, modelled with the coupled climate model CLIMBER-3α.

The palaeogeographical configuration of Scotese at 415 Ma83 is used, pre-industrial levels of other greenhouse gasses and idealized orbital parameters are assumed (circular orbit, obliquity of 23.5°; see supplementary Fig. 17 for two different orbital states). See supplement for details about modelled precipitation and humidity. Nine localities are marked where the studied fossil lycophyte flora have been found in 13 distinct outcrops: 1 – Green County, NY, USA and Schoharia Co, NY, USA; 2 – Traveller Mountain, ME, USA; 3 – North Shore and Seal Rock, Gaspe Qbc., Canada; 4 – Abitibi River, Ontario, Canada; 5 – Dalhousie, New Brunswick, Canada (close to Maple Green); 6 – Munchshecke, Siegburg and Mosel Valley, Germany; 7 – Rhynie Chert, Aberdeenshire, Scotland; 8– Yea district, Victoria, Australia; 9 – Xinjan, China. Six localities with palaeo-SST data from calcitic brachiopods are marked with white stars: a – Section Madène el Mrakib, Anti-Atlas Mountains, Morocco; b – Colle area, Cantabrian Mountains, Spain; c – Eifel Mountains, Germany; d – Buffalo Quarry, Iowa, USA; e– Mafeking Quarry, Manitoba, Canada; f– Gur’evsk town, Altai Mts./Salair area, Siberia. Seven localities with palaeo-SST data from phosphatic conodonts are marked with purple stars: g – Changputang, SE Yunnan, China and Nayi, Guangxi, China; h – Broken River, Queensland, Australia; i – Buchan, Boola Quarry, Tyers Quarry, Wellington, Victoria, Australia; j – Mungallala and Windellama, Wellington New South Wales, Australia; k – Barrandov and Na Skrabku, Prague basin, Czech Republic; m – Nevada, USA; n – CaiZiyan, Guangxi, China. Phosphatic SST constraints from d –Glory Quarry, Iowa, USA and b – Puech de la Suque, France are plotted along with the calcitic SST data from nearby outcrops.

Direct evidence for glaciation occurs in Late Devonian – Carboniferous deposits68, but is rare in Lower- and Middle Devonian strata. Indirect evidence of eustatic sea level change occurs in the Lower Devonian sedimentary successions from North Africa that could potentially have a glaciogenic origin69. Further, oxygen isotopic palaeotemperature records60 from the shallow tropical and subtropical oceans are compatible with our predicted palaeoclimate data, but also display considerable local temperature variations61. Also, Devonian plants are mostly found at lower-mid palaeolatitudes (~45°S), where the continent is predicted permanently ice-free. One exception might be the early Lockhovian Cooksonia-rich flora from the Paraná basin in Brazil positioned at high palaeolatitudes (~70°S)70. However, even at 500 ppm CO2, we find that the snow cover on Gondwana was not always perennial (supplementary Fig. 15). Although, the Paraná flora notably predates our record (~419 Ma) when atmospheric CO2 levels were presumably higher, plants could also have grown at high latitudes either during warmer Southern summers or because climate was warmer for other reasons than high CO2 (e.g. higher CH4, N2O).

In summary, the revised atmospheric CO2 record is compatible with all palaeoclimatic constraints from the geological record, with global palaeoclimate models, and with global models for the coupled biogeochemical cycles that include a larger weathering demand of the earliest terrestrial flora and where early afforestation played only a minor role on global CO2 sequestration.



Carbon isotopes, stomata density and stomata pore size from fossil leaves was compiled from the literature (see supplementary information) along with new data collected from the Baragwanathia flora in Victoria, Australia. Plant macrofossils from the fossil collection of Museums Victoria were originally collected by Isabel Cookson from Mt Pleasant Rd, Victoria, Australia71. The fossils are preserved as incrustations in fine-grained sandstone sandwiched in a 130 m thick stratigraphic section of mostly siltstone and shale (supplementary Fig. 1), interpreted as occasional bursts of high-energy turbidites carrying allochthonous fossils from shallower waters into very-low-energy marine depositional environment72. Index fossils (Uncinatograptus sp. cf. U. thomasi and Nowakia sp. ex gr. N. acuaria) confines the flora to the Pragian or earliest Emsian, corresponding to ~409.1 ± 1.5 Ma according to GTS202073.

Small fossil fragments found in four specimens (#15154, #15173, #15174, #15183) which also contain larger fragments of Baragwanathia longifolia and Zosterophyllum australianum, were selected for analyses based on the preservation of black organic matter contained within the brown mineralized fossils (supplementary Figs. 24). From each fragment (37 in total), 0.05–7.34 mg of material was extracted using a scalpel or a 0.8 mm Dremel drill. Only material visibly containing black organic matter was extracted.

Carbon isotope analyses

Carbon isotopic analyses for were performed in the Syracuse University GAPP Lab using an automated ‘nano-EA’ system adapted from that described in Polissar et al.74. The Syracuse University nano-EA comprises an Elementar Isotope Cube elemental analyzer coupled to an Isoprime 100 continuous-flow stable isotope mass spectrometer via an Isoprime Trace Gas analyzer. Though the presence of carbonate was not suspected, nor observed via testing of sample powders with 6 N HCl under a binocular scope, we decided to decarbonate the sample materials prior to analysis to ensure that sample materials were carbonate free. Sample powders were fumigated in ashed glass vials in the presence of neat hydrochloric acid within an evacuated glass bell jar for 24 hours and then were dried in an oven at 40 °C. For isotopic analysis, sample materials were transferred to silver or quartz cups (6 ×6 mm; EA Consumables) and nested in a small amount of quartz wool to ensure retention of sample materials within the cup. The cups and quartz wool were ashed at 480 °C for 8 hours. Sample materials were loaded into cups within a Class 100 laminar flow isolation cabinet with HEPA filtered air to minimize the potential of particulate contamination.

During isotopic analysis, sample cups were evacuated and purged with helium prior to introduction into the EA. Reaction conditions were as follows: oxidation and reduction reactor temperatures were 1100 °C and 650 °C, respectively; helium carrier gas flow was 158 ml/min and the O2 pulse was set for 45 seconds. Carbon dioxide generated during sample combustion was trapped within the EA in a molecular sieve trap. Following passage of the N2 peak, the primary EA trap was heated and carbon dioxide was released to a secondary, silica gel-filled cryotrap which was immersed in liquid nitrogen. Trapping duration was calibrated using the EA thermal conductivity detector data to ensure complete collection of the CO2 peak. Following collection of CO2, the cryotrap gas flow was switched to a lower-flow He carrier gas (~1 mL/min) via an automated Vici Valco 6-port valve. The trap was warmed, and sample gas was released to the IRMS through an Agilent CarboBond capillary chromatography column (25 m x 0.53 mm x 5 µm). The resulting raw carbon isotope data are blank-corrected using direct blank subtraction and normalized to the VPDB scale using the two-point correction scheme75 with the international reference materials NIST 1547-Peach Leaves (δ13C = 26.0 ± 0.2‰) and IAEA C6-sucrose (δ13C = −10.45 ± 0.03‰) which are run as solids. The reference materials USGS 61- Caffeine (−35.05‰) and USGS 62-Caffeine (−14.79‰) are dissolved in UV-treated MilliQ water, dispensed in known quantities, and have long-term laboratory reproducibility of ±0.47‰ for USGS 62 and ±0.39‰ for USGS 61 over a 20–90 nanomole range. Reproducibility of the δ13C values of reference materials with carbon contents greater than 50 nanomoles is ±0.3‰ (1 sd) and is equivalent to that previously reported74.

TOF-SIMS analyses

The morphology and chemical composition of the fragments were characterized non-destructively by SEM and TOF-SIMS to confirm the presence of organic tissue well suited for carbon isotope analysis. Time-Of-Flight Secondary Ion Mass Spectrometry (TOF-SIMS) was used to produce semi-quantitative maps of the elemental composition in sample #15153 (Supplementary Fig. 5).

This technique bombards the surface with Bi ions that causes a collision cascade in the uppermost atom layers of the specimen (~10 nm). This releases secondary ions that are accelerated in an electric field and their time of flight to the detector in vacuum is a function of their mass and sample depth76,77,78. Supplementary Fig. 5 shows elements bound to organic matter (incl. C, N, P) in the sample that produce polyatomic charged species, such as CN- and CNO- when emitted from the same sample depth. The correlations of organic-bound elements enable us to distinguish the presence of organic carbon from inorganic phases (e.g. carbonate minerals) in the sample.

Palaeoclimate modelling

The relatively fast coupled Earth-system model of intermediate complexity CLIMBER-3α was used to simulate the Devonian climate with an atmospheric CO2 level of 500 ppm. CLIMBER-3α encompasses a modified version of the ocean circulation model (MOM379,80) with a horizontal resolution of 3.75° x 3.75° and 24 vertical levels, a dynamic/thermodynamic sea-ice model81 with the same resolution and a fast atmospheric model82 of 22.5° longitudinal and 7.5° latitudinal resolution. The model does not explicitly model ice sheet growth on the continents, but snow cover on the continents is considered. The model was run for Lower Devonian boundary conditions (415 Ma) in terms of continental configuration, solar luminosity and vegetation cover83. Based on previous results67, three different orbital configurations were explored: the standard configuration (obliquity 23.5°, eccentricity e = 0) as well as cold (obliquity 22.0°, eccentricity e = 0) and warm (obliquity 24.5°, eccentricity e = 0.069, precession angle 0°) orbital configurations (supplementary Fig. 17). A sensitivity analysis considering seasonal surface air temperatures and sea-ice distribution for these different insolation patterns is shown in Supplementary table 6. For determining climate variables at specific proxy locations (see Supplementary Data 2) the simulated values were bilinearly interpolated on a 1°-by-1° grid and evaluated using coordinates transformed from the present-day values using GPlates.

Long-term global biogeochemical modelling

We used the Carbon-Oxygen-Phosphorous-Sulphur Evolution (COPSE) model to predict the histories of atmospheric pCO2pO2 and ocean composition over the Phanerozoic (550 Ma–today). This forward modelling approach enables hypothesis testing of mechanistic cause-effect relationships in the Earth system. A set of coupled differential equations describing the dynamic evolution of the C, O, P and S cycles were solved using an inbuilt variable timestep solver for ‘stiff’ Ordinary Differential Equation systems in Matlab®. The most recent version of the model (COPSE Reloaded, denoted ‘CR’; Supplementary Figs. 1718) were adapted and modified48,49. As input, we let the C/P ratio of buried terrestrial biomass increase with the colonization of land by non-vascular plants as in CR and adjusted the forcings on plant weathering (W), plant evolution (E), selective P weathering (F), volcanic outgassing (D) in a manner to simultaneously fit the effect of shallow vascular ecosystems on weathering processes and produce outputs consistent with the palaeorecords. Figure 3 shows the atmospheric pO2 and pCO2 trajectories predicted by the revised version of CR with forcings shown in supplementary Fig. 21. Further details on the revised COPSE modelling are found in the supplementary text section S5.

Data availability

All data are available in the main text or the supplementary materials.

Code availability

The models developed in this analysis are made available at University of Copenhagen’s Electronic Research Data Archive (ERD). https://doi.org/10.17894/ucph.214a6434-b7eb-4e62-aaac-afefc1247da4. The source code for the CLIMBER paleoclimate model used in this study is archived at the Potsdam Institute for Climate Impact Research and is made available upon request. The postprocessing scripts used to analyse the CLIMBER data and to generate the paleoclimate model figures presented in the study are stored on ERDA.


  1. Stein, W. E. et al. Mid-devonian archaeopteris roots signal revolutionary change in earliest fossil forests. Curr. Biol. 30, 421–431.e2 (2020).Article CAS Google Scholar 
  2. Foster, G. L., Royer, D. L. & Lunt, D. J. Future climate forcing potentially without precedent in the last 420 million years. Nat. Commun. 8, 14845 (2017).Article ADS CAS Google Scholar 
  3. Berner, R. A. Paleozoic atmospheric CO2: importance of solar radiation and plant evolution. Science 261, 68–70 (1993).Article ADS CAS Google Scholar 
  4. Walker, J. C. G., Hays, P. B. & Kasting, J. F. A negative feedback mechanism for the long-term stabilization of Earth’s surface temperature. J. Geophys. Res. 86, 9776 (1981).Article ADS CAS Google Scholar 
  5. D’Antonio, M., Ibarra, D. E. & Boyce, C. K. Land plant evolution decreased, rather than increased, weathering rates. Geology https://doi.org/10.1130/G46776.1 (2019).Article Google Scholar 
  6. Caves, J. K., Jost, A. B., Lau, K. V. & Maher, K. Cenozoic carbon cycle imbalances and a variable weathering feedback. Earth Planet. Sci. Lett. 450, 152–163 (2016).Article ADS CAS Google Scholar 
  7. Maher, K. & Chamberlain, C. P. Hydrologic regulation of chemical weathering and the geologic carbon cycle. Science 343, 1502–1504 (2014).Article ADS CAS Google Scholar 
  8. Goddéris, Y. et al. Onset and ending of the late Palaeozoic ice age triggered by tectonically paced rock weathering. Nat. Geosci. 10, 382–386 (2017).Article ADS Google Scholar 
  9. Ibarra, D. E. et al. Modeling the consequences of land plant evolution on silicate weathering. Am. J. Sci. 319, 1–43 (2019).Article ADS CAS Google Scholar 
  10. Berner, R. A. Atmospheric carbon dioxide levels over phanerozoic time. Science 249, 1382–1386 (1990).Article ADS CAS Google Scholar 
  11. Berner, R. A. GEOCARBSULF. A combined model for Phanerozoic atmospheric O2 and CO2Geochimica et Cosmochimica Acta 70, 5653–5664, (2006).
  12. Lenton, T. M. et al. Earliest land plants created modern levels of atmospheric oxygen. Proc. Natl. Acad. Sci. USA 113, 9704–9709 (2016).Article ADS CAS Google Scholar 
  13. Ghienne, J.-F. et al. A Cenozoic-style scenario for the end-Ordovician glaciation. Nat. Commun. 5, 4485 (2014).Article ADS CAS Google Scholar 
  14. Pohl, A. et al. Glacial onset predated Late Ordovician climate cooling. Paleoceanography 31, 800–821 (2016).Article ADS Google Scholar 
  15. Royer, D. L. Stomatal density and stomatal index as indicators of paleoatmospheric CO2 concentration. Rev. Palaeobot. Palynol. 114, 1–28 (2001).Article Google Scholar 
  16. Mora, C. I., Driese, S. G. & Colarusso, L. A.Middleto Late Paleozoic Atmospheric CO$_2$ Levels from Soil Carbonate and Organic Matter. Sci. New Ser. 271, 1105–1107 (1996).CAS Google Scholar 
  17. Ekart, D. D. A 400 million year carbon isotope record of pedogenic carbonate; implications for paleoatomospheric carbon dioxide. Am. J. Sci. 299, 805–827 (1999).Article ADS CAS Google Scholar 
  18. Witkowski, C. R., Weijers, J. W. H., Blais, B., Schouten, S. & Sinninghe Damsté, J. S. Molecular fossils from phytoplankton reveal secular P co 2 trend over the Phanerozoic. Sci. Adv. 4, eaat4556 (2018).Article ADS CAS Google Scholar 
  19. Yapp, C. J. & Poths, H. Ancient atmospheric CO2 pressures inferred from natural goethites. Nature 355, 342–344 (1992).Article ADS CAS Google Scholar 
  20. McElwain, J. & Chaloner, W. Stomatal density and index of fossil plants track atmospheric carbon dioxide in the palaeozoic. Ann. Bot. 76, 389–395 (1995).Article Google Scholar 
  21. Yapp, C. J. A possible goethite-iron(III) carbonate solid solution and the determination of CO2 partial pressures in low-temperature geologic systems. Chem. Geol. 64, 259–268 (1987).Article ADS CAS Google Scholar 
  22. Russell, J. D., Paterson, E., Fraser, A. R. & Farmer, V. C. Adsorption of carbon dioxide on goethite (α-FeOOH) surfaces, and its implications for anion adsorption. J. Chem. Soc. Faraday Trans. 1: Phys. Chem. Condens. Phases 71, 1623 (1975).Article CAS Google Scholar 
  23. Breecker, D. O., Sharp, Z. D. & McFadden, L. D. Atmospheric CO2 concentrations during ancient greenhouse climates were similar to those predicted for A.D. 2100. Proc. Natl. Acad. Sci. 107, 576–580 (2010).Article ADS CAS Google Scholar 
  24. Edwards, D. Stomata in early land plants: an anatomical and ecophysiological approach. J. Exp. Bot. 49, 255–278 (1998).Article Google Scholar 
  25. Witkowski, C. R., van der Meer, M. T. J., Smit, N. T., Sinninghe Damsté, J. S. & Schouten, S. Testing algal-based pCO2 proxies at a modern CO2 seep (Vulcano, Italy). Sci. Rep. 10, 10508 (2020).Article ADS CAS Google Scholar 
  26. Freeman, K. H. & Hayes, J. M. Fractionation of carbon isotopes by phytoplankton and estimates of ancient CO 2 levels. Glob. Biogeochem. Cycles 6, 185–198 (1992).Article ADS CAS Google Scholar 
  27. Franks, P. J. et al. New constraints on atmospheric CO2 concentration for the Phanerozoic. Geophys. Res. Lett. 41, 4685–4694 (2014).Article ADS CAS Google Scholar 
  28. Diefendorf, A. F., Mueller, K. E., Wing, S. L., Koch, P. L. & Freeman, K. H. Global patterns in leaf 13C discrimination and implications for studies of past and future climate. Proc. Natl Acad. Sci. 107, 5738–5743 (2010).Article ADS CAS Google Scholar 
  29. Schubert, B. A. & Jahren, A. H. Global increase in plant carbon isotope fractionation following the Last Glacial Maximum caused by increase in atmospheric pCO2Geology 43, 435–438 (2015).Article ADS CAS Google Scholar 
  30. Schubert, B. A. & Jahren, A. H. The effect of atmospheric CO2 concentration on carbon isotope fractionation in C3 land plants. Geochimica et. Cosmochimica Acta 96, 29–43 (2012).Article ADS CAS Google Scholar 
  31. Cui, Y. & Schubert, B. A. A 23 m.y. record of low atmospheric CO2Geology 48, 888–892 (2020).Article ADS CAS Google Scholar 
  32. Diefendorf, A. F., Leslie, A. B. & Wing, S. L. Leaf wax composition and carbon isotopes vary among major conifer groups. Geochimica et. Cosmochimica Acta 170, 145–156 (2015).Article ADS CAS Google Scholar 
  33. Lomax, B. H., Lake, J. A., Leng, M. J. & Jardine, P. E. An experimental evaluation of the use of Δ13C as a proxy for palaeoatmospheric CO2Geochimica et. Cosmochimica Acta 247, 162–174 (2019).Article ADS CAS Google Scholar 
  34. Stein, R. A., Sheldon, N. D. & Smith, S. Rapid response to anthropogenic climate change by Thuja occidentalis: implications for past climate reconstructions and future climate predictions. PeerJ 7, e7378 (2019).Article Google Scholar 
  35. Stein, R. A., Sheldon, N. D. & Smith, S. Y. C3 plant carbon isotope discrimination does not respond to CO2 concentration on decadal to centennial timescales. N. Phytologist 229, 2576–2585 (2021).Article CAS Google Scholar 
  36. Jardine, P. E. & Lomax, B. H. A 23 m.y. record of low atmospheric CO2Geology 49, e523 (2021).Article ADS Google Scholar 
  37. Franks, P. J. & Farquhar, G. D. The mechanical diversity of stomata and its significance in gas-exchange control. Plant Physiol. 143, 78–87 (2007).Article CAS Google Scholar 
  38. Kohn, M. J. Carbon isotope compositions of terrestrial C3 plants as indicators of (paleo)ecology and (paleo)climate. Proc. Natl Acad. Sci. 107, 19691–19695 (2010).Article ADS CAS Google Scholar 
  39. Fletcher, B. J., Brentnall, S. J., Quick, W. P. & Beerling, D. J. BRYOCARB: A process-based model of thallose liverwort carbon isotope fractionation in response to CO2, O2, light and temperature. Geochimica et. Cosmochimica Acta 70, 5676–5691 (2006).Article ADS CAS Google Scholar 
  40. Porter, A. S., Yiotis, C., Montañez, I. P. & McElwain, J. C. Evolutionary differences in Δ13C detected between spore and seed bearing plants following exposure to a range of atmospheric O2:CO2 ratios; implications for paleoatmosphere reconstruction. Geochimica et. Cosmochimica Acta 213, 517–533 (2017).Article ADS CAS Google Scholar 
  41. Porter, A. S. et al. Testing the accuracy of new paleoatmospheric CO2 proxies based on plant stable carbon isotopic composition and stomatal traits in a range of simulated paleoatmospheric O2:CO2 ratios. Geochimica et. Cosmochimica Acta 2019, 69–90 (2019).Article ADS Google Scholar 
  42. Boyce, C. K. How green was Cooksonia? The importance of size in understanding the early evolution of physiology in the vascular plant lineage. Paleobiology 34, 179–194 (2008).Article Google Scholar 
  43. Streel, M. Late Frasnian–Famennian climates based on palynomorph analyses and the question of the Late Devonian glaciations. Earth-Sci. Rev. 52, 121–173 (2000).Article ADS Google Scholar 
  44. Cleal, C. & Cascales-Miñana, B. Evolutionary Floras – revealing large-scale patterns in Palaeozoic vegetation history. J. Palaeosci. 28, 1–20 (2021).
  45. Dahl, T. W. & Arens, S. K. M. The impacts of land plant evolution on Earth’s climate and oxygenation state – An interdisciplinary review. Chem. Geol. 547, 119665 (2020).Article ADS CAS Google Scholar 
  46. McMahon, W. J. & Davies, N. S. Evolution of alluvial mudrock forced by early land plants. Science 359, 1022–1024 (2018).Article ADS CAS Google Scholar 
  47. Libertín, M., Kvaček, J., Bek, J., Žárský, V. & Štorch, P. Sporophytes of polysporangiate land plants from the early Silurian period may have been photosynthetically autonomous. Nat. Plants 4, 269–271 (2018).Article Google Scholar 
  48. Lenton, T. M., Daines, S. J. & Mills, B. J. W. COPSE reloaded: An improved model of biogeochemical cycling over Phanerozoic time. Earth-Sci. Rev. 178, 1–28 (2018).Article ADS CAS Google Scholar 
  49. Krause, A. J. et al. Stepwise oxygenation of the Paleozoic atmosphere. Nat. Commun. 9, 4081 (2018).Article ADS Google Scholar 
  50. Porada, P. et al. High potential for weathering and climate effects of non-vascular vegetation in the Late Ordovician. Nat. Commun. 7, 12113 (2016).Article ADS CAS Google Scholar 
  51. Halder, S., Arens, S. K. M., Jensen, K., Dahl, T. W. & Porada, P. A dynamic local-scale vegetation model for lycopsids (LYCOm v1.0). Geosci. Model Dev. 15, 2325–2343 (2022).Article ADS CAS Google Scholar 
  52. Bergman, N. M. COPSE: A new model of biogeochemical cycling over Phanerozoic time. Am. J. Sci. 304, 397–437 (2004).Article ADS CAS Google Scholar 
  53. McKenzie, N. R. et al. Continental arc volcanism as the principal driver of icehouse-greenhouse variability. Science 352, 444–447 (2016).Article ADS CAS Google Scholar 
  54. Glasspool, I. J., Edwards, D. & Axe, L. Charcoal in the Silurian as evidence for the earliest wildfire. Geol32, 381 (2004).Article ADS Google Scholar 
  55. Glasspool, I. J. & Scott, A. C. Phanerozoic concentrations of atmospheric oxygen reconstructed from sedimentary charcoal. Nat. Geosci. 3, 627–630 (2010).Article ADS CAS Google Scholar 
  56. Sønderholm, F. & Bjerrum, C. J. Minimum levels of atmospheric oxygen from fossil tree roots imply new plant−oxygen feedback. Geobiology 19, 250–260 (2021).Article Google Scholar 
  57. Algeo, T. J., Scheckler, S. E. & Scott, A. C. Terrestrial-marine teleconnections in the devonian: links between the evolution of land plants, weathering processes, and marine anoxic events [and discussion]. Philos. Trans.: Biol. Sci. 353, 113–130 (1998).Article Google Scholar 
  58. Montañez, I. P. et al. Climate, pCO2 and terrestrial carbon cycle linkages during late Palaeozoic glacial–interglacial cycles. Nat. Geosci. 9, 824–828 (2016).Article ADS Google Scholar 
  59. Arora, V. K. & Montenegro, A. Small temperature benefits provided by realistic afforestation efforts. Nat. Geosci. 4, 514–518 (2011).Article ADS CAS Google Scholar 
  60. Grossman, E. L. & Joachimski, M. M. Oxygen Isotope Stratigraphy. in Geologic Time Scale 2020 279–307 (Elsevier, 2020). https://doi.org/10.1016/B978-0-12-824360-2.00010-3.
  61. Chen, B. et al. Devonian paleoclimate and its drivers: A reassessment based on a new conodont δ18O record from South China. Earth-Sci. Rev. 222, 103814 (2021).Article CAS Google Scholar 
  62. van Geldern, R. et al. Carbon, oxygen and strontium isotope records of Devonian brachiopod shell calcite. Palaeogeogr., Palaeoclimatol., Palaeoecol. 240, 47–67 (2006).Article Google Scholar 
  63. Liu, Y., Peltier, W. R., Yang, J. & Vettoretti, G. The initiation of Neoproterozoic “snowball” climates in CCSM3: the influence of paleocontinental configuration. Clim. Past 23, 2555–2577 (2013).
  64. Feulner, G. Climate simulations of Neoproterozoic snowball Earth events: Similar critical carbon dioxide levels for the Sturtian and Marinoan glaciations. Earth Planet. Sci. Lett. 404, 200–205 (2014).Article ADS CAS Google Scholar 
  65. Feulner, G. Formation of most of our coal brought Earth close to global glaciation. Proc. Natl. Acad. Sci. USA 114, 11333–11337 (2017).Article ADS CAS Google Scholar 
  66. De Vleeschouwer, D., Crucifix, M., Bounceur, N. & Claeys, P. The impact of astronomical forcing on the Late Devonian greenhouse climate. Glob. Planet. Change 120, 65–80 (2014).Article ADS Google Scholar 
  67. Brugger, J., Hofmann, M., Petri, S. & Feulner, G. On the sensitivity of the devonian climate to continental configuration, vegetation cover, orbital configuration, CO2 concentration, and insolation. Paleoceanogr. Paleoclimatology 34, 1375–1398 (2019).Article ADS Google Scholar 
  68. Isbell, J. L. et al. Glacial paradoxes during the late Paleozoic ice age: Evaluating the equilibrium line altitude as a control on glaciation. Gondwana Res. 22, 1–19 (2012).Article ADS Google Scholar 
  69. Gindre, L., Le Heron, D. & Bjørnseth, H. M. High resolution facies analysis and sequence stratigraphy of the Siluro-Devonian succession of Al Kufrah basin (SE Libya). J. Afr. Earth Sci. 76, 8–26 (2012).Article ADS Google Scholar 
  70. Gerrienne, P., Bergamaschi, S., Pereira, E., Rodrigues, M.-A. C. & Steemans, P. An Early Devonian flora, including Cooksonia, from the Paraná Basin (Brazil). Rev. Palaeobot. Palynol. 116, 19–38 (2001).Article Google Scholar 
  71. Cookson, I. IV—On plant-remains from the Silurian of Victoria, Australia, that extend and connect floras hitherto described. Philos. Trans. R. Soc. Lond. B 225, 127–148 (1935).Article ADS Google Scholar 
  72. Earp, C. Costulatotheca schleigeri (Hyolitha: Orthothecida) from the Walhalla Group (Early Devonian) at Mount Pleasant, central Victoria, Australia. Alcheringa: Australas. J. Palaeontol. 43, 220–227 (2019).Article Google Scholar 
  73. Cohen, K. M., Finney, S. C., Gibbard, P. L. & Fan, J.-X. The ICS international chronostratigraphic chart. Episodes 36, 199–204 (2013).Article Google Scholar 
  74. Polissar, P. J., Fulton, J. M., Junium, C. K., Turich, C. C. & Freeman, K. H. Measurement of 13 C and 15 N Isotopic Composition on Nanomolar Quantities of C and N. Anal. Chem. 81, 755–763 (2009).Article CAS Google Scholar 
  75. Coplen, T. B. et al. New guidelines for δ 13 C measurements. Anal. Chem. 78, 2439–2441 (2006).Article CAS Google Scholar 
  76. Bertrand, P. & Lu-Tao, W. Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS). in Microbeam and Nanobeam Analysis (eds. Benoit, D., Bresse, J.-F., Van’t Dack, L., Werner, H. & Wernisch, J.) 167–182 (Springer Vienna, 1996). https://doi.org/10.1007/978-3-7091-6555-3_8.
  77. Belu, A. M., Graham, D. J. & Castner, D. G. Time-of-flight secondary ion mass spectrometry: techniques and applications for the characterization of biomaterial surfaces. Biomaterials 24, 3635–3653 (2003).Article CAS Google Scholar 
  78. Sodhi, R. N. S. Time-of-flight secondary ion mass spectrometry (TOF-SIMS):—versatility in chemical and imaging surface analysis. Analyst 129, 483–487 (2004).Article ADS CAS Google Scholar 
  79. Pacanowski, R. & Griffies, S. The MOM 3 manual, alpha version. NOAA. Geophysical Fluid Dynamics Laboratory (1999).
  80. Hofmann, M. & M.A., Morales Maqueda. Performance of a second-order moments advection scheme in an Ocean General Circulation Model. J. Geophys. Res. 111, C05006 (2006).ADS Google Scholar 
  81. Fichefet, T. & Morales Maqueda, M. A. Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics. J. Geophys. Res. 102, 12609–12646 (1997).Article ADS Google Scholar 
  82. Petoukhov, V. et al. CLIMBER-2: a climate system model of intermediate complexity. Part I: model description and performance for present climate. Clim. Dyn. 16, 1–17 (2000).Article Google Scholar 
  83. Scotese, C. R. Atlas of Devonian Paleogeographic Maps, PALEOMAP Atlas for ArcGIS, The Late Paleozoic, Maps 65-72Mollweide Projection. (2014).
  84. Belcher, C. M. & McElwain, J. C. Limits for combustion in low O2 redefine paleoatmospheric predictions for the mesozoic. Science 321, 1197–1200 (2008).Article ADS CAS Google Scholar 

Download references


We thank T. Ziegler and F. McSweeney for providing specimens from the palaeobotanical collection of Museum Victoria, O. Seberg and A. A. Pedersen for providing access to lycophyte cultures in the greenhouse of the Botanical Garden, University of Copenhagen, and to F. Sønderholm and C. J. Bjerrum for assistance and use of pCO2 metre. T.W.D. was funded by the Carlsberg Foundation through its Distinguished Associate Professsor program (grant no. CF16–0876) and the Danish Council for Independent Research (grant nos. 7014-00295B, 8102-00005B). C. J. was funded by the National Science Foundation (NSF EAR-1455258). The authors gratefully acknowledge the European Regional Development Fund (ERDF), the German Federal Ministry of Education and Research and the Land Brandenburg for supporting this project by providing resources on the high-performance computer system at the Potsdam Institute for Climate Impact Research. J. B. was funded through the VeWA consortium (Past Warm Periods as Natural Analogues of our high-CO2 Climate Future) by the LOEWE programme of the Hessen Ministry of Higher Education, Research and the Arts, Germany. B.H.L acknowledges funding from the NERC (grant nos NE/R001324/1, NE/T00392/1).

Author information

Authors and Affiliations

  1. Globe institute, University of Copenhagen; Øster Voldgade 5−7, Copenhagen, DenmarkTais W. Dahl & Magnus A. R. Harding
  2. State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan, ChinaTais W. Dahl
  3. Sino-Danish College (SDC), University of Chinese Academy of Sciences, Beijing, ChinaMagnus A. R. Harding
  4. Earth System Analysis, Potsdam Institute for Climate Impact Research, Member of the Leibniz Association, Potsdam, GermanyJulia Brugger & Georg Feulner
  5. Senckenberg Biodiversity and Climate Research Centre, Frankfurt am Main, GermanyJulia Brugger
  6. Center for Integrative Petroleum Research, King Fahd University of Petroleum and Minerals, Dhahran, Saudi ArabiaKion Norrman
  7. School of Biosciences, University of Nottingham, Sutton Bonington Campus, Leicestershire, UKBarry H. Lomax
  8. Department of Earth and Environmental Sciences, Syracuse University; Syracuse, New York, USAChristopher K. Junium


TWD designed the research with BHL’s assistance; C.K.J., M.A.R.H., K.N., and T.W.D collected data; G.F., J.B., B.H.L., and T.W.D. performed research; T.W.D., C.K.J., and G.F. acquired funding. T.W.D. wrote the paper with input from all authors.

Corresponding author

Correspondence to Tais W. Dahl.

Ethics declarations

Competing interests

The authors declare no competing interests

Peer review

Peer review information

Nature Communications thanks Alexander Krause and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary information

Peer Review File

Description of Additional Supplementary Files

Supplementary Dataset 1

Supplementary Dataset 2

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Cite this article

Dahl, T.W., Harding, M.A.R., Brugger, J. et al. Low atmospheric CO2 levels before the rise of forested ecosystems. Nat Commun 13, 7616 (2022). https://doi.org/10.1038/s41467-022-35085-9

Download citation

  • Received18 January 2022
  • Accepted18 November 2022
  • Published20 December 2022
  • DOIhttps://doi.org/10.1038/s41467-022-35085-9
2.2 13 votes
Article Rating
Newest Most Voted
Inline Feedbacks
View all comments
AGW is Not Science
December 22, 2022 6:19 am

The first sentence of the “Introduction” is utter nonsense. There is ZERO EMPIRICAL EVIDENCE that atmospheric CO2 levels have ever had an “effect” on the Earth’s climate.

Once again, hypothetical bullshit represented as supposed “fact.”

Curious George
Reply to  AGW is Not Science
December 22, 2022 11:07 am

Isn’t there a hypothesis that the original composition of the earth’s atmosphere was 80% N2, 20% CO2? Or maybe 75%-25%. At this level an effect of CO2 on the climate seems likely.

December 22, 2022 6:25 am

Models tell us little to nothing other than a sort of confirmation of what is or was thought to be by the scientists/modellers involved.

“At only 500 ppm CO2, however, our Early Devonian palaeoclimate model predicts a clement climate with global mean annual temperatures of 14.8–15.7C for different orbital configurations and a strong latitudinal temperature gradient not too dissimilar from today.”

At only 500ppm? Yet according to the media and the Church of Anthropological Global Warming and Latter Day Climate Change Saints…

Latest CO2 level: 418.2ppm
Safe Level: 350ppm

The Grauniad/Environment

Reply to  strativarius
December 22, 2022 7:52 am

I’m dreaming of a 15C Christmas. It’s very white right now.

(It’s already warmed up to -26C, on its way to -19C at peak warming this afternoon.)

Crispin in Val Quentin but really in Kigali
Reply to  Scissor
December 22, 2022 1:04 pm


We in Val Quentin broke the dawn below -40 and it has struggled up to -27 at 2 PM. The first day of winter was the same. The only advantage of this early Arctic freeze is that we can drive the pickup (if it will start) on Lac Ste Anne. The lake is being populated with ice fishing huts as I write, more each day.

The thin ice warnings have stopped and the government is texting us 3 or 4 times a day warning about the dangerous low temperatures. Some Low records in the region have been broken by as much as 6 deg C. It is warmer in Inuvik, NWT.

Reply to  Scissor
December 22, 2022 1:43 pm

It’s very white right now.

This is what “global warming” really looks like.

Formation of snow reduces OLR from around 280W/m^2 for freezing ground under clear sky to 180W/m^2 so ice forming is retaining up to 100W/m^2. That is a lot of warming.

Most “global warming” is occurring on land north of 40N in winter when the ground temperature is below freezing.

Snowfall is increasing. It is already overtaking melt on Greenland and Iceland, which have increasing permanent ice cover. Greenland has gained 170mm in elevation in the past decade.

Reply to  strativarius
December 22, 2022 1:37 pm

I agree. 500ppm and stable climate is heresy.

This paper starts off with all the political correct acknowledgement that CO2 can somehow change earth’s energy balance then comes and states that the climate was stable at 500ppm. How did it ever get past the gatekeepers?

Frank from NoVA
December 22, 2022 6:29 am

‘Atmospheric CO2 is a greenhouse gas that has affected Earth’s climate throughout geological history.’

I know ‘this’ is a requirement to obtain funding, but it’s demonstrably false based on an honest analysis of the CENOGRID data.

Joseph Zorzin
Reply to  Frank from NoVA
December 22, 2022 9:31 am

It may affect the climate, but only trivially.

Reply to  Joseph Zorzin
December 22, 2022 1:55 pm

The only impact it has is the contribution to total mass and the buoyancy of water vapour in the atmosphere. Negligible is the correct word.

Most land north of 40N is now covered in snow. More than 60% of the solar EMR available at the surface is not thermalised due to the high reflectivity of the snow. Where is the “greenhouse gas” in this energy balance.

And snow coverage extent is increasing year-on-year – where does the changing surface albedo come into the “greenhouse effect”?

Ice in the atmosphere, on water and on land controls earth’s energy balance. Most people living in the NH right now should appreciate that just by taking a look outside.

The tipping point comes when the permanent ice coverage expands. It is a positive feedback that accelerates summer cooling. So far only being observed on Greenland and Iceland.

Javier Vinós
December 22, 2022 7:21 am

Their proposed CO2 levels during the early Phanerozoic are startingly different from those proposed by other methods, like Geocarb III or proxies that show upward of 7000 ppm. It shows how little we know about distant past CO2 levels.

Javier Vinós
Reply to  Javier Vinós
December 22, 2022 7:32 am

comment image

The Early Devonian shows 4000 ppm in Geocarb III or COPSE models, and in Royer’s proxy compilation. This article defends about 500 ppm, more in line with the much-argued Rothman’s estimate.

Reply to  Javier Vinós
December 22, 2022 11:16 am

A fine figure, which used to be on Wikipedia. https://cliscep.com/2021/08/13/down-the-memory-hole/

Frank from NoVA
Reply to  Javier Vinós
December 22, 2022 11:25 am

This Rothman?

‘In 2017, Rothman made a dire prediction: By the end of this century, the planet is likely to reach a critical threshold, based on the rapid rate at which humans are adding carbon dioxide to the atmosphere. When we cross that threshold, we are likely to set in motion a freight train of consequences, potentially culminating in the Earth’s sixth mass extinction.‘

Oh my!

Reply to  Frank from NoVA
December 22, 2022 12:43 pm

According to the Rothman curve (purple) there are no coal measures. But purple is the outlier of the two models, which also match the (important) observations.

I also made an “Ed Hawkins” style graphic, only showing CO2 levels over the last 450,000,000 years:

Scroll down to see the darkmode version, which is better. Or I might paste it in below.

Crispin in Val Quentin but really in Kigali
Reply to  Frank from NoVA
December 22, 2022 1:31 pm

There will be no freight trains a century from now. They run on diesel which will be banned by then. Also by then people will own nothing, have nothing and be happy, or off to the re-education camps. “Vee haff vays of making you happy…”

The article sets up their conclusion about half way through – they are going to show that the low end of their uncertainty is similar to today’s concentration of CO2. That’s fine, though misleading because the central value is one only defensible one.

The off -the-wall claim that it indicates anything about the warming value of CO2 is predicated on the assumption that CO2 controls the temperature. Huh?? If the paleorecord shows us anything it is that a host of things are more important than the CO2 concentration when it comes to temperature.

I understand that the CO2 concentration was very high 1 billion years ago, and something brought it crashing down. What?

The paper concludes that prior to the spread of continental forests, the CO2 dropped to ~500 ppm. Then the massive absorption started, sequestering hundreds of billions of tons of it, all the while the concentration in air remained about the same. Huh. That is quite a claim.

Where did the replacement carbon dioxide gas come from? Further, they hint that as the forests spread and the coal deposits began forming, the CO2 stayed about the same. It is clear they imply that even after hundreds of millions of years during which the whole fossil fuel resourced was laid down, the CO2 concentration, again, stayed roughly the same. How is that possible?

It could have been hiding in that place where it went before the forests started spreading, I suppose.

To claim that, before the coal beds began forming (~360 m.y. ago), the CO2 had already dropped to 500 ppm, one must explain where all the carbon came from to supply the observable reserves of coal and biomass (at least) which is about 43.5 trillion tons, 0.002 of which is in the atmosphere at 400 ppm. How can the spread and sequestration of biomass drop the global concentration from 500 ppm to 400 ppm while creating 43.5 trillion tons of fossil fuels? That doesn’t add up by a factor of 500,

I think even Houdini couldn’t do that.

Javier Vinós
Reply to  Frank from NoVA
December 22, 2022 1:40 pm

Yes, Daniel H. Rothman.
Rothman, D.H., 2002. Atmospheric carbon dioxide levels for the last 500 million yearsProceedings of the National Academy of Sciences99(7), pp.4167-4171.

I’m not saying he is right or wrong about past CO2 levels. I think the lack of proper data prevents us from knowing distant past CO2 levels. Even the data for the past 50 million years is very weak. It is shameful to read Dana Royer that ancient CO2 data is consistent with its role in driving past temperature.

Jim Steele
December 22, 2022 9:10 am

Many of these paleoclimate studies are infested with numerous assumptions that cannot be proven. I get the strong impression researchers are simply trying to “adjust the numbers” to get estimates of CO2 concentrations and temperatures to match, and thus present a unified CO2 drive climate theory. All the divergent numbers weaken the current global warming theory. So create “revised” data estimates to fit the narratives.

As Javier points out estimates of CO2 have an exceedingly wide range. The use of stomata density and size to determine global CO2 are fraught with uncertainties as stomata are influenced by the pressure to allow CO2 in for photosynthesis while controlling transpiration and water out. Microhabitats play a huge role. I can collect leaves from plants within a 100 meter radius and find cottonwood leaves from trees near a stream with abundant stomata on both upper and lower leave surfaces and leaves fro adjacent stonecrop plants with very scare stomata that only open at night.

Furthermore Paleozoic plants were in the midst of an evolutionary development of various vascular systems for better water acquisition and stomatal control.

Jim Steele
Reply to  Jim Steele
December 22, 2022 9:41 am

Their own words reveal their intent, “We suggest from modern palaeoclimate model that the original conjecture that several thousands of ppm CO2 in the atmosphere were necessary to compensate for the ~3% weaker Palaeozoic Sun should be abandoned, in favour of much higher climate sensitivity of CO2 as supported by the new palaeo-CO2 proxy data”

Thus they perpetuate a myth that we are doomed by rising CO2 because after revising data they can argue world is so vulnerable because its sensitivity to CO2 is eeeven higher than first suggested.

Joseph Zorzin
Reply to  Jim Steele
December 22, 2022 9:46 am

but…but… it’s a “modern model” so it must be right! /sarc

Joseph Zorzin
December 22, 2022 9:30 am


Could someone please summarize that in a short paragraph?

Jim Steele
Reply to  Joseph Zorzin
December 22, 2022 9:44 am

To summarize in one sentence: after revising the data and eliminating inconvenient data, they can argue that the world is more vulnerable to global warming because sensitivity to just small increments of CO2 is much higher than ever suggested.”

Their own words reveal their intent, “We suggest from modern palaeoclimate model that the original conjecture that several thousands of ppm CO2 in the atmosphere were necessary to compensate for the ~3% weaker Palaeozoic Sun should be abandoned, in favour of much higher climate sensitivity of CO2 as supported by the new palaeo-CO2 proxy data”

Last edited 3 months ago by Jim Steele
Dean S
Reply to  Jim Steele
December 22, 2022 10:26 am

Which is distinctly at odds with the research showing actual climate sensitivity to CO2 is much lower than anything in the models.

Jim Steele
Reply to  Jim Steele
December 22, 2022 11:25 am

The alarmists are always trying to suggest higher CO2 sensitivity to counteract research finding lower sensitivity such as:

Climate Sensitivity Estimated From Temperature Reconstructions 2 of the Last Glacial Maximum Schmittner (2011)
“combining extensive sea and land surface temperature reconstructions from the Last Glacial Maximum with climate model simulations we estimate a lower median (2.3 K) and reduced uncertainty (1.7–2.6 K 66% probability). Assuming paleoclimatic constraints apply to the future as predicted by our model, these results imply lower probability of imminent extreme climatic change than previously thought.”

Last edited 3 months ago by Jim Steele
Ed Zuiderwijk
December 22, 2022 10:45 am

When climate researchers use the characterisation ‘robust’ then you can be practically sure that it is nothing of the kind.

Peta of Newark
December 22, 2022 4:40 pm

Quote:About this article
Contrived. Long winded, Minutia mining, Magical thinking. Appeals to a computer authority. Politically correct.
And as regards any actual science: Trying to get a glass slipper to fit Herman Munster via a monstrous cause/effect error

and that is what this article is about.
haha:’article‘ has more than one meaning in modern English
so yes, somebody got it right: It is ‘an article

December 22, 2022 6:39 pm

The science is settled
this is just a lot of BS

December 22, 2022 8:29 pm

In the absence of anthropogenic fossil fuel combustion, the dominant atmospheric CO2 source is volcanic outgassing, 

These people really have no idea about the role and movement of CO2 if they think this is how it worked.

Since the IPCC has published the fact in all its COP reports since 2002 that the amount of anthropogenic CO2 in the atmosphere has remained at a steady 3-4%, this claim already looks like someone who is quite clueless about the amount of CO2 in the atmosphere anyway. It may also come as some surprise to them that while CO2 outgassed from oceans or released by microbial respiration is identifiably different to anthropogenic CO2, that which comes from geothermal activity and anthropogenic activity cannot be differentiated.

It would benefit all if scientists looking into climate change (although it really should be physicists as that is the actual field they are investigating) went back to First Principles as their starting point as papers like this start at the wrong point, and just stay wrong.

December 22, 2022 10:45 pm

As soon as I read that CO2 was a green house gas I knew this was nothing but bologna.

%d bloggers like this:
Verified by MonsterInsights