Feeds:
Posts
Comments

Archive for 2010

In 1967 Journal of Atmospheric Sciences published the paper: Thermal Equilibrium of the Atmosphere with a Given Distribution of Relative Humidity by Manabe and Wetherald.

Here is one interesting model projection:

 

Model predictions 1967

Model predictions 1967

 

The corresponding note says:

 

Stratospheric cooling from increasing CO2

Stratospheric cooling from increasing CO2

 

Can this be true? How can “greenhouse” gases reduce temperature? Is this another “global warming causes more snow storms” type story?

First, a little about the stratosphere.

Stratospheric Basics

 

Atmospheric Pressure and Temperature, Bigg (2005)

Atmospheric Pressure and Temperature, Bigg (2005)

 

The stratosphere is the region of the atmosphere from around 10km to 50km. In pressure terms it’s the pressure between about 200mbar and 1mbar.

Ultraviolet radiation is almost completely absorbed in the stratosphere. The high energy photons of wavelength less than 0.24μm can break up molecular oxygen, O2, into atomic oxygen, O+O.

O2 and O combine to create O3, or ozone, which is again broken up with absorption of more ultraviolet.

Ozone production is greatest at a height around 25km. At higher levels, there are too few oxygen molecules to intercept all of the photons. At lower levels, there are few high energy photons left.

Here’s an interesting way of seeing how the absorption of solar energy at different wavelengths changes as thicker sections of the atmosphere,  especially the stratosphere, are traversed:

 

Absorption effects of different "amounts" of the atmosphere, Taylor (2005)

Absorption effects of different “amounts” of the atmosphere, Taylor (2005)

 

The reason why the troposphere (lower atmosphere) warms from the bottom is that once the UV is absorbed the atmosphere is mostly transparent to the rest of the solar radiation. Therefore, the radiation passes straight through and is absorbed by the earth’s surface, which warms up and consequently warms the atmosphere from beneath.

Air that warms expands, and so rises, causing convection to dominate the temperature profile of the lower atmosphere.

By contrast, the stratosphere is warmer at the top because of the effect of solar absorption by O2 and O3. If there was no absorption by O2 or O3 the stratosphere would be cooler at the top (as it would only be heated from underneath by the troposphere).

Just about everyone has heard about ozone depletion in the stratosphere due to CFCs (and other chemicals). Less ozone must also cause cooling in the stratosphere. This is easier to understand than the model results at the beginning (from increased “greenhouse” gases). Less ozone means less ability to absorb solar radiation. If less energy is absorbed, then the equilibrium stratospheric temperature must be lower.

Stratospheric Temperature Trends

Temperature measurements of the stratosphere are limited. We have satellite data since 1979 which doesn’t provide as much vertical resolution as we need. We have radiosonde data since the 1940s which is limited geographically and also is primary below 30hPa (around 25km).

Lots of painful work has gone into recreating temperature trends by height/pressure and by latitude. For example, in the 2001 review paper by Ramaswamy and many co-workers (reference below), the analysis/re-analysis of the data took 23 of the 52 pages.

Here is one temperature profile reconstruction from Thompson and Solomon:

 

Stratospheric Temperature Trends 1979-2003, Thompson (2005)

Stratospheric Temperature Trends 1979-2003, Thompson (2005)

 

From Thompson & Solomon (2005):

From 1979 to 1994, global-mean stratospheric temperatures dropped by 0.75 K / decade in the stratosphere below 35 km and 2.5 K / decade near 50 km
Another reconstruction from Randel (2008):

 

Stratospheric temperature trends by pressure, 1979-2007, Randel (2008)

Stratospheric temperature trends by pressure, 1979-2007, Randel (2008)

 

Before explaining why more CO2 and other trace gases could cause “stratospheric cooling”, it’s worth looking at the model results to understand the expected temperature effects of less ozone – and more CO2.

Observations and Recent Model Results

Notice that in the 1967 paper the predicted temperature drop was larger the higher up in the stratosphere. The effects of ozone are more complex and also there is more uncertainty in the ozone trends because ozone depletion has been more localized.

Here are model results for ozone – the best estimate of the observed temperature changes are in brown but aren’t expected to match the models because ozone is only one of the factors affecting stratospheric temperature:

 

Stratospheric observations and models, Shine (2003)

Stratospheric observations and models for ozone changes, Shine (2003)

 

Note that the effect of ozone depletion has a projected peak cooling around 1hPa (50km) and a second peak cooling around 80hPa.

Now the same paper reviews the latest model results for stratospheric temperature from changes in “greenhouse” gases:

 

Stratospheric observations and models for "greenhouse" gas changes, Shine (2003)

Stratospheric observations and models for “greenhouse” gas changes, Shine (2003)

 

The same paper reviews the model results for changes in stratospheric water vapor. This is a subject which deserves a separate post (watch this space):

 

Stratospheric observations and models for water vapor, Shine (2003)

Stratospheric observations and models for water vapor, Shine (2003)

 

Finally, the model results when all of the effects are combined together:

 

Stratospheric observations and models for ozone, GHG and water vapor changes, Shine (2003)

Stratospheric observations and models for ozone, GHG and water vapor changes, Shine (2003)

 

The model results are a reasonable match with the observed trends – but a long way off perfect. By “reasonable match” I mean that they reproduce the general trends of decadal cooling vs height.

There are many uncertainties in the observations, and there are many uncertainties in the changes in concentration of stratospheric ozone and stratospheric water vapor (but not so much uncertainty about changes in the well-mixed “greenhouse” gases).

A couple of comments from A comparison of model-simulated trends in stratospheric temperatures, by Shine et al, first on the upper stratosphere, reviewing possible explanations of the discrepancies:

None of these potential explanations is compelling and so the possibility remains that the discrepancy is real, which would indicate that there is a temperature trend mechanism missing from the models.

and then on the 20-70hPa region:

Nonetheless, assuming that at least some part of this discrepancy is real, one possible explanation is stratospheric water vapour changes. Figure 3 indicates that an extra cooling of a few tenths of a K/decade would result if the Boulder sonde-based water vapour trends were used rather than the HALOE water vapour trends. If this were one explanation for the model–observation difference, water vapour could dominate over ozone as the main cause of temperature trends in this altitude region.

Why Is the Stratosphere Expected to Cool from Increases in “Greenhouse” Gases?

This is a difficult one to answer with a 30-second soundbite. You can find a few “explanations” on the web which don’t really explain it, and others which appear to get the explanation wrong.

The simplest approach to explaining it is to say that the physics of absorption and emission in the atmosphere – when calculated over a vertical section through the atmosphere and across all wavelengths – produces this result. That is – the maths produces this result..

You can see an introduction to absorption and re-emission in CO2 – An Insignificant Trace Gas? Part Three.

[Note added to this article much later, the series Visualizing Atmospheric Radiation has an article Part Eleven – Stratospheric Cooling – from January 2013 on why the stratosphere is expected to cool as CO2 increases. It is quite involved but shows the detailed mechanism behind stratospheric cooling].

After all, this approach is what led Manabe and Wetherald to their results in 1967. But of course, we all want to understand conceptually how an increase in CO2 – which causes surface and troposphere warming – can lead to stratospheric cooling.

The great Ramanathan in his 1998 review paper Trace-Gas Greenhouse Effect and Global Warming (thanks to Gary Thompson of American Thinker for recommending this paper) says this:

As we mentioned earlier, in our explanation of the greenhouse effect, OLR reduces (with an increase in CO2) because of the decrease in temperature with altitude.

In the stratosphere, however, temperature increases with altitude and as a result the cooling to space is larger than the absorption from layers below. This is the fundamental reason for the CO2 induced cooling.

In Ramaswamy (2001):

For carbon dioxide the main 15-um band is saturated over quite short distances. Hence the upwelling radiation reaching the lower stratosphere originates from the cold upper troposphere. When the CO2 concentration is increased, the increase in absorbed radiation is quite small and the effect of the increased emission dominates, leading to a cooling at all heights in the stratosphere.

Are they saying the same thing? Yes (probably).

If these explanations help – wonderful. If they don’t, refer to the maths. That is, the mathematical result provides this solution and overall “hand waving” explanations are only ever a second-best “guide”. Also check out The Earth’s Energy Budget – Part Three for explanations about emissions from various levels in the atmosphere.

Conclusion

Understanding stratospheric temperature trends is a difficult challenge. Understanding the mechanisms behind this changes is much more of a conceptual challenge.

But over 40 years ago, it was predicted that the upper stratosphere would cool significantly from increases in CO2.

The depletion of ozone is also predicted to have an effect on stratospheric temperatures – in the upper stratosphere (where CO2 increases will also have the most effect) and again in the lower stratosphere where ozone is the dominant factor.

Stratospheric water vapor also has an effect in the lower stratosphere (where more water vapor leads to more warming and vice-versa), but more on this in a later post.

For some, who feel/believe that CO2 can’t really significantly affect anything in climate – this post isn’t for you – check out the CO2 – An Insignificant Trace Gas? series.

There will be others who will say “Ozone is the reason the upper stratosphere has cooled“. True, but increases in CO2 are also an important factor. The same calculations (maths and physics) that lead to the conclusion that less ozone will cool also lead to the conclusion that more CO2 will cool the upper stratosphere.

This subject also has two other possible consequences. One is about attribution. Global temperatures have increased over the last 40 years and many people want to understand the cause.

If solar heating was the direct cause (see Here Comes the Sun) the stratosphere would not be cooling. However, other effects could possibly also cause stratospheric cooling at the same time as tropospheric and surface heating. It’s a complex subject. But something to question for those other potential causes – would they also cause stratospheric cooling?

The other consequence is about GCMs. Some say that stratospheric cooling is a “vindication” of GCMs. In so far as we have covered the subject in this post we couldn’t reach that conclusion. The modeling of tropospheric and stratospheric temperature profiles can be done (and was by Manabe and Wetherald) with 1D radiative-convective models. Certainly 3d GCMs have also been used to calculate the effect by latitude but these results have more issues – well, the whole subject is much more complex because the change of ozone with height and latitude are not well understood.

But it is important to understand the difference between a GCM solving the general climate problem and a more constrained mathematical model solving the temperature profile against height through the atmosphere.

However, stratospheric cooling while the surface and troposphere are warming does indicate that CO2 and other “greenhouse” gases are likely influencers.

References

Thermal Equilibrium of the Atmosphere with a Given Distribution of Relative Humidity, Manabe and Wetherald, Journal of Atmospheric Sciences (1967)

Trace-Gas Greenhouse Effect and Global Warming, Ramanathan, Royal Swedish Academy of Sciences (1998)

Stratospheric Temperature Trends: Observations and Model Simulations, Ramaswamy et al, Review of Geophysics (2001)

A comparison of model-simulated trends in stratospheric temperatures, Shine et al, Q. J. R. Meteorol. Soc. (2003)

Recent Stratospheric Climate Trends as Evidenced in Radiosonde Data: Global Structure and Tropospheric Linkages, Thompson & Solomon, Journal of Climate (2005)

An update of observed stratospheric temperature trends, Randel, Journal of Geophysical Research (2008)

Read Full Post »

Many questions have recently been asked about the relative importance of various mechanisms for moving heat to and from the surface, so this article covers a few basics.

One Fine Day – the Radiation Components

 

Surface Radiation - clear day and cloudy day, from Robinson (1999)

Surface Radiation - clear day and cloudy day, from Robinson (1999)

 

I added some color to help pick out the different elements, note that temperature variation is also superimposed on the graph (on its own axis). The blue line is net longwave radiation.

Not so easy to see with the size of graphic, here they are expanded:

 

Clear sky

Clear sky

 

 

Cloudy sky

Cloudy sky

 

Note that the night-time is not shown, which is why the net radiation is almost always positive. You can see that the downward longwave radiation measured from the sky (in clear violation of the Imaginary Second Law of Thermodynamics) doesn’t change very much – equally so for the upwards longwave radiation from the ground. You can see the terrestrial (upwards longwave) radiation follows the temperature changes – as you would expect.

Sensible and Latent Heat

The energy change at the surface is the sum of:

  • Net radiation
  • “Sensible” heat
  • Latent heat
  • Heat flux into the ground

“Sensible” heat is that caused by conduction and convection. For example, with a warm surface and a cooler atmosphere, at the boundary layer heat will be conducted into the atmosphere and then convection will move the heat higher up into the atmosphere.

Latent heat is the heat moved by water evaporating and condensing higher up in the atmosphere. Heat is absorbed in evaporation and released by condensation – so the result is a movement of heat from the surface to higher levels in the atmosphere.

Heat flux into the ground is usually low, except into water.

 

Surface Heat Components in 3 Locations, Robinson (1999)

Surface Heat Components in 3 Locations, Robinson (1999)

 

All of these observations were made under clear skies in light to moderate wind conditions.

Note the low latent heat for the dry lake – of course.

The negative sensible heat in Arizona (2nd graphic) is because it is being drawn from the surface to evaporate water. It is more usual to see positive sensible heat during the daytime as the surface warms the lower levels of the atmosphere.

The latent heat is higher in Arizona than Wisconsin because of the drier air in Arizona (lower relative humidity).

The ratio of sensible heat to latent heat is called the Bowen ratio and the physics of the various processes mean that this ratio is kept to a minimum – a moist surface will hardly increase in temperature while evaporation is occurring, but once it has dried out there will be a rapid rise in temperature as the sensible heat flux takes over.

Heat into the Ground

 

Temperature at two depths in soil - annual variation, Robinson (1999)

Temperature at two depths in soil - annual variation, Robinson (1999)

 

We can see that heat doesn’t get very far into soil – because it is not a good conductor of heat.

Here is a useful table of properties of various substances:

The rate of heat penetration (e.g. into the soil) is dependent on the thermal diffusivity. This is a combination of two factors – the thermal conductivity (how well heat is conducted through the substance) divided by the heat capacity (how much heat it takes to increase the temperature of the substance).

The lower the value of the thermal diffusivity the lower the temperature rise further into the substance. So heat doesn’t get very far into dry sand, or still water. But it does get 10x further into wet soil (correction thanks to Nullius in Verba- really it gets 3x further into wet soil because “Thickness penetrated is proportional to the square root of diffusivity times time” – and I didn’t just take his word for it..)

Why is still water so similar to dry sand? Water has 4x the ability to conduct heat, but also it takes almost 4x as much heat to lift the temperature of water by 1°C.

Note that stirred water is a much better conductor of heat – due to convection. The same applies to air, even more so – “stirred” air (= moving air) conducts heat a million times more effectively than still air.

Temperature Profiles Throughout a 24-Hour Period

 

Temperature profiles throughout the day, Robinson (1999)

Temperature profiles throughout the day, Robinson (1999)

 

I’ll cover more about temperature profiles in a later article about why the troposphere has the temperature profile it does.

During the day the ground is being heated up by the sun and by the longwave radiation from the atmosphere. Once the sun sets, the ground cools faster and starts to take the lower levels of the atmosphere with it.

Conclusion

Just some basic measurements of the various components that affect the surface temperature to help establish their relative importance.

Note: All of the graphics were taken from Contemporary Climatology by Peter Robinson and Ann Henderson-Sellers (1999)

Read Full Post »

This post covers a dull subject. If you are new to Science of Doom, the subject matter here will quite possibly be the least interesting in the entire blog. At least, up until now. It’s possible that new questions will be asked in future which will compel me to write posts that climb to new heights of breath-taking dullness.

So commenters take note – you have a duty as well. And new readers, quickly jump to another post..

Recap

In an earlier post – Why Global Mean Surface Temperature Should be Relegated, Or Mostly Ignored – we looked at the many problems of trying to measure the surface of the earth by measuring the air temperature a few feet off the ground. And also the problems encountered in calculating the average temperature by an arithmetic mean. (An arithmetic mean for those not familiar with the subject is the “usual” and traditional averaging where you add up all the numbers and divide by how many values you had).

We looked at an example where the average temperature increased, but the amount of energy radiated went down. Energy radiated out would seem to be a more useful measure of “real temperature” so clearly arithmetic averages of temperature have issues. This is how GMST is calculated – well not exactly, as the values are area-weighted, but there is no factoring in of how surface temperature affects energy radiated.

But in the discussion someone brought up emissivity and what effect it has on the calculation of energy radiated. So in the interests of completeness we arrive here.

Emissivity of the Earth’s Surface

Our commenter asked:

So what are the non-black body corrections required for the initial calculation 396W/sqm? And what are the corrections for the equivalent temperature calculation? And do they cancel out (I think not due to the non-linearity issue) ?

What’s this about? (Of course, read the earlier post if you haven’t already).

Energy radiated from a body, E=εσT4

where T is absolute temperature (in K), σ=5.67×10-8 and ε is the emissivity.

ε is a value between 0 and 1, and 1 is the “blackbody”. The value – very important to note – is dependent on wavelength.

So the calculations I showed (in the thought experiment) where temperature went up but energy radiated went down need adjustment for this non-blackbody emissivity.

How Emissivity Changes

Here we consult the “page-turner”, Surface Emissivity Maps for use in Satellite Retrievals of Longwave Radiation by Wilber (1999).

Emissivity vs wavelength for various substances, Wilber (1999)

Emissivity vs wavelength for various substances, Wilber (1999)

And yet more graphs at the end of the post – spreading out the excitement..

Note the key point, in the wavelengths of interest emissivity is close to 1 – close to a blackbody.

For beginners to the subject, who somehow find this interesting and are therefore still reading, the wavelengths in question: 4-30μm are the wavelengths where most of the longwave radiation takes place from the earth’s surface. Check out CO2 – An Insignificant Trace Gas? for more on this.

I did wonder why the measurements weren’t carried on to 30μm and as far as I can determine it is less interesting for satellite measurements – because satellites can see the surface the best in the “atmospheric window” of 8-14μm.

So with the data we have we see that generally the value is close to unity – the earth’s surface is very close to a “blackbody”. Energy radiated in 4-16μm wavelengths only account for 50-60% of the typical energy radiated from the earth’s surface, so we don’t have the full answer. Still with my excitement already at fever pitch on this topic I think others should take on the task of tracking down emissivity of representative earth surface types at >16μm and report back.

So we have some ideas of emissivities, they are not 1, but generally very close. How does this affect the calculation of energy radiated?

Mostly Harmless

Not much effect.

I took the original example with 7 equal areas at particular temperatures for 1999 and show emissivities (these are arbitrarily chosen to see what happens):

  • Equatorial region: 30°C ;  ε = 0.99
  • Sub-tropics: 22°C, 22°C ;  ε = 0.99
  • Mid-latitude regions: 12°C, 12°C ;  ε = 0.80
  • Polar regions: 0°C, 0°C ;  ε = 0.80

The average temperature, or “global mean surface temperature” = 14°C.

And in 2009 (same temperatures as in the previous article):

  • Equatorial region: 26°C ;  ε = 0.99
  • Sub-tropics: 20°C, 20°C ;  ε = 0.99
  • Mid-latitude regions: 12°C, 12°C ;  ε = 0.80
  • Polar regions: 5°C, 5°C ;  ε = 0.80

The average temperature, or “global mean surface temperature” = 14.3°C.

The calculation of the energy radiated is done by simply taking each temperature and applying the equation above – E=εσT4

Because we are calculating the total energy we are simply adding up the energy value from each area. All the emissivity does is weight the energy from each location.

  • With the emissivity values as shown, the 1999 energy = 2426 W/ arbitrary area
  • With the emissivity values as shown, the 2009 energy = 2416 W/ same arbitrary area

So once again the energy radiated has gone down, even though the GMST has increased.

If we change around the emissivities, so that ε=0.8 for Equatorial & Sub-Tropics, while ε=0.99 for Mid-Latitude and Polar regions, the GMST values are the same.

  • With the new emissivity values, the 1999 energy = 2434 W/ arbitrary area
  • With the emissivity values as shown, the 2009 energy = 2442 W/ same arbitrary area

So the temperature has gone up and the energy radiated has also gone up.

Therefore, emissivity does change the situation a little. I chose more extreme values of emissivity than are typically found to see what the effect was.

The result is not complex or non-linear because emissivity simple “weights” the value of energy making it more or less important as the emissivity is higher or lower.

In the second example above, if the magnitude of temperature changes was slightly greater in the polar and equatorial regions this would be enough to still show a decrease in energy while “GMST” was increasing.

More Emissivity Graphs

Emissivity vs wavelength of various substances, Wilber (1999)

Emissivity vs wavelength of various substances, Wilber (1999)

Conclusion

Emissivity in the wavelengths of interest for the earth’s radiation is generally very close to 1. Assuming “blackbody” radiation is a reasonable assumption for most calculations of interest – as other unknowns are typically a higher source of error.

Because the earth’s surface has been mapped out and linked to the emissivities, if a particular calculation does need high level accuracy the emissivities can be used.

In the terms of how emissivity changes the “surprising” result that temperature can increase while energy radiated decreases – the answer is “not much”.

Read Full Post »

In On Having a Laugh – by Gerlich and Tscheuschner (2009) I commented that I had only got to page 50 and there were 115 pages in total.

Because there were so many errors already spotted, none central to the argument (the argument hadn’t started even at page 50), it seemed a pointless exercise to read it further. After all, many interesting papers await, on the thermohaline circulation, on models, on stratospheric cooling..

Perhaps most important of the criticisms was that Gerlich and Tscheuschner didn’t appear at all familiar with the climate science they were “debunking” – instead of commenting on encyclopedia references or throwaway comments in introductions to works unrelated to proving the inappropriately-named “greenhouse effect” they should be commenting on papers like Climate Modeling through Radiative-Convective Models by Ramanathan and Coakley (1978).

Clearly they were “having a laugh”

However, after noticing that a recent commenter actually cited Gerlich and Tscheuschner I went back and reviewed their paper. And in doing so I realized that many many misinformed comments by enthusiastic people on other popular blogs, and also this one, were included in the ground-breaking On Falsification Of The Atmospheric CO2 Greenhouse Effects by Gerlich and Tscheuschner.

It’s possible that rather than enthusiastic commenters obtaining misinformation from our duo that instead our duo have combined a knowledge of theoretical thermodynamics with climate science that they themselves obtained from blogs. The question of precedence is left as an exercise for the interested reader.

Miseducation

It is hard to know where to start with this paper because there is no logical flow.

Conductivity

The paper begins by reviewing the conductivity of various gases.

It is obvious that a doubling of the concentration of the trace gas CO2, whose thermal conductivity is approximately one half than that of nitrogen and oxygen, does change the thermal conductivity at the most by 0.03% and the isochoric thermal diffusivity at the most by 0.07 %. These numbers lie within the range of the measuring inaccuracy and other uncertainties such as rounding errors and therefore have no significance at all.

Clearly conductivity is the least important of means of heat transfer in the atmosphere. Radiation, convection and latent heat all get a decent treatment in studies of energy balance in the atmosphere.

If our duo had even read one book on atmospheric physics, or one central paper they would be aware of it.

Uninformed people might conclude from this exciting development that they have already demonstrated something of importance rather than just agreeing wholeheartedly with the work of atmospheric physicists.

Pseudo-Explanations to be Revealed in Part Two? Or Left as an Exercise for the Interested Student?

Following some demonstrations of their familiarity with mathematics and especially integration, they provide three conclusions, one of which refers to the Stefan-Boltzmann law, j=σT4:

The constant appearing in the T4 law is not a universal constant of physics. It strongly depends on the particular geometry of the problem considered.

and finish with (p21):

Many pseudo-explanations in the context of global climatology are already falsified by these three fundamental observations of mathematical physics.

Unfortunately they don’t explain which ones. The climate science world waits with baited breath..

The footnote to their comment on Stefan-Boltzmann:

For instance, to compute the radiative transfer in a multi-layer setup, the correct point of departure is the infinitesimal expression for the radiation intensity, not an integrated Stefan-Boltzmann expression already computed for an entirely divergent situation.

Sadly they are unfamiliar with the standard works in the field of the radiative-convective model.

Solar Energy Breakdown and A Huge Success in Miseducation

Solar Radiation Breakdown

Solar Radiation Breakdown

They followed up this table with the hugely popular comment:

In any case, a larger portion of the incoming sunlight lies in the infrared range than in the visible range. In most papers discussing the supposed greenhouse effect this important fact is completely ignored.

First, a comment on the “benefit” of this miseducation – being able to separate out solar radiation from terrestrial radiation is a huge benefit in climate understanding – it allows us to measure radiation at a particular wavelength and know its source. But many people are confused and say we can’t because 50% of the solar radiation is “infrared”. Infrared means >0.7μm. Conventionally, climate scientists use “shortwave” to mean radiation < 4μm and “longwave” to mean radiation > 4μm. As less than 1% of solar radiation is >4μm this is a very useful convention. Any radiation greater than 4μm is terrestrial (to 99% accuracy).

Many uninformed people who have become miseducated are certain that much solar radiation is >4μm – possibly due to confusing infrared with longwave.

We don’t speculate on motives on this blog so I’ll just point out that Gerlich and Tscheuschner know very little about any climate science, and from this comment probably don’t even understand the inappropriately-named “greenhouse” effect.

Why? Well, what has the visibility of the radiation have to do with the “greenhouse” effect? Of course it’s ignored. Our duo are just demonstrating their ignorance of the absolute basics.

Or they have some amazing insight into how the visibility or not of solar radiation affects the radiative transfer equations. All to be shared in part two probably..

The Core Question – the Radiative Transfer Equations

After a brief explanation of Kirchoff’s law, our duo discuss the core equations, the radiative transfer equations (RTE):

LTE [local thermodynamic equilibrium] does only bear a certain significance for the radiation transport calculations, if the absorption coefficients were not dependent on the temperature, which is not the case at low temperatures. Nevertheless, in modern climate model computations, this approach is used unscrupulously.

Absorption and emission coefficients get a very thorough treatment in the numerical solutions to the RTE, however, our duo are only familiar with work around the 1900’s and skip all modern work on the subject. Perhaps a more accurate statement would be:

We have no idea what anyone does but we read somewhere that stuff wasn’t done right..

Or they could actually show what effect that dependency actually had..

Then they decide to support the RTE:

Fantastic, 50 pages in we find the real RTE. This is what atmospheric physicists use to calculate the absorption and re-emission of radiation for each layer in the atmosphere. They follow this up with:

The integrations for the separate directions are independent of one another. In particular, the ones up have nothing to do with the ones down. It cannot be overemphasized, that differential equations only allow the calculation of changes on the basis of known parameters.

The initial values (or boundary conditions) cannot be derived from the differential equations to be solved. In particular, this even holds for this simple integral.

What do they mean? Of course you need boundary conditions to solve all real-world equations.

The separate directions are independent of one another? Yes, you find that in all treatments of radiative transfer.

So Gerlich and Tscheuschner agree that the RTE can be used to solve the problem? Or not? No one can tell from the comments here. If they do, the paper should be over now with support for the inappropriately-named “greenhouse effect”, unless they demonstrate that they can solve them for the atmosphere and get a different result from everyone else.

But they don’t.

Fortunately for those interested in what our duo really know and understand – they tell us..

The Modern Solution to the RTE – or How to Miss an Important 100 Years

After surveying works from more than 100 years ago, they conclude:

Callendar and Keeling, the founders of the modern greenhouse hypothesis, recycled Arrhenius’ discussion of yesterday and the day before yesterday by perpetuating the errors of the past and adding lots of new ones.

In the 70s and 80s two developments coincided: A accelerating progress in computer technology and an emergence of two contrary policy preferences, one supporting the development of civil nuclear technology, the other supporting Green Political movements. Suddenly the CO2 issue became on-topic, and so did computer simulations of the climate. The research results have been vague ever since.

No explanation of Callendar and Keeling’s mistakes – this is left as an exercise for the interested student.

And no mention of the critical work in the 1960s and 1970s which used the radiative transfer equations and the convective structure of the atmosphere to find the currently accepted solutions.

In fact, the research results haven’t been vague at all. Regular readers of this blog will know about Ramanathan and Coakley 1978, and there are many more specific papers which find solutions to the RTE – using boundary conditions and separation of upward and downward fluxes, as wonderfully endorsed by our comedic duo.

More recent work has of course refined and improved the work of the 1960s and 1970s. And the measurements match the calculations.

But what a great way to write off a huge area of research. Show some flaws in the formative work 100 or so years ago and then skip the modern work and pretend you have demonstrated that the modern theory is wrong.

As we saw in the last section, our duo appear to support the modern equations – although they are careful not to come out and say it. Luckily, they are blissfully ignorant of modern work in the field, which all helps in the miseducation of the uninformed.

The main work of the paper should now be over, but our duo haven’t realized it. So instead they move randomly to the radiative balance concept..

Radiative Balance and Mathematical Confusion

In every introduction to atmospheric physics you find the concept of radiative balance – solar energy absorbed = terrestrial radiation emitted from the top of the atmosphere. These concepts are used to demonstrate that the atmosphere must absorb longwave (terrestrial) radiation.

This concept can be found in CO2 – An Insignificant Trace Gas? Part One

After looking at the basics of the energy balance, they comment – on the right value for albedo (or ‘1-albedo’):

In summary, the factor 0.7 will enter the equations if one assumes that a grey body absorber is a black body radiator, contrary to the laws of physics. Other choices are possible, the result is arbitrary.

Being obscure impresses the uninformed. However, the informed will know that the earth’s emissivity and absorptivity will of course be different because the solar radiation is centered on 0.5μm while the terrestrial radiation is centered on 10μm. And the emissivity (and absorptivity) around 10um is very close to 1 (typically 0.98) while around 0.5μm the absorptivity is somewhat lower.

At this point, if we were to do a parody of our duo, we would write how their physics is extremely poor and do a three page derivation of absorptivity and emissivity as a function of wavelength.

Now follows many pages of maths explaining the impossibility of working out an average temperature for the earth during which they make the following interesting comment:

While it is incorrect to determine a temperature from a given radiation intensity, one is allowed to compute an effective radiation temperature Terad from T averages representing a mean radiation emitted from the Earth and to compare it with an assumed Earth’s average temperature Tmean

What they are saying is that for energy balance if we work out the radiation emitted from the earth we have dealt with the problem.

Fortunately for our intrepid duo, they are unacquainted with any contemporary climate science so the fact that someone has already done this work can be safely ignored. Earth’s Global Energy Budget by Trenberth, Fassulo and Kiehl (2008) covers this work.

To compute these effects more exactly, we have taken the surface skin temperature from the NRA at T62 resolution and 6-hour sampling and computed the correct global mean surface radiation from (1) as 396.4 W/m2. If we instead take the daily average values, thereby removing the diurnal cycle effects, the value drops to 396.1 W/m2 or a small negative bias. However, large changes occur if we first take the global mean temperature. In that case the answer is the same for 6 hourly, daily or climatological means at 389.2 W/m2. Hence the lack of resolution of the spatial structure leads to a low bias of about 7.2 W/m2. Indeed, when we compare the surface upward radiation from reanalyses that resolve the full spatial structure the values range from 393.4 to 396.0 W/m2.

The surface emissivity is not unity except perhaps in snow and ice regions, and it tends to be lowest in sand and desert regions, thereby slightly offsetting effects of the high temperatures on longwave (LW) upwelling radiation. It also varies with spectral band (see Chédin et al. 2004 for discussion). Wilber et al. (1999) estimate the broadband water emissivity as 0.9907 and compute emissions for their best estimated surface emissivity versus unity. Differences are up to 6 W/m2 in deserts, and can exceed 1.5 W/m2 in barren areas and shrublands.

So there is potential variation of a few W/m2 depending on the approach, and Trenberth et al settles on 396 W/m2 average – at least the values can be calculated, whereas our duo decided it was computationally impossible – perhaps as they saw the problem as requiring a totally accurate GCM.

With this information, the radiative balance problem can be resolved and we can see that there is a discrepancy between the solar energy absorbed and the terrestrial radiation emitted which requires explanation. The inappropriately-named “greenhouse effect”.

Without this information we can delight in much maths and pretend that nothing can be known about anything.

Why Conduction Can be Safely Ignored and Why We Just Demonstrated It

In many climatological texts it seems to be implicated that thermal radiation does not need to be taken into account when dealing with heat conduction, which is incorrect. Rather, always the entire heat flow density q must be taken into account..  It is inadmissible to separate the radiation transfer from the heat conduction, when balances are computed..

Unfortunately, the work on even the simplest examples of heat conduction problems needs techniques of mathematical physics, which are far beyond the undergraduate level.

In fact in many texts on atmospheric physics conduction is safely ignored due to the very low value of heat conduction through gases. Strictly speaking, if we write an equation then all terms should be included, including latent heat and convection. Why just radiation and conduction?

As Ramanathan and Coakley pointed out in their 1978 paper, convection is what determines the temperature gradient of the atmosphere but solving the equations for convection is a significant problem – so the radiative convective approach is to use the known temperature profile in the lower atmosphere to solve the radiative transfer equations.

Still, no thought of conduction as that term is so insignificant – as our intrepid duo go on to realize..

Commenting on the insolubility of heat flow via conduction they take a “typical example”:

If the radius of the Moon were used as the characteristic length and typical values for the other variables, the relaxation time would be equivalent to many times the age of the universe.

Therefore, an average ground temperature (over hundreds of years) is no indicator at all that the total irradiated solar energy is emitted. If there were a difference, it would be impossible to measure it, due to the large relaxation times. At long relaxation times, the heat flow from the Earth’s core is an important factor for the long term reactions of the average ground temperature; after all, according to certain hypotheses the surfaces of the planetary bodies are supposed to have been very hot and to have cooled down. These temperature changes can never be separated experimentally from those, which were caused by solar radiation.

So heat flow by conduction is so low that achieving balance by this method will take more than the age of the universe. Therefore, it is insignificant in comparison with convection and radiation.

Good so we can move on and climate scientists are right to ignore it. Was that the point that Gerlich and Tscheuschner were making? Yes, although possibly without realizing it..

Finally, the Imaginary Second Law of Thermodynamics

In their almost concluding section we see where countless climate enthusiasts have obtained their knowledge (or the reverse).

First, here’s an extract from a contemporary work on thermodynamics. This is from Fundamentals of Heat and Mass Transfer, 6th edition (2007), by Incropera & Dewitt:

As can be seen in the text, radiation can be absorbed by a higher temperature surface from a lower temperature surface and vice versa. Of course, the net result is a heat transfer from the hotter to the cooler.

The same uncontroversial description can be found in any standard thermodynamics work, unless they consider it too unimportant to mention. Certainly, none will have a warning sign up saying “this doesn’t happen”.

The explanation of the “greenhouse effect” is that the earth’s surface warms the lower atmosphere by radiation (as well as convection and latent heat transfer). And the atmosphere in turn radiates energy in all directions – one of which is back to the earth’s surface. Believers in the imaginary second law of thermodynamics don’t think this can happen. And this is possibly due to the miseducation by our intrepid duo. Or perhaps they learnt their thermodynamics from many “climate science” blogs.

The result of the actual climate situation is that the earth’s surface is warmer than it would have been without this atmospheric radiation. Pretty simple in concept.

Here’s how Gerlich and Tscheuschner explain things:

Everyone agrees.

Now the confusion. What are they saying? This isn’t what atmospheric physicists describe. The net heat transfer is from the earth’s surface (which was warmed by the sun) to the atmosphere.

Are they saying that it is impossible for any radiation to transfer heat from the atmosphere to the earth? It would appear so –

Following their diagram above, they comment, first quoting Rahmstorf:

Some `sceptics’ state that the greenhouse effect cannot work since (according to the second law of thermodynamics) no radiative energy can be transferred from a colder body (the atmosphere) to a warmer one (the surface). However, the second law is not violated by the greenhouse effect, of course, since, during the radiative exchange, in both directions the net energy flows from the warmth to the cold.

Rahmstorf’s reference to the second law of thermodynamics is plainly wrong. The second law is a statement about heat, not about energy. Furthermore the author introduces an obscure notion of “net energy flow”. The relevant quantity is the “net heat flow”, which, of course, is the sum of the upward and the downward heat flow within a fixed system, here the atmospheric system. It is inadmissible to apply the second law for the upward and downward heat separately redefining the thermodynamic system on the fly.

Our duo first attempt to confuse, as they frequently do in their opus by claiming that a clear explanation is obscure because precise enough terms aren’t used. It’s not obscure because they make the “correction” themselves.

Then add their masterstroke. It is inadmissible to apply the second law for the upward and downward heat separately redefining the thermodynamic system on the fly.

What on earth do they mean? Our comedic duo are the ones separating the system into upward and downward heat, followed by an enthusiastic army on the internet. Everyone else considers net heat flow.

As we saw in a standard work on thermodynamics, now in its 6th edition after two or three decades in print, there is no scientific problem with radiation from a colder to a hotter body – so long as there is a higher radiation from the hotter to the colder.

At this point I wonder – should I revisit the library and scan in 20 thermodynamic works? 50? What would it take to convince those who have been miseducated by our intrepid duo?

Perhaps Gerlich and Tscheuschner can now turn their attention to all of the unscientific text books like the one shown at the start of this section..

Conclusion

There is much to admire in Gerlich and Tscheuschner’s work. It can surely become a new standard for miseducation and we can expect its deconstruction by psychologists and those who study theories of learning.

From a scientific point of view, there is less to admire.

They have no understanding of modern climate science, content to dwell on works from over 100 years ago and ignoring any modern work. They appear to believe that the basis for the “greenhouse” effect is an actual greenhouse (as was covered in On Having a Laugh) even though no serious work on the subject relies on greenhouses. (Some don’t even mention it, some mention it to point out that the atmosphere doesn’t really work like a greenhouse).

In fact, the serious work of the last few decades relies on the radiative transfer equations – equations apparently endorsed by our duo, although their comments are “obscure”.

They take many other snipes at climate science by the approach of pointing out a term or dependency has been “neglected” (for example, like conduction through the atmosphere) without showing that the neglect has a significant impact – except in the case of conduction where (unwittingly?) they appear to show that conduction should definitely be ignored!

Someone could take issue with even modern work on climate science by the fact that they ignore relativistic effects.

Within the frame of modern physics, climate science is badly flawed to ignore relativity

And after 18 pages of unnecessary re-derivation of general relativity we find that “it’s therefore impossible to calculate this and the problem is insoluble“..

Well, although they haven’t read any modern climate science, it’s hard to see how they could be so confused about the application of the 2nd law of thermodynamics.

Perhaps in their follow up work they can explain why all the thermodynamics works are wrong, and especially where this 15μm (longwave) radiation comes from:

Measured downward longwave radiation at the earth's surface

Measured downward longwave radiation at the earth's surface

According to their interpretation of the 2nd law of thermodynamics this can’t happen. No heat can flow from the colder atmosphere to the warmer surface as that would be a “perpetuum mobile” and therefore impossible.

Where is it coming from Gerlich and Tscheuschner?

Read Full Post »

There are many misconceptions about how atmospheric processes work, and one that often seems to present a mental barrier is the idea of How much work can one molecule do?

This idea – presented in many ways – has been a regular occurence in comments here and it also appears in many blogs with eloquent essays on the “real role” of CO2 in the atmosphere, usually unencumbered by any actual knowledge of the scientific discipline known as physics.

Well, we all need mental images of how invisible or microscopic stuff really works.

When we consider CO2 (or any trace gas) absorbing longwave radiation the mental picture is first of trying to find a needle in a haystack.

And second, we found it, but it’s so tiny and insignificant it can’t possibly do all this work itself?

How much can one man or woman really do?

This article is really about the second mental picture, but a quick concept for the first mental picture for new readers of this blog..

Finding a Needle in a Haystack

Think of a beam of energy around 15.5μm. Here is the graph of CO2 absorption around this wavelength. It’s a linear plot so as not to confuse people less familiar with log plots. Water vapor is also plotted on this graph but you can’t see it because the absorption ability of water vapor in this band is so much lower than CO2.

CO2 absorption, 15.4-15.6um, linear, from spectralcalc.com

CO2 absorption, 15.4-15.6um, linear, from spectralcalc.com

The vertical axis down the side has some meaning but just think of it for now as a relative measure of how effective CO2 is at each specific wavelength.

Here’s the log plot of both water vapor and CO2. You can see some black vertical lines – water vapor – further down in the graph. Remember as you move down each black horizontal grid line on the graph the absorption ability is dropping by a factor of 100. Move down two black grid lines and the absorption ability has dropped by a factor of 10,000.

CO2 absorption - log graph - 15.4-15.6um, from spectralcalc.com

CO2 absorption - log graph - 15.4-15.6um, from spectralcalc.com

Now, I’ll add in the absorption ability of O2 and N2 – the gases that make up most of the atmosphere – check out the difference:

O2 and N2 added..

O2 and N2 added..

Spectralcalc wouldn’t churn anything out – nothing in the database.

15.5μm photons go right through O2 and N2 as if they didn’t exist. They are transparent at this wavelength.

So, on our needle in the haystack idea, picture a field – a very very long field. The haystacks are just one after the other going on for miles. Each haystack has one needle. You crouch down and look along the line of sight of all these haystacks – of course you can only see the hay right in front of you in the first one.

Some magic happens and suddenly you can see through hay.

Picture it..  Hay is now invisible.

Will you be able to see any needles?

That’s the world of a 15.5μm photon travelling up through the atmosphere. Even though CO2 is only 380ppm, or around 0.04% of the atmosphere, CO2 is all that exists for this photon and the chances of this 15.5μm photon being absorbed by a CO2 molecule, before leaving this world for a better place, is quite high.

In fact, there is a mathematical equation which tells us exactly the proportion of radiation of any wavelength being absorbed, but we’ll stay away from maths in this post. You can see the equation in CO2 – An Insignificant Trace Gas? Part Three. And if you see any “analysis” of the effectiveness of CO2 or any trace gas which concludes it’s insignificant, but doesn’t mention this equation, you will know that it is more of a poem than science. Nothing wrong with a bit of poetry, if it’s well written..

Anyway, it’s just a mental picture I wanted to create. It’s not a perfect mental picture and it’s just an analogy – a poem, if you will. If you want real science, check out the CO2 – An Insignificant Trace Gas Series.

CO2 – The Stakhanovite of the Atmospheric World?

Back in the heady days of Stalinist Russia a mythological figure was created (like most myths, probably from some grain of truth) when Aleksei Stakhanov allegedly mined 14 times his quote of coal in one shift. And so the rest of the workforce was called upon to make his or her real contribution to the movement. To become Stakhanovites.

This appears to be the picture of the atmospheric gases.

Most molecules are just hanging around doing little, perhaps like working for the _____ (mentally insert name of least favorite and laziest organization but don’t share – we try not to offend people here, except for poor science)

So there’s a large organization with little being done, and now we bring in the Stakhanovites – these champions of the work ethic. Well, even if they do 14x or 100x the work of their colleagues, how can it really make much difference?

After all, they only make up 0.04% of the workforce.

But this is not what the real atmosphere is like..

Let’s try and explain how the atmosphere really works, and to aid that process..

A Thought Experiment

For everyone thinking, “there’s only so much one molecule can do”, let’s consider a small “parcel” of the atmosphere at 0°C.

We shine 15.5μm radiation through this parcel of the atmosphere and gradually wind up the intensity. Because it’s a thought experiment all of the molecules involved just stay around and don’t drift off downwind.

The CO2 molecules are absorbing energy – more and more. The O2 and N2 molecules are just ignoring it, they don’t know why the CO2 molecules are getting so worked up.

What is your mental picture? What’s happening with these CO2 molecules?

a) they are just getting hotter and hotter? So the O2 and N2 molecules are still at 0°C and CO2 is at first 10°C, then 100°C, then 1000°C?

b) they get to a certain temperature and just put up a “time out” signal so the photons “back off”?

c) other suggestions?

The Real Atmosphere – From Each According to His Ability, To Each According to His Need

What is the everyday life of a molecule like?

It very much depends on temperature. The absolute temperature of a molecule (in K) is proportional to the kinetic energy of the molecule. Kinetic energy is all about speed and mass. Molecules zing around very fast if they are at any typical atmospheric temperature.

Here’s a nice illustration of the idea (from http://www.chem.ufl.edu/~itl/2045/lectures/lec_d.html).

At sea level, a typical molecule will experience around 1010 (10 billion) collisions with other molecules every second. The numbers vary with temperature and molecule.

Think of another way – at sea level 8×1023 molecules hit every cm2 of surface per second.

Every time molecules collide they effectively “share” energy.

Therefore, if a CO2 molecule starts getting a huge amount of energy from photons that “hit the spot” (are the right wavelength) then it will heat up, move even faster, and before it’s had time to say “¤” it will have collided with other molecules and shared out its energy.

This section of the atmosphere heats up together. CO2 can keep absorbing energy all day long even as a tiny proportion of the molecular population. It takes in the energy and it shares the energy.

If we can calculate how much energy CO2 absorbs in a given volume of the atmosphere we know that will be the energy absorbed by that whole volume of atmosphere. And therefore we can apply other well-known principles:

  • heating rates will be determined by the specific heat capacity of that whole volume of atmosphere
  • re-radiation of energy will be determined by the new temperature and ability of each molecule to radiate energy at wavelengths corresponding to those temperatures

Conclusion

The ability of a CO2 molecule to be “effective” in the atmosphere isn’t dependent on its specific heat capacity.

Molecules have embraced “communism” – they share totally, and extremely quickly.

Update – New post on the related topic of understanding the various heat transfer components at the earth’s surface – Sensible Heat, Latent Heat and Radiation

Read Full Post »

One commenter asked about CO2 absorption in the solar spectrum.

If CO2 absorbs incoming solar radiation then surely an increase in CO2 will reduce incoming radiation and balance any increase in longwave radiation.

The important factor is the usual question of quantifying the different effects.

Let’s take a look.

CO2 absorption in the 0.17-5um band, with solar spectrum overlaid

CO2 absorption in the 0.17-5um band, with solar spectrum overlaid

The CO2 absorption spectrum is from the line list browser of the recommended spectralcalc.com. The line list only goes down to 0.17μm (170nm), hence the reason for the graph not starting at 0.0μm.

The solar radiation is overlaid. Well, more accurately, the Planck function for 5780K is overlaid (simply drawn using Excel). Note that the CO2 absorption spectrum is on a log graph, while the radiation is on a linear graph. For those not so familiar with logarithmic graphs, the peak absorption around 4.3μm is 10-18, while the two peak absorptions just below 1μm are at 10-26 – which is 100,000,000 less.

The value of seeing the solar radiation spectrum overlaid is it enables you to see the relative importance of each absorption area of CO2. For example, the solar radiation between 2 – 4μm is only 5% of the solar radiation, so any absorption by CO2 will be quite limited.

Here’s the comparison with the important 15μm band of CO2. A 6μm width is shown, overlaid (blue line) with the 12-18um longwave radiation of a 288K (15°C) blackbody:

CO2 absorption in the 12-18um band, with terrestrial spectrum overlaid

CO2 absorption in the 12-18um band, with terrestrial spectrum overlaid

Just a little explanation of this graph and how to compare it to the solar version.

The average surface temperature of the earth is 15ºC, and it emits radiation very close to blackbody radiation (watch out for a dull post on Emissivity soon).

The proportion of radiation of a 288K blackbody between 12-18μm is 28%. What we want to do is enable a comparison between the CO2 absorption of solar radiation and terrestrial radiation.

Averaged across the globe and the year the incoming solar radiation at the top of atmosphere (TOA) is 239 W/m2 and the radiation from the earth’s surface is 396 W/m2. This works out to 65% higher, but so as not to upset people who don’t quite believe the earth’s surface radiation is higher than incoming solar radiation I simply assumed they were equal and scaled this section of the earth’s terrestrial radiation to about 28% of the solar radiation on the earlier graph. We are only eyeballing the two graphs anyway.

So with this information digested, the way to compare the two graphs is to think about the absorption spectra of CO2 simply being scaled by the amount of radiation shown overlaid in both cases.

As you can see the amount of absorption by CO2 of solar radiation is a lot less than the absorption of longwave radiation. Remember that we are looking at the log plot of absorption.

Is That the Complete Story?

Really, it’s more complicated, as always with atmospheric physics. There’s nothing wrong with taking a look at the approximate difference between the two absorption spectra, but luckily someone’s already done some heavy lifting with the complete solution to the radiative transfer equations using line by line calculations. For more on these equations, see the CO2 – An Insignificant Trace Gas series, especially Part Three, Four and Five.

The paper with the heavy lifting is Radiative forcing by well-mixed greenhouse gases: Estimates from climate models in the IPCC AR4 by W.D. Collins (2006). There’s a lot in this paper and aspects of it will show up in the long awaited Part Eight of the CO2 series and also in Models, On – and Off, the Catwalk.

Solving these equations is important because we can look at the absorption spectrum of CO2 in the 15μm band, but then we have to think about the absorption already taking place and what change in absorption we can expect from more CO2. Likewise for the solar spectrum.

Here are the two graphs, which include other important trace gases, as well as the impact of a change in water vapor. Note the difference in vertical axis values – the forcing effect of these gases on solar radiation has to be multiplied by a factor of 1000 to show up on the graph. The blue lines are CO2.

Net absorption of solar radiation by various "greenhouse" gases

Net absorption of solar radiation by various "greenhouse" gases

Longwave radiative forcing from increases in various "greenhouse" gases

Longwave radiative forcing from increases in various "greenhouse" gases

You can also see that the CO2 absorption in shortwave is across quite narrow bands (as well as being scaled a lot lower than terrestrial radiation) – therefore the total energy is less again. The vertical scale is energy per μm..

From these calculations we can see that with a doubling of CO2 there will be a very small impact on the radiation received at the surface, but a comparatively huge increase in longwave radiation retained – “radiative forcing” at the tropopause (the top of the troposphere at 200mbar).

So Is That the Complete Story?

Not quite. If trace gases in the atmosphere absorb solar radiation, is that so different from the surface absorbing solar radiation?

Or to put it another way, if the radiation doesn’t strike the ground, where does it go? It’s still absorbed into the climate system, but in a different location (somewhere in the atmosphere).

But as one commenter said:

The other point [this one] you make is simply not true and/or also not proven. There is only so much energy that can be taken up by a molecule.

This is a theme that has arrived in various comments from various posts. So the concept of How much work can one molecule do? is worth exploring in a separate post.

Hopefully, it’s clear from what is presented here that increases in CO2 absorption of the solar radiation are very small compared with absorption of longwave radiation.

Read Full Post »

Recap

This post is a follow on from my original article: American Thinker – the Difference between a Smoking Gun and a Science Paper.

Gary Thompson who wrote the article in American Thinker that prompted my article was kind enough to make some comments and clarifications. Instead of responding in the comments to the first article, I thought it was worth a new post, especially to get the extra formatting that is allowed in posts rather than comments.

I appreciate him commenting and I’m sure he has a thick skin, but just in case, any criticisms are aimed at getting to the heart of the matter, rather than at him – or anyone else who has different views.

For people who have landed at this post and haven’t read the original.. the heart of Gary’s article were 3 papers, of which I examined one (the first). The paper compared two 3-month periods 27 years apart in the East and West Pacific. Gary commented that the actual OLR (outgoing longwave radiation) was higher in the later period in the important CO2 band (or what we could see of it).

His claim – the theory says that more CO2 should lead to less emission at those wavelengths and therefore the theory has been disproved.

My point – Gary doesn’t understand the theory. The temperature was higher in the later period in this region and therefore the radiation leaving the earth’s surface would be higher. We’ll see this explained again, but did I mention that you should read Gary’s article and my article before moving forward? Also at the end of the Science of Doom post you can see Gary’s comments. Always worth reading what people actually wrote rather than what someone (me) with the opposite point of view highlighted from their words..

The Unqualified Statements in Papers

Gary started by saying:

I know why the authors of the papers were using climate models to simulate the removal of effect from surface temperatures and humidity and that the ‘theory’ says you must do that. But my problem lies in two peer reviewed papers that casts doubt on that theory and that method.

And cites two papers. The first, a 1998 paper: The Trace Gas Greenhouse Effect and Global Warming by the great V. Ramanathan (I will continue to call him ‘great’ even though he didn’t reply to my email about his 1978 paper.. possibly busy, but still..).

I recommend this 1998 paper to everyone reading this article. Even though it is 12 years old, it is all relevant and a very readable summary.

The Great Ramanathan

Gary pointed out page 3 where statements appeared to back up his (Gary’s) interpretation of the later OLR study. Here’s what Ramanathan said:

Why does the presence of gases reduce OLR? These gases absorb the longwave radiation emitted by the surface of the earth and re-emit to space at the colder atmospheric temperatures. Since the emission increases with temperature, the absorbed energy is much larger than the emitted energy, leading to a net trapping of longwave photons in the atmosphere. The fundamental cause for this trapping is that the surface is warmer than the atmosphere; by the same reasoning decrease of temperature with altitude also contributes to the trapping since radiation emitted by the warmer lower layers are trapped in the regions above.

By deduction.. an increase in a greenhouse gas such as CO2 will lead to a further reduction in OLR. If the solar absorption remains the same, there will be a net heating of the planet.

Gary commented on the last part of this:

Notice there is no clarifying statement about having to use model simulated graphs to ‘correct’ for surface temperatures and water vapor before seeing that OLR reduction.

And on the first part:

“since the emission increases with temperature, the absorbed energy is much larger than the emitted energy, leading to a net trapping of longwave photons in the atmosphere.” – here the author stated clearly that even taking into account higher emissions from warmer surfaces, the net will still be a reduction.

For half the readers here, they are shaking their heads.. But for Gary and the other half of the readership, let’s press on.

First of all, if we took any of 1000 papers on the “greenhouse” effect and observations, models, theoretical adjustments, impacts on GCMs – I bet you could find at least 700 – 900 of them at some point will make a statement that could be pulled out which has no “clarifying statement”. That “cast doubt” on the theory. Perhaps 1000 out of a 1000.

Context, context, context as they say in real estate.

Where to begin? Let’s look at “the theory” first. And then come back and examine Ramanathan’s statements.

The Theory

There are a few basics. For newcomers, you can take an extended look at the theory in CO2 – An Insignificant Trace Gas? It’s in seven parts! Actually it’s a compressed treatment.

This is itself is a clue.

In the books I have seen on Atmospheric Physics many tens of pages are devoted to radiation, including absorption and re-radiation – the “greenhouse” effect – and many tens of pages are devoted to convection. Understanding the basics is critical.

For Gary and his followers, the theory is on a precipice and these papers are giving us that clue. For people who’ve studied the subject the theory of the “greenhouse” effect is as solid as the theory of angular momentum, or the 2nd law of thermodynamics (the real one, not the imaginary one).

As Ramanthan says in the same paper Gary cites:

It is convenient to separate the greenhouse effect from global warming. The former is based on observations and physical laws such as Planck’s law for black-body emission. The concept of warming that results from the greenhouse effect, is based on deductions from sound physical principles. Numerous feedback processes, determine the magnitude of the warming; these feedbacks are treated with varying degrees of sophistication in GCMs and other climate models. As a result, predictions of the magnitude of the warming are not only model dependent but are subject to large uncertainties..

If I can paraphrase:

Greenhouse effect – dead solid. Global warming – lots of factors, need GCMs, pretty complicated.

Perhaps he has not realized that his words in the same paper combined with experiments demonstrate a flaw in the theory of the dead-solid “greenhouse” effect..

Back to the theory, but first..

A Quick, possibly Annoying, Diversion to the Theory of Gravity

But before we start, I thought it was worth an analogy. Analogies are illustrations, not proof of anything. They can often inflame an argument, but that’s not the intention here. Many people who are still undecided about the amazing theory of the inappropriately-named “greenhouse” effect might welcome a break from thinking about it.

The theory of planetary movements is my analogy. I picture a world where gravity – and its effects on the planets orbiting the sun – is strangely controversial. A concerned citizen, leafing through some fairly recent scientific papers notices that planets don’t really go around the sun in ellipses as the theory claims. In fact, there are some quite odd movements. And so, surprised that the scientists can’t see what’s in plain sight, this citizen draws attention to them.

When some detail-orientated commenters point out that the theory is actually (in part):

F = GMm/r2

where F= force between 2 bodies, G is a constant, M and m are the masses of the 2 bodies and r is the distance between them

And the ellipse idea is just a handy generalization of the results of the laws of graviational attraction.

And the reason why some planetary movements recently measured don’t follow an ellipse is because a few planets are a bit closer together, there’s a large asteroid flying between them and so when we do the maths it all works out pretty well.

The original concerned citizen then pulls out a few papers where, in the introduction, gravity is explained as that force that produces elliptical movements in the planets.. with no disclaimers about F=GMm/r2 and claims that this theory is, therefore, under question.

Why this annoying analogy? Most “theories” in physics are developed because of some observations which get analyzed to death – and finally someone produces a “comprehensive theory” that satisfies most of the relevant scientific community. The paper with the comprehensive theory usually contains some equations, some observations, some matching of the two – but often in common everyday usage the shorthand version of the theory is used.

“The theory of gravity tells us that planets orbit the sun in an elliptical manner, and ..”

So much more tedious to keep saying F=GMm/r2, and the other formulae..

(I know, I haven’t proved anything, but maybe a few readers can take a moment and see a parallel..)

Back to the Radiative-Convective Theory

First, bodies radiate energy according to Planck’s formula, which looks complicated when written down. The idea is simplified by the total energy radiated according to the Stefan-Boltzmann law which says that total energy is proportional to the 4th power of temperature.

Temperature goes up, energy radiated goes up (quite a bit more) – and the peak wavelength is a little lower. Here’s a graphical look of the Planck formula which takes away the mathematical pain:

Blackbody Radiation at 288K and 289K (15'C and 16'C)

Blackbody Radiation at 288K and 289K (15'C and 16'C)

Two curves – 288K and 289K. Now zoomed in a little where most of the energy is:

Close up of the peak energy of 288K and 289K

Close up of the peak energy of 288K and 289K

Total energy in the top curve (289K) is 396W/m^2, and in the bottom curve (288K) is 390W/m^2

Second, trace gases absorb energy according to a formula, which when simplifed is the Beer-Lambert law.

Absorption of Radiation as "optical thickness" increases, Iz=I0.exp (-x)

Absorption of Radiation as "optical thickness" increases

Without going into a lot of maths, as the concentration of a “greenhouse” gas increases, the absorption graph falls off more steeply – more energy is absorbed.

Third, re-radiation of this energy takes place according to an energy balance equation that you can see in CO2 – Part Three. The energy balance or radiative transfer equations rely on knowledge of the temperature profile in the lower atmosphere (the troposphere), because the actual temperature here is dominated by convection not radiation. (Radiation still takes place and the temperature profile is a major factor in the radiation – but this effect is not the primary determinant of the temperature profile).

When Gary says:

I know why the authors of the papers were using climate models to simulate the removal of effect from surface temperatures and humidity and that the ‘theory’ says you must do that. But my problem lies in two peer reviewed papers that casts doubt on that theory and that method.

It sounds like Gary believes this “model” is some suspicious extra that tries to deal with problems between the theory and the real world. But it’s the foundation. Anyone who had read an introduction to atmospheric physics would understand that. Someone who had tried hard to understand a few papers without a proper foundation would easily miss it.

  • Higher temperatures increase OLR
  • More trace gases reduce OLR in certain wavelengths

Which effect dominates in a particular situation?

It’s simple conceptually. But, if we want to find out exact results – such as, which effect dominates in a particular situation – we need a “model” = an equation or set of equations. Because if we want to quantify the effects we have to solve some tricky equations which you can see in a paper by.. the Great Ramanathan. Well, Ramanathan & Coakley 1978  – a seminal paper on Climate Modeling through Radiative-Convective Models, here’s an extract from p7:

A few equations from Ramanathan and Coakley, p7

A few equations from Ramanathan and Coakley, p7

Lots of maths. The rest of the paper is similar. Let’s move on to Ramanathan’s much later paper and what he said and meant.

Has Ramanathan given up on The Theory?

Let’s review the words Gary pulled out of the paper:

Why does the presence of gases reduce OLR? These gases absorb the longwave radiation emitted by the surface of the earth and re-emit to space at the colder atmospheric temperatures. Since the emission increases with temperature, the absorbed energy is much larger than the emitted energy, leading to a net trapping of longwave photons in the atmosphere.

Gary reads into this “here the author stated clearly that even taking into account higher emissions from warmer surfaces, the net will still be a reduction“. By which Gary thinks Ramanathan is saying something like:

..measurements from a specific location at a later date when CO2 has increased will always lead to a reduction in OLR..

-my paraphrase.

But no, he has totally misunderstood what the author is saying. Ramanathan is doing a quick drive through of the basics and explaining how the greenhouse effect works.

Lower temperatures in the atmosphere mean that the radiation to space from this lower temperature atmosphere is lower than the radiation from the surface. (This is the “net trapping”). Therefore – the “greenhouse” effect. There is no conclusion here that “at all times in all situations increases in “greenhouse” gases will lead to a reduction in OLR in these bands“. The conclusion is just that “greenhouse” gases mean that the surface is warmer than it would be without these gases

And the first claim, Ramanathan said:

By deduction.. an increase in a greenhouse gas such as CO2 will lead to a further reduction in OLR. If the solar absorption remains the same, there will be a net heating of the planet.

Gary said:

Notice there is no clarifying statement about having to use model simulated graphs to ‘correct’ for surface temperatures and water vapor before seeing that OLR reduction.

That’s because the basics of the theory are the solid foundation and don’t need to be restated as qualifiers to every statement. Ramanathan helped write the theory! For people who think that this stuff is just some added extra, read his 1978 paper and see all the maths and the explanations. This is the theory.

In the earlier part of his earlier statement he said “Since the emission increases with temperature” but didn’t qualify it with “emission increases in proportion to the 4th power of temperature“.

Is Ramanathan losing confidence in the Stefan-Boltzmann formula? Or Planck? Are these rocks crumbling?

It’s only because Gary has decided that these points are somehow rocky that he would reach these conclusions. (And I could pick 10’s of other statements in the paper which, without qualifiers, could be taken to be the beginning of the end of a specific theory).

As a general point – when we look at the hypothetical global annual average after “new equilibrium” from increased CO2 is reached – if this new equilibrium exists – the theory (1st law of thermodynamics) predicts that the “new” OLR will match the “old” OLR (global annual average). And in that case the OLR in “greenhouse” gas bands will be reduced a little, while the OLR outside of those bands will be increased a little.

But when we look at one local situation we need the theory, also known as “the model”, also known as equations, to tell us what exactly the result will be.

Ramanathan hasn’t cast any doubt on it. He believes it. His papers from the 1970s to today work it all out. He just doesn’t write qualifiers to each statement as if every line will be read by people who don’t understand the theory..

Gary’s Maths

In Gary’s comment he reproduced his back of envelope calculations of how much surface temperature should have changed over this 27 year period.

He takes the charitable approach of considering a net reduction in OLR from one of the three papers he originally reviewed – if I understood this step correctly. The figure he takes to work with is a 1K reduction. And then tries to work out how much this 1K reduction in OLR in the CO2 band would have on surface temperatures.

Like all good science this starts on napkins and the back of envelopes, because everyone studying a problem first has to attempt to quantify it using available data and available formulae. Then when the first results are worked out and it seems like something new is discovered – or something old overturned – then the scientist, patent clerk, writer now has to turn to more serious methods.

In Gary’s preliminary results he shows that a reduction in OLR due to CO2 might have contributed something like 0.3°C to surface temperature change over 30 years or so – where the GISS temperature increase for the period is something like 0.7°C. The essence of the calculation was comparing the Planck function (see the first and second graph in this post) of two temperatures 1K apart, then considering how much is in the CO2 band and so calculating the approximate change in W/m2. Then by applying a “climate sensitivity” number from realclimate.org, converting that into a temperature change.

Radiative physics has the potential to confuse everyone. Perhaps if we consider the “equilibrium case”, one problem with Gary’s calculation above will become clearer.

What’s the equilibrium case? In fact, it’s one generalized result of the real theory. And because it’s easier to understand than dI = -Inσdz + Bnσdz = (I – B)dχ people start to think the generalized result under equilibrium is the theory..

Under equilibrium, energy out = energy in. This is for the whole climate system. So we calculate the incoming solar radiation that is absorbed, averaged across the complete surface area of the earth = 239 W/m2.

So if the planet is not heating up or cooling down, energy out (OLR) must also = 239 W/m2.

So we consider the golden age of equilibrium in 1800 or thereabouts. We’ll assume the above numbers are true for that case. Then lots of CO2 (and other stuff) was added. Let’s suppose CO2 has reached the new current CO2 level of 380ppm – and stays constant from now on. The downward radiative forcing as calculated at “top of atmosphere” from the increase in CO2 is 1.7 W/m2. And nothing else happens in the climate (“all other things being equal” as we say).

Eventually we will reach the new equilibrium. Doesn’t matter how long it takes, but let’s pretend it’s 2020. Before equilibrium the planet will be heating up, which means OLR < 239 W/m2 (because more energy must come into the planet than leave the planet for it to warm up). At equilibrium, in 2020, OLR = 239 W/m2 – once again. (Of course, at wavelengths around 15μm the energy will be lower than in 1800 and other wavelengths the energy will be higher).

At this new equilibrium point we still have a “radiative forcing” of 1.7W/m2, which is why the surface temperature is higher, but no change in OLR when measured at 1800 and again at 2020.

We’ll assume, like Gary, the realclimate.org climate sensitivity – how do we calculate the new equilibrium temperature?

The change in OLR is 0.0 W/m2. Therefore, change in temperature = 0’C. Ha, we’ve proved realclimate wrong. Climate sensitivity cannot exist.

Or, it was wrongly applied.

As trace gases like CO2 increase in concentration they absorb energy. The atmosphere warms up and re-radiates the energy in all direction. Simplifying, we say it re-radiates both up and down. The extra downward radiation is what is used to work out changes in surface temperature. Not the immediate or eventual change in OLR.

The extra downward part is usually called “radiative forcing” and comes with a number of definitions you can see in Part Seven of the CO2 series (along with my mistaken attempt to do a “back of envelope” calculation of surface temperature changes without feedback).

How do we work out the change in surface temperature?

In Gary’s case, what he will need to know is the change in downward longwave radiation. Not upwards.

Conclusion

The theory of radiative transfer in the atmosphere, at its heart, is a relatively simple one – but in application is a very challenging one computationally speaking. Therefore, it’s hard to grasp it in its details intuitively.

The theory isn’t the generalized idea of what will happen when moving from one equilibrium state to another equilibrium state with more “greenhouse” gases. That is a consequence of the theory under specific (and idealized) circumstances.

But because everyone likes shorthand, to many people this has become “the theory”. So when someone applies the maths to a specific situation it is “suspicious”. Maths becomes equated with “models” which to many people means “GCMs” which means “make-belief”. I think that if Gary had read Climate Modeling through Radiative-Convective Models by the Great Ramanthan, he might have a different opinion about the later paper and whether Ramanathan is “bailing out” on the extremely solid theory of “greenhouse” gases.

Anyone who has read a book on atmospheric physics would know that Gary’s claim is – no nice way to say it – “unsupported”. Yes – cruel, harsh, wounding words – but it had to be said.

However, most of the readers here and most on American Thinker haven’t read an undergraduate book on the topic. So it makes it worth trying to explain in some detail.

It’s also great to see people trying to validate or falsify theories with a little maths. Working out some numbers is an essential step in proving or disproving our own theories and those of others. In the example Gary provided he didn’t apply the correct calculation. (The subsequent pages of Ramanathan’s 1998 paper also run through these basics).

For those convinced that the idea that more CO2 will warm the surface of the planet is some crazy theory – these words are in vain.

But for the many more people who want to understand climate science, hopefully this article provides some perspective. The claims in American Thinker might at first sight seem to be a major problem to an important theory. But they aren’t.

Note – I haven’t yet opened the Philipona paper, but will do so in coming weeks and probably add another article about it. I didn’t want to leave Gary’s comments unanswered, and this post is already long enough..

Note on a Few Technical Matters

A few clarifications are in order, for people who like to get their teeth into the technical details.

Strictly speaking, at the very top of atmosphere there is no downward longwave forcing at all. That’s because there’s no atmosphere to radiate.

The IPCC definition of “radiative forcing” at “top of atmosphere” is a handy comparison tool of extra downward radiation before feedbacks and before equilibrium of the surface or troposphere is reached. In reality the increase in downward radiation doesn’t occur in just one location, extra downward longwave radiation occurs all throughout the troposphere and stratosphere.

In the example above, from 2010 to 2020 as surface and troposphere temperatures increased, radiative forcing would also increase slightly so it wouldn’t necessarily be constant at 1.7W/m2.

Climate sensitivity is calculated by GCMs to work out the resulting long term (“equilibrium”) temperature change, with climate feedbacks, from increased radiative forcing. No warranty express or implied as to their accuracy or usefulness, just explaining how to apply the climate sensitivity value correctly.

Read Full Post »

In Part One, we introduced some climate model basics, including uses of climate models (not all of which are about “projecting” the future).

And we took at a look at them in their best light – on the catwalk, as it were.

Well, really, we took a look at the ensemble of climate models. We didn’t actually see a climate model at all..

Ensembles

The overall evaluation in Part One was the presentation of a “multi-model mean” or an ensemble. An ensemble can be the average of many models, or the average of one model run many times, or both combined.

We will return to more discussion about the curious nature of ensembles in a later post. Just as a starter, two observations from the IPCC.

IPCC AR4 in Chapter 8, Climate Models and their Evaluation, comments:

There is some evidence that the multi-model mean field is often in better agreement with observations than any of the fields simulated by the individual models (see Section 8.3.1.1.2), which supports continued reliance on a diversity of modelling approaches in projecting future climate change and provides some further interest in evaluating the multi-model mean results.

and a little later:

Why the multi-model mean field turns out to be closer to the observed than the fields in any of the individual models is the subject of ongoing research; a superficial explanation is that at each location and for each month, the model estimates tend to scatter around the correct value (more or less symmetrically), with no single model consistently closest to the observations. This, however, does not explain why the results should scatter in this way.

One interpretation of this would be:

We like ensembles because they give more accurate results, but we don’t really understand why..

A subject to come back to, now it’s time for a real model..

Step Forward Climate Model “Cici” – CCSM3

CCSM3, “Cici”, is the model from NCAR (National Center for Atmospheric Research) in the USA. Out of all the GCMs discussed in the IPCC AR4, Cici has the “best curves” – the highest resolution grid. Well, she comes from the prestigious NCAR..

The model’s vital statistics – first the atmosphere:

  • top of atmosphere = 2.2 hPa (=2.2mbar), this is pretty much the top of the stratosphere, around 50km
  • grid size = 1.4° x 1.4° (T85)
  • number of layers vertically = 26 (L26)

second, the oceans:

  • grid size = 0.3°–1° x 1°
  • number of vertical layers = 40 (L40)

The vital statistics give a quick indication of the level of resolution in the model. And there are also model components for sea ice and land. The model doesn’t need the infamous “flux adjustment” which is the balancing term for energy, momentum and water between the atmosphere and oceans required in most models to keep the two parts of the model working correctly.

The CCSM3 model is described in the paper: The Community Climate System Model Version 3 (CCSM3) by W.D. Collins et al, Journal of Climate (2006). The source code and information about the model is accessible at http://www.ccsm.ucar.edu/models/.

And for those who love equations, especially lots of vector calculus, take a look at the 220 page technical document on CAM3, the atmospheric component.

It will be surprising for many to learn that just about everything on this model is out in the open.

CCSM3 Off the Catwalk – Hindcast Results

As with the multi-model means results in Part One we will take a look through a similar set of results for CCSM3.

Annual temperature

CCSM3 Annual Land & Sea Temperature Actual (top) vs Model (bottom)

CCSM3 Annual Land & Sea Temperature Actual (top) vs Model (bottom)

Cici looks pretty good.

Details – The HadISST (Rayner et al., 2003) climatology of SST for 1980-1999 and the CRU (Jones et al., 1999) climatology of surface air tempeature over land for 1961–1990 are shown here. The model results are for the same period of the CMIP3 20th Century simulations. In the presence of sea ice, the SST is assumed to be at the approximate freezing point of sea water (–1.8 °C).

However, it’s hard to tell looking at two sets of absolute values, so of course we turn to the difference between model and reality.

Annual Temperature – Model Error

Model simulations of annual average temperature less observed values for Cici and for the “ensemble” or multi-model mean:

Annual Temperatures - Simulated minus observed for CCSM3 and the ensemble

Annual Temperatures - Simulated minus observed for CCSM3 and the ensemble

In terms of absolute error around the globe, Cici and the ensemble are very close (using the Anglotzen statistical method).

We could note that even though the values are “close”, there are areas where Cici – and the ensemble – don’t do so well. In Cici’s case southern Greenland and the Labrador Sea, which might be very important for predicting the future of the thermohaline circulation. And both are particularly bad for Antarctica, a general problem for models.

To give an idea of the variation of models, here are all of the models reviewed by the IPCC in AR4 (2007):

Annual Temperature - Model less Actual - All models

Annual Temperature - Simulated minus observed - All models

The top right is Cici (red circle). It’s clear that Cici is a supermodel..

Standard Deviation of Temperature

The standard deviation of temperature – “over the climatological monthly mean annual cycle ” – simulated less observed for Cici and the ensemble. We could describe it as how good is the model at working out how much temperature actually varies over the year in each location?

First, however, to make sense of the “error” of model less actual, we need to know what actual values look like:

Standard Deviation of Temperature over the climatological monthly mean annual cycle

Standard Deviation of Temperature over the climatological monthly mean annual cycle

As we would expect, the oceans show a lot less temperature variation than the land and around the tropics and sub-tropics the variation is close to zero.

Now let’s take a look at the model less actual, or “model error”:

Std Deviation of Temperature - Simulated minus observed for CCSM3 and the ensemble

Std Deviation of Temperature - Simulated minus observed for CCSM3 and the ensemble

We can see that Cici has some problems in modeling temperature variation especially under-estimating the actual variation around northern Russia and Canada and over-estimating the variation in the Middle East and Brazil. The ensemble appears to be in slightly better shape here.

Of course, these areas are where the largest temperature variation takes place.

Diurnal Range of Land Temperature

As before, first the actual values:

Annual average of diurnal temperature range over land

Annual average of diurnal temperature range over land

And now the model less actual, or “model error”:

Diurnal temperature range over land - Actual less Model for CCSM3 and ensemble

Diurnal temperature range over land - Actual less Model for CCSM3 and ensemble

We can see a lot of areas where the model error is quite large, usually corresponding to larger measured values. In the case of Greenland, for example, the annual average diurnal temperature range is over 20°C, while the model under-estimates this by more than 10°C. Given the legend the error might be as big as the actual value..

We can also see that on average Cici under-estimates the diurnal temperature range, and the ensemble is closer to neutral but still appears to under-estimate.

Here’s another comparison which demonstrates the problem of all the models vs observation:

Diurnal temperature range vs latitude - Observed compared with all models

Diurnal temperature range vs latitude - Observed compared with all models

The black line is the observed value. We can see that all of the models except for one are definitely under-estimating, and none of the models are particularly close to the observed values.

Now we can get to see more fundamental values.

Reflected Solar Radiation

This value is essential for calculating the basic radiation budget for the earth.

First the actual values as measured by ERBE (1985-1989):

Average Reflected Solar Radiation, ERBE

Average Reflected Solar Radiation, ERBE

And now the model error – model less actual:

Reflected Solar Radiation - Actual less Model for CCSM3 and ensemble

Reflected Solar Radiation - Actual less Model for CCSM3 and ensemble

The ensemble is definitely better than Cici. Cici has some large errors, for example, North Africa, Pacific Ocean and the Western Indian Ocean where the model error seems to be up to half of the actual value.

If we look at the values averaged by latitude the results appear a little better:

Reflected Solar Radiation vs latitude - Observed compared with all models

Reflected Solar Radiation vs latitude - Observed compared with all models

But the deviations give us a better view:

Reflected Solar Radiation vs latitude - Model error for all models

Reflected Solar Radiation vs latitude - Model error for all models

Note Cici in the solid blue line. The ensemble is proving to be the pick of the bunch..

So the model’s ability to simulate reflected solar radiation is much better by latitude than by location. But most or all of the models have significant discrepancies even when averaged over each latitude.

Outgoing Longwave Radiation

The other side of the radiation budget, first the actual ERBE measurement (1985-1989):

Outgoing Longwave Radiation, OLR, ERBE

Measured Outgoing Longwave Radiation, ERBE

And now the model error – model less actual:

OLR - Actual less Model for CCSM3 and ensemble

OLR - Actual less Model for CCSM3 and ensemble

As with reflected SW radiation, the ensemble performs better than Cici. So while measured values are in the range of 200-300 W/m2, Cici has some areas where the (absolute) error is in excess of 30W/m2.

Looking at the OLR values averaged by latitude, the results appear a little better:

OLR vs latitude - Observed compared with all models

OLR vs latitude - Observed compared with all models

And the deviations, or model error:

OLR vs latitude - Model error for all models

OLR vs latitude - Model error for all models

Rainfall

Measured from CMAP, 1980-1999:

Rainfall 1980-1999

Rainfall 1980-1999

Units are in cm of rainfall per year. And now the model error – model less actual:

Rainfall - Actual less Model for CCSM3 and ensemble

Rainfall - Actual less Model for CCSM3 and ensemble

Once again the ensemble outshines Cici. There are some substantial errors in the areas where rainfall is high.

As with some of the previous model results, if we look at the model vs observed by latitude the picture is somewhat better:

Rainfall vs latitude - Observed vs all models

Rainfall vs latitude - Observed vs all models

Humidity

Lastly, we will take a look at specific humidity. First the “measured”, as recalculated by ERA-40:

Specific Humidity vs Latitude and Altitude, from ERA40

Specific Humidity in g/kg vs Latitude and Altitude, from ERA40

Observed annual mean specific humidity in g/kg, averaged zonally, 1980-1999. Note that the vertical axis is pressure on the left in mbar and km in height on the right.

And now the model error – but this time in % = (model – actual)/actual x 100:

Specific Humidity - % Error for CCSM3 and ensemble

Specific Humidity - % Error for CCSM3 and ensemble

Once again the ensemble appears to outperform Cici. And both, but especially Cici, have problems in the top half of the troposphere (around 500-200mbar) with 20-50% error in some regions in Cici’s case.

Conclusion

This has been a quick survey of model results for different parameters across the globe, but averaged annually, compared with observations.

In Part One, we saw the ensemble in its best light. But when we take a look at a real model, the supermodel Cici, we can see that she has a lot of areas for improvement.

There’s lots more to investigate about models, all to come in future parts of this series.

As always, comments and questions are welcome, but remember the etiquette.

Read Full Post »

In the previous article in this series, The Earth’s Energy Budget – Part Two we looked at outgoing longwave radiation (OLR) and energy imbalance. At the end of the article I promised that we would look at problems of measuring things and albedo but much time has passed, promises have been forgotten and the fascinating subject of how the earth really radiates energy needs to be looked at.

If you are new to the idea of incoming (absorbed) solar radiation being balanced by OLR, or wonder how the solar “constant” of 1367W/m2 can be balanced by the earth’s OLR of 239W/m2 then take a look at Part One and Part Two.

Introduction

If you’ve read more in depth discussions about energy balance or CO2 “saturation” you might have read statements like:

More absorption by CO2 causes emission of radiation to move to higher, colder layers of the atmosphere

If these kind of comments confuse you, sound plain wrong, or cause you to furrow your brow because “it sounds like it’s probably right but what does it actually mean?” – well, hopefully some enlightenment can be found.

Effective Radiation

The sun’s core temperature is millions of degrees but we see a radiation from the sun that matches 5780K – its surface temperature:

Solar Radiation, top of atmosphere and at earth's surface, Taylor (2005)

Solar Radiation, top of atmosphere and at earth's surface, Taylor (2005)

In this figure there are two spectra: the top one is how the sun’s radiation looks before it reaches the top of the earth’s atmosphere – contrasted with the dotted line of a “blackbody” – or perfect radiator – at 5780K (5507°C for people new to Kelvin or absolute temperature).

The bottom one – of less interest for this article – is how the sun’s radiation looks at the earth’s surface after the atmosphere has absorbed at various wavelengths.

Why don’t we see a radiation spectrum from the sun that matches millions of degrees?

If we measure the upward longwave radiation from the earth’s surface at 15°C we see an effective “blackbody” radiator of 288K (15°C). But why don’t we see a radiation spectrum of 5000K – the temperature somewhere near the core?

The answer to both questions is that radiation from the hotter inner areas of these bodies gets completely absorbed by outer layers, which in turn heat up and radiate at lower temperatures. In the case of the sun, the radiation spectrum includes hotter areas below the surface that are not absorbed at some wavelengths, as well as the surface itself.

In the case of the earth it’s really the top skin layer that emits longwave radiation.

So when we measure the radiation from the earth with a surface temperature of 15°C (288K) we know we will see a longwave radiation that matches this 288K. This will be a total energy radiated of 390W/m2 with the peak wavelength of 10.1μm. The temperature below the surface is irrelevant.

(Well, it’s not really irrelevant. The hotter layers below warm up the layers above – through conduction and radiation).

This is what the radiation looks like:

Blackbody Radiation at 15'C or 288K

Blackbody Radiation at 15'C or 288K

This assumes an emissivity of 1. The emissivity of the surface of the earth varies slightly but is close to 1, typically around 0.98. Watch out for a dull post on emissivity at some stage..

At the top of atmosphere, as many know, the OLR is around 239W/m2. For those confused by how it can be 390W/m2 at the surface and 239W/m2 at the top, the answer is due to absorption and re-radiation of longwave radiation by trace gases – the “greenhouse” effect. See the CO2 – An Insignificant Trace Gas? series, and especially Part Six – Visualization and CO2 Can’t Have that Effect Because.. if you don’t understand or agree with these well-proven ideas.

If the earth’s atmosphere was completely transparent to longwave radiation this spectrum would look exactly the same at the earth’s surface and at the top of atmosphere (TOA).

Here’s what it does look like with some typical blackbody radiation curves overlaid:

Outgoing longwave radiation at TOA, Taylor (2005)

Outgoing longwave radiation at TOA, Taylor (2005)

(Note that the spectrum is shown in wavenumber in cm-1. For convenience I added wavelength in μm under the wavenumber axis. Wavelength in μm = 10,000/wavenumber).

For energy balance – if the earth is not warming up or cooling down – we would expect the earth to radiate out the same amount of energy that it absorbs from the sun. That amount is 239W/m2, which equates to an average temperature of 255K (-18°C).

As the text for this graphic shows, when the energy under the curve is integrated this is what it comes to! But as you can see the actual spectrum is not a “blackbody curve” for 255K. So let’s take a closer look.

Everything Gets Through or Nothing Gets Through – a Few Thought Experiments

Imagine a world where the upwards longwave radiation from the earth’s surface didn’t get absorbed by any gases in the atmosphere.

Most people are familiar with that thought experiment – it’s a staple of the most basic radiation model in climate science. The radiation at the top of the atmosphere would look like this (the top graph):

255K radiation, 100% transmittance

255K radiation, 100% transmittance

This is the blackbody radiation at 255K (-18°C) with 100% “transmittance” through the atmosphere. The area under the curve, if we extend it out to infinity, is 239W/m2.

And of course, because the radiation hadn’t been absorbed or attenuated in any way, the temperature at the earth’s surface would also be 255K. Chilly.

Now let’s think about what would happen if the atmosphere allowed radiation only through the “atmospheric window” and everywhere else the transmittance was zero:

323K radiation through a perfect "atmospheric window", 8-14um

323K radiation through a perfect "atmospheric window", 8-14um

The bottom graph shows how the transmittance of the atmosphere varies with wavelength in this thought experiment.

The top graph in this case is the blackbody radiation from 323K (50°C) only allowed through between 8-14μm. The energy under the curve is 239W/m2. (Note the higher values on the vertical scale compared with the earlier graphs).

So if the atmosphere absorbed all of the surface radiation below 8μm and above 14μm the earth’s surface would heat up until it reached 50°C (323K). Why? Because if the temperature was only 15°C the amount of energy radiated out would only be 141W/m2. More energy coming in than going out = earth heats up. The surface temperature would keep heating up until eventually 239W/m2 made it out through the atmospheric window – which is 50°C.

Closer to The Real World – Illustration of Radiation from Multiple Layers in the Atmosphere

Even in the atmospheric window some radiation is absorbed, i.e. the transmittance is not 1. But let’s assume for sake of argument it is 1. So energy in the 8-14μm band just passes straight through the atmosphere. It’s still a thought experiment.

Lots of gases absorb at lots of wavelengths – which makes thinking about it as a whole very difficult. So we’ll just assume that the rest of the atmosphere outside the atmospheric window all shares the same absorption characteristics – that is, every wavelength is identical in terms of absorption of radiation.

Now let’s try and consider what really happens in the atmosphere. Each “layer” of the atmosphere radiates out energy according to the temperature in that layer. For reference, here is the temperature (and pressure) at different heights:

Atmospheric Temperature & Pressure Profile, Bigg (2005)

Atmospheric Temperature & Pressure Profile, Bigg (2005)

The highlighted area at the bottom – the troposphere – is the area of interest. This is where most of the atmosphere (by molecules and mass) actually resides.

In our thought experiment radiation from the surface (outside the atmospheric window) gets completely absorbed by the atmosphere, or at least the amount that gets through is very small. Taken to the extreme we would get the result shown a few graphs earlier where the surface temperature rises up to 50’C.

But just because surface radiation doesn’t get out doesn’t mean that radiation from the atmosphere can’t get out.

Each layer of the atmosphere radiates according to its temperature. Even if the atmosphere’s transmittance is zero when considering the entire thickness of the atmosphere, there will be some layer where radiation starts to get through.

This is partly because there is less atmosphere to absorb the closer we get to the “top”. And also because as we get higher in the atmosphere it gets thinner. Less molecules to absorb radiation. Even if some gas is a fantastically good absorber of energy, there must be a point where radiation is hardly absorbed. For example, at the top of the stratosphere, about 50km, the pressure is around 1mbar – 1000x less than at the surface. At the top of the troposphere (the tropopause) the pressure is around 200mbar – 5x less than at the surface.

The challenge in thinking about the atmosphere radiating is that unlike the surface of the earth where all radiation is emitted from the very surface, instead radiation is emitted from lots of different layers:

  • Higher up – less absorption, more radiation makes it through
  • Lower down – more absorption, less radiation makes it through

But let’s still keep it simple and think about the surface temperature being our standard 15°C (288K) and the atmospheric window letting through everything between 8-14μm. This means 141W/m2 makes it out through this window.

If we have energy balance, the OLR = 239W/m2 in total = 141 (through the atmospheric window) + 98 (radiated from the atmosphere at some height, wavelengths outside 8-14μm).

What temperature equates to this layer in the atmosphere? Well, assuming no absorption above this radiating layer (not really the case), and only radiation outside 8-14μm, the temperature of the atmosphere would have to be 219K, or -54°C. Take a look back at the temperature profile above – this is pretty much the top of the troposphere, around 11km.

Remember that this isn’t exactly how radiation gets radiated out to space – it doesn’t come from one “skin layer”. We might consider that if the transmittance of the atmosphere is 1 at this height, then maybe at 10km the transmittance is 0.8 and at 9km the transmittance is 0.5, and at 8km the transmittance is 0.1..

So each layer is radiating energy, with higher layers being colder but more of their radiation getting through, and lower layers being warmer – so radiating a higher amount – but less of their radiation getting through.

For many people reading, this is a straightforward concept, why so long.. for others it might still seem tough to grasp..

So here is a sample radiation diagram with illustrative values only (and values a little different from above):

Radiation from different heights in the atmosphere, illustrative values only

Radiation from different heights in the atmosphere, illustrative values only

What the diagram shows is the radiation outside the 8-14μm band. That’s because in our thought experiment the 8-14μm band doesn’t absorb any radiation (and therefore can’t radiate in this band either).

Take the top layer at 11km. If we calculate the blackbody radiation of 219K (-54°C) and exclude radiation in the 8-14μm band the radiation is 98W/m2. Then the grey block above with “0.6” is the atmosphere above with transmittance of 0.6, so the radiation actually getting through from this layer to the top of atmosphere is 59W/m2. Similarly for the two other layers (with different values).

In total, the energy leaving the top of atmosphere (outside of the atmospheric window) is 98W/m2. (It’s just a coincidence that this is the value of the top layer before any absorption). And inside the atmospheric window was the number we already calculated of 141W/m2, so the total OLR is 239W/m2.

Of course, we all know the real atmosphere is much more complex with lots of different absorption at different wavelengths. But hopefully this “intermediate” example help to explain how the atmosphere radiates out energy.

So finally, onto the real point..

What Happens with More Absorbing Gases?

Remember how this long post started..

If you’ve read more in depth discussions about energy balance or CO2 “saturation” you might have read statements like:

More absorption by CO2 causes emission of radiation to move to higher, colder layers of the atmosphere

Now, maybe this kind of statement will make more sense.

In our model – our thought experiment – above, we had a uniform absorber of radiation outside the atmospheric window. Suppose we increase the amount of this absorber – the skies open and someone pours some more in and stirs it around. Let’s say the amount increases by 10%.

Well, take a look back at the last diagram. See the transmittance values for each layer in the atmosphere – 0.6 at 11km high, 0.25 at 10km high and 0.1 at 9km high.

Regardless of how realistic these actual numbers are, increasing the amount of absorbing gas by 10% will automatically mean that each of the transmittance numbers is reduced by 10%. And so less radiation makes it out to the TOA (top of atmosphere).

Effectively because lower layers are contributing less energy out through TOA the effective radiating height has moved up. It’s not because some directive has been passed down from a higher authority. And it’s not because one layer has stopped and another layer has taken over.

It’s just that lower layers contribute less, so the “average radiating height” is now higher and colder.

(Note: it might look at first sight that the average height is still the same even though the amount of radiation has reduced. This is not really the case, see the note at end).

In our particular example what would happen is that the OLR would reduce from 239W/m2 down to 141+98*0.9=229 W/m2. So the surface would warm up and this would warm up each layer of the atmosphere until eventually a new hotter steady state was reached.

Conclusion

This has been a long post to try and create more of an understanding of how the earth actually radiates energy, and why more of any trace gas increases the “greenhouse” effect.

It does it because more “absorbing gases” reduce the amount of radiation that can make it out from lower layers in the atmosphere. These lower layers are hotter and radiate much more energy. Proportionately more energy will then be radiated from higher layers which are colder, and therefore these radiate less energy.

It’s a not a mystical force that raises the “effective radiating height” in the atmosphere. But the effective radiating height does increase.

Note

In the example above, the three layers together contributed 98W/m2 at TOA. That is an “effective temperature” of 219K – remembering that we are excluding radiation from the 8-14μm window. If we reduce the radiation from these three layers by 10%, we now have 89W/m2 which is about 212K – effectively radiating from a colder level in the atmosphere.

Read Full Post »

Many readers of this blog would like progress towards the solution of the great questions in climate science. Other readers have stopped by still pondering the basics.

Some of those pondering the basics might have read many of the exciting claims on the internet that the “greenhouse” effect can’t exist because it would violate the 2nd law of thermodynamics.

It’s not a claim that you find in any books on atmospheric physics by the way, it’s strictly an “internet phenomenon”.

Just a few basics in case this is the first post you have read from this site.

The inappropriately named “greenhouse” effect can be summed up in a few sentences:

Longwave radiation from the earth’s surface is absorbed by many trace gases, including water vapor and CO2. The absorption causes these gases to heat up and energy is radiated back out – both up and down. The upward radiation is effectively “no change”. The downward radiation adds to the energy received from the sun and heats up the surface of the earth more than if this downward radiation did not occur.

If there was no absorption of radiation by “greenhouse” gases the surface of the earth would be a lot colder. Here is a very simplified graphic to draw people’s attention to the fact that “something big” is going on:

Upward Longwave Radiation, Numbers from Kiehl & Trenberth

Upward Longwave Radiation, Numbers from Kiehl & Trenberth (1997)

(TOA = top of atmosphere). If there was no absorption and re-radiation back down the two numbers would be the same. (Note that the downward radiation is not shown to crystallize the issue)

These numbers are global annual averages under a clear sky. Under a cloudy sky the numbers are different but similar – and still the radiation from the surface of the earth is a lot greater than that leaving through the top of atmosphere. For more on this take a look at CO2 – An Insignificant Trace Gas? – Part Six – Visualization and the followup CO2 Can’t Have That Effect Because.. as well as the start of the series on CO2.

Many people have said that the numbers are obviously wrong, I’m mixing up solar radiation and longwave radiation, it can’t happen, the temperature varies a lot from equator to poles so that’s why the radiation numbers are wrong..

As one person said on another blog, possibly commenting on one of these earlier posts:

I even saw one where the guy had 100 watts/m2 going into the atmosphere MORE than was coming out FOREVER.

Sharp-eyed readers might notice that I haven’t drawn in the downward radiation. Energy is balanced in the atmosphere because of the downward radiation from the atmosphere (not drawn). This is the “greenhouse” effect.

Onto the Imaginary Second Law of Thermodynamics

How can a colder atmosphere add heat to a warmer surface?

Can a candle warm the sun?

There are many popular restatements of the imaginary 2nd law. These two should be a representative sample. And so follows the Q.E.D. claim that the “greenhouse” effect plainly contradicts the second law of thermodynamics.

What is the second law?

The Real Second Law of Thermodynamics

My boring thermodynamics books and I have long since had a parting of the ways, so I looked it up on Wikipedia. Not a 100% reliable source, but the (real) second law is just as I remember it so I looked no further.

It’s possible that the imaginary second law has taken a strong hold because anyone who does look it up finds statements like dS/dt>=0, where S is entropy. Wow. Clever people. What’s entropy? How does this relate to candles? Candles can’t warm the sun, so I guess the second law has just proved the “greenhouse” effect wrong..

According to Wikipedia, Clausius expressed the second law (validly) like this:

Heat generally cannot flow spontaneously from a material at lower temperature to a material at higher temperature

Again, that seems right and it doesn’t have any entropy involved in the description. I never did like entropy. It never seemed real.

Perhaps this formulation has been the inspiration for the imaginary second law. As a not very precise definition many people might read this and think no energy at all can flow from a cold body to a hot body.

In fact, no net energy can flow from a cold body to a hot body.

In the case of the real “greenhouse” effect and the real 2nd law of thermodynamics, net energy is flowing from the earth to the atmosphere. But this doesn’t mean no energy can flow from the colder atmosphere to the warmer ground.

It simply means more energy flows from the warmer surface to the colder atmosphere than in the reverse direction.

Another likely reason the imaginary second law has become popular is most people are much more familiar with conduction of heat than radiation. Conduction of heat only appears to flow one way.

A Thought Experiment

We’ll do a thought experiment to demonstrate why the imaginary second law of thermodynamics is wrong. It’s simpler, safer, cheaper AND more reliable than assembling equipment. After all, we are going to look at radiation and if we do an experiment we would need to ensure that no convection or conduction was taking place.

And the thought experiment will, I hope, be more powerful. Plus it will have the added benefit for those already convinced by the beguiling imaginary second law that they can say “you haven’t proven anything, it’s all in your head” and so the popular imaginary law can live on.

In our thought experiment we will consider the sun. It’s hot. It doesn’t conduct or convect any heat outside its immediate surface because space is a vacuum and heat can only travel by radiation through a vacuum.

So energy is radiated out from the sun equally in all directions. At a 1000km distance from the sun, we get our measuring instrument out and find that energy radiated is 10,000W/m2 (because I can’t be bothered to work out the actual number).

Now we fly in a cold large rock and park it at 1000km from the sun. Energy from the sun is absorbed on this cold rock and it heats up to some equilibrium value where it is also radiating out what it is receiving.

The temperature of this large once-cold rock is now a toasty 648K (375’C). All is well with both the real and imaginary formulations of the 2nd law, so far.

Now, from a galaxy far far away, we fly in a new star. Before we started moving it we checked the radiation 1000km away from the star and found that it was 11,000W/m2. We are careful in our relocation of this star that nothing changes in its inner generation of radiation. The new star is parked 1000km away from the once-cold rock and 1000km away from the sun.

It’s a love triangle. Due to the new star’s welcome appearance, the rock heats up further. It now receives 21,000W/m2. Its new equilibrium temperature is 780K. All is still well with both the real and imaginary formulations of the 2nd law.

Trouble in Paradise

But now a problem.. the new star is radiating out in all directions. Believers in the imaginary second law have no problem with the idea that the sun receives energy from the new star. After all the sun was a little colder.

But what about the sun? It is also radiating out in all directions. Or it was before the new star arrived.

Now that the new star is parked 1000km from the sun, squarely in the path of some portion of the sun’s radiation we have to ask ourselves what actually happens?

Believers in the real second law of thermodynamics are quite happy. No cognitive dissonance there. The energy from the sun which is incident on the new star’s surface actually increases the new star’s surface temperature compared with what it was before.

11,000W/m2 are flowing from the new star to the sun, and 10,000W/m2 are flowing from the sun to the new star. Some kind of new equilibrium might be reached, but for real second law believers there is no angst. The net flow of energy is from the hotter to the colder.

Believers in the imaginary second law, what happens?

One obvious suggestion is that the sun’s paltry 10,000W/m2 which was flowing through that exact spot now divert around the new star as if it had some kind of force field. Perhaps all the energy lines completely redistribute so that (depending on the diameter of this new star) about 10,015W/m2 flow in all directions except through the location of the new star.

Another obvious suggestion is that the sun “realizes” the new star is there and energy is flowing from the new star to it so just stops radiating in that exact direction. I put “realizes” in quotes of course because we all know the sun is not sentient. It’s just terminology. Some process that drives the imaginary second law will no doubt make this happen.

And the most likely suggestion of all is that this radiation from the sun, when it strikes the surface of the new star, just bounces off. Or is absorbed but doesn’t actually heat up the surface of the new star (unlike the inner radiation of this new star which does warm the surface from the inside).

Conclusion

I can’t help thinking that all my explanations for the imaginary second law have their own problems. And so I welcome explanations from promoters of the theory for the physical processes that take place near the surface of the new star.

Perhaps the problem is in the thought experiment itself. After all, you can’t just fly a star in from another galaxy and park it close to the sun. Barking mad!

And so the Imaginary Second Law of Thermodynamics lives on!

Update – the imaginary law also covered (possibly created by) On the Miseducation of the Uninformed by Gerlich and Tscheuschner (2009)

Update – a worked example with the maths, Radiation Basics and the Imaginary Second Law of Thermodynamics

Update – and more explanation with reference to one advocates explanation of this imaginary law, Intelligent Materials and the Imaginary Second Law of Thermodynamics

Read Full Post »

« Newer Posts - Older Posts »