The next logical step is to allow spatial resolution into our model. We will do this by breaking the northern hemisphere into 9 latitude regions of 10o each. The general geometry is shown in the following diagram:
In the model we will consider, the solar flux and albedo
both have a latitude dependence (Si and
i, where i runs from 1 to 9 indicating the latitude band).
Also, if a latitudinal band is
colder or warmer than the average global temperature, heat will flow
into/out-of the region. We will assume that this heat flow depends linearly
on the temperature difference between the region and the global temperature,
in other words it is: F×(Ti-TAve).
Putting this additional term into
from the Global Climate Modeling Project
results in an energy balance
|Si||Solar Flux in Latitude Band i|
|This is the product of S/4 (the average global solar constant) times the insolation si|
|The fraction of solar flux incident on each latitudinal band.|
|α i||Albedo in Latitude Band i|
|The albedo of ice is much larger than the albedo of land/water. We can do a crude model of the temperature dependence of the albedo by using α i=0.3 for Ti > Tc or α i=0.6 for Ti <= Tc|
|Tc||Below this Temperature, we assume a Permanent Ice Pack (Tc = -10oC)|
|F||Heat Transport Coefficient (F=3.80 W m-2 oC-1)|
|TAve||Global Weighted Average Temperature|
|This temperature is the weighted average of the temperature in all of the latitude zones on the previous iteration. The weighting factors fi are just the relative fraction of the surface area of the sphere in each latitude zone.|
|A & B||Coefficients expressing Infrared Radiation Loss (A=204 W m-2 and B=2.17 W m-2 oC-1)|
|CE||Heat Capacity (CE = 2.08 × 108 J/m2 oC)|
The Rough Flow of the
iterates Eqn. 5) is:
keyboardcommand to allow you to easily change any variables
MaxTempDiff=1E6- we will be checking for a small value in the next step
MaxTempDiff) is smaller than the tolerance specified by
Run the Matlab script using the parameters given in the program
and this handout.
Try adjusting the factor
(which is multiplied by the Solar Constant) to find out
how sensitive the model is to variations in the solar constant.
By what factor does the solar
constant need to decrease before the earth is completely
glaciated? Before all the permanent ice pack melts?
Modify your program so instead of finding a steady-state solution
it will find a time dependent solution (using Eqn. 4).
|Zone||Annual Mean |
|Solar Insolation |
Solar Constant) b