# Simple Climate Models that Predict Climate Change

I have a physics/maths background, and to me the most convincing evidence of a scientific fact is often a very simple model that gives reasonable results. A simplified model, which you can fully understand, is often more convincing than the numerical output of a complex model.

So I'm looking for a simple climate model, which could be used to predict how increasing $\ce{CO_2}$ concentrations would lead to an increase in temperature. For instance, this model need not account for the specific geography of the earth: it might just consider the sun's radiation passing through an atmosphere and striking a flat, homogeneous surface with some constant absorptivity. Then we could just do an energy balance to find equilibrium temperature as a function of $\ce{CO_2}$ concentration.

If I could do a back-of-the-envelope calculation to see that doubling $\ce{CO_2}$ concentrations would be expected to increase the Earth's temperature by a few degrees (I expect nothing more than this order-of-magnitude evidence) then I think that would be convincing evidence that climate change is worth worrying about. You can argue with satellite data and complex models, but you can't argue with basic physics.

So far, the only source I've found that gives such a simplified model is Arrhenius himself, in his 1896 paper. I actually find his model quite convincing, and its prediction ($\sim5^oC$ increase from doubling $\ce{CO_2}$) gives a reasonable prediction too. However, his data on the absorptivity of the atmosphere was (rather quaintly) obtained by early measurements of the intensity of moonlight at different angles in the sky. While this is ingenious, I'm sure with modern technology we could get much better estimates.

So I'm basically looking for either a modernised version of Arrhenius' calculations, or alternatively another simple model which convincingly demonstrates that increasing $\ce{CO_2}$ has the power to alter the earth's temperature by as much as a few degrees.

• My suggestion would be to look at the "climate sensitivity" parameters which are meant to simplify the calculation needed to estimate the impact of a GHG. en.wikipedia.org/wiki/Climate_sensitivity – farrenthorpe May 16 '15 at 3:17
• The thing with "a very simple model that gives reasonable results" for very complex systems, is that you can create two such models that give diametrically opposite answers; a pair of "plausible" models for climate might tell you that there both is, and is not, sufficient grounds to urgently act on greenhouse gases. I'm afraid that this heuristic isn't always the right heuristic, and this is one case where it's not. – 410 gone May 21 '15 at 16:22

It sounds like you are asking for a simple climate model that treats the globe as a homogeneous system. A Zero-dimensional energy balance model may be of interest (see this Penn State site ) which uses a starting point that is also basically described in this post: Temperature as a function of luminosity and greenhouse gas concentrations

If you solve for the temperature of the Earth using a blackbody approximation and thermal equilibrium (with an albedo of about 0.3), you are only taking into account the temperature of the Earth without a greenhouse effect (about 33 °C cooler than actual temperature). It is very difficult to calculate the influence of infrared radiative forcing components (e.g. oxygen and carbon dioxide) to an empty atmosphere, since there are many biogeochemical/physical feedbacks in the climate system which are nonlinear.

However, if you start with current temperature/concentrations, you can add terms that represent particular increases to the GHG concentrations. Those terms are constantly being refined by the scientific community and are summarized in this Wikipedia article about Radiative Forcing (see the link to climate sensitivity at the bottom of the page). This is a kind of "back of the envelope" modeling approach that most climate scientists use to discuss past and future influence of climate forcers, which can be very helpful when discussing the expected climate effects of various emissions/mitigation scenarios.

The Radiative Forcing article in Wikipedia says:

Radiative forcing can be used to estimate a subsequent change in equilibrium surface temperature (ΔTs) arising from that radiative forcing via the equation:

$$\Delta T_s =~ \lambda~\Delta F$$

where λ is the climate sensitivity, usually with units in K/(W/m2), and ΔF is the radiative forcing. A typical value of λ for CO2 is 0.8 K/(W/m2), which gives a warming of 3K for a doubling of CO2 [from pre-industrial times, assuming RF of CO2 is ~ 3.7)].

It's important to note that most climate sensitivity parameters and radiative forcing values are derived with respect to some time period (e.g. 1750 - 2011) and GHG feedbacks are set for a particular climate regime. Revised radiative forcing and climate sensitivity parameters are needed to make calculations for other climate regimes.

The IPCC AR5 figure of Radiative Forcing components attributes 1.65 W/m2 to the actual increase in CO2 concentrations from 1750 to 2011 (or ~1.3K in temperature rise, if lambda is 0.8). The values of Radiative Forcing that IPCC provides are likely close enough for doing realistic scenario testing. However, something in an entirely different feedback regime would need very different parameters. For instance, halving the CO2 in the atmosphere would have a much greater absolute effect on temperature than doubling it; temperature increases logarithmically as CO2 concentrations increase. This is because strong CO2 absorption bands are nearly saturated under present conditions, and increasing CO2 from present conditions mainly affects the flanks of the absorption bands. However, water precipitates out of the system easily under low-temperature conditions, so CO2 has a much greater forcing value in those types of climate regimes.

The most difficult problems with testing simple climate scenarios are related to the oceans and the biosphere, which could have countless effects on the average climate and GHG regime.

The Meteorological Institute of the University of Hamburg, Germany, provides and simple model called "Planet Simulator" (PlaSim).

The Planet Simulator is mainly made for educational purposes and runs on personal computers. However, it is also parallalized and can be run on HPC clusters. There exists an Documentation and the source code is said to be well commented. The whole source code is open source. The programming language is Fortran. One should be able to compile it with gfortran under Linux.

The Planet Simulator offers a GUI (after compiling successfully) to set up da model simulation and run it. There is code for postprocessing included. Thus, there is not need for additional postprocessing software. If you write out the output as netCDF, which is the quasi-standard output format of climate models, it might be practical to install ncview or Panoply because these are common netCDF viewers with which you can look directly into the output files.