Guest Post by Willis Eschenbach
I had one of my investigations take a curious turn recently. I was going to use as my springboard the most interesting 2015 paper entitled The Albedo Of Earth, by Graeme L. Stephens et al. However, a strange thing happened along the way. I got to thinking about their Figure 5, in particular Panel (a).:
Figure 1. This shows Figure 5(a) of Stephens2015. ORIGINAL CAPTION: Annual cycles of (a) the globally averaged albedo … The solid curves are for all-sky fluxes, and the dashed curves are the clear-sky fluxes. The error bars represent the interannual variability. … Annual means of all quantities have been subtracted.
To start with, a few definitions. The “albedo” of a planet, moon, or other celestial object is a number from 0.0 to 1.0 that measures the fraction of solar radiation that is reflected from the surface of the object. It’s often given as a decimal fraction (e.g. 0.30), although I prefer it as a percentage (e.g. 30%). The albedo of the earth is about 0.29, meaning 29% of the sunlight is reflected back to space.
But somehow, their graph didn’t look right to me. I’ve looked at a lot of graphs of the planetary albedo. It’s a curious curve. The maximum albedo occurs on the summer and winter solstices in December and June, when one of the poles are maximally pointed at the sun. At that point, the most ice-covered area is exposed, which makes for a high albedo. In addition, because the snow- and ice-covered area is much larger in the northern hemisphere than in the southern, the December albedo is higher than the June albedo.
And on the other hand, on the equinoxes the sun sees mostly ice-free zones, and so the albedo is lowest around March and September.
The problem is that this means that the two peaks in albedo will be somewhere around the winter and summer solstices (December 21 and June 21) with the two minimum albedo measurements between the two … and their graph shows nothing of the sort. So to see if their graph was correct, I made my own graph using the CERES data. Here’s that result.
Figure 2. Plot of the albedo from the CERES dataset, March 2000 to February 2014 (14 years).
Note that as expected, the peaks are around the solstices (June and December) with the high peak in December.
To determine which version is correct, the Stephens2015 version or my version, let me offer the following graph from the Encyclopedia of Climate, by Gerald North et al. It shows, not the average, but a single year of the albedo variations of the CERES dataset. However, the albedo varies little from year to year.
Figure 3. The month-by-month albedo cycle as shown in the Encyclopedia of Climate.
Note that the form of the albedo record in the North book is identical to that of my average shown in Figure 2—it peaks around the summer and winter solstices, is highest in December, and is lowest around the two equinoxes. This is in complete contradiction to the Stephens2015 results shown in Figure 1.
This becomes a significant issue because Stephens et al. use the same graph in a later section of their paper to show how different their results are from the output of many climate models. This is shown in Figure 4 below.
Figure 4. This shows Figure 11(a) of the Stephens 2015 paper. ORIGINAL CAPTION: The global mean annual cycle of (a) TOA albedo … The solid lines are CERES observations taken from Figure 5, and the colored lines are 10 year average seasonal cycle of individual CMIP5 models, and the dashed lines are the multimodel mean seasonal cycle.
Based on this graph, the authors of Stephens2015 reasonably say:
From the comparisons presented in these figures, it is evident that models and measurements differ in potentially important ways.
Unfortunately, when we put in the correct values for the albedo variations, a very different picture emerges:
Figure 5. Actual average albedo variations from the CERES data (thick red line) overlaid on Figure 4 (Stephens2015 Fig. 11).
As you can see, while the model average (dotted line) still has problems compared to the red line that shows the CERES albedo variations, at least they generally get the peaks around the solstices (21st of June & December), the low points at the equinoxes, and the larger peak in December.
Now, I noticed this important error back at the end of May, and I was in a bit of a quandary regarding what I should do. In the event, I noted that Peter Webster, an associate of Judith Curry, was one of the authors. Given my respect for Dr. Curry, I didn’t want to blow the whistle on Peter, and I thought I might actually be able to take another path. I didn’t have his email address, so I wrote to Dr. Curry and sent her the above analysis, and I asked her to pass it on to Peter Webster, which she kindly and promptly did.
Dr. Webster was very good about the matter. He responded to me immediately, and said that he had passed my email on to Dr. Stephens because he is the lead author of the paper. So I waited.
After a short while I asked what was happening. Peter said that Dr. Stephens was in England, it would be a week or so. When nothing transpired, I got back in touch with Peter three weeks later. He said that like me, he had heard nothing from Dr. Stephens. Finally, after almost a month and a half had expired with no answer from Dr. Stephens, Peter said that I might as well go ahead and publish, and that I should also send a formal note to Review of Geophysics.
So I am going to take him up on the first half of his suggestion, and any fault in my taking that step is mine alone, not his. I gotta say, this is very frustrating. I tried to do the right thing, even had the backing of one of the co-authors, and I got exactly zip in return. Note that I do not fault Peter Webster in this matter in any way. He has been most responsive and supportive throughout, but his lead author is not replying to the question, which leaves Peter with no options. Ah, well. I’d hoped that a nudge would be as good as a wink to a blind horse, but I guess sometimes that’s not enough, you need a nudge plus a baseball bat.
However, I don’t want to write a dang letter to Rev Geo. I always feel like I have to give myself a lobotomy to write in the black-letter long-paragraph obscurantist style favored by the journals …
So I’ve chosen to make my first (and perhaps only) move by reporting the matter here. At least now, all of the modelers that Stephens claimed were wrong in important ways will know that a) Dr. Stephens was wrong on this particular point, and b) he declined my offer to reveal his mistake himself and correct it in his preferred manner.
We’ll see what comes after that.
Let me add that I think I have acted in this matter with what passes for my best manners. For example, here is what I said to Peter in my earliest emails:
So please take as much time as it needs (within reason) for you to both verify the error, and then decide how you wish to proceed. I am well aware of how difficult it is to be publicly found in error. I’m one of the few bloggers out there with a post entitled Wrong Again.
Let me assume for a moment that your graph is in error, based on my looking at lots and lots of albedo graphs, and based on physics (higher albedo at the solstices when the poles are most exposed to the sun), and based on the Encyclopedia of Climate as I showed, and based on my own calculations shown in Figure 2, and based on the vague but visible similarity to the models.
If that is indeed the case, I’m more than willing to have you and Graeme make the initial announcement, assuming of course that I’m given appropriate credit. I’ll reserve the right to write it up for WUWT if that comes to pass, but you are free to pick the time, place, manner, and content of the announcement.
Like I say, I’ve been there under the spotlight. So I’ll do whatever I can to make it work as best it can for you.
I don’t think I could have been more supportive than that … and I got nothing. People often say something like “Hey, why don’t you write to the authors instead of simply posting your results up at WUWT”. This is one of the many examples of why I’ve generally given up on that approach … it often doesn’t work, and when it doesn’t, it is the source of much delay and frustration.
Despite this error, Stephens2015 is still an interesting paper, don’t let me put you off reading it. I’m sorry that it’s marred by this single problem.
In any case, TGIF, work week’s over. I’ve been up on a scaffolding working on second story windows the last three days, and although I don’t mind the heights, it is still tiring to have to spend the day with the pucker factor up somewhere in the low seventies, safety line or not … I’ll be glad when this part of the job is done.
Regards to all, and my wish is that you get to spend the weekend doing something other than working high up on a scaffolding with your life depending on some overgrown piece of string tying you to safety …
w.
Don’cha Know: If you disagree with someone please have the courtesy to quote the exact words they used. That way we can all be clear exactly what it is that you are objecting to.
Discover more from Watts Up With That?
Subscribe to get the latest posts sent to your email.
I wonder what this prediction will do to Ceres calculations of albedo in 15 years time
http://www.telegraph.co.uk/news/science/11733369/Earth-heading-for-mini-ice-age-within-15-years.html
“Earth heading for ‘mini ice age’ within 15 years
River Thames could freeze over in 2030s when Northern Hemisphere faces bitterly cold winters, scientists say”
Living in Scotland I might live long enough to walk on the Firths of Clyde and Forth.
I wonder what Leif Svalgaard thinks about this. Should be up his street.
Leif, hvad mener du om denne forudsigelse?
Could it be he used the right graph for the Northern hemisphere but mislabelled it as “Global”?
Willis – For us lay persons. How does the CERES Albedo calculation relate to total solar energy received by earth? Isn’t most of the energy in the infrared spectrum?
No it’s not it’s mostly visible.
Look for the bright yellow white object in the sky; That’s where it’s coming from.
There were long discussions about this paper in March on Judith Curry’s and ATTP’s blogs
http://judithcurry.com/2015/03/10/the-albedo-of-earth/
and
https://andthentheresphysics.wordpress.com/2015/03/11/new-albedo-paper/
Willis,
There are two different definitions of albedo.
The “Bond albedo” is weighted according to intensity of insolation. That is the albedo you have plotted.
The “geometric albedo” is weighted by surface area, without regard to the angle of the sun. That should peak around the equinoxes, for reasons given above by Richard Verney.
Stephens et al. do not seem to say which one they use, but it appears that they use the latter. If that is standard for the type of work in their paper, and in reporting albedo from models, then there might well be nothing wrong. They still *ought* to be clear about which definition they use, but from my experience such lack of clarity is not surprising.
If the difference is only a matter of definition, then I am not surprised that Stephens did not reply to you. He is under no obligation to explain things to you and he is obviously very busy (look at the title page: 3 institutions separated by 8 time zones).
If the issue is just a matter of definition, it is disconcerting that Webster could not point out the misunderstanding. Probably due to overspecialization and diffusion of responsibility.
I should add that I am not claiming that Stephens et al. have not made an embarrassing mistake. That may well have inadvertently calculated a geometric albedo an incorrectly compared it to a Bond albedo. I only point out that it *might* not be a mistake.
Certainly the context of the paper implies that “albedo” means “Bond albedo”.
Mike,
When talking about climate, albedo means ‘bond albedo’. I can imagine no circumstances where ‘geometric albedo’ would be relevant, but then maybe someone will surprise me by suggesting one.
MikeB,
“When talking about climate, albedo means ‘bond albedo’.”
I would think so too.
“I can imagine no circumstances where ‘geometric albedo’ would be relevant”
I can, but only as a step towards another result. You can find numbers for the albedo of the ocean, or vegetation, or ice; those are geometric albedo.
Let’s say that you want to compute global albedo from either satellite data or surface properties. You would divide the surface into a grid, obtain or estimate a value of geometric albedo for each grid cell, then average over the surface to get the albedo. A simple average would give you global geometric albedo and an insolation weighted average would give Bond albedo.
Off hand, I don’t see any use for globally averaged geometric albedo. But that seems to be what Stephens et al. provide. They refer to the values in Figure 5 (a), (b), and (c) as “globally averaged albedo” and “global mean albedo”, suggesting simple averages. And if those were Bond albedo, then the reflected fluxes in 5 (d), (e), and (f) would be just the albedo times insolation. But it does not look like that is the case. But the purpose in plotting geometric albedo is unclear.
What is more significant is whether the CMIP5 model output in Figure 11 is also geometric albedo and, if so, why they compared that rather than reflected flux.
Well, you could do something complicated like that, but it would not be geometric albedo.
Geometric albedo is simply the ratio of how much a planet reflects back towards a light source compared to a Lambertian reflector. What use it is, I don’t know (something to do with astronomy).
Bond albedo is the proportion of insolation a planet and its atmosphere ‘reflect’ back to space, and that is the one we want.
Oh…and geometric albedo can be greater than 1.
Thanks guys.
I had not picked up on the distinction and seperation of bond albedo and geometric albedo.
I had assumed that when one assesses the albedo of the Earth and how it varies throughout the year one looks at changes in geometric albedo (changes in the area of ice & snow cover/extent, changes in vegetation, changes in patterns of cloudiness) and overlay on top of that changes in solar insolation/changes in the face of the planet being presented to the sun (the latter will of course include changes brought about by changes in reflectivity of the oceans as more or less ocean is being presented to the sun).
MikeB,
“Well, you could do something complicated like that, but it would not be geometric albedo.”
But there is nothing complicated about calculating geometric albedo. It is Bond albedo that is complicated. The latter could be measured directly if you are far enough away. LEO satellites are too close for that.
Well I’ve always taken albedo as the total fraction (percentage) of the solar spectrum energy arriving at earth that is redirected back to space. No need to weight anything for anything. And note it is what fraction goes back out with no wavelength shift.
In addition to refractive scattering albedo by clouds (it’s not reflection) there is also a good deal of absorption by the H2O and other GHGs in the atmosphere, which lowers surface level insolation. The blue sky Raleigh scattering is of course a part of the albedo contribution.
That radiant energy is not re-emitted at solar spectrum wavelengths but as a LWIR thermal spectrum dependent on the cloud (water) temperature.
Less than half of that can ever reach the surface to contribute to surface heating.
I always feel like I have to give myself a lobotomy to write in the black-letter long-paragraph obscurantist style favored by the journals …
“it is always felt.”
Oh, very good, cracked me up entirely. I forgot about the required use of the passive case.
w.
Don’t you mean the required use of the passive case was forgotten by the writer? HAHAHAHA
Good luck.
Joe Born wrote above.
“Objecting to the mode of the criticism in order to evade addressing its substance is the last refuge of a charlatan.”
IMO that should be somewhere in the WUWT header.
Fig. 02 and 03 indicate a higher albedo at perihelion and at the highest solar insolation mainly for the SH due to the eccentricity. The total energy received from the sun is at minimum versus the reversed case, e.g. NH.
Does this means, combination with an actual decline in obliquity, a cooling trend is enhanced?
Aha! This explains it:
http://www.fnal.gov/pub/today/images11/ROW042111figure1.jpg
The chart in Figure one of the article appears to follow sea ice extent. Sea ice extent peaks in March in NH, and in September in the SH.
I am surprised the All Sky and Clear Sky curves are so close – this implies that albedo changes due to clouds is relatively low.
This is the explanation from the study. It looks fine to me.
“This seasonal cycle behavior is set in part by the seasonal cycle of surface albedo
which has a boreal spring-time maximum resulting from the reflection off the brighter snow-covered land
surfaces between 30°N and 60°N (Figure 5c). The second weaker maximum in the boreal fall season is
influenced by reflection from midlatitude clouds of the SH (Figure 5b)”
Maybe I’m missing something willis, but the graph is clearly labeled as all-sky and clear-sky albedo… it doesn’t sound like they’re using ground-features at all. (no, haven’t read the actual study.)
“Annual means of all quantities have been subtracted.”
Would subtract out things like Antarctic/Arctic ice, right?
That would be the static ice. Ice that grows/shrinks would still be included.
ie. This is a graph of albedo anomoly, not albedo itself
Be polite and get no response. Send the note off to a peer reviewed journal, and have the believers circle the wagons and allow something pertinent to get buried until the opposition can safely take credit for it themselves. Pick your poison. Etiquette and manners are extra-scientific values, and as such are a burden to the scientific enterprise. Send the impolite note to the journal–the response often exposes “feet of clay.”
From an answer to richard verney:
“The albedo peaks represent some kind of mix of winter snow, winter ice, winter clouds, and solar angle”
I wonder what the solar angle has to do with albedo?
Svend Ferdinandsen replying to richard verney)
Snow and ice fields do not generally show any effect of solar angle on measured albedo, but then again, almost all of the snow and ice albedo studies I read for ice and snow fields on land deliberately measure near noon each day over a series of days deliberately to prevent “solar angle” from affecting their results.
Assume their is little solar angle effect because the reflection randomly bounce around under the first few cm of ice and snow then bounce up and out of the ice. The number of bounces under the surface and the random crystal structure of the near-surface ice in a natural snow field don’t generate predictable angle-of-independence reflection changes from the sun rays.
Water, on the other hand, has a very, very large albedo changes for direct solar radiation that is controlled affected by both wind speed and solar elevation angles under 33 degrees. Water’s (open ocean measurements) are as high as 0.38 – 0.52 at solar elevation angles under 4-7 degrees.
Above 33 degrees from the horizon, water has a near-constant albedo = 0.066 for both direct and indirect radiation.
Anyone who has looked down a crevasse knows that bluegreen radiation propagates for meters down into the ice.
After a few hours in the sun the snow surface has melted and refrozen, and is now a transmitting surface that lest solar radiation in and traps it in the ice by TIR. As in the oceans, the red radiation is absorbed first, and the blue green last, hence the crevasse illumination.
After about 72 hours fresh snow is no longer highly reflective, hardly more so than green grass.
Snow looks very “bright” to the eyes, because it is essentially a Lambertian reflecting surface so the Radiance and Luminance are the same in all directions. The cosine fall off in Intensity is compensated by the cosine reduction in perceived emitting area so the radiance / Luminance is independent of viewing angle.
But it is all a mirage. It is brighter than other surfaces around it, and the eye is not a good “light meter” since it has a very powerful automatic gain control.
Very nice that that “bright” Antarctic ice patch is right next to the black outer space patch. Very nice trick and a very untypical earth photo with unusually small northern cloud cover.
g
I see the point, to discriminate between general albedo and albedo relative to the Sun, but then it should be stated what kind of albedo we talk about.
Just think about the albedo of a mirror? Unless you are perpendicular to the mirror the albedo might be 0.
Albedo is not reflection coefficient.
Nothing. Albedo is a single number for the entire earth, and knows nothing about what direction the sun is looking or pointing.
Well Albedo is NOT reflection coefficient.
The latter does depend on incidence angles and such like.
Albedo simply means what percentage of incident solar spectrum energy is returned to space without alteration by any earth effects. That is alteration in the sense of any change in photon energy or frequency.
1-albedo is the total solar energy that is left to interact with the earth system.
It matters not by what process the incident solar photons are redirected to space, whether by Raleigh scattering or refractive scattering from water droplets or ice crystal, or reflection from solid or liquid surfaces.
However, once those solar photons get absorbed in some process by some physical material, that energy ceases to be a part of albedo.
Talking about geometry with relation to albedo is just nonsense.