Hansen Scenarios A and B – Revised

This is a somewhat restated version of an earlier post seeking to understand the differences between Hansen Scenarios A and B. Rather than trying to clarify matters here, Gavin Schmidt posted over at Tim Lambert’s. In this morning’s post, I correctly identified that the difference between Scenarios A and B for periods up to the present pertained primarily to the handling of CFCs and had nothing to do with one being “exponential” and one being “linear”, as Gavin had stated at realclimate. However, I got wrongfooted somewhat in my interpretation of the data as archived at realclimate.

In Hansen et al 1988, they stated that they dealt with other CFCs and trace gases by doubling the effect of CFC11 and CFC12, a point that I noted in my post yesterday. They said:

Potential effects of several other trace gases are approximated by multiplying the CFC11 and CFC12 amounts by 2.

In my post yesterday, I incorrectly surmised that this would not be a substantial effect. When I analysed the data as archived at realclimate, I assumed that the effect of other trace gases would be counted by a separate line item in which they set the impact equal to the combined CFC11 and CFC12 impact. I did not consider the possibility that they would actually insert incorrect CFC11 and CFC12 values into their Schedule A physical inventories, making these values not directly reconcilable to physical measurements. However, this turns out to be what they’ve done. Schmidt snickered at my failing to consider the possibility of them altering their physical inventories to incorrect values. I see now what they’ve done and can follow through their reasoning, but I’d have been inclined to keep the Other Trace Gases as a distinct entry so that people can keep track of it. Climate scientists can be tricky accountants.

The more substantive issue is that this ad hoc handling of Other Trace Gases is what accounts for the near time differences between Scenarios A and B – a point that I made this morning, even if I didn’t quite connect to the doubling for OTGs.

Here is a diagram showing CO2 concentrations for Scenarios A,B and C (per realclimate version), as compared with observations to 2006 (GISS). Obviously, there is no material difference between Scenario A and Scenario B concentrations or, for that matter, with observations up to 2006 or so. CO2 concentrations start to diverge between the two scenarios in the next decade, but, for the near time analysis, calling one graph “exponential” and another “linear” is not really salient to the quantitative analysis (contra Schmidt’s characterization at realclimate.)

cfc_ha45.gif
Figure 1. CO2 concentrations. A,B, C from realclimate; observed from NASA GISS.

Next methane, the second most important GHG. The left panel shows methane levels in Hansen et al 1988 as compared with present estimates; the right panel shows radiative forcing estimates. There is a slight difference between Scenario A and Scenario B up to 1998, when the Hansen-Michaels dispute occurred, but the difference is insufficient to yield the difference in Scenario A and B reported outcomes. A couple of interesting points though: the levels in Hansen et al 1988 reported for 1987 are noticeably higher than levels presently believed to have existed in 1987. Is this because of changing ways of measuring methane content? At present, I don’t know, but the differences are surprisingly high, as compared to CO2. As noted yesterday, methane concentrations have not increased nearly as much as Hansen projected, following a course parallel to Scenario C at a lower concentration. The forcing estimates obviously parallel the concentration estimates.

 cfc_ha46.gif  cfc_ha54.gif

There is negligible difference between Scenario A and Scenario B N2O.

The main action is in CFCs and, in particular, the Trace Gases other than CFC11 and CFC12. Here’s how Hansen et al 1988 described their CFC scenarios:

In Scenario A, … CCl3F (F-11) and CCl2F2 (F-12) emissions are from reported rates (Chemical Manufacturers Association (CMA) 1982] and assume 3% yr-1 increased emission in the future with atmospheric lifetimes for the gases of 75 and 150 years respectively. ,, In Scenario B… the annual growth of CFC11 and CFC12 is reduced from 3% yr-1 today to 2% yr-1 in 1990, 1% yr-1 in 2000 and 0 in 2010…. In Scenario C, CFC11 and CFC12 abundances are the same as scenarios A and B until 1990; thereafter CFC11 and CFC12 emissions decrease linearly to zero in 2000.

As noted above, they also said:

Potential effects of several other trace gases are approximated by multiplying the CFC11 and CFC12 amounts by 2.

The Figure below shows CFC11 and CFC12 concentrations, using the data from realclimate only pertaining to CFC11 and CFC12. In this calculation, the arbitrary doubling of CFC11 and CFC12 values in Scenario A (representing OTGs) has been backed out and will be shown separately under OTGs. Scenario A and B projections for CFC11 and CFC12 are identical ; and thus this is not a reason for Scenario A and B differences. Both scenarios over-estimated actual concentrations, which were like Scenario C.

 otgsh69.gif  otgsh70.gif

In my first try at this post yesterday, when I noted the doubling as an ad hoc method for Other Trace Gases, I hypothesized that the effect would not be material, but this does not prove to be the case. The figure below show how the primary near-time difference between Scenarios A and B arises. The solid black line shows the Scenario A radiative forcing (equal by assumption to the sum of Scenario A-B forcing from CFC11 and CFC12) – Hansen’s “Business as Usual” case; the dotted line (zero for all years) shows the Scenario B forcing (Hansen’s “more plausible” case). The red points show the radiative forcing for OTGs from realclimate; green lines are a first try at IPCC A1B and A2.

otgsh68.gif

So where does this leave us? To clarify one point, I don’t have any problem with properly articulated scenarios as a way of analyzing projections. How else are you going to do it? With respect to projections of CO2 concentration made in the 1980s, the projections are pretty much on the money. But even if they were off to some degree, all it would do is defer or accelerate the date at which CO2 doubled somewhat. And I agree with concerned climate scientists that doubling CO2 is an unplanned real-life experiment that could have very serious consequences. If the impact is 2.5 deg C, then we should probably be glad that that’s all it is, as the results could have been much worse.

If Hansen’s 1988 model didn’t work very well, then that does not mean for me that some other model mightn’t work, even some other model by Hansen. Many readers here take much more hard-edged positions on these things than I believe to be justified.

Having said that, I’m also trying to figure out how one gets from point A to point B in all of this, a path made more difficult because the field is all too often characterized by sloppy practices.

The different handling of Other Trace Gases was material to Scenarios A and B outcomes, accounting for virtually the entire difference between scenarios. Scenario B hypothesized that there would be no contribution from the Other Trace Gases. Why did Hansen think in 1987 that this was “more plausible” than the doubling in Scenario A. With the benefit of hindsight, it looks like OTG forcing was about halfway between Scenarios B and Scenario A (before any counting from tropospheric ozone etc. which would move the results closer to Scenario A.

What then is the basis, if any, for concluding that Scenario B was the more “plausible”? I doubt whether Hansen was really thinking about differences between near-time Other Trace Gas concentrations when he said that Scenario B was “more plausible” for the near-time dispute with Michaels. The 1988 paper strongly suggests that he was thinking about “deep time” differences out past 2025 where exponential versus linear actually matters. However the term was “convenient” when Hansen re-visited the matter in his 1998 debate.

However, trying to put oneself in Hansen’s shoes back in 1988, there’s nothing to suggest that Hansen thought that it was “less realistic” to account for OTGs by doubling than to ignore them (Scenario B.) Of course, this ad hoc method appears to have over-estimated the impact of OTGs, but, once again, there’s nothing to suggest that Hansen in 1988 thought that it was less “realistic” to ignore the contribution of OTGs.

As to how Hansen’s model is faring, I need to do some more analysis. But it looks to me like forcings are coming in below even Scenario B projections. So I agree that it’s unfair for Hansen critics to compare Scenario A temperature results to actual outcomes as a test of the model mechanics. On the other hand, Hansen’s supporters have also been far too quick to claim vindication given the hodgepodge of GHG concentration results. If it’s unfair to blame the blame the model for differences between actual and projected if the GHG projections are wrong, then it i equally unfair to credit the model with “success” if it gets a “right” answer using wrong GHG projections. One would really have to re-run the 1988 model with observed GHG concentrations to make an assessment. Given that GISS have changed their models since 1988, GISS would presumably argue that the run is pointless, but the cost of doing the run doesn’t appear to be large and it seems like a reasonable exercise for someone to do. It would be interesting to obtain a listing of the 1988 model to that end.

Pielke Comparison of Trends 1990-2007

Roger Pielke Jr reports that he has posted up a comparison of Hansen, IPCC, and observed trends 1990-2007 here:

http://sciencepolicy.colorado.edu/prometheus/archives/climate_change/001330temperature_trends_1.html

and suggested that I start a thread (done!).

Publishing NASA Data at Realclimate

In 2006, Willis Eschenbach digitized Hansen A,B and C scenarios to 2005 (see here ). I used these values together with my own visual digitization extension to 2010 in my recent post. Subsequently Lucia drew my attention to a digital version of the NASA data placed online at realclimate here , referenced in a realclimate post here.

I carried out a routine comparison of the two versions. The realclimate version of Hansen Scenario C was 0.166 deg C warmer than Willis’ digitized version. In reviewing the data, realclimate Scenario C was higher than realclimate Scenario B, so an error has been introduced somewhere in the process.

Which raises the question: how did this error get introduced into the NASA data published digitally for the first time at realclimate? Did it get introduced in digital copying? Or did NASA itself digitize the Hansen scenarios from print media and introduce the error then? In fact, exactly what is the provenance of the digital version presently archived at realclimate? Gavin did not say in his post. Did Gavin digitize the print media? Did someone else digitize it? Or is it digital data? [Update: these matters have been resolved]

Steve: For some reason, the 1994 value of Scenario C is lower than the 1994 value of Scenario B in the original graphic, notwithstanding the lower forcing of C versus B. It appears that Willis’ digitization is from a later and muddier version of the graphic and is incorrect at a few points though the average difference is not material. This leaves the question of why C is higher than B in 1994.

Update Jan 21: Willis digitized an image from 1999. Gavin said in an RC comment on Dec 22, 2007 (but not in the note itself) that the scenario data was digitized (presumably from the Hansen et al 1988 graphic). However,as far as I can presently determine, the other data sets e.g. radiative forcing in wm-2, also archived by Gavin, does not correspond to a Hansen et al 1988 graphic and therefore could not be digitized from a Hansen et al 1988 figure. Maybe it’s a digitization from another publication or maybe it was calculated; hard to say at present. I’ll email Schmidt and ask him.

Hansen GHG Concentration Projections

In the various disputes over Hansen et al 1988, Roger Pielke Jr and NASA apologist Eli Rabett (who has been said to be occasional NASA contractor Josh Halpern) have each attempted to disentangle the forcing projections implied by Hansen et al 1988 – Pielke here and Rabett here for CO2 and here for other gases.

Commenters on the previous thread have naturally puzzled over the actual differences in GHG projections between Scenario A, B and C and, in an effort to clarify matters, I’ve attempted to calculate GHG concentrations through to 2030 implementing the assumptions described in Hansen et al 1988. (Roger Jr posted up a chart similar to the ones below for CO2). In another post, I’ll apply these results to calculate radiative forcings. I’ve saved my collations of Hansen projected GHG concentrations for the three scenarios at:
http://data.climateaudit.org/data/hansen/hansenscenario_A.dat
http://data.climateaudit.org/data/hansen/hansenscenario_B.dat
http://data.climateaudit.org/data/hansen/hansenscenario_C.dat

The collation of these tables is shown in the script http://data.climateaudit.org/scripts/hansen/collation.hansen_ghg.txt. I’ll review the results some more when I get to radiative forcings, but the differences between major GHG concentrations in Scenario A and B are very slight and it’s a little puzzling how the differences arise between the two scenarios. I’ll look at each GHG contribution below. Continue reading

Hansen Scenarios A and B – Original

I’m making some changes to this post in light of comments from Gavin Schmift. Continue reading

Thoughts on Hansen et al 1988

Update (Jul 28, 2008): On Jan 18, 2008, two days after this article was posted, RSS issued a revised version of their data set. The graphics below are based on RSS versions as of Jan 16, 2008, the date of this article, and, contrary to some allegations on the internet, I did not “erroneously” use an obsolete data set. I used a then current data set, which was later adjusted, slightly reducing the downtick in observations. On Jan 23, 2008, I updated the graphic comparing Hansen projections using the revised RSS version. Today I re-visited this data, posting a further update of this data including the most recent months. While some commentators have criticized this post because the RSS adjustment reduced the downtick slightly, the downtick based on the most recent data as of July 28, 2008 is larger than the RSS adjustment as of Jan 2008.)

In 1988, Hansen made a famous presentation to Congress, including predictions from then current Hansen et al (JGR 1988) online here . This presentation has provoked a small industry of commentary. Lucia has recently re-visited the topic in an interesting post ; Willis discussed it in 2006 on CA here .

Such discussions have a long pedigree. In 1998, it came up in a debate between Hansen and Pat Michaels (here); Hansen purported to rebut Crichton here, NASA employee Gavin Schmidt on his “private time” supported his NASA supervisor, Jim Hansen here , NASA apologist Eli Rabett believed to be NASA contractor Josh Halpern here . Doubtless others.

It seems like every step of the calculation is disputed – which scenario was the “main” scenario? whether Hansen’s projections were a success or a failure? even how to set reference periods for the results. I thought it would be worthwhile collating some of the data, doing chores like actually constructing collated versions of Hansen’s A, B and C forcings so that others can check things – all the little things that are the typical gauntlet in climate science.

Here is my best interpretation of how Hansen’s 1988 projections compare to recent temperature histories.
hansen20.gif

I’ll compare this graphic with some other versions. On another occasion, I’ll discuss the forcings in Hansen et al 1988. First, I’m going to review the prior history of this and related images.
Continue reading

Unthreaded #30

Continuation of Unthreaded #29

Can I recommend however, that instead of using this thread, commenters should register and use the new CA Forum under the relevant headings. The same rules of conduct and topic apply on the forum of course, but it should be a lot quicker and easier to follow individual discussions than on a single thread of a multi-purpose blog post.

User registration is now e-mail activation rather than Steve or I having to activate new users every few hours, as is password recovery.

I have also installed LaTeX on the board, so that the brilliance of your mathematical reasoning is not obscured by poor formatting. 😉

I think it’s likely that Unthreaded might be retired in the future in favour of the more elegant message board solution, but we shall see.

John A

PS I’m still retired. I just temporarily had some time on my hands. Really.

Pielke Jr on IPCC Predictions

Excellent post by Pielke Jr here

From Lacis et al 1981 to Archer Modtran

In several recent posts, we’ve been reviewing the provenance of the radiative forcing estimate of about 4 wm-2 and the logarithmic form for estimating climate response to increased CO2 levels. This has led on the one hand to several primary references, including Lacis et al 1981 and, on the other hand, to references to realclimate collaborator David Archer’s MODTRAN calculator.

I want to focus today on the structure of these calculations – or, at least, to do so to the extent possible given that the methodologies are poorly described. In our discussions, many readers have volunteered observations relating to formulas for simple absorbance of radiation by gas. However, what these observations fail to address is that the relevant calculation here requires a detailed specification of the atmosphere, which includes not just CO2 levels, but a temperature varying from surface to atmosphere exit, other radiative gases including water vapor, clouds, etc.

Given this atmospheric profile, the calculations yield downwelling (or upwelling) infrared radiation in wm-2. The Archer calculator provides Java results for a variety of conditions, permitting the calculation of upwelling and downwelling radiation at specified altitudes. (It would be nice to have this algorithm available in a Matlab or R function and probably it wouldn’t be too hard to convert.)

What catches my eye about these calculations is that atmospheric conditions are held rigid – as though it were a coincidence that the CO2 radiation-to-space maximum and tropopause are at the same elevation. My own intuition is that the atmospheric profile is itself affected by presence of CO2. In Luboš short note on the matter, he took a similar perspective, hypothesizing that additional CO2 would raise the tropopause – a change in atmospheric profile that doesn’t occur in the rigid calculations of Myhre et al and Archer, where the entire impact is concentrated on downwelling wm-2. The form of the problem strikes me as an interesting type of calculus of variations problem and it would be interesting to see what this type of mathematician could do with it.

On to a review of the sources.
Continue reading

Gore Scientific "Adviser" says that he has no "responsibility" for AIT errors

In earlier posts, we observed that Al Gore’s An Inconvenient Truth claimed that “Dr Thompson’s thermometer” confirmed Michael Mann’s hockey stick, but, when analysed, what Gore described as “Dr Thompson’s thermometer” merely proved to be Michael Mann’s hockey stick mis-identified. No wonder it resembled Mann’s hockey stick – or, to use the phrase more common in climate science, no wonder there was a “remarkable” resemblance.

Recently we’ve seen Pierrehumbert’s hysterics over at RC about Courtillot’s misidentification of a solar series. I’m sure that all of you recall similar hysterics from Pierrehumbert calling Al Gore out on this error. You don’t?

Hu McCulloch of Ohio State University now writes about a recent encounter with Lonnie Thompson, the serial ice core non-archiver:

On January 11, Lonnie Thompson gave a talk on Climate change at Ohio State. After his talk, I asked him if the graph identified by Al Gore as “Dr. Thompson’s Thermometer” in his book and film was really based on his ice core research.

Thompson admitted that an error had been made, and even had a slide ready that showed the data of the Mann Hockey Stick plus Jones instrumental data that Gore’s figure was based on, alongside an average of dO18 z-scores from 6 of his Andean and Himalayan ice cores, similar to the 7-series graph that appeared in his 2006 PNAS article. He stated that he recognized the error right away, and even sent Gore (and Mann, as I recall) an e-mail pointing out the mistake.

When I pressed him if it wouldn’t be appropriate to make a more public announcement, given the high-profile nature of the error, Ellen Mosley-Thompson, his wife and co-author, stood up and offered that it was Gore’s error, not theirs, so that they had no responsibility for it, and that in any event there was no forum in which to make a correction.

I suggested that since OSU’s Byrd Polar Research Center has a website with a News page, it would be trivial and virtually costless to post a press release clarifying the matter for the millions of readers and viewers of Gore’s book and film who are not on Thompson’s e-mail list.

“Gore’s error”. “No responsibility.” “No forum”.

Here’s a picture from the OSU website. See if you can find Waldo. Ohio State University press releases have also stated that Thompson was an adviser to Gore for the documentary. Thompson’s online CV says that he was on the “Science Advisory Board” for Inconvenient Truth prior to its release in April 2006. So he was on the Board but he didn’t bear any responsibility. Sure, Lonnie. Sure, Ellen.

Gore used the term “my friend Lonnie Thompson” and Thompson doesn’t know how to correct the error. Sure, Lonnie.

“No responsibility.” “No forum”. No shame.