Measurements of outgoing longwave radiation (OLR) are essential for understanding many aspects of climate. Many people are confused about the factors that affect OLR. And its rich variability is often not appreciated.
There have been a number of satellite projects since the late 1970’s, with the highlight (prior to 2001) being the five year period of ERBE.
AIRS & CERES were launched on the NASA AQUA satellite in May 2002. These provide much better quality data, with much better accuracy and resolution.
CERES has three instruments:
- Solar Reflected Radiation (Shortwave): 0.3 – 5.0 μm
- Window: 8 – 12 μm
- Total: 0.3 to > 100 μm
AIRS is an infrared spectrometer/radiometer that covers the 3.7–15.4 μm spectral range with 2378 spectral channels. It runs alongside two microwave instruments (better viewing through clouds): AMSU is a 15-channel microwave radiometer operating between 23 and 89 GHz; HSB is a four-channel microwave radiometer that makes measurements between 150 and 190 GHz.
From Aumann et al (2003):
The simultaneous use of the data from the three instruments provides both new and improved measurements of cloud properties, atmospheric temperature and humidity, and land and ocean skin temperatures, with the accuracy, resolution, and coverage required by numerical weather prediction and climate models.
Among the important datasets that AIRS will contribute to climate studies are as follows:
- atmospheric temperature profiles;
- sea-surface temperature;
- land-surface temperature and emissivity;
- relative humidity profiles and total precipitable water vapor;
- fractional cloud cover;
- cloud spectral IR emissivity;
- cloud-top pressure and temperature;
- total ozone burden of the atmosphere;
- column abundances of minor atmospheric gases such as CO, CH, CO, and N2O;
- outgoing longwave radiation and longwave cloud radiative forcing;
- precipitation rate
More about AIRS = Atmospheric Infrared Sounder, at Wikipedia, plus the AIRS website.
More about CERES = Clouds and the Earth’s Radiant Energy System, at Wikipedia, plus the CERES website – where you can select and view or download your own data.
How do CERES & AIRS compare?
CERES and AIRS have different jobs. CERES directly measures OLR. AIRS measures lots of spectral channels that don’t cover the complete range needed to just “add up” OLR. Instead, OLR can be calculated from AIRS data by deriving surface temperature, water vapour concentration vs height, CO2 concentration, etc and using a radiative transfer algorithm to determine OLR.
Here is a comparison of the two measurement systems from Susskind et al (2012) over almost a decade:
Figure 1
The second thing to observe is that the measurements have a bias between the two datasets. But because we have two high accuracy measurement systems on the same satellite we do have a reasonable opportunity to identify the source of the bias (total OLR as shown in the graph is made of many components). If we only had one satellite, and then a new satellite took over with a small time overlap any biases would be much more difficult to identify. Of course, that doesn’t stop many people from trying but success would be much harder to judge.
In this paper, as we might expect, the error sources between the two datasets get considerable discussion. One important point is that version 6 AIRS data (prototyped at the time the paper was written) is much closer to CERES. The second point, probably more interesting, is that once we look at anomaly data the results are very close. We’ll see a number of comparisons as we review what the paper shows.
The authors comment:
Behavior of OLR over this short time period should not be taken in any way as being indicative of what long-term trends might be. The ability to begin to draw potential conclusions as to whether there are long-term drifts with regard to the Earth’s OLR, beyond the effects of normal interannual variability, would require consistent calibrated global observations for a time period of at least 20 years, if not longer. Nevertheless, a very close agreement of the 8-year, 10-month OLR anomaly time series derived using two different instruments in two very different manners is an encouraging result.
It demonstrates that one can have confidence in the 1° x 1° OLR anomaly time series as observed by each instrument over the same time period. The second objective of the paper is to explain why recent values of global mean, and especially tropical mean, OLR have been strongly correlated with El Niño/La Niña variability and why both have decreased over the time period under study.
Why Has OLR Varied?
The authors define the average rate of change (ARC) of an anomaly time series as “the slope of the linear least squares fit of the anomaly time series”.
We can see excellent correlation between the two datasets and we can see that OLR has, on average, decreased over this time period.
Below is a comparison with the El Nino index.
We define the term El Niño Index as the difference of the NOAA monthly mean oceanic Sea Surface Temperature (SST), averaged over the NOAA Niño-4 spatial area 5°N to 5°S latitude and 150°W westward to 160°E longitude, from an 8-year NOAA Niño-4 SST monthly mean climatology which we generated based on use of the same 8 years that we used in the generation of the OLR climatologies.
Figure 2
It gets interesting when we look at the geographical distribution of the OLR changes over this time period:
Figure 3 – Click to Enlarge
We see that the tropics have the larger changes (also seen clearly in figure 2) but that some regions of the tropics have strong positive values and other regions have strong negative values. The grey square square centered on 180 longitude is the Nino-4 region. Values as large as +4 W/m²/decade are found in this region. And values as large as -3 W/m²/decade are found over Indonesia (WPMC region).
Let’s look at the time series to see how these changes in OLR took place:
Figure 4 – Click to Enlarge
The main parameters which affect changes in OLR month to month and year to year are a) surface temperatures b) humidity c) clouds. As temperature increases, OLR increases. As humidity and clouds increase, OLR decreases.
Here are the changes in surface temperature, specific humidity at 500mbar and cloud fraction:
Figure 5 – Click to Enlarge
So, focusing again on the Nino-4 region, we might expect to find that OLR has decreased because of the surface temperature decrease (lower emission of surface radiation) – or we might expect to find that the OLR has increased because the specific humidity and cloud fraction have decreased (thus allowing more surface and lower atmosphere radiation to make it through to TOA). These are mechanisms pulling in opposite directions.
In fact we see that the reduced specific humidity and cloud fraction have outweighed the effect of the surface temperature decrease. So the physics should be clear (still considering the Nino-4 region) – if surface temperature has decreased and OLR has increased then the explanation is the reduction in “greenhouse” gases (in this case water vapor) and clouds, which contain water.
Correlations
We can see similar relationships through correlations.
The term ENC in the graphs stands for El Nino Correlation. This is essentially the correlation of the time-series data with time-series temperature change in the Nino-4 region (more specifically the Nino-4 temperature less the global temperature).
As the Nino-4 temperature declined over the period in question, a positive correlation means the value declined, while a negative correlation means the value increased.
The first graph below is the geographical distribution of rate of change of surface temperature. Of course we see that the Nino-4 region has been declining in temperature (as already seen in figure 2). The second graph shows this as well, but also indicates that the regions west and east of the Nino-4 region have a stronger (negative) correlation than other areas of larger temperature change (like the arctic region).
The third graph shows that 500 mb humidity has been decreasing in the Nino-4 region, and increasing to the west and east of this region. Likewise for the cloud fraction. And all of these are strongly correlated to the Nino-4 time-series temperature:
Figure 6 – Click to expand
For OLR correlations with Nino-4 temperature we find a strong negative correlation, meaning the OLR has increased in the Nino-4 region. And the opposite – a strong positive correlation – in the highlighted regions to east and west of Nino-4:
Figure 7 – Click to expand
Note the two highlighted regions
- to the west: WPMC, Warm Pool Maritime Continent;
- and to the east: EEMA, Equatorial Eastern Pacific and Atlantic Region
We can see the correlations between the global & tropical OLR and the OLR changes in these regions:
Figure 8 – Click to expand
Both WPMC and EEPA regions together explain the reduction over 10 years in OLR. Without these two regions the change is indistinguishable from zero.
Conclusion
This article is interesting for a number of reasons.
It shows the amazing variability of climate – we can see adjacent regions in the tropics with completely opposite changes over 10 years.
It shows that CERES gets almost identical anomaly results (changes in OLR) to AIRS. CERES directly measures OLR, while AIRS retrieves surface temperature, humidity profiles, cloud fractions and “greenhouse” gas concentrations and uses these to calculate OLR.
AIRS results demonstrate how surface temperature, humidity and cloud fraction affect OLR.
OLR has – over the globe – decreased over 10 years. This is a result of the El-Nino phase – at the start of the measurement period we were coming out of a large El-Nino event, and at the end of the measurement period we were in a La Nina event.
The reduction in OLR is explained by the change in the two regions identified, which are themselves strongly correlated to the Nino-4 region.
References
Interannual variability of outgoing longwave radiation as observed by AIRS and CERES, Susskind et al, Journal of Geophysical Research (2012) – paywall paper
AIRS/AMSU/HSB on the Aqua Mission: Design, Science Objectives, Data Products, and Processing Systems, Aumann et al, IEEE Transactions on Geoscience and Remote Sensing (2003) – free paper
Four things, SoD:
First of all, thanks for the hat tip on the Susskind paper, much of the basis for this whole post.
Second, the ERBE project did not last for 5 years, but 14-15 (1985-99).
Third, you compare CERES and AIRS, noting that the former actually measures OLR and that the latter calculate it based on other ‘weather-based’ measurements, finding assurance in their striking similarity/correlation. But you forget to acknowledge the fact that this same assurance can be found when comparing the ISCCP FD with ERBE/ERBS and HIRS in the years between 1983-85 and 2000-02. (I know the post is not about that period, but just saying …)
Fourth and last, good to see that ENSO is given its rightful place in this whole issue. That process, fed by energy from the Sun, is after all what drives global temperatures. Not CO2.
Kristian,
Oh, please. ENSO drives cyclic variation in global temperatures. That’s cyclic variation around a long term mean. CO2 is a driver of the long term mean.
Kristian,
First – I assume a tongue in cheek thanks. Which I can understand from your point of view. Although from my point of view I had that paper and had started working on a post on that topic. Afterwards I found it was one of your papers.. well, hat tip to you.
Second – ERBE – I’ll come back to this one. As you haven’t done your own research I will provide mine.
Third – there is a difference between a principle and a value. You are comparing principles and implying that this means the values also are equivalent.
CERES and AIRS find similar but not exact OLR.
The reason that datasets like ISCCP FD exist is precisely because the measured datasets from that earlier era are not accurate enough from a number of points of view, just one of which is long term drift.
The same principle governs both parallel activities. My question is, what is the accuracy of the pre-CERES era?
If we look at the CERES era, a big step forward in measurement we find interesting commentary..
From Advances in Understanding Top-of-Atmosphere Radiation Variability from Satellite Observations, Loeb et al (2012):
In the (much-improved) CERES era, with 2 CERES instruments on 2 satellites and an AIRS instrument all in direct comparison we can identify 0.5 Wm-2 per decade.
Which is precisely the reason for my asking the questions I did ask and you didn’t answer. What do you believe is the accuracy (absolute and wrt drift) of the pre-CERES satellite observations?
On item 2 above:
The ERBE scanning radiometer data was available from 1984- 1989:
ERBE/ERBS Nov 1984 – Jan 1989
ERBE/NOAA 9 Jan 1985 – Jan 1987
ERBE/NOAA 10 Oct 1986 – May 1989
The ERBE non-scanner data was available from 1984-1999.
ERBE/ERBS Nov 1984 – Oct 1999
ERBE/NOAA 9 Jan 1985 – Jul 1990
ERBE/NOAA 10 Oct 1986 – Oct 1994
The scanner data has been the most useful for understanding climate because it provides OLR and Reflected SW as a function of location. The non-scanner data, by comparison, provides the total OLR via the “WFOV” = wide field of view, with no spatial breakdown.
The stated long-term drift accuracy for ERBE non-scanner OLR:
This comes from a source I can’t verify – ‘Appendix A’ which appears to be by Wong and Wielicki (two well-known people in this field).
This appendix does reference the later data correction from Reexamination of the Observed Decadal Variability of the Earth Radiation Budget Using Altitude-Corrected ERBE/ERBS Nonscanner WFOV Data, Wong et al, Journal of Climate (2006).
The point of referencing this paper and the quote is that stated accuracy is usually a best-case result. Then reality causes either a revised dataset, or a widened accuracy. It’s rare to find it the other way round, although I believe CERES provided better than expected accuracy.
And from The ScaRaB Earth Radiation Budget Dataset, Kandel et al, Bulletin of the American Meteorological Society (1998):
[Emphasis added]
And Comparison of measurements from satellite radiation budget instruments, Louis Smith et al, JGR (2006) has a lot to say about inter-satellite comparison, but in a way that I find hard to turn into decadal drift numbers.
Kristian,
I’m looking forward to your comment fifth and last, in response to my comment explaining why increases in GHGs, all other things being equal, would not be expected to cause a decadal reduction in OLR.
As far as I can understand your hypothesis: “..increases in GHGs will lead to a decadal reduction in OLR, which we haven’t seen..”
That’s my best estimate of your hypothesis. I probably got it wrong though.
One can obviously bias trends by picking the difference point in the ENSO cycle as starting and ending dates. Why isn’t it possible to pick starting and ending dates which represent similar states in the ENSO cycle?
Frank,
That’s a good question. In the case of this paper it is just using the longest time-series available.
The correlation of OLR with El-Niño seems pretty clear – or do you disagree?
For matching up the same states in the El Niño cycle I wonder whether this is possible. After all, it is a chaotic system. How do you define similar states in a chaotic system? We have El Ninos and La Ninas, but they don’t follow a rigid alternative pattern, they have very uneven strengths and durations..
So, the broadband measurements ( back to 1974 ) are probably not as accurate as subsequent systems. Still, it’s interesting to look at the series.
In that series, the two most notable El Nino years (’82-’83 and ’97-’98 ) don’t stand out at all in OLR:
http://www.climate4you.com/GlobalTemperatures.htm
Borne out better by this comparison:
And this plot of the tropics only:
Climate Weenie,
What is the source of the NOAA graphs of OLR?
Climate4you seems to get the data from knmi:
http://climexp.knmi.nl/selectfield_obs.cgi?someone@somewhere
With these references:
http://journals.ametsoc.org/doi/pdf/10.1175/1520-0477%281978%29059%3C1570%3AEARHBO%3E2.0.CO%3B2
http://www.esrl.noaa.gov/psd/data/gridded/data.interp_OLR.html#detail
Click to access OLRinterpolated.pdf
Problems likely abound with this dataset, but it’s probably the only longer term estimate of OLR.
SOD: The correlation with ENSO seems clear, but not how to address it. By choosing an appropriate time in the ENSO cycle, one can get all different kinds of answers about whether surface temperature has “significantly increased” over roughly the last decade, so I would hope they would have avoided that problem here.
To try to match states, one could try one of the various indexes of ENSO. Look at the value of the index for the last month of data and find the earliest month with a similar value. Then look at the value of the index for the first month of data and find the latest month with a similar value. Calculate the trend over both periods and report the range.
[…] 2013/02/07: TSoD: CERES, AIRS, Outgoing Longwave Radiation & El Nino […]
A random question for anyone:
Looking at graph b of figure 2 (tropical all sky OLR anomaly), what would be a good way to model this as some kind of stochastic process?
What question are you trying to answer?
Off the top of my head, an ARMA process with drift. That is an ARIMA (p,d,q) process with d constrained to be zero. and all of the AR coefficients must be less than one with a possible linear trend. I think you can make a very strong argument that no weather related process is a random walk. It could be fractionally integrated, but you don’t have a long enough time series to tell.
I’m building an interesting model. I’m not trying to answer any question about El Nino.
I’m adding OLR noise “in the style of” El Nino.
Right now I’ve got a couple of sine waves of different frequencies (8 yrs and 2.5 years) added plus some monthly noise and that’s actually not too bad.
But it’s too regular.
I think I need some frequency variation as well – but if you go from sin(2πft) to sin(2πt (f . random_f_factor)) you get a feature that changes as time increases – higher frequency components appearing later in the time series. There’s an obvious answer but I’ve only just arrived at this point.
Probably should be phase noise instead on the main component. I don’t think frequency noise makes sense.
Here’s a sample, with time in years on x-axis and OLR on the y-axis:
It’s not ideal. It will do for now. But if anyone has a suggestion for tropical OLR vs time “in the style of” El Nino then it would be appreciated.
You don’t want frequency noise, you want the period to vary. There’s a difference. It might be interesting to compare the model with El Nino data in the frequency domain.
SOD: You could check you output for realism by looking at the Fourier transform.
You might be able to take FT of the earth’s temperature record and create an alternative records using random phase shifts at each frequency
[…] CERES and AIRS – satellite instruments – are explained in CERES, AIRS, Outgoing Longwave Radiation & El Nino. […]