Thursday, July 30, 2020

26. The temperature trend in Australia since 1853 - PARABOLIC

In Posts 18-24 I calculated the temperature trend for each Australian state using all the available instrumental temperature records with more than 40 years of data, and then compared the results with those claimed by Berkeley Earth. These results were summarized in my last post.

The overall picture, based on the actual raw data, was that most states had experienced significant warming since 1980 but that the current temperatures were generally the same as those seen in 1880 or earlier. In most cases, the temperature records exhibited large swings in temperature between 1890 and 2000, typically a large decline in temperature up to 1950 and then a recovery; only Queensland exhibited a temperature trend that was anything like that advanced by the IPCC in their global instrumental temperature record, but Queensland has no data before 1887. And only in the case of Queensland did the temperature construction using the raw data qualitatively agree with the equivalent Berkeley Earth version. In the case of all other states the two sets of data (i.e. the one based on raw anomalies and the Berkeley Earth adjusted anomaly version) diverged significantly the further back in time before 1950 that you look, and generally the divergence was even greater in the 19th century.

What also became abundantly clear was that due to the large fluctuations seen over time intervals of more than 50 years, only temperature records that are at least 150 years in length allow the observer to put those fluctuations into their true context and to draw any meaningful conclusions. The temperature record is not as most climate scientists and the IPCC appear to claim: namely stable and constant before 1900 and exponentially rising thereafter. What has become clear to me is that the data is largely random but on multiple timescales. There are short-term fluctuations in the monthly means of up to ±5 °C, and long-term fluctuations on timescales of many decades or centuries that are often in excess of ±0.5 °C. In previous posts (see Post 9 and Post 17 in particular) I have speculated on the nature of these fluctuations, their scaling properties and whether they are fractal in nature.

In this post I will complete the analysis for Australia by using the regional temperature trends for each state to construct an average for the whole of Australia. The method will be straightforward: to weight each state's overall trend by the area of that state relative to that of the total area of Australia. Based on this methodology the weightings for each state are as follows:

NSW
Victoria
Tasmania
South Australia
Western Australia
Northern Territory
Queensland
     
0.1004
0.0284
0.0085
0.1229
0.3306
0.1776
0.2316

In addition, the trend from each state is renormalized with a fixed temperature offset in order to ensure that the mean anomaly of each state for the period 1961-1990 is zero. After the scaling by weighting and the renormalization, the contributions from each state are added. The resulting trend in the monthly mean temperature for the whole of Australia is shown below in Fig. 26.1 together with the 5-year moving average. As a result of the offsetting process, the temperature changes in Fig. 26.1 are all defined relative to the 1961-1990 average.



Fig. 26.1: The temperature trend for Australia since 1853. The best fit is applied to the interval 1871-2010 and has a gradient of 0.24 ± 0.04 °C per century. The temperature changes are relative to the 1961-1990 average.


The data in Fig. 26.1 shows that the five year average of the regional temperature for Australia, relative to the average for 1961-1990, decreased from a peak of +0.37 °C in 1879 to a minimum of -0.46 °C in 1909. It then recovered from a low of -0.38 °C in 1948 to reach another peak of +0.47 °C in 2007. So while it is true that temperatures rose by 0.85 °C from 1948 to 2007, it is also true that there was an almost equal but opposite fall of -0.75 °C during the 70 years prior to 1948. So, when put into context, the temperature rise since 1948 does not appear that cataclysmic. Yet when you look at the Berkeley Earth adjusted data for the same period, the picture is very different.




Fig. 26.2: The temperature trend for Australia since 1853 based on Berkeley Earth adjusted temperature data. The best fit is applied to the interval 1951-2008 and has a gradient of 1.51 ± 0.06 °C per century. The temperature changes are relative to the 1961-1990 average.


Fig. 26.2 shows the result of the same averaging process for the temperature anomalies in Australia as previously outlined for the data in Fig. 26.1. The only difference is that Fig. 26.2 uses Berkeley Earth adjusted anomaly data, and Fig. 26.1 uses the original data with the correct anomalies calculated for each station. Yet the two datasets look completely different. Whereas the raw data in Fig. 26.1 combines to form a U-shaped curve or parabola, the adjusted data forms a highly asymmetric curve, with little or no change in temperature before 1950 and a steep rise afterwards. That rise has a gradient of 1.51 ± 0.06 °C per century, which equates to a total rise from 1951 to 2008 of 0.88 °C. This is similar to that seen in Fig. 26.1 for the same time period, but without the preceding decline in temperature. However, while the adjusted data in Fig. 26.2 differs significantly from the raw data in Fig. 26.1, it agrees quite well with the trend published by Berkeley Earth and shown below in Fig. 26.3.




Fig. 26.3: The temperature trend for Australia since 1840 according to the Berkeley Earth website.


I believe that the level of agreement between the 12-month and 10-year moving average curves in Fig. 26.3 with those in Fig. 26.2 validates my approach, and therefore also validates the data in Fig. 26.1. The other point to note about the data in Fig. 26.1 is the gradient of the best fit line. This is 0.24 ± 0.04 °C per century and therefore equates to a total temperature rise of about 0.34 °C from 1871 to 2010. This is far more than the actual change of about 0.10 °C and illustrates the caution that needs to be applied to linear trends that are applied to non-linear data. 




Fig. 26.4: A comparison of the actual ten year average temperature trend for Australia since 1853 and the Berkeley Earth adjusted temperature version. The best fit is applied to the actual raw data over the interval 1859-2008 and has a gradient of 0.184 ± 0.056 °C per century. The temperature changes are defined relative to the 1991-2000 average.


If we contrast the raw data in Fig. 26.1 with the adjusted data in Fig. 26.2 by comparing the 10-year moving average of each (see Fig. 26.4 above) we see that after 1980 there is relatively close agreement between the two. However, as we go back in time before 1980, we see the curves diverge. If we compare the temperature change from 1880 to 2008 in each case, we see that the raw data shows evidence of a rise of only about 0.2 °C, while for the adjusted data the rise is over 0.6 °C. This difference is principally due to adjustments made by Berkeley Earth to the data. These adjustments are shown below in Fig. 26.5.




Fig. 26.5: The contribution of Berkeley Earth adjustments to the anomaly data after smoothing with a 12-month moving average. The linear best fit to the data is for the period 1901-2010 (red line) and the gradient is +0.295 ± 0.016 °C per century. The orange curve represents the contribution made to the blue adjustment curve by breakpoint adjustments only.


The total adjustments in Fig. 26.5 (blue curve) were determined by subtracting the raw monthly data in Fig. 26.1 from the equivalent curve for the sum of the adjusted data. The data was then smoothed with a 12-month moving average which removed over 90% of the noise. The breakpoint adjustment curve in Fig. 26.5 (the orange curve) is just the weighted sum of the breakpoint adjustment curves from the different states (see Posts 18-24). The breakpoints adjustments are determined from the Berkeley Earth station data by subtracting the Berkeley Earth raw anomaly (which is different from my raw anomaly because it uses homogenization) from the Berkeley Earth adjusted anomaly.

The curves in Fig. 26.5 above indicate that summing the Berkeley Earth adjustments to the anomaly data results in a curve that has a trend of gradient +0.295 ± 0.016 °C per century for the period 1901-2010. This in turn amounts to a contribution of over 0.4 °C between 1870 and 2010 as the same trend appears to extend back until at least 1870.


Concluding points

1) According to the raw data (shown in Fig. 26.1), temperatures in Australia may have risen by up to 0.8 °C over the last 60 years, but if so, they are still, at worst, no more than about 0.2 °C above some of the peak values seen in previous centuries. In all likelihood, the temperature rise is probably even less (the data in Fig. 26.1 suggests it may be less than 0.1 °C). Unfortunately a lack of data prior to 1853 precludes any more definitive conclusions than this.

2) The large changes in mean temperature seen over time in the instrumental temperature record for Australia appear to be natural and reversible, and therefore probably occur on a regular basis.

3) The adjustments made to individual temperature records through a combination of breakpoints and homogenization by Berkeley Earth do not appear to completely cancel when multiple station records are averaged. The analysis above suggests that such adjustments could actually add more than 0.4 °C to the warming trend for Australia, which would be more than double the amount that could be attributed to the raw data, as noted in Point 1 above.


4 comments:

  1. A few questions.

    Do you have any reason to prefer raw data over the data that has been cut wherever there are clear discontinuities. I understand Berkeley doesn't homogenise data in the way that other groups do, so if you do prefer the raw data with all its problems (like station moves, changes in instrumentation, time of observation changes etc), could you explain why?

    When using the sparser data set pre 1910, is the sparsity factored in to the uncertainty estimates of trends (eg, 1859 to 2008), or not?

    How many stations are available continuously from 1859 to 1885? Are they well spread out, or grouped together? I choose this period as QLD has no data for then, so I'm wondering what there actually IS.

    What has the IPCC said about Australian temperatures prior to 1910? Specific to what you said in the article, did the IPCC suggest the temps were stable in Australia? I ask this because of you are comparing global with regional, there is no reason expect significant correlation, and especially if the pre-1910 record is sparse or located in an even smaller region (of Australia).

    It seems likely to me that the large swings we see early in the longer record are due to small sample size, and/or poor coverage, and that if coverage then matched that now, the amplitude of the national anomalies would be less. We see the same in the instrumental record, and we see it in the modern record when we sub-sample the global data. Fewer stations means wider swings in temp anomalies.

    ReplyDelete
    Replies
    1. Thanks for your comments, Barry. You pose some good questions, and I will try to answer them.

      In terms of the pre-1910 data, yes it is much sparser than recent data and this is a big problem. But rather than look at Australia as a whole, it is probably better to look at individual states first. I have analysed each state separately (see Posts 18 - 24) and in each case I have indicated the locations of all the stations with more than 40 years and 100 years of data on a map. As you will see, there is very little clustering of stations except in SA and WA, and perhaps around Brisbane in QLD. When it comes to data before 1910, NSW has over 20 records that go back to 1881 or earlier, VCT has seven, SA has two, NT has two, WA has one, Tasmania has one, and QLD has none. While this is not a lot of data, most of it behaves in the same way. The observation of higher temperatures in 1880 and earlier is almost universal. The only exception to this appears to be for data in the middle of the major cities of Melbourne and Sydney, and maybe QLD. Otherwise, there is little evidence of large random swings that need to be averaged out, although before 1870 I would say that the data is more problematic and less reliable.

      I have not compared my results with IPCC, but instead used the results from Berkeley Earth (BE). There are three reasons for this. Firstly, the methodology used by BE is more transparent than other climate groups. Secondly, BE has published trends for every region and country that are easily accessible on its website; the IPCC, as far as I know, hasn't. And thirdly, the BE global results broadly match those of the IPCC and other groups (HadCRUT4, NOAA and GISS). So, by comparing my results with BE, I am by proxy (and climate scientists love proxy data) comparing with the IPCC.

      Finally, I return to your first question: why do I use the raw data and not adjusted data? Firstly, because this should always be the starting point for any data analysis. The raw data is the true data, even if it is corrupted. If it is corrupted, then you need to prove it and explain why. In my view climate scientists have failed to prove it. The implicit assumptions in their analysis seem to be that all temperature trends should vary slowly and smoothly over time, that there should be no abrupt jumps, and that before 1900 all temperatures were fairly stable. Yet the data contradicts this. Most pre-1900 data shows significant natural variation, and many regions exhibit sudden temperature changes. For example, in the 1980s in Botswana there was an abrupt rise of almost 2 °C in many of the temperature records (see Post 38). Similar, but smaller, effects were also seen in data across the border in South Africa and Namibia. If different neighbouring records agree, then you have to trust the data as the datasets are effectively corroborating and endorsing each other. This is why I trust the data more than I trust the attempts of climate scientist to "correct" the data. I am not saying that the raw data is perfect, but I see no evidence so far to indicate that the corrections are necessary or make the data better, and I have analysed every temperature record in the Southern Hemisphere so far, plus many in Europe and North America.

      As for cutting data at discontinuities, here the BE data is very informative. BE uses a statistical test to identify such breakpoints. But if this test is accurate, then I would expect most breakpoints to correspond to a large temperature discontinuity and adjustment. The problem is that many don’t. Whether the proportion that don’t is significant is something that I need to investigate further. Until then, the plan is to study the raw data, and look to see if the pattern of temperature trends is consistent or not from territory to territory. If it is broadly consistent, then I think that throws some doubt over the need to adjust the data. If it isn’t, then I’ll need to think again.

      Delete
  2. What happened to the spike in teh 1930s? newspapers from around the world confirm this was hotter than any time since?

    ReplyDelete
    Replies
    1. What spike? From what I have seen so far, global temperatures in the 1930s were generally cooler than now, and cooler than in the 19th century.

      What the data above shows is how the average monthly temperatures for the whole of Australia have deviated from their long term average. It is possible that a few days in the 1930s at a few locations were the hottest on record, but the average for those months could still have been fairly normal.

      Delete