**- JoNova - http://joannenova.com.au -**

New Science 3: The Conventional Basic Climate Model — In Full

Posted By __David Evans__ On September 25, 2015 @ 1:57 am In __Global Warming__ | __Comments Disabled__

A feast. A feast! For those who want the meat, the math and the diagrams (*don’t miss the diagram*s). As far as we know, this is the first time the architecture of the basic climate model has been laid out in one place on the Net. This is the most math heavy post this series, but it has to be done, and properly. This is where the 1.2 °C direct effect of doubling CO2 gets amplified to 2.5 °C with fairly basic physics. If the equations are not your forte, look at the “the Establishment Case” below the equations to get some idea why establishment scientists find it mind-bendingly hard to imagine how climate sensitivity could possibly be much different.

For commenters who know there are problems with this model (don’t we all), one of the points of doing this is to get through to the establishment leaders — to speak their language instead of having separate conversations. Of course, for some minds it will not matter what skeptical scientists say, but for other, key people, *it will. *We would expect seeing the flaws laid out so clearly will undercut the implacable confidence of leaders, though they may not say so.

Again, this ties for Most Uncontroversial Post on this blog. Everything here, the IPCC would agree with (except maybe the last sentence).

Do admire the diagram –Figure 2. It’s no accident it is similar to an electrical circuit diagram. Modeling and feedbacks have been tested to the nth by thousands of electrical engineers around the world building things we use every day. Thinking of the climate model this way is a useful technique to figure out where it goes wrong.

In a sense this is where silicon chip wisdom starts to scrutinize the climate maze.

– Jo

Dr David Evans, 24 September 2015, David Evans’ Basic Climate Models Home, Intro, Previous, Next, Nomenclature.

Here is the conventional basic climate model, in full. It builds on the previous post in the series, which explained how the model worked in the case of no feedbacks and only a CO2 input. This post uses the same terminology, notation, and ideas, without necessarily explaining them again. This is “heaviest” post in the series.

Readers who don’t want to see details of the mathematical modeling might want to skip straight to the diagram and the calculation of ECS at the end — a general understanding of the previous post and of the model diagram below is sufficient to make sense of the ensuing posts in this series.

Consider the hypothetical situation where only the only things that can change are:

- Temperature. While some variables might alter the relativities between various mean atmospheric temperatures, we will take the temperature to be the one we are most interested in, the surface temperature
*T*_{S}. - The
*n*driver variables,*V*_{1},…,*V*, are the climate influences considered in the model, such as the CO2 concentration, the total solar irradiance (TSI), or the externally-driven albedo (EDA). The driver variables are not directly dependent on the temperature, each other, or the feedback variables._{n} - The
*m*feedback variables,*U*_{1},…,*U*, are the feedbacks that the model considers, such as the average lapse rate, the surface-driven albedo, or the average height of the water vapor emissions layer (WVEL). Each feedback variable is dependent on the temperature, but is not directly dependent on the driver variables or on the other feedback variables. Because they depend on temperature, as the Earth warms these variables further affect absorbed solar radiation (ASR) or outgoing longwave radiation (OLR), which affects warming, which further affects these feedback variables, and so on in a feedback loop._{m}

The previous post is the special case where *n* equals one, *m* is zero, and *V*_{1} is the base-2 logarithm of the CO2 concentration *L*.

We are going to focus on the balance of radiation entering and leaving the planet, so it is helpful to explicitly focus on the net top-of-atmosphere (TOA) downward flux, namely

If the temperature, the drivers, and the feedbacks all change, how much does the OLR change?

We generalize the process that led to Eq. (3) of the previous post. *G* depends on *T*_{S} as well as on each of the driver and feedback variables, but on nothing else because everything else is held constant. The arguments of *G* are treated as 1+*n*+*m* independent variables, each of which varies by an incremental amount (using the “*d*” notation). We apply the well-known calculus rule for partial derivatives: the total change in *G* is the sum of the changes due to each variable, and the change due to a variable is the product of the partial derivative of *G* in that variable with the actual change in that variable –

(The Σ symbols means to sum what follows, one term for each value of the index *i*, whose values are in turn 1, 2,…, *n* or *m*. Thus each Σ is the sum of *n* or *m* terms, where each of the *n* or *m* terms in that sum is like what follows the Σ but with the appropriate value substituted for *i*.)

As in the simple case, we convert this description of how OLR changes into a description of a balance between increases and decreases in net flux, by considering a move between two steady states. Assuming the planet stays close to steady state as the climate influences do their work, the incremental changes *dX* become small changes Δ*X*, so Eq. (2) becomes

The radiation is balanced in a steady state (i.e. *G* is zero), so there is no change in net downward flux between steady states, so Δ*G* is zero. Also, because ASR is not directly dependent on surface warming, by Eq. (4) of the previous post,

Applying these two relationships to Eq. (3) gives the equation describing the balance:

Ignoring the feedbacks (by pretending all the *U _{i}* and their derivatives are zero), this says that the increase in OLR due to surface warming is balanced by an increase in net downward flux (i.e. increases in ASR plus decreases in OLR) due to changes in the driver variables

Changes in the drivers can be viewed as a process with *n*+1 stages:

- In stage
*i*, where*i*is in {1,2,…,*n*}, the climate driver variable*V*changes by Δ_{i}*V*and everything else is held constant, which increases the net TOA downward flux_{i}*G*by

- In stage
*n*+1, the surface warms by Δ*T*_{S}holding all the drivers constant, which increases OLR by exactly enough to restore the net change in*G*to zero.

Or these *n*+1 stages can be viewed as *n*+1 processes that occur simultaneously, without significantly interfering with each other, and we add the OLR changes that each causes (the climate is linear for small changes, so processes and changes superpose).

Feedbacks complicate the balance in Eq. (5) by adding yet more downward flux (which can be negative) due to changing values of the feedback variables *U*_{1},…,*U _{m}*, which change in response to surface warming. These extra downward fluxes have to be balanced by yet more OLR due to surface warming, requiring more surface warming, which further changes the values of the feedback variables, and so on in a feedback loop.

Solving Eq. (5) for Δ*T*_{S}, the surface warming is

where the total feedback *f* is defined as

This expression for Δ*T*_{S} is the same as it is without feedbacks (i.e. with *f* equals zero) except that it is multiplied by

which is the distinctive form of a feedback *f* around a multiplier *λ*_{0} — apply Fig. 1 with *a* equal to *λ*_{0} and *b* equal to *f*.

Figure 1: The closed-loop and open-loop forms of a feedback loop are equivalent.

AR5 reports the individual and total feedbacks from the CMIP5 as (Table 9.5 and Fig. 9.43, and p. 591), in W m^{−2} °C^{−1}:

- water vapor +1.6±0.3
- lapse rate −0.6±0.4
- water vapor and lapse rate combined +1.1±0.2
- surface albedo +0.3±0.1,
- cloud +0.3±0.7,

and the total feedback is

(Incidentally, the loop gain, the total amplification going once around the feedback loop, is

It is less than one. If it were greater than or equal to one then the infinite sum implicit in the denominator in Eq.s (7) and (9) would not converge and the computed surface warming would be infinite.)

In Eq. (7) the contribution to the surface warming by every driver shares a part in common: the Planck sensitivity *λ*_{0}, modified by the feedbacks *f*. This inspires the concept of the “forcing” of a driver, as the part peculiar to the driver variable: the (instantaneous radiative) forcing of a change Δ*V _{i}* in the climate driver variable

(see AR5, section 8.1.1.1 and Fig. 8.1; this is the appropriate form of forcing for changes between steady states). The contribution the change in climate driver *V _{i}* makes to surface warming, before feedbacks, is the product of its forcing and the Planck sensitivity:

The (total) radiation imbalance is the total change in downward flux before feedbacks, or the sum of forcings due to the drivers:

The situation is illustrated in Fig. 2:

- The first driver,
*V*_{1}, is the no-feedbacks ASR*A*_{NF}, the portion of ASR that is unaffected by feedbacks in the period under consideration. (Moving between steady states, changes in ASR are due either to Δ*A*_{NF}or to albedo feedbacks. Changes in*A*_{NF}are caused by something external to Earth, such as a fluctuation in TSI or a change in externally-driven albedo.) - The second driver,
*V*_{2}, is the logarithm of the CO2 concentration,*L*. - The feedbacks are all combined into the total feedback
*f*, see Eq. (10) above. - The computational paths from the drivers to the surface warming Δ
*T*_{S}are provided by Eq. (7), noting that (see Eq. (5) of the previous post for the last)

Figure 2: Conventional basic climate model, for changes from one steady state to another. The radiation imbalances (or forcings) due to the drivers are added (the purple circle), then balanced by the responding increase in feedbacks-adjusted OLR due to surface warming that would occur if all else was constant.

Adding the forcings and applying the feedbacks in Fig. 2, the surface warming is

where

Each driver adds its forcing to the input of the Planck sensitivity.

The equilibrium climate sensitivity (ECS) is the surface warming Δ*T*_{S} when the CO2 concentration doubles (i.e. Δ*L* is one) and the other drivers are unchanged, which by Eq.s (16) and (17) is

This agrees with AR5 (p. 1033), which finds the ECS as likely to be 1.5°C to 4.5°C (for all greenhouse gases, not just for CO2, because it is for a doubling of equivalent CO2 concentration — AR5 p.1451, p.1452). Setting *f* to zero in Eq. (18) gives the no-feedbacks-ECS found in the previous post.

This calculation in Eq. (18) shows that the sensitivity of carbon dioxide is about 2.5 °C (though in AR4 (2007) the central value was higher, about 3.3 °C), though there is a pretty wide uncertainty. Note that the three parameters in Eq. (18) were obtained mainly by basic physics:

*D*_{R,2X}, the decrease in OLR if CO2 doubles holding all else constant, is based mainly on spectroscopy that can be replicated in laboratories.*λ*_{0}, the Planck sensitivity, comes from Stefan-Boltzmann’s law (though it gets modified by 20% or so because we are concerned here with the surface temperature rather than the radiating temperature of the Earth). Stefan-Boltzmann’s law is based on replicated laboratory physics, and is almost certainly beyond reproach.*f*, the total feedback, comes from the large computer models, but its main components are arguably based largely on known physical properties of water vapor and so are surely about right. Although the possibility exists that there is a large and unknown negative feedback that has inadvertently been omitted, it hasn’t been found despite decades of searching.

The first doubling of CO2 from pre-industrial levels level might occur by around 2070, and it is widely thought that the CO2 level will top out at most after about three doublings, in the next century or two. If three doublings from pre-industrial occur, the expected surface warming by Eq. (18) will be 7.5 [4.5,16.5] °C. Even 4.5 °C would be pretty disruptive, while the upper end possibilities would be disastrous. Hence the alarm.

Skeptics who are familiar with the conventional model have mainly accepted Eq. (18) but queried the value of *f*, arguing that because the model’s implications don’t match a lot of empirical evidence there must be a large negative feedback missing from *f*.

In this blog series we accept the parameters, including *f*, as roughly correct, but argue that the logic in arriving at Eq. (18) is faulty, and that when the modeling is done correctly the ECS is an order of magnitude lower.

Article printed from JoNova: **http://joannenova.com.au**

URL to article: **http://joannenova.com.au/2015/09/new-science-3-the-conventional-basic-climate-model-in-full/**

Click here to print.

Copyright © 2008 JoNova. All rights reserved.