Coriolis Forces in Hopkins and Simmons Vorticity Equation
In the thread vector operations in Hopkins and Simmons, I compute the components of the curl as:
In my post Coriolis forces in Hopkins and Simmons I compute the coriolis force as:
And the partial derivatives are given by:
I’ll derive the rest of this later but this doesn’t seem to be the form of the prognostic equation used by Hopkins and Simmons.
Divergence Free Flow
In the post Vector Operations in Hopkins and Simmons I derived the divergence operator as follows:
If the divergence of the velocity equals zero then:
Which implies:
Lagrangian Mechanics and The Heat Equation
I have been noticing a lot of discussion on Lucia’s blog about an energy Balance Climate Model.
Two Box Model: Algebra for ver. 1 test.
Two Box Model: Now assuming surface temperature are ‘mixed’ values.
Two Box Model: Rough idea how to obtain parameters.
Two Box Models & The 2nd Law of Thermodynamics
Which started based on a post by Tamino, called “Two Box”
The biggest criticism I’ve scene about these models is that there is very little physics which is used to form these models and consequently they won’t capture much of the nonlinear dynamics.
The basic premises of the model is that the time lags are associated with the heat capacity of the ocean. Consequently, the ocean and some of the atmosphere is partitioned into boxes and the model tries to determine how these boxes are coupled.
The end result though is the time legs end up being associated with eignvectors which, project onto both boxes. Consequently the modes of the system do not coincide with the boxes. Moreover, it is desirable to minimize the number of model parameters which need to be fit and therefore principles of physics, should be used to derive a more realistic “Two box model” or should I say more generally “Two mode model” as the modes do not coincide with the boxes.
The most simple equation with regards to the transfer of energy is the heat equation:
1)
The heat equation is essentially a diffusion equation based on Brownian Motion. But we can adjust the constant to try and account for other heat exchange processes.
In one dimension the heat equation can be written as:
2)
In this equation the time derivative only depends on the spatial derivatives. Keeping a mind on the dynamics, the higher the frequency that the climate forcing is, the more shallow it will penetrate into the ocean because the more ocean you include the greater the thermal inertia (heat capacity) and consequently the greater the damping. Thus a good fit for the temperature response at a given frequency as a function of depth, may be an exponential function of the depth.
Thus, from this observation and because the mathematics are fairly simple consider a set of basis functions which the special component is exponentially decreasing from the surface. That is let:
3)
Plugging this into the heat equation one gets:
4)
Where:
The constant represents how quickly heat diffuses in the ocean at a given depth. Equating terms on the heat equation gives:
5)
Now let there be a function , and let the quanity being diffused in the heat equation be the temperature. Then the energy constraint is given by:
6)
This constraint makes one of the generalized coordinates redundant.
Differentiating (5) with respect to time, gives:
7)
is determined by the forcing. Now some notes, with regards to the coordinate system:
1)I suggest taking the as the redundant basis function, the one that has the smallest value .
2)Since this model is a linearization (or approximation), I suggest using as the temperature variable at a given depth, the amount by which it exceeds the average temperature at that depth. That is model the temperature anomaly instead of the actual temperature.
3)Given two then we are looking at the change in the heat transfer from the mean instead of the actual heat transfer.
4)These are only suggestions.
Now the Lagrangian, is essentially a cost function. In our cost function, we want to minimize the error in the derivatives (the square of equation (5) integrated over y) and as an additional constraint, I want to minimize the change in entropy with respect to y.
When a parcel of water rises because it is hotter then the surrounding and therefore less dense it will expand with constant entropy (adiabatic expansion) if it exchanges no heat with the surroundings. Therefore, the heat induced ocean currents will act in order to try and minimize the entropy change between adjacent regions. (Note the actual effect of this process is to increase entropy).
From thermodynamics:
8 )
Therefore:
9)
from which we can get:
10)
Substituting into (10) the expression for temperature (equation (3)) equation (10) gives:
11)
The entropy cost function is:
12)
The Lagrangian is given by:
13)
Where is a stiffness parameter for the entropy.
Equation 13) is subject to the constraints for some (These are the energy constraints mentioned above.
14)
15)
Now the Euler Lagrange Equation is used to obtain the dynamics (differential equations) from the Lagrangian
16)
This will give a second order differential equation and it might be necessary to use some linear algebra to rearrange the equation. You can convert this into a first order differential equation by using the Hamiltonian form. If only three basis functions are used then one gets “two mode model”. There are no restrictions on the number of bais functions used.
Coriolis Forces
A derivation for coriolis forces can be found be found on the wikipedia page for fictitious forces.
In general for an accelerating reference frame in rectangular coordinates the factious forces are given by:
Where:
he first term is the Coriolis force, the second term is the centrifugal force, and the third term is the Euler force. When the rate of rotation doesn’t change, as is typically the case for a planet, the Euler force is zero.
Looking specifically at the Coriolis force:
which gives in (eastwest,northsourth, height) coordinates:
Where is the latitudinal coordinate (equator=zero latitude).
In general the cross product for a coordinate with orthonormal direction vectors is given by:
http://en.wikipedia.org/wiki/Coriolis_effect#Formula
since the basis direction vectors are orthogonal in hopkins and simmons coordinates write write:
(note with regards to weather the system is right handed we can choose the direction of the logitudanal cordinate to make it right handed.)
Just to recall from the post (Hoskins and Simmons (1974) Coordinate System):
where theta is the latitude.
Where is the surface pressure and is the vertical coordinate.
is the longitude.
Additionally:
U is the longitudinal component of the velocity
V is the latitudinal component of the velocity
W is the vertical component of the veolicty (not used in Hopkins and Simmons 1974)
Now the angular velocity of the earth in Hopkins and Simmons is given by:
Where the sign of is positive for the northern hemisphere and negative for the southern hemisphere.
Therefore:
Some comments:
The result obtained is essentialy the same result that one would get if, they took the(east west,north south, altitude) coordinate system and replaced with .
The only differences are the order and sign of the components. These are the only differences because both coordinate contain the same unit vectors. In my example of a Hopkins and Simon’s like coordinate system I used a different order for the components then was used in my example for the (eastwest, north south altitude) coordinate system. This will effect the sign in the cross product.
I wrote the z component of the angular velocity as to emphasize that the positive direction for the z component in Simpons coordinate system is downward. However, the actual angular rotation of the earth in simons coordinate system still have a postive component depending on the which direction is defined as positive for the longitudinal coordinate.
The order which we specify the coordinates determines the right handedness of the coordinate system. Therefore, righthandedness is not inherently a geometric property because it depends on the order of the coordinates. For instance, in standard Cartesian coordinates
In our case the first coordinat,e , was specified in the downward direction, our second coordinate, , points south, now using the right hand rule means that gives the positive direction for the third coordinate in the east direction.
It is for these reasons that differences can arrise, and therefore it is very important when doing cross products to clearly express the postive direction of the coordinate unit vectors and the order of the coordinates.
Hoskins and Simmons (1974) Coordinate System
sHoskins and Simmons (1974) Coordinate System
Hoskin and Simon used the following coordinate system:
where theta is the latitude.
Where is the surface pressure and is the verticle coordinate.
is the longitude.
This coordinate system is also used in the following model which was recomended to me as a simple model:
As the complexity of general circulation models has been and still is growing considerably, it is not surprising that, for both education and research, models simpler than those comprehensive GCMs at the cutting edge of the development, are becoming more and more attractive. These medium complexity models do not simply enhance the climate model hierarchy. They support understanding atmospheric or climate phenomena by simplifying the system gradually to reveal
the key mechanisms. They also provide an ideal tool kit for students to be educated and to teach themselves, gaining practice in model building or modeling. Our aim is to provide such a model of intermediate complexity for the university environment: the PlanetSimulator. It
can be used for training the next GCM developers, to support scientists to understand climate processes, and to do fundamental research.
http://www.mi.unihamburg.de/216.0.html?&L=1
In order to use the coordinates given by Hoskins and Simmons in the Navier Stokes equations it is necissary to find the operations for Div Grad and Curl in this coordinate system. You can get a free PDF here:
http://www.scribd.com/doc/2590597/LecturesonTransformationofCoordinates
Which shows how to derive these operations for any coordinate system. These derivations express the apove operations in terms of scaling quantities. The scaling quantitys are the magnitude of the derivative of the change in postion, with respect to a change in one of the coordinates in the new coordinate system.
In our case for the coordinate the scalling quanity is:
because when you change latitude the distance moved is proportional to the radius multiplied by the change in latitudinal (measured in radians).
Now expressing the latitude in terms of
Now for the which is the longitude coordinate, depends on how you define zero latitude. The distance changed for a given change in longitude is the product of the change in logitude (measured in radians) multiplied by the distance to the earths axis of rotation. If zero lattitude is the center of the earth then this distance is if zero latitude is the north pole (typical in spherical coordinates) then this distance is .
We will see that taking zero latitude as the north pole will be simpler.
But Therefore:
The vertical coordinate is the ratio of the hydrostatic pressure to the hydrostatic pressure at sea level.
In Hoskins & Simmons the hydrostatic pressure used in the coordinate ignores the lapse rate. In some models variations on this are used. For instance in a model used by Judith Curry this coordinate is somewhat adjusted based on the surface topology of the earth. (see my post: hybird Coordinate System)
In Curries model this is done to try and keep wind flow on the horizontal plane. Anyway, the equation for hydrostatic equilibrium is:
Now:
Hmmmmmmm: From what I’ve done the only way I can see to make the constants disappear is to set
Here is what Hoskins and Simmons have to say about temperature scaling and such.
all nondimensionalized using as length scale the radius of the planet, a, ; as time scale the reciprocal of its angular velocity ; temperature scale (R as the gas constant) and pressure scale Bar.
In http://www.mi.unihamburg.de/fileadmin/files/forschung/theomet/planet_simulator/downloads/PS_ReferenceGuide.pdf was taken to be 250K. Some numbers of rellevence:
Gas constant of air: 287.058 J kg−1 K−1
and for the molar mass of air:
Standard temperature =273.16k
Source for the molar mass of air:
Standard pressure=1 Atm=1013.250*10^2 Pa.
The molar mass of air at standard temperature and pressure is 28.964*10^3 kg/mole. Molar
http://www.blurtit.com/q139778.html
Standard Gravity: 9.80665 m/s2
A Comparison of GCMs and Some History (Two Links)
I found two links which seem to be a good place to start with understanding the differences between the current GCM, and the history of GCM.
The GCMReality Intercomparison
Project for SPARC (GRIPS):
Scientific Issues and Initial Results
S. Pawson,a,b,c K. Kodera,d K. Hamilton,e T. G. Shepherd,f S. R. Beagley,g B. A. Boville,h
J. D. Farrara,i T. D. A. Fairlie,j A. Kitoh,d W. A. Lahoz,k U. Langematz,c E. Manzini,l
D. H. Rind,m A. A. Scaife,n K. Shibata,e P. Simon,o R. Swinbank,n L. Takacs,p
R. J. Wilson,e J. A. AlSaadi,q M. Amodei,o M. Chiba,r L. Coy,p J. de Grandpré,g
R. S. Eckman,q M. Fiorino,s,t W. L. Grose,q H. Koide,d J. N. Koshyk,f D. Li,k,n
J. Lerner,m J. D. Mahlman,e N. A. McFarlane,u C. R. Mechoso,i A. Molod,p
A. O’Neill,k R. B. Pierce,q W. J. Randel,h R. B. Rood,b and F. Wuh
ABSTRACT
To investigate the effects of the middle atmosphere on climate, the World Climate Research Programme is supporting the project “Stratospheric Processes and their Role in Climate” (SPARC). A entral theme of SPARC, to examine model simulations of the coupled troposphere–middle atmosphere system, is being erformed through the initiative called GRIPS (GCMReality Intercomparison Project for SPARC). In this paper, an overview of the objectives of GRIPS is given. Initial activities include an assessment of the performance of middle atmosphere climate models, and preliminary results from this evaluation are presented here. It is shown that although all 13 models evaluated represent most major features of the mean atmospheric state, there are deficiencies in the magnitude and location of the features, which cannot easily be traced to the formulation (resolution or the parameterizations included) of the models. Most models show a cold bias in all locations, apart from the tropical tropopause region where they can be either too warm or too cold. The strengths and locations of the major jets are often misrepresented in the models. Looking at threedimensional fields reveals, for some models, more severe deficiencies in the magnitude and positioning of the dominant structures (such as the Aleutian high in the stratosphere), although undersampling might explain some of these differences from observations. All the models have shortcomings in their simulations of the presentday climate, which might limit the accuracy of predictions of the climate response to ozone change and other anomalous forcing.
http://www.gfdl.noaa.gov/bibliography/related_files/pawson0001.pdf
196575: Spread of AGCMs
 Modeling Groups Proliferate
 Modeling Innovations
 Research on Carbon Dioxide and Climate
 Early Climate Politics and AGCMs
 References
By 1965, three groups in the United States had established ongoing efforts in general circulation modeling:
 Geophysical Fluid Dynamics Laboratory
 UCLA Dept. of Meteorology
 National Center for Atmospheric Research
At this point, AGCMs and modeling techniques began to spread by a variety of means. Commonly, new modeling groups began with some version of another group’s model. Some new groups were started by postdocs or graduate students from one of the three original AGCM groups. Others built new models from scratch. An AGCM family tree offers a visual map of these relationships, with links to some of the models and modeling groups.
http://www.aip.org/history/sloan/gcm/1965_75.html
I plan to have more to say on this later.

Recent
 Laplace Transform Via Limits
 log(CO2) and Scary Graphs
 Numeric Solutions to The Heat Equation
 Coriolis Forces in Hopkins and Simmons Vorticity Equation
 The Cross Product in Non Orthogonal Coordinate Systems
 Lagrangian Mechanics and The Heat Equation
 Laplace Transform of f(t) Related to smoothed f(t)?
 Coriolis Forces
 Vector Operations in Hoskins and Simmons Coordinates
 API/Object Viewers/Memory Mapping/
 Defining a Microsoft access Datasource
 Fractal Modeling of Turbulence

Links

Archives
 July 2012 (1)
 September 2009 (5)
 August 2009 (19)
 March 2009 (2)

Categories

RSS
Entries RSS
Comments RSS