Guest Post by Willis Eschenbach
I came across a curious graph and claim today in a peer-reviewed scientific paper. Here’s the graph relating sunspots and the change in sea level:
And here is the claim about the graph:
Sea level change and solar activity
A stronger effect related to solar cycles is seen in Fig. 2, where the yearly averaged sunspot numbers are plotted together with the yearly change in coastal sea level (Holgate, 2007). The sea level rates are calculated from nine distributed tidal gauges with long records, which were compared with a larger set of data from 177 stations available in the last part of the century. In most of the century the sea level varied in phase with the solar activity, with the Sun leading the ocean, but in the beginning of the century they were in opposite phases, and during SC17 and 19 the sea level increased before the solar activity.
Let me see if I have this straight. At the start of the record, sunspots and sea level moved in opposite directions. Then for most of the time they were in phase. In both those cases, sunspots were leading sea level, suggesting the possibility that sunspots might affect sea level … except in opposite directions at different times. And in addition, in about 20% of the data, the sea level moved first, followed by the sunspots, suggesting the possibility that at times, the sea level might affect the number of sunspots …
Now, when I see a claim like that, after I get done laughing, I look around for some numerical measure of how similar the two series actually are. This is usually the “R2” (R squared) value, which varies from zero (no relationship) to 1 (they always move proportionately). Accompanying this R2 measure there is usually a “p-value”. The p-value measures how likely it is that we’re just seeing random variations. In other words, the p-value is the odds that the outcome has occurred by chance. A p-value of 0.05, for example, means that the odds are one in twenty that it’s a random occurrence.
So … what did the author of the paper put forwards as the R2 and p-value for this relationship?
Sad to relate, that part of the analysis seems to have slipped his mind. He doesn’t give us any guess as to how correlated the two series are, or whether we’re just looking at a random relationship.
So I thought, well, I’ll just get his data and measure the relationship myself. However, despite the journal’s policy requiring public archiving of the data necessary for replication, as is too common these days there was no public data, no code, and not even a Supplementary Online Information.
However, years of messing around with recalcitrant climate scientists has shown me that digitizing data is both fast and easy, so I simply digitized the graph of the data so I could analyze it. It’s quite accurate when done carefully.
And what did I find? Well, the R2 between sunspots and sea level is a mere 0.13, very little relationship. And even worse, the p-value of the relationship is 0.08 … sorry, no cigar. There is no statistically significant relationship between the two. In part this is because both datasets are so highly auto-correlated (~0.8 for both), and in part it’s because … well, it’s because as near as we can tell, sunspots [or whatever sunspots are a proxy for] don’t affect the sea level.
My conclusions from this, in no particular order, are:
• If this is the author’s “stronger effect related to solar cycles”, I’m not gonna worry about his weaker effect.
• This is not science in any sense of the word. There is no data. There is no code. There is no mathematical analysis of any kind, just bald assertions of a “stronger” relationship.
• Seems to me the idea that sunspots rule sea level would be pretty much scuttled by sunspot cycles 17 and 19 where the sea level moves first and sunspots follow … as well as by the phase reversal in the early data. At a minimum, you’d have to explain those large anomalies to make the case for a relationship. However, the author makes no effort to do so.
• The reviewers, as is far too often the case these days, were asleep at the switch. This study needs serious revision and buttressing to meet even the most minimal scientific standards.
• The editor bears responsibility as well, because the study is not replicable without the data as used, and the editor has not required the author to archive the data.
So … why am I bothering with a case of pseudo-science that is so easy to refute?
Because it is one of the papers in the Special Issue of the Copernicus journal, Pattern Recognition in Physics … and by no means the worst of the lot. There has been much disturbance in the farce lately regarding the journal being shut down, with many people saying that it was closed for political reasons. And perhaps that is the case.
However, if I ran Copernicus, I would have shut the journal down myself, but not for political reasons. I’d have closed it as soon as possible, for both scientific and business reasons.
I’d have shut it for scientific reasons because as we see in this example, peer-review was absent, the editorial actions were laughable, the authors reviewed each others papers, and the result was lots of handwaving and very little science.
And I’d have shut it for business reasons because Copernicus, as a publisher of scientific journals, cannot afford to become known as a place where reviewers don’t review and editors don’t edit. It would make them the laughing stock of the journal world, and being the butt of that kind of joke is something that no journal publisher can survive.
To me, it’s a huge tragedy, for two reasons. One is that I and other skeptical researchers get tarred with the same brush. The media commentary never says “a bunch of fringe pseudo-scientists” brought the journal down. No, it’s “climate skeptics” who get the blame, with no distinctions made despite the fact that we’ve falsified some of the claims of the Special Issue authors here on WUWT.
The other reason it’s a tragedy is that they were offered an unparalleled opportunity, the control of special issue of a reputable journal. I would give much to have the chance that they had. And they simply threw that away with nepotistic reviewing, inept editorship, wildly overblown claims, and a wholesale lack of science.
It’s a tragedy because you can be sure that if I, or many other skeptical researchers, got the chance to shape such a special issue, we wouldn’t give the publisher any reason to be unhappy with the quality of the peer-review, the strength of the editorship, or the scientific quality of the papers. The Copernicus folks might not like the conclusions, but they would be well researched, cited, and supported, with all data and code made public.
Ah, well … sic transit gloria monday, it’s already tuesday, and the struggle continues …
w.
PS—Based on … well, I’m not exactly sure what he’s basing it on, but the author says in the abstract:
The recent global warming may be interpreted as a rising branch of a millennium cycle, identified in ice cores and sediments and also recorded in history. This cycle peaks in the second half of this century, and then a 500 yr cooling trend will start.
Glad that’s settled. I was concerned about the next half millennium … you see what I mean about the absence of science in the Special Edition.
PPS—The usual request. I can defend my own words. I can’t defend your interpretation of my words. If you disagree with something I or anyone has written, please quote the exact words that you object to, and then tell us your objections. It prevents a host of misunderstandings, and it makes it clear just what you think is wrong, and why.

vukcevik says:
January 29, 2014 at 8:57 am
I have superimposed GeoPolar magnetic field data and the by science ‘accepted’ dipole graph … agreement between two appears to be reasonable
The field at any point of the Earth’s surface would be a good fit too. There is nothing special about the ‘GeoPolar’ magnetic field, except perhaps that that field is affected by the largest uncertainties [as the are no paleomagnetic measurements at all up there] and so is the least well known field of all.
lsvalgaard says:
January 29, 2014 at 9:11 am
The field at any point of the Earth’s surface would be a good fit too. There is nothing special about the ‘GeoPolar’ magnetic field,
Dr. S
Let’s data speak for it, here is what Korte data show for your little California town
http://www.vukcevic.talktalk.net/DrS-town.gif
On that bet you would have lost your shirt.
GeoPolar MF is possibly most accurate data for Earth’s dipole during last 7,000 years.
I could email you file if you wish.
vukcevic says:
January 29, 2014 at 9:39 am
Let’s data speak for it, here is what Korte data show for your little California town
You do not say what you plot. The total field strength would be best.
GeoPolar MF is possibly most accurate data for Earth’s dipole during last 7,000 years.
‘possibly’? There are no measurements up there. If you would try to get a measure of the dipole strength you should plot the field at the corrected geomagnetic latitude of 90 degrees, not the meaningless geographic pole.
Have you a file of numbers, which you think is more accurate?
“There are no measurements up there.”, Doesn’t matter much, see http://www.ngdc.noaa.gov/geomag/WMM/data/WMM2010/WMM2010_Report.pdf
Plotted are values from CALSK7K (I assume radial F) Br N & S absolute values.
Current 90 degrees geomagnetic latitude may not be a good representation due to poles movements. For the last 2Ky the movements were more or less cantered around 90N (see http://dourbes.meteo.be/aarch.net/linford.pdf page5 ).
I have a program that will scan gumf file for a set Lat/Long box for max values, but from Korte’s file one might not be able to work it out by scanning data, since GMF pole is not same as the max Br. I suppose by spherical triangulation of declination data from two locations would it be possible.
If you can suggest alternative numbers (possibly for both N & S ) I’ll have a go at it and compare to the accepted dipole reconstruction.
vukcevik says:
January 29, 2014 at 11:55 am
Current 90 degrees geomagnetic latitude may not be a good representation due to poles movements
No location on the surface is a good representation of the dipole movements. All locations are bad for this, although some might be worse than others. Your whole exercise is meaningless as the surface field is not what cosmic rays and the solar wind see.
9:11 am The field at any point of the Earth’s surface would be a good fit too.
11:55 am All locations are bad for this, although some might be worse than others.
Well, well, well!
A scientist has file of raw C14 or 10Be data going back 6-7 millennia. To work out the heliospheric mf he has to subtract the Earth MF dipole.
Where the a file for the Earth MF dipole data is to be found?
‘Sine wave’ like curve from your link
http://www.leif.org/research/CosmicRays-GeoDipole.jpg demonstrably is no good !
If there is a raw C14 or 10Be data going back 2 or more millennia ?
let’s do the science properly and compare
I’ve got the data, do you?
As Willis implies : No data, no contest !
vukcevic says:
January 29, 2014 at 2:11 pm
Jeez, Vuk, lose the attitude, it just makes you look arrogant.

There’s a file here of five different estimates of the changes in the Earth’s geomagnetic dipole based on a variety of proxies proxies. I found them by actually following up the links that Leif posted upthread … you might consider doing the same.
I’ve graphed them up:
Happy now?
w.
vukcevic says:
January 24, 2014 at 1:49 pm
Usoskin & Korte:
We conclude that changes of the regional tropospheric ionization at midlatitudes are defined by both geomagnetic changes and solar activity, and none of the two processes can be neglected.
———————————————–
Nope can’t forget the regional tropospheric ionization at midlatitudes. Another way for GCR to enter the system. Just give them a little push out their homestead..pulses and waves or changes in solar winds.
I wish this movie would go a little slower though. You easily miss the polar, cusp and equator action ..
The constant in motion radiation belts remind me of magnetopause simulations..and the constant motion that creates WAVES………………………………………………………………………………………………………
Solar Anomalous and Magnetospheric Particle Explorer (SAMPEX)
Published on May 18, 2013
Movie of the changing radiation belts as measured by SAMPEX/LICA from January 1, 1998 to March 1, 2005
Thanks Vuks and Dr. S., for links and discussion.
Those pesky multi linked videos over there.
Last video good but not the intended video..
vukcevik says:
January 29, 2014 at 2:11 pm
If there is a raw C14 or 10Be data going back 2 or more millennia ?
Going back 10 millennia. And easy to find
let’s do the science properly and compare
You can’t do science and knowledgeable people have already done that. One problem is that the data is not good enough to settle some of the questions, e.g. how much of the variation is due to the Sun, the earth’s magnetic field, and to climate. These questions [on which I am sort of an expert] are active areas of research, see e.g. http://www.leif.org/research/Svalgaard_ISSI_Proposal_Base.pdf and http://www.leif.org/research/Long-term-Variation-Solar-Activity.pdf
Carla says:
January 29, 2014 at 6:16 pm
Nope can’t forget the regional tropospheric ionization at midlatitudes. Another way for GCR to enter the system.
You [as does also Vuk] misunderstand the science. The tropospheric ionization has nothing to do with the Radiation belts or waves or any of what else you rave about..
Usoskin & Korte: We conclude that changes of the regional tropospheric ionization at midlatitudes are defined by both geomagnetic changes and solar activity
All they are saying is the well-known fact that the GCRs that ionize the troposphere are modulated both by the slow, long-term changes of the main magnetic field of the Earth and by the activity of the solar magnetic field in the [outer] heliosphere. Nothing new there.
@ur momisugly Willis
I knew I’ll get myself into trouble. Me “arrogant”, the man who claims to have ‘corrected’ all known and yet unknown science, surely you’re joking …
Dr. S and I, or our ‘conversations’ go long way back, I accept his put downs, he accepts my ignorance, and even managed to improve on the spelling of my name. We are kind of friends and occasionally exchange emails. Unfortunately, if it was not for Dr. S you wouldn’t ever heard of me.
Thanks for the graph, I actually did look at the files:
Korte’s geomagnetic dipole ( http://earthref.org/ERDA/973/ ) is interesting
http://wattsupwiththat.com/2014/01/21/sunspots-and-sea-level/#comment-1551530
and plotted data against an earlier Korte file, as shown here, but didn’t link to it!
http://www.vukcevic.talktalk.net/MF-6D.htm
Differences are of similar order, except that in my set of calculated data (Korte’s radial geomagnetic model) the data go further back.
What I was actually looking for is the raw numerical data filesfor C14.and or 10Be and data for dipole which were used before Korte’s compilation (2009) for eliminating Earth’s field.
It does help if one is more precise with language, but being vague it does help one dig himself out of a deep hole, and some are in it more often than the others.
@ur momisugly Dr. Svalgaard
Thanks doc, for the links, have seen one or two before but I’ll check it out again.
vukcevic says:
January 30, 2014 at 1:28 am
What I was actually looking for is the raw numerical data files for C14.and or 10Be and data for dipole which were used before Korte’s compilation (2009) for eliminating Earth’s field.
The raw data comes from a variety of sources as detailed in the reports in PNAS http://www.leif.org/EOS/PNAS-2012-Steinhilber.pdf and http://www.leif.org/EOS/PNAS-2012-Steinhilber-Appendix.pdf
Figures S1 to S7 shows ‘semi-raw’ data [top panels] which you can digitize yourself. The data presented has already been corrected for the varying strength of the Earth’s dipole. You cannot duplicate that by simple correlation as it is necessary to apply a complicated physical model [as described in the papers] and actually compute the correction. Using the surface field is meaningless as that is not what the cosmic rays see. You must calculate from the spherical harmonic coefficients and a model of the magnetosphere what the magnetic field in space looks like and then calculate the orbits followed by the cosmic rays for all energies and integrate the modulation in order to correct the observed flux.
Some of the ‘raw’ data can be found here: http://www.ncdc.noaa.gov/data-access/paleoclimatology-data/datasets/climate-forcing
Real science is not just comparing time series. There is physics in between and that physics must be understood and applied.
vukcevik says:
January 30, 2014 at 1:28 am
even managed to improve on the spelling of my name.
Your name is on WUWT’s ‘black list’ that causes a comment to go into ‘moderation’. Misspelling it [as done in this comment] avoids that…
vukcevik says:
January 30, 2014 at 1:28 am
What I was actually looking for is the raw numerical data files for C14.and or 10Be and data for dipole which were used before Korte’s compilation (2009) for eliminating Earth’s field.
The raw data comes from a variety of sources as detailed in the reports in PNAS http://www.leif.org/EOS/PNAS-2012-Steinhilber.pdf and http://www.leif.org/EOS/PNAS-2012-Steinhilber-Appendix.pdf
Figures S1 to S7 shows ‘semi-raw’ data [top panels] which you can digitize yourself. The data presented has already been corrected for the varying strength of the Earth’s dipole. You cannot duplicate that by simple correlation as it is necessary to apply a complicated physical model [as described in the papers] and actually compute the correction. Using the surface field is meaningless as that is not what the cosmic rays see. You must calculate from the spherical harmonic coefficients and a model of the magnetosphere what the magnetic field in space looks like and then calculate the orbits followed by the cosmic rays for all energies and integrate the modulation in order to correct the observed flux.
Some of the ‘raw’ data can be found here: http://www.ncdc.noaa.gov/data-access/paleoclimatology-data/datasets/climate-forcing
Real science is not just comparing time series. There is physics in between and that physics must be understood and applied.
Thanks doc, If I wasn’t ‘block-head ignorant’, wouldn’t know about latest Steinhilber’s paper link.
In my inimitable manner, I like to interfere with the ‘settled science’ and found these two oddities
http://www.vukcevic.talktalk.net/Stein-Vuk.htm
I expect you might say: So what? but would appreciate more elaborate comment
I will look into some of the ‘paleoclimatology-data/datasets/climate-forcing’ data, and then I will be off to my winter bolt-hole, and write a bit more about http://www.vukcevic.talktalk.net/Ap-NHT.htm
p.s
you are welcome to use
any spelling you choose.
vukcevik says:
January 30, 2014 at 6:06 am
In my inimitable manner, I like to interfere with the ‘settled science’ and found these two oddities
Since it is meaningless to compare with the surface field the ‘so what’ reply seems appropriate.
I am sure you can do better than that, perhaps Steinhilber needs to go back to his data.
If he does so, Steinhilber et al would find that the ‘paleo’ TSI is as flat as pancake (except tiny 11 year mod) as certain Dr. Svalgaard of Stanford and of the WUWT, has been telling us for some time now; ahh… but some people will never learn, let alone do the science.
vukcevic says:
January 30, 2014 at 1:28 am
Now you’re not only being arrogant, you are accusing me without quoting my words, which is the action of a willful, spiteful, spoiled child who thinks he can just ignore polite requests, viz:
w.
Willis Eschenbach says:
January 30, 2014 at 11:33 am
…………
Jeez, Vuk, lose the attitude, it just makes you look arrogant.
I really didn’t object to it, I heard a lot worse from elsewhere.
Sorry, it appears to be a misunderstanding, my response was meant solely as a bit of a ‘self deprecating humour’ with no intention of malice, etc .
Nevertheless, I offer an apology.
All the best.
Willis,
Now I see where the misunderstanding comes from, after reading your comment for the third time.
I knew I’ll get myself into trouble. Me “arrogant”, [from?] the man who claims to have ‘corrected’ all known and yet unknown science, surely you’re joking …
I was describing my ‘scientific endowers’ not yours.
Word ‘from’ it was not and never meant to be there, it was ‘ me’ and not you.
Oh, never mind.
Have a good day.
vukcevic says:
January 30, 2014 at 12:03 pm
my response was meant solely as a bit of a ‘self deprecating humour’ with no intention of malice,
If you were sincere there would have been no need for quote marks around self deprecating humour. There is but a thin line between the humorous and the ridiculous and it is not clear in each case on which side of the line you are at that time. Better simply to stick to your continuing science education.
I often put quote marks for something that is a colloquial phrase, which -self deprecating humour- is . Well you can make of it anything you like, At first I didn’t read italics since I know what I written and what I meant, only at third reading I did realise [from?] was inserted. . If you read whole post again
http://wattsupwiththat.com/2014/01/21/sunspots-and-sea-level/#comment-1554317
you can see that I was talking about myself, and that is how it was meant. It is not my habit to comment on anyone personally, regardless what others may or may not say or do. I shall leave it at that.
vukcevik says:
January 30, 2014 at 1:28 am
Me “arrogant”, the man who claims to have ‘corrected’ all known and yet unknown science
I think that claim is unfounded and plain silly. I think you actually [in your heart] mean what you claim which is quite sad.
vukcevic says:
January 30, 2014 at 12:03 pm
My apologies as well if I misunderstood your words.
w