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:
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:
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.
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.
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:
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]: