New Computer Fund

Sunday, May 5, 2013

Land Temperature and Ocean Heat Content and Baseline Selection


This is one of those in progress posts.  Since there was such a good correlation between North Atlantic Sea Surface temperatures and "Global" land surface temperatures, I was curious how well the Ocean Heat Content (OHC) data agreed with the land surface temperature data.  The data is the Climate Research Unit CRUtemp4 data (global) and the 0-700 meter OHC data from NOAA.

The data is massaged a bit with a common baseline anomaly from 1955 to 2012 to match the OHC data period, then normalized by dividing by the standard deviation of the baseline period. As you can see there is a very good fit between the CRUtemp4 data and the North Atlantic OHC data.  Near the end, about 2005, there appears to be a divergence between CRU4 and Global OHC, but that is a pretty short period and it includes the ARGO data.  Time will tell if that is real or not.

Pretty spiffy eh?  Well, since I used the full length of the OHC data for the anomaly period, that pretty much has to fit.  I have beat the data into submission.

Here I have added the CMAR Sea Level Rise data using the same baseline, almost, 1955 to 2006 so you can see how well I can whip data.  If you extent the trend lines for each of the series the OHC and GMSL will have about the same trend while the CRU4 series will have a much different trend.  That appears to make Global Mean Land Surface Temperature the odd data out using this anomaly baseline. 

But is it Land Surface Temperature or just the Hemispheres since most of the land surface is in the Northern Hemisphere and has a lot of influence from the North Atlantic and the Gulf Stream?



Here is the Hadley Center SST version 3 with the GMSL with the same massaging using 1880 to 2006 for the baseline period.  Using the GMSL as a reference, in this case it appears that the Northern Hemisphere has much larger variations around the mean sea level "mean". 




So let's replace the NH SST with the "Global Land Surface Temperature".  This is using the same 1880 to 2006 baseline for anomaly with the same normalization.  Temperature of the Southern Oceans, Sea Level Rise and Global land surface temperature seem to make some sense with the MSL "mean".   Base on this plot, temperatures and sea level have been rising together since at least 1880 with some "wandering" of both the SST and GLST.  I could select a different baseline period and get radically different results.

Finding the right "reference" to avoid chasing wild geese should be priority one in paleo climate reconstructions. 


Friday, May 3, 2013

The Golden Ratio Control Range

A. M. Selvam, is an Indian professor that is into chaos theory and all the associated things that go along with it.  Most of her papers include the Penrose pattern which is an example of fractals based on the Golden Ratio which is 1-SQRT(5) divided by two with is an infinitely repeating number.

While the Golden ratio or Phi is neat, I found that the ratio part 0.6180_  appears to be the heart of the situation.  The range 0.6160_ with its counter 1-0.6180_=0.3819_ makes for a very stable control range.

Since entropy in a system that exists over eons of time must not continuously gain/lose significant energy or mass to exist in a stable state must be equal to 50%, i.e. energy in must equal energy out and mass in must equal mass out, the range 0.6180_ to 0.50 to 0.3819_ is a common natural control range.  Just consider the difference, 0.6180_ - 0.5 = 0.1180_ and 0.5 - 0.3819_ = 0.1180_

This is likely not a "law" of nature, more like a convenience of nature.  If you stray too far from "normal" or Ein=Eout, you cease to exist.  The things we see in nature are still here, the Golden Ratio appears to be a reasonable limit for a "stable" system whether it is biological, mechanical, thermal or chemical.

The only reason I am posting this is because the range of natural variability is a bit of a question mark.  When I have mentioned that there are some parts of chaotic systems that are simple, eyes tend to glaze over because I am not a certified expert on chaotic systems.  I have worked on chaotic systems and am capable of making observations, but they are chaotic systems so you should take any observation with a grain of salt.  Recurrent patterns though are typical of chaotic systems and the Golden Ratio Control Range appears to be a convenient reference for many chaotic systems.  

So this post is just to stimulate some that wonder about order out of chaos to consider simple recurrent patterns in nature as references.  Well, also because of the conservation of water/energy part of the ocean puzzle.

In order for water to exist on Earth, Ein the ocean has to equal Eout of the oceans and evaporation has to equal condensation, over some reasonable time frame.  Since the clouds in the atmosphere are there because of ocean evaporation, there should be a "simple" mechanism for the clouds to regulate the energy.

With better satellite and surface data, it appears that the latent heat of evaporation from the oceans happens to equal the solar energy reflected by the clouds. 

The data in this table if from the Trenberth, Fausto and Keihl 2009 energy budget which I have criticied on several occasions.  Even in their study, solar reflected energy closely match the ocean latent heat loss.

Depending on what values you use for averages, total solar energy reflected with 1361Wm-2 TSI at the top of the atmosphere with 30% albedo is 408 Wm-2.  That is roughly the equivalent energy of the average surface of the Earth.  If you use TSI/4 as your "Average" TSI, then your "average" energy reflected would be 108 Wm-2.  As a general rule though, reflected energy is equal to the stored energy.  Adding CO2 to the atmosphere would increase the potential to store energy and the system should respond by reflecting more energy.  Since clouds are the main source of the albedo, increased temperature should increase cloud cover which will offset some portion of the additional stored energy.   

That is a fairly logical train of thought that required a bit of creativity to step around in order to get the higher end of climate sensitivity estimates.  If cloud albedo is directly related to conservation of energy and water in the oceans, more than just a feed back, an actual control setting, that would be interesting.  It appears to be even more interesting than just clouds. 

A while back I was comparing the TSI and Albedo of various planets in our solar system.  NASA has Planetary Fact Sheets with the information I was using.  While setting up a spread sheet to compare the inner planet I made a little mistake that was comical enough to do a post on.  My mistake indicated that the surface temperature of a planet can be estimate by just the TSI times the Bond albedo.  When both the Bond and Geometric albedo are considered, even a more interesting estimate is possible. 

 

Mercury Venus Earth Mars
Bond 0.07 0.90 0.31 0.25
Geometric 0.14 0.67 0.37 0.17
Ratio 2.09 0.74 1.20 0.68
TSI 9126.60 2613.90 1367.60 589.20
Bond Flux 620.61 2352.51 418.49 147.30
Geo Flux 1295.98 1751.31 501.91 100.16
Diff 675.37 -601.20 83.42 -47.14
Sub surface -54.76 2953.71 335.06 194.44
Sub K 0.00 477.74 277.26 241.99
Bond K 323.45 451.32 293.11 225.76
Geo K 388.82 419.22 306.73 205.01


The Bond and Geometric Albedo values and Total Solar Irradiance (TSI) are from the NASA Planetary fact sheets in the earlier link. The Ratio row is just Geometric divided by Bond albedo values.  The Golden Ratio range is 1.26 to 0.76.  Venus and Earth fit in that range and both have relatively stable atmospheres.  Mars has a thin atmosphere and Mercury virtually no atmosphere. 

The Flux values are just the respective albedo times the TSI.  The Diff or difference is provided and below that the estimated surface temperatures are calculated using the Stephan-Boltzmann relationship.  Sub surface is the difference between the Bond and the Diff value. 

Since the Geometric flux is based on the maximum albedo value, that should be roughly the peak equatorial flux.  For Earth the values 501.91, 418.49 and 83.42 are fairly close to the current equatorial, mean ocean and meridional flux values.  The sub-surface value of 335.06 is also very close to the "average" energy of the oceans on Earth.  For Venus, all the values are lower than the "true" surface values on Venus, but there is some indication that geothermal energy is a large portion of the energy near Venus' hard surface.  Mars' sub surface value based on this theory indicate there may be some geothermal energy impact, though not much.  Mercury has a very large discrepancy that may be due to a number of factors.  Then there is of course estimate errors and my theory may be crap.   There is some circular reasoning involved.

The surface temperature dependence on albedo is nothing new.  Surface albedo responding to surface temperature is nothing new.  The combination with Bond and Geometric albedo though appears more than just a numeric oddity.   I can't help but think something has been missed with all the climate Paradoxes that seem to crop up so frequently.

This is a more interesting puzzle it seems every day.  I am going to leave this as food for thought since the data for other planets is so coarse and focus on the one we have the most data, Earth.


 
This chart shows the approximate Sea Surface Energy (SSE) based on the ERSST data set using the Reynolds Oiv2 satellite data from 1981 to 2013 for "calibration".  30 degree latitude zones 60S to 60N were used then averaged as 323.6 Average which is assuming a uniform ocean area distribution or simple average and Areal which considers the actual surface area.  The simple average has a mean of ~415 Wm-2 and the Areal Average has a mean of ~420.5 Wm-2.  Both compare well with the 418.5 Wm-2 based on reflectance using Bond Albedo. 

Using the 30S-equator and 30N to equator bands, the average SSE is 449 Wm-2 and 457Wm-2 respectively, lower than the 501.9 estimated with the geometric albedo, but the bands are some what wider than what would be considered direct peak insolation.  The meridional difference in the SH is 80 Wm-2 and in the NH 73 Wm-2, some what lower than the 83 Wm-2 difference estimated with the geometric minus Bond flux but reasonable given the difference between Areal and simple averaging. 

The reason for using the satellite SST for this comparison is because the actual SST and SSE is very uncertain.  According to Stephens et al., the true surface energy of the entire Earth is close to 400 Wm-2 with a +/- 17 Wm-2 margin of error.  Not having a better estimate of the absolute surface temperature appears to matter much more that most will admit. 

Using more current and likely more realistic data, the albedo reflect energy appears to provide a reasonable gut check of the overall system state.

As a reminder, to exist, energy and matter have to be conserved or Ein must equal Eout over some reasonable time period.  That requires 50% entropy, energy lost must equal work produced in an open system that is to remain stable.  The Golden Ratio appears to provide a reasonable estimate of the range of possible imbalance allowed during that reasonable time period for an open system to return to a stable condition.  If albedo can adjust to maintain that range it appears to be quite adept at doing it. 

This is an observation not an attempt at a proof, but interesting none the less. 

Thursday, May 2, 2013

Monty Hall or Not.

The Monty Hall is a statistical puzzle that I try to use to explain how decision making can be skewed toward simpler logic or off into a tizzy.

Monty Hall's game show had three doors.  One with a prize and two with trash or goats.  With three doors you have a one in three chance of picking the right door.  After picking a door, Monte would give clues or offer bribes to get you to change.  Without anymore information, each door still had a one in three chance of holding the prize.  If Monte open one of the doors that did not hold the prize, then the odds changed.  Your door was still one in three, but since the other two had a 2 in three chance, with one opened that remaing door still had a 2 in three chance of holding the prize.  You would always be better off switching doors in that case.

The confusion is that the initial odds can't be forgotten.  Even though there are two doors and one prize, you only had a one in three chance to start, you will do better switching because Monte wouldn't open the prize door until the end of the show.  That biases the odds.

I mention this again because the discussion invariably turns to the economics of climate mitigation.  The original odds were 4C proposed by James Hansen  and 2 C proposed by S. Manabe.  Because of a disagreement between those two scientists, Jules Charney proposed a compromise and used the average of the two estimates.  The estimate for sensitivity became 3C +/- 1.5C. 

Now if you are an economist trying to figure the cost and benefits of CO2 policy, you have choices.  Since the 3C is not a true estimate, you can ignore that value and use 2 C or 4C.  With 2C considered "safe" or at least managable, you have an act/don't act decision.  That makes wait a while and see a valid option for policy action and no regrets policy action with cost effective options with incremental steps as new information and technology becomes available the best option.

Using 3C, the odds are stacked toward imprudent policy action choices with less flexibility.  The game is biased not by science, but by a political compromise.  A goat door was added that was not part of the initial odds.  With three options and two deemed "bad", the activists win.  You lose the game. 

You can't even point out this sleight of hand, whether it was intentional or not, because it is supposedly based on science and doubting science makes you some kind of Neanderthal idiot.  This is the legacy of the Merchants of Doubt.  Flawed but scientific sounding statistics are used to separate you from your wallet.

Just a thought for a rainy day.