Future climate forcing potentially without precedent in the last 420 million years — TEI-tagged view

Gavin L. FosterpersName, Dana L. RoyerpersName, Daniel J. LuntpersName (2017date). “Future climate forcing potentially without precedent in the last 420 million years.” Nature Communicationstitle·j.

Object of study — atmospheric CO2 and combined CO2-solar radiative forcingterm·key over the last 420 million years

Research question — whether Earth's stable surface temperature reflects stable climate forcing over geological time

Contribution — a CO2 compilation showing CO2 decline cancelled rising solar forcing, placing future forcing beyond precedent

Full text — tagged

The reference list is not part of the text we hold for this paper. Our copy was extracted from the publisher's page and the bibliography did not survive that extraction, so the article below ends before its references. Nothing has been removed by us.

Future climate forcing potentially without precedent in the last 420 million years Abstract The evolution of Earth's climate on geological timescales is largely driven by variations in the magnitude of total solar irradianceterm·dom (TSIterm·dom) and changes in the greenhouse gas content of the atmosphere. Here we show that the slow ∼50 Wm −2 increase in TSIterm·dom over the last ∼420 million years (an increase of ∼9 Wm −2 of radiative forcingterm·key) was almost completely negated by a long-term decline in atmospheric CO 2term·key . This was likely due to the silicate weathering-negative feedbackterm·dom and the expansion of land plants that together ensured Earth's long-term habitability. Humanity's fossil-fuel use, if unabated, risks taking us, by the middle of the twenty-first century, to values of CO 2 not seen since the early Eocene (50 million years ago). If CO 2 continues to rise further into the twenty-third century, then the associated large increase in radiative forcingterm·key, and how the Earth system would respond, would likely be without geological precedent in the last half a billion years.

The primary source of energy in the Earth's climate system is incoming solar radiation, termed total solar irradianceterm·dom, or TSIterm·dom. At equilibrium, because of the requirement of energy conservation, the Earth's radiative budget must balance such that TSIterm·dom is equal to outgoing longwave radiation at the top of the atmosphere (a state known as radiative equilibrium). Since the radiation emitted by a body is a function of surface temperature, the Earth's ‘effective temperature' ( T E ) is the temperature at which radiative equilibrium is achieved assuming the Earth acts like a blackbody, and can be calculated (in K) using the following expression: where F s is the TSIterm·dom (currently ∼1,368 Wm −2 ), σ is the Stefan–Boltzmann constant (5.67 × 10 −8 W m −2 K −4 ), A =0.29 (ref. 1 ) is the Earth's average planetary albedo (the fraction of incoming radiation scattered or reflected back out of the atmosphere by clouds, particulates in the atmosphere and the Earth's surface), and the factor of 4 accounts for the spherical and rotating nature of the Earth. The 31 K difference between T E and the observed surface temperature of the Earth (+14.0 °C or 287.1 K is the 1961–1990 mean 2 ) is almost entirely due to the action of the greenhouse effectterm·key (for example, ref. 3 ).

The majority (∼75%) of the greenhouse effectterm·key is due to the warming effects of water vapour and clouds, with the non-condensing greenhouse gasses (predominantly CO 2 and CH 4 ) accounting for the remaining 25% (ref. 3 ). However, at the temperatures and pressures typical of the Earth's surface, water vapour and clouds act as feedbacksterm·dom rather than driversterm·dom of the greenhouse effectterm·key, with CO 2 and CH 4 , and the other non-condensing GHGs (for example, N 2 O) determining the overall strength of the greenhouse effectterm·key 3 . Given this understanding, and that summarized in equation (1) , the climatic evolution of the Earth over geological time is largely a function of the concentration of the non-condensing greenhouse gases, planetary albedo (A) and the TSIterm·dom ( F s ; for example, ref. 4 ).

Owing to the way the Sun generates energy though the nuclear fusion of hydrogen into helium, over time its luminosity has increased 5 . This is a relatively well-understood process and the TSIterm·dom at any time can be approximated by the following function 5 : where F s is the present day TSIterm·dom (1,368 Wm −2 ), t is the time of interest since the formation of the Earth (Myrs) and t 0 is the age of the Earth (4,567 Myrs (ref. 6 ), although here following ref. 5 we use a t 0 =4,700 Myrs). As a result, there has been an increase in TSIterm·dom of ∼400 Wm −2 since the formation of the Earth. The following equation relates this change in TSIterm·dom relative to today to radiative forcingterm·key (ΔF sol ) assuming a constant planetary albedo ( A ) through time: The global mean surface temperature change (Δ T , in K) for a given change in radiative forcingterm·key (Δ F , in Wm −2 ) can be described by: Where S is the sensitivity parameter 7 , in K W −1 m 2 . Following ref. 4 , and assuming an effective emissivity, ɛ , of the Earth of 0.6 (that is, 40% of longwave radiation is absorbed by greenhouse gases in the atmosphere), equation (1) can be used to define the following equation that describes S for today's climate (where the surface temperature=287.1 K) in the absence of any climate feedbacksterm·dom (for example, water vapour, sea-ice etc.), also known as the Planck response (hence the subscript ‘P') 8 : S P depends on the overall strength of the greenhouse effectterm·key and the surface temperature of the Earth, and so this value is not necessarily applicable throughout Earth's history; nonetheless, this treatment provides a first-order constraint that the long-term secular increase in TSIterm·dom of 400 Wm −2 would be associated with a secular warming of at least ∼20 K over the last 4.5 billion years, if albedo and emissivity remained constant. Since the combined effect of the greenhouse gas and other climate feedbacksterm·dom (for example, water vapour, lapse rate, sea-ice etc.) is positive 8 , this is a minimum estimate of S and when all climate feedbacksterm·dom are considered, S a (where ‘a' denotes actuo after ref. 7 ) is likely in the range of 0.8– 1.6 K W −1 m 2 (for example, refs 7 , 9 , 10 , 11 ).

The actual evolution of Earth's temperature through geological time is a subject of considerable debate (for example, ref. 12 ), yet there is a longstanding view that, despite the increase in solar output, Earth's surface temperature was, for much of geological time, warmer, not colder, than today (for example, refs 13 , 14 ). This apparent inconsistency is known as the ‘Faint Young Sunterm·dom' paradox and to reconcile the observation of a relatively stable climate in the face of increasing solar output through time requires a parallel change in some other factors that influence Earth's radiative budget. First and foremost among these factors is thought to be a concomitant reduction in the strength of the greenhouse effectterm·key, and in atmospheric CO 2term·key concentration in particular (for example, refs 15 , 16 ). The radiative forcingterm·key from changing CO 2 can be estimated using the following formulation 17 : where Δ F CO 2 is the radiative forcingterm·key from CO 2 change (in Wm −2 ), C represents the concentration of CO 2 at the time of interest and C 0 is the pre-industrial concentration of CO 2 =278 p.p.m. For the decrease in atmospheric CO 2term·key to balance the increase in solar output over time, the following equation is required to be true (or at least approximately so): In this formulation we are ignoring the other non-condensing greenhouse gases (for example, CH 4 , N 2 O). This is partly out of practical necessity since there are currently no proxiesterm·key for these gases, but is justifiable to a first order because variations in CO 2 accounted for ∼80% of the greenhouse gas forcing on glacial–interglacial timescales 4 and CO 2 was likely the dominant greenhouse gas for the past 500 Myrs (refs 18 , 19 , 20 ).

Testing whether the apparent stability of surface temperature is a consequence of the relative stability of climate forcing, therefore, requires high-density, high-quality data on past atmospheric CO 2term·key levels. Consequently, in this study we compile the available CO 2 data for the last 420 Myrs where, although the magnitude in the change of solar output is reduced relative to changes on longer timescales, data availability is sufficient to provide a detailed picture of the evolution of atmospheric CO 2term·key . The relative stability of climate over this interval is constrained by three lines of evidence: (i) the continued presence of complex life, given the thermal limits of metazoans (<40–50 °C; ref. 21 ) and despite a number of mass extinctions 22 ; (ii) the continuous presence of liquid water and the occurrence of climate-sensitive sedimentary deposits (for example, coal versus glacial deposits) that indicate warm and cold periods but no evidence of an overall long-term secular trend ( Fig. 1 ) 23 ; and (iii) moderate tropical sea surface temperatures (25–40 °C) as recorded by oxygen and clumped isotope analysis of fossil marine carbonates 24 25 26 27 . Our new compilation shows that this climate stability was the result of a long-term decline in atmospheric CO 2term·key that, in terms of radiative forcingterm·key, approximately cancelled out the increase in solar output. This long-term view of climate forcing provides a valuable geological context for potential levels of CO 2 in our warming future.

Results A new CO 2 compilation In order to better understand the role of the greenhouse effectterm·key in Earth's climate evolution, we have compiled ∼1,500 discrete estimates of atmospheric CO 2term·key from five independent techniques drawn from 112 published studies covering the last 420 Myrs ( Supplementary Data 1 ). We uniformly apply a set of criteria and the latest understanding, described in the supplement, to screen the available CO 2 records. By following these criteria, the ages and CO 2 values associated with some records are revised (see Methods), and ∼1/5 of the published estimates have been excluded (leaving n =1,241 in our final compilation). This standardization process helps ensure the highest-quality compilation whose individual records can be more cleanly compared to one another.

Our new compilation is shown in Fig. 1 and, compared to older compilations, there is better agreement between the different methods of CO 2 reconstruction 28 . In this case, this is largely due to the refinement of the pedogenic carbonate proxy following ref. 29 (see Methods). However, despite our efforts to improve the compilation, there are still relatively large uncertainties at times and disagreements remain between techniques for some time intervals ( Fig. 1 ). It is also apparent that the early parts of the record rely on fewer observations and are characterized by a reduced diversity in proxy type ( Supplementary Fig. 1 ). In order to gain a better appreciation of the multi-million-year evolution of CO 2 , and in light of these uncertainties and limitations we have followed a probabilistic approach where Monte Carlo resampling is used to generate 1,000 artificial time series of CO 2 with each data point randomly perturbed within its age ( X ) and CO 2 ( Y ) uncertainty. Each realization was then interpolated to a regular 0.5 Myr spacing and a LOESSterm·key fit was performed with the degree of smoothing optimized by generalized cross-validation 30 . At each time step the distribution of LOESSterm·key fits was evaluated and the maximum probability, hence the most likely value for long-term CO 2 , and the associated 68 and 95 percentile ranges were determined ( Fig. 1 and Supplementary Data 2 ).

Our new compilation and probabilistic treatment reveals that for the last ∼420 Myrs, CO 2 , on the whole, has been elevated compared to pre-industrial values (278 p.p.m.). The highest CO 2 values of ∼2,000 p.p.m. were reached during the Devonian (∼400 Myrs ago) and Triassic (220–200 Myrs ago), with individual estimates ranging up to a maximum of ∼3,700±1,600 p.p.m. at 215 Myrs. In contrast, values close to pre-industrial are found during much of the Carboniferous (∼300 Myrs ago) and late Cretaceous (∼80 Myrs ago). A linear fit to either the entire CO 2 compilation, or a resampling of our LOESSterm·key fit to reflect the original data density, reveals that long-term average CO 2 has declined over the last 420 Myrs by ∼3.4 p.p.m. per Myrs ( Fig. 1 ).

Radiative forcingterm·key through the last 420 million years Radiative forcingterm·key from CO 2 change can be calculated using equation (6) . This can then be combined with the radiative forcingterm·key from the increase in solar output defined by equation (3) ( Fig. 2a ); the temporal evolution of the simple sum of these terms (Δ F CO2,sol ) is shown in Fig. 2b . As has been noted elsewhere with prior Phanerozoicterm·key CO 2 compilations (for example, refs 31 , 32 ), there is a good first-order agreement between CO 2 levels and the occurrence of greenhouse/icehouse states ( Fig. 1 ): CO 2 is high during greenhouse climate states and low during icehouse climates. Our data support previous work suggesting that, when Δ F CO2,sol is considered, icehouse states generally occur when climate forcing drops below approximately +1–2 Wm −2 (refs 31 , 32 ). This treatment also reveals that for the last 420 Myrs, despite a gradual 4% increase in solar luminosity (equivalent to +9 Wm −2 of radiative forcingterm·key), there has been very little long-term change in Δ F CO2,sol (despite shorter-term fluctuations of up to 10 Wm −2 ), with a linear fit to all the data suggesting only a slight decrease of 0.004 Wm −2 per Myrs (±0.001 1 s.e.m.; P =0.0003; red line on Fig. 2b ); a linear fit to our LOESSterm·key fit resampled with the original data density yields a slightly higher negative slope (−0.008±0.001 1 s.e.m.; P <0.0001; black line Fig. 2b ). This indicates that the long-term decrease in CO 2 over the last 420 Myrs has largely compensated for the increase in solar output over the interval. In terms of shorter-term radiative imbalances, 95% of the Δ F CO2,sol values are within ±7 Wm −2 and 68% within +5/−3 Wm −2 ( Fig. 3 ). Thus, even on multi-million-year timescales (and sometimes shorter) the external radiative forcingterm·key (not considering the action of additional feedbacksterm·dom internal to the climate system such as the ice-albedo feedbackterm·dom) has not varied beyond ±2.0% of the total incoming radiation (∼340 Wm −2 ).

Discussion The concentration of atmospheric CO 2term·key on multi-million-year timescales depends largely on the balance between the input of CO 2 from volcanism, metamorphism and organic carbon weathering and the output of CO 2 from silicate weatheringterm·dom and organic carbon burial 15 . Since the magnitude of silicate weatheringterm·dom is climatically sensitive (that is, a function of temperature and precipitation; for example, ref. 33 ), silicate weatheringterm·dom also buffers the CO 2 content of the atmosphere by acting as a negative feedbackterm·dom 15 . However, we note that in many modern settings silicate weatheringterm·dom is supply-limited 33 . Thus, for silicate weatheringterm·dom to operate as a climate feedbackterm·dom, the Earth must be sufficiently tectonically active to ensure an adequate supply of fresh minerals 34 .

The marked ‘double-hump' pattern of the CO 2 reconstruction ( Fig. 1 ) that is also common to other compilations (for example, refs 28 , 35 ) is likely caused by changes in the inputs and outputs of CO 2 in response to the supercontinent cycle 36 . For example, the low CO 2 during the Carboniferous (∼300 Myrs ago) and during the later parts of the Cenozoic (last 65 Myrs) was likely a result of a reduced volcanic flux and/or enhanced silicate weatheringterm·dom, at least in part due to higher continental relief during supercontinent construction phase. During the intervening greenhouse intervals the reverse was likely true (high volcanic flux and/or low silicate weatheringterm·dom due to low relief) 36 37 . An additional factor in the decline in CO 2 was the overall enhancement of silicate weatheringterm·dom due to the expansion of the terrestrial biosphereterm·dom over the last 400 million years or so (refs 16 , 38 ). It is probably this expansion, along with the silicate weathering-negative feedbackterm·dom (facilitated by sufficiently active plate tectonic regime 34 ), that was the key in keeping Δ F CO2,Sol relatively constant over the long term (for example, ref. 16 ). Indeed, as with previous compilations 35 , our new CO 2 record largely agrees with the GEOCARB carbon cycle model that incorporates these long-term changes in silicate weatheringterm·dom ( Supplementary Fig. 2 ), especially when the uncertainties in both records are considered 35 . This serves to further underscore the importance of the rise of the terrestrial biosphere, along with plate tectonics and silicate weatheringterm·dom, in ensuring climate stability through time 16 and provides additional support for our broad understanding of the Earth's long-term carbon cycle as encapsulated by the GEOCARB model 35 38 ( Supplementary Fig. 2 ).

Regardless of the ultimate cause for the observed relative stability in Δ F CO2,sol over the last 420 million years, business-as-usual emission scenarios (for example, representative concentration pathway RCP8.5term·key) 39 for fossil fuel emissions suggest that atmospheric CO 2term·key could peak in 2,250 AD at ∼2,000 p.p.m. CO 2 values as high as this were last seen in the Triassic around 220–200 Myrs ago ( Figs 3 and 4 ). However, because of the steady increase in solar output over time, in terms of radiative forcingterm·key by the end of this century RCP8.5term·key is similar to the early Eocene, and by 2,250 AD exceeds what is recorded in the geological record for at least 99.9% of the last 420 Myrs ( Figs 3 and 4 ). A recent study suggested that if both conventional and non-conventional fossil fuel reserves (amounting to ∼12,000 Pg C; ref. 40 ) were exhausted in such a business-as-usual scenario, atmospheric CO 2term·key could rise to ∼5,000 p.p.m. by 2,400 AD (refs 41 , 42 ), which is clearly higher, in terms of both forcing and absolute CO 2 , than at any time captured by our compilation ( Figs 3 and 4 , Wink12K scenario). Such a scenario therefore risks subjecting the Earth to a climate forcing that has no apparent geological precedent, for at least the last 420 Myrs. We should be aware of course of the limitations of the geological record, and it is debatable whether an extreme climate change event analogous to the Anthropocene, if it existed at all, would leave a detectable signal, given our current CO 2 proxiesterm·key and records 43 . Nonetheless, prolonged warm greenhouse climate states have occurred in the past, but the rates of climate change in the geological record are on the whole very likely slower than what we are currently experiencing (for example, ref. 44 ). Unabated fossil fuel use therefore has the potential to push the climate system into a state that has not been seen on Earth in at least the last 420 Myrs.

Methods Assembling the CO 2 compilation The new CO 2 compilation here consists of 1,241 independent estimates coming from five proxy methods and 112 published studies. In putting this compilation together, we applied a uniform set of criteria described next to screen potential CO 2 records. By following these criteria, the CO 2 values associated with some records are revised, while other records are excluded all together. This standardization process helps ensure a high-quality compilation whose individual records can be more cleanly compared to one another. Of course, some of our criteria may be judged at some point to be incorrect. It is in this light, and in the spirit of full transparency, that we describe our criteria. We consider our compilation to be a ‘living document', not only because in the future new records will be added but because existing records will be interpreted in new ways. The compilation can be found in Supplementary Data 1 .

First, we excluded all goethite-based CO 2 estimates 45 46 47 48 due to uncertainties in modelling some of the isotopic fractionation factors 49 . We also excluded all B/Ca-based estimates 50 because the environmental controls on B/Ca are currently not well understood 51 . The nahcolite proxy is based on well-described mineral phase equilibria 52 , but we exclude the published estimates because it is not possible to assign mean values (all values within a defined CO 2 interval are equally probable). It is important to note, however, that the range described by these nahcolite estimates agree well with the other data in our compilation. We do not include the boron-based estimates of Pearson and Palmer 53 owing to problems related to potential diagenesis, potential analytical issues, vital effects of extinct species and the evolution of seawater δ 11 B and alkalinity (as discussed in ref. 19 ). For the alkenone-based records of Pagani and colleagues 54 55 56 57 58 , we adopt the compilation of ref. 58 who imposed a uniform set of quality controls on the original data sets and used a TEX 86 temperature record (versus the δ 18 O of ref. 54 ). We update all Cenozoic liverwort-based estimates 59 with the atmospheric δ 13 C model of ref. 60 .

For estimates calculated from stomatal ratios (SR, where SR=ratio of stomatal frequency in nearest living equivalent to fossil), two transfer functions to CO 2 have been proposed: 1 SR=1 RCO 2 and 1 SR=2 RCO 2 , where RCO 2 =ratio of atmospheric CO 2term·key in the past to pre-industrial conditions (300 p.p.m.; refs 61 , 62 ). Following Beerling and RoyerpersName 63 , we use these two functions to establish the uncertainty bounds in estimated CO 2 . Owing to the ad hoc nature of these functions, CO 2 estimates from stomatal ratios should be considered semiquantitative only. We modify many of the stomatal-based estimates of Retallack 64 . First, we only use estimates associated with five or more cuticle fragments, as this is the minimum sampling required in most cases for robust estimates 65 . We use the quantitative transfer function of RoyerpersName 65 for estimates coming from Ginkgo adiantoides , the fossil species considered conspecific with extant G. biloba 66 . For estimates coming from other species within the Ginkgo genus we use the stomatal ratio approach, and we exclude estimates not coming from Ginkgo as these other groups are too distantly related to Ginkgo and have not been calibrated correctly 67 .

A recent development that we have incorporated is the lowering of CO 2 estimates from the pedogenic carbonate method 29 . This is because one of the key input parameters, the concentration of biologically derived CO 2 in the soil ( S (z)), had been overestimated in most paleosols by a factor of two or more 29 68 . Here we adjust most CO 2 estimates assuming a S (z) of 2,000 p.p.m. (ref. 68 ). This is a simple correction that masks variability in S (z) 69 , but it is nonetheless a useful place to start and in most cases represents an improvement over previous estimates 70 . This is an area where we anticipate substantial revision in the near future. Along these lines, several proxiesterm·key for S (z) have been proposed 69 70 71 72 73 74 75 76 , and we use estimates of S (z) derived from these methods when available. We also replace some of the pedogenic carbonate-based estimates of ref. 73 with those based on the same sediments 74 75 because these newer studies provide better constraints on organic matter δ 13 C. Lastly, we collapse the high-resolution stomatal-based record of ref. 76 into one estimate calculated from the mean of the individual estimates.

Uncertainties in paleo-CO 2 estimates have been assessed in different ways. Most earlier uncertainties associated with the pedogenic carbonate method, for example, only incorporate uncertainty in S (z); potential uncertainties in other input variables are ignored. In addition, errors associated with stomatal-based methods that involve quantitative regression equations usually only reflect uncertainty in the regressions; uncertainties in the fossil measurements are ignored. This means that error bars can convey different meanings, which hinders comparisons across methods. In recent years, Monte Carlo simulations for propagating uncertainty in multiple input variables have been developed for all five proxy systems 11 55 59 77 . So long as the reported percentile ranges are identical, these simulations facilitate better cross-method comparisons.

We attempt here to standardize reported uncertainties so that they incorporate uncertainty in all input variables; we also standardize the percentile levels (16 and 84; equivalent to ±1 s.d. for a normally distributed population). For uncertainties derived from Monte Carlo simulations that propagate uncertainties in most or all input variables (this covers most estimates from the alkenone, boron and liverwort proxiesterm·key, as well as most recent estimates from pedogenic carbonates and stomata), we adjust if needed to the 16 and 84 percentile levels: for example, if the 2.5 and 97.5 percentiles are reported, we assume the population has a normal distribution and cut the uncertainties in half. Because most paleo-CO 2 estimates have a right-skewed distribution, this means that we are overestimating the magnitude of the upper uncertainty and underestimating the magnitude of the lower uncertainty. For uncertainties with the pedogenic carbonate method derived by varying only S (z), we adopt for the 16 and 84 percentiles the generic recommendation of +100%/−50% of the median value calculated from multifactorial Monte Carlo simulations 78 ; these uncertainties are larger and probably more representative than the reported uncertainties. Similarly, for stomatal estimates whose uncertainties only take into account variance in the transfer function, we assume 16 and 84 percentiles equal to +100%/−40% of the median value, a range that takes into account representative uncertainty in both the transfer function and the fossil measurements 77 .

All ages have been updated to the 2012 Geologic Time Scale 79 . Many CO 2 estimates from terrestrial proxiesterm·key have biostratigraphic (for example, ‘Maastrichtian'), magnetostratigraphic or chronostratigraphic constraints. For records with biostratigraphic or magnetostratigraphic constraints, we update their ages following the new Time Scale. Many CO 2 estimates have no reported age uncertainties. In almost all cases, these estimates have excellent age control, the most common example being estimates from marine sediment cores (for example, all boron- and alkenone-based estimates). Some terrestrial-based data fall under this category too, for example, the Triassic–Jurassic sediments in the Newark Supergroup 80 . When no age uncertainty is reported we take the 1 s.d. to be 0.001 Myr for marine records and 4% for terrestrial records.

Monte Carlo sampling and LOESSterm·key fitting Monte Carlo approach was used to generate 1,000 realizations of the multiproxy CO 2 time series with each data point randomly perturbed within its age and CO 2 uncertainty as given in Supplementary Data 1 . For CO 2 estimates with asymmetric uncertainties, we force symmetry by assuming ±1 s.d. is equal to the average offset between the median and the 16 and 84 percentiles. Each realization of the time series was then interpolated to a regular 0.5 Myr spacing and LOESSterm·key curves were fitted to each interpolated realization using the programme R 81 . The optimal degree of smoothing was determined for each interpolated realization using generalized cross-validation 30 . At each 0.5 Myr time step the distribution of LOESSterm·key curves was evaluated and the probability distribution was determined. From this probability distribution the most likely value (or probability maximum) and the upper and lower limits corresponding to 68 and 95% confidence limits were identified. These variables are given in Supplementary Data 2 .

Owing to the nature of the proxiesterm·key used the make-up and coverage of our new CO 2 compilation varies considerably as a function of time. By grouping the data into 2.5 million year bins, Supplementary Fig. 1 shows that the number of observations in each 2.5 million year bin generally increases towards the present, with a notable exception of a peak at ∼200 million years ago at the Triassic–Jurassic boundary. Supplementary Fig. 1 also shows that the boron isotope and alkenone δ 13 C methods are restricted to the Cenozoic and liverwort δ 13 C proxy are concentrated in the middle part of the last 420 million years. Estimates based on fossil plant stomata appear focused in the last 200 million years, with only sporadic coverage over the early part of the interval. Only CO 2 estimates from pedogenic carbonate δ 13 C method are relatively evenly spread across the last 420 million years. To examine the influence of this uneven distribution of proxy types we performed the LOESSterm·key fitting a further two times but without one of the dominant proxy systems (no-pedogenic carbonate, no-stomata), these alternatives are shown with the complete data set in Supplementary Fig. 3 . While subtle variations exist, overall the large-scale structure remains intact, although if the pedogenic carbonate estimates are discarded the early part of the Phanerozoicterm·key is largely under sampled. In terms of long-term CO 2 change, however, discarding all the estimates from either the pedogenic carbonate method or the stomatal method does not change our principal conclusions. For instance, the long-term CO 2 decline over the last 420 million years without the pedogenic carbonates is 2 p.p.m. per million years, and without the estimates from the plant stomata method it is 4 p.p.m. per million years. Both of which are similar to the CO 2 decline we determine for the complete data set (3.4 p.p.m. per million years), indicating that a long-term decline in CO 2 over the last 420 million years is a robust observation. In terms of Δ F CO2,sol the long-term trends are also similar being 0.008±0.001 Wm −2 ( R 2 =0.03 and P <0.001) and 0.0001±0.001 Wm −2 ( R 2 =<0.001, P =0.94) for no-stomata and no-pedogenic carbonate, respectively.

Data availability The authors declare that all data supporting the findings of this study are available within the Supplementary Information files associated with this manuscript.

Additional information How to cite this article: FosterpersName, G. L. et al . Future climate forcing potentially without precedent in the last 420 million years. Nat. Commun. 8, 14845 doi: 10.1038/ncomms14845 (2017).

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

Supplementary Material Supplementary Figures Phanerozoicterm·key CO2 compilation used in analyses LOESSterm·key fit to the CO2 data set in Sup. Data 1.