ORIGINAL RESEARCH article

Front. Complex Syst., 12 August 2025

Sec. Multi- and Cross-Disciplinary Complexity

Volume 3 - 2025 | https://doi.org/10.3389/fcpxs.2025.1617092

Unsettling the settled: simple musings on the complex climatic system

  • Department of Water Resources and Environmental Engineering, School of Civil Engineering, National Technical University of Athens, Athens, Greece

Abstract

Our revisit of fundamental issues of climate challenges the notion and term of the “greenhouse effect”, and attempts a scientific reevaluation using minimal assumptions, such as Newton’s laws, maximum entropy and gas spectroscopy. It replaces terms like “greenhouse gas” with “radiatively active gas” (RAG) and “greenhouse effect” with “atmospheric radiative effect” (ARE). While ARE exists in several planets’ atmospheres, on Earth it is primarily driven by water vapor and clouds, with CO2 playing a minor role (especially anthropogenic CO2 which represents 4% of total emissions). Equilibrium thermodynamics, via entropy maximization or molecular collision simulation, leads to an isothermal atmosphere at about 250 K (the average temperature of the troposphere and stratosphere) irrespective of RAG presence or not. It is the troposphere’s 6.5 K/km temperature gradient (lapse rate), partly shaped by moist adiabatic processes, that drives the atmosphere away from this equilibrium and warms the surface to about 288 K on average, with ARE (mainly water vapor and clouds) contributing to the warming, but only when this gradient exists. The temperature gradient varies spatially and temporally and, since 1950, has weakened in the tropics and grown in the polar areas, resulting in a decrease of the surface equator-to-pole gradient, as expected in global warming conditions.

1 Introduction

Do we live in a greenhouse? Most people would give an affirmative reply given the enormous campaign to promote the climate initiatives that are founded on the concepts of greenhouse gases (GHG) and greenhouse effect (GHE). Systematic searches in Google Books (https://books.google.com/ngrams/) and Internet Archive (https://archive.org/) show that, while the term ‘greenhouse’ is old (already appearing in the 18th century with its literal meaning), the terms GHG and GHE are newer. Following Poynting (1907) who used the term GHE for planetary atmospheres (even though he was criticized by ; see Supplementary Appendix SA), the term GHE as it relates to planetary atmospheres was popularized in 1960s by NASA (; ). Since the late 1970s, however, the usage of GHE has changed to closely associate it with carbon dioxide (CO2) and its emission into the atmosphere from the burning of fossil fuels, but also to imply that this causes disastrous effects on climate, economy and every aspect of life (; ; ; ; ; ; ; ).

From a scientific perspective, clarity in terminology is essential (cf. the Aristotelian notion of sapheneia; ). In strict scientific terms, the reply to the question whether we live in a greenhouse is negative: The Earth’s atmosphere does not function like a greenhouse. A simple revisit of the definition of a greenhouse (a structure that is designed to regulate the temperature and humidity of the environment inside) suffices to see this, even though effects such as those that the term GHE purports to describe are present on Earth, as well as on other planets in the solar system with atmospheres (see section 3.2). Hence, the term GHE can be misleading as further explained in Supplementary Appendix SA.

For scientific clarity and accuracy, we propose replacing the term GHG with radiatively active gas (RAG), comprising water vapor (WV) and non-condensing radiatively active gases (NC-RAGs). We also propose replacing the term GHE with atmospheric radiative effect (ARE). ARE is the influence of RAGs on the energy balance in the atmosphere by absorbing, emitting, and scattering radiation across both shortwave (SW) and longwave (LW) spectra. These processes influence the flow of energy between the Earth’s surface, the atmosphere, and space. While the term GHE is often associated primarily with heat-trapping, ARE explicitly recognizes all radiative interactions in the atmosphere. It avoids the misleading analogy of a greenhouse and is more applicable to a scientific framework that considers both SW and LW radiation. For the reader’s convenience, these acronyms are listed and explained in Table 1, along with all other abbreviations used in this paper.

TABLE 1

AcronymExplanation and comments
ACSAmerican Chemical Society
AIPAmerican Institute of Physics
AREAtmospheric Radiative Effect (proposed to replace GHE)
CERESClouds and the Earth’s Radiant Energy System
CLIMEXPClimate Explorer (by the Koninklijk Nederlands Meteorologisch Instituut)
ECMWFEuropean Centre for Medium-Range Weather Forecasts
ERA5Fifth generation ECMWF ReAnalysis
GHCN-DGlobal Historical Climatology Network – Daily
GHEGreenhouse Effect (proposed to be replaced by ARE)
GHGGreenhouse Gas (proposed to be replaced by RAG)
HITRANHigh-Resolution Transmission (a molecular spectroscopic database)
ICAOInternational Civil Aviation Organization
IPCCIntergovernmental Panel on Climate Change
LWLongwave (radiation)
NASANational Aeronautics and Space Administration (United States)
NC-RAGNon-Condensing Radiatively Active Gas (e.g., CO₂)
NOAANational Oceanic and Atmospheric Administration (United States)
RAGRadiatively active gas (proposed to replace GHG)
REISRadiative Energy Imbalance Slope
RRTMRapid Radiative Transfer Model
SWShortwave (radiation)
TOATop-of-Atmosphere
WMOWorld Meteorological Organization
WRITWeb-based Reanalyses Intercomparison Tools (by NOAA)
WVWater Vapor

List of acronyms used in this paper.

Notably, WV and clouds (for which WV is responsible) dominate ARE, while CO2 contributes only 4%–5% to it (). Also, anthropogenic CO2 emissions are only 4% of the total, with the vast majority (96%) being natural (). Additionally, evidence suggests that changes in temperature precede those in CO2 concentration, thus challenging the assumption that CO2 drives temperature (; ). For this reason, here we avoid using popular terms such as “radiative forcing” (e.g., ) which imply that CO2 is the temperature driver.

These recent developments, summarized in and in a review paper led by a chatbot (large language model), , question the mainstream view that climate science is settled, particularly with respect to the CO2 and the human contribution to the increase of its atmospheric concentration in the last decades. Arguably, science can never be settled (cf. ), particularly that dealing with one of the most complex systems on Earth, the climatic system, which remains under active investigation. It is recalled that the climatic system consists of the atmosphere, the hydrosphere (including its solid phase—the cryosphere), the lithosphere and the biosphere, which mutually interact and respond to external influences (system inputs).

It is thus healthy to question even the fundamental ideas about climate, leaving aside the fact that they are established or settled. In this respect, this paper examines (in

Section 3

) the following fundamental questions, using empirical data and minimal assumptions, and aiming to clarify the roles of ARE and atmospheric dynamics without preconceived notions:

  • 1. Is there an empirically verified ARE in the atmosphere?

  • 2. What would be the atmospheric temperature profile in equilibrium?

  • 3. What is the empirically observed temperature profile, and does it differ from the equilibrium profile?

  • 4. Is the ARE responsible for the observed temperature profile in the atmosphere?

  • 5. What factors can explain the observed temperature profile?

  • 6. What is the relative importance of ARE and temperature gradient in climate dynamics?

  • 7. What are the recent changes in the atmospheric temperature gradient?

Note that the terms “temperature gradient” or simply “gradient” are used above and below as shortcuts for “(minus) vertical gradient of temperature in the troposphere” (also known as “lapse rate”). Whenever a different gradient direction is assumed, we specify it (e.g., surface equator-to-pole temperature gradient).

2 Data and software

The temperature data used here are taken from the ERA5 Reanalysis on a monthly scale. ERA5 stands for the fifth generation atmospheric reanalysis of the European Centre for Medium-Range Weather Forecasts (ECMWF; ECMWF ReAnalysis). Its data are publicly available for the period 1940 onwards at a spatial resolution of 0.5°. The data sets used here were retrieved from the Physical Sciences Laboratory platform of the US National Oceanic and Atmospheric Administration (NOAA) (WRIT: Monthly Timeseries; NOAA Physical Sciences Laboratory, https://psl.noaa.gov/cgi-bin/data/atmoswrit/timeseries.pl). For the period 1950 onwards, they are also accessible from the Climate Explorer (CLIMEXP) platform (https://climexp.knmi.nl/) which in addition allows extended processing of the data.

Radiation data from satellites were retrieved from NASA’s ongoing project Clouds and the Earth’s Radiant Energy System () from the Terra platform (operational since January 2001). The specific product used here is the CERES SSF1deg monthly averaged TOA LW radiative fluxes at a 1°-regional grid, constant-meteorology-temporally-interpolated. The top-of-atmosphere (TOA) fluxes are provided for clear-sky and all-sky conditions and are available online.

The main software tool used here is RRTM (standing for rapid radiative transfer model), a hybrid physical/statistical approximation of the full line-by-line models, developed for climate studies. Its results have been extensively evaluated and rigorously validated (). The model implementation used here is an interactive web application hosted by the University of Chicago (https://climatemodels.uchicago.edu/rrtm/) that simulates the radiation flux through Earth’s atmosphere at a range of altitudes, calculating both SW (incoming and reflected sunlight) and LW (emitted by the ground and atmosphere), both upward and downward. Users can adjust parameters such as solar input, albedo, and atmospheric composition in order to study their effects on Earth’s energy balance, revealing whether the planet gains or loses energy and enabling the parameters that lead to energy balance to be found.

Additional software was developed in-house to simulate the thermal equilibrium dynamics of a gas and to determine the altitude-dependent distribution of air density and temperature under gravity. The simulation is based solely on molecular collisions treated as perfectly elastic hard-sphere interactions, without any assumptions beyond Newton’s laws. A detailed description is provided in Supplementary Appendix SD, and the software is available online as Supplementary Software, allowing interested readers to review the code, replicate the calculations, or simply view the simulation videos (see Supplementary Video).

3 Analysis of questions

3.1 Is there an empirically verified ARE in the atmosphere?

The reply to this question is definitely affirmative and is based on observations. As noted by , a pioneer of the measurement and modelling of radiation, the first observations relating to the problem of Earth’s (LW) radiation to space were made between the years 1780 and 1850. Ångström (1915, p. 16) cites some sporadic measurements by several researchers for the period 1887-1912. From these measurements and from experiments in the 19th century, it was understood that the major constituents of the atmosphere, i.e., nitrogen (N2) and oxygen (O2), are transparent to LW radiation. In contrast, some minor constituents, particularly WV, carbon dioxide (CO2) and ozone (O3), absorb and reemit LW radiation, thus being RAG. pioneered our understanding of the dominance of WV in this process, of the importance of its presence for climate and life, and of the minor contribution of CO2 in the ARE.

Ångström (1915, p. 16) provided his own systematic measurements, as well as a rough quantification of ARE in the following way:

a surface at 15°C temperature ought to radiate 0.526 cal. If the observed effective radiation does not amount to more, for instance, than 0.15 cal, this must depend upon the fact that 0.376 cal is radiated to the surface from some other source of radiation. In the case of the earth this other source of radiation is probably to a large extent its own atmosphere.

Since the introduction of celebrated equation for natural evaporation from an open water surface, the quantification of the ARE has become part of the routine hydrological calculations for real-world problems, such as those of the hydrological balance (of which evaporation represents a substantial component) and the irrigation needs in agricultural applications. equation is based on the presence of WV in the atmosphere and disregards that of NC-RAGs such as CO2. Notably, while the role of CO2 in photosynthesis is important in biochemical terms, it becomes negligible in terms of its contribution to the surface energy balance. A recent study by analyzed a large set of LW radiation measurements distributed in time across a century, in which CO2 has escalated from 300 to about 420 ppm. They concluded that the observed increase of the atmospheric CO2 concentration has not altered the ARE in any discernible way. Thus, the ARE continues to be dominated by the quantity of WV in the atmosphere and CO2 remains insignificant in the ARE.

Since the 1960s, detailed spectroscopic studies of RAGs were conducted, initially developed for military purposes and subsequently expanded to cover broader scientific use, evolving into an essential tool for atmospheric and astrophysical research. The results of these studies were compiled into HITRAN (High-Resolution Transmission), a molecular spectroscopic database designed to support the study of electromagnetic transmission and emission in gaseous media (; ; ; ). The database includes detailed line-by-line spectroscopic parameters for various molecular species.

Experimental verification of the theory at laboratory has been provided by and (see Supplementary Appendix SA). Empirical atmospheric data for testing the applicability of the theory had been provided early by based on a Michaelson interferometer in a satellite. Recently, and using these data showed that theory is impressively consistent with the observations, as shown in Figure 1, reproduced from the latter publications.

FIGURE 1

(or Figure 10 in ): Vertical intensities Ī(0) at the top of the atmosphere observed with a Michaelson interferometer in a satellite (), and modeled with radiation transfer theory for the Sahara desert, Mediterranean and Antarctica. The “frequency” referred to in the horizontal axis should better read as “wavenumber”. The intensity unit in the vertical axis is 1 i. u. = 1 mW m−1 cm sr−1. The labels (a-f) are explained in each of the panels.

Some discrepancy appears in the Antarctica case, as the observed radiation is higher than that emitted by the Earth’s surface. The explanation given by and is this:

Radiative forcing is negative over wintertime Antarctica since the relatively warm greenhouse gases in the troposphere, mostly CO2, O3 and H2O radiate more to space than the cold ice surface at a temperature of T = 190 K, could radiate through a transparent atmosphere.

We note though that the T = 190 K, as well as the entire atmospheric profile, are estimated, rather than observed. mention that

no radiosonde data are available for comparison […] The profile can be anticipated intuitively from an examination of the spectrum […] The cold surface temperature results in the relatively low radiances found in the window region and between the water vapour lines.

Further examination of systematic ground temperature observations in the area shows that the assumed value of T = 190 K is unlikely. Specifically in the South Pole (station Amundsen-Scot; coordinates: 90°S, 0°E, prob: 2,770 m), the lowest daily temperature in April 1970 (the month of the experiment) did not fall below 204K. The Vostok station (coordinates: 78.45°S, 106.87°E, 3,420.0 m or prob: 3,468 m) is the location with the lowest temperatures among 41 GHCN-D stations in Antarctica, and its temperatures in April are, on average, 7 K lower than in Amundsen-Scot. Furthermore, the reconstructed vertical temperature profile given by does not seem plausible as it has an abrupt increase of >30 K at the level of 700 hPa.

On the other hand, the ERA5 Reanalysis data suggest that there is a permanent state of temperature inversion (increasing temperature with increased altitude) between the atmospheric levels of 700 and 600 hPa (Figure 2), albeit not that big. This explains the discrepancy seen in Figure 1. The presence of the temperature inversion over the Antarctic Plateau has been studied by , who found that, along with the scarcity of free tropospheric water vapor, it causes a negative greenhouse effect. In addition, stated that in the case of temperature inversion, RAGs actually cool, rather than warm, the atmosphere. (See additional explanation in Section 3.4 below.).

FIGURE 2

3.2 What would be the atmospheric temperature profile in equilibrium?

In their recent work, stated: “In the absence of greenhouse gases, the isothermal atmosphere will not change with time, since there is no thermal gradient to drive heat flow”. They also asserted that “without greenhouse gases, a thermally isolated adiabatic atmosphere would slowly evolve to an isothermal atmosphere because of conductive heat flow from the warmer lower atmosphere to the colder upper atmosphere”. Here we examine these statements with a slightly different formulation as seen in the question in the heading of this subsection.

It is well known that a gas in equilibrium tends to become isothermal, but here we examine the question also including gravity in the analysis. In fact, this problem is very old and is known as the Loschmidt paradox after who asserted that gravity could counteract the temperature uniformity implied by the second law as the kinetic energy of the molecules of a gas should be lower at higher elevation in compensation for the higher potential energy. The paradox puzzled many for years, yet both James Clerk Maxwell and Ludwig Boltzmann reacted negatively, supporting the isothermal hypothesis, as detailed in a recent analysis by .

Nonetheless, we deemed it useful to perform our own investigation in reply to this question independently of past arguments and using two methods. First, we maximize the entropy of a single molecule which is in motion in a vertical column under gravity. The single-molecule method has been successfully used before to derive the Clausius-Clapeyron equation (). The development and application of the method to a spherical monoatomic and to a diatomic molecule are presented in Supplementary Appendices SB, SC, respectively. In brief, the thermal (internal) energy per molecule turns out to be the same for all altitudes, equal to , where is the number of degrees of freedom of a molecule and ε is its total energy. This results in for monoatomic molecules and for diatomic molecules. The remaining energy is the average dynamic energy due to gravity. The temperature is also constant, independent of the altitude, equal towhere k is the Boltzmann constant. In other words, entropy maximization results in an isothermal atmosphere.

The same result emerges when we approach the problem with our second method, without invoking the principle of maximum entropy—by instead simulating the motion of molecules, modeled as perfectly elastic hard spheres, using only Newton’s laws. Once again, the temperature is found to remain constant with altitude, as shown in Figure 3 and discussed in detail in Supplementary Appendix SD. In other words, none of the approaches predicts a decrease in temperature with altitude, which is the observed behavior as will be explained in Section 3.3.

FIGURE 3

3.3 What is the empirically observed temperature profile, and does it differ from the equilibrium profile?

It is common knowledge that Earth’s atmosphere is not isothermal. The temperature decreases with altitude in the troposphere, then it remains constant in the lowest layer of the stratosphere and continues to change at even higher altitudes. The typical atmospheric temperature profile is represented by the standard atmosphere adopted by the International Civil Aviation Organization (). According to the , the standard atmosphere is “A hypothetical vertical distribution of atmospheric temperature, pressure and density which, by international agreement, is roughly representative of year-round, midlatitude conditions.” It is based on large inventories of observational data and is widely used for meteorological and engineering applications. The vertical profile of temperature vs. altitude and atmospheric pressure is shown in Figure 4, with the temperature being 15°C or 288.15 K at the surface, linearly decreasing in the troposphere with the (minus) gradient being 6.5 K/km up to the altitude of 11 km, and taking a constant value of 216.65 K above this up to 20 km.

FIGURE 4

standard atmosphere: (left) temperature and pressure as functions of altitude; (right) temperature and altitude as functions of pressure.

Clearly, the standard atmosphere is far from isothermal. Here we extensively use the standard atmosphere as satisfactorily representing reality, even though newer data and research () support slightly lower values of Γ, particularly at high latitudes. This will be verified below (Section 3.7). It is interesting to note that the average temperature in the standard atmosphere, calculated (by numerical integration) from the coordinates of Figure 4, iswhere, by the ideal gas law, the air density is proportional to , with denoting air pressure.

Figure 5 shows that temperature gradients appear in other planets and satellites in our solar system, and even in extrasolar planets (for additional rocky exoplanets see also ). In other words, isothermal atmospheres never seem to actually exist in planets.

FIGURE 5

, (Figure 3) with the permission of AIP Publishing. Saturation vapor pressure curves for the indicated condensates (substances forming clouds) for each planet are also plotted (dashed lines), details of which are given in the original publication.

Another important case, additional to the isothermal, is the isentropic atmosphere, where the per-molecule entropy is constant, independent of altitude. Following here we use the adjectives adiabatic and isentropic interchangeably, even though under some conditions, adiabatic processes can differ from isentropic processes. It is well known that in the absence of water vapor (dry conditions) the atmosphere has a constant temperature gradient equal towhere m/s is the gravity acceleration (typical value for the troposphere) and J kg–1 K−1 is the specific heat capacity of the dry air for constant pressure. This is depicted in Figure 6 for the troposphere and lower stratosphere.

FIGURE 6

Another limiting case is the moist (pseudo-)isotropic profile, also shown in Figure 6. To calculate it we have used the following equation (Beers, 1945, p. 359; ) which gives a satisfactory approximation for both the isentropic change of saturated air (in which the liquid water remains in the space considered) and saturated pseudo-adiabatic change (in which the liquid water precipitates immediately):

In this, J kg–1 K−1 is the gas constant of dry air, J kg–1 is the latent heat of vaporization, and and are the saturation water vapor pressure and mixing ratio, respectively, for temperature T. By taking derivatives in Equation 4, setting , and performing algebraic manipulations we readily find as a special case Equation 3 for dry conditions.

The saturation water vapor pressure is given by the Clausius-Clapeyron equation, as reformulated for water vapor by :with T0 and e0 representing the coordinates of the triple point of water. The mixing ratio iswhere ε = 0.622 is the ratio of molar masses of water vapor and dry air.

It is seen in Figure 6 that the moist adiabatic atmosphere does not have a constant temperature gradient. Rather the gradient is low (4.9 K/km) at the surface and becomes equal to the dry adiabatic rate (9.8 K/km) at high altitudes. This is a consequence of the fact that at low temperatures the saturation vapor pressure is close to zero and, as already explained, Equation 4 switches to Equation 3. The standard atmosphere has a profile close, but not identical, to that of the moist adiabatic one. Its slope is constant (6.5 K/km) for the entire troposphere. The differences become large at high altitudes near the tropopause and above.

It is useful to compare these profiles with observational data. To this aim, we enroll in the ERA5 reanalysis over the period 1940 – 2024 and find the global average temperature over this period at the available pressure levels, as well as the global average geopotential heights for these pressure levels. The results, areally averaged over the torrid and the north frigid zones, are shown in Figure 7. To adapt the standard profile to the conditions of Figure 7 we take the following steps: (a) we replace the standard ground temperature of 288.15 K with the average temperature of the zone at 1,000 hPa, (b) we apply the standard gradient Γ = 6.5 K/km up to the altitude where the temperature takes the standard stratospheric value of 216.65 K, and (c) we keep a constant value of 216.65 K above this altitude up to 20 km.

FIGURE 7

We observe that in the tropics the actual average temperature profile is located between those of the moist adiabatic and the standard atmosphere, but above 12.5 km it is closer to the former than the latter. However, in the polar area, observations are totally irrelevant to the moist adiabatic profile. Notice that the moist adiabatic conditions would suggest greater gradients in colder conditions (right panel) than in hotter (left panel) because of the Clausius-Clapeyron equation which results in lower WV quantities for lower temperatures (compare the gradients at the surface level, 4.0 and 7.7 K/km in the two panels). The reality is the exact opposite as seen in Figure 7 (see also Figures 12, 13 below). This means that the moist adiabatic process, despite often being invoked in the literature to explain atmospheric conditions, is not sufficient a tool to describe and explain what happens in reality. All these support the conclusion that the standard atmosphere, despite its empirical basis and its low explanatory power, is more representative for the entire range of conditions on Earth, while the moist adiabatic profile is representative only for the tropics. The dry adiabatic profile is always very different from average conditions. For this reason, the analyses that follow are based on the standard atmosphere.

3.4 Is the ARE responsible for the observed temperature profile in the atmosphere?

To address the question in the heading, we apply the established greenhouse theory, without considering doubts that have been cast on the validity of the theory or alternative hypotheses (e.g., ; ). Rather, we enroll the RRTM software, which calculates at each atmospheric level z the SW and LW radiation going up and down, namely, the four quantities . The algebraic sum of the four quantities for equal to the ground level is the net radiation balance at the surface, which is used in evaporation calculations. While, due to its “rapid” characteristic this software is not the most accurate, we deem it satisfactory for the exploratory character of this research. Taking the differences at two adjacent levels we calculate the quantities , the total radiative energy imbalance , and finally the radiative energy imbalance slope, REIS, based on the following equations.

By convention, quantities directed up are taken as positive and those directed down as negative. If REIS turns out to be positive at a certain layer, between the levels , this means that the layer receives energy of another type, such as sensible or latent thermal energy. This should be equal to REIS, so that the total energy balance be reinstated.

Figure 8 depicts examples of these quantities between the elevations z = 0 and z = 426 m (the fourth altitude given as output by RRTM) and for conditions specified in the figure caption, which include both RAGs and clouds. Two cases are examined, an isothermal atmosphere (Γ = 0) and the standard atmosphere (Γ = 6.5 K/km). In the former case, REIS is almost zero, which means that no forcing appears to drive the atmospheric state away of the equilibrium. Therefore, the atmosphere would remain at the equilibrium (isothermal) state, despite the presence of RAGs and clouds.

FIGURE 8

In contrast, for the standard atmosphere gradient a positive REIS appears, which means that the total energy balance needs to be reinstated by other energy forms. These forms are indeed present and are the sensible heat due to warm air parcels moving from the Earth’s surface up, and the latent heat due to evaporation, again moving from the surface up.

A more complete picture is provided in Figure 9, where REIS is calculated and plotted for the entire range of altitudes modeled by RRTM, for four different cases of ARE (RAG with clouds, RAG without clouds, no NC-RAG with clouds, no NC-RAG without clouds) and for two cases of temperature gradient, Γ = 6.5 K/km (standard atmosphere; upper row of panels) and Γ = 0 (isothermal atmosphere, lower row of panels). What was discussed with respect to Figure 8 is confirmed in Figure 9 for all cases examined and for all altitudes in the troposphere. The isothermal case yields almost zero REIS in the entire troposphere (z < 11 km) in all cases, which means that the presence of RAGs, whether NC or WV, is not enough for the observed gradient to emerge, nor for the ARE to appear.

FIGURE 9

The last observation may sound counterintuitive to many who regard the presence of RAGs as a sufficient reason for the temperature gradient. However, our finding is consistent with experimental evidence reported by , who stated:

the greenhouse effect is mainly the result of a temperature difference over the propagation path of the radiation and thus the lapse rate in the atmosphere.

It further agrees with the following two statements by and , respectively):

It should be emphasized that if dynamical mixing were to have led to an isothermal atmosphere, then there would be no warming due to added greenhouse gases. In the counterfactual case that mixing were to have lead to increasing temperature with altitude, then added greenhouse gases would actually cool the atmosphere. In brief, greenhouse warming depends crucially on the existence and properties of dynamic mixing within the troposphere, and not simply on the radiative picture.

It is an interesting curiosity that had convection produced a uniform temperature, there wouldn’t be a greenhouse effect.

On the other hand, with the standard atmosphere gradient (Γ = 6.5 K/km) there is a positive REIS in the entire troposphere, which leaves room for sensible and latent heat to act, reinstating the total energy balance and stabilizing the temperature gradient.

3.5 What factors can explain the observed temperature profile?

Given the analysis of the previous subsection, it becomes clear that it is not the ARE that creates and maintains the temperature gradient. Rather the mechanisms responsible for the temperature gradient are:

  • • the warming of the soil and liquid water by the sunshine during the day and their cooling during the night;

  • • the water evaporation and transpiration at the surface level and condensation aloft;

  • • the convection, and the implied vertical transfer of sensible and latent heat;

  • • the winds caused by spatial temperature differences and influenced by Coriolis forces.

These are not static forcings, but processes, i.e., perpetual changes in the climatic system. The processes occur on different time scales, some of which are too small to drive the atmosphere to the equilibrium (isothermal) state. It is noted that the processes occur at a macroscopic level, with the motion of masses of air, typically referred to as parcels. And as noted by , because of the very small thermal conductivity of air, it takes a very long time for appreciable heat to flow into or out of a parcel of reasonable size.

The driving mechanisms of these processes are the following:

  • 1. Clouds form and disappear, strongly affecting the SW and LW radiation processes.

  • 2. The Earth’s surface is not homogeneous in terms of radiation absorption and reflection (varying albedo). The massive presence of water in liquid and solid phases with different albedo values, both on the Earth’s surface (hydrosphere, cryosphere) and in the atmosphere (clouds) is responsible for spatial and temporal variations in the albedo, as well as in the thermodynamic properties of different parts of the Earth.

  • 3. Earth is round (not flat) and the sunrays come with different slopes at different places.

  • 4. Earth rotates around its axis on a daily basis.

  • 5. Earth rotates around the Sun on an annual basis.

  • 6. Earth’s orbit around the Sun is elliptical, resulting in changes in the distance between the two bodies.

  • 7. The climatic system is complex and is subject to irregular changes.

A rough quantification of changes produced by these mechanisms is shown in Table 2, where it is seen that some of them are of the order of 100%, reaching the absolute maximum of 200% on a daily scale. This quantification is just indicative and is unable to quantitatively predict the observed temperature profiles as seen in Figure 7. It is also unable to produce the standard temperature gradient of 6.5 K/km, which here we regard as an empirical model with satisfactory accuracy in representing the average conditions over the entire globe. It incorporates all mechanisms discussed here, but it can hardly be inferred from them by deduction. Sometimes the literature connects it to the moist adiabatic profile but the data above (particularly the right panel in Figure 7) do not support this idea.

TABLE 2

#ReasonTime scaleAssumptionsCondition 1Condition 2 (W/m2) (W/m2)% changea
1Change from clear to cloudy skyHoursφ = 45°, α = 0.3, JanuaryClear sky (C = 0)Cloudy sky (C = 1)7425100
2Spatial change of albedobHoursφ = 45°, C = 0, JanuaryLand (α = 0.3)Sea (α = 0.05)7410030
3Change of latitudebHoursα = 0.3, C = 0, Januaryφ = 45°φ = 46°74705
4Change from daylight to nightDailyφ = 45°, α = 0.3, C = 0, JanuaryDaylightNight1970200
5Change in the angle of incidence of the sun’s raysAnnualφ = 45°, α = 0.3, C = 0JanuaryJuly74245107
6Changing distance between Earth and SunAnnualφ = 0°, α = 0.3, C = 0JanuaryJuly2202067
7Recent global warmingDecadesGlobal average conditionsBeforeAfter160.61610.2

Typical changes that cause departure from the isothermal atmosphere and their quantified effect on SW radiation processes. R1 and R2 denote the downwelling solar radiation at the surface level for the specified conditions 1 and 2, respectively; φ, α and C denote the latitude, albedo and cloudiness, respectively. The calculations are based on and .

a

.

b

Spatial changes are translated into temporal changes assuming wind velocity of 10 m/s.

3.6 What is the relative importance of ARE and temperature gradient in climate dynamics?

The information provided in the different panels of Figure 9 includes the surface temperature at each of the cases examined (A1-A4, B1-B4). A summary of the results is given in Figure 10, along with some additional hypothetical cases. The first imaginary-world case is a planet shaped as a disk perpendicular to the Sun’s rays, without an atmosphere and receiving the average radiative energy that the Earth does. In this case, the so-called effective temperature is easy to calculate from the Stefan-Boltzmann law, i.e.,where I = 1360 W/m2 is the solar irradiance with I/4 = 340 W/m2 representing the average solar energy received by (the spherical) Earth’s surface, while σ = 5.67 × 10−8 W m−2 K−4 is the Stefan–Boltzmann constant, ε is the emissivity of the radiating body and α is its albedo. Taking ε = 1 and using the Earth’s current albedo , so that W/m2, the result is 254.5 K.

FIGURE 10

About the same (T = 252 K) would be the temperature in the imaginary-world case where there is an atmosphere but no ARE (no RAG and no clouds; second bar in Figure 10) or even with NC-RAGs, but without WV and clouds and without a temperature gradient. If the gradient of 6.5 K/km is present (third bar in Figure 10) the temperature increases to 262 K. Therefore, the effect of the NC-RAGs is zero for an isothermal atmosphere and 10 K for an atmosphere with temperature gradient of 6.5 K/km. By no means is it close to about 30 K as typically implied in the literature.

In general, in an isothermal atmosphere, the effect of RAGs, including NC-RAGs, WV and clouds, is practically zero, as shown by the continuous (blue) line in Figure 10, which fluctuates only slightly, from 249 to 256 K. In contrast, if there is a gradient of 6.5 K/km, the temperature increases up to ∼288 K in the realistic case A1, in which the atmosphere contains all RAGs and clouds.

Hence, it is the temperature gradient that makes the surface-level temperature increase from about 252 K (which is close to the average 250 K of Equation 2) to about 288 K (i.e., by 36 K). This increase is usually attributed to the “greenhouse effect”, but it is mainly the result of the temperature gradient, whose origin is not the ARE but the processes described in section 3.5. This is another reason why it is necessary to stop using the term “greenhouse effect” when talking about the climatic system.

We stress that, among the various cases shown in Figure 10, only cases A1 and A2 are realistic, while all the others are hypothetical. A final hypothetical case shown at the rightmost end of the figure is where the NC-RAGs are doubled. This results in a temperature increase of zero in comparison to case B1 or 1.5 K in comparison to case A1.

While the temperature gradient is the determinant factor that maintains the Earth’s temperature much higher than the effective temperature of 254.5 K, with the mechanisms causing this gradient being also causes of maintaining warm conditions, the RAGs also play a significant role as can be inferred from Figure 10. An atmosphere free of RAGs (both WV and NC) and clouds would keep the temperature close to the average, i.e. 252 K (second bar in Figure 10). However, this case is equally imaginary as that of no atmosphere (first bar in Figure 10) and does not represent planet Earth. For Earth is massively (71%) covered by water, which evaporates, thus giving rise to WV and clouds in the atmosphere. Even without NC-RAGs but with WV and clouds, the surface temperature would be 277 K (case A3)—a substantial increase with respect to the effective temperature.

While Figure 10 provides information on the importance of different agents of Earth’s warming, it should be stressed that the figure compares the realistic cases (A1, A2) with several imaginary-world cases (all others). Hence, it cannot be a basis for understanding the importance of each factor in real-world conditions. The scientific approach for the latter task is to determine the partial derivatives of a suitable multivariate function representing the ARE, such as the downwelling or outgoing LW radiation, at the point of the current conditions and compare them. Such an analysis has been done in with the following resulting percent contributions: (a) for downwelling LW radiation, WV and clouds 95%, CO2 4%, all other 1%; (b) for outgoing LW radiation, WV and clouds 87%, CO2 5%, all other 8% (primarily due to the influence of the stratospheric ozone).

Water vapor, besides being the determinant RAG and the cause of clouds, also has overwhelmingly more impact than CO2 because its quantity in the atmosphere varies substantially over time and space. This makes WV an agent of perpetual change on all time scales. Relevant changes, mostly on hourly to annual time scales, have already been described in Section 3.5. However, changes also occur at decadal, centennial scales and beyond. These long-term changes are inherent in all geophysical processes, including hydroclimatic ones, and are described in stochastic terms by the so-called Hurst-Kolmogorov dynamics (; ; ; ; ).

As far as long-term changes are concerned, the CERES satellite data show a decline in the albedo of about 0.004 for the entire observation period (post 2000), which translates to 1.4 W/m2 (higher than the average imbalance of 0.4 W/m2 shown in case 7 of Table 2). This is consistent with an observed decline in cloud area fraction (; Appendix A.2, ; Appendix B; see also ; ). Interestingly, a recent study by identifies a record-low planetary albedo in 2023, apparently caused largely by a reduced low-cloud cover in the northern mid-latitudes and tropics, in continuation of a multi-annual trend, all of which contributing to the recent global temperature surge.

Coming back to the results shown in Figure 10, a relevant issue to stress is that the analysis made is unidimensional. Is that representative for the average condition of Earth, which is spherical? A negative reply has been suggested by , who showed that the mean physical temperature of a spherical body is always lower than its effective radiating temperature computed from the globally integrated absorbed solar flux. This happens because of Hölder’s inequality between integrals and the fact that Stefan-Boltzmann law is a power law with an exponent 4.

This suggestion indeed makes sense for the condition of a planet without an atmosphere (first bar in Figure 10), where an integration over the Earth’s sphere would result in a lower value than shown in the figure. However, if we consider the realistic conditions of Earth, in particular cases A1 and A2 of Figure 10, rather than the imaginary-world conditions, the cautionary note by Nikolov and Zeller does not apply. For the relationship between outgoing radiation and temperature is no longer a power law but a virtually linear relationship, as shown in Figure 11, constructed from CERES radiation data and ERA5 Reanalysis near-surface temperature data. Additional information on this linear relationship is provided by , Figure 1 and , Figure 6, which in essence confirm early observations and analyses by . Hence, the quantities depicted in Figure 10 are representative of average conditions for the (spherical) Earth, at least for the realistic cases A1 and A2, albeit not for the imaginary case of an atmosphere-free Earth.

FIGURE 11

3.7 What are the recent changes in the atmospheric temperature gradient?

Once we have recognized that the temperature gradient is a result of macroscopic changes on Earth (Section 3.5) and that it is a critical factor determining the climate (Section 3.6), it is useful to monitor and analyze its changes, which are certainly much more influential than the famed changes in CO2 concentration. To this end, Figure 12 (upper) shows the variation of the atmospheric temperature averaged over the five geographical zones at two pressure levels, 1,000 and 200 hPa, as given by the ERA5 reanalysis.

FIGURE 12

By using the geopotential height at these pressure levels, also available in the ERA5 dataset, we determined and plotted in Figure 12 (lower) the evolution of the vertical temperature gradient per geographical zone, calculated as the (minus) ratio of differences of temperature and geopotential heights between these two levels.

It is clear from Figure 12 that, as we move from the equator to the poles and the surface temperatures decrease, so do the temperature gradients. This is fully compatible with the analyses of Sections 3.4–3.6, according to which higher gradients lead to higher surface temperatures (see also Figure 13, left).

FIGURE 13

Surface temperatures present increasing trends over time in all zones. In contrast, the temperatures at the level of 200 hPa do not show geographical or temporal changes. It is interesting to examine the temporal trends of temperature gradients, visualized in Figure 12 (lower) and summarized in Figure 13 (right). At the torrid zone we observe a falling trend, which can be viewed as negative feedback to the increasing surface temperatures. In contrast, at the frigid zones the trends are increasing, and hence we would expect greater warming. In the temperate zones no notable trends appear.

Indeed, the zonal distribution of the climatic surface temperature differences between the latest and oldest 30-year periods (for which ERA5 Reanalysis data are available in the CLIMEXP platform), seen in Figure 14, confirm these observations. The figure plots show that for the greatest part of Earth, between latitudes 50°S and 50°N there has been a warming trend of about 1.4 K/century, slightly increasing as we move north, as the percentage of land increases. Between 50°N and the Arctic Circle there is a further warming trend, while between 50°S and the Antarctic Circle the trend becomes cooling. The greatest increases of surface temperature are observed in the frigid zones, where we have increasing temperature gradient trends.

FIGURE 14

All the above observations are consistent with the evidence from paleoclimatic studies, which suggests that the surface equator-to-pole temperature gradients were weaker during warm periods, while tropical temperatures changed very little (, p. 17; , p. 289; ; ; ). They are also consistent with the recent study by , which identified wind stilling (compatible with the decrease in the surface equator-to-pole temperature gradient), resulting in weaker evaporation in two-thirds of the ocean and a slight decreasing trend in global-averaged ocean evaporation during 2008–2017. In addition, as already mentioned, the falling multiyear trend in the planetary albedo, largely caused by a reduced low-cloud cover, has been a main driver of the recent global temperature.

4 Discussion and conclusions

The research presented here challenges the ideas that we live in a greenhouse and that science can be settled, and revisits the most fundamental topics related to climate. This is attempted in the simplest possible way and using the fewest premises, such as Newton’s laws, the principle of maximum entropy and the spectroscopic properties of gases. Additionally, the study proposes that commonly used vocabulary should be replaced by rigorous scientific terminology, the main examples being “greenhouse gas” and “greenhouse effect”, which could be replaced by “radiatively active gas” (RAG, comprising water vapor—WV—and non-condensing radiatively active gases—NC-RAGs) and “atmospheric radiative effect” (ARE), respectively.

The conclusions of the analyses presented here can be summarized as follows:

  • • There is an empirically verified ARE in the atmosphere, not only on Earth but also on the other planets. On Earth, ARE is dominated by WV and clouds, with CO2 playing a very minor role—let alone human added CO2 which represents only 4% of the total emissions to the atmosphere.

  • • Equilibrium thermodynamics clearly show (either using the principle of maximum entropy, or stochastic simulation of molecule collisions) that Earth’s atmosphere would be isothermal at the equilibrium, with or without RAGs. In an isothermal atmosphere the temperature would be slightly higher than 250 K, a value which represents the vertically average temperature of the standard atmosphere over the troposphere and stratosphere.

  • • The fact is that the atmosphere is not isothermal. Rather, the troposphere has a vertical temperature gradient of about 6.5 K/km, which is imprinted in the standard atmosphere. The gradient is resultant of macroscopic changes that drive the atmosphere out of equilibrium. While the moist adiabatic changes play a role in shaping this gradient, they cannot fully predict real atmospheric conditions.

  • • The mean surface temperature of 288 K, also imprinted in the standard atmosphere, is much higher than the equilibrium temperature. While RAGs (mostly WV and clouds) play a role in yielding this temperature, the critical factor is the vertical temperature gradient, without which the ARE alone would not be able to increase the equilibrium temperature.

  • • Given the importance of the atmospheric temperature gradient, which is described by large-scale atmospheric thermodynamics, rather than radiative physics, it is puzzling that the emphasis in climate research has been on the latter. This gradient is not a universal constant, but it varies with space and time. It is therefore useful to monitor and analyze its changes. The data show that since 1950 the gradient has weakened in the tropics and grown in the polar areas resulting in decrease of the surface equator-to-pole gradient, as expected in global warming conditions.

A final point worth stressing is that in complex systems, such as the climatic system, observational data are the only scientific test bed for making hypotheses and assessing their validity. Focusing on one of the factors affecting the climatic system, namely, the anthropogenic CO2 emissions, and basing on models that emphasize this factor, may distort our perception of the big picture and be detrimental to science, whose objective is to pursue the truth.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories, namely in the links given in section 2 and the references cited there.

Author contributions

DK: Data curation, Supervision, Conceptualization, Software, Investigation, Methodology, Writing – original draft, Writing – review and editing, Resources, Project administration, Validation, Formal Analysis. GT: Validation, Visualization, Writing – review and editing, Methodology, Formal Analysis, Investigation, Software.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

We are grateful to the colleagues and organizations who have put their huge data sets online along with the data processing and computational systems they have developed. These include the CERES data, the ERA5 Reanalysis, the CLIMEXP and WRIT data and software platforms, and the RRTM software system. We thank two reviewers for the positive assessment of the paper and their detailed and constructive comments. We also thank the Guest Editors and Chief Editors for the processing of the paper and their suggestions, and for approving the paper for publication. We acknowledge useful comments (including points of disagreement) on a preprint of the paper (doi: 10.13140/RG.2.2.14752.49923) which we received from Ned Nikolov, William Happer, Willis Eschenbach, Hermann Harde and Richard Cina. All the above helped us improve the final version of the paper.

In memoriam

Dedicated to the memory of Anna (Annouska) Patrikiou–Koutsoyiannis, who left this world while this research was conducted.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcpxs.2025.1617092/full#supplementary-material

References

Summary

Keywords

climate, climatic system, atmosphere, thermodynamics, greenhouse

Citation

Koutsoyiannis D and Tsakalias G (2025) Unsettling the settled: simple musings on the complex climatic system. Front. Complex Syst. 3:1617092. doi: 10.3389/fcpxs.2025.1617092

Received

23 April 2025

Accepted

23 July 2025

Published

12 August 2025

Volume

3 - 2025

Edited by

Ning Wang, Ministry of Emergency Management, China

Reviewed by

Stavros Alexandris, Agricultural University of Athens, Greece

Patrice Poyet, Independent Researcher, Moorea, French Polynesia

Updates

Copyright

*Correspondence: Demetris Koutsoyiannis,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics