Energy balance models

"Previous studies identified instabilities for a shrinking ice cover in two types of idealized climate models: (i) annual-mean latitudinally varying diffusive energy balance models (EBMs) and (ii) seasonally varying single-column models (SCMs). The instabilities in these low-order models stand in contrast with results from compre- hensive global climate models (GCMs), which typically do not simulate any such instability. To help bridge the gap between low-order models and GCMs, an idealized model is developed that includes both latitudinal and seasonal variations. The following EBM is a idealized representation of sea ice and climate with seasonal and latitudinal variations in a global domain." (Extraced from: Till J.W Wagner and Ian Eisenman, "How Climate Model Complexity Influences Sea Ice Stability" (2015) p.2 Abstract) For simplicity, I will call the paper "paper" so that I can refer to it more easily on this page.

  ☞ Paper and Documentation

  ☞ Python source code: simple EBM, complex EBM

Temperature EBM

Input Section
\(D^* = Wm^{-2} K^{-1}\) Diffusivity for heat transport:
0.6
\(A = Wm^{-2}\) OLR when \(T = T_m\):
193
\(B = Wm^{-2}K^{-1}\) OLR temperature dependence:
2.1
\(c_w = W yr m^{-2}K^{-1}\) Ocean mixed layer heat capacity:
9.8
\(S_0 = Wm^{-2}\) Insolation at equator:
420
\(S_2 = Wm^{-2}\) Insolation spatial dependence:
240
\(A_0\) Ice-free coalbedo at equator:
0.7
\(A_2\) Ice-free coalbedo spatial dependence:
0.1
\(\alpha_i\) Coalbedo when there is sea ice:
0.4
\(Wm^{-2}\) Radiative forcing:
0
\(\gamma\) Gamma:
1

Aquaplanet EBM with seasonal cycle

"This model is an idealized representation of sea ice and climate with seasonal and latitudinal variations in a global domain. The surface is an aquaplanet with an ocean mixed layer that includes sea ice when conditions are sufficiently cold." (paper p.3; Chapter 2)

"Horizontal diffusion occurs in a ghost layer with heat capacity \(cg\), all other processes occur in the main layer, and the temperature of the ghost layer is relaxed toward the temperature of the main layer with time scale \(tg\). [...] the ghost layer does not represent a separate physical layer such as the atmosphere, which would add physical complexity to the mode." (paper p.7; Chapter 2.e)



\(D^* = Wm^{-2} K^{-1}\) \(A = Wm^{-2}\) \(\alpha_i\)
\(S_0 = Wm^{-2}\) \(B = Wm^{-2}K^{-1}\) \(c_g = W*yr*m^{-2}K^{-1}\)
\( S^{*}_{1} = Wm^{-2} \) \(c_w = W*yr*m^{-2}K^{-1}\) \(\tau_{g}\)
\( S_2 = Wm^{-2} \) \(F_b = Wm^{-2}\)
\(A_0\) \( L_f = W*yr*m^{-3}\) (0-99)
\(A_2\) \(K = Wm^{-1}K^{-1}\) (0-99)

With default parameters the seasonal cycle of the equilibrated climate is plotted.

(a) shows the seasonal cycle of \(E(t, x)\), which fully represents the model state since \(E\) is the only prognostic variable and the forcing varies seasonally. The associated surface temperature (b) and ice thickness (c) are roughly consistent with present- day climate observations in the Northern Hemisphere.

The red curve in Fig. (a) - (c) indicates the ice edge. The blue line in Fig. \(c\) indicates the time of coldest (winter) and the red line the time of warmest (summer). \(x = 0\) represents the latitude at the equator and \(x = 1\) at the North Pole. Coalbedo is the fraction of incident solar radiation \(S\) that is absorbed and not reflected to space (\(1-\alpha\)).

\(T_m\) - melting temperature
OLR - outgoing longwave radiation
\(S_0\) - annual mean insolation at the equator
\(S^{*}_1\) - determines the amplitude of seasonal insolation variations (annual frequency: \(\omega = 2\pi yr^{-1}\) )
\(S^{*}_1\) - determines the equator-to-pole insolation gradient

NOTE: Cloud cover and water vapor are not included in this idealized sea ice model.