Category Archives: Climate

Climate studies.

Global Temperature Reanalysis Baseline Comparisons

I have been downloading the daily average global surface air temperature data initial condition output for the Sigma 0.995 level from the National Centers for Environmental Prediction and National Center for Atmospheric Research Reanalysis 1 (NCAR R1) cooperative effort.  This output is still being updated daily about 2 days behind the current day and extends back to 1948.  So far I have downloaded and processed the daily temperature output back to 1979.  The Sigma 0.995 level corresponds to the pressure altitude at 99.5 % of the surface air pressure, which is roughly about 50 meters above ground level.  The actual height above ground level varies somewhat depending on atmospheric conditions.  The NCAR R1 model uses a 2.5 degree latitude by 2.5 degree longitude grid.  In this post I am comparing the resulting global mean surface air temperature anomalies (GMSATA) for several time periods and different reference baseline periods.  See the Methods section at the end for calculation details and links to sources.

The first two graphs, Figures 1 and 2, display NCAR R1 daily average GMSATA time series output for two different reference baselines.  I was expecting the older 1979-2000 baseline to show higher temperature anomalies than the much more recent 2011-2015 baseline as can be seen in the graphs, but I was surprised that the baseline shift also changed the seasonal pattern across each year.  Apparently the seasonal patterns have shifted from one baseline period to another.  Figure 3 is provided to compare the daily average 2-meter above ground level surface air temperature anomalies from the more modern Climate Forecast System Reanalysis (CFSR)  to the NCAR R1 output.  The CFSR has a higher resolution of 0.5 degree latitude by 0.5 degree longitude.  See the Daily Updates page, which can be accessed from the menu bar at the top of this page, for more information about the daily CFSR output.

Figure 1. NCAR R1 GMSATA daily 2018 January through 2019 February 13, referenced to 2011-2015 baseline

Figure 2. NCAR R1 GMSATA daily 2018 January through 2019 February 13, referenced to 1979-2000 baseline

Figure 3. CFSR GMSATA daily 2018 January through 2019 February 14, referenced to 1979-2000 baseline

I also calculated GMSATA for two other reference baseline periods, 1994-2013 and 1981-2010.  The time series results for all four reference baseline periods for 2018-2019 are plotted together in Figure 4.  They all converge around January-February and diverge the most around September-October.  I’m not sure why.

I included the 1994-2013 reference period used by Nick Stokes for reporting daily averages here and the 1979-2000 period used by the University of Maine Climate Reanalyzer here.  The 1981-2010 period is the most recent three decade (30-year) period commonly used for climatological data reporting.  I previously used the 2011-2015 period for comparing monthly NCAR R1 versus CFSR temperature anomalies here.

Figure 4. NCAR R1 GMSATA 2018 through 2019 February 13 reference baseline comparisons

The next three graphs, Figures 5 through 7, are like the first three graphs, but covering a longer time period, from 2014 to 2019 so far.  Again the general patterns are similar, but the details differ.

Figure 5. NCAR R1 GMSATA daily 2014 through 2019 February 13, referenced to 2011-2015 baseline

Figure 6. NCAR R1 GMSATA daily 2014 through 2019 February 13, referenced to 1979-2000 baseline

Figure 7. CFSR GMSATA daily 2014 through 2019 February 14, referenced to 1979-2000 baseline

The last two graphs in Figures 8 and 9 cover a longer time period, for the current century so far, beginning in 2001.  I do not yet have CFSR daily averages for all of this period, so only the NCAR R1 results for two different reference periods are presented.

Figure 8. NCAR R1 GMSATA daily 2001 through 2019 February 13, referenced to 2011-2015 baseline

Figure 9. NCAR R1 GMSATA daily 2001 through 2019 February 13, referenced to 1979-2000 baseline

Overall, this exercise goes to show that changing reference baseline periods for daily GMSATA does cause quite a bit of variation in the results – more than I expected.  However, the general trends as indicated by the running 365-day averages did not appear to be affected by changing baselines.


For the NCAR R1 daily averages I downloaded the gridded Sigma 0.995 level temperature output which is provided in compacted annual files (thanks to Nick Stokes for providing the link below on his blog).  I used the National Aeronautic and Space Administration (NASA) Panopoly program to extract the temperature grids from the compacted data files and then loaded the daily temperature grid data into spreadsheets by year.  For each day I calculated daily averages by latitude zone, weighted by area, to compile zonal and global averages.

I calculated the reference period baseline averages for each day and used centered running 5-day averages to smooth the baseline results.  Once the reference baseline temperatures were calculated, the temperature anomalies were calculated for each day by subtracting the reference baseline value for the day from the actual daily temperature average for that day.  For the temperature anomalies I also calculated running 365-day and 91-day averages to show annual and seasonal scale tendencies.

NCAR R1 annual compacted files with daily global temperature grids:

NASA Panoply program:




GOES-17 Is Now Operational

“In its new role, GOES-17 will serve as NOAA’s primary geostationary satellite for detecting and monitoring Pacific storm systems, fog, wildfires, and other weather phenomena that affect the western United States, Alaska, and Hawaii.”

“The latest milestone for GOES-17 comes exactly eleven months after the satellite first reached its geostationary orbit 22,000 above Earth. Launched March 1, 2018, GOES-17 is NOAA’s second advanced geostationary weather satellite and the sister satellite to GOES-16 (also known as GOES East). Together the two satellites provide high-resolution visible and infrared imagery as well as lightning observations of more than half the globe – from the west coast of Africa to New Zealand, and from near the Arctic Circle to the Antarctic Circle.”

Large winter storm system moving into California 2019 February 2:

High resolution animation of clouds around Hawaii 2019 January 15:

High resolution animation of brown smoke plume blowing toward the west from the large from the Woolsey Fire west of Los Angeles  in California, as high level cirrus clouds above blow to the east on 2018 November 13:

Read more here.

Real-time images here.

The climate is changing!

The climate is changing!
The climate is changing!
The climate is changing!

And its all YOUR FAULT!!!

Say the Chicken Littles of today.

Mostly for pride, prestige, fame, and funding.
And to promote potentially catastrophic policies.
All in the false pretense of “green” and “renewable” energy.
Endeavors that will enrich the Chicken Littles promoting them.
Ironically these policies are not likely to stop climate change.
But will instead impoverish and endanger most of us.
Wasting huge amounts of hard earned tax payer dollars.

And to top it all off …
The hidden Chicken Little motto:
Do as I say, not as I do.
Hypocrisy at its best.

Chicken Littles of the world are in a class of their own.
Not subject to the rules they espouse for the rest of us.
They say we must quit using fossil fuels.
But that is not what they do:
Flying in private jets. Living in huge walled mansions.
Driving expensive cars. Living lavish lifestyles.
All provided by fossil fuels.

Climate change is now blamed for almost everything bad.
Never mind digging for real causes that are ignored.

Consequently, “climate change” has become a religion.
We must have faith in the sacred climate models.
Never mind that they have not been validated.
The unfaithful are persecuted as “deniers”.

When will the gullible public learn?
When lies are told often enough.
They are perceived as truth.
Buyer beware.

How the climate models fare

The graph below provided by Clive Best compares global surface air temperature anomaly projections from a large number of climate model runs for different Representative Concentration Pathways (RCPs) for carbon dioxide versus two estimates of global surface air temperature anomalies based on actual temperature measurements from Cowtan and Way and the UK Hadley Climate Research Unit Temperatures (HadCRUT4.6).  Notice that the measurement based estimates are indicating global temperatures at the low end of the model projections.

Not at all alarming.

Actual climate model temperature projections

Another graph from Clive Best below shows direct climate model output of projected global average surface air temperatures, indicating a very wide range in temperatures that is hidden by using temperature anomalies in the graph above.  Another sign that the climate models are not representing the real world very well.

Graphs sourced here, thanks to Clive Best.

And if you think climate is changing now, take a look at how much climate has changed in the last 3 million years without any help from humans:
Paleo Climate

Earth Orbital Eccentricity Effect on Global Mean Surface Temperature

I thought it would be an interesting exercise to look at the annual change in global mean surface temperature (GMST) in the atmosphere at 2 meters above ground level versus the annual change in solar radiation incoming (SRI) at the top of the atmosphere (TOA).

The earth’s orbit is slightly elliptical.  Eccentricity is a measure of the departure of an ellipse from a true circle.  The earth’s eccentricity varies over long time scales and is estimated to have ranged from a low of 0.000055 to a high of 0.0679, with a geometric mean of 0.0019.  The present eccentricity is 0.017 and decreasing very slowly.  That doesn’t seem like much, but the earth’s distance from the sun currently increases by 3.4% over the course of a year from the minimum distance in early January to the maximum distance in early July.  That variation in distance causes about a 6.8% increase in TOA SRI relative to the minimum amount at the greatest distance in early July in order to reach the maximum amount at closest distance in early January, assuming a constant output of radiation from the sun.

The US National Aeronautical and Space Administration (NASA) Clouds and the Earth’s Radiant Energy System (CERES) measurements from satellites include TOA SRI.  I downloaded monthly CERES TOA SRI for the period from March 2000 through July 2018.  I also have compiled monthly estimates of GMST from the Climate Forecast System Reanalysis (CFSR) output of initial conditions four times per day, with monthly data prior to 2010 provided by the University of Maine Climate Change Institute (UM CCI) and data from 2010 through 2018 calculated from the Climate Data Assimilation System (CDAS) output of 2-meter surface air temperature.

A time series graph comparing TOA SRI with GMST is provided in Figure 1 below (click to enlarge).  It shows the annual cycle of GMST with a peak in July and minimum in January, whereas the TOA SRI peaks in July with a minimum in January – the complete opposite.  At first, this seems counter-intuitive since increased incoming solar radiation should cause increased temperature and yet the data show the opposite.  My guess is that there is a time lag mediated primarily by the oceans and possibly also affected by the much greater percentage of land in the northern hemisphere, but I have not read up on the subject.  This graph also includes centered running 12-month averages (Run 12) for both TOA SRI and GMST.

Figure 1. Comparison of CERES TOA SRI with CFSR GMST
(click to enlarge)

Figure 2 below compares the CERES TOA SRI with the CFSR Northern Hemisphere (NH) mean surface temperature.  The phase of the annual cycle in temperature is similar to that of the global temperature cycle and is nearly opposite the phase of the TOA SRI cycle.  Notice that the amplitude of the annual temperature cycle is much larger than the global cycle.

Figure 2. Comparison of CERES TOA SRI with CFSR Northern Hemisphere Mean Surface Temperature (click to enlarge)

Figure 3 below compares the CERES TOA SRI with the CFSR Southern Hemisphere (SH) mean surface temperature.  The phase of the annual cycle in temperature is very close to the phase of the TOA SRI cycle, but advanced forward by about a month.  Notice that the annual temperature cycle is larger than the global cycle, but not near as large as the NH cycle.

Figure 3. Comparison of CERES TOA SRI with CFSR Southern Hemisphere Mean Surface Temperature (click to enlarge)

Figure 4 below compares the CERES TOA SRI with the CFSR tropics (30N-30S) mean surface temperature, which covers the 50% of the global surface that receives most of the incoming solar energy.

Figure 4. Comparison of CERES TOA SRI with CFSR Tropics Mean Surface Temperature (click to enlarge)

The not quite linear cyclical nature of the relationship between SRI and GMST is illustrated by the scatter plot in Figure 5 below.

Figure 5. Scatter plot of CFSR GMST versus CERES TOA SRI
(click to enlarge)

By lagging the GMST by 6 months, the relationship to TOA SRI looks more meaningful, as seen in Figure 6 below.

Figure 6. Scatter plot of CFSR GMST lagged 6 months versus CERES TOA SRI
(click to enlarge)

I calculated and compiled annual statistics for each year.  During the study period, the 2001-2017 TOA SRI average of the 17 annual averages was 340.023 watts per meter squared (watts/m2) with a standard deviation of only 0.087 watts/m2.  The monthly averages ranged from 329.113 to 351.516 watts/m2.  The annual range in monthly TOA SRI averaged 21.998 watts/m2 with a standard deviation of only 0.106 watts/m2.  The highest annual range was 21.838 and lowest 22.261 watts/m2.  The average percentage of the annual range relative to the lowest month each year was 6.682%.  Note that I did not bother to weight the annual averages based on number of days per month because I do not expect a significant difference.

Similarly, the 2001-2017 average of the 17 annual averages of GMST was 287.810 Kelvin (K) with a standard deviation of 0.133 K and a range in monthly averages from 285.505 K to 289.864 K.  The lowest annual average was 287.735 K in 2008 and the highest was 288.051 K in 2016.  The annual range in monthly GMST averaged 3.775 K with a standard deviation of 0.166 K and varied from 3.511 K to 4.128 K.

For 2001-2017, the CFSR monthly surface temperature data shows an average annual range of 12.589 K for the northern hemisphere and an average annual range of 5.227 K for the southern hemisphere. The southern hemisphere cycle is only one month delayed from the TOA SRI cycle. Both of the hemispheric cyclical temperature swings are larger than the global swing and include both seasonal earth tilt plus eccentricity effects. The global swing should be the net remainder resulting from the eccentricity induced TOA SRI effect.

Initial Thoughts

So, if we assume that the annual range in TOA SRI is the primary driver for the annual range in GMST, then we can calculate that GMST rises (or falls) by 3.775/21.998 = 0.172 K per 1 watt/m2 of TOA SRI change.  The implication is that if a doubling of CO2 causes a radiative forcing of 3.7 watts/m2, the corresponding rise in GMST would be 3.7 x 0.172 = 0.635 K and this would include any feedbacks that normally occur in the earth system over an annual cycle.  The main difference is that the ramp up in solar input from orbital eccentricity occurs over only a 6 month period each year, whereas the doubling of CO2 might take on the order of a century.  However, both are responses to changes to the earth’s radiation balance.  So far, I have not been able to think of any compelling reasons why the difference in time scales would make much difference in the resulting effect on GMST.  Possibly the short period cycling of the solar input may not allow enough time to reach full impact on GMST in either direction, leading to a quasi-steady state oscillating result?

I’m afraid this assessment may be too simplistic or that I may have overlooked some important influences.  Thus, I’m not at all certain this approach is a valid method for estimating the effect on GMST from a doubling of CO2.  I will be interested hear what readers have to say and I would be especially interested in learning how well the global climate models handle this annual cycling of both TOA SRI and GMST.

After Further Review

I’m now seeing that the annual cycle in global temperature is dominated by the seasonal cycle related to the earth’s axial tilt.  I thought it might still be possible to estimate the effects of eccentricity on the seasonal cycle in each hemisphere, and from that result, estimate the net effect on the global temperature cycle from eccentricity.  However, I found that backing out a rough estimate of the effect of eccentricity actually increased the global annual temperature range, because it increases the NH annual temperature range more than it decreases the SH annual temperature range.  Consequently it does not appear possible to determine the effect of eccentricity on global temperature without using a climate model to test the effect of varying degrees of eccentricity.  Like many things I’m learning about climate … it is very complex.


Eccentricity discussion:

CERES project description:’s_Radiant_Energy_System

CERES data download:

UM CCI Reanalyzer:

CDAS monthly average downloads:

Climate sensitivity discussion:

Earth’s energy budget discussion:

Arctic Buoy Decline

There has been a sharp decline in the number of Arctic buoys reporting temperature measurements to the global synoptic weather network the last couple of years, as compared to the previous several years. Below are comparisons of synoptic temperature observations for 1200 Universal Time Coordinates (UTC) on January 1 for each year from 2019 back to 2015 for the Arctic Ocean area. The first set of maps below, Figures 1 through 5,  are standard plots of weather data, including temperature (upper left from station circle), from the National Oceanographic and Atmospheric Administration (NOAA) Weather Prediction Center (WPC) for the southwestern portion of the Arctic Ocean.  Notice the increasing number of open ocean buoys reporting going back each year.  Click on any of the maps to see the full size image.

Figure 1. Synoptic weather observations for 2019 January 1 at 1200 UTC.

Figure 2. Synoptic weather observations for 2018 January 1 at 1200 UTC.

Figure 3. Synoptic weather observations for 2017 January 1 at 1200 UTC.

Figure 4. Synoptic weather observations for 2016 January 1 at 1200 UTC.

Figure 5. Synoptic weather observations for 2015 January 1 at 1200 UTC.

The next set of maps, Figures 6 through 10, show plots of synoptic temperature observations from OGIMET for the entire Arctic Ocean area for January 1 at 1200 UTC for each year from  2019 back to 2015.  Again notice the increasing number of temperature observations from the open Arctic Ocean going back each year.    Click on any of the maps to see the full size image.

Figure 6. Synoptic temperature observations for 2019 January 1 at 1200 UTC for the Arctic Ocean.

Figure 7. Synoptic temperature observations for 2018 January 1 at 1200 UTC for the Arctic Ocean.

Figure 8. Synoptic temperature observations for 2017 January 1 at 1200 UTC for the Arctic Ocean.

Figure 9. Synoptic temperature observations for 2016 January 1 at 1200 UTC for the Arctic Ocean.

Figure 10. Synoptic temperature observations for 2015 January 1 at 1200 UTC for the Arctic Ocean.

Oddly, the most recent International Arctic Buoy Program (IABP) map for 2018 December 27 shows numerous buoys reporting air temperature from the Arctic Ocean area as seen in Figure 11 below.

Figure 11. IABP map of air temperature observations for 2018 December 27.

Why the temperature measurements from these buoys are not being reported to the global synoptic weather network is puzzling, especially considering that much of the recent global warming has been occurring primarily in the winter night-time Arctic area.  I am also not certain whether any of these buoy measurements are being ingested into global weather models and associated reanalyses, separately from synoptic weather data.  It would be a shame if they are not.

The IABP also provides data and graphs of data from the Arctic Buoys here.   In the past I have compared the detailed buoy measurements from IABP to OGIMET plotted synoptic observations and code that were reported from some of the buoys.  I found that in some cases the IABP reported “surface temperature” from underneath the hull of the buoy was being erroneously reported in the synoptic data as “air temperature”.  Graphs of the IABP data indicate that the “surface temperature” (which is typically from ice or water beneath the buoy) is often substantially different than the “air temperature” when both are reported.  I also found flat-lined air temperatures and air temperatures that did not match well with those from nearby buoys in both data sets.  But these are different problems to investigate for another day.

Happy New Year!

Global Temperature January 2018

Notice:  the graph in this post has been moved to the Monthly Trends page accessible in the menu bar at the top of this page, where it will be updated monthly rather than by monthly posts.

Climate Forecast System Reanalysis (CFSR) monthly global surface temperature anomaly estimates for 2014 through January 2018 from the University of Maine Climate Change Institute (UM CCI) and from WeatherBELL (WxBELL) are graphed below along with monthly global temperature anomaly estimates for the lower troposphere derived from satellite measurements provided by the University of Alabama at Huntsville (UAH).  The UM CCI CFSR estimates have been adjusted (UM adj), while the WxBELL CFSR estimates have been left unadjusted to show the difference.  The UM CCI CFSR adjusted monthly estimates for June through December 2017 are based on final daily averages and the estimate for January 2018 is based on preliminary daily averages.  Most likely the UM CCI January estimate will drop by about 0.10C once final daily averages for January are released. The UM CCI January preliminary estimate at +0.36C was down by 0.04C from December while the WxBELL estimate at +0.26C was down by 0.13C and the UAH estimate also at +0.26C was down by 0.15C. Click on the graph below to see a larger copy.

Also shown for comparison are monthly global temperature anomaly estimates from six other major sources, including lower tropospheric estimates from the Remote Sensing Systems (RSS), and surface estimates from the European Centre for Medium-Range Weather Forecast (ECMWF) Reanalysis Interim adjusted (ERAI adj), US National Center for Environmental Information (NCEI), US National Aeronautics and Space Administration (NASA) Goddard Institute of Space Studies (GISS), the UK Hadley Climate Research Unit Temperature (CRUT), and the Berkeley Earth Surface Temperature (BEST), all final through December 2017. All estimates have been shifted to the latest climatological reference period 1981-2010.

Update 2018 February 6

Final January 2018 global temperature anomaly estimates from RSS for lower troposphere and ERAI adj for surface have been added to the graph.   Both showed decreases from December 2017 to January 2018, with RSS down 0.04C,  and ERAI adj down 0.14C.

GOES-16 Preview

The first next generation US Geostationary Operational Environmental Satellite (GOES) was recently launched by the US National Aeronautic and Space Administration (NASA) on November 19th of 2016, designated as GOES-R.  This new series of satellites will provide 34 meteorological, solar and space weather products.  They orbit above a fixed point at the earth’s equator at a distance of about 22,300 miles out in space.   As with all US meteorological satellites, the US National Oceanic and Atmospheric Administration (NOAA) has taken over operation of the satellite and designated it GOES-16.  Below is a link to a visible composite color high resolution full-disk test image from midday January 15th of 2017 provided by NOAA.  To see the image in full resolution, click on the reduced image below which will take you to NOAA’s web site to view the full resolution image (use scroll bars or browser magnification tool to navigate) and you can return here by using your browser back button.

Additional test images can be seen here:
GOES-16 Image Gallery

Here is an animation showing the 16 different imagery channels available:

Below is a description of the satellite and its uses.

In May 2017, NOAA will announce the new location for GOES-16.  It will replace either GOES-East or GOES-West and will become operational in November 2017.  The next satellite in the series, GOES-S, is scheduled for launch in spring 2018 and should be operational by a year later.

Information about data access can be found here:
GOES-R User Systems

NASA also has a useful web page for viewing real-time and archived high resolution imagery from three polar orbiting satellites here:
NASA Worldview