Saturday, August 31, 2013

Significant warming

I'm seeing more frequently conflation of "no recent warming" with "no statistically significant warming". At WUWT I saw a comment by Willis Eschenbach:
"the UAH record shows no trend since August 1994, a total of 18 years 9 months."
This surprised me, because I knew UAH had risen quite a lot, so I noted that the trend was actually 1.38°C/century over that period. I was reproved
"Not statistically significant is the same as no trend far as we’re concerned."
and indeed Willis came in later to say
"I assumed you knew that everyone was talking about statistically significant trends, so I didn’t mention that part."

Phil Jones was famously hassled about whether there had been statistically significant warming since 1995, which lead to lots of "no warming" stuff. And there were even Keenan-inspired questions of the Met Office in the HOuse of Lords.

You'll notice this talk is coming from skeptics. Scientists don't spend a lot of time worrying about whether trends are "significant". There's good reason for that.

Significance tests can't prove anything. They seek to disprove a null hypothesis. And the question of whether the temperature trend was zero is of little interest. No-one expects that it would have been. So the question of disproving that is not important to them. And I can't imagine why skeptics think it is important to disprove that the trend was zero.

Lucia writes a lot about significance tests, but she always tests trends which, rightly or not, she associates with AGW sources. That makes sense - disproving those could mean something.

Anyway, Willis' example above illustrates what is wrong with equating no warming with no significant warming. It's quite possible for warming exactly in line with AGW predictions to be sytatistically insignificant, because of the number of observations and their noisiness. Now a theory can't do better than get it right. So in this post, I'll show how, under various measures of significance, various trends do not become significant until sustained for quite a lot of years. That's not because there's any doubt about whether they are happening.

I'll also look at the effect of measures. I'm planning a few posts on significance and autocorrelation. This arises partly because I'm now maintaining current data in my triangle trend plots. SkS has a trend calculator, and they give significantly different trend levels. The reason is that they use a different time series model. They give wider confidence intervals which lead to longer quoted periods of "no significant warming". I'll discuss the implications, and possible other methods.


Models of statistical significance

I'll say more about the maths of this in an upcoming post. For the moment, I'll just identify three models of linear regression residuals:
  • Naive OLS
    Assumes that after linear regression, the residuals are independently distributed. Markedly untrue for monthly surface temperature data. Not bad for annual averages.
  • Autocorrelated residuals.
    With an AR(1) model, the one most commonly used (if not naive). Often in conjunction with a Quenouille approximation to get the uncertainty. This is quite good for moderate autocorrelation, and if the model is inadequate, it is probably not because of that approx. My trend viewer uses it.
  • ARMA(1,1) model, as described by Forster and Rahmstorf. This has an extra degree of freedom, which allows the autocorrelation model to flatten when the lag-1 correlation indicates a steep slope. It is used by the SkS calculator.
Each extra parameter tends to increase the uncertainty of the trend, because it implies greater relations between the data points and hence fewer degrees of freedom.

Methods

I'm using the R arima function. The actual call for a time series T is:
 h=arima(T,k,xreg=time(T)/100)
where k is c(0,0,0) for OLS, c(1,0,0) for AR(1) and c(1,0,1) for ARMA(1,1).
h is a structure with a vector of coefficients h$coefs and a covariance matrix h$var.coef. The last coef is the trend, and the last diagonal in the cov matrix is taken to be the variance. In fact there's no simple variance for a multiple parameter model, but the matrix is diagonally dominant.

Significant trends

I'll plot the 95% uncertainty levels of trend (1.96σ), for periods ending Jul 2013 and starting in years going back in time, for each of the five common datasets, HADCRUT4, GISS Land/Ocean, NOAA Land/Ocean. UAH T2LT and MSU RSS T2LT. For any given uncertainty level, trends starting to the right of the line are not significant. I'll also plot the actual trends, in the same color; where the two curves cross is where the period of significantly positive trend starts. I'll do this for each of the three model types.

But first, because it's easier to explain, I'll show just one data set, UAH T2LT, with the three models. The black curve is the trend for the period from the starting date on the x axis to now (July 2013). The colored curves are the 1,96σ levels for each of the models for the residuals. Where the black curve is above the color, the trend from that year is significant. Where they cross is the most recent month for which it could be said that there has been no significant trend. A complication is that it may cross more than once.

Because the se for the ARMA(1,1) model is larger than the others, it gives a longer period of "no statistically significant trend". However, it's worth noting that the crossing point is at quite a high level; nearly 1 °C/century. That's a general pattern, because generally the trend curves have a negative slope in this region. It illustrates the point above - the trend may not be significant, but may also be fairly close to what was predicted.

It's interesting that the CI does not monotonically decrease with increasing period, but has a small spike around 1997. The error of the trend actually increases here (going back) because of the big 1998 spike with dips on each side. We;ll see that this is mainly a feature of the satellite measures.

I've marked a horizontal at 2 °C/century, which is often said to be an AGW predicted level for this time. You can see that the multi-parameter models which give a long period of no significance, would say that even this level, agreeing with prediction, was insignificant for nearly a decade.

Now here are a series of plots showing the trends, CI's and crossing points for each data set - one plot for each stochastic model:





The same pattern - as parameters are added the CI's widen and the "no warming" period extends back, but the actual degree of warming becomes increasingly positive.

Conclusion

I don't think scouting around for a period free of significant trend is a useful activity, because it doesn't actually prove that the theory made a bad prediction. For that you have to test the deviation from the prediction.

But if you really want to do it, you have to deal with the fact that there are different models which give substantially different answers. And models that allow long periods mark as insignificant periods that are increasingly in line with what was predicted.

In my next post, I'll look at the Akaike Information Criterion, which gives at least some objective criterion for choosing models. I'm also hoping to present a model-free approach.

Friday, August 23, 2013

More maintained data vizualization pages

Most of the data vizualization pages that you see on the right now have a scheme for automatic updating. I have added this capability to two of them:

Monthly GHCN plots, with WebGL and shading. These will be updated approx weekly.
Temperature trend viewer, also updated weekly.

These pages are updated daily

Latest Ice and Temperature.
Collection of High Resolution NOAA SST with WebGL.
Regional Hi-Res SST movies..



Friday, August 16, 2013

July GISS Temp down by 0.12°C



GISS LOTI went from 0.66°C in June to 0.54°C in July. This just about cancels the previous month's rise. It's also the second month where GISS moved much more than TempLS. Satellite lower troposphere temperatures fell by comparable amounts.


Here is the GISS map for July 2013:



And here, with the same scale and color scheme, is the earlier TempLS map for July:



Previous Months

June
May
April
March
February
January
December 2012
November
October
September
August
July
June
May
April
March
February
January
December 2011
November
October
September
August 2011

More data and plots


Wednesday, August 14, 2013

Temperature drives CO2?


For years now, we have been seeing graphs of CO2 and temperature records, with an observation that the pattern of CO2 lags temperature change. The claim is then that temperature drives CO2, not the other way around. Often this is expanded to the proposition that the rapid recent increase in air CO2 has nothing to do with the huge amount that we have emitted.

In the past, indeed CO2 has responded to sea temperature changes. That is no surprise. And CO2 has not been driving temperature. The reason is that nothing has been driving CO2. Unless something pushes up (or down) CO2, it can't force anything. And for millenia, the total amount of carbon cycling through various forms on the Earth's surface has been fairly constant.

None of this has anything to do with our present situation. We are looking at the response to burning 370 Gigatons of carbon (so far). There's no doubt what is driving that CO2 production. And it's a new introduction to the surface total. It is about 2/3 of what the atmosphere held 150 years ago, so it's big. We are about to find out how CO2 drives temperature.

Below the break, I'll look at some recent expressions of the Temp drives CO2 meme.

Ice ages

This is the traditional form. You see a plot like this:

with the observation that CO2 peaks and rises lag temperature rises by a few hundred years. That article starts:
"There are still people who insist that changes in CO2 can explain the pattern of glacial and interglacial periods."
but doesn't say who they are. Know anyone? It's not the IPCC! In the AR3, in 2001 with the evidence very recent, they weren't sure, but had no trouble with the idea that CO2 lags temperature:
" From a detailed study of the last three glacial terminations in the Vostok ice core, Fischer et al. (1999) conclude that CO2 increases started 600 ± 400 years after the Antarctic warming. However, considering the large uncertainty in the ages of the CO2 and ice (1,000 years or more if we consider the ice accumulation rate uncertainty), Petit et al. (1999) felt it premature to ascertain the sign of the phase relationship between CO2 and Antarctic temperature at the initiation of the terminations."

Fluctuations in recent CO2

For a while now, variants of graphs matching short-term changes in CO2 with temperature have been circulated, popularized particularly by Murry Salby in various lectures (though not written documents). An apostle is a commenter called Bart, or Bartemis. In a typical version, the annual difference of CO2 is said to be predicted by the difference between average temperature anomaly and an "equilibrium value", which is a fudge parameter. Plots like this are shown:

It shows Southern Hemisphere temperature compared with the annual difference of the CO2 series, with an offset to get them to match. It shows a good correspondence of the El Nino peaks, and some ability to trach the other short term ups and downs.

Why the SH, you might ask? Bart says it's not a cherry pick, but is because the SH has most ocean. And yes, the global temp is not so good:

The peaks still line up, but the trends drift apart. NH is worse. The reason is that the SH had a particular ability to match the fairly modest rate of temp rise with the longer termCO2 difference trend, while keeping the scaling right. Other plots don't. It's hard to see that as an ocean effect.

It illustrates the main issue. Diffencing the CO2 curve emphasises the short term fluctuations in what is otherwise a rather boring CO2 history. But it is the long term AGW rise that is of interest, and differencing suppresses that. The trend would have been there as a constant, except that the two curves are arbitrarily displaced to match, so even that has gone. The wiggle matching in no way explains the main effect - the 30% or so increase in CO2 since we've been emitting.

This is clearly shown in the equivalent AR3 plot of differenced CO2:

Now the differences are not displaced, and are clearly positive. It's that steady positive value that is to be explained. And the AR3 shows the explanation - the emissions, which are much steadier, but rise with the air increases, and always exceed them. The question is not how the CO2 got into the atmosphere, but where some of it went.

Predicting CO2

That wiggle matching has been going on for a while, but Bart also has a plot which does try to do the right thing - to integrate his differential form and predict CO2 rise from temperature (his plots here). He calls it "very high fidelity":


But is it? I've reproduced the last plot, the differences in red, here, along with a quadratic regression in black:

Here the CO2 curve is just 24-month smoothed, and the temperatures were modified as integrated 0.184*(GISSLO-0.424). These are slightly different to Bart's. I used his numbers first, but the curve was rather different; with these parameters it looks the same. There may have been some rounding. In any case, the point is that the quadratic regression fits much better than the GISS derived version.

But the quadratic version could have been expressed as integrated 8.832E-5*(t-429) where t is time in months from the start (Feb 1959). It's of the same form as the GISSLO fit, but using t instead of GISSLO. So simple date is a better predictor than temperature!

Well, not quite. In fairness, I should optimise the coefficients by regressing against GISS as I did against time. That gives the green curve in the above plot. It's very close to the time regression. So a better conclusion is that GISS is no better that time. On the longer term, there is absolutely no basis for saying that temperature drives CO2.

Conclusion

Short term temperature fluctuations, like ENSO, cause short term CO2 fluctuations of a few ppmv. And very large longer term fluctuations, like Ice Age terminations, cause larger fluctuations of up to 100 ppmv. All dwarfed by the current increase, which is clearly driven by our emissions.

Differencing the recent CO2 history shows up short term fluctuations, which show this matching. But this says nothing about the effect of our steady emissions, which differencing and translating removes. When the differenced version is integrated as a predictor, it offers no benefit.




IPCC derivative

Monday, August 12, 2013

TempLS global temp down 0.03°C in July



The TempLS monthly anomaly for July 2013 was 0.469°C, compared with 0.499° in June. This continues a sequences of small changes. Satellite indices decreased more.

Here is the spherical harmonics plot of the temperature distribution:



Warm spots in East US and N Europe. Australia also continued its long warm spell. Arctic rather cool.


And here is the map of stations reporting:





Monday, August 5, 2013

Buffers, pH and ocean acidification

This post is a followup to some recent discussion of ocean acidification at Climate etc, which provoked scorn from Stoat, and in turn snark from WUWT.

Inevitably, the first comments are along the lines of "can't be acidification until pH7" (never with a reference). I blogged on that here. I'd like to take up two points from that - Lewis acidity and buffers. The general aim of this post is to put pH in context. It's over-used.

Lewis Acid

Early definitions of acidity relied on the properties of hydrogen ions, and of course pH is central. But in 1923, Lewis generalised the notion of acidity to the sharing of electron pairs, which needn't involve hydrogen.

That's an important simplification. In fact, H+ has very little final role in the seawater reactions. You can express them as:
CO2 + CO3-- + H2O ⇌ 2HCO3-
Ca++ + CO3-- ⇌ CaCO3
No protons, but the Lewis acid CO2 reacts with the base CO3--, shifting the solubility equilibrium.

Buffers

Buffering is often thought of as pH buffering, but again I'll aim to generalise. In fact, it's just a consequence of ternary equilibrium:
H + A ⇌ HA
The Law of Mass Action gives the equilibrium relation:
[H][A]/[HA] = K
If K is of moderate size, and [H] is small relative to [A] and [HA], then whether the reaction moves back or forward, [H] remains relatively fixed because the dominant species don't change much in proportion.

The familiar pH version arises when HA is a weak acid and [A] a weak base. Then [H+] is small, and relatively fixed (buffered) by the equilibrium.

A common sceptic scoff at acidification is, how could so few H+ ions make a difference? At one level it's right - they can't directly. But what buffering means is that if [H+] changes, then something else has changed, by a lot. And that is what is likely to cause changes.

Sea-water reactions

These were listed above:
CO2 + CO3-- + H2O ⇌ 2HCO3-
Ca++ + CO3-- ⇌ CaCO3(aragonite)
I'll ignore the effects of hydration of both H+ and CO2

There is an excellent review article by Zeebe which sets this out in much more detail, and I'll refer to it later.

The effect is that when CO2 is added, the first equilibrium shifts to the right, CO3-- is diminished, shifting the second equilibrium to the left. The equations can be combined:
CO2 + CaCO3 + H2O ⇌ 2HCO3- + Ca++
This emphasises that ultimately each extra molecule of CO2 tends to dissolve a molecule of CaCO3, if equilibrium shifts to the right.

pH Equilibrium relations

The equilibria can also be expressed as the more traditional diprotic acid-base expressions

CO2 + H2O ⇌ HCO3-+ + H+; pKa=5.94

HCO3- ⇌ CO3--+ + H+: pKa=9.13
The pKa values are the log10 of the equilibrium constants.
The equilibrium can be shown against pH as Bjerrum plot. This shows relative concentrations normalised against DIC = dissolved inorganic carbon. This is just the sum:
DIC=[CO2]+[HCO3-]+[CO3--] Here is the plot:

RTA stands for Relative (to DIC) Total Alkalinity - see next.

Equilibrium relations - some numbers

I'm following the Zeebe review article here. Typical values for dissolved species:
  • DIC 2.05E-3 M (M = 1 mol/l ~ mol/kg)
  • Total Alkalinity=[HCO3-]+2*[CO3--] = 2.35E-3 M
  • ratio %: [CO2]:[HCO3-]"[CO3--]=0.5:89:10.5 at pH 8.2
  • So [CO2]=1E-5 M, [HCO3-]=1.82E-3 M, [CO3--]=2.1E-4 M
  • At pH 8,2, [H+]=6.3E-9 M; [OH-]=1.6E-6 M;

Sea Water measurement

Eli had a post arising from these threads on pH measurement. It's hard. However, there is a new method, probably less accurate, but suitable for continuous monitoring. It uses an ion-selective field FET, SeaFET, which can be set in place for months at a time. But, as insisted here, pH is over-emphasised. It follows the major species, and they can be used to measure it. Two quantities easily measured by flask analysis are all you need. One is DIC. This can be measured gravimetrically. The other is Total Alkalinity (TA) TA=HCO3-]+2*[CO3--] This can be determined by acid titration with an indicator that changes about pH=4.5, which removes almost all bicarb and carb. The titration also picks up borate and hydroxide, and [H+] needs to be subtracted, but these are small effects which I will ignore here. All you need are the values of DIC and TA - both easily obtained by flask analysis. The point is they are present in much larger concentrations, and are more stable.

Use of equilibrium relations

Given DIC for normalisation and any other concentration (in particular TA), you can read a pH value from the Bjerrum plot, and deduce everything else. But here are some more convenient log plots. Firstly with pH again on the x-axis, but log10 on the y-axis. Note that both log([CO3--]) and log([CO2]) vary close to linearly
Now shifting the relative total alkalinity RTA to the x-axis. This allows the other values to be read explicitly:
 
RTA is easily obtained, but doesn't preserve the linear relations. Here is log([CO2]) on the x-axis:

Solubility equilibrium

The end problem is dissolution of CaCO3 structures. Two reasons are sometimes given why this may not matter

  • Biological deposition is complex - not just a solubility product
  • The oceans are supersaturated relative to aragonite formation (the CaCO3 form used by most organisms).
Deposition is complex, though I have trouble believing that reducing [CO3--] won't make it harder. But in any case, the structure once formed has to resist solution, at all times, when conditions are fluctuating. Another reason advanced for unconcern is that pH continuously monitored shows fluctuations large relative to the long term changes wrought by CO2 increase. It's not clear whether these track changes in carbonate, but if they do, then clearly the structures are intermittently exposed to dissolving conditions much worse than average. And dissolving once is enough.

Discussion

What I hope these plots will show is that pH is just one of a number of co-varying quantities, and because of the small number of ions represented, is far from the most important. You can use it to track the reaction if convenient, but you don't have to. It's an intermediate - the nett reaction is CO2 neutralising CO3--, producing HCO3-. This then affects the CaCO3 solubility equilibrium. The secondary role of pH is relevant to the following things that have been said in the acidification discussions:
  • "Acidification requires pH<7
    pH is just an indicator of the carbonates equilibrium. pH 7 has no relevance.
  • The number of H+, even with acidification, is too small to matter
    True, H+ are not the problem
  • pH measurements are sparse and inaccurate
    But not TA/DIC, which is enough.