Graphing Out Loud: curves and lines

I love graphs – my eyes quickly glaze over at a table of numeric data, but a graph, used correctly, can quickly and easily tell the whole story.

‘Used correctly’ is the key phrase – for all their power, graphs are infamously easy to bungle, and when used incorrectly they can misinform – or lie outright.

I’m going to look at an example that touches on a few graphical and statistical concepts near and dear to my heart, as well as carbon geochemistry.

Fig. 1: An image from C3Headlines; the 3 C's are "Climate, Conservative, Consumer". Oh, and the article is titled "The Left/Liberal Bizarro, Anti-Science Hyperbole Continues". It sure would be tragic if they made obvious n00b mistakes after using such language. Click for link!

Coming from an article on the website C3Headlines, this image claims that carbon dioxide concentrations have ‘Linear, Not Exponential Growth’. thereby ‘expos[ing] the lunacy of typical left/liberal/progressive/Democrat anti-science’, The author has reached this conclusion by graphing January CO2 levels* and fitting a linear trendline to them.

Already this is a warning sign – the comparisons the author makes are entirely qualitative, apparently  based up on eyeballing the graph. However, trend lines are created by a statistical process called a linear regression, which comes with a caveat: it will fit a trend line to ANY data given to it, linear or nonlinear. Fortunately, there are also ways of evaluating how good a trend line is.

temperature aNOMalies

If you are new to climate science, you might be wondering what, exactly, this ‘temperature anomaly’ thing is that you keep hearing about. I know I was a bit confused at first! This post explains the concept, using a real-world example.

Absolute temperatures (yearly averaged) from two sites in the UK: one urban (St. James Park, green) and one rural (Rothamsted, red). Although the urban site is consistently warmer, the two sites show the same warming trend. But is there a way to compare them directly? Data from Jones et al. 2008, kindly provided by Dr. Jones.

Cities tend to be warmer than their surrounding countrysides, a fact known as the urban heat island effect (UHI). This occasionally is offered as an alternative explanation for greenhouse warming, but it fails on closer inspection. We can use data from Jones et al. (2008) [PDF] to see one reason UHI can’t explain observed warming. One time series is from St. James Park, in the city of London; the other is from nearby Rothamsted, a rural site some tens of miles away. As you can see, the urban location is consistently about 2 C warmer; however, the warming is nearly identical at both sites (a strongly significant 0.03 deg C/year). Jones et al. note:

“… the evolution of the time series is almost identical. As for trends since 1961 all sites give similar values …  in terms of anomalies from a common base period, all sites would give similar values.”

This gives us a hint about what a temperature anomaly is:

Did chaos theory kill the climatology star?

Last time, we saw that some mathematical systems are so sensitive to initial conditions that even very small uncertainties in their initial state can snowball, causing even very similar states to evolve very differently. The equations describing fluid turbulence are examples of such a system; Lorenz’s discovery of extreme sensitivity to initial conditions ended hopes for long term weather forecasting. Because the state of the weather can only be known so well, the small errors and uncertainties will quickly build up, rendering weather simulations useless for looking more than a few days ahead of time.

But Lorenz’s discovery doesn’t have much impact on climate modelling, contrary to the claims of some climate skuptix. Climate is not weather, and modelling is not forecasting.

Weather refers to the state of the atmosphere at a particular time and place: What temperature is it? Is it raining? How hard is the wind blowing, and in which direction? Climate, on the other hand, is defined in terms of the statistical behavior of these quantities:

“Climate in a narrow sense is usually defined as the average weather, or more rigorously, as the statistical description in terms of the mean and variability of relevant quantities over a period of time ranging from months to thousands or millions of years. [...] Climate change refers to a change in the state of the climate that can be identified (e.g., by using statistical tests) by changes in the mean and/or the variability of its properties, and that persists for an extended period, typically decades or longer. ” IPCC

Many climate skuptik talking points derive from confusing these two quantities, in much the same way that a gambler might win a few hands of poker and decide that they are on a roll.

Although it is generally not possible to predict a specific future state of the weather (there is no telling what temperature it will be in Oregon on December 21 2012), it is still possible to make statistical claims about the climate (it is very likely that Oregon’s December 2012 temperatures will be colder than its July 2012 temperatures). It is very likely that the reverse will be true in New Zealand. It is safe to conclude that precipitation will be more frequent in the Amazon than in the Sahara, even if you can’t tell exactly when and where that rain will fall.

In fact, Lorenz’s groundbreaking paper, ‘Deterministic Nonperiodic Flow’, would seem to endorse this sort of statistical approach to understanding fluid dynamics:

“Because instantaneous turbulent flow patterns are so irregular, attention is often confined to the statistics of turbulence, which, in contrast to the details of turbulence, often behave in a regular well-organized manner.” (Lorenz 1963)

Let’s take a closer look.

Fig. 1. Three solutions of the Lorenz equations, starting at virtually identical points. Although the solutions are similar at first, they rapidly decouple around T=12.

The Lorenz equations consist of three variables describing turbulent fluid flow (X,Y, and Z), and three controlling parameters (r, b, and s). The equations are differential equations, meaning that a variable is described in terms of how it changes over time- saying ‘Johnny is driving west at 60 miles per hour’ is a simple differential equation. In order to solve a DiffEq, you need an initial condition – “Johnny started in Chicago” is an initial condition; without knowing that, you can’t say where she will be after driving for three hours.

John Everett, part II.5: Significance Junkies

A part of my John Everett series – read more: 0/I - II.0 - II.5 - II.75 -  III.0 - III.3 - IV.0 - IV.4 - IV.8 - V - VII - VIII - Full Report 

The CO2 scenarios are literally falling flat and need revision. The observational trend line shows monotonic growth – pretty much a straight line as in the chart below of global marine CO2 measurements (NOAA data)4, while the IPCC scenarios used in most research rely on an accelerating growth. Certainly the predicted rapid acceleration of the IS92a model (see solid black line in middle of the figure on the right) is missing from the NOAA data plotted below. In fact, if the last 8 or 12 years are representative of the future, we might imagine a downward slope in the growth rate.

Last time, we looked at one claim Dr. Everett makes in this paragraph: that the measured rate of change in atmospheric carbon dioxide is inconsistent with the emissions scenarios used to predict future ocean acidification. To do this, he plays fast and loose with quantities and their derivatives (the rates at which they change.) The imprecision extends even to his quoted numbers: in pidginthe previous paragraph, he gives a growth rate as “3.05 ppm”. That’s a not a growth rate; it’s a concentration. He means 3.05 ppm per year. His projection is an extrapolation of “the average rate of increase for the past 10 years (1.87/year)…” 1.87 WHAT per year? I know that he means 1.87 ppm/year, but a lot of people wouldn’t, and I shouldn’t have to make assumptions. If Everett is being sloppy with his units, he’s being sloppy with his science.

The other claim that Dr. Everett draws from the rate of change in CO2 is that “the growth rate seems to be leveling off, if not declining [...] In fact, if the last 8 or 12 years are representative of the future, we might imagine a downward slope in the growth rate. ” Look at the graph of the growth rate again. It goes up and down- a lot.

The record of changes in atmospheric carbon dioxide since 1980. It's got it's ups and downs. Click to see the full record back to 1959.