first record of thermospheric ir cooling rates over a complete solar cycleir cooling in current solar maximum conditions much weaker than prior maximumvariability in thermospheric ir cooling observed on scale of days to 11 years 
  first record of thermospheric ir cooling rates over a complete solar cycle ir cooling in current solar maximum conditions much weaker than prior maximum variability in thermospheric ir cooling observed on scale of days to 11 years
the sounding of the atmosphere using broadband emission radiometry ( saber ) instrument [ russell et al .
, 1999 ; mlynczak , 1997 ] on the nasa thermosphere - ionosphere - mesosphere energetics and dynamics ( timed ) satellite in late january 2014 marked 12  years of continuous observations of the atmosphere from the tropopause up to the middle thermosphere .
in particular , saber has provided continuous observations of the infrared cooling rates in the lower and middle thermosphere ( 100 to as high as 300  km altitude ) due to emission from the vibration - rotation bands of carbon dioxide ( co2 , 15  m ) and nitric oxide ( no , 5.3  m ) .
these data represent the first continuous , global observations of the two major radiative cooling mechanisms of the thermosphere over a nominal complete solar cycle .
radiative cooling by co2 and no depends primarily on three factors : the kinetic temperature , the atomic oxygen abundance , and the abundance of no or co2 . in the thermosphere , excitation of the vibrational modes of each molecule
occurs largely by collisions with atomic oxygen ( o ) [ crutzen , 1970 ; kockarts , 1980 ] .
the vibrationally excited molecules may then radiate spontaneously or may be quenched , again largely due to collisions with o. the rate of collisional excitation depends on the product of the o abundance and the no or co2 abundance .
collisional excitation is also strongly temperature dependent and is proportional to exp(e / kt ) , where e is the upper state energy , k is boltzmann 's constant , and t is the kinetic temperature .
variations in cooling are thus linear in o , no , or co2 , and nonlinear in temperature .
the o and no abundances are sensitive to variations in solar irradiance and to geomagnetic effects , the latter of which can be quite substantial for no . in addition , co2 is expected to be steadily increasing due to anthropogenic activity . above 100  km the radiative transfer for both co2 and no is effectively optically thin and any emitted radiation escapes either to space or to the lower atmosphere .
the fine structure line of atomic oxygen at 63  m also contributes to the radiative cooling of the thermosphere [ bates , 1951 ] .
the temperature dependence of cooling by o at 63  m is proportional to exp(228/t ) .
consequently , infrared cooling by o changes very little even with large changes in temperature because of the relatively small argument ( i.e. , small energy at 63  m ) in the exponential [ mlynczak et al . ,
this letter therefore concentrates on the substantial observed variations in cooling by no and co2 made by the saber instrument .
saber commenced scientific operations in january 2002 after a successful launch into orbit on the timed spacecraft in december 2001 . during the last 12  years saber
has chronicled the effects of solar variability on the infrared radiative cooling of the thermosphere on a variety of time scales .
the large increases in radiative cooling of the thermosphere in response to coronal mass ejections and particle precipitation , now referred to as the natural thermostat effect , were documented early in the saber mission [ mlynczak et al .
prior to saber observations , the response of the thermosphere to particle precipitation has been extensively studied [ e.g. , grard and barth , 1977 ; siskind et al . , 1989 ] .
subsequent to these initial saber observations , mlynczak et al .   documented the decrease in the radiative cooling by no associated with the start of the declining phase of solar cycle 23 .
short - term periodic features corresponding to the first four harmonics of the 27 day solar rotation period were observed in the radiative cooling by no during solar minimum [ mlynczak et al . , 2008 , 2010a ] .
high - speed solar wind streams emanating from multiple coronal holes in the sun 's atmosphere caused these features .
although not shown in this paper , these short - term periodic features have not reappeared from 2009 to 2013 .
[ 2010b ] further documented the occurrence of solar minimum and demonstrated that cooling by co2 dominates the radiative cooling of the thermosphere at that time .
noted the onset of solar cycle 24 as evidenced by increasing rates of radiative cooling by both no and co2 simultaneously with increasing rates of solar activity . in this paper
we present the continuous 12 year saber record of radiative cooling of the thermosphere and discuss the effects of solar variability now fully evident in this unique data set .
our intent is to document , for the first time , a complete record of thermospheric infrared cooling by no and co2 over a nominal complete solar cycle .
solar cycle 24 is now being studied for its relative weakness compared to prior solar cycles .
concurrent observations of changes in solar ultraviolet irradiance and in geomagnetic activity also correlate with the observed changes in cooling and are discussed below .
the next section briefly discusses the methodology used to derive radiative cooling rates from the saber limb radiance measurements at 5.3  m and 15  m .
results are presented in section 3 and a summary and discussion are given in section 4 .
the saber instrument is an infrared radiometer that records vertical profiles of limb radiance in 10 distinct spectral channels . from these data , vertical profiles of kinetic temperature , atmospheric composition , and rates of radiative cooling by no and
the thermospheric cooling rate profiles are in units of nw / m ( as opposed to k / d ) because saber does not measure atmospheric density over the depth of the thermosphere .
these vertical cooling rate profiles can be integrated with respect to altitude and then with respect to latitude and longitude to compute the total infrared power ( w ) radiated by no and by co2 over the entire earth on a given day . to illustrate the influence of the solar cycle
, we will use time series of both the daily global infrared power and the global annual mean cooling rate profiles over the last 12  years .
we begin by showing the daily time series of global infrared power ( w ) radiated by no and co2 for the past 12  years in figure  1 .
the time series covers 22 january 2002 through 11 march 2014 , over 4400  days of data .
saber observes all local times over 60  days as a consequence of the 74 inclined orbit of the timed spacecraft .
the 11 year solar cycle is visually evident in the data , with maximum emission in 2002 , minimum in 2008/2009 , and rising into the current maximum extending from 2011 to the present ( march 2014 ) .
the deep minimum in cooling by no is evident as the radiated power decreases by nearly an order of magnitude from 2002 to 2008 .
the increase in both no and co2 cooling at the onset of solar cycle 24 is evident early in 2011 [ hunt et al . ,
both the no and co2 data exhibit high frequency  spikes  in cooling which are the result of geomagnetic activity [ mlynczak et al . , 2003 ] and possibly solar flare events . the co2 cooling exhibits periodicities corresponding to the first three harmonics of the annual cycle [ mlynczak et al . ,
2008 ] . these longer periods in co2 are suggestive of forcing from the lower atmosphere weather .
the 60 day running mean curves for both no and co2 clearly indicate that the current maxima conditions of solar cycle 24 are substantially weaker than the prior maximum conditions in 20022003 .
peak co2 emission in 2002 is about 15% larger than in 20112014 , while peak no emission is about 1.75 times larger .
time series of saber daily global infrared power for ( top ) no and ( bottom ) co2 , from 22 january 2002 to 11 march 2014 .
data from more than 4400  days of saber observations are in each series . as discussed in mlynczak et al . [ 2010b ] , the primary driver of changes in co2 cooling over the solar cycle is changes in temperature , although co2 increases over the 12  years likely contribute as well .
changes in both temperature and no abundance drive the changes in no cooling over the solar cycle . variability or changes in atomic oxygen may also influence both no and co2 cooling rates .
temperature and composition changes in turn are driven by changes in solar irradiance and by geomagnetic activity that enhances auroral energy deposition and that drives composition changes in the thermosphere [ e.g. , richards , 2004 ] .
this time series of daily global radiated power has proven to be a very powerful tool to examine the response of the atmosphere to short - term variability of the sun .
in particular , the data are now being used to diagnose problems with forecasting geomagnetic storm intensity [ knipp et al . , 2013 ] and to study effects of corotating interaction regions on the ionosphere [ verkhoglyadova et al . , 2013 ] . to further investigate the influence of solar variability , we have computed the deseasonalized no and co2 power since the beginning of the mission .
every 60  days , the timed spacecraft is rotated 180 about its yaw axis to maintain thermal stability .
the mean global power is computed for no and co2 for each spacecraft yaw cycle ( which is also the time required to precess through all local times ) , with 6 yaw cycles in each  yaw year . 
the yaw cycle repeats annually from mid - january in one year to mid - january in the next .
this results in 72 ( 12  years by 6  cycles per year ) 60 day global means . next
the average power of the first yaw cycle in each year is computed , followed by the average power of the second yaw cycle , then the average of the third yaw cycle , etc .
this gives six mean powers ( over 12  years ) , one for each of 6 annual yaw cycles .
the mean power is subtracted from the corresponding individual yaw cycle average power for each of the 12  years , i.e. , the mean of yaw cycle 1 from each of the 12 first yaw cycles , and so forth , resulting in the deseasonalized no and co2 powers shown in figures  2 and 3 .
for each 60 day period we also computed the difference between the mean ultraviolet solar irradiance measured by the solar euv experiment ( see ) instrument [ woods et al . ,
the spectral interval chosen for the see data is 1190  nm that includes the lyman - alpha line and the schumann - runge continuum that heat the lower thermosphere .
( top ) deseasonalized ( over period of 60  days ) co2 thermospheric infrared power and ( bottom ) see 60 day radiative anomalies ( relative to the see long - term mean ) .
( top ) deseasonalized ( over period of 60  days ) no thermospheric infrared power and ( bottom ) anomalies in the kp geomagnetic index ( relative to the kp long - term mean ) .
correlations with geomagnetic activity appear stronger during the 20022003 maximum than during the current weaker maximum period .
the results are shown in figure  2 for co2 and the see irradiance . in this figure
we clearly see the progression of the solar cycle in the co2 power ( green bars ) from january 2002 to january 2014 .
a peak in deseasonalized power occurred in late 2011 , corresponding to a peak in the deseasonalized see irradiance .
additional , larger peaks in co2 emission occurred in early and late 2013 . in a similar way figure  3
shows the deseasonalized no power and the anomalies in the 60 day mean kp index relative to the long - term mean of the kp index .
the deseasonalized no power also shows a peak in late 2011 and two peaks in 2013 that correspond to the peaks in co2 power . the visual correlation between see and co2 emission and kp and no emissions is much higher during the stronger solar maximum conditions of 20022003 than those observed between 2011 and 2014 .
these data can also be used to compute the difference in infrared energy emitted by no and co2 at different times . in 2002
no and co2 combined to radiate 7    10 more joules of infrared energy than in 2013 , with each gas accounting for about one half of the total difference .
the daily infrared radiative power results above have proven extremely useful in studying the influence of the variable sun on the earth 's thermosphere .
however , the power data offer no information on the variability as a function of altitude . for this purpose
we present in figures  4 and 5 the time series of the difference between the global annual mean cooling rate profile in 2002 and that observed each year from 2003 to 2013 .
vertical profiles of the differences in cooling ( nw / m ) are shown as a function of yaw year .
figure  4 indicates that the no cooling rates decreased until the middle or latter part of 2008 and then began increasing until 2012 , when they appear to level off through 2013 .
figure  5 shows the difference in the vertical profile of global annual mean radiative cooling by co2 from that observed in 2002 .
again , negative values indicate lower radiative cooling than in 2002 . as with the no cooling ,
the maximum change occurs in the middle to latter part of 2008 but at a much lower altitude of about 105  km . the cooling rates by co2 decreased by 25%30% at the peak altitude from those observed in 2002 .
in contrast , the no cooling rates decreased by more than 85% at the altitude of peak change .
time series of the difference in the vertical profile of global annual average no cooling rate by year from 2002 .
time series of the difference in the vertical profile of global annual average co2 cooling rate by year from 2002 .
we begin by showing the daily time series of global infrared power ( w ) radiated by no and co2 for the past 12  years in figure  1 .
the time series covers 22 january 2002 through 11 march 2014 , over 4400  days of data .
saber observes all local times over 60  days as a consequence of the 74 inclined orbit of the timed spacecraft .
the 11 year solar cycle is visually evident in the data , with maximum emission in 2002 , minimum in 2008/2009 , and rising into the current maximum extending from 2011 to the present ( march 2014 ) .
the deep minimum in cooling by no is evident as the radiated power decreases by nearly an order of magnitude from 2002 to 2008 .
the increase in both no and co2 cooling at the onset of solar cycle 24 is evident early in 2011 [ hunt et al . ,
both the no and co2 data exhibit high frequency  spikes  in cooling which are the result of geomagnetic activity [ mlynczak et al . , 2003 ] and possibly solar flare events . the co2 cooling exhibits periodicities corresponding to the first three harmonics of the annual cycle [ mlynczak et al . ,
2008 ] . these longer periods in co2 are suggestive of forcing from the lower atmosphere weather .
the 60 day running mean curves for both no and co2 clearly indicate that the current maxima conditions of solar cycle 24 are substantially weaker than the prior maximum conditions in 20022003 .
peak co2 emission in 2002 is about 15% larger than in 20112014 , while peak no emission is about 1.75 times larger .
time series of saber daily global infrared power for ( top ) no and ( bottom ) co2 , from 22 january 2002 to 11 march 2014 .
data from more than 4400  days of saber observations are in each series . as discussed in mlynczak et al . [ 2010b ] , the primary driver of changes in co2 cooling over the solar cycle is changes in temperature , although co2 increases over the 12  years likely contribute as well .
changes in both temperature and no abundance drive the changes in no cooling over the solar cycle . variability or changes in atomic oxygen may also influence both no and co2 cooling rates .
temperature and composition changes in turn are driven by changes in solar irradiance and by geomagnetic activity that enhances auroral energy deposition and that drives composition changes in the thermosphere [ e.g. , richards , 2004 ] .
this time series of daily global radiated power has proven to be a very powerful tool to examine the response of the atmosphere to short - term variability of the sun .
in particular , the data are now being used to diagnose problems with forecasting geomagnetic storm intensity [ knipp et al . , 2013 ] and to study effects of corotating interaction regions on the ionosphere [ verkhoglyadova et al . , 2013 ] . to further investigate the influence of solar variability , we have computed the deseasonalized no and co2 power since the beginning of the mission .
every 60  days , the timed spacecraft is rotated 180 about its yaw axis to maintain thermal stability .
the mean global power is computed for no and co2 for each spacecraft yaw cycle ( which is also the time required to precess through all local times ) , with 6 yaw cycles in each  yaw year . 
the yaw cycle repeats annually from mid - january in one year to mid - january in the next .
this results in 72 ( 12  years by 6  cycles per year ) 60 day global means . next
the average power of the first yaw cycle in each year is computed , followed by the average power of the second yaw cycle , then the average of the third yaw cycle , etc .
this gives six mean powers ( over 12  years ) , one for each of 6 annual yaw cycles .
the mean power is subtracted from the corresponding individual yaw cycle average power for each of the 12  years , i.e. , the mean of yaw cycle 1 from each of the 12 first yaw cycles , and so forth , resulting in the deseasonalized no and co2 powers shown in figures  2 and 3 .
for each 60 day period we also computed the difference between the mean ultraviolet solar irradiance measured by the solar euv experiment ( see ) instrument [ woods et al . ,
the spectral interval chosen for the see data is 1190  nm that includes the lyman - alpha line and the schumann - runge continuum that heat the lower thermosphere .
( top ) deseasonalized ( over period of 60  days ) co2 thermospheric infrared power and ( bottom ) see 60 day radiative anomalies ( relative to the see long - term mean ) .
( top ) deseasonalized ( over period of 60  days ) no thermospheric infrared power and ( bottom ) anomalies in the kp geomagnetic index ( relative to the kp long - term mean ) .
correlations with geomagnetic activity appear stronger during the 20022003 maximum than during the current weaker maximum period .
the results are shown in figure  2 for co2 and the see irradiance . in this figure
we clearly see the progression of the solar cycle in the co2 power ( green bars ) from january 2002 to january 2014 .
a peak in deseasonalized power occurred in late 2011 , corresponding to a peak in the deseasonalized see irradiance .
additional , larger peaks in co2 emission occurred in early and late 2013 . in a similar way figure  3
shows the deseasonalized no power and the anomalies in the 60 day mean kp index relative to the long - term mean of the kp index .
the deseasonalized no power also shows a peak in late 2011 and two peaks in 2013 that correspond to the peaks in co2 power . the visual correlation between see and co2 emission and kp and no emissions is much higher during the stronger solar maximum conditions of 20022003 than those observed between 2011 and 2014 .
these data can also be used to compute the difference in infrared energy emitted by no and co2 at different times . in 2002
no and co2 combined to radiate 7    10 more joules of infrared energy than in 2013 , with each gas accounting for about one half of the total difference .
the daily infrared radiative power results above have proven extremely useful in studying the influence of the variable sun on the earth 's thermosphere . however , the power data offer no information on the variability as a function of altitude . for this purpose
we present in figures  4 and 5 the time series of the difference between the global annual mean cooling rate profile in 2002 and that observed each year from 2003 to 2013 .
vertical profiles of the differences in cooling ( nw / m ) are shown as a function of yaw year .
figure  4 indicates that the no cooling rates decreased until the middle or latter part of 2008 and then began increasing until 2012 , when they appear to level off through 2013 .
figure  5 shows the difference in the vertical profile of global annual mean radiative cooling by co2 from that observed in 2002 .
again , negative values indicate lower radiative cooling than in 2002 . as with the no cooling ,
the maximum change occurs in the middle to latter part of 2008 but at a much lower altitude of about 105  km . the cooling rates by co2 decreased by 25%30% at the peak altitude from those observed in 2002 .
in contrast , the no cooling rates decreased by more than 85% at the altitude of peak change .
time series of the difference in the vertical profile of global annual average no cooling rate by year from 2002 .
time series of the difference in the vertical profile of global annual average co2 cooling rate by year from 2002 .
we have presented 12  years of radiative cooling data associated with no and co2 in the earth 's lower and middle thermosphere .
the influence of the variability of the sun during solar cycle 23 and the first part of solar cycle 24 are evident in these data . a long and deep solar minimum previously noted
is now observed in its entirety , and solar maximum conditions that are substantially weaker than the prior maximum conditions are now occurring .
the deseasonalized vertically integrated cooling rates ( radiated power ) appear to be more strongly correlated with changes in solar ultraviolet irradiance and geomagnetic activity during the prior maximum than during the current maximum conditions .
the influence of solar variability is evident in both no and co2 cooling throughout the entire depth of the lower and middle thermosphere .
the saber radiative cooling data for no and co2 constitute a unique climate data record for testing upper atmosphere models in terms of both radiative / chemical physics and response to solar variability .
one key application would be to the ongoing debate [ e.g. , solomon et al . ,
2013 ] regarding the causes and consequences of the last solar minimum in 20082009 during which record low thermospheric densities occurred [ emmert et al . , 2010 ] .
solomon et al .   attributed the low density to changes in solar ultraviolet irradiance ,
while emmert et al .   conclude that other factors including composition changes must play a role .
the saber radiative cooling data can be used to test different scenarios on composition ( atomic oxygen ) changes , temperature changes , and co2 changes to further resolve this enigma .
 infrared radiative cooling of the thermosphere by carbon dioxide ( co2 , 15  m ) and by nitric oxide ( no , 5.3  m ) has been observed for 12  years by the sounding of the atmosphere using broadband emission radiometry ( saber ) instrument on the thermosphere - ionosphere - mesosphere energetics and dynamics satellite . 
 for the first time we present a record of the two most important thermospheric infrared cooling agents over a complete solar cycle . 
 saber has documented dramatic variability in the radiative cooling on time scales ranging from days to the 11 year solar cycle . 
 deep minima in global mean vertical profiles of radiative cooling are observed in 20082009 . 
 current solar maximum conditions , evidenced in the rates of radiative cooling , are substantially weaker than prior maximum conditions in 20022003 . 
 the observed changes in thermospheric cooling correlate well with changes in solar ultraviolet irradiance and geomagnetic activity during the prior maximum conditions . 
 no and co2 combine to emit 7    1018 more joules annually at solar maximum than at solar minimum.key pointsfirst record of thermospheric ir cooling rates over a complete solar cycleir cooling in current solar maximum conditions much weaker than prior maximumvariability in thermospheric ir cooling observed on scale of days to 11 years 
 
Key Points
1. Introduction
2. Methodology
3. Results
3.1. Daily Global Power
3.2. Global Annual Mean Cooling Rate Profiles
4. Summary and Discussion