Annoying Fortran Formats

Earth to climate scientists: no one uses Fortran punchcards anymore. When you put data onto the web, it doesn’t have to fit onto 80 columns as though it were a punchcard. Also when home computers have 100 GB of memory, you don’t have to squish multiple records onto one line. Hurricane track data as archived – Atlantic here is archived on the web so that it fits into 80-column punch cards just as in the 1960s. (Dendro data is the same.) The first 3 storms are shown below. There are header “cards”, trailer “cards” and fixed format fields. A Fortran program is provided at the site to read the data. To read this into R requires a fiddly little program – it’s easy enough to do, but it’s a pointless nuisance. I’m in the process of making R-tables out of the data for the 6 ocean basins and will save the results as flat files. Aside from the nuisance, the continued use of Fortran makes everything much more work than it really is with modern languages like R.

00005 06/25/1851 M= 4 1 SNBR= 1 NOT NAMED XING=1 SSS=1
00010 06/25*280 948 80 0*280 954 80 0*280 960 80 0*281 965 80 0*
00015 06/26*282 970 70 0*283 976 60 0*284 983 60 0*286 989 50 0*
00020 06/27*290 994 50 0*295 998 40 0*3001000 40 0*3051001 40 0*
00025 06/28*3101002 40 0* 0 0 0 0* 0 0 0 0* 0 0 0 0*
00030 HRBTX1
00035 07/05/1851 M= 1 2 SNBR= 2 NOT NAMED XING=0
00040 07/05* 0 0 0 0* 0 0 0 0*222 976 80 0* 0 0 0 0*
00045 HR
00050 07/10/1851 M= 1 3 SNBR= 3 NOT NAMED XING=0
00055 07/10* 0 0 0 0* 0 0 0 0*120 600 50 0* 0 0 0 0*
00060 TS
00065 08/16/1851 M=12 4 SNBR= 4 NOT NAMED XING=1 SSS=3
00070 08/16*134 480 40 0*137 495 40 0*140 510 50 0*144 528 50 0*
00075 08/17*149 546 60 0*154 565 60 0*159 585 70 0*161 604 70 0*
00080 08/18*166 625 80 0*169 641 80 0*172 660 90 0*176 676 90 0*
00085 08/19*180 693 90 0*184 711 70 0*189 726 60 0*194 743 60 0*
00090 08/20*199 759 70 0*205 776 70 0*212 790 70 0*219 804 70 0*
00095 08/21*226 814 60 0*232 825 60 0*239 836 70 0*244 843 70 0*
00100 08/22*250 849 80 0*256 855 80 0*262 860 90 0*268 863 90 0*
00105 08/23*274 865 100 0*280 866 100 0*285 866 100 0*296 861 100 0*
00110 08/24*307 851 90 0*316 841 70 0*325 830 60 0*334 814 50 0*
00115 08/25*340 800 40 0*348 786 40 0*358 770 40 0*368 751 40 0*
00120 08/26*378 736 40 0*389 718 40 0*400 700 40 0*413 668 40 0*
00125 08/27*428 633 40 0*445 602 40 0*464 572 40 0*485 542 40 0*
00130 HRAFL3 GA1

PDI: Elsner

Willis writes:
Well, looking at these studies is giving me a headache. My latest one is High frequency variability in hurricane power dissipation and its relationship to global temperature, James B. Elsner et al. Continue reading

Curry on the Wegman Reports

Here are Judith Curry’s Comments on the Wegman Report. I appreciate these sorts of contributions and am obviously relying on such contributions (Willis, bender, etc.) more and more. Continue reading

Curry's Comments on Klotzbach

I’m quite happy to publish other people’s criticisms of topical and semi-topical papers. Judith Curry has sent in the following comments on Klotzbach et al, (mentioned on roadmap), abridged by her from a previous post on the tropical listserv, responding to specific talking points on the Klotzbach paper that were posted on Gray’s website.
Continue reading

Re-post: To Browsing Undergraduates

Climate Audit has been considered – at least as a phenomenon – in a couple of courses. Kenneth Blumenfeld’s students had a different reaction than the Georgia Tech students. Earlier this year, I wrote a short comment about a post that Kenneth had made at realclimate about this, which I am re-posting in its entirety. It’s amazing how perceptions differ. The students had to represent various viewpoints and one even chose to represent CA. Kenneth recently asked me to write another into piece of a similar nature for the fall term which I’ve appended below – following the reprint. Continue reading

Hurricane Data Compilation Thread

UKweatherworld has a thread where people don’t comment but simply contribute references on the Holocene Optimum. I’m going to start doing things like this as a way of keeping track of things. People have contributed data but it’s impossible to keep track of. This thread is going to consist only of scripts and data relevant to replicating the hurricane papers. Any general commentary will be deleted. I’ll put some of the data sets into *.txt files so that they can be called. Please provide references to some of the data sets.

Prior discussion http://www.climateaudit.org/p?790, http://www.climateaudit.org/p?803, http://www.climateaudit.org/p?822.

Emanuel, Nature 2005 url Website
Landsea, Nature 2005 url

Webster publications
Webster, P.J., G. J. Holland, J. A. Curry, and H.-R. Chang, 2005: Changes in Tropical Cyclone Number, Duration, and Intensity in a Warming Environment, Science, 309 (5742), 1844-1846. url
Webster P.J., J.A. Curry , J. Liu, and G. J. Holland, 2006: Response to comment on "Changes in tropical cyclone number, duration, and intensity in a warming environment", Science, 311 (5768), 1713c url
Mann and Emanuel, EOS, 2006. url
Curry, J.A., P.J. Webster and G.J. Holland, 2006: Mixing Politics and Science in Testing the Hypothesis That Greenhouse Warming Is Causing a Global Increase in Hurricane Intensity. Bull. Amer. Met. Soc., 87 (8), 1025-1037 url

Moberg’s G. Bulloides

Last year, when Moberg was published, I pointed out witha slightly arched eyebrow that one of the two most important contributors to any 20th century HS-ness in Moberg was the increasing percentage of subpolar foraminifera (G. Bulloides) in the Arabian Sea – intuitively not a direct indicator of warming. Having visited the foraminifera literature in more detail in my reading of Hansen (who is more interested in G. ruber), I noticed some pertinent discussion of G. Bulloides in the Arabian Sea, which confirmed my initial reading. G bulloides percentage directly indicates upwelling cold water, which in turn is linked to a stronger winds (i.e. stronger monsoon), which is said to be associated with warming onshore. Continue reading

Underground Problems with Mann-Holes

Willis writes in:

While researching ocean drill cores at the WCDC, I stumbled across Mann’s borehole data. One of the proxies used for historical temperature reconstruction is "borehole temperature", the temperature down in the ground. In 2002, Michael Mann et al. published a study called Optimal surface temperature reconstructions using terrestrial borehole data. It is available here, for $9.00. In it they use all of the available weapons to construct the temperature proxy “¢’‚¬? EOFs, PCA, and of course that perennial favorite, "optimal fingerprinting", viz:

We employ a spatial signal detection approach that bears a loose relationship with “‹Å“”‹Å“optimal detection” approaches used in anthropogenic climate signal fingerprinting [Mitchell et al., 2001]. In such “‹Å“”‹Å“optimal detection” approaches, one seeks to identify, through generalized linear regression, the estimate of a target signal (as predicted by a model) in empirical data. Detection is accomplished through rotation of the empirical data, in EOF state-space, away from the direction of maximal noise (as estimated from, e.g., a control model simulation). Continue reading

The Georgia Tech Report Card

Thanks to Judith Curry for sending along these candid comments from a couple of her students about climateaudit. There has been discussion at the other thread wshich I’d prefer move to this thread.

Here is the report from the Georgia Tech hurricane class discussion on the climateaudit hurricane threads. Two students were assigned to make presentations: Student #1 is a 2nd year graduate student, slightly older and with a mature and broad perspective; student #2 is a recent Ph.D. awardee with good knowledge of statistics.

Student #1 gave an overview of the blogosphere and climate-related blogging activities, and some history of the climateaudit site. He described climateaudit’s practice as:
1. attacking a paper on global warming, before reading it very carefully or understanding the context of the paper, assuming that the author is either dumb or has an “agenda”
2. a plethora of statistical activity of a fairly rudimentary nature
3. realization that the issues are complex
4. some attempts at trying to gain physical understanding of what is going on
5. realization that the issues are even more complex
6. give up and move onto something else

Student #1 then asked the following questions (which I answered):
1. How influential is climateaudit?
2. What items have they raised that we should pay attention to?
3. What can we learn and avoid the next time?
4. Was Dr. Curry’s blogging time well spent, or did it legitimize and prolong a discussion that in the end hasn’t really accomplished anything?

Student #2 focused on the statistical issues surrounding WHCC and Emanuel papers. He raised the following main points:
1. The climateauditors do not seem to understand parameteric vs nonparametric tests. The Kendall test (rank based test) used by WHCC does not require a normal distribution and is also fairly insensitive to serial correlation, so the emphasis on autocorrelation and distributions did not add anything.
2. The climateauditors show a general lack of physical interpretation and a lack of appreciation of the fundamentally Bayesian approach (if not explicitly, then implicitly) to climate science statistics, whereby physics and prior knowledge suggests your predictors.
3. ARMA (Spanish for weapon) is a brute force method used (not very productively) when nothing is known about the physics.
4. WHCC statistics were robust and appropriate; the Curry et al. BAMS article was unfairly criticized since the readers did not go back to the original paper cited in Figure 1, which explained what went into Fig 1 and how the trend was determined.
5. There were problems with Emanuel’s statistical analysis that should have been caught in the review process
6. Student #2 was pretty hot under the collar about the whole thing
7. “A lot of personal attacks. Not using bad manners… but still personal attacks. An example? Their opening lines on the hurricane thread: There are statistical issues in fitting trend lines to spiky data like this, which bender is well aware of and pointed out in the predecessor thread. If Curry is unaware of these issues, what does that say? If she is aware of these issues and ignored them, what does that say?”
8. “A biased blog that pretends it is not. In terms of most of the statistics they seem to know what they are talking about, but they should. Most of the stuff is part of basic statistical training. While they appear to be curious about some physics, there is a general lack of good physical interpretation.”

Topics raised in the discussion:

People reading only the thread leader and first few posts get the impression that the paper is wrong, when further down the thread the paper gets vindicated. This gives the casual visitor to the site a negatively biased impression of climate science.

One student raised the issue that statistical mistakes such as made by Emanuel (2005) should have been weeded out in the review process; suggested that a “statistical editor” was needed for climate journals to review the papers for basic sound statistical practices.

The students thought that the fact that the climateauditors did not have “external funding” to do this work diminished their credibility

The students agreed that statistics should be done correctly, data should be made publicly available (but extra work should not be done to make the data and programs convenient for the skeptics), and funding sources should be disclosed.

The “biases” of the climateauditors were discussed. Bender was perceived as a hardcore anti- warmer. SteveM and Willis were perceived has hardcore statistical skeptics, assuming that all analyses done be climate people are suspect. Steve Bloom was viewed as a somewhat heroic glutton for punishment. David Smith was viewed as the voice of reason.

I then went on to describe what I thought was useful and interesting about the site and about the hurricane threads, and the blogospheric approach to science. Everyone agreed that the climateauditors spotted things in the Emanuel paper that none of us had spotted.

Overall, the students were pretty negative about the site. I suggested that the two students post their comments; they did not want to, and I agreed to summarize the discussion (I was asked not to mention their names). They viewed blogging on climateaudit as entering a black hole of trying defend yourself against a prejudged guilty verdict. Well, I am not exactly sure what I expected from this discussion, but it doesn’t sound like the younger generation of scientists are very keen to enter the blogospheric discussions on climate science.

Student #2 ended with 3 quotes and a joke:

Bayesian statistics is difficult in the sense that thinking is difficult. Donald A. Berry

Some people use statistics as a drunken man uses lamp-posts”¢’‚¬?for support rather than illumination. Andrew Lang

Facts do not “speak for themselves.” They speak for or against competing theories. Facts divorced from theories or visions are mere isolated curiosities. Thomas Sowell

Two statisticians were traveling in an airplane from LA to New York. About an hour into the flight, the pilot announced that they had lost an engine, but don’t worry, there are three left. However, instead of 5 hours it would take 7 hours to get to New York. A little later, he announced that a second engine failed, and they still had two left, but it would take 10 hours to get to New York. Somewhat later, the pilot again came on the intercom and announced that a third engine had died. Never fear, he announced, because the plane could fly on a single engine. However, it would now take 18 hours to get to New York. At this point, one statistician turned to the other and said, "Gee, I hope we don’t lose that last engine, or we’ll be up here forever!"

Truth Machines

Here is a sociological study of how GCM modelers relate to their models. [Lahsen, 2005. Seductive Simulations?
Uncertainty Distribution Around Climate Models, Social Studies of Science, 35, 895-922.] I kid you not. Lahsen:

“At least at the time of my fieldwork, close users and potential close users at NCAR (mostly synoptically trained meteorologists who would like to have a chance to validate the models) complained that modelers had a “fortress mentality”. In the words of one such user I interviewed, the model developers had “built themselves into a shell into which external ideas do not enter’. His criticism suggests that users who were more removed from the sites of GCM development sometimes have knowledge of model limitations that modelers themselves are unwilling, and perhaps unable, to countenance.”

Sound anything like Wegman’s size-up of the Team? Another comment:

Recognition of this tendency may be reflected in modelers’ jokes among themselves. For example, one group joked about a “dream button” allowing them “Star Wars style” to blow up a satellite when its data did not support their model output. They then jokingly discussed a second best option of inserting their model’s output straight into the satellite data output.

Afterwards re-read the executive summary of the U.S. CCSP in which they discuss the inconsistency between models predicting that the trend in troposphere temperatures would exceed surface trends and observations, which show the opposite. Faced with this inconsistency, Wigley and other senior climate scientists blamed their lying eyes:

A potentially serious inconsistency [between model results and observations] has been identified in the tropics. The favored explanation for this is residual error in the observations, but the issue is still open.

It is quite a remarkable article. Enjoy.