No, David Stockwell, it’s not the Sun

A few days ago I posted a video that, in less than 4 minutes, rebutted 10 common climate myths. The most amusing was the emphatic No to the claim that it’s because the Sun is getting hotter. To a certain extent, this is all that one need say with respect to this claim. It really can’t just be the Sun. Since then, however, I’ve encountered David Stockwell’s work suggesting that it is the Sun because the rise in surface temperature correlates well with the cumulative Total Solar Irradiance (TSI).

Now, I haven’t actually read David’s papers so can’t really comment about what he’s proposed. You may think that it is somewhat unreasonable to criticise someone’s work without actually reading it. There is, however, a good reason for this. It’s because it really can’t be the Sun. Why is this? A very obvious reason relates to climate sensitivity. The surface temperature has increased by 0.85oC since pre-industrial times. There is still a 0.7 Wm-2 radiative imbalance. If you’re suggesting that this is a response to a change in solar forcing, then you’re suggesting that the change in solar forcing since pre-industrial times has resulted in an increase in equilibrium surface temperature of around 1oC.

The problem, though, is that the solar forcing today is only about 0.1 Wm-2 greater than it was in pre-industrial times. An increase in surface temperature of 1oC produces an increase in surface flux of around 5.4 Wm-2. That would imply an incredibly large amplification factor. There is no evidence to suggest that our climate is that sensitive to changes in solar forcing. Furthermore, if it were that sensitive to changes in solar forcing, it would also be very sensitive to changes in other external forcings, such as anthropogenic forcings. The climate can’t be incredibly sensitive to changes in solar forcing while having no sensitivity at all to changes in other forcings. It really doesn’t make sense.

One of the reasons I wanted to write this post, though, was that there was something I’d been meaning to work out. An alternative to the system being very sensitive to changes in solar forcing would be that, for some unknown reason, the surface temperature in the mid-1800s was 1oC below the actual equilibrium value. Since then, therefore, we’ve simply been recovering towards equilibrium. If the equilibrium temperature is To and the current surface temperature is T, then the amount of energy, dE, that the system will accrue in a time dt is
TempTime
where A is the surface area of the Earth.

However, as the Earth accrues energy the surface temperature, T, will rise and so we need some estimate of how this will change with energy. The land and atmosphere have a mass of around 1019kg and a specific heat capacity of 1000 J kg-1 K-1. This means that it would take 1022J of energy to increase the temperature of the land and atmosphere by 1oC. The oceans, however, have a heat capacity 100 times greater than that of the land and atmosphere. So, increasing the temperature by 1oC would require 1024 J. As a check of this, the ocean heat content has increased by about 2.5 x 1023 J since 1960, while surface temperatures have risen by about 0.6oC. That would imply 5 x 1023 J per degree, but 1024J per degree is close enough for what I’m trying to illustrate here. So, the other condition we need then is
dTdE
So, I wrote a short computer code that would work out how the temperature would change with time if initially T = 288 K, To = 289 K, the energy changes according to the first equation above, and the temperature changes according to the second equation. The result is shown in the figure below. Basically, if we were simply recovering towards equilibrium from being 1oC below equilibrium, it would take only 20 years for temperatures to rise by 0.85oC and within 60 years we’d be virtually at equilibrium. Also, the profile is clearly not linear. The dashed-line simply shows what would happen if we happened to end up 1oC degree above equilibrium. Bear in mind that my estimate for the relationship between energy and temperature is possibly on the high side, so it may actually be faster. Also, this is quite different to what would happen if a pulse of energy (from an ENSO event for example) were to suddenly heat the land and atmosphere. That would decay very quickly – probably only a few months.
TempTime
This post was partly just a bit of fun. I had really just wanted to check how long it would take for the temperature to rise if it just happened to start 1oC below equilibrium. Essentially, far faster than the timescale it has actually taken for surface temperatures to rise by 0.85oC, and with a temporal profile that is very different to what’s been observed. Ultimately, however, it wasn’t really necessary to do that calculation. Basically we know that it can’t be the Sun because that would imply that the climate is much more sensitive to changes in external forcings than other evidence suggests, and this would apply both to changes in solar forcings and to changes in anthropogenic forcings. You can’t suggest that it is both much more sensitive to changes in solar forcings than expected while, at the same time, being much less sensitive to anthropogenic forcings.

This entry was posted in Climate change, Climate sensitivity, ENSO, Global warming, IPCC and tagged , , , , . Bookmark the permalink.

18 Responses to No, David Stockwell, it’s not the Sun

  1. I was leaving this to see if anyone noticed, but I think should probably have used T = 254 K and To = 255 K, rather than T = 288 K and To = 289 K. It doesn’t, however, make much difference. Temperature increases by 0.85oC in about 30 years and is virtually at equilibrium by about 80 years.

  2. Rachel says:

    I have to confess that I tend to skim over the bits with numbers.

  3. As I suspect most do. That’s why I get away with so many mistakes 🙂

  4. OPatrick says:

    Talking of numbers could you clarify this value: “The land and atmosphere have a mass of around 1019kg” – how do you decide which part of ‘the land’ to include in the mass calculation?

  5. OPatrick, it is an approximation, but what I did there was the following. The atmosphere has a mass of 5 x 1018 kg. If you consider the land to be the top 1m, then you get a mass of about 7 x 1017kg. So, I don’t actually know how deep one should consider the land in this context but assuming that radiative heating doesn’t extend below 10m, then 1019 kg seems like a reasonable upper limit for the land and atmosphere. If someone knows better than that, feel free to point it out.

  6. In fact, this figure here suggests the continents and atmosphere absorb about as much heat as each other. Since they have roughly the same specific heat capacity that would seem to indicate a similar mass is involved, so a total of about 1019 kg seems reasonable. Of course, as I suggest in the post, my factor of 100 greater for oceans is almost certainly too high, but that would just mean that it is quite a bit faster than I estimate here.

  7. Also talking of numbers, I finally have to ask the following (perhaps rather dumb) question regarding your flux estimate (as you’ve mentioned it in several posts now). You write: “An increase in surface temperature of 1K produces an increase in surface flux of around 5.4 W/m2.”
    I was always under the impression that you are talking about the Planck response here. However, it is usually specified as -3.2 W/m2/K, which leaves me with a change in surface flux of 3.2 W/m2. Am I missing something fundamental here? Thanks for enlightening me, Wotts 🙂

  8. Ooops, italics not properly closed … always reminds me that preview option would be great 😉

    Also talking of numbers, I finally have to ask the following (perhaps rather dumb) question regarding your flux estimate (as you’ve mentioned it in several posts now). You write: “An increase in surface temperature of 1K produces an increase in surface flux of around 5.4 W/m2.
    I was always under the impression that you are talking about the Planck response here. However, it is usually specified as -3.2 W/m2/K, which leaves me with a change in surface flux of 3.2 W/m2. Am I missing something fundamental here? Thanks for enlightening me, Wotts

  9. Karsten, I was talking of the Planck response. Indeed I simply gave the magnitude and not the sign, which I should have done. I’ve also ignored the emissivity (around 0.6 I believe). So, yes, it should really be -3.2 Wm-2. However, to be fair, I did specifically say “surface flux”, but you’re correct that including the emissivity would be more appropriate (assuming that’s what you’re suggesting).

  10. OPatrick says:

    Thanks. I ask partly because I was recently involved in a conversation with someone who claimed that as the Earth has a mass of 5.98 x1024 kg we are only experiencing a temperature increase of 1.67 x 10-6 degrees per year.

    Incidentally, can I use superscript in comments? It ignored my ‘sup’ request.

  11. You should be able to use sup. I’ve fixed your comment by doing exactly that.

    The Guardian comment is quite remarkable. I think your response to it was probably all it deserved.

  12. Okay, thanks Wotts. Just double-checked, the feedback/response convention seems to be Earth’s radiating (average) temperature, which makes it -3.2W/m2/K in GCM’s. For the surface, your blackbody estimate of -5.4W/m2/K is about right. The difference determines Earth’s emissivity. So surface flux change and Planck response doesn’t seem to be the same thing after all. Guess it was me who messed things up a bit here 😉

  13. BBD says:

    Oh look. Two scientists exploring a problem admit errors, refine analysis, seek improved understanding.

    Look, “sceptics”. Look.

  14. Karsten, to be honest I’ve become a little confused about which is more appropriate for determining the amplification/feedback. I was reading Lockwood (2010) which seems to use the change in surface flux compared to the change in forcings to estimate an amplification factor (or, feedbacks I guess). Having said that, it’s been a long day and I may be just get myself more confused if I continue to think about this any longer 🙂

  15. DocMartyn says:

    “If the equilibrium temperature is To and the current surface temperature is T, then the amount of energy, dE, that the system will accrue in a time dt is………………”

    Do you think that the term ‘equilibrium temperature’ is a correct usage for the average of (Tmax + Tmin)/2, taken at points of a rotating spheres surface?

    If I were to pace a meter long steel rod into a fire, would the ‘equilibrium temperature’ be (Hot Tip + Cool Tip)/2 ?

  16. Wotts, I’ve just checked Lockwood 2010 and you’re absolutely spot on. 5.4W/m2/K is the Planck response of the surface temperature change. The confusion arises from the fact that the Planck response in GCM’s is expressed in terms of flux changes at the top of the atmosphere. Therefore, it’s a function of Earth’s total emissivity ( = 1/4*S(1-A) / (sigma*Ts^4) = 0.61 ; S=SW in ; A=albedo ; sigma=Boltzm.C. ; Ts=sfc Temperature ). With this number, we get straight to the model estimate of approx. 3.2W/m2/K. Guess that’s about it.

    A minor nit though: It should probably be sigma*T^4 twice in your first equation (rather than T^2 for the first term).

  17. DocMartyn, I wasn’t meaning the average though. I was simply meaning the temperature at which the system would be in radiative equilibrium.

    Karsten, indeed – thanks – a typo.

  18. Pingback: Nigel Lawson and the Royal Society | Wotts Up With That Blog

Comments are closed.