Chapter7 rota frms atmosphere

7-1Chapter VII. Rotating Coordinate Systems 7.1. Frames of References In order to really look at particle dynamics in the context of the atmosphere, we must now deal with the fact that we live and observe the weather in a non-inertial reference frame . Speci ally, we will look at a rotating coordinate system and introduce the Coriolis and centrifugal force . Rotating and Non -rotating Frames of Reference First, start by recognizing that presence of 2 coordinate systems when dealing with problems related to the earth: a) one fixed to the earth that rotates and is thus accelerating (non -inertial), our real life frame of reference b) one fixed with respect to the remote "star" , i.e., an inertial frame where the Newton's laws are valid. Apparent Force In order to apply Newton's laws in our earth reference frame, we must take into account the acceleration of our coordinate system (earth). This leads to the apparent forces that get added to Newton's 2nd law of motion: (inertial)realF dV dtm=rr (7.1) (non-inertial)apparent realF F dV dtmm=+rrr (7.2) where apparentF mr are due solely to the fact that we operate / observe in a non -inertial reference frame. We need to relate to inertialnoninertialdAdA dtdt−rr . 7-27.2. Time Derivatives in Fixed and Rotating Coord inates Reading Section 7.1 of Symon. Using concepts developed in the last chapter, let's introduce some notion that will prove to be useful. Ωr z y TVr rr Latitude circle a=colatitude latitudef= x Equatorial Plane Let rr = position vector from the origin (center of the earth). Then, we can find out that the instantaneous tangential velocity of a particle at any distance rr ( az+rr where a = mean radius of earth and z = altitude above ground) is given by TVr =Ω×rrr . (7.1) Here Ωr is the vector angular velocity of the earth rotation. We can see this by first looking at the magnitude of TVr . According to the above figure, ||||||||sin||TVRrr a =Ω=Ω=Ω×rrrr rr (7.2) Rr 7-3Now, let's try to relate dA dtr in 2 reference frames: a) Fixed (absolute) system ˆˆˆ(,,)ijk b) Rotating ( Ωr ) system ˆ ˆˆ(',',')ijk z z' y' Ωr y x x' Now, a vector Ar is the same vector no matter what coordinate system it is viewed from. This ˆˆ ˆˆˆˆ '''''' xyzxyzAAiAjAkAiAjAk=++=++r . Given that the primed system is rotating, howev er, the time derivative of Ar will be different if viewed from the 2 systems. Mathematically, we have 'a' for absolute ˆ ˆˆ ' ' 'ˆ ˆˆ''' ˆ ˆˆ''''''y ax z y x z xyzdA dAdA dAijkdtdtdtdt dA dA dAijkdtdtdt didjdkAAAdtdtdt=++ =++ +++r (7.3) Note that the primed unit vectors vary with time! (Remember our earlier derivation of the velocity in plane polar coordinates?) 7-4 ˆ'()itt+Δ ˆ'iΔ due to rotation of unit vector ˆ'()it Recall earlier that TVr =Ω×rrr for a particle rotating at angular velocity Ωr . Since Tdr Vdt=rr therefore dr dtr = rΩ×rr (7.4) The above is correct only for a vector whose length does not change. The same is t rue for the unit vectors, therefore ˆ'ˆ' ˆ'ˆ' ˆ'ˆ'diidt djjdt dkkdt=Ω× =Ω× =Ω×r r r (7.5) If we now loot at dA dtr in the reference frame of the rotating (primed) system, then ˆ ˆˆ',' and 'ijk appear to be fixed! ' ' 'ˆ ˆˆ'''y x rzdA dA dAdAijkdtdtdtdt=++r (7.6) which is a time rate of change in Ar seen from the rotating reference frame and is analogous to adA dtr earlier. Thus, from (7.3), we can write 7-5adA dtr = ˆ ˆˆ '(')'(')'(')r xyydAAiAjAkdt+Ω×+Ω×+Ω×rrrr =ˆ ˆˆ ('''''')rr xyydAdAAiAjAkAdtdt+Ω×++=+Ω×rrr rr à adA dtr =rdAAdt=+Ω×rrr (7.7) Let's look at a physical interpretation . If Ar=rr, then we can write adr dtr rdrrdt=+Ω×rrr (7.8) In another word, absrelcoordVVV=+rrr . Think of the earlier example of throwing a baseball towards a person on the ground from a moving railway car. If you can remember this simple analogy, you will be able to remember the formula in the box. Absolute velocity in the inertial frame aV Velocity relative to moving (rotating) coordinates ()rV Velocity of the coordinate system itself (remember Vr=Ω×rr r) 7-67.3 Equation of Motion in Absolute and Rotating Coordinates Let's now apply concepts from the previous section to any arbitrary velocity vector of a particle assuming that the origins of the rotatin g and fixed (inertial) frames are the same. In this case, the position vector to any location is the same in both systems. Our goal is to find the equations of motion in the absolute and rotating coordinates. Now, we just showed that aVVr =+Ω×rrrr (7.9) where Vr is the relative velocity (we shall drop subscript r from now on). Apply Eq.(7.7) to aVr in (7.9) à 2adV dVrVrdtdt dVdrVrdtdt dVddrrVrdtdtdt dVVrdt=+Ω×+Ω×+Ω×  =+Ω×+Ω×+Ω×Ω× Ω =+×+Ω×+Ω×+Ω×Ω×  =+Ω×+Ω×Ω×rrrrrr rr rrrrrrrr rr rrrrrrrr rrrrr r (7.10) Apply the Newton's second law à 2iaiF dV dVVrdtdtm =+Ω×+Ω×Ω×=∑rr rrrrr r (7.11) The part in the box is the equation of motion in the rotating coordinate system! It describes the change of (relative) velocity in time subjecting the net force. The forces on the right hand side are real forces, and t he second and third term on the left arises because of the coordinate rotation, and there are apparent (not real) forces. We will discuss them in more details in the following. 7-7We will look at the last term one the left with respect to an earth oriented coordinate system: Ωr z y TVr 2R−Ωr rr Latitude circle a=colatitude latitudef= x Equatorial Plane Note that ||||cosRr f =rr. Look at the last term in (7.11): r Ω×Ω×rr r. SHOW IT FOR YOUSELF that 2rR Ω×Ω×=−Ωrrr r! (7.12) (Hint – use the right ha nd rule to find out the direction of this vector, then use the definition of cross product to determine its magnitude) Using this we can rewrite equation (7.11) as 22iaiF dV dVVRdtdtm=+Ω×−Ω=∑rr rrrr (7.13) Abs. Accel Rel. Accel. Coriolis ac cel Centripetal accel Net real force This is a VERY IMPORTANT EQUATION . Rr 7-8How do we make use of this equation? Let's look at the forces on the right hand side of (7.13) in the real atmosphere. They are Gravity = gr (forc e per unit mass) PGF (per unit mass) = 1pr−∇ (you have seen this before and you will derive this term in Dynamics I) Friction = …. (neglected for now) These are real forces (those seen and felt in an inertial / fixed reference f rame), so let's put them into the equation: 212dVVRgpdt r+Ω×−Ω=−∇rrrr r (7.14) this is the acceleration that we see and measure in our earth coordinate system – not an inertial frame, but the next 2 terms take that into account. Rearranging (7. 14) à Coriolis force Centrifugal force 212dVpgVRdt r=−∇+−Ω×+Ωrrrr r (7.15) Equation (7.15) is "THE" EQUATION OF MOTION used in meteorology – it is the backbone of meteorology! Relative acceleration Real forces Apparent forces due solely to the rotation of the coordinate system 7-97.4 Forces in Equation of Motion We are now g oing to spend more time looking in details at each of these terms. 7.4.1. Gravity Let's look first at gravity . Recall that ˆeGMm gr a=r (7.16) where Me = mass of earth and a = mean radius of the earth. It turns out it is convenient to combine andgR2Ωr r, the centrifugal force. Due to the presence of the centrifugal force that's directed away from the axis of rotation, the earth is actually elliptic, with the radius at the equator slightly larger than that at the poles. If one measures that force that's pulling a mass downward, it is the sum of true gravity (the one due to gravitational effect) and the centrifigual force, i.e., 2 netggR =+Ωrrr . (7.17) Note that gr points t owards the center of the earth while netgr is slightly equatorward. Since netgr is what we really measure at the surface of the earth, from now on we drop the subscript net, and assume gr includes both effects. gr points in the direction of plumbline and is perpendicular to a static water surface – the horizontal plane. Equation (7.15) becomes: 12dVpgVdt r=−∇+−Ω×rrr r (7.17a) here gr is the netgr in (7.17). See figure below. 7-10 Read Symon Sections 7.2 and 7.3, Riegal Pages 144 -153! 7.4.2. Coriolis Force Let's now look at the other apparent force – the Coriolis force. In contrast to the centrifugal force, the Coriolis force is acting only when the object moves relative to the rotating / non -inertial coordinate system. Consider the case of a merry -go-around. 2 If person #1 rolls a ball toward 2 on a frictionless surface at a un iform speed in a straight line, it will appear to deflect to the right (in this case) – or in a direction opposite the rotation. To a person on the ground, the ball appears to travel in a straight line (non -1 7-11accelerating) – but from the coordinate system of the merry -go-around, the ball accelerates (velocity changes). Let's look physically at the mathematical expression for the Coriolis force : 2 corFV =−Ω×rrr (7.18) (1) The minus sign is important, just like in advection. Don't forget it! (2) The Coriolis force acts perpendicular to the direction of the velocity vector and in direct proportion to its amplitude: Ωr Vr VΩ×rr - VΩ×rr In the northern hemisphere , where Ωr is counterclockwise, the Coriolis force acts to the right of the velocity vector , i.e., it tends to deflect air parcel to the right of their dire ction of motion (see merry -go-around example). (3) The Coriolis force changes only the direction of the particle's motion, not it's speed, because the force is always perpendicular to the direction of motion, just like the centripetal force in uniform c ircular motion. Components of Coriolis Force Now, because the Coriolis force is a vector, we need to look at its components. It turns out, as we will see, that it affects motion in all three directions: x y zduCotherdt dvCotherdt dwCotherdt=+ =+ =+ 7-12How do we obtain those three components? One is purely mathematical – we just need to fine the components of 2 corFV =−Ω×rrr in a chosen Cartesian coordinate system. We can also obtain them in a more physical way. Let's do it in the physical way first. Consider fi rst the situation where we have a particle of unit mass moving freely on the frictionless surface of the rotating earth if the particle is initially at rest in the rotating frame , then only the centrifugal force is acting, along with the gravitational forc e. We combined the sum of these two into a new gravity that is perpendicular to the local horizontal. Suppose now that the particle is suddenly set in motion toward the east . It is now rotating faster than the earth, and thus the centrifugal force will b e stronger . If Ω = magnitude of angular velocity of the earth, Rr the position vector from the axis of rotation , as shown earlier, and u the eastward speed relative to the ground (i.e., the speed relative to the rotating coordinate system), then the total centrifugal force is that from earth's rotation plus that due to the eastward motion, i.e., 2 2 2 22 uuuRRRRRRRΩ Ω+=Ω++rrrr (7.19) For synoptic -scale motion, ||||uuRorRΩΩ== , the last term can be negle cted to a first approximation. How do we show this? 2 222212uuuuuRRRRRRRRR ΩΩΩ+=+ Ωrrrr ; Thus, the only term left is the second one, which is the Coriolis force (recall that it's due to the relative motion only) due to motion along a latitude circle . Let's look at the component of this vector! Original rotation Due to relative motion Centrifugal force due t o earth rotation as before. It is added to gr Deflecting force that acts outward, in the direction of Rr 7-13 Ωr ˆ2cos()uk f Ω Rr f 2ˆ 2uRuRRΩ=Ωr ˆ 2sin()uj f −Ω f We see that, due to the Coriolis force alone, 2sin() 2cos()cor cordvudt dwudtf f=−Ω  =Ω  (7.20) Thus, if a particle is moving eastward in the horizontal plane in the Northern Hemisphere, it is deflecting Southward and Upward! The opposite is true for a particle moving from east to west. If you plug in numbers, you will find that the vertical component of the Coriolis force is much smaller than gravity , and thus it usually does not cause big change in the vertical position (height) of an object, but the horizontal components can be significant, compared to other horizontal forces (horizontal pressure gradient force is the other dominant force). Now, let's repeat this exercise for a particle displaced toward the equator . What's different here? Now, R, the distance from the particle from the axis of earth rotation changes (increases). 7-14 Ωr North pole R1 start R 2 stop Because the forces acting on the particle are central, the torque is zero – there is no torque in the east -west direction, therefore, the particle conserves angular momentum . So, if R increases, what will happen to the velocity? Initially, it has no motion tangential to the direction of rotation. The angular momentum is Ω R2, so as R increases, the Ω of the particle has to decrease à the absolute eastward velocity decreases , the particle starts to deflect toward west. So, the particle moving south is deflected to the west. Let's try to derive mathematically the speed at which the deflection occurs. Let δR be the change in R and δu the change in zonal wind speed. Note that δφ <0 for an equator -ward displacement but δR>0. δu is to be determined. 22 startstopddRRdtdtqq=  At the start, the particle has no speed in the east -west direction, so it's A.M. is Ω R2 7-15where Ω is earth's angular velocity. At the ending poi nt, it has acquired a relative angular velocity of u RRd d+, so the final A.M. is 22()()()uRRRRuRRRRddddddΩ++=Ω++++ . Equaling the A.M at start and ending point à 2222()() RRRRuRRR dddd Ω+Ω+Ω++=Ω Consider small displacements therefore small δu and δR, the 2 δ terms can be neglected (since they are even smaller), and solve for δu, we get 2 uRdd=−Ω (7.21) Now, we can rewrite δR in terms of δφ as follows: R = a cos( φ) à δR = - a sin( φ) δφ Thus, δu = 2 Ω a δφ sin(φ) à 2sin() corduad dtdtff=Ω  Furthermore, a d φ = dy, therefore 2sin() corduvdtf=Ω  (7.22) (7.22) gives the acceleration (deflection) in east -west direction due to the earth -relative motion in the north -south direction. For northward m otion, v>0 à du/dt >0 à eastward deflection, deflection to the right. If you repeat the same analysis for a particle launched vertically, then the conservation of absolute angular moment shows that you will see a 7-162cos() corduwdtf=−Ω  (7.23) i.e., the particle will accelerate toward the east. Adding this to the result just obtained gives the component of the Coriolis force in the zonal direction! 2sin()2cos() corduvwdtff=Ω−Ω  (7.24) Having derived the components of the Coliolis force, 2 V −Ω×rr , the equation of motion, (7.15) can now be written as (neglecting friction): 12sin()2cos()dupvwdtxffr∂=−+Ω−Ω∂ 12sin()dvpudtyfr∂=−−Ω∂ (7.25) 12cos()dwpgudtzfr∂=−−+Ω∂ If we let 2sin()f f =Ω and 2cos()f f =Ω% (they are called Coriolis parameters) then 1 dupfvfwdtx r∂=−+−∂% 1 dvpfudty r∂=−−∂ (7.26) 1 dwpgfudtz r∂=−−+∂% The above equations are EXTREMELY important for meteorology! In the above, we derived the co mponents of the Coriolis force using physical methods. Let's now derive these components using a purely mathematical method . 7-17We choose a Cartesian coordinate system whose x -y coordinate surface is tangent to the surface of the surface (see figure below) . The coordinate origin is at the earth surface of a given latitude and longitude. The Coriolis force per unit mass is given by 2 V −Ω×rr . We can expand this mathematically for the Cartesian -coordinate as follows: ˆ ˆˆ 22 xyzijk V uvw−Ω×=−ΩΩΩrr Of course, we need the components of the earth's rotation vector (,,)xyzΩΩΩ , and we can most easily obtain them by defining a coordinate system with its origin at a point on 7-18the earth's surface with the z -axis parallel to the radial vector ()rr with z increasing away from the earth's surface (local vertical -- see figure above). The x and y axes together form a plane that is tangent to the earth, with x increasing toward the east and y increasing toward t he north. The figure shows the earth's rotation vector at the defined point on the earth, along with its components (projected onto the three Cartesian axes). Note that the x -component is zero (x -axis into the page) because the rotation vector is perpendic ular to the x -axis. Thus, based upon the geometry shown, we find that 0,cos,sinxyzff Ω=Ω=ΩΩ=Ω . Using these in the above determinant, we have ˆ ˆˆ 2(2sin2cos)(2sin)(2cos) Vivwjuku ffff −Ω×=Ω−Ω+−Ω+Ωrr (7.27) Plug (7.27) into (7.17a), we obtain exactly the same equations for vel ocity components as given in (7.25)! The mathematical way is apparently simpler and more elegant. The drawback is that the physical meaning of each individual term is lost in the derivation. Comments: You can see from the equations that the accelerati on due to Coriolis effect is proportional to the velocity – and deflection it causes has to be considered when the force acts on an object for a long time, i.e., when the object travels over a long distance, such as the upper -level air stream does! Theref ore for large -scale atmospheric flows, Coriolis force can not be neglected. Example 1: Suppose a ballistic missile is first due eastward at 43 degrees N latitude. If the missile travels 1000km at a speed of 500 m/s, by how much is it deflected due to the Coriolis force? Given: φ = 43ºN, u of missle = 500m/s distance traveled =1000km = 106 m To find: lateral deflection, The equations to use is (how do we know)? 7-192sin()dvudtf =−Ω . Note. Missile starts with u ? 0 and v =0, but soon is deflected to the right (sourth ) à u ? 0 and v < 0. So we should get v < 0 and Δy < 0, i.e., a southward deflection. Δ y Δy is what we ultimately need. To be Δy, we first have to find v. 51 12sin()2(7.29210)(500/)sin(43) 0.0497dvusmsdt sf−− −=−Ω=−×° =− 0()1 000.0497 ()0.0497vtt vdvsdt dyvttdt− ==− =−=∫∫ () 00 2(0.0497) ()0.04972yTt dytdt TyT=− =−∫∫ What is T? 10000002000500/mTsms== . Therefore 220000.049799.42ykmΔ=−=− . 7-20Example 2: Suppose we drop an object from a tower of height h at the equator. When will it land relative to a "plumb line" at the release point? Let's first try to guess the answer – will it fall to the east or west? In this problem, gravitational and centrifugal forces are in the radial direction, i.e., no projection onto the longitudinal circle (c.f. Figure in page 7 -13). So, our picture is relatively simple: Looking from above: East N.P. Object UP Ω West The Coriolis force is in the equatorial plane (no y deflection). What equation governs the motion? The object falls at a speed w (we use info from c hapter 4 on falling body), and we expect the deflection to be east -west. So, we need an equation that, given w, predicts the evolution of u. Thus, we use equation (7.24) 2sin()2cos()duvwdtff =Ω−Ω . Initially v = 0, so the first term is zero à 2cos(0)2duwwdt=−Ω=−Ω . Now, dwmmgdt=− à 0() 0()wtt wdwgdt =−∫∫ ()wtgt =− 7-21Plug w(t) into du/dt equation and integrating à () 00 22 ()utt dugtdt utgt=Ω =Ω∫∫ The deflection is found from dx/dt = u(t) = Ωgt2 or Δx = Ωgt3/3. Since Δx > 0 , the object falls to the EAST of the drop point. If the height to fall is h, then 0 2/2zt hdzwgtdt dzgtdt zhgt==− =− =−∫∫ For z =0, the time T is 2gTh= à 3/212 3gxghΔ=Ω  For a 500 m tower, Δx ~ 0.77 cm at the equator. This is no different from the merry -go- around example! 7-227.5. Foucault Pendulum We live on the rotating earth. How do we know that the earth is rotating? It's easy to tell if you observe the earth from the space, but not so if you stand at the surface of the earth. There is actually a simple way to demonstrate that the earth is rotating, even inside a closed building. It is called Foucault pendulum. It was invented by French physics, Foucault*, by accident. *Foucaul t, Jean Bernard Leon . 1819 -1868. French physicist who estimated the speed of light and determined that it travels more slowly in water than in air (1850). He also employed a pendulum to prove the rotation of the earth (1851) and invented the gyroscope (185 2). An interesting article on foucault pendulum can be found at http://www.abc.net.au/surf/pendulum/default.htm . See illustration on the next page and a photo of a real thing. A picture of Foucault pendulum at Smithsonian Institution's (http://www.si.edu) National Museum of American History (http://americanhistory.si.edu) in Washington, D.C. Taken in 1996. Courtesy of April Strout. As the pendulum moves, it knocks over the tiny red thin gs around the circumference of the circle. 7-23 7-24It is basically a large pendulum hanging from a tall building / church hall that is fixed to the rotating earth. Let's look at the simple pendulum first. We know for a regular pendulum made up of a mass m hanging from a massless string of length l, θ l m mg sin( q) mg the equation describing the motion is ()sin()ddlmmgdtdtqq=− i.e., sin() lgqq=−&& . Since we assume that the swing is mild, i.e., q remains relatively small, we can approximate sin( q) with q, i.e., sin()qq; (because 2sin()/2.... qqq=++ ) à 0g lqq+=&& . (7.28) This is ex actly the SHO equation we solved earlier in the semester. It has a general solution of ()cos()sin()tAtBtqww =+ 7-25where g lw= (you can show it yourself). It has the unit of s-1 and is the frequency of oscillation. If the pendulum starts from rest at angle θ0 (dθ/dt = 0 at t =0 … from rest), then A = θ0 and B=0 à 0()cos()ttqqw = . (7.29) In the absence of friction, the motion is perpetual and periodic. The period is 22lTgppw== . (7.30) Note that the period is independ ent of the amplitude θ0, and this is why pendulum are good for use in clocks! They are very regular! Now let's go back to the Foucault pendulum. The pendulum is set in motion, if we neglect the effect of friction and observe it from a fixed star , what wi ll we see? We should see that the pendulum behaves like a simple pendulum, doing its periodic motion perpetually, and within a fixed vertical plane whose angle does not change relative to the star , but because the earth is rotating underneath it, the angl e of the plane changes relative to the earth. This is easier to understand if the pendulum is located at the north pole. The earth rotates underneath the pendulum while the hanging point of the pendulum does not change. As a result, we see relative motion between the pendulum and the earth. What if we observe the pendulum standing on the earth ? Will we still see the simple periodic motion of a simple pendulum? Does the earth rotation have any effect? Yes – remember the Coriolis force in a rotating frame? Because the plane in which the pendulum swings is fixed in the absolute coordinate while the earth we stand on rotates counterclockwise relative to that plane, we (who rotates together with the earth) observe that the pendulum plane rotates clockwise rela tive to the us and the earth! In the earth coordinate, such 'mysterious' rotation of the pendulum plane is due to the Coriolis force, which we know always acts perpendicular to the direction of motion. If we view the pendulum from the top, the trajector y is a fixed straight line if there is no Coriolis force: When Coriolis force is present, it deflects the 'bob' of the pendulum to its right, the trajectory becomes something like 7-26 7-27 Coriolis force α Coriolis force Applying the equation of motion in the direction normal to the direction of bob movement, i.e., F cor = ma, we can find (details not show) that the change in the angle α of pendulum plane relative to the ear th is governed by equation sin()d dtaf =−Ω . (7.31) The right hand side is actually the (local) vertical component of the angular velocity of the earth, except for the negative sign. Since the earth rotates counterclockwise, the pendulum p lane rotates clockwise (indicated by the negative sign). So we observe a trajectory depicted in the previous page. The rate of rotation of the pendulum plane is equal to Ω at the north pole ( -Ω at the south pole) and zero at the equator (the Foucault pendulum will not work if placed at the equator). The period of rotation of the pendulum plane is 22 sin()Tpp af==Ω&. (7.32) At the north pole, it is exactly 24 hours, and at 45 degree latitude, it is about 34 hours, and at the equator, there is not rotation and the period is infinitely large. At the equator, we will not see the effect of earth rotation – the Coriolis force is zero there. It is the projection of earth's angular velocity to the local vertical direction [ sin() zf Ω=Ω ] that matters! 7-287.6. Coriolis Force and Weather The Coriolis force is extremely important for large -scale atmospheric flow. It is actually one of the largest force acting on air parcels in the horizontal direction. The balance between the Coriolis force and horizontal pressure gradient force give rise to the so - called geostrophic wind that is the dominant component of winds at the large scales. Earlier, we presente d the equations of motion 1 dupfvfwdtx r∂=−+−∂% 1 dvpfudty r∂=−−∂ (7.33) 1 dwpgfudtz r∂=−−+∂% where 2sin()f f =Ω and 2cos()f f =Ω%. At the large scales (scales on the order of 1000 km), the vertic al wind speed is much smaller than the horizontal speed. Vertical motion becomes large only in small -scale systems such as the thunderstorm, hurricanes and tornadoes. Therefore we focus our attention on horizontal motion (u,v). We also neglect the Coriolis force related to w as appearing in u equation. We therefore have 1 dupfvdtx r∂=−+∂ (7.34a) 1 dvpfudty r∂=−−∂. (7.34b) 7.6.1. Geostrophic Flow The geostrophic winds exist when there is a perfect balance between the pressure gradient force and the Coriolis force. This balance is called geostrophic balance . As a result, the acceleration is zero (assuming friction can be neglected). Scale analysis (you will learn it next semester) shows that this is nearly true for large -scale flows, hence 1 gpfvxr∂−+∂ = 0 (7.35a) 1 gpfuyr∂−−∂ = 0 (7.35b) We denote winds that satisfy geostrophic balance, ˆˆ (,)gggggVuviujv==+r . 7-29 Low pressure p PGFF =r gVr Coriolis cF=r H pressure Equation (7.35) can be written in a vector form: 1ˆ 0gpfkVr−∇−×=r (7.36) (becaus e ˆˆ ˆˆˆˆ() gggggkVkiujvjuiv×=×+=−r ), or 1ˆg pfkVr−∇=×r (7.37) ˆk×(7.37) à 1ˆˆˆ () gg kpfkkVfVr−×∇=××=−rr therefore 1ˆgVkp fr=×∇r (7.38) (7.38) is the vector form of (7.35) (it can also be obta ined directly from 7.35). It defines the geostrophic wind. It says that, in the northern (southern) hemisphere, if you walk in the direction of geostrophic winds, high pressure is always on your right (left). The southern hemisphere is opposite because th e Coriolis parameter f is negative there. P - gVr P+ p∇ 7-30The large -scale (synoptic scale and above) upper -level (where friction can be neglected) atmosph ere flow is nearly geostrophic. Therefore the wind is essentially parallel to the geopotential height (a representation of pressure field on a constant pressure surface) contours or perpendicular to the pressure gradient force. See figure below. Let’s look closer at the geostrophic wind equation 1 gp ufyr∂=−∂. 1) If 0p y∂<∂ (pressure decreases to the north), then u > 0 à high pressure on the right. 2) At the equator, f = 0, the geostrophic wind is undefined – i.e., there is no geostrophic balance to talk about at the equator. 3) f = 2 Wsin(φ) is largest at the poles. Given the same pressure gradient force, the velocity is smallest at the poles, there one needs a smaller u g to obtain the same Coriolis force to ba lance the pressure gradient force. 7-314) If we know that the wind is in a geostrophic balance, we only need to observe the pressure field in order to know the wind field. The geostrophic balance can be used to quality control wind and pressure observations. 5) Near the surface, the geostrophic balance is usually not as well satisfied because of friction, one more force that needs to be considered in the force balance. Also for flows with large curvature, centrifugal force (arising from centripetal accelera tion, not from earth rotation) also needs to be considered. 7.6.2. Curved Flow, Natural Coordinates For straight airflow along straight parallel isobaric contours, when a geostrophic balanced is reached, acceleration due to either speed and direction ch ange is zero. In reality, airflow is rarely perfectly straight, its trajectory is often curved . To be study curved flows, it is convenient to define and use the natural coordinates. Trajectory of particle / air parcel ˆn ˆt Nature coordinates consist of 2 unit vectors, one along the direction of motion ( ˆt) and one normal to and directed to the left of the velocity vector ( ˆn). The coordinate directions are s and n, respectively. Cartesian ˆ,xi ˆ,yj Natural ˆ,st ˆ,nn Here we assume we only deal with 2D motion. In the natural coordinates, because ˆt is always in the direction of velocity, ˆ VVt=r . (7.39) 7-32If s is length along the trajectory, dsVdt=. V is always positive. The acceleration in the natural coordinates is 2ˆˆ()ˆˆ ˆdVdVtdVdtdVVtVtndtdtdtdtdtR==+=+r (7.40) In obtaining (7.40), we used ˆˆdtVndtR= (7.41) where R = radius of curvature following the parcel's motion. It is analogous to the time derivative of the unit ve ctor in radial direction of the plane polar coordinate and we will not derive it rigorously here. Unlike the radius of a circle, the radius of curvature R can have both positive and negative signs. It is defined as center of curvature Vr Vr ˆn ˆn center of curvature The center of curvature on the left side The center of curvature on the right side (+ˆn side) of trajectory R > 0 (-ˆn side) of trajectory, R < 0 Note that R is locally calculated and can change long the trajectory. 7-33In (7.40), we have 2 ˆ ˆdVdVVtndtdtR=+r (7.42) In the equation of motion, 1ˆdVpfkVdt r=−∇−×rr , the pressure gradient force can be written as 111ˆ ˆppptnsn rrr∂∂−∇=−−∂∂. (7.43) What about the Coriolis force ˆfkV−×r ? We know it always points to the right of trajectory (assuming northern hemisphere where f>0. The result we obtain will still be valid when f < 0), therefore ˆ ˆ fkVfVn−×=−r (7.44) Putting all three expressions together à 2 ˆ ˆdVVtndtR+=11ˆ ˆpptnsnrr∂∂−−∂∂ˆfVn− Equating the ˆ ˆ and nt components, we obtain the equations of motion for the respective components: dV dt=1p sr∂−∂ (7.45a) 2V R=1p nr∂−∂fV− (7.45b) Notes: Total acceleration Acceleration due to speed change Acceleration due to direction change – centri petal acceleration 7-34(1) The equations are coupled through V and p. (2) For motion parallel to the isobars, p s∂ ∂= 0 à V = constant following the motion. (3) For purely straig ht flow, R = ∞ à centripetal acceleration and the corresponding centrifugal force is zero. We get the geostrophic relation 1p nr∂−∂fV− = 0 à 1 gp Vfnr∂=−∂. (4) We know tha t if we choose a coordinate system that follows the air parcel undergoing a curved motion, we experience a centrifugal force that is related to the centripetal acceleration see in a fixed coordinate. (7.45b) can be rewritten as 210pVfVnRr∂=−−−∂ (7.46) where 2V R− is the centrifugal force. We will look at two additional types of flow in the following sections. 7.6.3. Cyclostrophic Flow (Read 4 -page handout of today. From Holton – Introduction to Dynamic Meteorology). For flows with large curvature (small R. e.g., tornadoes, dust devils, water sp routs), the centrifugal force >> Coriolis force, therefore the latter can be neglected. We then have the following equation 210pV nRr∂=−−∂ which is a balance between the pressure gradient force and the centrifugal force. The velocity given i s: RpVnr∂=−∂. (7.47) What kind of pressure pattern is possible? Physically, we know that the centrifugal force always directs away from the center of rotation (center of curvature in general), therefore to balance it, the pressure gradient force has to point to the center of curvature. This 7-35means for a closed circulation, the circulation center always has to be a center of low pressure ! Does (7.47) fit our physical reasoning above? There are two possible cases. CASE I: 0and0,is realpRVn∂<>∂. Vr Centuuuuur ˆn L PGFuuuuur ˆn Vr 0and0,is realpRVn∂<>∂ CASE II: 0and0,is realpRVn∂><∂. Vr L Centuuuuur PGFuuuuur ˆn ˆn Vr 0and0,is realpRVn∂><∂ 7-36Therefore , for cyclostrophic circulation, the pressure is always lower at the center and the circulation can be either cyclonic (counterclockwise) and anticyclonic (clockwise). Anticyclonic low-pressure systems rarely occur at the large scale s, however, as the Coriolis comes into play . The balance between pressure gradient force and the Coriolis force becomes dominant and the centrif ugal force is small relative to those two (partly due to large R). 7.6.4. Gradient Wind Gradient flow/wind is a three -way balance between the Coriolis force, pressure grad ient force and the centrifugal force. Because of the balance, the acceleration in the direction of trajectory dV/dt = 0. The equation of force balance is (7.46) 210pVfVnRr∂=−−−∂ solving for V à 1/2 22 24fRfRRpVnr ∂=−±− ∂  . (7.48) The ph ysically meaningful solution ( V is real and positive) exists only in certain situations. The four possible scenarios are illustrated in the following figure s. CASE I : 0and0,pRn∂>>∂then V will be either complex or negative – neither is permitte d. Vr H cFr ˆn PGFuuuuur All three forces point in the same direction – impossible to have a balance. Centuuuuur 7-37 CASE II : 0and0,is realpRVn∂><∂ and positive (with + sign in front of square root) Vr L ˆn cFr PGFuuuuur This is a case of Normal or R egular Low with baric flow . A flow is called baric when Coriolis f orce and PGF are opposite in directions. Such a flow and pressure pattern is typical of the low -pressure systems, mid -latitude cyclones obs erved in the atmosphere. CASE III: 0and0.pRn∂<>∂ Vr PGFuuuuur Centuuuuur ˆn L cFr In this case, the PGF and F c are in the same direction. Such flow is called antibaric . PGF and F c are balanced by the centrifugal force . Obviously the centrifugal force has to be large enough – which requires large speed, large curvature (small radius of curvature) or both. This type of circulation rarely occurs at the synoptic or large scales in the atmosphere where centrifugal force is typically not large enough to balance the other two force s. Centuuuuur 7-38Another reason is that such a circulation is hard to get established. If we have a pressure pattern with low center , an air parcel starting from rest will be driven towards low pressure center and turned to the ri ght by Coriolis force – resulting in a cyclonic rather than anti-cyclonic circulation . Such force balance is possible at the small scale s where F c plays a much smaller role. That circulation is then si milar to the cyclostrophic anti cyclonic circulation we discussed in the p revious section. This is the case of Anomalous low. CASE IV : 0and0.pRn∂<<∂ Vr Centuuuuur ˆn H cFr PGFuuuuur In this case, the radical could be negative, therefore we require that 22 4RpfR nr∂<∂. à get two positive real roots for V. 7-39Positive sign case , 1/222 0 242fRfRRpfRVnr ∂=−+−>−> ∂  , V is larger than the negative sign case. Vr L H CF PGF Cor Since Cor (Coriolis force) is proportional to V and CF (centrifugal forc e) proportional to V squared, V2/R (large R) , the PGF (which is in the same direction as CF) needs to be only relatively s mall for the balance to occur. Even through the above case (with positive sign) is t heoretically possible, it almost never occu r in the atmosphere. At about 43 °N (f = 10-4 s-1), with density ρ =1 km m-3, a radius of curvature R = - 1000 km, and a pre ssure gradient of /1/100pnmbkm∂∂=− , V = 11 m/s with negative sign and about 89 m/s with positive sign. Strong winds are rarely observed in high -pressure systems. Therefore, this is the case of Anomalous High . The negative sign case , 1/222 0 242fRfRRpfRVnr ∂=−−−<−> ∂ , V is smaller than the positive sign case à Cor and CF small er à PGF larger for b alance Vr L H CF PGF Cor 7-40According to earlier discussions, this is the c ase of Normal or Regular High . Subgeostrophic and Supergeostrophic Flows Within sharp troughs in the middle -latitude westerlies the observed velocities are often subgeostrophic by as much as 50%, even though the streamlines still tend to be oriented parallel to t he isobars. These large departures from geostrophic balance are a consequence of the large centripetal acceleration associated with the sharply curved flow within such regions. In order to illustrate how these accelerations alter the balance between the Co riolis force and the pressure gradient force it is convenient to represent them in terms of an apparent centrifugal force. The three -way balance between pressure gradient and Coriolis and centrifugal forces is shown in the following figure for cyclonic and anticyclonically curved trajectories. In both cases, the centrifugal force is directed outward from the center of curvature of the air trajectories (denoted by the dashed lines) and has a magnitude given by 2/ VR where R, is the loca l radius of curvature. In the cyclonic case, the centrifugal force reinforces the Coriolis force so that, in effect, a balance of forces can be achieved with a wind velocity smaller than would be required if the Coriolis force were acting alone. Thus, in t his case, it is possible to maintain a sub - geostrophic flow parallel to the isobars. For the anticyclonically curved trajectory the situation is just the opposite: the centrifugal force opposes the Coriolis force and, in effect, necessitates a super -geostr oghic wind velocity in order to bring about a three - way balance of forces. For the combination of sharp anticyclonic curvature and a strong pressure gradient it is not possible to achieve a three -way balance of forces; that is to say, for all possible val ues of V, 2 nV PfVR>− It is no coincidence, then, that the combination of sharp troughs and tight pressure gradients is not uncommon, whereas sharp ridges are rarely if ever observed in conjunction with tight pressure gradients 7-41 The three -way balance between the horizontal pressure gradient force, the Coriolis force and the centrifugal force, in flow along curved trajectories ( - -) in the Northern Hemisphere. We can also look at the above issues by comparing the gradient wind speed with the geostrophic wind speed correspondin g to the same pressure gradient. Since 1 gp fVnr∂=−∂ where gV is the geostrophic wind speed, equation (7.46) can rewritten as 2 0 gV fVfVR=−− (7.49) and V is the gradient wind speed. Therefore 2 1gV V VfR=+ (7.50) For c yclonic circulation in norther n hemisphere, fR>0, therefore V<Vg, i.e., in flows with cyclonic curvat ure, gradient wind speed (whi ch is a better appr oximation to the rea l wind ) is smaller than the geostr ophic wind speed . One gets sub geostrophic fl ow. Simil arly, V>Vg with antic yclonic curvature.subgeostrophic supergeostrophic 7-427.6.5. Frictional Effect Up to this point, we have neglected the influe nce of friction. Friction actually plays a key role in setting the structure of large -scale atmosphere for the development of deep convection. For now, let's take a more simplistic view of the quantitative effect of friction. In a geostrophic flow, the P GF and the Coriolis force are in balance. When friction is present, friction acts in the opposite direction as the velocity vector and is therefore also perpendicular to the Coriolis force. The balance between the three for ces are illustrated as follows: P - PGFuuuuur P 0 f Vr rFr Coruuuur P+ As can be seen above, the wind vect or is no longer parallel to the pressure contours, the flow is crossing the isobars to the lower pressure side! Based on this balance, it can be shown that this velocity is less than the geostrophic velocity. φ, the angle of departure from geo strophy, can be calculated from the force balance. Friction is usually non -negligible near the ground (in the boundary layer) , where the air is more turbulent and is subject to surface drag. The crossing of air to the low -pressure side at the lower levels contributes significant ly to the low -level conv ergence in low -pressure systems and the convergence induces (through mass continuity – remember the kinematic method for calculating vertical motion from horizontal divergence?) ascending motion that destabilizes the atmosphere so that convection to occur (bad weather) . Sim iliarly, friction induces low-level divergence in high pressure systems à descending motion à good weather. L conve rgence 7-43 Mean -sea level pressure and surface winds, showing the surface winds mostly follow the isobaric contours but the cross -contour (from high to low pressure) component is significant. At a result, low-pressure centers are associated with convergence and high -pressure centers associated with divergence. The convergence forces air in low -pressure system to rise, prod ucing weather.