Present and Last Glacial Maximum climates as states of maximum entropy production

Present and Last Glacial Maximum climates as states of maximum entropy production


The Earth, like other planets with a relatively thick atmosphere, is not locally in radiative equilibrium and the transport of energy by the geophysical fluids (atmosphere and ocean) plays a fundamental role in determining its climate. Using simple energy-balance models, it was suggested a few decades ago that the meridional energy fluxes might follow a thermodynamic Maximum Entropy Production (MEP) principle. In the present study, we assess the MEP hypothesis in the framework of a minimal climate model based solely on a robust radiative scheme and the MEP principle, with no extra assumptions. Specifically, we show that by choosing an adequate radiative exchange formulation, the Net Exchange Formulation, a rigorous derivation of all the physical parameters can be performed. The MEP principle is also extended to surface energy fluxes, in addition to meridional energy fluxes. The climate model presented here is extremely fast, needs very little empirical data and does not rely on ad hoc parameterizations. We investigate its range of validity by comparing its performances for pre-industrial climate and Last Glacial Maximum climate with corresponding simulations with the IPSL coupled atmosphere-ocean General Circulation Model IPSL_CM4, finding reasonable agreement. Beyond the practical interest of this result for climate modelling, it supports the idea that, to a certain extent, climate can be characterized with macroscale features with no need to compute the underlying microscale dynamics.

1 Introduction

The Earth receives a certain amount of energy from the Sun, in the form of visible light, which it has to radiate back to space, in the form of infrared light, to maintain a steady state. Most Earth System processes, including weather and climate, can be regarded as little more than steps in this process of energy conversion from one form to the other, going through various other forms of energy (potential energy, kinetic energy, heat…). Although in each of these steps the overall quantity of energy has to be conserved, its quality can vary [41, 40]. The field of physics that deals with such energy conversions in different forms and the corresponding evolution of its quality is called thermodynamics, and the quality of energy is measured by a state function called entropy. The constraint of conservation of the exchanged quantity of energy is well-known and commonly used in climate sciences. Less attention is paid to the additional information one can gain by monitoringcreation and exchanges of entropy, either as a diagnostic tool or as a prediction principle relying on an extremum property. This information can be discarded with little consequences for some processes for which physical laws naturally allow for a macroscopic description (and for the subsequent computation of the exchanged energy due to the process), but in some other cases it might bring in a fundamental constraint that should not be ignored.

A typical example is atmospheric heat transport: considering the observational fact that the net top-of-atmosphere radiative budget is positive in the Tropics and negative near the Poles, one has to invoke a poleward atmospheric transport to ensure global energy balance in the steady state (see e.g. [50] for a quantitative discussion based on observations and [28] for a modelling counterpart). An important difficulty comes from the turbulent nature of the laws of atmospheric motion and the lack of a satisfactory theory of turbulence. We are thus forced to integrate numerically the equations of motion at a very high computational cost. On the other hand, the laws of radiation allow for a macroscopic description more easily. To cope with this difficulty, the options up to now have been either to deal with the overdetailed (from the minimalist standpoint we are considering) microscopic equations of motion of the atmosphere, or to use empirical parameterizations.

An alternative was suggested by Paltridge [35, 36] (see [33] for a nice reformulation). He was looking for a large-scale description of heat transport and tested a variational principle stating that atmospheric heat transport adjusts so as to maximize the production of entropy. Several climate studies based on similar models have been published since this seminal work: [13, 11, 54, 25] for zonally-averaged climate, [34, 43] for vertical convection. The Maximum Entropy Production can be interpreted as a maximum energy transport efficiency requirement. The general idea was first expressed by Lorenz [24]. See for instance [40, 23, 52, 26, 27] and references therein for discussions in terms of available potential energy and Carnot engines, and [39, 38] for the competing effect of water vapour removal by the atmosphere.

The purpose of this paper is not to discuss the justification of the Maximum Entropy Production Principle. For a rather comprehensive review of different approaches and applications, the reader is refered to [31] and [20]. In this paper, we present a new formulation of Paltridge’s model: we have developed a different approach for the treatment of radiation, which we believe to be more rigorous and physically sound. In particular we provide a derivation of the analytical expression of the radiative coefficients which are introduced as ad hoc in Paltridge’s model. We also avoid the questionable maximum convection hypothesis, critical in Paltridge’s model [33, 44], by applying the MEP principle on the vertical dimension as well, and discard Paltridge’s cloud cover variable and oceanic transport for the sake of cogency. Formulated this way, the model is as free of tunable coefficient, empirical parameterization and spurious assumption as possible, therefore constituting a clean basis to evaluate the validity of the MEP hypothesis as applied to climate. The only information needed a priori to compute the coefficients are the surface albedo, the vertically-integrated water vapour and ozone density and the concentration. In particular, we are now able to run a sensitivity experiment with respect to the surface albedo parameter. Here, we investigate the effect of the surface albedo distribution corresponding to the presence of large ice-sheets in the Northern Hemisphere during the Last Glacial Maximum. The results for pre-industrial and LGM climates are compared with standard simulations with the IPSL_CM4 atmosphere-ocean general circulation model using similar boundary conditions.

2 Model description

We use a purely thermodynamic two dimensional (horizontal) energy balance model. Only energy transfers are represented. The exchanges of mass, momentum, angular momentum, etc are discarded here. Each grid cell consists of two layers: one atmospheric layer (at temperature ) exchanging energy with the surface, space and adjacents cells in the form of longwave and shortwave radiation, surface heat flux, and atmospheric heat transport, and one surface layer (at temperature ), only exchanging energy with the overlying atmospheric layer and with space, in the form of shortwave and longwave radiation and surface sensible heat (see figure 1). We note the area of cell , . All the simulations presented in this paper are carried out on a regular grid with and for easier comparison with the GCM results using the same grid.

2.1 Net-exchange formulation for longwave radiative transfer

Radiative exchanges are described in the Net Exchange Formulation (NEF) formalism - originally introduced by [14] (see [3, 7] for details). In the NEF, radiative fluxes are replaced by energy exchange rates between two layers as basic objects describing the radiative transfer. Instead of solving the radiative transfer equation, we compute the net energy exchanged per time unit between layers and (resp. layer and boundary surface or boundary surfaces and ) at frequency :


where represents the volume of layer (resp. the area of boundary surface ), the set of optical paths from one point of layer to one point of layer , is if is a boundary surface and the absorption coefficient otherwise, is the transmitivity of the path and finally is the Planck function. Depending on the nature of (surface or volume), is either an infinitesimal surface or volume element. Clearly .

The NEF was originally used in engineering rather than atmospheric sciences. In the few existing cases of atmospheric transfer studies, the main point is that NEF allows for a fine analysis of the greenhouse effect and how it evolves when one parameter changes: it is clear from the matrix which exchanges prevail in the atmosphere (generally, these are the exchanges between neighbouring layers, between the surface and any layer and the cooling to space of any layer and the surface). This description contains more information on the radiative transfer than the traditional flux approach insofar as one not only knows how much energy is radiated through a surface but also where the radiation originates from (computationally, for a layer atmosphere, NEF requires numbers while the flux approach only needs ). Other major advantages of this formulation are the handiness it provides in coupling radiative exchanges to other modes of energy transfer, and the fact that it automatically enforces physical laws such as energy conservation or the direction of energy exchanges depending on the temperature gradient, whatever numerical approximations are made. Clearly these points make it suitable for simple thermodynamic models such as ours.

In the very simple case of our two level cell, we only have to compute the infrared exchange between the surface and the atmosphere , the surface and space and the atmosphere and space . Applying the diffusive approximation ( [4], justified in [22]), angular and vertical integration ( being the height of the atmosphere) yields, for instance, for :


From now on, let us consider the exchanged energy per surface area, keeping the same symbols . Integrating over the whole spectrum, the equations for the exchanged energy can be recast into the following form:


where represents the emissivity of the atmosphere and . Of course the transmission functions of the atmosphere depend on the vertical profiles of absorbing gases, pressure and temperature. In order to preserve the simplicity of the analytical expressions and to allow for comparison with Paltridge, we impose here the heavy constraint that the coefficients are constant in each cell. Given a characteristic temperature for the local atmosphere and a humidity profile, we can evaluate these coefficients using Goody’s statistical model [12] and the spectral data from [45]. The temperature and humidity profiles used here come from a linear interpolation of McClatchey’s reference atmosphere profiles [32]. Hence we prescribe here the coefficients as functions of latitude only. In fact we do not need the full profiles here but only the vertically integrated values. Sensitivity tests were carried out to assess the dependence of the coefficient on the water vapour amount and characteristic temperature and its impact on the MEP climate (see section 3.3).

2.2 Shortwave radiation parameterization

At latitude , the annual mean insolation is given by


where is the solar constant, the eccentricity of the Earth’s orbit and the obliquity. The incoming top-of-atmosphere solar radiation in each cell is given by the latitudinal mean of over the cell, which we write , thus representing the annual mean of the cosine of the zenith angle of the Sun for the latitude zone (weighed by a distance factor) and is therefore a function of the orbital parameters. In all the simulations considered here, we use the present-day values for the orbital parameters: and . The Last Glacial Maximum orbital parameters do not differ much from those values.

Solar radiation exchange rates are parameterized as in [21]. With their notations, let us define


where is the surface albedo, (resp. ) the column-integrated water vapour (resp. ozone) concentration. is the magnification factor accounting for the slant path, (resp. ) are numeric functions corresponding to water-vapour (resp. ozone) absorption, and and account for Rayleigh scattering in the atmosphere (see [21] or [46] for details). Incoming solar radiation at surface and solar radiation absorbed in the atmosphere are then expressed in the simple form:


The coefficient can be interpreted as the proportion of solar radiation directly absorbed in the atmosphere while accounts for absorption after reflection at surface and represents the proportion of solar radiation absorbed by the ground.

To compute the value of the three coefficients , we use humidity and ozone profiles from McClatchey [32] again (as for the longwave coefficients, this is not a critical parameter here), with pressure scaling for water vapour. Surface albedo is prescribed from IPSL simulations (see section 2.4).

2.3 Entropy Production Maximization

For each grid box, we can compute the parameters as explained above and radiative energy exchange rates are given by equations 3 to 11 as functions of these parameters and temperatures and .

Since we are looking for steady states of the system, we impose the following energy balance constraints in each box:


where is the surface heat flux and the horizontal convergence (see figure 1). Formally these equations can be solved to express and as functions of the two unknown variables, and .

Similarly to the energy budget, we can write an entropy budget, separating the rate of change of the entropy of a parcel of the system into two contributions:


where stands for the entropy exchanged with the surroundings per unit time, while is the entropy created per unit time. The second law of thermodynamics states that . In the steady-state, we have .

Assuming local thermodynamic equilibrium, the variation of the entropy for a parcel of the system exchanging an amount of heat is given by . Hence, for a single cell of the model, the entropy production rate (per unit surface area) associated to the horizontal convergence is given by and the entropy production associated to the surface heat flux is .

We thus consider the total material entropy production rate:


Let us search for the energy fluxes distribution maximizing this function subject to the following global constraint: the total heat transport divergence over the globe must vanish, i.e.


We enforce this physical condition by introducing a Lagrange multiplier in the previous equation.


and we call Maximum Entropy Production (MEP) state the solution of the system (12)-(13) maximizing the entropy production rate (15) with the imposed constraint (16). In practice, solving the system reduces to finding the maximum of a function of two variables in each box. One could equivalently solve the optimization problem in terms of the temperature variables, replacing and from equations 12-13 into .

In the MEP procedure, it is only the entropy production rate due to turbulent processes which is maximized, motivated by the supposed connection of the MEP hypothesis with principles of maximum statistical entropy. The entropy production rate due to radiative processes, however, can provide interesting information as a diagnostic tool (see section 3.5).

2.4 IPSL simulations description and albedo forcing

The IPSL model is a fully coupled Ocean-Atmosphere General Circulation Model (OAGCM). The version used here is IPSL_CM4 [30] i.e. the version which has been used for the CMIP3/IPCC AR4 exercise. The model is composed of the LMDZ atmosphere model [16] run at resolution 96x72x19 (longitude-latitude-vertical levels), the ORCA2 ocean model [29], the LIM2 see-ice model [8, 9], all coupled via the OASIS coupler [51].

The present work uses two simulations performed with this model : (1) a pre-industrial simulation, for which the model is forced with constant, pre-industrial greenhouse gas concentration ( ppm, ppb, ppb), insolation, coastlines, topography, land-ice extent and (2) a pseudo LGM run. Since the simple climate model described above does not account for topography, nor for variations, we have chosen to use an OAGCM simulation in which the extent of the LGM ice-sheets is taken into account but not their height. The ice-sheet extent is prescribed from the Peltier ICE-5G reconstruction [42]. Greenhouse gases atmospheric concentrations are maintained at their pre-industrial level, as well as orbital parameters and topography. Note that none of the simulations uses an aerosol forcing, which explains that the simulated pre-industrial climate is slightly warmer than the standard IPCC simulation.

For each of these simulations, the annual mean surface albedo is computed (figure 2) and used as an input for the corresponding MEP simulation.

Hence the boundary counditions for the IPSL_CM4 simulations correspond as much as possible to those for the MEP simulations.

3 Results and discussion

3.1 Present day climate

Figure 3 shows the surface temperature at MEP state for present albedo forcing and a comparison with a control simulation from the IPSL model.

The global mean temperature simulated by our simple MEP model is C, which is C warmer than the IPSL model results. The reason for this strong bias is probably the lack of any clouds in the model. Another reason could be the crude representation of the vertical structure of the atmosphere in the model and the subsequent flaws in vertical radiative transfer calculations. Comparison with the IPSL model reveals that some of the regions of largest discrepancies are regions of high elevation (Antarctica, Tibetan plateau, and to a smaller extent Greenland, cf Fig. 3). Indeed, topography effects are not included in our model. We show in section 3.2 that the combined effects of clouds and elevation account for the bulk of the global difference with the IPSL model.

In spite of this global bias, considering the simplicity of the model and the absence of any degree of freedom to adjust the results, the model gives a strikingly accurate view of the global picture of the climate. Some important specific features, like the presence of the main deserts (Sahara, Arabian peninsula, Kalahari, Australia,..), are already visible taking into account only the albedo effect. Since the only non zonal parameter in the model is the albedo, and its value is roughly uniform over the oceans, we were expecting the surface temperature to be zonally homogeneous over the oceans, which is indeed the case. The positive anomaly of temperature over mid-latitudes regions (compared to the IPSL results, Fig. 3, right) is certainly imputable to the important cloud cover over those areas, lacking in the model.

Figure 4 provides the curves of zonally averaged atmospheric heat transport and net top of atmosphere radiative budget. The net top of atmosphere radiation budget switches from positive in the Tropics to negative near the Poles at roughly 40S and 35N, while the overall shape of atmospheric heat transport qualitatively fits the IPSL model results. It peaks approximatively at 35S and 40N with respective values of 3.8 PW and 4.4 PW. The location of these maximum values is consistent with [48] but the absolute value is slightly underestimated here: [48] gives the peak transport around 6 PW. The absence of any oceanic heat transport in our model prevents us to push the comparison forward. Comparison with the atmospheric transport from the IPSL simulation reveals that the MEP transport is indeed a little low, especially in the Southern Hemisphere.

It might also explain why surface energy fluxes are underestimated here: since all the heat flowing laterally is assumed to come from the atmosphere, the need for surface heat flux is reduced. Indeed, in the MEP state, the global mean surface flux is only (figure 5), compared to [49]. The maximum value of the surface energy flux is reached in the Tropics while it is negative in the polar regions. Comparison with the IPSL simulation (not shown) confirms that the largest discrepancy arises over the oceans, where the lack of oceanic transport in our MEP model is crucial.

These results can also be interpreted in a different manner, considering that the surface energy flux computed by MEP is actually closer to the sensible heat component rather than the sum of the sensible and latent heat flux. Formally, the physical nature of the surface flux is not specified in the equations of the model. Yet, we believe that since it appears in the entropy production rate as the product between the flux and the gradient of the inverse of temperature, and given that the sensible heat flux also depends directly on the temperature gradient at the surface but not the latent heat flux, it makes more sense to compare it to the sensible heat flux (figure 5). In that case, the global average surface flux is approximately twice the value given by [49]. The fact that the global mean average surface flux lies in between the value of the global mean sensible heat flux and the sum of the sensible and latent heat flux can be seen as a trade-off between the need for a large surface flux to ensure the energy balance and the need for a large enough temperature gradient from the entropy production rate point of view.

As discussed in [19], since a MEP state is supposed to be the most efficient heat conducting state given the constraints, one should expect more gentle equator-pole temperature gradients in MEP states than with more conventional models. Here, adjusting the overall surface temperature field by setting to eliminate the global bias, we observe that the Tropics are roughly at the same temperature in the IPSL model and the adjusted MEP model (actually slightly colder in the MEP model) while the mid-latitudes are substantially warmer (around 5 degrees) in the MEP state. However the temperature anomaly is not symmetric for the poles: the South Pole is colder by around C in the adjusted MEP model whereas the North Pole is warmer by around C. Therefore the equator-pole gradient is indeed reduced in the MEP state for the Northern Hemisphere, but not for the Southern Hemisphere, as expected from the relatively low poleward heat transport in the Southern Hemisphere (figure 4).

3.2 A representation of clouds in the MEP climate

As suggested in section 3.1, the major part of the difference between the MEP pre-industrial climate and the corresponding IPSL simulation is likely to be due to the absence of clouds in the MEP model. Inclusion of fixed or variable clouds in MEP models is possible, but unavoidably requires parameterization and extra assumptions, which is not the philosophy of this study. Nevertheless, to evaluate the effect of clouds on the MEP model presented here, it is possible to modify the radiation equations (10)-(11) to include an ad hoc cloud parameter :


The parameter is computed in each grid box from the cloud radiative forcing at the top of the atmosphere in the corresponding area in the IPSL model. On the global mean, the cloud radiative forcing at the top of the atmosphere in the IPSL model is roughly (in accordance with [18]), which corresponds to a mean value of , but the geographical distribution of the radiative forcing is very far from being uniform.

For pre-industrial climate, this parameterization of clouds in the MEP model reduces the temperature by 4.7C on global average. Thus a 2.5C difference remains with the corresponding IPSL run. The temperature anomaly map (not shown) reveals that this difference is mainly due to the effect of topography (Antarctica and Tibet are in particular very warm in the MEP model), not accounted for in the MEP model. As an approximate correction, we apply a constant lapse rate, using the elevation data from the IPSL run output. Figure 6 shows the combined effect of clouds and topography on the MEP pre-industrial climate: the global mean difference with the IPSL pre-industrial run becomes 0.9C, which is comparable to the uncertainty due to the water vapour content of the atmosphere (see section 3.3). Thus a crude parameterization of the main components lacking in the rigorous version of the MEP model suffices to recover the correct global picture of present climate. Note in particular that the warm bands in the mid-latitudes in figure 3 have completely disappeared with the inclusion of clouds.

3.3 Effect of the radiative parameters on the MEP state

As mentioned above, the only parameters in the model are the radiative parameters and the surface albedo (forced in each box with the surface albedo from the IPSL run). Computation of the radiative parameters only requires the value of the column-integrated water vapour density , and a characteristic temperature for .

In this study, we have fixed once and for all the value of the radiative parameters as a function of latitude only, using a linear interpolation of standard atmospheric profiles for and . To assess the sensitivity of the model with respect to this choice, we have also computed the MEP state for pre-industrial surface albedo using a constant value for , independent of latitude. Three choices were made, corresponding to the three standard profiles: Sub-Arctic, Mid-Latitude and Tropical. Table 1 shows the global mean temperature resulting from these uniform choices.

The values for the column-integrated water vapour density and characteristic temperature span a realistic interval: as a comparison, in the IPSL pre-industrial run, ranges from around near the poles to around at the equator. Thus we can estimate the sensitivity of the model to realistic choices in and to be on the order of 1K.

3.4 Last Glacial Maximum

Figure 7 shows the surface temperature difference between the Last Glacial Maximum and pre-industrial control run, at MEP and with the IPSL model. The global mean difference is K colder for the LGM at MEP, compared to K colder for the IPSL model. This remarkable overall agreement means that the sensitivity of the MEP model to surface albedo is comparable to that of the GCM. The major part of the cooling occurs over Canada and Northern Europe for both models, but it is much stronger in the MEP model: between 20 and 30 K for these areas, compared to K for the IPSL model. The cooling is also overestimated in Patagonia - by a factor 2 approximatively - and Tibet - by a factor 3. These extreme values are compensated for by a small underestimate of the cooling almost everywhere else: around K for most of the oceans, South America, Africa, Australia, Antarctica and Southern and Eastern Asia vs to K for the IPSL model. This milder cooling becomes even more dramatic over Russia - K for MEP vs K for the IPSL model - and Greenland - K vs K. Surprisingly the MEP model predicts an intense warming - K - in the Baffin Bay. We suspect this warming to be due to seasonality effects of the albedo.

Beyond the general relative agreement, the MEP model response to surface albedo is a little too strong over land (Canada, Northern Europe, Patagonia and Tibet) but not over sea: the Barents sea is slightly cooler in the IPSL-LGM than in the MEP-LGM.

Figure 8 shows the meridional heat transport for the Last Glacial Maximum in MEP state. There is no major change of regime compared to present climate. Poleward transport is slightly more important (roughly more energy at the maximum) during the LGM, as can be expected from colder high latitudes. This enhancement of the MEP-LGM poleward transport, larger than that obtained by the IPSL simulation, can probably be attributed to the exaggerated cooling of the Northern high latitudes in MEP-LGM. During the Last Glacial Maximum, surface fluxes (not shown) become negative over ice-covered areas, the remaining of the globe showing little difference with present surface fluxes at MEP state.

The effect of clouds for the LGM at MEP state, with the cloud parameter being calculated from local cloud radiative forcing in the IPSL LGM simulation, is less than for pre-industrial climate. On the global mean, clouds cool the LGM climate by C in the MEP model. As a result, the sensitivity of the MEP model with clouds is reduced to C, leading to a sensitivity difference with the IPSL model of C, as compared to C for the model without clouds. Topography has no influence on the LGM-PI surface temperature difference as it is the same for both periods in both models.

Since the radiative coefficients are set once and for all in the MEP model, there is no representation of the water-vapour feedback, contrary to the IPSL model. In the IPSL model, the globally averaged water-vapour density changes from for pre-industrial climate to at the LGM, corresponding to a variation in of around . From section 3.3, one can estimate that the effect of water-vapour feedback on the MEP-LGM climate must be on the order of a few tenth of degrees.

3.5 Entropy budgets

We can decompose the total entropy production rate of section 2.3 into contributions from various processes as follows:


where are the entropy production rates due to, respectively, absorption of solar radiation in the atmosphere, absorption of solar radiation at the surface and surface-atmosphere longwave interaction:


The entropy production rate of a thermodynamic process can be interpreted as a measure of its irreversibility. Therefore, the inspection of the relative values of the different terms above, gives valuable insight with respect to the role played by the different processes in driving the Earth out-of-equilibrium.

The expressions for the contributions to entropy production involving interaction between radiation and matter are only approximate, as noticed by [5], but we shall not delve further into this here, both for simplicity and to compare with published results. The reader is referred to [47] or [53], for example, for computations taking into account the correct entropy of the radiation field. Nevertheless, the nonlocality of the NEF is already an improvement, in principle, as compared to the usual flux formulation, which is a local expression of radiative transfer [6].

Table 2 shows the entropy budget of the Earth under pre-industrial and Last Glacial Maximum conditions at the MEP state. These values are compared with estimations for present climate from observations [41], from the intermediate complexity model PLASIM [10] and from the Hadley Centre GCM [37]. In both cases, the entropy production due to absorption of solar radiation by the atmosphere is a little lower than other estimates. More importantly, the entropy production due to absorption of solar radiation at the surface is significantly higher than other estimates: this is due to the high solar absorption at the surface (in the absence of clouds) which is not compensated for by a warmer surface. The entropy production due to longwave interaction between the surface and the atmosphere (the greenhouse effect) is comparable to the values for PLASIM and HadCM3, and the material entropy production is considerably lower than other estimates. This last point is explained by the fact that latent heat fluxes are not explicitly represented in the MEP model, while they account for the major part of the material entropy production in all the other studies. Note that the total entropy production at MEP state is slightly lower for the LGM as compared to the pre-industrial, in spite of a higher material entropy production.

4 Conclusion

The aim of this study was to provide a physically reliable, more consistent and more easily generalized version of the Paltridge model. The radiative part has been entirely reformulated using the Net-Exchange Formulation in which the physical meaning of all the coefficients is clearly identified. Therefore the calculation of these coefficients, to the desired degree of approximation, only involves standard radiative data. Besides, all the energy fluxes other than radiation (horizontal convergence and surface fluxes) are subjected to the same principle: Maximum Entropy Production.

We have shown here that a thermodynamic model based only on a simple but robust radiative scheme and the principle of maximum entropy production yields results comparable to those of the IPSL model, for surface temperature, for both the pre-industrial and the Last Glacial Maximum climates, with a negligible computational cost. Obviously, the description adopted in the MEP model discards an important number of climate features. The inclusion of some key processes, like the seasonal or the hydrological cycles, in this thermodynamic framework would be of great interest. But the model in its present form already serves the purpose of providing a compromise between a minimal model with very little a priori data and information needed, and acceptable efficiency.

In particular, one of the strongest points in the model is certainly that it does not include any adjustable parameter. This ability to get rid of the usual, varyingly important, parameter calibration makes it a good candidate to investigate climates where little is known or where some phenomena are likely to be different from the usual parameterization validity range, like for instance climates of other planets, inside or outside the solar system, or paleoclimates.

Further developments also include coupling the model with a state-of-the-art radiative code to come up with a full three dimensional model, which would allow for more realistic paleoclimate simulations.

The basic idea of variational thermodynamic principles such as the principle of Maximum Entropy Production, which still remains to be proved, is that microscopic details are irrelevant, to a certain extent and depending on the question we are trying to answer, to macroscopic behaviour. The general process of not making any unjustified assumption on those microscopic details finds a nice formulation at equilibrium with the principle of Maximum Entropy and the usual apparatus of equilibrium thermodynamics can be derived easily from the application of this principle to the classical ensembles of statistical mechanics [17]. But in spite of the conceptual similarity, the link between Maximum Entropy Production Principle and Maximum Entropy Principle still is not clear (see [1, 2] for an attempt to derive the former from the latter and [15] for a comment). This study expresses the view that regardless of the veracity of the principle in full generality, it can certainly be useful in application in the field of climate sciences, as well as in many others.


The authors would like to thank Jean-Louis Dufresne for fruitful discussions and advice, as well as Gilles Ramstein and Didier Roche for their stimulating remarks and constant encouragements.

Profile () T (K) (K)
Sub-Artic 1.5 272 0.432 295.6
Mid-Latitude 2.3 283 0.436 295.9
Tropical 5.1 300. 0.444 296.5
Table 1: Column-integrated water vapour density (), characteristic temperature (T) and resulting longwave atmospheric emissivity () for the three standard atmospheric profiles. The last column () gives the global mean temperature in the MEP state for pre-industrial surface albedo, with uniform atmospheric emissivity.
Peixoto et al PLASIM HadCM3 MEP Pre-industrial MEP LGM
258 255 216 217
561 557 812 667 650
24 7 11 7 7
- 28 39 - -
32 29 38 16 18
17 7 12 0 0
892 883 912 906 892
Table 2: Comparison of the entropy budget of the Earth computed from observations ([41], some values have been updated after [20]) and two models of different complexity (PLASIM: [10], HadCM3: [37]), in . designates all the other contributions to the total entropy production rate, including eventual imbalance.
Figure 1: Energy exchanges in one model box : incoming solar energy (thin solid lines) absorbed by atmosphere and ground , longwave radiative exchanges (dashed) ; exchange between the surface and the atmosphere , cooling to space of the surface and of the atmosphere ) and large scale heat transport convergence () and surface sensible heat () (thick solid arrows).
Figure 2: Surface albedo forcing for the MEP model, computed from the output of the IPSL_CM4 simulations. Left : annual mean surface albedo values from the pre-industrial simulation with the IPSL model. Right : annual mean surface albedo values from the LGM ice extent simulation with the IPSL model. The surface albedo is higher over continents than over the open ocean, and it is higher in high-latitude areas than in low-latitude areas because of the presence of snow and ice. The main difference between pre-industrial and LGM surface albedo distributions is the presence of a large ice-sheet over North America and Scandinavia during the LGM.
Figure 3: Results for surface temperature. Left: Surface temperature at MEP state for present conditions, in C. Right: Surface temperature difference between the MEP state and the IPSL simulation. Contour lines interval is C, positive contours are drawn in solid lines, negative contours in dashed lines and the null contour as a dotted line.
Figure 4: Left: Total Meridional energy transport predicted by MEP for present conditions (solid) and for the IPSL model (dashed), in PW. Right: Net radiative budget at the top of atmosphere as a function of latitude (zonally averaged) for the MEP state (solid) and for the IPSL model (dashed) at present conditions, in .
Figure 5: Left: Surface heat flux predicted by MEP for present conditions, in . The mean value is . Right: Difference between the surface energy flux at MEP state and the surface sensible heat flux for the IPSL simulation for present conditions. Contour lines interval is , positive contours are drawn in solid lines, negative contours in dashed lines and the null contour as a dotted line.
Figure 6: Pre-industrial surface temperature difference between the MEP model with clouds, including a topography correction, and the IPSL model. Contour lines represent the -5 isoline (dashed), the null isoline (dotted) and the +5 isoline (solid). This figure is to be compared with figure 3, right.
Figure 7: Left: Surface temperature difference between Last Glacial Maximum conditions and present conditions at MEP states. Right: Surface temperature difference between Last Glacial Maximum conditions and present conditions for the IPSL model. Contour lines space is C, positive contours are drawn in solid lines, negative contours in dashed lines and the null contour as a dotted line. The global mean LGM cooling is C for the MEP model and C for the IPSL model.
Figure 8: Total Meridional energy transport predicted by MEP for LGM conditions (solid line) and present conditions (dashed line) at MEP state (left) and for the IPSL simulation (right), in PW.


  1. R. Dewar. Information theory explanation of the fluctuation theorem, maximum entropy production and self-organized criticality in non-equilibrium stationary states. Journal of Physics A: Mathematical and General, 36:631–641, 2003.
  2. R. Dewar. Maximum entropy production and the fluctuation theorem. Journal of Physics A: Mathematical and General, 38:371–381, 2005.
  3. J.-L. Dufresne, R. Fournier, C. Hourdin, and F. Hourdin. Net exchange reformulation of radiative transfer in the m band on Mars. J. Atmos. Sci, 62:3303–3319, 2005.
  4. W. M. Elsasser. Heat Transfer by Infrared Radiation in the Atmosphere., volume 6 of Harvard Meteorological Studies. Harvard University Press, 1942.
  5. C. Essex. Radiation and the irreversible thermodynamics of climate. J. Atmos. Sci, 41(12):1985–1991, 1984.
  6. C. Essex. Radiation and the violation of bilinearity in the thermodynamics of irreversible processes. Planetary and Space Science, 32(8):1035–1043, 1984.
  7. V. Eymet, J.-L. Dufresne, P. Ricchiazzi, R. Fournier, and S. Blanco. Long-wave radiative analysis of cloudy scattering atmospheres using a net exchange formulation. Atmospheric Research, 72:239–261, 2004.
  8. T. Fichefet and M. Maqueda. Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics. Journal of Geophysical Research, 102:12609, 1997.
  9. T. Fichefet and M. Maqueda. Modelling the influence of snow accumulation and snow-ice formation on the seasonal cycle of the antarctic sea-ice cover. Climate Dynamics, 15:251–268, 1999.
  10. K. Fraedrich and F. Lunkeit. Diagnosing the entropy budget of a climate model. Tellus A, 60(5):921–931, 2008.
  11. J. Gerard, D. Delcourt, and L. Francois. The maximum entropy production principle in climate models: application to the faint young sun paradox. Q. J. R. Meteorol. Soc., 116:1123–1132, 1990.
  12. R. Goody. A statistical model for water vapour absorption. Q. J. R. Meteorol. Soc., 78:165–169, 1952.
  13. H. Grassl. The climate at maximum entropy production by meridional atmospheric and oceanic heat fluxes. Q. J. R. Meteorol. Soc., 107:153–166, 1981.
  14. J. Green. Division of radiative streams into internal transfer and cooling to space. Q. J. R. Meteorol. Soc., 93:371–372, 1967.
  15. G. Grinstein and R. Linsker. Comments on a derivation and application of the ‘maximum entropy production’ principle. Journal of Physics A: Mathematical and Theoretical, 40:9717–9720, 2007.
  16. F. Hourdin, I. Musat, S. Bony, P. Braconnot, F. Codron, J.-L. Dufresne, L. Fairhead, M.-A. Filiberti, P. Friedlingstein, J.-Y. Grandpeix, G. Krinner, P. LeVan, Z.-X. Li, and F. Lott. The LMDZ4 general circulation model: climate performance and sensitivity to parametrized physics with emphasis on tropical convection. Climate Dynamics, 27:787–813, 2006.
  17. E. Jaynes. Information theory and statistical mechanics. Physical Review, 106:620–630, 1957.
  18. J. Kiehl and K. Trenberth. Earth’s annual global mean energy budget. Bull. Amer. Meteor. Soc, 78:197–208, 1997.
  19. A. Kleidon, K. Fraedrich, T. Kunz, and F. Lunkeit. The atmospheric circulation and states of maximum entropy production. Geophys. Res. Lett, 30(23):2223, 2003.
  20. A. Kleidon and R. Lorenz, editors. Non-equilibrium Thermodynamics and the Production of Entropy: Life, Earth, and Beyond. Springer, Berlin, 2005.
  21. A. Lacis and J. Hansen. A parameterization for the absorption of solar radiation in the Earth’s atmosphere. J. Atmos. Sci, 31:118–133, 1974.
  22. J. Li. Gaussian quadrature and its application to infrared radiation. J. Atmos. Sci, 57:753–765, 2000.
  23. C. Lin. An extremal principle for a one-dimensional climate model. Geophys. Res. Lett, 9:716–718, 1982.
  24. E. N. Lorenz. Generation of available potential energy and the intensity of the general circulation. Dynamics of Climate, pages 86–92, 1960.
  25. R. Lorenz, J. Lunine, P. Withers, and C. McKay. Titan, Mars and Earth: Entropy production by latitudinal heat transport. Geophys. Res. Lett, 28(3):415–418, 2001.
  26. V. Lucarini. Thermodynamic efficiency and entropy production in the climate system. Phys. Rev. E, 80:021118, 2009.
  27. V. Lucarini, K. Fraedrich, and F. Lunkeit. Thermodynamic analysis of snowball Earth hysteresis experiment: Efficiency, entropy production and irreversibility. Q. J. R. Meteorol. Soc., 136:2–11, 2010.
  28. V. Lucarini and F. Ragone. Energetics of PCMDI/CMIP3 climate models: Net energy balance and meridional enthalpy transports. Rev. Geophys., 49:RG1001, 2011.
  29. G. Madec, P. Delecluse, M. Imbard, and C. Lévy. OPA version 8.1 Ocean general circulation model reference manual, 3. Technical report, LODYC, 1997.
  30. O. Marti, P. Braconnot, J.-L. Dufresne, J. Bellier, R. Benshila, S. Bony, P. Brockmann, P. Cadule, A. Caubel, F. Codron, N. D. Noblet-Decoudré, S. Denvil, L. Fairhead, T. Fichefet, M.-A. Foujols, P. Friedlingstein, H. Goosse, J.-Y. Grandpeix, E. Guilyardi, F. Hourdin, A. Idelkadi, M. Kageyama, G. Krinner, C. Lévy, G. Madec, J. Mignot, I. Musat, D. Swingedouw, and C. Talandier. Key features of the IPSL ocean atmosphere model and its sensitivity to atmospheric resolution. Climate Dynamics, 34:1–26, 2010.
  31. L. Martyushev and V. Seleznev. Maximum entropy production principle in physics, chemistry and biology. Physics Reports, 426:1–45, 2006.
  32. R. McClatchey, J. Selby, F. Volz, R. Fenn, and J. Garing. Optical properties of the atmosphere. Air Force Cambridge Res. Lab., 1972.
  33. D. O’Brien and G. Stephens. Entropy and climate. ii: Simple models. Q. J. R. Meteorol. Soc., 121:1773–1796, 1995.
  34. H. Ozawa and A. Ohmura. Thermodynamics of a global-mean state of the atmosphere—a state of maximum entropy increase. J. Climate, 10:441–445, 1997.
  35. G. Paltridge. Global dynamics and climate-a system of minimum entropy exchange. Q. J. R. Meteorol. Soc., 101:475–484, 1975.
  36. G. Paltridge. The steady-state format of global climate. Q. J. R. Meteorol. Soc., 104:927–945, 1978.
  37. S. Pascale, J. Gregory, M. Ambaum, and R. Tailleux. Climate entropy budget of the hadcm3 atmosphere–ocean general circulation model and of famous, its low-resolution version. Climate Dynamics, pages 1–18, 2009.
  38. O. Pauluis and I. Held. Entropy budget of an atmosphere in radiative–convective equilibrium. part i: Maximum work and frictional dissipation. J. Atmos. Sci, 59(2):125–139, 2002.
  39. O. Pauluis and I. Held. Entropy budget of an atmosphere in radiative–convective equilibrium. part ii: Latent heat transport and moist processes. J. Atmos. Sci, 59(2):140–149, 2002.
  40. J. P. Peixoto and A. H. Oort. Physics of Climate. Springer: New-York, 1992.
  41. J. P. Peixoto, A. H. Oort, M. D. Almeida, and A. Tome. Entropy budget of the atmosphere. J. Geophys. Res, 96:10981–10988, 1991.
  42. W. Peltier. Global glacial isostasy and the surface of the ice-age Earth: the ICE-5G (VM2) model and GRACE. Annual Reviews of Earth and Planetary Science, 32:111–149, 2004.
  43. T. Pujol and J. Fort. States of maximum entropy production in a one-dimensional vertical model with convective adjustment. Tellus, 54:363–369, 2002.
  44. T. Pujol and J. Llebot. Extremal climatic states simulated by a 2-dimensional model: Part i: Sensitivity of the model and present state. Tellus A, 52:422–439, 2000.
  45. C. Rodgers and C. Walshaw. The computation of infra-red cooling rate in planetary atmospheres. Q. J. R. Meteorol. Soc., 92:57–92, 1966.
  46. G. Stephens. The parameterization of radiation for numerical weather prediction and climate models. Mon. Wea. Rev, 112:826–867, 1984.
  47. G. Stephens and D. O’Brien. Entropy and climate. i: Erbe observations of the entropy production of the earth. Q. J. R. Meteorol. Soc., 119:121–152, 1993.
  48. K. Trenberth and J. Caron. Estimates of meridional atmosphere and ocean heat transports. J. Climate, 14:3433–3443, 2001.
  49. K. Trenberth, J. Fasullo, and J. Kiehl. Earth’s global energy budget. Bull. Amer. Meteor. Soc, 90(3):311–323, 2009.
  50. K. Trenberth and A. Solomon. The global heat balance: Heat transports in the atmosphere and ocean. Climate Dynamics, 10:107–134, 1994.
  51. S. Valcke. OASIS3 user guide (prism_2-5). Technical report, PRISM support initiative, 2006.
  52. A. Vos and P. Wel. The efficiency of the conversion of solar energy into wind energy by means of Hadley cells. Theoretical and Applied Climatology, 46:193–202, 1993.
  53. W. Wu and Y. Liu. Radiation entropy flux and entropy production of the earth system. Rev. Geophys., 48:RG2003, 2010.
  54. P. Wyant, A. Mongroo, and S. Hameed. Determination of the heat-transport coefficient in energy-balance climate models by extremization of entropy production. J. Atmos. Sci, 45(2):189–193, 1988.
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Comments 0
Request comment
The feedback must be of minumum 40 characters
Add comment
Loading ...

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description