Acta Polytechnica doi:10.14311/AP.2013.53.0832 Acta Polytechnica 53(Supplement):832–838, 2013 © Czech Technical University in Prague, 2013 available online at http://ojs.cvut.cz/ojs/index.php/ap THE SENSITIVITY OF THE GREENHOUSE EFFECT TO CHANGES IN THE CONCENTRATION OF GASES IN PLANETARY ATMOSPHERES Smadar Bresslera,∗, Giora Shaviva, Nir J Shavivb a Dept. of Physics, Israel Institute of Technology, Haifa 32000, Israel b Racah Institute of Theoretical Physics, The Hebrew University, Jerusalem 91904, Israel ∗ corresponding author: smadar@physics.technion.ac.il Abstract. We present a radiative transfer model for Earth-Like-Planets (ELP). The model allows the assessment of the effect of a change in the concentration of an atmospheric component, especially a greenhouse gas (GHG), on the surface temperature of a planet. The model is based on the separation between the contribution of the short wavelength molecular absorption and the long wavelength one. A unique feature of the model is the condition of energy conservation at every point in the atmosphere. The radiative transfer equation is solved in the two stream approximation without assuming the existence of an LTE in any wavelength range. The model allows us to solve the Simpson paradox, whereby the greenhouse effect (GHE) has no temperature limit. On the contrary, we show that the temperature saturates, and its value depends primarily on the distance of the planet from the central star. We also show how the relative humidity affects the surface temperature of a planet and explain why the effect is smaller than the one derived when the above assumptions are neglected. Keywords: greenhouse effect, radiative transfer, earth like planets, surface temperature. 1. Introduction The influence of concentration changes of atmospheric gases on the surface temperature of planetary atmo- spheres is a leading thread in current planetary re- search [7]. Due to the importance of the problem, it is desirable to have a model which can predict cor- rectly the greenhouse effect and its dependence on various changes, while at the same time being suf- ficiently simple to provide a tool to understand the details of the radiative transfer physical processes affecting this problem. We devised such a model. The model consists of two parts – radiative transfer and the molecular absorption dependent optical depth. We find that the minimum number of bands needed for a model to be faithful to the underlying physics, is a two-band semi-grey model. Consequently, we first solve the radiative transfer problem in terms of two optical depths in the two chosen bands, chosen in such a way as to provide a faithful representation of the underlying radiative transfer. We denote the two bands by “vis” and “fir” and we will specify them shortly. The radiative transfer equation is then solved in terms of the optical depths τvis and τfir to yield a universal function Tsurf (τvis,τfir). As the solution is found in terms of dimensionless quantities, it is uni- versal and does not depend on many of the planetary parameters such as its mass or specific composition of the atmosphere. Once we obtain the universal solution to the radia- tive transfer problem, we calculate the optical depths τvis and τfir from the basic molecular absorption data. As the molecular absorption is predominantly line absorption with wide windows, special care must be exercised in devising the algorithm which converts the molecular absorption coefficients κ(λ) into the above optical depths. With the universal radiative transfer solution and the optical depths, a change ΔX in the concentration of a certain gas yields a change ΔT according to: ∂Tsurf ∂X = ∂Tsurf ∂τvis ∂τvis ∂X + ∂Tsurf ∂τfir ∂τfir ∂X , (1) where the optical depths τvis and τfir will be defined shortly. Our study is unique is several ways. In particu- lar, we define the two wavelength bands according to physical properties. Since the original treatment by Simpson [10, 11], the semi-grey model has only been applied to RT problems by Thomas and Stamnes [12] and practically no other similar distinction was made. It is also crucial to note the fundamental difference be- tween stellar and planetary atmospheres. The molecu- lar absorption, with its large variation in wavelength, and the total optical depth of planetary atmospheres, are such that the atmosphere contains spectral win- dows through which the planetary radiation can leak to space almost freely. Such a phenomenon does not exist in stellar atmospheres which are hotter and in which the absorption is due to ions. Consequently, the assumption of LTE, frequently implemented in stellar atmospheres, is not really justified in planetary 832 http://dx.doi.org/10.14311/AP.2013.53.0832 http://ojs.cvut.cz/ojs/index.php/ap vol. 53 supplement/2013 The Sensitivity of the Greenhouse Effect atmospheres. Note that the assumption of LTE is usually referred to the distribution of the electrons in the various levels and to the distribution of the radi- ation field with wavelength. Many treatments treat the visible as a heat source and the IR in the diffusion approximation. We assume that the electrons satisfy the Boltzmann distribution but we do not assume that the radiation field is Planckian. In this sense, the radi- ation field in illuminated planetary atmospheres is not in LTE and it is improper to assume that F = ∇B(T) where F is the radiation flux. Separate treatment of the radiative transfer in the visible and infrared was already carried out in the past (cf. [12, section 12.3]). Our treatment differs in several points: We impose the energy conservation at each layer, we redefine the absorption coefficients, we apply the radiative transfer and not the diffusion approximation, and we find the temperature feedback directly from the radiative equation and not from the diffusion approximation for the far infrared part of the spectrum. 2. Basic assumptions We show in Fig. 1 the specific intensities of the in- solating star and the thermal emission of the planet. We define λrad as the wavelength at which the two intensities are equal, namely Ip(λrad) = (1 −a) 4 ( R∗ d )2 I∗(λrad), (2) where Ip and I∗ are the specific intensities of the planet and the insolating star at the top of the planetary at- mosphere. a is the mean albedo. Energy conservation implies that 1 CV dQ(z) dt = ∫ ∞ 0 κ(λ,z) [J(λ,z) −B(λ,z)] dλ, (3) where J is the mean specific intensity. A further requirement of steady state implies that dQ(z)/dt = 0. As apparent in the figure, J(λ) � B(λ(Tatm)) for λ < λrad. Consequently, we can split the energy integral, and write ∫ λrad 0 κ(λ,z)J(λ,z)dλ = ∫ ∞ λrad κ(λ,z) [J(λ,z) −B(λ,z)] dλ (4) The first term is positive definite and hence always represents heating. The second term must therefore be negative and thus represents cooling. Clearly, the two wavelength ranges describe different phenomena. At this point, we can also state that J = B, which is the condition for LTE, only if the first integral vanishes, namely there is no heating of the atmosphere by absorbed radiation. Our semi-grey model therefore has two bands, the vis band in the range up to λrad, and the fir band for wavelength In te ns ity /A bs or pt io n Central star emission Planet emission !(") "cut "rad Figure 1. The definition of λrad, where the specific in- tensity of insolation equals the specific intensity of the planetary thermal emission. Also shown is the defini- tion of λcut, the wavelength above which the molecular absorption (denoted by κ) becomes significant. λ > λrad. The radiative transfer equation is solved in the two stream approximation. However, while the optical depths τvis and τfir are constant over the respec- tive wavelengths, we allow I(λ) to change with λ. We assumed in the calculations reported here: Nvis = 200 wavelengths in the range (103Å÷λrad) and Nfir = 400 wavelengths in the range (λrad ÷ 8 × 105Å). The at- mosphere was divided into 50 slabs. Each slab has an optical depth of τvis/50 for λ < λrad and an optical depth τfir/50 for λ > λrad. The Nvis and Nfir wave- lengths were distributed logarithmically. The energy condition was used to calculate the temperature of slab i. The energy condition in our calculation is given by: τvis Natm j=Nvis∑ j=1 [Ji,j −B(Ti,j)] (5) + τfir Natm j=Nfir∑ j=1 [Ji,j −B(Ti,j)] = 0 where Natm is the number of layers in the atmosphere (50 in our case) and the index i runs over all layers in the atmosphere. The above condition must be satisfied for every i. Nvis and Nfir are the number of wavelengths we use in the respective wavelength range. Ji,j is J at atmospheric layer i and wavelength j and B(Ti,j) is the Planck function for temperature Ti at layer i and wavelength λj. Although it can be eliminated, Natm is kept for clarity. There are Natm such equations to solve for the temperature Ti in each layer i. Note that we kept the Planck function in the λ < λrad range despite the fact that it is very small. We assume that the total optical depth is divided equally in all the layers of the atmosphere. This implies that the physical width of the atmosphere varies with height. Moreover, we work with the two τ’s as the independent variables, not κ. 833 Smadar Bressler, Giora Shaviv, Nir J Shaviv Acta Polytechnica Figure 2. The surface temperature as a function of τfir for a fixed τvis. Also shown is the sky temperature Tsky. 0.01 0.1 1 10 100 1000 300 400 500 600 rad Åλ =20,000 rad Åλ =10,000 rad Åλ =50,000 Su rf ac e T em pe ra tu re K τfir>>>1 τvis Figure 3. The effect of τvis on Tsurf for three cases. 3. Greenhouse and Anti-Greenhouse models Figure 2 illustrates a typical case, where τfir increases indefinitely while keeping τvis fixed. The surface tem- perature increases slowly for small τfir’s (i.e., per given logarithmic increase of τfir), but when τfir ∼ 1, the rate of increase becomes larger, up to about τfir ∼ 100, where the surface temperature saturates and levels off. The reason is that as the temperature rises, the peak of the Planck spectrum progressively moves to- wards shorter wavelengths and thermal radiation leaks through the vis range (cf. [9]). Thus, the greenhouse effect does not experience a runaway when the con- centration of any gas, and its corresponding optical depth, increase indefinitely. Clearly, Simpson’s para- dox does not exist when the problem is solved properly. Previous attempts to solve the paradox assumed the existence of windows in the absorption coefficient. The saturation is obtained, however, even if no windows are assumed. Moreover, the saturation appears even without assuming various feedbacks like increased albedo etc. The figure also depicts the sky temperature Tsky, defined such that σT 4sky is equal to the radiation flux from the atmosphere to the surface of the planet. As Tsky < Tsurf in this equilibrium model, the surface cools by exchanging energy with the atmosphere. Figure 3 demonstrates the effect of a changing τvis on the surface temperature, for extremely large τfir for which the saturation temperature is reached. We show that τvis is sufficiently powerful to create an anti-GHE even under the most adverse condition. As long as τvis ≤ 1 the effect of τvis is negligible, but for larger values, the effect is very noticeable. In the limit of τvis � 100, the decrease in temperature approaches an asymptotic value. It is interesting to note that irrespective of λrad, the minimal temperature reached is the same. Finally, the saturation temperature is not a monotonic function of λrad. 4. The resolution of the Simpson paradox In 1927 Simpson treated the radiative transfer prob- lem in planetary atmospheres, and tacitly assumed: (a) that λcut ≡ λrad. (b) the existence of LTE for the long wavelength range. Under these assumptions, he found that the temperature of the atmosphere is given by: T 4p (τ) = ( 1 + 3τ 4 ) (1 −〈a〉) 4 ( R∗ d )2 T 4∗ , (6) where τ is the mean optical depth for λ ≥ λrad, which is measured from the top of the atmosphere down- ward. 〈a〉 is the mean albedo. Simpson did not specify how τ is evaluated, presumably it was the Planck mean. In particular, the temperature near the surface is obtained by substituting τ = τtot. Obvi- ously, as τtot →∞, so does Tp. For sufficiently large τtot (τtot ≈ 3.8 × 105), the temperature of the planet reaches the temperature of the central star and can even surpass it. We note that there are particular wavelength ranges in the far IR, for which the total optical depth is as high as 104. The possibility of a temperature runaway, as predicted by the simple Simpson solution, was coined the Simpson paradox. There were attempts to resolve the paradox by assum- ing the existence of windows in the fir. However, it is obvious that as long as the vis band is ignored and the radiative transfer is solved with a single band, there is no way to eliminate the paradox. Our treatment (see [9]) solves the problem, as we demonstrate that for moderate optical depths of τfir, the temperature saturates. 5. The (τvis,τfir) plane The main result is shown in Fig. 4 where the con- tour lines of constant Tsurf in the (τvis,τfir) plane are plotted. The classical picture of the greenhouse ef- fect, where τvis does not exist or is merged with τfir, is along the τfir axis, where the surface temperature rises monotonically and saturates. As long as τvis < 1, its effect is small, but as τvis increases, its effect be- comes more prominent. For τvis ∼ 10, it reduces the saturation temperature. 834 vol. 53 supplement/2013 The Sensitivity of the Greenhouse Effect 200 250 300 350 400 Greenhouse Saturation Anti-Greenhouse - 4 - 3 - 2 - 1 0 1 2 - 1 0 1 2 3 Log10 Τvis L og 10 Τ fi r Figure 4. Curves of constant Tsurf in the (τvis,τfir) plane. The regions of saturation, greenhouse and anti- greenhouse are marked. The (τvis,τfir) plane allows us to determine the effect that the change in concentration of an atmospheric constituent has, once we evaluate the two optical depths τvis and τfir and their changes with concentra- tion. In the next section we show how to evaluate the mean optical depth of the relevant bands. 6. The algorithm for the optical depth Both the Planck and the Rosseland means are poor averages of the absorption when it comes to planetary atmospheres, where molecular absorption dominates. Two factors play here, the existence of spectral win- dows and the large variation as a function of wave- length. In any averaging of this sort, the optimal weight function is the one which is the closest to the actual solution, which is of course not known. Never- theless, it is imperative to have a good guess for the weight function, or else the results may be completely skewed. The failure of the Rosseland mean is a good example for a poor weight function which yields an unacceptable result. Consider first the vis range. Since the temperature of the radiation is that of the central star and hence very high relative to that of the planetary atmosphere, the zeroth solution for the transmission of specific intensity I(z,ν) is given by I(z,ν) = ITOA∗,ν e −κ(ν)z. (7) Here ITOA∗,ν is the stellar specific intensity at the top of the atmosphere. To secure the energy flux transfer through the atmosphere, we therefore write that:∫ ν2 ν1 ITOA∗,ν e −κ(ν)zdν = e−〈τvis〉 ∫ ν2 ν1 ITOA∗,ν dν. (8) Next we note that the radiation interacts with the atmosphere, which is at temperature Tatm, and the stellar radiation is to a good approximation that of a black body at temperature T∗, so we have: 〈τvis〉 = − log (∫ν2 ν1 B(T∗,ν)e−τ(ν,Tatm)dν∫ν2 ν1 B(T∗,ν)dν ) , (9) where ν1 and ν2 correspond to the boundaries of the vis range, Tatm is the temperature of the atmosphere and τvis is the total average optical depth for this range. It is important to note that the temperature in the weighting function is not that of the atmosphere through which the radiation passes, but that of the insolating star – the sun in the particular case of the Earth or the star in the general case. The optical depth, ∫Z 0 κ(λ,T,P)dλ, however, is calculated with the temperature of the atmosphere. The expression so obtained reduces to the trivial results in various limits. If there is no absorption in the vis range, then τvis = 0. When the optical depth is constant, then τvis is equal to this constant as well. 6.1. Transition in the far infrared domain Consider now the radiative transfer in the fir range. Let τ be measured from the top of the atmosphere downwards. If we write I+(τ) as the thermal flux towards larger optical depths (downwards) and I−(τ) as the flux towards smaller optical depths, then the solutions for I±(τ) under the above approximations are: F(τ) ≡ I−(τ) − I+(τ) = const. (10) E(τ) ≡ I−(τ) + I+(τ) (11) = [I−(τ) − I+(τ)] τ + const. By comparing the conditions at the top (τ = 0) to the bottom (τ = τtot), we obtain I−(τtot) − I+(τtot) = I−(0) − I+(0), (12) I−(τtot) + I+(τtot) = [I−(0) − I+(0)] τtot + [I−(0) + I+(0)] . The boundary conditions we have are: I+(0) = I∗,fir and I−(τtot) = Ip,fir, (13) where I∗,fir is the insolation for λ > λrad and Ip,fir is the planet’s emission at λ > λrad. It is generally assumed that I∗,fir = 0. However, if λrad decreases significantly, this assumption may no longer be justi- fied. Next we consider the thermal equilibrium of the surface, i.e., total absorption equals the total emission: (1−a(λ))I∗,vis +Iatm,↓ = (1−a(λ))Ip,vis +Ip,fir, (14) where I∗,vis is the insolation for λ < λrad, Iatm,↓ is the emission of the atmosphere towards the surface 835 Smadar Bressler, Giora Shaviv, Nir J Shaviv Acta Polytechnica (at λ > λrad), Ip,vis the planet’s emission at λ < λrad and a is the albedo at the short wavelengths. Our main point is that Ip,vis must be included at relatively high surface temperatures. Using the two sets of Eqs. 12 and 13, the thermal equilibrium becomes: (1 −a(λ)) (I∗,vis − Ip,vis) = 2(Ip,fir − I∗,fir) 2 + τfir,tot . (15) From the above set of equations, we can derive an expression for the average net fir flux (per unit fre- quency) over a finite band Δν, and define an effective opacity through the following: ΔIfir = 1 Δν ∫ ν2 ν1 [Ip,fir(Tp) − I∗,fir] 1 + 3τ(ν)/4 dν (16) ≈ [ Ip,fir(Tp) − I∗,fir ] Δν ∫ ν2 ν1 dν 1 + 3τ(ν)/4 ≡ [ Ip,fir(Tp) − I∗,fir ] Δν 1 1 + 3τfir,tot/4 , that is, τfir,tot = 4 3  ∫ν2ν1 B(Tatm,ν)dν∫ν2 ν1 B(Tatm,ν)dν 1+3τtot(ν)/4 − 1   . (17) This expression for the grey absorption is useful be- cause it encapsulates the different behaviors in the fir. In particular, wavelength regions for which the optical depth is small allow for a larger flux and there- fore receive a larger weight in the averaging. This is present in the Rosseland mean as well. However, unlike the Rosseland mean, the flux does not diverge if the wavelength region becomes optically thin. In such a case, the emission saturates at its surface emission. Thus, the mean can adequately describe the effect of spectral windows. 7. Actual calculation The data of molecular absorption was taken from the HITRAN compilation [6]. The procedure for calculat- ing the line absorption is described in the manual of this compilation. 8. The effect of water vapor We calculated the optical depths for increasing degrees of column densities of water molecules. The results are shown in Fig. 5. The effect of increasing the column density of water vapor is shown by the green line. The interesting phenomenon is that the curve starts in the region for greenhouse effect and continues for sufficiently high amounts of water, towards the anti- greenhouse domain. The water curve is not vertical but has a slope. The pink arrow denotes the result that a model which does not distinguish between the two domains, the vis and fir, would yield. As we can see, the inclusion of the effect of τvis lowers the predicted 5 2.0 1.5 1.0 0.5 0.0 0.5 1.0 0.4 0.2 0.0 0.2 0.4 0.6 0.8 1.0 330 320 310 300 290 280 270 260 250 230 360 vis fi r 350 340 Figure 5. The effect of gaseous water molecules in the (τfir,τvis) plane. The points are increasing column density starting from zero to the maximum amount of water vapor (relative humidity of 100 %). The green line depicts the increase in Tsurf while the pink arrow depicts the result that would have been obtained from the classical approach, which neglects τvis. The arrow marks the location of the Earth according to the Earth mean column density of water vapor as given by Crisp [2]. surface temperature. The arrow marks the location of the Earth for a column density of 8.12 × 1022#/cm2, as given by Crisp [2] for the mean Earth. 8.1. Line-by-Line models Of all the methods applied to describe the greenhouse effect of any gas on the atmospheric temperature [3] the Line-by-line (LBL) models are the most elaborate. These are not full radiative transfer models but radi- ation transmittance models. These models calculate the absorbed energy by each absorption line. The LBL model treats only one downward stream of radia- tion. In the vis, it is pure insolation, whereas beyond λrad, it is the downward self emission from the top of the atmosphere (TOA) and subsequent lower layers, in response to the long wavelength radiation emitted from the planetary surface. The absorbed energy so calculated (in terms of W/m2) is transferred to a global circulation model (GCM), where the effect on the temperature is calcu- lated. Since in this way the increase in concentration of any absorber always leads to increased energy ab- sorbed, the results are always positive, namely heating. Less detailed methods, like the correlated-k distribu- tion are even less accurate. 836 vol. 53 supplement/2013 The Sensitivity of the Greenhouse Effect 9. A maximum temperature for a planet The fact that the greenhouse saturates implies that a planet at a given distance from the central star has a maximal temperature, the temperature of saturation. This is a strict limit which does not depend on the pa- rameters of the planet like atmospheric composition or the mass or structure of the atmosphere (like pressure and density). The saturation temperature depends only on the distance of the planet from the central star and the mean albedo. Hence, this limit can serve to distinguish between brown dwarfs revolving around a central star and a planet. This model should be further developed for jovian planets, to distinguish them from brown dwarfs. This should be done by changing the lower boundary conditions. 10. Summary The radiative transfer model for planetary atmo- spheres presented here enjoys simplicity yet it does not compromise the fundamental physics of the green- house effect. The prediction of the greenhouse effect is correct and reliable, as we saw with the prediction of the location of the Earth in the (τvis,τfir) plane. In general, having the universal solution Tsurf (τvis,τfir) allows an easy determination of the surface tempera- ture of the atmosphere, irrespective of the particular composition. We demonstrated the solution in the case of water vapor, and have shown that water va- por do not lead to a runway greenhouse effect, and even does not drive an Earth-like atmosphere into the saturation region. The model generalizes the greenhouse and the anti- greenhouse effects and generates a comprehensive pic- ture of the phenomenon. Radiative models which are essentially transmit- tance models like the LBL, cannot predict the result- ing temperature changes, as they do not impose an energy balance equation. Such models can predict how much energy is absorbed by a given atmospheric struc- ture but they do not have the feedback to evaluate the resulting temperature changes. For this purpose, one has to resort to a General Circulation Model, which is a dynamic model and not static. The saturation of the greenhouse effect leads to the existence of a strict limit to the temperature of a planet. This limit can help to distinguish between planets and brown dwarfs. Acknowledgements Special thanks to Dr. Eyal Bressler and to Jenia and Ram Oren. References [1] Bressler,S., Shaviv, G., Shaviv, N.J. 2012, in preparation [2] Crisp, D. chap. 11 in Allen Astrophysical Quantities, 2000 [3] IPCC Climate Change 2007: Synthesis Report [4] Kasting, J.F., 1988 Runaway and moist greenhouse atmospheres and the evolution of Earth and Venus Icarus 74, 472 [5] Mihalas, D. 1970 Stellar atmospheres, W.H. Freeman and Company, San Francisco [6] Rothman, L. S. et al. The HITRAN 2008 molecular spectroscopic database, 2009 JQSRT 110, 533 [7] Seager, S. 2010 Exoplanet Atmospheres, Princeton University Press, New Jersey [8] Shaviv, G. & Wehrse, R., 1991, Continuous energy distributions of accretion discs, A&A, 251, 117 [9] Shaviv, N.J., Shaviv, G. & Wehrse, R., 2011, The maximal runaway temperature of Earth-like planets, Icarus, 216, 403 [10] Simpson, G.C., 1927. Some studies in terrestrial radiation. Mem. R. Meteorol. Soc. II 1015 (16), 69–95. 1016 [11] Simpson, G.C., 1928. Further studies in terrestrial radiation. Mem. R. Meteorol. Soc. 1017 III (21), 1–26 [12] Thomas, G.E. and Stamnes, K. (1999) Radiative Transfer in the Atmosphere and Ocean, section 12.3, Cambridge University Discussion Jim Beall — There is a correlation between tempera- ture and CO2 in the industrial period between 1880–2004. The CO2 goes from 280 ppm to 380 ppm. But the temper- ature only increased by 3/4 °C during the interval from 1880–2004. Smadar — Different analyses give a wide range for the anthropogenic contribution to the temperature rise, and even the measured temperature rise itself has a large uncertainty, such that our results all fall within these limits. It is only with better and more detailed models that we will be able to rule out (through this methodology) possible misconceptions about the extent that anthropogenic CO2 is the sole cause for heating of the atmosphere. Maurice van Putten — Is it known how much CH4 is released when there is a 1 °C increase? The result, I expect, can be readily included in your model. Smadar — You are correct that CH4 is the next runner up in our model after CO2, being the third important greenhouse gas in Earth’s atmosphere. The answer to your question can be obtained using our model in a similar manner to that of CO2 and water. On the one hand, we expect a smaller effect due to the still relatively low column density of CH4 (about 200 times less than CO2). Its strong absorption, on the other hand, might compensate for this, and it is hard to predict what will be its final effect without calculating it. We have developed a greenhouse indicator aimed exactly at comparing different greenhouse gases by their slope over the τfir and τvis plane, due to changes in concentration. For example, from the presented data, it seems that CO2 has a stronger greenhouse effect than 837 Smadar Bressler, Giora Shaviv, Nir J Shaviv Acta Polytechnica water, at least for some humidity regions. So we should wait and see. Bozena Czerny — 1. What about cloud coverage? This varies across the surface and it is also coupled vertically with high cloud reflectivity. 2. What about mechanic effects like convection or winds? This likely predominantly cools the surface. Smadar — 1. Cloud coverage has not yet been consid- ered, as it is a complicated matter: clouds at different heights have different albedos, and also the complication of combining the surface albedo with partial cloud coverage, changing the effective albedo. Initial experiments with changing surface albedo, though, show very low sensitivity of the model to even 30 % changes, but this is yet to be determined carefully in further calculations. 2. Convection was still not considered in this preliminary model, only pure radiative transfer. As convection will onset only at the steepest temperature gradient, which is the adiabatic limit, our calculation yields the upper limit of the temper- ature. We can parameterize convection from our model by translating the dT/dτfir into dT/dz through a linear scale factor l = (d ln τ dz )−1 . (18) This will be carried out in further studies. As to winds, our model is not suitable for calculating the effect of global winds and wind jets (in gravitationally locked planets especially), and it requires coupling with climate global circulation models (GCM), which are used today to calcu- late the greenhouse temperatures from LBL-models, which yield radiation fluxes in W/m2. We hope to match our radiative transfer software in the future with a good GCM, in order to feed our results to climatic models. As you see, there is still much work to be done! We are only at the beginning . . . 838 Acta Polytechnica 53(Supplement):832–838, 2013 1 Introduction 2 Basic assumptions 3 Greenhouse and Anti-Greenhouse models 4 The resolution of the Simpson paradox 5 The (tau_vis, tau_fir) plane 6 The algorithm for the optical depth 6.1 Transition in the far infrared domain 7 Actual calculation 8 The effect of water vapor 8.1 Line-by-Line models 9 A maximum temperature for a planet 10 Summary Acknowledgements References