1 Introduction
Carbon dioxide is the major infrared cooler of the atmosphere from the lower stratosphere up to the lower thermosphere, where emission by nitric oxide becomes important . However, the CO infrared emissions in the bands near 15 that are responsible for the cooling are in non-local thermodynamic equilibrium (non-LTE) above around 70 km
The most used parameterization of the CO 15 cooling for Earth's middle and upper atmosphere was developed by . That parameterization is applicable for a limited range of CO abundances, up to double the pre-industrial CO concentration. Nowadays, however, with the rapid increase in the CO concentration in the atmosphere and its expected increase in the coming decades, climate model projections are being carried out in much higher CO scenarios that even quadruple the pre-industrial CO abundance . For example, such scenarios are considered in the Coupled Climate Model Intercomparison Projects. Therefore, parameterizations coping with such large CO concentrations are highly in demand, and that is precisely the purpose of our work. In addition to that of , other parameterizations of the CO 15 cooling rates have been developed in the past. In the case of Earth's atmosphere, it worth mentioning the comprehensive review of the early works reported by , the summary presented in Sect. 5.8 of and the more recent work of . Further, just before this work was submitted, a new parameterization was developed that uses the accelerated lambda iteration and opacity distribution function techniques . For the Mars and Venus atmospheres, where CO is the most abundant species, the problem has been tackled in several studies, e.g. , , and . In our case, we have the option of developing a completely new parameterization to adapt other CO parameterizations (like those cited above) or to extend and improve the parameterization of . Attending mainly to practical reasons of promptness, we opted for the latter.
The paper is structured as follows. A very basic description of the parameterization is presented in Sect. . Section describes the input atmospheric parameters used in the parameterization and required for calculating the reference cooling rates. In Sect. we describe the calculations of the reference LTE and non-LTE cooling rates. A detailed description of the parameterization is presented in Sect. . The testing and accuracy of the parameterization against (i) the reference cooling rates, (ii) the measured temperature fields of the middle atmosphere of MIPAS (Michelson Interferometer for Passive Atmospheric Spectroscopy) and (iii) the temperature profiles obtained by a high-resolution version of the Whole Atmosphere Community Climate Model (WACCM-X) are discussed in Sects. , and , respectively. The previous cooling rate parameterization was commonly used together with a parameterization of the CO near-infrared heating rates in GCMs. As we have extended the former to higher CO volume mixing ratios (VMRs) and we do not plan to extend the latter to higher CO VMRs in the near future, we assess in Sect. the validity of the current near-infrared heating parameterization for high CO VMRs. In Sect. , we summarize the main conclusions of the study.
2 Framework of the parameterization
As discussed above, this parameterization is essentially based on that of . To compute the CO cooling rate, the atmosphere is divided into five regions: one LTE and four different non-LTE regions. The method and approximations for computing the cooling rates in those regions are described in detail in Sect. . The new parameterization has a finer grid and, because it has been developed to cover a larger range of CO VMRs, the boundaries of the non-LTE regions were revised and, in general, their upper boundaries were shifted to higher altitudes. The scheme consists of 83 levels in ), covering to spaced by 0.25. The parameterization, however, can also be used for , where it uses the same scheme as for the NLTE4 region (see Sect. ). The relationship between pressure and geometrical altitude for the reference temperature profiles is shown in Fig. S1 in the Supplement. To a first approximation, the geometric altitude below km is related to by (km) .
Figure 1
Input data used in the reference calculations. (a) The six temperature profiles used in the reference calculations up to the lower thermosphere (Fig. shows them for the full range of pressures). MLE: mid-latitude equinox (September, 40° N); TRO: tropics (June, Equator); MLS: mid-latitude summer (June, 40° N); MLW: mid-latitude winter (June, 40° S); SAS: subarctic summer (June, 70° N); SAW: subarctic winter (June, 70° S). A profile typical of an elevated stratopause event (ES, mean of MIPAS temperature measured for latitudes of 70–90° N for 15 February 2009) is also shown for comparison and is discussed in Sect. . (b) CO volume mixing ratio profiles used in this work for the entire altitude range. Solid lines: the profiles used in the reference calculations; dashed lines: intermediate profiles used to test the parameterization. (c) The O( volume mixing ratio profiles. Here and in the following figures, the geometric altitude is approximate, corresponding to the pressure–altitude relationship of the MLE reference atmosphere.
[Figure omitted. See PDF]
The parameterization computes cooling rates for given inputs of temperature and concentrations of CO, O(, O and N as a function of pressure. No specific grid is required, and it can be irregular. The routine interpolates the given parameters to its internal pressure grid. Possible cooling effects caused by temperature disturbances at vertical scales smaller than the internal grid of the parameterization, km
Further, the collisional (de)activation of CO) levels by the main atmospheric molecules (N and O) and by O( can also be prescribed. To compute the different coefficients employed by the parameterization (see Sect. ), reference LTE and non-LTE cooling rates are required (see Sect. and ). These are calculated for selected reference atmospheres and are described in the next section.
3 The reference atmospheres3.1 Temperature
We used the same six pressure–temperature reference atmospheres as in for altitudes below km. Above this altitude, they were extended up to km with the empirical US Naval Research Laboratory Mass Spectrometer Incoherent Scatter Radar version 2.0 (MSIS2) model for medium conditions of solar activity: sfu (June 2011) for all atmospheres except for MLE, which was sfu (September 2011). These six – profiles cover the envelope of the climatological zonal mean temperatures of the current middle atmosphere very well, e.g. as measured by MIPAS from 2007 to 2012 (see Sect. ). However, they do not cover the small-scale temporal and spatial temperature variability (see Fig. ). The performance of the parameterization for such variability is addressed in Sect. . Further, the range of the six temperature profiles does not cover the episodes of stratospheric warming with an elevated stratopause well. During these events, the altitude region of the typical stratopause, at about 50 km, is much colder, being even 50 K colder than during normal conditions, and the altitude of the typical mesopause, near 85–90 km, is warmer by a similar amount (see Fig. a, profile ES). For these conditions the temperature profile is nearly isothermal from the tropopause up to about 0.1 hPa and exhibits an inversion above, with a peak near the mesopause. We anticipate that, for these rare conditions, the error incurred by this parameterization can be significant (see Sect. ).
We should also mention that the envelope of these reference atmospheres does not fully cover the predicted temperatures for the end of this century for projections with high CO emissions. In particular, WACCM simulations for this century under the RCP6.0 scenario yield zonal mean temperatures which are colder in the middle atmosphere. In order to cover such predictions, the envelope of the six – profiles assumed here would have to be widened by about K in the upper stratosphere and by about K in the mesosphere. The parameterization accuracy for such predicted temperatures has not been fully assessed in this work as we have considered only the projection of high CO VMR profiles but not the corresponding predicted temperature fields. This will be the subject of future work.
3.2
CO, O(, O and N abundances
The valid range of the parameterization of with respect to CO volume mixing ratios is exceeded by the CO concentration of several scenarios considered in the 6th Coupled Climate Model Intercomparison Project (CMIP6), in particular for the 4CO experiment. Several CO scenarios have been proposed for the future. Thus, proposed the scenario RCP2.6, which reaches tropospheric CO values near 1000 ppmv by the end of the century. Likewise, suggested the high CO scenario of RCP8.5 (CMIP5), which has CO concentrations of 2000 ppmv in the second half of the 23rd century or even higher than 2000 ppmv (e.g. SSP5-8.5 in CMIP6; see ). Here we used a wide range of tropospheric CO values ranging from about half of the pre-industrial (1851) value of 285 ppmv to about 10 times this value (see Fig. b). The specific profiles were built from a WACCM run under the CMIP6 SSP5-8.5 scenario . Global annual mean profiles of CO were taken from WACCM simulations for years: 1851 (pre-industrial), CO profile no. 2; 2014, CO profile no. 3; 2050, CO profile no. 4 (pre-industrial); and 2099, CO profile no. 6 (pre-industrial). In addition, we set up the low CO profile (no. 1) by halving pre-industrial profile no. 2, the intermediate CO profile no. 5 (pre-industrial) from the mean of WACCM outputs for 2050 and 2099, the high CO profile no. 7 (pre-industrial) by multiplying WACCM output for 2099 by a factor of 1.25, and the highest CO profile no. 8 ( pre-industrial) by multiplying WACCM output for 2099 by a factor of 2.7. In addition to those CO VMR profiles, we also composed the intermediate profile nos. 9, 10 and 11 for testing the parameterization (see Sect. ), which are shown in Fig. b with dashed lines. Profile nos. 9 and 11 were obtained by multiplying WACCM outputs for the years 2050 and 2099 by factors of 0.979 and 1.8, respectively. Profile no. 10 was calculated by weighting the WACCM annual mean for the years 2050 and 2099 by 0.76875 and 0.256250, respectively. WACCM provides the CO VMR profiles up to about 130 km. Above that altitude, they were calculated by using a WACCM-X run for 2008, which provides a CO VMR up to near 500 km, and scaling them, in pressure levels, by the CO value of the corresponding CO profile at a pressure of hPa.
As discussed above, the parameterization requires the N, O and O( volume mixing ratio profiles for the six – reference atmospheres. They were taken from the MSIS2 model and are shown in Fig. S2.
4 Cooling rates for the reference atmospheresWe describe in this section the non-LTE cooling rates used as a reference. To compute the coefficients of the parameterization and the boundaries of the different layers, they also require the calculations of the cooling rates in LTE, which are also described in this section. Further, we have assessed the accuracy of the LTE cooling rates by comparing them with those calculated by an independent code, the Reference Forward Model
Figure 2
LTE cooling rates for the US standard temperature profile and the CO VMR of computed by the GRANADA algorithm and the RFM code.
[Figure omitted. See PDF]
4.1 Reference LTE cooling ratesThe LTE cooling rates have been computed using a modified Curtis matrix formulation . In that computation we used as the basis for the radiative transfer calculations (e.g. the optical depths, the transmittances and their differences) the Karlsruhe Optimised and Precise Radiative Transfer Algorithm
In order to ensure the accuracy of these LTE cooling rates, we have compared them with those obtained with another very well-tested and widely used radiative transfer code, RFM . This code has been used in many studies relevant to the MIPAS instrument and for the retrieval of MIPAS level-2 data obtained by the University of Oxford. It worth mentioning that RFM uses a classical Curtis matrix method (double-flux transmittance differences), while we use the modified Curtis matrix method. Figure shows the results of the comparison for the US standard temperature profile and the CO VMR of . Here we used a common fine-altitude grid of 0.5 km. We see that the agreement between both codes is very good, with differences at most of the altitudes smaller than 0.1–0.2 K d. Note that some of the major differences appear to be associated with small oscillations in the RFM results.
The same formulation has been used to calculate the Curtis matrices of all the CO bands which are required to compute the coefficients of the parameterization in the LTE region (see Sect. ).
4.2 Reference non-LTE cooling rates
The reference line-by-line non-LTE cooling rates have been computed by using the GRANADA non-LTE code. The details of the method for solving the system of equations for CO are given in . In addition to the solution of the statistical and radiative transfer equations described in that work for the calculation of the non-LTE populations of the CO levels, here, in order to compute accurate non-LTE cooling rates and to account for the overlap between the different CO bands, we included an additional final iteration computing the radiation fields in all the bands by using the lambda iteration method. This algorithm shares the radiative transfer algorithm with KOPRA . Thus, the details of the radiative transfer calculation related to KOPRA, e.g. line shape, spectroscopic data, wavenumber grid, or line mixing, given in LTE Sect. above, also apply to the non-LTE calculations described here.
For this case of non-LTE cooling rates, each ro-vibrational band contributes according to the non-LTE populations of their upper and lower levels. The non-LTE cooling rates calculated here comprise 16 vibrational bands emitting and absorbing in the 15 region, i.e. the fundamental band, three first hot bands, seven second hot bands of the major CO isotopologue and fundamental bands of isotopologues OCO, OCO, OCO, OCO and OCO. The contributions of other weaker bands arising from higher levels, e.g. , 5 or 6, are included in the calculation but have negligible contributions for the conditions of Earth's atmosphere.
Table 1
The main collisional processes affecting the CO vibrational levels included in the calculations of the non-LTE cooling rates
No. | Process | Rate coefficient |
---|---|---|
1a | CO N CO N | |
1b | CO O CO + O | |
(, 3, 4.5, 6.75, 9, 12 for , …, 6) | ||
1c | CO + O CO + O | |
2a | CO+N CO(, 3 or 4,) N | + |
2b | CO O CO(, 3 or 4,) O | |
2c | CO O CO(+2, 3 or 4,) O | |
3 | CO(,) CO CO(+1,) CO() | |
4 | CO(,) O CO(+1,) O(1) | |
5a | CO+CO CO) CO() | See Table 6 in . |
5b | CO+CO CO) CO+1) | ; –4; |
6 | CO N CO) N(1) | for , 1, 3, 5, 7 |
for , 4, 6 | ||
7 | CO N CO N | for or & |
for all other cases |
Rate coefficient for the forward sense of the process (cm s). This rate is taken as 10 cm s for temperatures lower than 150 K
For the calculations of the non-LTE cooling rates, a collisional scheme and collisional rates are required. Although the collisional rates affecting the CO levels are an input parameter for the parameterization, here we have used, for the calculations of the reference cooling rates and for testing the parameterization, the collisional rates described in . They have been recently revised and used in the non-LTE retrieval of temperature from SABER and MIPAS measurements . The most relevant collisional processes concerning the populations of the levels emitting in the different bands described above, and their rates, are listed in Table for easier reference. We should note that these rates and their temperature dependencies are different from those used in the previous parameterization of . The values are in general of very similar magnitudes, except for the rate (process 1c in Table ) that has been considered here with its upper limit, i.e. about a factor of 2 larger than in the parameterization of . This rate coefficient is not well known, with uncertainties of the order of a factor of 2
The cooling rates near 15 change very little with the illumination conditions. However, those cooling rates (or more strictly speaking, the flux divergence) of the CO bands computed by GRANADA under daytime conditions are affected by some emission from the relaxation and/or redistribution of the solar energy absorbed in the near-infrared bands
Figure 3
The non-LTE cooling rates for four reference atmospheres shown for altitudes up to the lower thermosphere. The cooling rates extended to the thermosphere are shown in Fig. S5. Note the different scales.
[Figure omitted. See PDF]
Figure 4
Contributions of the different CO bands to the cooling rates for four – reference atmospheres and the contemporary CO VMR no. 3 profile. Note the different scales.
[Figure omitted. See PDF]
The results for the accurate, line-by-line non-LTE cooling rates computed for the six reference – profiles and the eight CO VMRs are shown in Fig. from the stratosphere up to the lower thermosphere and in Fig. S5 for the upper part of the parameterization, i.e. from 80 to 200 km. We observe that the altitude distribution of the cooling rates depends very much on the temperature profile. This is the major difficulty in building the parameterization. A general common feature is the maximum near the stratopause, because at these altitudes the non-LTE cooling rates do not differ significantly from those in LTE, and these are mainly driven by the high temperature of this region. Above the stratopause, the total non-LTE cooling rates depend very much on the contributions of the different bands, e.g. the fundamental band of the major isotopologue (FB), the contributions of the first and hot bands (Hots) and those of the fundamental bands of the five minor isotopologues. These contributions are shown in Fig. for the contemporary CO VMR profile (no. 3) and four – profiles. The non-LTE cooling rates generally decrease with altitude above the stratopause, reaching a minimum near the mesopause for several – profiles; see e.g. the SAS and MLE atmospheres in Figs. and . The cooling can even be negative, e.g. heating for the very cold sub-arctic summer (SAS) mesopause, where heating can be several Kelvin per day (see the bottom-left panel of Fig. ). Exceptional cases are the winter atmospheres (mid-latitude winter, MLW (not shown), and sub-arctic winter, SAW), where the mesopause is warmer and the cooling rates are high in this region. Above the mesopause, the cooling rate rapidly increases following the enhancement of the kinetic temperature. Above about 130 km, however, the cooling rates decline (see Fig. ). At these altitudes, cooling to space is a very good approximation to the non-LTE cooling rate, which, when expressed in Kelvin per day, is proportional to the CO VMR, to the atomic oxygen density [O] and to the temperature through
Figure 5
Non-LTE–LTE cooling rate differences for four – reference atmospheres and the contemporary CO VMR no. 3 profile. Differences for the higher CO VMR no. 6 profile ( the pre-industrial profile) are shown in Fig. S6. The “*” symbol indicates the pressure level (hPa) where the non-LTE–LTE difference reaches 5 %. Note the different scales.
[Figure omitted. See PDF]
The dependence of the non-LTE cooling rates on the CO abundances is illustrated in Fig. . We observe that, in general, the cooling rate correlates very well with the CO abundance, although that correlation is not always linear and generally depends on altitude. This is also true for the cases where we have net heating for low CO VMR, e.g. the SAS atmosphere between about 75 km and 95 km. For the MLE and MLS atmospheres, the cooling rate near km changes from net cooling to net heating for the largest CO VMRs. Further, the very low cooling for the tropical (TRO) – profile near 70 km remains very low even when the CO VMR varies by a factor of 20. In the lower thermosphere, e.g. above around km, however, the dependency of the cooling rate on the CO is very close to being linear (see Fig. S5).
Figure 6
Effect of the collisional rate (or, equivalently, the O( concentration) on the non-LTE cooling rates for four – reference atmospheres and the contemporary CO VMR no. 3 profile. Note the different scales.
[Figure omitted. See PDF]
A comparison of the non-LTE and LTE cooling rates for the six – reference atmospheres for CO VMR nos. 3 (current VMR) and 6 (4 times the pre-industrial value) is shown in Figs. and S6. This comparison is necessary in order to establish the boundaries of the different atmospheric regions of the parameterization (see Sect. ). We first observe (Fig. ) that the altitude of the departure of the cooling rate from LTE (considered the altitude where the non-LTE–LTE difference is larger than 5 %) depends on the temperature profile and ranges from hPa ( km) for the SAS atmosphere to hPa ( km) for the TRO atmosphere. A similar plot but for the higher CO profile no. 6 ( pre-industrial) is shown in Fig. S6. An overview of the altitude or pressure level of the deviation of the cooling rate from LTE is shown in Fig. for the six – profiles and the eight CO VMR profiles. We see that the lower altitudes (higher pressures) occur for the SAS and SAW reference atmospheres. It is also evident that this altitude increases with the CO VMR, except for the SAS and SAW cases, for which it is nearly independent of the CO VMR. That is expected as, for a more abundant CO atmosphere, the 15 bands become optically thicker and fewer collisions are sufficient for keeping the emitting levels in LTE. Figure suggests that, for higher CO VMRs, the LTE–non-LTE transition region could be placed at higher altitudes. However, as the parameterization is intended to cover the full range of CO VMR profiles, we have to be conservative, and we placed it at the lowest altitude for any – or CO VMR profile. Thus, it has been taken at ( hPa, km), which, except for the SAW atmosphere with the lowest CO profile, fulfils the LTE–non-LTE transition region for all the – and CO VMR profiles.
For completeness, Fig. shows an example of the comparison of non-LTE and LTE cooling rates, including the thermosphere for the six – profiles and the current CO values. This shows the enormous difference between LTE and non-LTE cooling rates (non-LTE values being much smaller) in the thermosphere.
The atomic oxygen concentration is an input to the parameterization and plays a crucial role in determining the CO infrared cooling. As a consequence, it is very important in establishing the different non-LTE regions of the parameterization . To identify the atmospheric regions where it is important, we have performed a calculation by dividing the collisional rate by a factor of 2, which is almost equivalent to changing the O( concentration by the same factor. Figure shows this effect for four – profiles and considering the current CO. Generally, it is most important above around 10 hPa ( km). However, for the SAS and SAW atmospheres, it is also important down to hPa ( km). The fact that its importance starts being significant at different atmospheric levels for the different – profiles poses an additional difficulty in the development of the parameterization.
Figure 7
Atmospheric regions considered in the parameterization. represents the boundaries of the layers: , , and . The temperature profiles used in the reference calculations are also shown as a reference.
[Figure omitted. See PDF]
5 The parameterizationEssentially, here we follow the parameterization developed by . A brief description of the method, including the most important features and equations, is given in this section.
The atmosphere is divided into five different regions (see Fig. ) where different approaches are used for calculating the cooling rates. These regions are qualitatively the same as those defined by , but their altitude extensions (except for the LTE region) have been significantly revised, mainly as a consequence of the ample range of CO abundances for which this parameterization is developed. In fact, their upper boundaries have been moved upwards (except for LTE), resulting in the following ranges: LTE –9.875 (– km), NLTE1 –12.625 (–87 km), NLTE2 –16.375 (–109 km), NLTE3 –19.875 (–180 km) and NLTE4 ( km).
The lowermost (LTE) and uppermost (NLTE4) regions are the most straightforward and also the regions where the errors are in general smaller. The most difficult parts are the transition regions from LTE to non-LTE, where (i) several bands contribute to the cooling with different source functions and their relative contributions depend very much on the actual temperature profiles (see Fig. ), and (ii) the exchange of radiation between the layers is significant and different for the considered bands. Further, although most of the radiative excitation at a given layer is produced by the absorption of photons travelling from below, the absorption of photons travelling downwards can also contribute significantly. This is the case, for example, for the stronger fundamental band near the mesopause. Further, the cooling above around 90 km also depends on the collisions with atomic oxygen. This effect can be accurately taken into account in the upper non-LTE regions where all the bands become optically thin. However, it is very difficult to represent it properly between around the mesopause and a few tens of kilometres above, where the atomic oxygen concentration varies largely and the exchange of radiation between the layers is still important.
5.1 The LTE region
The parameterization in the LTE region is based on the Curtis matrix method. The cooling rate at a given pressure level , in a spectral region and for a particular temperature profile is given by
1 where the indices refer to pressure levels and , and the sum is extended over pressure levels ranging from the lower boundary, , until . This upper boundary of the Curtis matrix has been selected to minimize the error in the lowest non-LTE region, NLTE1 (see more details in Sect. below). is the modified Curtis matrix (slightly different from its usual definition; see e.g. ). The factor in Eq. () represents the exponential part of the Planck function and is given by 2 where is the Planck constant, is the Boltzmann constant, and is the temperature of the – profile at level . Similarly, corresponds to the exponential part of the Planck function for the surface temperature . The term accounts for the absorption at level times the transmission from the surface up to level at , so that accounts for the heating rate due to the absorption of the radiation from the surface (or lower boundary). The cooling rate is calculated in the spectral range from 540 to 800 cm divided into frequency intervals, , 10 cm wide. Those LTE cooling rate profiles have been calculated for each of the six – reference atmospheres and the eight CO VMR profiles.
In the parameterization, the Curtis matrix is expressed with an explicit temperature dependence by where the matrix coefficients and are given by and , and are the band strengths of the fundamental, first hot and second hot bands, respectively, in the interval. In this way, the temperature dependence, mainly caused by the band strength of the first and second hot bands, is carried out in . Those matrix coefficients are calculated for each spectral interval . We obtain the coefficients for the entire spectral region of the CO 15 bands by summing over all the intervals and weighting with the dependency of the factor, e.g. with cm being the frequency of the fundamental band of the major isotopologue.
Next, we define global and matrix coefficients, to be used for any input temperature profile, as weighted averages of and for the six reference – profiles. We introduce a set of normalized weights , altitude-dependent, for each temperature profile so that 3 Analogously to the matrix coefficients and , we define the corresponding vector coefficients for the surface flux, and , so that with 4 In this way, the cooling rate , at a pressure level , for a given input temperature profile is calculated in the parameterization by 5 The weights are obtained by minimizing the cost function at each pressure level , given by the sum of the square of the differences of the reference LTE cooling rates, , and those computed by the parameterization by using Eq. (), , for each – profile , e.g. or, in more detail, by The normalized coefficients were originally introduced to consider different fractions of the area of Earth ascribed to each – reference profile. Thus, in the previous parameterization they were taken to be equal to 0.05 for the subarctic (winter and summer) profiles, 0.1 for the mid-latitude (winter and summer) profiles, 0.4 for the tropical profile and 0.3 for the mid-latitude equinox – profile. In this study, we have explored different options including the original coefficients and a uniform weighting for the six – profiles, and we found a smaller for the latter, e.g. for all the profiles. Hence, that was included in this version.
In that way, we have parameterized the cooling rates as a function of temperature. The cooling rates also depend on the CO VMR profiles (see Fig. ). The parameterization incorporates the dependence on the CO abundance by calculating and for a generic CO profile by assuming a linear interpolation in and from the adjacent CO VMR profiles. Thus, the and coefficients of Eq. () have been calculated (and are provided) for the eight CO VMRs shown in Fig. b.
5.2 The NLTE1 region: the transition from LTE to non-LTEThis region is difficult to parameterize because we have several bands contributing to the cooling (see Fig. ), and their relative contributions depend significantly on both the temperature structure and the CO VMR profile. Note that, at certain levels, the cooling induced by the weaker hot bands is larger than that of the stronger fundamental band. We should also note that the contribution of the first hot bands at high altitudes, –150 km, is not negligible (5 %–10 %, Fig. ). This contribution is accounted for in the parameterization by implicitly assuming that it is produced by the fundamental band of the main isotopologue (see below and Sect. ).
The lower boundary of this region, i.e. the LTE–NLTE1 transition, occurs, depending on the temperature profile, at altitudes from km up to km (0.08 to 0.004 hPa; see Fig. ), taking place a few kilometres lower for the subarctic summer atmosphere and the lowest CO VMR. This transition region occurs at higher altitudes for larger CO VMRs. That is, the atmosphere becomes optically thicker and fewer collisions are enough to keep the levels in LTE. However, since we also need to represent the low CO VMRs, we decided to conservatively set up this region at rather low altitudes, ( km), the same level used in the previous parameterization.
The upper limit of this region was set up in the previous parameterization at the pressure levels where collisions with O( start affecting the cooling rates significantly. Again, that pressure level depends on the temperature profile (and also on the O( concentration), being lower at hPa (), for the subarctic summer and winter conditions (see Fig. ). Here, we have taken the upper boundary of ( km), slightly higher than the 12.5 value assumed in the original parameterization.
In this region we followed, as in , the matrix approach discussed in Sect. above. Thus, Eq. () was used but with modified and coefficients that account for the non-LTE corrections. For each – profile, , we define and where and are the reference LTE and non-LTE cooling rates, respectively. Then, the general and coefficients were calculated by following the same procedure as for the LTE region, i.e. by weighting the –-specific and coefficients with a set of altitude-dependent weights and minimizing the total cost function (see Sect. ). In this way, we obtain and and the cooling rates are computed by using Eq. () but replacing and by and , respectively, i.e.
6 This procedure, while producing a perfect match for a single atmosphere by construction, generates irregularities for other atmospheres at some levels, e.g. when using for atmosphere at points where is close to zero. We observed that the irregularities were significantly mitigated by reducing the dimensions of the Curtis matrix from to , where corresponds to ( hPa, km), i.e. by placing slightly above the boundary between the NLTE1 and NLTE2 regions. Errors induced in the LTE cooling rates by the matrix reduction are negligible (smaller than 0.05 K d at the upper boundary).
5.3 The NLTE2 and NLTE3 regions: the recurrence formula with and without correctionThe parameterization in the NLTE2, NLTE3 and NLTE4 regions is based on the recurrence formula proposed by . This approach is valid when the cooling rate is dominated by the fundamental band, and the absorption of radiation coming from the layers above the layer at work can also be neglected . The boundaries of these regions are then adapted to the applicability of that approach. The NLTE3 boundaries were chosen to embrace the region where those conditions are fulfilled to a large degree. In the layers below, i.e. in the NLTE2 region, that formula is not accurate and requires a correction term to account for the absorption of radiation coming from the layers above and for the cooling of bands other than the fundamental one of the main isotopologue. The recurrence formula is also the basis for the calculation of the cooling rate in the NLTE4 region (see Sect. ), but it is simplified because the exchange of photons within the layers of this region can be neglected. Further, we should emphasize that the dependence of the cooling rate on the CO VMR in these regions is mainly twofold: on the one hand, its direct dependence (see Eq. below) and, on the other hand, the escape function which depends on the CO column above a given layer (see Figs. S7 and S8). We discuss below the boundaries of the NLTE2 and NLTE3 layers and the expressions used for the cooling rates, i.e. the recurrence formulation.
The lower boundary of the NLTE2 region is set up in the layer where the cooling rate obtained by the corrected recurrence formula is more accurate than that given by the non-LTE-corrected Curtis matrix approach (used in NLTE1). This has been set up at ( km), which is very similar to the value in the original parameterization of 12.5 ( km). The upper boundary of the NLTE2 region is set up at the pressure level where the recurrence formula does not need to be corrected to yield an accurate estimation of the cooling rate. In this work, it has been set up at ( km), which is significantly higher than the value of 14 ( km) set up in the previous parameterization. The main reason for this is that, for the higher CO VMRs used here, the atmosphere becomes optically thicker; hence, the absorption of radiation of the layer above needs to be taken into account, also at lower pressures.
Thus, the cooling rates in the NLTE2 and NLTE3 regions are calculated by
7 where is a constant that depends on the Einstein coefficient of the fundamental band (), on and on the units of
Note that this constant has been changed from its value of in to the value used here of .
. VMR() is the CO VMR; is the mean molecular weight, , where , , and are the collisional rate constants with N, O and O( (see Table ); and [N], [O] and [O(] are the concentrations of the respective species. Note that the collisional rates depend on through their temperature dependencies.The at level , , is obtained by the recurrence formula 8 starting from the lower boundary at , where, using Eq. (), 9 and is obtained by Eq. (). The coefficients above are given by 10 where 11 is the escape function which mainly depends on the CO column, , above a given level . The temperature of those layers affects this function as it influences the line shape of the CO lines and hence the probability of photons escaping into space. This is reflected in Fig. S7a, which shows as a function of the CO column for the six – reference atmospheres and a single CO VMR profile (no. 3, current VMR). The dependence of with the CO VMR profiles is shown in Fig. S7. In our calculations, we have used for the average of this function for the six – reference atmospheres.
The parameter entering Eq. () and needed in the NLTE2 region has been computed by minimizing the following cost function at each point : After performing some sensitivity tests, we used uniform weighting for the different reference atmospheres ( for all atmospheres) rather than the area weighting used in the previous parameterization. Other tests were performed to determine the optimal upper boundary for the correction: extending the region upwards reduces the error in the –19 region but results in a spurious jump at the uppermost boundary, which is avoided when using a lower of 16.375. It is worth noting that above takes values below unity, thus decreasing the escape in the region. For the fit of the optimal , the parameterized value of is considered a starting point rather than the reference value.
In the NLTE3 region, we used the same method as in region NLTE2, except that no correction for the function is applied, i.e. .
5.4 The NLTE4 regionThe recurrence formula described above is also valid in the uppermost NLTE4 region, but, as the CO bands are so optically thin here, the exchange of radiation within the layers of this region can be neglected, and the recurrence formula is reduced to a simpler expression, i.e. a cooling-to-space term and an additional term that accounts for the absorption of the radiation emitted by the layers below its boundary. Thus, the cooling rate for this region is computed by using Eq. () but with a simple expression for , , that gives a smooth transition to the cooling-to-space approximation, e.g.
12 where is obtained from the boundary condition 13 and uses the recurrence formula in Eq. ().
6 Testing the parameterization for the reference atmospheresThe parameterization has been tested against the reference cooling rates calculated for the reference atmospheres (the six – profiles and the eight CO VMR profiles) (see the next section) and for intermediate CO VMRs and the collisional rate (Sect. ). Further, it has been verified for measured temperature profiles that exhibit a large variability (Sect. ) and for the temperature profiles obtained by a high-resolution version of WACCM-X (Sect. ).
Figure 8
Comparison of the cooling rates of the current and previous parameterizations with respect to reference non-LTE cooling rates for the present-day CO VMR profile no. 3 and the six – reference atmospheres. Note that the latter cooling rates are hardly visible in the left panels of the figures. See Figs. S9 and S10 for lower CO concentrations and Figs. S11–S15 for higher CO VMRs.
[Figure omitted. See PDF]
6.1 Accuracy of the parameterization for the reference atmospheresIn this section, we discuss the accuracy of the current parameterization for the assumed reference atmospheres. The non-LTE models used in the original and current parameterizations are different. Hence we expect some differences caused not just by the parameterization itself, but possibly also by the different non-LTE models. Figure shows the cooling rates of this parameterization compared to those of the previous parameterization and the reference ones for a contemporary CO VMR profile (no. 3) and the six – profiles. The comparisons for the lower CO profiles are shown in Figs. S9 and S10 and in Figs. S11–S15 for high CO VMRs. We should clarify that, to make the comparison meaningful, the three sets of cooling rates shown here include the same updated collisional rates (Table ). Note that these rates are different from those used in the previous parameterization . The new parameterization also supports the previous collisional rates, but it has been optimized for the new ones in Table . As expected, larger differences are obtained in the region between 10 hPa ( km) and hPa ( km) and are more marked for the SAS and SAW atmospheres.
Figure 9
Mean of the cooling rate differences of the current and previous parameterizations with respect to the reference cooling rates for the four lowest CO VMR profiles. The shaded areas show the standard deviation of the differences between the six – atmospheres. Detailed comparisons for each – profile are shown in Fig. above and in Figs. S9, S10 and S11 in the Supplement.
[Figure omitted. See PDF]
The differences are more clearly illustrated in Figs. and S16, where we show the mean and standard deviation of the differences for the four lowest and four highest CO VMR profiles, respectively. The improvement of the new parameterization is noticeable (compare the blue and red lines). In general, the cooling rates of the current parameterization are more accurate than in the previous one for most of the regions and temperature structures. We observe that the errors (e.g. the differences with respect to the reference non-LTE cooling rates) of the new parameterization (red curves) are very small overall. They are below K d for the current and lower CO abundances (see Fig. ). For higher CO concentrations, between about 2 and 3 times the pre-industrial values, the largest errors are –2 K d and are located near 110–120 km (see Fig. and the top-left panel in Fig. S16). The quoted values refer to the mean of the differences, although they are larger for the individual – atmospheres. The spread of these values is larger in the region of 10 hPa ( km) to 10 hPa ( km), where the root mean square (rms) reaches values between and K d (Fig. ).
For the very high CO concentrations (4, 5 and 10 times the pre-industrial abundances), the errors are also very small, below K d for most regions and conditions, except in the 107–135 km region, where we found maximum positive biases of , and K d for , and the pre-industrial CO VMR profiles (see Fig. S16). Those maximum errors are however comparable when expressed in relative terms, all about 1.2 %. The significant rms in the region of –120 km is also notable; clearly, this region is more difficult to parameterize, particularly for such a large range of CO abundances.
That increase in the differences of the new parameterization with respect to the reference calculations for the very high CO VMRs near 110 km seems to be related to the transition region from NLTE2 to NLTE3 (see Fig. ). It looks like the cooling in the lower part of the NLTE3 region also requires correction by the factor for high CO VMRs. This suggests that, for higher CO VMRs, the parameterization would be more accurate if this transition altitude has risen. Such a rise, however, would worsen the cooling below this boundary. This manifests the difficulty in obtaining very accurate cooling rates for a large range of CO VMRs with this method.
6.2Assessment of the cooling rates for intermediate CO VMRs and for the collisional rate
The aim of the parameterization is to be used for any CO VMR input profile in the range of profile nos. 1 and 8 of Fig. b and any plausible value for the rates discussed in Sect. . In this section we demonstrate that the parameterization is also very accurate for CO VMRs that fall between the reference profiles used for its development and also when using different values. In particular, we show results for the intermediate CO VMR profile nos. 9, 10 and 11 (see Fig. b) and the collisional rate used in the reference calculations divided by a factor of 2.
Figure 10
Mean of the cooling rate differences of the current parameterization with respect to the reference cooling rates for the intermediate CO VMR profiles (nos. 9–11). The shaded areas show the difference spread (standard deviation) between the six – atmospheres. Note the different scales of the axis. Detailed comparisons for each – profile are shown in Figs. , S17 and S18.
[Figure omitted. See PDF]
Figure shows the results of the calculation for the intermediate CO VMR profile no. 9, which is between the current CO VMR value and that projected for 2050 (2 times the pre-industrial value). We can observe similar features in the calculations for the contemporary CO VMR profile (no. 3) (see Fig. ), although the differences are slightly larger because the CO VMR profile is larger. The distinctions are more clearly seen in Fig. , where we show the mean and standard deviation of the differences for the six – profiles. The patterns in the differences, as well as their values and spreads, are very similar to those described above in Sect. for the CO reference profiles. The major differences appear between 105 and 135 km, reaching maximum values of 1, 2 and 9 K d for VMR profile nos. 9, 10 and 11, respectively. Again, we observe that the new parameterization is more accurate at practically all altitude levels. Further, the maximum values of the standard deviations of the differences for the various – profiles also resemble very much those discussed before, reaching maximum values of about 2 K d, 3 K d and 15 K d for the respective CO VMR profiles. Note that, although these values are larger for higher CO VMRs, they are very similar when expressed in percentage.
As the CO()–O( collisional rate, , is still uncertain nowadays by about a factor of 2
Figure 11
Mean of the cooling rate differences of the current parameterization with respect to the reference cooling rates for the results in Fig. . The shaded areas show the differences spread (standard deviation) between the six – atmospheres.
[Figure omitted. See PDF]
6.3 Performance of the parameterizationIn Table we list some examples of the time taken to execute the parameterization for two processors, two compilers and three atmospheric intervals. Better performance is noticeable (up to a factor of 5 faster) when using the Intel Fortran Compiler Classic (
Table 2
Performance of the parameterization.
Processor | Compiler | Input atmosphere | Time (ms) |
---|---|---|---|
type, specs | (no. of levels, altitude range) | ||
1, x86_64 64-bit AMD EPYC 7742 2.245 GHz | 273, 0–550 km | 0.074 | |
1, x86_64 64-bit AMD EPYC 7742 2.245 GHz | 185, 50–550 km | 0.042 | |
1, x86_64 64-bit AMD EPYC 7742 2.245 GHz | 162, 50–270 km | 0.041 | |
1, x86_64 64-bit AMD EPYC 7742 2.245 GHz | 273, 0–550 km | 0.360 | |
1, x86_64 64-bit AMD EPYC 7742 2.245 GHz | 185, 50–550 km | 0.204 | |
1, x86_64 64-bit AMD EPYC 7742 2.245 GHz | 162, 50–270 km | 0.197 | |
2, x86_64 64-bit Intel Core i7 4.2 GHz | 273, 0–550 km | 0.133 | |
2, x86_64 64-bit Intel Core i7 4.2 GHz | 185, 50–550 km | 0.077 | |
2, x86_64 64-bit Intel Core i7 4.2 GHz | 162, 50–270 km | 0.073 |
Figure 12
Zonal mean of the differences in the cooling rates of the old parameterization (a, c) and the new one (b, d) with respect to the reference cooling rates obtained for MIPAS temperatures for 14 January 2010 (solstice) and 25 March 2010 (equinox). Similar figures for 13 June 2010 (northern summer hemisphere) and 21 September 2010 (northern autumn hemisphere) are shown in Fig. .
[Figure omitted. See PDF]
7 Testing the parameterization for the MIPAS-measured temperatures7.1 Solstice and equinox conditions
We have compared the cooling rates estimated by the parameterization with the reference ones for realistic, e.g. measured, temperature profiles that present a large variability and very variable vertical structure (see e.g. Fig. ). Specifically, we compared them for the – profiles measured by the MIPAS instrument for 5 full days of measurements (about 2500 profiles) with global latitude coverage and covering 2 d for solstice conditions (14 January and 13 June) and 2 d for equinox conditions (25 March and 21 September) for 2010. Further, we compare the results for the temperatures of 15 February 2009 when strong stratospheric warming followed by an elevated stratopause event occurred in the northern polar hemisphere (see Sect. ). The comparison is carried out for the MIPAS measurements taken only under nighttime conditions, as the MIPAS non-LTE cooling rates for daytime, obtained simultaneously with the temperature inversion, also include the fraction of the 15 cooling which is produced by the relaxation of the solar energy absorbed by CO NIR bands, which is not accounted for in this parameterization (see Sect. ). The zonal means of the temperatures, CO VMRs and O( abundances for those conditions are shown in Figs. , S19 and S20, respectively.
Figure 13
Mean (solid lines) and rms (dash lines) of the differences in the cooling rates of the old parameterization (blue) and the new one (red) with respect to the reference accurate cooling rates obtained for MIPAS temperatures for 14 January and 13 June 2010 (solstice, a, c) and 25 March and 21 September 2010 (equinox, b, d).
[Figure omitted. See PDF]
The results are presented in Fig. for the zonal mean of the differences for 1 d of solstice and 1 d of equinox conditions and in Fig. as the global mean difference for all latitudes for each of the 4 individual days. In general, the new parameterization is slightly more accurate. For example, the deviations of the cooling rates from the reference calculations in the altitude range of 105–115 km are larger in the old parameterization (about 2 K d) than in the new one and are negligible in this region. Also, the differences with respect to the reference calculations are larger in the altitude range of 80–95 km for solstice conditions and at altitudes of 80–100 km for equinox conditions (see Fig. ).
Overall, the errors in the mean profiles of the cooling rates of the new parameterization for 1 d of measurements are below 0.5 K d, except in the region between and hPa (–95 km), where they can reach values of 1–2 K d. This region is the most difficult to parameterize because several bands contribute to the cooling rate, and they are very sensitive to the temperature structure of the middle atmosphere (e.g. even outside this region). Note also that this is precisely the region where the rms of the differences of the cooling rate with respect to the reference ones are largest, reaching values of up to 6 K d (see Fig. ).
Figure 14
As in Fig. but for the MIPAS temperatures taken on 15 February 2009, when a major elevated stratopause event occurred.
[Figure omitted. See PDF]
Figure 15
Mean (solid) and rms (dash) of the differences in the cooling rates of the old parameterization (in blue) and the new one (red) with respect to the reference cooling rates obtained for MIPAS temperatures for 15 February 2009 for latitudes north of 50° N.
[Figure omitted. See PDF]
7.2 Elevated stratopause conditionsThe comparison of the cooling rates estimated by the old and new parameterizations with respect to the reference calculations for 15 February 2009, a day with a pronounced and unusual elevated stratopause event (see the zonal mean temperatures in Fig. ), is shown in Fig. . Similar features to the other conditions shown above can be appreciated, except in the polar winter region. The mean of the differences and the standard deviations for all the profiles at latitudes north of 50° N are shown in Fig. . The differences are significantly larger than for other latitudes in the 80–95 km altitude region. Both parameterizations underestimate the cooling in that atmospheric region. The new parameterization has, however, better performance above about 80 km, but in the elevated stratopause region (80–100 km), it still underestimates the cooling by 3–7 K d (10 %).
It seems clear that part of this underestimation is caused by the fact that such atypical temperature profiles (see Sect. ) were not considered in the parameterization. However, its inclusion would not solve the problem, as in the calculations of the coefficients a trade-off of the weighting of the different – reference atmospheres has to be chosen (see Sect. and ). Thus, it might ameliorate the inaccuracy for these elevated stratopause events but would worsen the accuracy for other general situations. This manifests the difficulty or limitation of this method in providing accurate non-LTE cooling rates for all the temperature structures (gradients) that we might find in the real atmosphere. Nevertheless, we have to keep in mind that these situations are sporadic and limited to high polar regions. Hence they should not impact significantly the accuracy of the cooling rates of this parameterization in global multiyear GCM simulations.
8 Testing the parameterization for WACCM-X high-resolution temperatures
In addition to the tests above, we have also tested the parameterizations for the temperature structure obtained with a high-resolution version of the WACCM-X model . This version of WACCM-X has a fine grid of 0.25° 0.25° in latitude longitude and a vertical resolution of 0.1 scale heights ( km) in most of the middle and upper atmosphere, transitioning to 0.25 scale heights in the top three scale heights
We recall that the grid of the parameterization is 0.25.
. With such a fine grid, the model itself can internally generate gravity waves, thus providing temperature profiles with a vertical structure very similar to that measured by high-vertical-resolution lidars of the mesosphere and lower thermosphere. Some examples of – profiles of the model exhibiting those vertical features and the latitudinal and longitudinal variabilities are shown in Fig. S23 of the Supplement. Further, the model spans from the surface up to nearly 600 km ( hPa), which is ideal for testing the parameterization. In addition to the pressure–temperature profiles of the model, their O(, O and N VMR profiles have been used. A contemporary CO VMR (profile no. 3) was included in these calculations.The parameterization has been tested for a total of 225 temperature profiles. They have been selected from the model output for January conditions at four latitudes, 20, 40, 60 and 70° N of the northern winter hemisphere, and two additional latitudes, 60 and 70° S of the southern summer hemisphere. For each latitude, 36 profiles corresponding to longitudes from 0 to 360° every 10° were selected. A few – profiles are shown in the left column of Fig. , and all the profiles for latitudes 20° N, 60° N, 70° N and 70° S are shown in Fig. S23 in the Supplement. We should note that those temperature structures very much resemble those measured by lidar instruments.
Figure 16
Comparison of cooling rates of the parameterization and the reference calculations for high-resolution WACCM-X temperature profiles. The left column shows the panels for the pressure–temperature profiles, with two profiles (black and red) in each panel. The middle column shows the reference cooling rates (solid) and the cooling rate from the parameterization (dotted) for the corresponding – profiles. The right column shows the differences of the parameterized cooling and the reference calculation. The different rows are results for the – profiles at different latitudes.
[Figure omitted. See PDF]
Figure 17
Mean and rms of the differences in the cooling rates of the parameterization with respect to the reference cooling rates obtained for the WACCM-X high-resolution temperature profiles.
[Figure omitted. See PDF]
Figure 18
Solar NIR heating rates for the tropical atmosphere and a SZA of 44.5° computed by the solar NIR heating parameterization of and those obtained in the updated GRANADA model (see the text). The right panel shows the differences in the heating rates of the parameterization minus those computed in GRANADA.
[Figure omitted. See PDF]
The results for a few representative – profiles are shown in Fig. . A few more examples are shown in Figs. S24–S26 in the Supplement. In general, the results are very similar to those obtained for the MIPAS temperatures. The parameterization works very well below hPa ( km), with differences generally smaller than 1–2 K d. In the upper regions, above hPa ( km, not fully shown in Fig. because of the small scale chosen to highlight the differences in the region of larger differences), it also works very well. In this region, the cooling rate differences are slightly larger than near hPa ( km), generally below 5 K d, but they are much smaller in relative terms since the cooling rates at high altitudes are very large (of the order of 100–300 K d). In the intermediate region, between hPa and hPa, the algorithm still reproduces the reference calculations rather well but not as well as in the other regions. The cooling rate differences can reach up to 10 K d at a few isolated levels of some individual profiles which can represent up to about 20 %. It is noticeable though that, while the differences can be significant, the profile shapes of the reference and parameterized calculations are very similar (see the middle column of Fig. ).
To have a global perspective, we have plotted in Fig. the mean of the differences for all the – profiles together with their rms. We see that the mean (bias) of the parameterization is very small, practically below 1.5 K d anywhere and below 0.5 K d for most of the atmospheric layers. The rms, a representative error of individual profiles, is also small at levels below hPa ( km) with values of 1–2 K d ( %) and above 10 hPa ( km) with values smaller than 4 K d ( %). In the intermediate region, between and hPa, the rms is however significant, with most values in the range of 5–12 K d. While these values are significant in percentage at – hPa, they are very small above hPa.
9Discussion: the use of this parameterization with a previous CO solar NIR heating rate parameterization
Some of the GCM models use the parameterization of the CO 15 cooling together with that of the CO NIR heating of . Hence, as we are updating the former for larger CO abundances and as the update of the NIR heating parameterization for the large CO abundances is beyond the scope of this work, we investigate whether the latter is still valid for the large CO abundances. For that purpose, we compute CO NIR heating rates with the parameterization of and with the GRANADA model for the large CO concentrations for the six – reference atmospheres. We should note that the non-LTE models used in both parameterizations are different, and hence we expect some difference caused not just by the parameterization itself but also by the differences between the underlying non-LTE model in the NIR heating parameterization and GRANADA. The CO NIR heating rates of GRANADA were calculated with the rate coefficients and photolysis rates described in but updated with those described in and also with those described below. In particular, the rate used in these calculations is % smaller than in below 100 km and thus leads to an [O(] of % smaller below 90 km but that is very similar near 100 km. Above km, the coefficient used in the present calculations is about 40 % smaller than in , leading to a similar reduction in [O(]. Further, we updated the following collisional rates. The rate coefficient of N O( N(1) O has been increased by a factor of 1.08, and the collisional deactivation of N(1) with atomic oxygen (which has an important role in the heating rates; see e.g. ) has been updated from values of to cm s.
The results of the comparison are shown in Fig. for the tropical atmosphere and an intermediate solar zenith angle (SZA) of 44.5°. The region of most importance for the CO NIR heating rates is that comprised between 0.1 and 0.01 hPa . In this region, the differences between the algorithm of and GRANADA are in the range of to K d for CO VMRs up to 5 times the pre-industrial CO profile, e.g. about 10 % to 15 %. Hence, given that they have been computed with very different non-LTE models and the significant effect that parameters like the CO VMR above km, the collisional rate between N(1) and O(, the O( concentration itself and the rate of exchange of CO quanta with N(1) have on this solar heating (see e.g. ), these differences are reasonable. Hence the new CO cooling rate parameterization reported here can be safely used together with the CO solar NIR heating parameterization of for CO VMRs up to 5 times the pre-industrial CO profile.
10 Summary and conclusionsAn improved and extended parameterization of the CO 15 cooling rates of Earth's middle and upper atmosphere has been developed. It essentially follows the same method of the parameterization of . The major novelty is its extended range of CO abundances, ranging from CO profiles with tropospheric values close to half of the pre-industrial value to 10 times that value. This extension of CO profiles can still be safely applied to the parameterization of the CO near-infrared heating of up to at least 5 times the pre-industrial CO values and is normally combined with this cooling rate parameterization.
Other improvements or updates are as follows. They have an extended and finer vertical grid, increasing the number of levels from 8 to 83. The CO line list has been updated, from HITRAN 1992 to HITRAN 2016. Although the collisional rate coefficients affecting the CO and levels are input parameters for the parameterization, in this version we have used more contemporary values, e.g. as currently used in the non-LTE retrieval of temperature from CO 15 emissions of SABER and MIPAS measurements . The rate coefficients are in general of a very similar magnitude, except for the collisional deactivation of CO(,) levels by atomic oxygen, which is now larger by approximately a factor of 2, e.g. close to its accepted upper limit. As a consequence of the larger range of CO VMR profiles, the different NLTE layers for computing the cooling rates have been significantly revised. For example, it is worth mentioning that the lowermost altitude of the cooling-to-space approximation (the uppermost NLTE layer) has risen from km to 160–170 km.
The new parameterization has been thoroughly tested against line-by-line LTE and non-LTE cooling rates for (i) the six – reference atmospheres; (ii) the two most important input parameters (besides temperature), the CO VMR profiles and the collisional rate of CO(,) by atomic oxygen; (iii) realistic measured temperature fields of the middle atmosphere (about 2500 profiles), including an episode of strong stratospheric warming with a very elevated stratopause; and (iv) the temperature profiles (225 profiles) obtained by a high-resolution version of WACCM-X capable of generating internally gravity waves and hence with temperatures showing a large variability and pronounced vertical wave structures. Further, to illustrate the improvements, the comparisons of points (i) to (iii) have also been performed for the previous parameterization.
For the reference temperature profiles, the errors of the new parameterization (mean of the differences in the cooling rates with respect to the reference calculations for the six – atmospheres) are below 0.5 K d for the current and lower CO VMRs. For higher CO concentrations, between about 2 and 3 times the pre-industrial values, the largest errors are –2 K d and are located near 110–120 km. For the very high CO concentrations (from 4 to 10 times the pre-industrial abundances) the errors are also very small, below K d, for most regions and conditions, except in the 107–135 km region, where the parameterization overestimates them in a few Kelvin per day ( %). For these reference atmospheres, the new parameterization has a better performance for most of the atmospheric layers and temperature structures.
From the testing of the parameterization for realistic current temperature fields of the middle atmosphere as measured by MIPAS, we found that, in general, the new parameterization is slightly more accurate. In particular, in the 105–115 km range, the previous parameterization overestimates the cooling rate by 1.5 K d, while the new one is very accurate. However, in the other height regions the difference is not so important. The new parameterization has a better performance in the 80–95 km altitude region. Overall, the errors in the mean profiles (bias) of the cooling rates of the new parameterization, calculated for four different atmospheric conditions with about 500 profiles in each of them, are below 0.5 K d, except between and hPa (–95 km), where they can reach biases of 1–2 K d. That region is the most challenging to parameterize because several CO 15 bands contribute to the cooling rate, and they depend very heavily on the temperature structure of the whole middle atmosphere (e.g. even outside this region). For single-temperature profiles, the cooling rate error (characterized by the rms of the difference between the reference and the parameterized cooling rates) is about 1–2 K d below hPa ( km) and above hPa ( km). In the intermediate region, however, it is significant, between 2 and 7 K d. We have further tested the parameterization against very rare and demanding situations, such as the temperature structures of stratospheric warming events with an elevated stratopause. In these situations, however, the parameterization underestimates the cooling rates by 3–7 K d ( %) at altitudes of 80–100 km, and the individual cooling rates show a significant rms (5–15 K d).
In addition, we have tested the parameterization for the temperature structure obtained by a high-resolution version of WACCM-X, with the temperatures showing a large variability and pronounced vertical wave structure. The mean (bias) error of the parameterization is very small, smaller than 0.5 K d for most atmospheric layers, and below 1.5 K d for almost any altitude from the surface up to 200 km. The rms of the differences in the cooling rates from the parameterization and the reference model is similar to that obtained for MIPAS temperatures, with values of 1–2 K d ( %) below hPa ( km) and smaller than 4 K d ( %) above 10 hPa ( km). In the intermediate region, between and hPa, they are slightly larger than for MIPAS, with values in the range of 5–12 K d. These values, while they are significant in relative terms at – hPa, are very small in percentage above hPa.
As has been shown, this parameterization has some limitations (see Sects. , and ). In order to be able to apply specific approximations for the cooling rates, it has been designed for fixed atmospheric regions where specific radiative transfer regimes prevail. Thus, its extension to a very large range of CO abundances inevitably causes a loss of accuracy for extreme cases in specific atmospheric layers. A possible solution for future updates could be to use different extensions of the non-LTE regions (i.e. Fig. ) for different abundances of CO. Likewise, this parameterization (like the original one) was devised for use in GCMs, i.e. to produce accurate cooling rates globally, e.g. when considering all expected temperature profiles covering the different latitudinal and seasonal conditions. Thus, the ability of the parameterization to compute accurate cooling rates for individual temperature profiles with large temperature gradients in the hPa ( km) to hPa ( km) region is limited. On the contrary, it is extremely fast. The routine takes only 15 of CPU time to calculate a profile in the range of 50 to 270 km on a machine with an Intel Core i7 4.2 GHz processor when compiled with
Appendix A Notes and recommendations for using the parameterization
The routine source code is written in Fortran 90 and is available at 10.5281/zenodo.10849970 . It has been devised for implementation in general circulation models, although it can also be used for other purposes, e.g. to compute the CO 15 cooling rate for a given reference atmosphere.
The code is organized in a library (in the directory source/modules/) that can be included in a more complex GCM model. The subroutine to be called is
The following inputs are required (in order) by
-
Atmospheric profiles as a function of pressure for temperature and four VMRs of CO, O, O and N
-
lev0 : the index of the given pressures so that(lev0) is the maximum pressure level (lower boundary) to be considered for calculating the heating rate. Heating rates will be calculated from that pressure up to the minimum pressure specified in the pressure array. For example, if is given in the range of 10 to 10 hPa (or 10 to 10 hPa) and(lev0) 1 hPa, the heating rate will be calculated in the range of 1 to 10 hPa. -
surf_temp : surface temperature – if set to a negative value, the temperature of the maximum pressure level will be used. -
hr : heating rate. This is an input–output array with the same dimension of pressure. It will only be calculated at pressures in the range ofp(lev0) (the maximum pressure considered) to the minimum specified pressure (minimum(pressure)). Note that throughout this paper we have used the term “cooling rates”, e.g. thehr values with changed signs. -
The values are the temperature (K), pressure (hPa), VMRs (mol mol, not ppm) and heating rate (K d).
-
Input profiles can run either from the ground to the top of the atmosphere (decreasing pressures) or reverse (top to ground with increasing pressures). The pressure grid can be irregular.
-
Important notes: calculations in the LTE region.
-
Pressure levels should include the surface pressure (near 10 hPa), even if the 15 cooling is to be calculated only at lower pressure levels (higher altitudes), i.e.
p(lev0) 10 hPa. -
If 15 cooling is only calculated in the non-LTE regime, it is recommended to set up the lower boundary,
p(lev0) , close to the limit of the LTE–non-LTE transition, e.g. near 1 hPa. In this way, more time-consuming calculations in the LTE region will be avoided.
-
To compile the routine, follow these steps.
-
Edit
Makefile and change the Fortran compiler to your preferred choice (e.g.gfortran orifort ). -
From this folder, run
make .
A test program,
-
First input at line 9:
n_lev ,lev0 andT_surf -
Start from line 12.
-
Six atmospheric profiles are read (
n_lev rows are expected).
The output heating rates are written in the
-
cp input_test.dat input.dat -
./run_cool -
Check that the results in
output.dat are consistent withoutput_test.dat . -
The same procedure can be done for test 2.
-
The rates are defined in the form .
-
for CO–O:
a_zo, b_zo, g_zo (default: , , 76.75) -
for CO–O:
a_zo2, b_zo2, g_zo2 (default: , , 83.8) -
for CO–N:
a_zn2, b_zn2, g_zn2 (default: , , 83.8)
Although the parameterization is specifically developed for the CO 15 non-LTE region, it also works for the LTE region, but the user should be cautious that other important cooling rates in the LTE region, such as those of O and HO, are not included. We recommend that GCM users utilize their radiation scheme in the LTE region and this parameterization in the non-LTE region (e.g. above or 60 km).
Boundaries of the parameterization. About the lower boundary (maximum pressure), see the notes above. About the upper boundary (minimum pressure), there is in principle no limitation, but we recommend setting it as high as the upper lid of your model. There is a large number of GC and CC models with an upper lid at hPa (or km). This parameterization can be used for such models to compute the CO non-LTE cooling rates between and km. We note that, under these circumstances, the cooling rates near the upper lid might not be accurate, as the contribution of the layers above the upper lid is not considered. This, however, is not a limitation of the parameterization itself but an intrinsic limitation of this kind of model. There is no restriction either on the upper limit of the upper boundary, provided it is physically meaningful. That is, it can be placed at altitudes as high as 500 km or higher.
Appendix B Additional figures
Figure B1
Contributions of the different CO bands to the cooling rates for the MLE – atmosphere. That is, like the left panel of Fig. but for the entire altitude range.
[Figure omitted. See PDF]
Figure B2
Altitude (a) and pressure level (b) of the deviation ( %) of the non-LTE cooling rates from LTE values for the six – reference atmospheres and the eight CO VMR profiles (see Fig. b): no. 1, half of the pre-industrial; no. 2, the pre-industrial; no. 3, present day; no. 4, 2pre-industrial; no. 5, pre-industrial; no. 6, pre-industrial; no. 7, pre-industrial; and no. 8, pre-industrial.
[Figure omitted. See PDF]
Figure B3
Non-LTE–LTE cooling rate differences for the six – reference atmospheres and the present-day CO VMR profile no. 3 (as in Fig. ) but including the thermosphere.
[Figure omitted. See PDF]
Figure B4
Comparison of the cooling rates of the current and previous parameterizations with respect to accurate cooling rates for the intermediate CO VMR profile no. 9 (see Fig. b) and the six – reference atmospheres.
[Figure omitted. See PDF]
Figure B5
Testing the effect of the CO()–O( collisional rate on the parameterization. Comparison of the cooling rates of the current and previous parameterizations with respect to the reference cooling rates for the present-day CO VMR abundance (profile no. 3; see Fig. b) and the six – reference atmospheres when reducing the nominal collisional rate in Table by a factor of 2.
[Figure omitted. See PDF]
Figure B6
An example of the MIPAS nighttime temperature profiles (15 February 2009) used for verifying the parameterization accuracy. Note the large variability of the temperature profiles.
[Figure omitted. See PDF]
Figure B7
MIPAS zonal mean nighttime temperatures for 14 January 2010 (northern winter hemisphere, a), 25 March 2010 (northern spring hemisphere, b), 13 June 2010 (northern summer hemisphere, c) and 21 September 2010 (northern autumn hemisphere, d).
[Figure omitted. See PDF]
Figure B8
Zonal mean of the differences in the cooling rates of the old parameterization (a, c) and the new one (b, d) with respect to the reference cooling rates obtained for MIPAS temperatures for 13 June 2010 (northern summer hemisphere) and 21 September 2010 (northern autumn hemisphere).
[Figure omitted. See PDF]
Figure B9
The MIPAS nighttime temperature zonal mean for 15 February 2009 (see also Fig. ), used for verifying the parameterization accuracy under stratospheric sudden warming and elevated stratopause conditions. Note the high-altitude location of the stratopause in the northern polar region.
[Figure omitted. See PDF]
Code and data availability
The code is available at 10.5281/zenodo.10849970 . The parameterization is also available as a Python routine for calculating cooling rates for specific purposes at 10.5281/zenodo.10567258 . Note that the Python version is much slower than the Fortran version, and it is not recommended for use in GCMs.
The supplement related to this article is available online at:
Author contributions
MLP performed the LTE and non-LTE reference calculations, participated in the adaptation of the original parameterization, wrote the manuscript and had the final editorial responsibility for this paper. FF led (together with BF) the adaptation of the original parameterization of , wrote the code of the new parameterization, performed all the tests of the parameterization and calculated the cooling rates of the previous and new parameterizations presented here. VF made a critical contribution to the adaptation of its original algorithm and to its update. BF co-led the adaptation of the original parameterization and designed the accuracy tests. DRM provided the CO abundance and advice on the development of the parameterization for its use and implementation in climate models. All the authors contributed to the discussions and provided text and comments.
Competing interests
The authors declare that they have no conflict of interest.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
The IAA team acknowledges financial support from the Agencia Estatal de Investigación, MCIN/AEI/10.13039/501100011033, through grant nos. PID2019-110689RB-I00, PID2022-141216NB-I00 and CEX2021-001131-S. We thank two anonymous referees for their very valuable suggestions leading to an improvement of this work.
Financial support
This research has been supported by the Agencia Estatal de Investigación (grant nos. PID2019-110689RB-I00, PID2022-141216NB-I00 and CEX2021-001131-S).The article processing charges for this open-access publication were covered by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
Review statement
This paper was edited by Tatiana Egorova and reviewed by two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The radiative infrared cooling of CO
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details




1 Instituto de Astrofísica de Andalucía, CSIC, Granada, Spain
2 Institute of Atmospheric Sciences and Climate, CNR-ISAC, Bologna, Italy
3 York University, Centre for Research in Earth and Space Science, Toronto, Canada; retired
4 School of Physics & Astronomy, University of Leeds, Leeds, UK