ATM 623: Climate Modeling

Assignment 2: Albedo feedback in the zero-dimensional energy balance model

In lecture we defined a zero-dimensional energy balance model for the global mean surface temperature $T$ as follows

$$C \frac{dT}{dt} =(1-α)Q- OLR(T)$$$$OLR= \sigma (\beta T)^4$$

where we defined these terms:

  • $C = 4\times 10^8$ J m$^{-2}$ K$^{-1}$ is a heat capacity for the atmosphere-ocean column
  • $\alpha$ is the global mean planetary albedo
  • $\sigma = 5.67 \times 10^{-8}$ W m$^{-2}$ K$^{-4}$ is the Stefan-Boltzmann constant
  • $\beta=0.885$ is our parameter for the proportionality between surface temperature and emission temperature
  • $Q = 341.3$ W m$^{-2}$ is the global-mean incoming solar radiation.

For this exercise, we will introduce a new physical process into our model by letting the planetary albedo depend on temperature. The idea is that a warmer planet has less ice and snow at the surface, and thus a lower planetary albedo.

Represent the ice-albedo feedback through the following formula:

$$ \alpha(T) = \left\{ \begin{array}{ccc} \alpha_i & & T \le T_i \\ \alpha_o + (\alpha_i-\alpha_o) \frac{(T-T_o)^2}{(T_i-T_o)^2} & & T_i < T < T_o \\ \alpha_o & & T \ge T_o \end{array} \right\}$$

with the following parameter values:

  • $\alpha_o = 0.289$ is the albedo of a warm, ice-free planet
  • $\alpha_i = 0.7$ is the albedo of a very cold, completely ice-covered planet
  • $T_o = 293$ K is the threshold temperature above which our model assumes the planet is ice-free
  • $T_i = 260$ K is the threshold temperature below which our model assumes the planet is completely ice covered.

For intermediate temperature, this formula gives a smooth variation in albedo with global mean temperature. It is tuned to reproduce the observed albedo $\alpha = 0.299$ for $T = 288$ K.

Your assignment:

Complete the three problems below. Submit your solutions in a single IPython notebook that contains your code, your figures, and your text.

I suggest starting by making a copy of this file, and adding your answers in additional code and text cells following each question.

Submit your completed notebook to me by email before class next Tuesday February 3.

Problem 1

  • Define a Python function that implements the above albedo formula. There is definitely more than one way to do it. It doesn't matter how you do it as long as it works!
  • Use your function to calculate albedos for a wide range on planetary temperature (e.g. from $T=250$ K to $T=300$ K.)
  • Present your results (albedo as a function of global mean temperature) in a nicely labeled figure.

Problem 2

The climate system will be at equilibrium whenever the net shortwave flux into the system balances the longwave flux out of the system, or $ASR = OLR$. Here, you will use a graphical approach to look for equilibrium solutions. On a single graph, plot both $ASR(T)$ and $OLR(T)$ as functions of global mean surface temperature $T$. Use the same temperature range you used in Problem 1.

Intersections of the two graphs indicate temperature(s) at which the energy budget is balanced. According to our model, if the Earth was at this temperature, it should remain at this temperature.

Based on your graphs, try to answer the following questions:

  • How many equilibria do you find in your graph?
  • At approximately what temperatures?
  • If you found more than one equilibrium, try to explain in words why this is possible in the model.

Problem 3

Use a numerical time-stepping method to integrate the full model from an arbitrary initial temperature out to (near) equilibrium. Experiment with different initial temperatures to see what happens over time. Is it possible to get the model to equilibrate in all the different solutions you found in Problem 2?

Present your results in a single graph with temperature on the y axis and time on the x axis. Each curve on the graph should indicate the temperature timeseries from a single model integration. Different curves should show the solutions beginning from different initial temperatures. Use as many different initial temperatures as you think is necessary to fully document the characteristics of the model.

In words, comment on what this exercise has revealed (if anything) about the physics of this system.


In [0]: