Bern Model when emissions stabilised at 2013 levels

The Bern model has been used by IPCC for emission scenarios since SAR in 1995. We will now use it to extrapolate CO2 levels 500 years into the future with annual emissions fixed at 2013 levels (~10 GT C/year). All values are expressed as an equivalent atmospheric CO2 concentration in ppm for convenience.

To do this we integrate emissions over time and and calculate atmospheric CO2 concentrations.

CO_2(t) = CO_2(1750) + fac \times \int_{1750}^{t} Em(t') \times ( a_0 + \sum_{1}^{3} a_i e^{-\frac{t-t'}{\tau_i}} ) dt'  

where fac converts GTC to ppm, Em = emissions in year t’ and the a’s and are \tau's are the parameters of the Bern model. CO_2(1750) = 280ppm . For the integration I use this dataset for the emissions data.

The Bern model simply describes the time decay of an annual pulse of CO2 added to the atmosphere. It is parameterised as follows.

\Delta C = a_0 + \sum_{1}^{3} a_i e^{-\frac{t}{\tau_i}}

where for AR4:  a_0 = 0.217,  a_1 = 0.259,  a_2 = 0.338,  a_3 = 0.186,  \tau_1 = 172.9 y ,  \tau_2 = 18.51 y , \tau_3 = 1.186 y

and for TAR:   a_0 = 0.152,  a_1 = 0.253,  a_2 = 0.279,  a_3 = 0.316

Here are  the results.

CO2 levels in the atmosphere following a stabilisation at 2013 emission rates held constant for 500 years.

Levels do not reach an equilibrium value in the Bern model because a fixed fraction (a0) in any year are assumed to remain for ever in the atmosphere. However even in this case levels rise (just) to ~1270 ppm with AR4 parameters and ~950 ppm using TAR values. The airborne fraction shows how it falls asymptotically to the fixed retention rate.

Airborne fraction for the Bern model 500 years into the future for fixed 2013 emissions.

We can also ask the question as to how well the Bern model describes the measured CO2 levels between 1950 to 2016 based only on the emissions data. Here is that comparison.

Detail comparison. The dark green trace is the Moana Loa CO2 measurements .

The bern model agrees with the current value of CO2 (400ppm), but it does not give a good description of the trend. In fact the actual CO2 growth is slower than that of the model. This supports the hypothesis that a0 is actually very small and that CO2 levels will approach stability much faster.

Let’s consider two estimates of climate sensitivity (ECS) a median value of 2.5C (preferred by Gavin Schmidt) and a low value of 1.5C (preferred by say Nic Lewis).

ECS Model Version Net Warming in 2620
2.5C AR4 5.5C
2.5C TAR 4.2C
1.5C AR4 3.2C
1.5C TAR 2.5C

When you consider that we are assuming annual emissions held constant at ~10 GtC/y for 500y into the future, then these final result are not really that scary at all.

About Clive Best

PhD High Energy Physics Worked at CERN, Rutherford Lab, JET, JRC, OSVision
This entry was posted in AGW, Climate Change, climate science, Science and tagged , . Bookmark the permalink.

8 Responses to Bern Model when emissions stabilised at 2013 levels

  1. Lance Wallace says:

    You have two green lines that seem to overlap–Mauna Loa CO2 measurements and then something called Av-Measurements (undefined). Is that the average of all other CO2 measuring locations?

  2. Ron Graf says:

    Clive, this is very informative. It looks from your chart that Bern Model over-estimates the trend for CO2 accumulation (or underestimates the CO2 natural sinks). I would also be interested to see a chart of a simulation of emissions diminishing over the next 83 years to zero, to reflect an assumption that we would be an entirely nuclear fusion economy by then, able to create synthetic carbon fuels cheaply be then by reduction of CO2.

    On the quoting of Gavin’s ECS, I don’t dispute it but I would love to have a neat reference if you have, and also for other leading IPCC figures.

    • Clive Best says:

      It would be easy to put in a gradual reduction into the emissions trend. However, the physical assumption built in to the Bern model is that 20% of emissions from any year remain in the atmosphere forever. This cannot make physical sense. The only thing that does make sense is that the balance of PCO2 is reached first with the surface layer which then ‘diffuses’ down to the deep ocean over a longer time scale, allowing more absorption. Eventually all the ocean is ‘balanced’ so supposedly the rest has to wait for geological burial, but why 20% and not say 2%? Why would it even be constant?

      Gavin quote is from
      OK he says 2.5C-3.0C, but I seem to remember somewhere else when he was forced to give a preference and he said 2.5C

      On Twitter
      @ClimateOfGavin @Balinteractive When you say climate sensitivity is probably around 2.5-3C I assume you mean ECS. Correct?
      6:14 PM – 16 Oct 2015

      Gavin Schmidt ?@ClimateOfGavin 16 Oct 2015
      @clivehbest Yes. @Balinteractive
      1 reply 0 retweets 0 likes

  3. Randall says:

    None of this has any relevance because the basic premise that carbon dioxide has any affect on climate is wrong. Solar energy controls climate.

  4. Jim Bouldin says:

    Hi Clive,

    The form of the atm decay function is an interesting question w/ uncertainty. But I’m wondering where you got your model parameters because they don’t match the ones given by Joos et al 2013, from which I got mine. Thanks.

    • Clive Best says:

      I programmed up the Joos et al. 2013 values and although they are very different the end result is more or less the same as those quoted and displayed above in AR4 (labelled here as AR5)

      • Jim Bouldin says:

        Thanks for so doing Clive. I have to believe the high similarity of result with very different parameters, means simply a difference in time step between the two.

        I’m thinking about this issue a lot and have a post partially written but I’m moving at a glacial pace right now. Based on some fairly simple math, I cannot in any way see how they come up with such long, slow relaxation rates to a high asymptote. Strongly skeptical on that concept. And this in turn has implications for estimated GWPs (global warming potentials), via direct effects on time-integrated RF.

        Your statement in the other post, to a commenter, regarding the idea that continually increasing emissions is masking detection of faster relaxation times, is both interesting and quite important. I’m thinking a lot about that.

        This is quite a big deal.

Leave a Reply