Guest Post by Willis Eschenbach
I had an idea a couple days ago about how to estimate cloud feedback from observations, and it appears to have panned out well. You tell me.
Figure 1. Month-to-month change in 5° gridcell actual temperature ∆T, versus gridcell change in net cloud forcing ∆F. Curved green lines are for illustration only, to highlight how many of the datapoints fall outside those lines in each of the four quadrants. Results have been area-weighted, giving a slightly smaller slope (-1.7 W/m2 per degree) than initally reported (- 1.9 W/m2 per degree). Data colors indicate the location of the gridcell, with the Northern hemisphere starting with blue at the far north, slowly changing to yellow and to red at the equator. From there, purple is southern tropic, through pink to green for the farthest south latitudes. Updated.
Cloud feedback is what effect the changing clouds have if the earth warms. Will the clouds act to increase a warming, or to diminish it? The actual value of the cloud feedback is one of the big unknowns in our current understanding of the climate.
The climate models used by the IPCC all say that as the earth warms, the clouds will act to increase that warming. They all have a strong positive cloud feedback. My thunderstorm and cloud thermostat hypothesis, on the other hand, requires that the cloud feedback be strongly negative, that clouds act to decrease the warming.
My idea involved the use of what are called “gridded monthly climatologies”. A monthly climatology is a long-term month-by-month average of some climate variable of interest. “Gridded” means that the values are given for each, say, 5° latitude by 5° longitude gridbox on the surface of the planet.
My thought was to obtain the monthly actual temperature gridded climatology. This is the real temperature “T” as measured, not the anomaly. In addition, I would need the gridded net cloud forcing “F” from the ERBE (Earth Radiation Budget Experiment) data. Net cloud forcing is the balance of how much solar energy the clouds reflect away from the earth on the one hand, and on the other, how much the same clouds increase the “greenhouse” downwelling longwave radiation (DLR). Net cloud forcing varies depending on the type, thickness, altitude, droplet size, and color of a given cloud. Both positive and negative cloud forcing are common. By convention, positive net cloud forcing (e.g. winter night-time cloud) is warming, while a negative net cloud forcing (e.g. thick afternoon cumulus) is cooling.
Remembering that a cloud feedback is a change in net forcing in reference to a change in temperature, I took the month to month differences of each of the two climatologies . I did this in a circular fashion, each month minus the previous month, starting from February minus January, around to January minus December. That gave me the change in temperature (∆T) and the change in forcing (∆F) for each of the twelve months.
The ERBE satellite only covers between the Arctic and Antarctic circles, the poles aren’t covered. So I trimmed the polar regions from the HadCRUT absolute temperature to match. Then, the HadCRUT3 absolute temperature data are on a 5° grid size, while the ERBE satellite data is on a 2.5° grid. Since the grid midpoints coincided, I was able to use simple averaging to “downsample” the satellite cloud forcing data to correspond with the larger temperature gridcell size.
The results of the investigation are shown in Figure 1. The globally averaged cloud feedback is on the order of -1.9 watts per square metre for every one degree of monthly warming.
This result, if confirmed, strongly supports my hypothesis that the clouds act as a very powerful brake on any warming. At typical Earth surface temperatures, the Stefan-Boltzmann equation gives about five watts per square metre (W/m2) of additional radiation per degree. That is to say, to warm the surface by 1°C, the amount of incoming energy has to increase by about 5 W/m2. This, of course, means that if there were no feedbacks, a doubling of CO2 (+3.7 watts per square metre per the IPCC) would only cause about 3.7/5 or about three-quarters of a degree of warming. The models jack this three-quarters of a degree up to three degrees of warming by, among things, their large positive cloud feedback.
But this analysis says that the cloud feedback is strongly negative, not positive at all. As a result, a doubling of CO2 could easily cause less than eight-tenths of a degree of warming. If the cloud negative feedback is actually -1.9 W/m2 per degree as shown above, and it were the only feedback, a doubling of CO2 would only cause half a degree of warming …
If confirmed, I think that this is a significant result, so I put it up here for people to check my math and my logic. I’ve fooled myself with simple mistakes before …
Code for the procedures and data is appended below.
All the best,
w.
PS – please, no claims that the “greenhouse effect” is a myth or that DLR doesn’t exist or that DLR can’t transfer energy to the ocean. I’m beyond that, whether you are or not, and more to the point, there are plenty of other places to have that debate. This is a scientific thread with a specific subject, and if necessary I may snip such claims (and responses) to avoid thread drift. If so, I will indicate such excisions.
NOTE: The slope of the trend line in Figure 1 is now properly area-adjusted, making the following section andFigure 2 superfluous. .[UPDATE] I’ve gone back and forth about whether to area-average. The problem is that the gridcells are not the same size everywhere. The usual way to area-average is to multiply the data by the cosine of the mid latitude, so I have done that.
Figure 2 shows the area-adjusted version. Still a significant negative feedback from clouds, but smaller than in the non-adjusted version.
FIGURE 2 REMOVED
Figure 2. Area adjusted cloud feedback. Note the lower estimate of the cloud feedback, a bit smaller than my initial estimate. Color of the dots indicates latitude, ranging from blue at the furthest north through cyan to the equator, then in the southern hemisphere through yellow to red at the furthest south.
Note that we still see the same form in the four quadrants. It is still rare for a large temperature drop to be associated with anything but a rise in the cloud forcing.
I’m still not completely happy with this method of area-adjusting, because it adjusts the data itself. But I think it’s better than no area-adjusting at all. The best way would be to convert both of the datasets to equal-area cells … but that’s a large undertaking and I think the final result won’t be much different from this one.
[UPDATE] Here’s the two hemispheres:
Figure 3. Northern Hemisphere Cloud Feedback. Color of the dots indicates latitude, ranging from blue at the furthest north through yellow in the subtropics, to red at the equator.
Figure 4. Southern Hemisphere Cloud Feedback. Color of the dots indicates latitude, ranging from green at the furthest south through pink in the subtropics, to purple at the equator.
[UPDATE] To better inform the discussion, I have made up the following maps of the variables of interest, month by month. These are the monthly absolute temperatures T, the monthly net cloud forcings F, and the month by month changes (deltas) of those variables, ∆T and ∆F.
Figure 5. Absolute temperature (T)
Figure 6. Net Cloud Forcing (F)
Figure 7. Change in absolute temperature (∆T)
Figure 8. Change in net cloud forcing (∆F)
APPENDIX: R code to read and process the data (not including the updated charts). I've tried to keep wordpress from munging the code, but it likes to either put in or not put in carriage returns.
===================================
# data is read into a three dimentional array [longitude, latitude, month]
diffannual=function(x){# returns month(t+1) minus month(t)
x[,,c(2:12,1)]-x
}
# rotates the circle of months by n
rotannual=function(x,n){
if (n!=0) {
if (n>=0){
x[,,c((n+1):12,1:n)]
} else {
x[,,c((13+n):12,1:(12+n))]
}
} else
x
}
#_averages_2.5°_gridcells_into_5°_gridcells,_for_[long,lat,mon]_array
downsample=function(x){
dx=dim(x)
if (length(dx)==3){
reply=array(NA,c(dx[1]/2,dx[2]/2,dx[3]))
for (i in 1:dx[3]){
reply[,,i]=downsample2d(x[,,i])
}
} else {
reply=downsample2d(x)
}
reply
}
# averages 2.5° gridcells into 5° gridcells for [long, lat] 2D array
downsample2d=function(x){
width=ncol(x)
height=nrow(x)
smallforcing=matrix(NA,height/2,width/2)
for (i in seq(1,height-1,2)){
for (j in seq(1,width-1,2)){
smallforcing[(i+1)/2,(j+1)/2]=mean(c(x[i,j],x[i+1,j],x[i,j+1],x[i+1,j+1]),na.rm=T)
}
}
as.matrix(smallforcing)
}
# EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEE End Functions
# LLLLLLLLLLLLLLLLLLLLLLLLLL LOAD DATA ----- gets the files from the web
# HadCRUT absolute temperature data
absurl="http://www.cru.uea.ac.uk/cru/data/temperature/absolute.nc"
download.file(absurl,"HadCRUT absolute.nc")
absnc=open.ncdf("HadCRUT absolute.nc")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5jan/data.txt","albedojan.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5feb/data.txt","albedofeb.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5mar/data.txt","albedomar.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5apr/data.txt","albedoapr.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5may/data.txt","albedomay.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5jun/data.txt","albedojun.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5jul/data.txt","albedojul.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5aug/data.txt","albedoaug.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5sep/data.txt","albedosep.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5oct/data.txt","albedooct.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5nov/data.txt","albedonov.txt")
download.file("http://badc.nerc.ac.uk/browse/badc/CDs/erbe/erbedata/erbs/mean5dec/data.txt","albedodec.txt")
albnames=c("albedojan.txt","albedofeb.txt","albedomar.txt","albedoapr.txt","albedomay.txt","albedojun.txt","albedojul.txt","albedoaug.txt","albedosep.txt","albedooct.txt","albedonov.txt","albedodec.txt")
# read data into array
forcingblock=array(NA,c(52,144,12))
for (i in 1:12){
erbelist=read.fwf(albnames[i],skip=19,widths=rep(7,13))
erbelist[erbelist==999.99]=NA
erbeout=erbelist[,13][which((erbelist[,1]>-65) & (erbelist[,1]
length(erbeout)
forcingblock[,,i]=matrix(erbeout,52,144,byrow=T)
}
# DOWNSAMPLE FORCING DATA TO MATCH TEMPERATURE DATA,
# and swap lat and long to match HadCRUT data
smallforcing=aperm(downsample(forcingblock),c(2,1,3))
smallforcing[1:72,,]=smallforcing[c(37:72,1:36),,]# adjust start point
# GET ABSOLUTE DATA, TRIM POLAR REGIONS
absblock= get.var.ncdf(absnc,"tem")
smallabs=absblock[,6:31,]
#dim(absblock)
# GET MONTH-TO-MONTH DIFFERENCES
dabs=diffannual(smallabs)
dforcing=diffannual(smallforcing)
dim(dforcing)
#SAVE DATA
save(forcingblock,smallforcing,smallabs,dabs,dforcing,file="erbe_cloud_forcing.tab")
# make cosine weight array
cosarray=array(NA,c(72,26,12))
cosmatrix=matrix(rep(cos(seq(-62.5,62.5,by=5)*2*3.14159/360),72),72,26,byrow=T)
cosmatrix=cosmatrix/mean(cosmatrix[1,])
cosarray[,,1:12]=cosmatrix
cosarray[,,2]
# GET CORRELATION, SLOPE, AND INTERCEPT
#cor(dabs,dforcing,use="pairwise.complete.obs")
module=lm(dforcing~dabs)
m=module$coefficients[2]
b=module$coefficients[1]
#Plot Results
par(mgp=c(2,1,0))
plot(dforcing~dabs,pch=".",main="Cloud Feedback, 65°N to 65°S", col="deepskyblue3",xlab="∆ Temperature (°C)",ylab="∆ Cloud Forcing (W/m2)")
lines(c(m*(-20:15)+b)~c(-20:15),col="blue",lwd=2)
textcolor="lightgoldenrod4"
text(-20,-60,"N = 18,444",adj=c(0,0),col= textcolor)
text(-20,-70,paste("Slope =",round(m,1),"W/m2 per degree C of warming"),adj=c(0,0),col= textcolor)
text(-20,-80,paste("p = ","2E-16"),adj=c(0,0),col= textcolor)
P. Solar says:
October 9, 2011 at 12:03 pm
These are gridcell monthly averages, about twenty thousand of them. Which among them are you proposing to connect with dots, and on what basis? Please think this stuff through before tossing out un-finished ideas.
Please re-read what I have done. ∆T is the change in temperature with time, which is the discrete form of the derivative you are using, dT(t)/dt. No way to put this politely. You’re suggesting that I do what I just did, you’re not paying attention.
Nick Stokes above already pointed out the proper way to deal with the weighting. Read the thread.
Lots of possible ways to do it, but since Nick solved the problem already, not necessary … plus the fact that there’s no way I’d do a smooth as you suggest, that’s just losing information.
w.
P. Solar says:
October 9, 2011 at 12:11 pm
Report back with your results …
w.
P. Solar says:
October 9, 2011 at 1:03 pm
Brilliant, I love unconventional solutions like that. I’ll take a look.
w.
Willis Eschenbach says:
October 9, 2011 at 11:31 am
Steve from Rockwood says:
October 9, 2011 at 10:19 am (Edit)
Willis,
I made an error in my earlier post. Of the radiation reaching the earth, only about 1% is longwave (> 4 um). The other 99% is shortwave. Sorry to display my ignorance but that’s how I learn.
———————-
So virtually all incoming solar radiation is shortwave. That explains Dave’s comments.
Which leads me to a new question: Are clouds better at blocking shortwave radiation (negative) or trapping outgoing longwave radiation (positive)?
Your question is the great unknown in climate science, and is what I’m attempting to answer in this thread.
w.
How much would this change what happens to the clouds?
aaron
“Have you ever thought of doing aerosol masks? ”
OT. willis wants a tight thread, lets give him that seeing as he has devoted a bunch of time to the work
I have updated the head post to include the area-weighting as suggested by Nick Stokes, who has my thanks.
w.
FYI
Longwave infrared (LWIR) begins at 5 micrometers. For all practical purposes sunlight contains no LWIR energy. What little it does have never makes it to the surface as greenhouse gases block it. Also for all practical purposes the earth emits nothing but longwave infrared. There is a wide stark divide in the frequency domains of sunlight and thermal emission from the earth.
Moreover, different materials, not just just greenhouse gases but surface materials such as ice, water, rocks, soil, and vegetation have radicially different absorption and emission characteristics across all frequency ranges. It borders on ludicrous to lump all electromagnetic energy into “Watts” and all matter into “surface” or “atmosphere”. When you make simplifications that border on ludicrous you get answers that border on ludicrous. Garbage in, garbage out. All the fancy statistical math in the world won’t restore highly critical details of how EMR interacts with different forms of matter once they are left out.
Richard G says:
October 9, 2011 at 12:24 pm
“I would like to point you to a musing from the Chiefio. It is about clouds.”
“Gives Us This Day Our Daily Enthalpy”
http://chiefio.wordpress.com/2011/09/27/gives-us-this-day-our-daily-enthalpy/
Excellent. Thank you. E.M. Smith is one sharp cookie…
Peter Miller says:
The piece you are missing is that infinite series sometimes converge to finite values. For example, the geometric series 1 + (1/2) + (1/4) + (1/8) + (1/16) + … converges to the value 2. Hence, if a positive feedback responded to a warming by CO2 of 1 C by producing an additional half C of warming, and then the feedback responded to that half C of warming by producing an additional quarter C of warming, and so on…then what you would end up with is not exponential heating but simply the multiplication of the original heating by a factor of 2.
Steve from Rockwood says:
October 9, 2011 at 9:37 am
Because the sun is its only significant source of energy and for all practical purposes there is zero longwave infrared energy in sunlight. But I see you already discovered for yourself that I was correct.
Now you need to discriminate between the ways by which the ocean absorbs and gives up energy. Radiation from the sun is the only thing that penetrates it. That is it’s only source of energy. Longwave infrared is totally blocked by a surface layer that is only a few microns thick. This skin doesn’t mix downwards. Viscosity overwhelms all other forces in skin layer. The skin layer evaporates. Evaporation is the primary way that the ocean gives up absorbed solar energy – 70% of it. 10% is given up by conduction and the remaining 20% is via radiation. This is RADICALLY different from the way land and ice gain and lose energy. There is no comparison between the two. Greenhouse gases are quite effective over land because their mode of operation is via radiation and radiation IS the primary way that land surfaces give up solar heat. Radiative emission from the ocean on the other hand is not the primary means of giving up energy so greenhouse gases play commensurately smaller role there.
Eschenbach refuses to acknowledge the difference between land and water and how they are quite different in the way they absorb and release solar energy. He’s decided that bringing up actual, well established physics that are counter to his anedotally acquired beliefs is “thread hijacking” and has promised to deprive you of this information so that he may pursue his ridiculous wild geese flying through the clouds.
It’s really kind of amusing. An ostrich with his head stuck in the sand chasing after wild geese…
I’ll clip this comment and post it to tallbloke’s website. Tallbloke has forgotten more science than Willis ever knew.
http://tallbloke.wordpress.com/2011/08/25/konrad-empirical-test-of-ocean-cooling-and-back-radiation-theory/
@Myrrh
> How much would this change what happens to the clouds?
> http://journals.ametsoc.org/doi/abs/10.1175/2010JAS3282.1
It doesn’t change anything. Total solar insolation is 1200 Watts at top of atmosphere. Energy in longwave infrared range (>5um) is 1% of that or 12 Watts. Virtually none of that 12 Watts makes it to the cloud deck as it is absorbed by O3, CO2, CH4, H2O, sulfate aerosols, nitrogen compounds, and a host of other molecules between the cloud deck and TOA.
Dave, how is it possible that shortwave radiation is the only source of ocean heating?
Dave, how is it possible that longwave radiation from the sun has no effect on ocean temperature?
Thanks.
[UPDATE – See, this is the problem. Dave starts off with his claims, and right away someone (reasonably) says “Hey, that makes no sense at all”, and then we’re off the rails. I’ll leave this here to illustrate the point, but I will snip further responses along the same line. If you want to dispute the DLR/ocean issue with Dave, please do so elsewhere. w.]
UPDATE TO UPDATE No Willis, if you had the first clue about the physics you’d know that the response was not reasonable. There is practically no longwave infrared energy in sunlight. The responder apologized and corrected himself later in this thread for his error. ow if only Willis had the means, motive, and self-esteem to correct himself we could make some progress and people could actually learn some science from one of Willis’ misadventures into statistical analysis which, by all indications, is another anecdotal skill he picked up via osmosis while either surfing, sailing, fishing, trying to get a dishonarble mental discharge from the Army, or whatever his source of knowledge which decidedly does not include going to school or studying any science textbooks.
[Reply: You are invited to submit an article of your own if you wish. You might be surprised at the comments. ~dbs, mod.]
Willis,
Good work. Applause, applause.
you wrote: That’s why we invented mathematics, all right. Part of the problem that there are over eighteen thousand data points, so the underlying pattern is obscured. That doesn’t mean it’s not there, just that this particular visualization breaks down under the weight.
Usually, the more information you have on a graph, the better. This might (repeat might) be one of the exceptions: I recommend a separate panel for each region. Also, I recommend a separate analysis for each region, with a test of whether the regression coefficient is the same in all regions. Lastly, I recommend a region BY month breakdown, and a test of whether the slope varies across the year, between regions, of if there is a region BY month interaction. These are not criticisms, these are requests for follow-up analyses.
Mike McMillan wrote this: Perhaps a contour map of point density would help.
Think of our recommendations as complimentary, not conflicting.
It’s true that the R^2 value is small, but small effects over large regions and long times make non-trivial differences. Besides, the sign is opposite to the “conventional wisdom”.
Kevan Hashemi wrote: Clouds reflect sunlight back into space. Whatever effect they have upon long-wave radiation is entirely secondary to the fact that they block the sun’s heat from reaching the surface of the planet. Here’s a simulation that shows the effect:
It is important to have actual data backing up important claims. Willis has actual data at the “explantory level” that counts: actual clouds and actual temperatures.
Dave Springer wrote: R^2 = 0.11 means there’s no significant pattern.
And yet, the AGW promoter Andy Dessler got an even smaller R^2, probably with the wrong sign, published in Science.
Dave Springer wrote: This generalization does not hunt. The type of energy and type of surface makes a difference. For instance if the surface is ice and the energy is shortwave it will raise the temperature very little. If the surface is water and the energy is shortwave it will raise the temperature very much.
…
More to the point, if the energy is longwave and the surface is water it will not raise the temperature very much. This is predicted by the physics and experimentally confirmed.
That’s pretty incidental to the empirical relationship that Willis has reported. Exactly how much of this is longwave, and exactly how much is shortwave, and exactly how much is on land and exactly how much is at sea, and exactly how much is am and how much pm (see Willis’ work on the TAO/TRITON data) can be subjects for subsequent research and analysis.
Willis, I see that you got requests for more work, not just mine. Have fun.
And keep up the good work.
Willis says: “These are gridcell monthly averages, about twenty thousand of them. Which among them are you proposing to connect with dots, and on what basis? Please think this stuff through before tossing out un-finished ideas.”
First of all Willis, don’t get pissy. I am trying to make useful and constructive criticism. I’ve already done this kind of analysis and more and can probably point out some errors and make useful suggestions. If anyone is ” tossing out un-finished ideas” it’s you , as you readily admit in your post.
You specifically asked for people to check over what you’d done so don’t get offended when they do.
Just looking at the pretty colors 🙂 it would appear that some of the commenters are correct: you need to divide it up between landmass and ocean.
As a preliminary proxy for that, do latitude bands. 0 – 65 at 5 degree cells gives 13 cell bands. Do it 5/4/4, that is, 0-25, 25-45, and 45-65. I predict that the 0-25 band will look very much the same whether Northern or Southern hemisphere, but the others will be quite different because the land distribution is quite different. You really want to separate out land and ocean, but bands might be simpler as a first cut.
Regards,
Ric
Willis says: These are gridcell monthly averages, about twenty thousand of them. Which among them are you proposing to connect with dots, and on what basis? Please think this stuff through before tossing out un-finished ideas.
Why don’t you think it through before dismissing it? Part of your problem is dumping that much incompatible data it one bucket. Read my posts , you may learn something. Don’t let it go to your head that Anthony let you post on his blog.
Please re-read what I have done. ∆T is the change in temperature with time, which is the discrete form of the derivative you are using, dT(t)/dt. No way to put this politely. You’re suggesting that I do what I just did, you’re not paying attention.
Yeah, but you are plotting against dR/dt , so once again you have not even though about what I said before rejecting it.
I don’t see your objection to the cosine weighting . ERBS data is W/m2 so you must relate it to an area related term.
Nick Stokes above already pointed out the proper way to deal with the weighting. Read the thread.
This thread was already >110 posts by the time I got to read it. Most of that is the usual garbage , bitching and inane comments. Life’s too short. I did scan the whole thread for anything you’d posted and picked up on help Nick had provided you with.
Your post said you were still unhappy about the idea of weighting, I was saying it was correct to do it and why. Again I see no reason for your off-hand response.
As a means of reducing the 2.5 gridded data I’d suggest a 1:2:1 binomially weighted average in lat , then the same in longitude. Alternatively interpolate the 5 degree data with a Catmull-romm spline (quite easy actually)
Lots of possible ways to do it, but since Nick solved the problem already, not necessary … plus the fact that there’s no way I’d do a smooth as you suggest, that’s just losing information.
If you want to move from a 2.5 grid to a 5 deg grid, how do you imagine you can do that without “just losing information”? The fact that you later throw away the whole time dependency of BOTH datasets makes you comment about losing information rather comical.
Don’t get your back and get narky, just read the posts, you may learn something. If you are unable to learn anything I’m sorry, I tried let you benefit from some of the things I’d learnt in doing a similar analysis, I don’t intend to spend my days arguing about it.
best.
PS, one way to test what I say about regression attenuation is to invert the axes and redo the ols fit. If the result is good one slope will be the inverse of the other. I think you will get a shock how different they are. Clearly the result should not depend on which way you plot the data … unless I’m right about OLS.
Willis says: “Report back with your results …”
My result IS just this, that doing OLS on this kind of data is fundamentally wrong and I’ve said why. I don’t need to reproduce your errors to prove it to you. Do it yourself. [hint reverse the arguments you pass to lm() ]
I did post a link to a graph where I have done this on similar data and and link to a thread where I covered this with TroyCA. He did listen and agreed.
If you want to dismiss the “positive feedback” results of Dessler, F&G et al this is the easiest and most incontrovertible route. Since the results of this flawed method are the only empirical justification for positive feedback in models and that is the only reason the models produce alarming warming this seems to me to be the crux of the whole AWG battle.
To spell this out: Dessler (and others) do a similar scatter plot method and get an incorrect slope for similar reasons. It is a mathematical result that the ols estimator WILL be less than the true “slope”. How much less depends on the amount of noise. The data has more noise than signal. The regression attenuation will be very significant.
Since the ols value is less than the neutral 3.3 W/m2/K of the plank response he comes to the erroneous conclusion that he has demonstrated positive feedback.
“”””” Steve from Rockwood says:
October 9, 2011 at 10:19 am (Edit)
Willis,
I made an error in my earlier post. Of the radiation reaching the earth, only about 1% is longwave (> 4 um). The other 99% is shortwave. Sorry to display my ignorance but that’s how I learn.
———————-
So virtually all incoming solar radiation is shortwave. That explains Dave’s comments. “””””
Well take heart Steve; ignorance is NOT a disease; we are all born with it. The Library of Congress is one of the world’s largest repositories of information about most of which, most of us are totally ignorant.
But I’m happy to see that you found some correction for your earlier statement that 53% of the solar spectrum was “longwave”.
It would be useful, if persons who are interested in this whole climate energy business acquainted themselves with certain rather basic pieces of information. The first I would suggest is the specific shape of the Planck Black Body radiation curve.
The version I would most strongly recommend is the universal normalized BB curve as published in the standard Optics Text; Modern Optical Engineering, by the late Warren J Smith, who worked for many years at Infra_Red Industries in Santa Barbara California.
He plots the Spectral Radiant Emittance (W.cm^-2.micron^-1) as a fraction of the peak value at the spectrum peak, as a function of the wavelength also scaled to the spectral peak wavelength.
So the relative wavelength goes from 0.15 to 50, and the vertical axis has two scales ; a linear one from 0.0 to 1.0, and a logarithmic one from 1.0 down to 10^-5. From this one can deduce the behavior of almost any thermal radiator, at least qualitatively and often quite quantitatively.
Some important general properties (of BB radiation) are.
1/ (almost) exactly 25% of the total energy is radiated at wavelengths shorter than the peak wavelength, and 75% at longer wavelengths.
2/ Only 1% of the total energy survives at wavelengths shorter than 0.50 times the peak wavelength.
3/ 50% of the total energy is emitted below about 1.5 times the peak wavelength (closer to 1.45).
4/ Only 1% of the total energy is emitted at wavelengths longer than 8 times the peak wavelength.
5/ The Spectral Radiant Emittance is down to 10^-5 at 0.20 times the peak wavelength and 40 times the peak wavelength.
The other most important relationship is the Wien Displacement Law, which says that the Temperature (Kelvins) times the peak wavelength, is a constant (2897.8 Kelvin.microns.
Apologies to the chemists and molecular spectroscopists who think in terms of Wave numbers (frequency) rather than wavelength. Despite the rationale for using frequency; given that
E = h.(nu) , I have yet to see my first plot of the solar spectrum plotted on a frequency scale.
Translated to the solar spectrum TSI, in its near BB guise at roughly 6,000 K. Other sources give 5900 K as the best fit BB curve fit, or 5770 as the value best meeting the TSI value which was 1353 Wm^-2, when that was published, so likely 5800 K at today’s 1366 Wm^-2.
So Wien’s law gives about 0.5 microns for the TSI spectral Peak, so 98% of the energy lies between 0.25 microns, and 4.0 microns. and the 50% division is at about 740 nm.
As for “visibility”, The commission on Colorimetry of the Optical Society of America cites Standard Luminosity values for the Photopic Luminosity curve of the eye from about 385 nm to 760 nm. at which points the value is 10^-4. This is not the extreme of the eye sensitive range. The dark adapted eye is capable of seeing single photons, at least in the blue-green range, so 10^-4 is not invisible. Warren Smith gives eye sensitivity values in the near infra-red out to 1.1 microns wavelength, where the sensitivity is 10^-11..
I tend to place the solar energy “short wave” radiation as being the 200 nm to about 2.0 microns range even though about 7% of solar energy survives past that. But the “long wave” part of the infra-red spectrum has some effect from about 2.0 and certainly by 4.0 microns.
A final key plot to have is the air mass one ground level solar spectrum. This plot corresponds roughly to the clear air surface maximum irradiance; before any albedo effects are considered.
Well not all, because the first major return of solar spectrum energy to space, in elastic scattering processes, is the Raleigh scattering of the short wavelengths which result in the blue sky. The blue sky looks the same looking down on the blue planet, as it does looking up, so the blue sky you see looking ap remote from the solar disk, represents exactly the amount of solar energy in the shorter wavelengths that is immediately lost to space, and is a large part of the loss from the 1366 Wm^-2.
The big mystery to me, is this almost childish fascination with the major climate calamity that is purportedly created from a source that is modelld quite closely by the LWIR emission from an ordinary bottle of drinking water at about room Temperature. It emits a near BB spectrum at about 400 Wm^-2 with a spectral peak of about 10.1 microns; corresponding to the average surface Temperature of the planet, as immortalized By Dr Trenberth’s amazing cartoon drawing.
Meanwhile, there is virtually total neglect of the effect of water both as a vapor or as clouds, on the real energy source driving our climate; the energy supplied by the sun.
FEEDBACK networks normally operate, by sensing the system OUTPUT, and modifying the system INPUT. One does not normally connect a feedback network to some largely irrelevant noisy output, which is what the surface LWIR actually is.
On any average mid-summer day in the north, the earth can enjoy a total surface Temperature range from as low as -90 C at places around Vostock in Antarctica, up to about +60 deg C in North African and middle eastern deserts. The surface LWIR emittance covers at least an 11 to 1 range, from nearly double the sacred average of 390 W/m^2, to about 1/6th of that value.
And that entire range can exist simultaneously; so we worry about whether or if, the average may have changed by 0.6 deg C in the last 150 years, and that from a quite corrupted and improper global sampling methodology.
When you connect the major climate feedback (water) to the INPUT (the sun) you find, that the feedback can never ever be positive. (more) Water can only lower the solar energy that reaches the earth surface; it can never increase it. CO2 also acts to reduce incoming solar energy, but not to the extent that water does.
Interesting and informative post George. However you seem to have misunderstood what ‘positive’ is all about.
“Water can only lower the solar energy that reaches the earth surface; it can never increase it.”
I’m not aware of anyone saying that is not the case. That part is negative. However, when you look at the similar effect of outgoing IR a similar argument would lead to saying this can only be positive. The big question is the overall, net effect on ERB.
Unfortunately, for the reasons I’ve stated above, this cannot be resolved by doing invalid OLS regression on scatter plots.
That Willis makes this error in a blog post is excusable. That the likes of Dr. Dessler do this in the published literature is rather less admirable.
Troy’s post I linked, looking at water vapour vs CERES radiation budget had an r2 of 0.03
He initially posted a “slope” of -0.19 , when he did it the other way around he got -6.4
We’re not talking +/-10% here.
Least squares does not work unless your x variable is a CONTROLLED variable.
In Willis’ case here I expect it will give a much more negative slope the other way around. So it may establish the sign of whatever it is measuring.
In the Dessler / F&G analysis the fundamental question is which side of 3.3 is it. This means the attenuation can change the “net feedback” from negative to positive.
Figure 1 states: “Data colors indicate the location of the gridcell, with the Northern hemisphere starting with red at the far north, slowly changing to yellow through the tropics, and to red at the far southern latitudes.”
Did you mean to mention blue dots in there somewhere?
Jim
Willis Eschenbach wrote: Report back with your results …
Too often you get angry instead of reading. You have a nice analysis. Why ruin it by ignoring its shortcomings? P. Solar is correct: you can not estimate the slope that you want to estimate. It’s well known.
Willis wrote: These are gridcell monthly averages, about twenty thousand of them. Which among them are you proposing to connect with dots, and on what basis? Please think this stuff through before tossing out un-finished ideas.
OK, so you don’t have time. Let it go at that instead of the pseudo-intellectual dismissal.
I may be completely misunderstanding what you’re doing here, but if you are subtracting the temp of the month ahead from the current month I can see why Nick Stokes implies there could be a seasonal bias – and one which would presumably apply in one direction for one half of the year in each hemisphere (days lengthening or shortening) and then in the reverse direction for the rest of the year (for that hemisphere). The hemispheres would be out of phase and you’d probably expect a 4 quadrant plot which is what you show. If you treated the hemispheres separately AND split the year into 2 parts for each hemisphere you could run the subtraction both forward and backward for each of the 4 hemisphere/season combinations. IF each forward/backward run showed the same trends you could then rule out the seasonal signal (which may be much of your ‘noise’). If there’s still a negative trend in temp with increasing cloud, you’ve got a very interesting result.