
Studies of Carbon 14 in the atmosphere emitted by nuclear tests indicate that the Bern model used by the IPCC is inconsistent with virtually all reported experimental results.
Guest essay by Gösta Pettersson
The Keeling curve establishes that the atmospheric carbon dioxide level has shown a steady long-term increase since 1958. Proponents of the antropogenic global warming (AGW) hypothesis have attributed the increasing carbon dioxide level to human activities such as combustion of fossil fuels and land-use changes. Opponents of the AGW hypothesis have argued that this would require that the turnover time for atmospheric carbon dioxide is about 100 years, which is inconsistent with a multitude of experimental studies indicating that the turnover time is of the order of 10 years.
Since its constitution in 1988, the United Nation’s Intergovernmental Panel on Climate Change (IPCC) has disregarded the empirically determined turnover times, claiming that they lack bearing on the rate at which anthropogenic carbon dioxide emissions are removed from the atmosphere. Instead, the fourth IPCC assessment report argues that the removal of carbon dioxide emissions is adequately described by the ‘Bern model‘, a carbon cycle model designed by prominent climatologists at the Bern University. The Bern model is based on the presumption that the increasing levels of atmospheric carbon dioxide derive exclusively from anthropogenic emissions. Tuned to fit the Keeling curve, the model prescribes that the relaxation of an emission pulse of carbon dioxide is multiphasic with slow components reflecting slow transfer of carbon dioxide from the oceanic surface to the deep-sea regions. The problem is that empirical observations tell us an entirely different story.
The nuclear weapon tests in the early 1960s have initiated a scientifically ideal tracer experiment describing the kinetics of removal of an excess of airborne carbon dioxide. When the atmospheric bomb tests ceased in 1963, they had raised the air level of C14-carbon dioxide to almost twice its original background value. The relaxation of this pulse of excess C14-carbon dioxide has now been monitored for fifty years. Representative results providing direct experimental records of more than 95% of the relaxation process are shown in Fig.1.
Figure 1. Relaxation of the excess of airborne C14-carbon dioxide produced by atmospheric tests of nuclear weapons before the tests ceased in 1963
The IPCC has disregarded the bombtest data in Fig. 1 (which refer to the C14/C12 ratio), arguing that “an atmospheric perturbation in the isotopic ratio disappears much faster than the perturbation in the number of C14 atoms”. That argument cannot be followed and certainly is incorrect. Fig. 2 shows the data in Fig. 1 after rescaling and correction for the minor dilution effects caused by the increased atmospheric concentration of C12-carbon dioxide during the examined period of time.
Figure 2. The bombtest curve. Experimentally observed relaxation of C14-carbon dioxide (black) compared with model descriptions of the process.
The resulting series of experimental points (black data i Fig. 2) describes the disappearance of “the perturbation in the number of C14 atoms”, is almost indistinguishable from the data in Fig. 1, and will be referred to as the ‘bombtest curve’.
To draw attention to the bombtest curve and its important implications, I have made public a trilogy of strict reaction kinetic analyses addressing the controversial views expressed on the interpretation of the Keeling curve by proponents and opponents of the AGW hypothesis.
(Note: links to all three papers are below also)
Paper 1 in the trilogy clarifies that
a. The bombtest curve provides an empirical record of more than 95% of the relaxation of airborne C14-carbon dioxide. Since kinetic carbon isotope effects are small, the bombtest curve can be taken to be representative for the relaxation of emission pulses of carbon dioxide in general.
b. The relaxation process conforms to a monoexponential relationship (red curve in Fig. 2) and hence can be described in terms of a single relaxation time (turnover time). There is no kinetically valid reason to disregard reported experimental estimates (5–14 years) of this relaxation time.
c. The exponential character of the relaxation implies that the rate of removal of C14 has been proportional to the amount of C14. This means that the observed 95% of the relaxation process have been governed by the atmospheric concentration of C14-carbon dioxide according to the law of mass action, without any detectable contributions from slow oceanic events.
d. The Bern model prescriptions (blue curve in Fig. 2) are inconsistent with the observations that have been made, and gravely underestimate both the rate and the extent of removal of anthropogenic carbon dioxide emissions. On basis of the Bern model predictions, the IPCC states that it takes a few hundreds of years before the first 80% of anthropogenic carbon dioxide emissions are removed from the air. The bombtest curve shows that it takes less than 25 years.
Paper 2 in the trilogy uses the kinetic relationships derived from the bombtest curve to calculate how much the atmospheric carbon dioxide level has been affected by emissions of anthropogenic carbon dioxide since 1850. The results show that only half of the Keeling curve’s longterm trend towards increased carbon dioxide levels originates from anthropogenic emissions.
The Bern model and other carbon cycle models tuned to fit the Keeling curve are routinely used by climate modellers to obtain input estimates of future carbon dioxide levels for postulated emissions scenarios. Paper 2 shows that estimates thus obtained exaggerate man-made contributions to future carbon dioxide levels (and consequent global temperatures) by factors of 3–14 for representative emission scenarios and time periods extending to year 2100 or longer. For empirically supported parameter values, the climate model projections actually provide evidence that global warming due to emissions of fossil carbon dioxide will remain within acceptable limits.
Paper 3 in the trilogy draws attention to the fact that hot water holds less dissolved carbon dioxide than cold water. This means that global warming during the 2000th century by necessity has led to a thermal out-gassing of carbon dioxide from the hydrosphere. Using a kinetic air-ocean model, the strength of this thermal effect can be estimated by analysis of the temperature dependence of the multiannual fluctuations of the Keeling curve and be described in terms of the activation energy for the out-gassing process.
For the empirically estimated parameter values obtained according to Paper 1 and Paper 3, the model shows that thermal out-gassing and anthropogenic emissions have provided approximately equal contributions to the increasing carbon dioxide levels over the examined period 1850–2010. During the last two decades, contributions from thermal out-gassing have been almost 40% larger than those from anthropogenic emissions. This is illustrated by the model data in Fig. 3, which also indicate that the Keeling curve can be quantitatively accounted for in terms of the combined effects of thermal out-gassing and anthropogenic emissions.
Figure 3. Variation of the atmospheric carbon dioxide level, as indicated by empirical data (green) and by the model described in Paper 3 (red). Blue and black curves show the contributions provided by thermal out-gassing and emissions, respectively.
The results in Fig. 3 call for a drastic revision of the carbon cycle budget presented by the IPCC. In particular, the extensively discussed ‘missing sink’ (called ‘residual terrestrial sink´ in the fourth IPCC report) can be identified as the hydrosphere; the amount of emissions taken up by the oceans has been gravely underestimated by the IPCC due to neglect of thermal out-gassing. Furthermore, the strength of the thermal out-gassing effect places climate modellers in the delicate situation that they have to know what the future temperatures will be before they can predict them by consideration of the greenhouse effect caused by future carbon dioxide levels.
By supporting the Bern model and similar carbon cycle models, the IPCC and climate modellers have taken the stand that the Keeling curve can be presumed to reflect only anthropogenic carbon dioxide emissions. The results in Paper 1–3 show that this presumption is inconsistent with virtually all reported experimental results that have a direct bearing on the relaxation kinetics of atmospheric carbon dioxide. As long as climate modellers continue to disregard the available empirical information on thermal out-gassing and on the relaxation kinetics of airborne carbon dioxide, their model predictions will remain too biased to provide any inferences of significant scientific or political interest.
References:
Climate Change 2007: IPCC Working Group I: The Physical Science Basis section 10.4 – Changes Associated with Biogeochemical Feedbacks and Ocean Acidification
http://www.ipcc.ch/publications_and_data/ar4/wg1/en/ch10s10-4.html
Climate Change 2007: IPCC Working Group I: The Physical Science Basis section 2.10.2 Direct Global Warming Potentials
http://www.ipcc.ch/publications_and_data/ar4/wg1/en/ch2s2-10-2.html
GLOBAL BIOGEOCHEMICAL CYCLES, VOL. 15, NO. 4, PAGES 891–907, DECEMBER 2001 Joos et al. Global warming feedbacks on terrestrial carbon uptake under the Intergovernmental Panel on Climate Change (IPCC) emission scenarios
ftp://ftp.elet.polimi.it/users/Giorgio.Guariso/papers/joos01gbc[1]-1.pdf
Click below for a free download of the three papers referenced in the essay as PDF files.
Paper 1 Relaxation kinetics of atmospheric carbon dioxide
Paper 2 Anthropogenic contributions to the atmospheric content of carbon dioxide during the industrial era
Paper 3 Temperature effects on the atmospheric carbon dioxide level
================================================================
Gösta Pettersson is a retired professor in biochemistry at the University of Lund (Sweden) and a previous editor of the European Journal of Biochemistry as an expert on reaction kinetics and mathematical modelling. My scientific reasearch has focused on the fixation of carbon dioxide by plants, which has made me familiar with the carbon cycle research carried out by climatologists and others.
Greg Goodman says:
July 2, 2013 at 9:46 am
Your example of the hottest water is not the one that is most sensitive to change. Since the North pole is an ocean with large expanses of exposed water for a large part of the year, it seems to be more important than the continental south pole
The reaction to temperature changes at the sink places is similar as at the release places, with as difference that in general the temperature change is not that much (the main sinks are at the edge of the ice), but the changes are in the exact place where the sinks are and the area involved. I have no idea what drives the connection between AO and the CO2 rate of change. Probably the AO drives ocean temperature changes and that drives the CO2 rate of change variability… Similar connections can be made between ENSO and ocean temperatures and CO2 rate of change…
Re: jkanders
I think you have misunderstood the Bern Formula TerryS. Where have you got the idea that the model assumes that the CO2 does not mix?
I got the idea from the Berne formula. The Berne formula represents a system whereby CO2 does not mix. It doesn’t matter what the designers of the Berne model specified (or the IPCC) what matters is what the model actually ends up representing. It represents a system without CO2 mixing. Obviously the IPCC did not specify this but that is what they got.
The other error I think you make is that you mix your formula for residence time based on the ratio between CO2 content in the atmosphere and the annual flow of CO2 in and out of the atmosphere.
No, I don’t. See AR4 2.10.2 footnote 1. There they present a simplified version of the formula with 4 terms instead of 6. They describe the formula as “The decay of a pulse of CO2 with time t”.
If two systems can be represented by the same mathematical formula then they are functionally the same. The bucket with the water divided into sections and never mixing has the same mathematical formula as the Berne model. Therefore the Berne model is representing an atmosphere where CO2 never mixes.
Phil,
Mother Nature “fails to do” what? The tonnage of new vegetation/plankton/animalea every year has no direct relationship with the tonnage of rot. The oceans drink or spew CO2 as they will. We search for a formula to deduce the future of a random walk here, no proof it is anything else. Search away…
Ferdi: “That makes that the partial pressure in seawater at equilibrium increases with ~16 microatm. That means that X increases:
Xi = X/(750-400)*(766-400) = 1.046 X
…
After some time, usually 1-3 years, the CO2 level increased to 416 ppmv. That gives:
Xi = X/(750-400)*(766-416) = 1.000 X”
How do you get the 1-3 years here? At current rates 16ppmv will take about 8 years based on the 750 figure for hottest water. A more typical value you make that more. Since the process seems to be dominated by AO we should probably be looking at colder waters.
This needs looking at with some proper numbers but may be a means to apportion the ratio of outgassing and residual aGHG.
Clive Best says:
July 2, 2013 at 10:41 am
Carbon sinks may increase or decrease with rising CO2 levels. However for now lets assume that carbon sinks remain constant so we can take a single effective “pulse” e-folding time – tau years. Man made emissions of fossil fuels are currently running at 5.5 Gtons per year, and the atmosphere currently contains 750g tons CO2.
So you accept that the carbon dioxide sinks could be variable – but let’s disregard reality and do some maths? Are you a climate modeler by any chance? 😉
tumetuestumefaisdubien1 says:
July 2, 2013 at 9:14 am
If we would assume the temperature increase and CO2 outgassing only in the upper 100m epipelagic zone of the ocean* then we can calculate the temperature dependent CO2 outgassing in absolute numbers as:
You make the fundamental omission that outgassing of the oceans increases the CO2 content (= partial pressure) of the atmosphere. An increase of 16 ppmv is sufficient to fully compensate for the increase of 1 K in temperature of the entire ocean. It doesn’t matter how much CO2 the oceans contain, it only matters what the pressure difference between CO2 in the oceans at equilibrium and in the atmosphere is.
Take a few bottles of coke: 0.5, 1 and 1.5 liters. Shake well (with the still closed bottles!). Measure temperature and pressure. For the same fill of coke and at the same temperature, the pressure in the different bottles will be equal (except that the relative loss of CO2 from the smaller amount of liquid is somewhat higher). Thus the total amount of CO2 doesn’t play any role, only pressure is important, which is influenced by concentrations and temperature.
Ferdi says: ” I have no idea what drives the connection between AO and the CO2 rate of change”
Err, pressure perhaps? AO in pressure based index. If atmospheric pressure drops , so does 400 ppmv worth of partial pressure. That will affect out-gassing, requiring more ppmv CO2 to maintain the same balance at a lower pressure.
http://climategrog.wordpress.com/?attachment_id=231
The graph shows AO rose between 1960 and 1995 and has ‘plateaued’ since. So has d/dt(CO2) so has temperature.
Human emissions on the other hand …
The biggest mystery is why CO2 levels in the atmosphere are naturally so low. Why have levels been around ~300ppm for the last 3 million years, while higher levels existed in the past ? What mechanism determines the “equilibrium” level for CO2? Does life determine how much CO2 remains in the atmosphere ?
I think one clue to this mystery is the remarkable fact that atmospheric CO2 radiation to space is currently maximised at ~300ppm. Under today’s climate conditions with Tavg ~ 288K maximum cooling of the atmosphere is ensured with 300ppm. This is the second law of thermodynamics at work in the atmosphere. Much below 200ppm and plants will stop growing and CO2 levels rise – so I suspect life really does determine the climate on Earth !
P.S. Humans are also part of life.
Michael Moon says:
July 2, 2013 at 11:05 am
Phil,
Mother Nature “fails to do” what? The tonnage of new vegetation/plankton/animalea every year has no direct relationship with the tonnage of rot. The oceans drink or spew CO2 as they will. We search for a formula to deduce the future of a random walk here, no proof it is anything else.
Except in your ‘random walk’ it’s always going the same way! Every year sources exceed sinks by about the same amount, where are the years when sinks exceed sources as you’d expect with a random walk? That’s not a random walk!
Greg Goodman says:
July 2, 2013 at 11:07 am
How do you get the 1-3 years here? At current rates 16ppmv will take about 8 years based on the 750 figure for hottest water.
Sorry, I was not clear: the 1-3 years is the half life time for equilibrium between the ocean surface and the atmosphere, not an absolute figure. For huge changes like an instant 1 K increase, that indeed will take more years, for smaller changes that will be shorter as most changes are reversed within the 1-3 years time frame…
Greg Goodman says:
July 2, 2013 at 7:15 am
Ferdi says: “Further, any extra natural release from the oceans or vegetation would show up in the 13C/12C ratio’s of the atmosphere. The biosphere is a net sink of CO2, preferably of 12CO2, thus increasing the 13C/12C ratio in the atmosphere. The oceans have way higher 13C/12C ratio’s than the atmosphere, thus should increase the ratio in the atmosphere. But we see a firm decline…”
Then some of our trivial assumptions about the carbon cycle are wrong….
>>>>>>>>>>>>>>>>>>>>
Yes, Take a look at The Trouble With C12 C13 Ratios
And speaking of atom bomb testing….
Plastics allow migration too which is why soda in an plastic bottle will go flat as it ages.
@Clive Best
This graphic makes it hard for me to believe that tau can be measured in years:
http://ds.data.jma.go.jp/ghg/kanshi/co2map/gmapplot_e.html
It comes up displaying wintertime China. Change the month from 12 to 7. The anomalies are now localized and mostly minimal. Rotate the globe to check other northern hemisphere locations. The US is a net sink for CO2 in summer along with other continental areas with active plant life. Kingdom Plantae eagerly devours all of the CO2 we can produce in a matter of days or weeks when they are not dormant.
I expected paper #1 from Gösta Pettersson to say something about isotope fractination during absorbtion and outgassing in the oceans, but I can’t seem to find it.
As 14CO2 molecules are heavier and slower than nornal 12CO2 it’s easier captured by the water and more difficult to outgas, therefore one would expect that the overturning of 14CO2 is quicker than normal CO2. The question is if this is quantitatively significant or not.
Phil,
Always? Really? The Earth is old, the Scripps Institute is young. There is a lot of evidence that CO2 has varied in the past. The farther back we look, the more variation it has, as high as 7000 ppm according to some sources. Ice cores represent an average of the concentration in the 80 to several hundred years it took the firn to close up, and there is always some diffusion within even solid ice.
“Always” is a dangerous word, back it up for us?
Author’s paper 1. :
Among proponents of the AGW hypothesis, the ability of the Bern model to simulate the Keeling
curve has lent credence to the model and its ability to predict future levels of airborn carbon dioxide for presumed emission scenarios. The Bern model (or closely related carbon cycle models tuned to the Keeling curve) is routinely used by climate modellers to obtain the carbon dioxide input data they require to arrive at predictions of anthropogenic effects on the future climate.
===
The idea that the Bern model reproduces Keeling curve is frivolous. It is easy to “reproduce” a short section of a cumulative integral like that by any number of models. Whether is truely reflects the system behaviour is better seen by plotting the rate of change.
One thing is clear, feed the continually increasing human emission totals into the Bern model will not produce a plateau after 1995.
http://climategrog.wordpress.com/?attachment_id=259
Gösta,
Please study carefully the difference between tracer mixing and pulse diffusion uptake.
Bomb test data only measure tracer mixing. Peter Dietze published already in 1997 on the problems with the Bern Model. He found a pulse diffusion half life time of 38 years (e-folding time 55 years).
http://www.john-daly.com/carbon.htm
In the following graph this difference between the Bern Model and a constant rate diffusion (as observed) is shown.
http://members.casema.nl/errenwijlens/co2/co2afname.gif
leftturnandre says:
I expected paper #1 from Gösta Pettersson to say something about isotope fractination during absorbtion and outgassing in the oceans, but I can’t seem to find it.
As 14CO2 molecules are heavier and slower than nornal 12CO2 it’s easier captured by the water and more difficult to outgas, therefore one would expect that the overturning of 14CO2 is quicker than normal CO2. The question is if this is quantitatively significant or not.
===
So any out-gassing by the oceans and the massive annual too and fro will deplete the atmospheric dC13 that Ferdinand is telling me must be a sign of anthropogenic CO2. Like I said earlier, some of our assumptions about the carbon cycle seem badly flawed.
“I expected paper #1 from Gösta Pettersson to say something about isotope fractination during absorbtion and outgassing in the oceans, but I can’t seem to find it.”
Did you read paper 1 ?
see Page 7.
Greg Goodman says:
July 2, 2013 at 11:20 am
Ferdi says: ” I have no idea what drives the connection between AO and the CO2 rate of change”
Err, pressure perhaps? AO in pressure based index. If atmospheric pressure drops , so does 400 ppmv worth of partial pressure. That will affect out-gassing, requiring more ppmv CO2 to maintain the same balance at a lower pressure.
As far as I have found some real figures, the Arctic pressure at sealevel seems to change some +/- 20 mbar with the AO, that means a variability of +/- 8 microatm of pCO2 pressure. Not really a huge difference… A few years increase at the current rate will already exceed the variability in uptake.
I look up the formula in AR4 page 213 and I totally disagree that this represents a system where CO2 does not mix.
The formula is:
P(t) = 0.217 + 0.259*e(-t/172.9) + 0.338*e(-t/18.51) + 0.186*e(-t/1.186)
Where t is in years, and P(t) is the pulse.
This formula is not very uncommon or exotic in any way. This is a perfectly normal way of modeling a system with several different sink rates.
In words it says that if you add one unit of CO2 to the atmosphere then:
21,7 % of it will remain there indefinitely
25,9% will have a lifetime of 172.9 years
33.8% will have a lifetime of 18.51 years
18.6% will have a lifetime of 1.186 years
I have not gone further to investigate the justification for each parameter, but one can assume that the long lifetime of 172.9 years represent the sink rate into the deep ocean, the middle lifetime of 18.51 years may represent surface oceans and the low lifetime of 1.186 years to represent the biosphere exchange. However, whether this assumption is correct or not is not important, the important thing is that they represent different sink rates all acting on one mixed content of atmosphere.
You may also use this formula for calculating the depletion of the water level in a bucket with three holes, and the water can mix freely.
Stephen Wilde says:
July 2, 2013 at 2:39 am
….Murry Salby also suggests soil moisture on land as a significant player.
>>>>>>>>>>>>>
I would agree with him. Think of how limestone caves are formed. link
Too bad the author of that article got the chemistry wrong. (My thesis topic) You get differential dissolving of limestone beds based on the amount of clastic (sand and clay particles) in the limestone. The more clastic the longer the dissolving rate. The amount of surface area available to be dissolved by the H2CO3 determines the rate of reaction. The clastics effectively ‘hide’ or mask some of the limestone and protect it from dissolving. Erosion does not explain the differences in how the different bedding dissolves that is seen in cave formation but differences in dissolving rate do.
This photo sort of shows how the width of the cave changes as the groundwater dissolves its way through different bedding planes.
Greg Goodman says:
July 2, 2013 at 12:17 pm
So any out-gassing by the oceans and the massive annual too and fro will deplete the atmospheric dC13 that Ferdinand is telling me must be a sign of anthropogenic CO2.
Please reread my former comment: the back and forth release of CO2 from the oceans drops the d13C level with about 8 per mil. That is sufficient to maintain the difference between the ocean surface d13C at +1 to +5 per mil and the pre-industrial d13C level in the atmosphere of -6.4 +/- 0.2 per mil. But in no way that can explain the drop of 1.6 per mil since 1850, in lockstep with human emissions…
@Ian W.
Actually no !
I am proposing that if it were still possible to ignore all the doom mongers and simply continue on “business as usual” (thereby improving the lives of nearly everyone on Earth) my simple “maths” says that at worst CO2 levels should stabilize at ~1000ppm for the indefinite future leading to about 2C warming (using latest measurements).
Instead it looks likely we will suffer a much more alarming fate chasing phantoms !
fhhaynie says:
July 2, 2013 at 8:30 am
Ferdinand,
The biosphere is not a net sink for C12.
Yes it is, the earth is greening, thanks to all that extra CO2 in the atmosphere. And it is calculated from the oxygen balance: there is about 1 GtC more CO2 uptake by the whole biosphere than release:
http://www.bowdoin.edu/~mbattle/papers_posters_and_talks/BenderGBC2005.pdf
Richard M says:
July 2, 2013 at 9:55 am
One explanation for lower CO2 levels since the Eocene (down from 1000 ppmv to ~300) has been the putative Azolla Event.
http://en.wikipedia.org/wiki/Azolla_event
But IMO a cooling Earth is sufficient explanation. In other Ice House phases, as during the late Carboniferous & early Permian Periods, CO2 appears to have dropped into the 300s ppmv of dry air.
CACCAs use the AE to argue for lowered CO2 causing Cenozoic cooling, but IMO the lower CO2 is a result of cooling from other “forcings”, such as the arrangement of continents & orbital mechanics. There could of course be feedback effects.