2011/01/16

Reconstructing HADCRUT from Cosine functions

A "for fun" project - not scientific!

Obtaining the repetitive frequencies:
Using bandpass filters set to a bandwidth of f/150 (where f is the lower frequency of the BP filter) the filter is swept from 0.5year to 150 years looking for amplitude peaks.
As each peak is found the output from the filter is compared to the output of a cosine wave set to the filter output amplitude and to the centre frequency of the filter. The cosine wave is is then phase shifted to align the two waveforms.
This is repeated for as many peaks as are found in the HADCRUT3V record. Outputs from the BP filters that fall into likely TSI frequencies 10, 11, 12, 14, 21 years have their corresponding synthesized outputs nulled. 0.5 year output is also nulled.

TSI is phase shifted, amplitude modified, then it added or subtracted from sum of cosines (shift is +48 months and is subtracted). However from the start of records to a round 1900 a lower error is provided by adding the TSI!

The remaining synthesized outputs are then summed and modified by a multiplier constant (approx 2.69) set by minimising the error between synthesized and original HADCRUT.

No low frequency peaks were found that could fit the rise in temperature. The rise is obtained by curve fitting the HADCRUT data with a 3rd order polynomial:
y = 2.40389E-07x^3 - 1.34093E-03x^2 + 2.49320E+00x - 1.54550E+03

This is added into the synthesised data to provide the output data.

Because the reconstructed data is just the sum of cosines and TSI is relatively unimportant the data can be extended to the future in this case to about 2017.
The plots show 2011 to be cooling to mid year.

The following plots show the effects of
No TSI
TSI subtracted over whole period
TSI added until 1904 and then subtracted to current date








This is a plot of the frequencies and amplitudes used to obtain these plots. Disabled frequencies are shown as Zero amplitude. The output at 11 years represents TSI amplitude







Period yearsamplitudePhase
0.5011650 - artificially set to zero3977.413703
0.8419638890.00675024954.17800684
0.8900402780.010820097107.3740716
1.7683750.00609514863.60020753
1.9372694440.00696353446.86597827
1.9923776040.0050929035.83418163
2.0383916670.00410546540.65715905
2.1070.00465193923.05126836
2.2648577780.00741553929.18466231
2.31393750.00917761269.3504266
2.5417777780.00600900123.33534812
2.6755555560.00786754452.76779754
2.7457888890.00489206410.93806807
2.874550.01041719240.32797096
3.1554833330.00524658275.96412155
3.2691944440.0046675686.630682169
3.5793916670.00830903140.98807579
3.76250.0080482139.685707583
3.7750416670.00821947329.02716352
3.9907583330.00551225755.81368265
4.7031250.0056539496.33244882
5.1270333330.0049196056.188780868
5.2256944440.00482918734.08701089
5.7273611110.0028169156.746840105
5.9865555560.00483971833.42814581
6.2975888890.0039478623.39636614
6.5985888890.0022942663.380314824
7.5860361110.00283688134.9110203
8.4413777780.00273765334.79750828
9.1219722220.00621127617.3232396
10.078483330 - artificially set to zero6.888853744
TSI  WAVEFORM  set at 11 yrs0.006495  THE TSI DATA48 Months
11.452213890- artificially set to zero28.41974504
12.382805560- artificially set to zero22.96039729
14.799166670- artificially set to zero16.96130996
21.178694440- artificially set to zero23.39280935
35.885888890.00178587110.89922244
64.196611110.00278583516.71200197
64.213333330.0027857764.134103705
2508.3333331.56572E-060.302469618
2508.3333331.56572E-060.30372208
2508.3333331.56572E-060.462315105


Problems:
Just what is the filter capability of the BandPass excel plugin? A better version is required.
Need a better trend profile than the 3rd order polynomial
Bandpass filter is not good doeas not detect longer periods - insufficient length of data
=======
Yes - I know that with enough parameters you can match anything with a sum of sines. But this plot gives a reasonable match for 160 years so a few years into the future may also be well matched!
It should also be noted that there is an underlying trend (the polynomial == GHG forcing?) which is always positive. But there are many negative slopes and zero warming areas synthesized.



4 comments:

  1. Cool post, TFP.

    A couple of questions:

    1. What did you use for TSI?

    2. Are there any albedo adjustments?

    3. What frequencies did you detect?

    4. Have you tried fitting the TSI against just the residuals of a linear or exponential trend?

    More to do!

    ReplyDelete
  2. The table (now) above contains all the frequencies seen using this method. Obviously with only 150 years of data detecting long periods is not possible.

    What did I use for TSI. It should have been onthe spreadsheet but isn't. I now cannot find the data I kow I was toying with the idea of interpolating Leif Svalgaards data to month periods. But I then found the data. I will continue to look.

    There are no albedo adjustments

    I have not tried linear+TSi or exponential+tsi. I may try this but I'm not sure what this will proove???

    ReplyDelete
  3. Do the various possible periods between 11:35 and 35:64 then come up as 0 amplitude?

    I am intrigued by the 11 (solar cycle) and the 64 (which is a near miss to my best fit for sines with nls)

    ReplyDelete
  4. Enjoy!
    http://www.fel.duke.edu/~scafetta/pdf/

    ReplyDelete