Next Article in Journal
Characteristics of the Transport of a Typical Pollution Event in the Chengdu Area Based on Remote Sensing Data and Numerical Simulations
Next Article in Special Issue
Satellite Assessments of Tropopause Dry Intrusions Correlated to Mid-Latitude Storms
Previous Article in Journal
Seasonality of Precipitation over Himalayan Watersheds in CORDEX South Asia and their Driving CMIP5 Experiments
Previous Article in Special Issue
Sensitivity of the Reaction Mechanism of the Ozone Depletion Events during the Arctic Spring on the Initial Atmospheric Composition of the Troposphere
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phases of the Isobaric Surface Shapes in the Geostrophic State of the Atmosphere and Connection to the Polar Vortices

Department of General and Theoretical Physics, Institute of Mathematics and Natural Sciences, North Caucasus Federal University, 1 Pushkin Street, 355009 Stavropol, Russia
*
Author to whom correspondence should be addressed.
Atmosphere 2016, 7(10), 126; https://doi.org/10.3390/atmos7100126
Submission received: 3 August 2016 / Accepted: 30 September 2016 / Published: 5 October 2016

Abstract

:
This paper presents a theoretical study of the disturbed isobaric surface shape in the geostrophic state of the atmosphere. It has been shown that, depending on the overheat sign at the equator, the isobaric surface has the shape of an oblate or prolate geoid. If the geostrophic wind velocity is nonzero at the poles, the local pressure extrema (minima for oblate geoid and maxima for prolate geoid) appear at the poles in the geostrophic state. This result correlates with the well-known polar vortex phenomenon and possibly can refine our understanding and interpretation of the phenomenon. In other words, the existence of polar minima and maxima of the pressure field can be the peculiarity of the geostrophic state of the atmosphere. It has been found that air must be colder than the surrounding atmosphere for initiation of the zonal eastward transport. For warm air mass, only easterly winds will be observed.

Graphical Abstract

1. Introduction

As is generally known, the geostrophic state plays a significant part in atmosphere dynamics. It is important for understanding global climate system and local geophysical processes to know the characteristic properties of geostrophic wind [1,2,3,4]. Scale analysis of variables involved in the atmosphere dynamics equation leads to the following expressions for the geostrophic wind velocity components [5,6,7,8]:
u g = 1 2 ω 0 ρ s sin φ p g y
v g = 1 2 ω 0 ρ s sin φ p g x
where ug is the geostrophic wind along the parallel (x-axis) velocity expression; v g the geostrophic wind along the meridian (y-axis) velocity expression; ω0 = 7.29 × 10−5 rad/s the Earth’s angular velocity; ρs the air density of static atmosphere; φ the latitude; and pg the geostrophic pressure disturbance relative to the static state. White and Bromley [9] show that this analysis is only applicable to a shallow atmosphere. Thus, the geostrophic wind velocity can be obtained if the disturbed isobaric surface is known. However, the question of the disturbed isobaric surface shape is still an open question [6,7]. Note that Equations (1) and (2) have been obtained assuming that w u ,   v (where w is the vertical velocity projection) and that the vertical projection of the motion equation comes to the atmosphere static state equation.
The knowledge of the disturbed isobaric surface shape can improve our fundamental understanding of a number of significant atmospheric phenomena such as the large-scale atmospheric circulation, the planetary waves’ dynamics, etc. In particular, in this paper, we place special emphasis on the polar vortex phenomenon. This phenomenon represents a persistent, large-scale low-pressure area located near either of the Earth’s poles [10,11,12]. Distinct polar vortices have also been observed on other planetary bodies of the solar system [13,14,15,16,17,18]. The nature of the polar vortex is still not clear. The question is raised here as to whether the polar vortex results from the atmosphere geostrophic state disturbance owing to the non-uniform heating of the Earth from the equator to the pole or if it is inherent in the geostrophic state itself. At present, a number of various approaches to study the polar vortex formation and development mechanisms exist [19,20,21,22,23,24,25]. Most of them consider the polar extrema of the pressure field as an effect of the geostrophic state disturbance but not as the geostrophic state feature. They associate the polar pressure minima with a geostrophic response to thermal forcing and radiative processes. According to the modern understanding of the polar vortex formation mechanism, it results from the large-scale atmospheric circulation driven by the difference of the Earth surface heating at the poles and at the middle latitudes. Warm air at 60 degrees of north and south latitudes goes up and moves towards the poles in the upper troposphere. Near the poles, the air begins to cool and goes down, making a high-pressure zone. The near-surface air moves between the high-pressure zone and the low-pressure zone of the polar front deflecting to the west under the Coriolis force. As a result, the easterly winds develop at the surface surrounding the pole in a vortex form. At the same time, the polar cyclone forms at the pole.
Contrary to the described scenario, in the present work, we offer a hypothesis that the polar vortex corresponds to the basic state of the atmosphere. We believe that such a persistent, long-lived state should be one of the attributes of the atmosphere geostrophic state, i.e., this is a steady state. Any deviation from the geostrophic state will lead to the ageostrophic oscillations. Particularly, the circulatory mechanism described above will disturb the state of the polar vortex, but it is not the reason for the polar vortex appearance. It can be expected that the analysis of the isobaric surface geometry should demonstrate correlation with the fact of existence of a pronounced low-pressure area near the poles. The goal of the present work is to determine (even if qualitatively) the disturbed isobaric surface shape in the geostrophic state of the atmosphere. This will allow solving the above-mentioned problems.

2. Main Equations

The expressions for the geostrophic wind components can be obtained from the analysis of the atmospheric dynamics equation in the vector form [5,6,7,8,26]:
v t + ( v ) v = g 0 1 ρ i p + 2 [ v ω 0 ] + ω 0 2 R + f fr
where g0 is the gravity acceleration; p the pressure gradient; 2[0] the Coriolis acceleration; ω 0 2 R the centrifugal acceleration; R the distance between the air parcel location and the rotation axis; ffr the frictional force per unit mass; and ρi the density of moving air (in the static state, the density will be denoted as ρs). A system of coordinates is related to the Earth’s surface; abscissa axis (x-axis) is directed along the parallel; ordinate axis (y-axis) is directed along meridian; and applicate axis (z-axis) is perpendicular to the Earth’s surface (Figure 1). Equation (3) is Newton’s second law for the rotating Earth’s atmosphere.
In the atmosphere static condition, when v = 0, Equation (3) will have the form:
0 = g 0 1 ρ s p s + ω 0 2 R
The static condition corresponds to the balance of the gravity force, the pressure gradient and the centrifugal force. It is more convenient to introduce the free fall acceleration vector equal to the gravity acceleration g0 and centrifugal acceleration vectorial sum:
g = g 0 + ω 0 2 R
Hence, the Earth’s geoidal surface is perpendicular to the free fall acceleration g.
Then, the atmosphere static state equation has the form:
0 = g 1 ρ s p s
Under this form, the static state corresponds to the balance of the gravity force and the pressure gradient. In such notation, it follows from here that, in the static state, the isobaric surfaces are perpendicular to the free fall acceleration vector, i.e., parallel to the Earth’s geoidal surface. These isobaric surfaces are assumed as the atmosphere undisturbed state. Note that, in the static state, the undisturbed isobaric surface is usually approximated by a sphere [5,6,27,28,29]. Actually, the undisturbed isobaric surfaces are geoid-like. Thus, the static state is the basic state of the atmosphere in the so-called zero approximation. In this case, the baric and velocity fields are exactly determined.
At the steady motion dv/dt = 0, the isobaric surfaces of geoidal shape become disturbed; hence, the pressure may be written as [7]:
p = p s + p g
Similarly, the air density in Boussinesq approximation may be presented as [5,6]:
ρ g = ρ s ( 1 α Δ T )
where Δ T = T g T s is the overheat function which is the air temperature disturbance relative to the static state: α = 1/T0, T0 = 273 K. The overheat function has a positive sign at warm air mass motion and a negative sign at cold air mass motion. The geostrophic state is the next basic state of the atmosphere in the so-called first approximation. In this case, the baric and velocity fields also should be determined exactly.
The equation of frictionless steady motion ffr = 0 will have the following form:
0 = g 1 ρ g ( p s + p g ) + 2 [ v , ω 0 ] α Δ T ρ s p s 1 ρ s p g + 2 [ v , ω 0 ] = α Δ T g 1 ρ s p g + 2 [ v , ω 0 ]
The geostrophic state of the atmosphere is determined by the balance of the gravity force, the pressure gradient and the Coriolis force. Here, we determined that, according to Equation (8), ρ g ρ s .
The following expressions define the projections of the Earth rotation angular velocity [5,6,30]:
ω 0 x = 0 ,   ω 0 y = ω 0 cos φ ,   ω 0 z = ω 0 sin φ
Consider the projections of the stationary state motion Equation (9) in the coordinates system where the x-y-plane is tangent to geoid:
0 = 1 ρ s p g x + 2 v ω 0 sin φ 2 w ω 0 cos φ
0 = 1 ρ s p g y 2 u ω 0 sin φ
0 = 1 ρ s p g z + α Δ T g   + 2 u ω 0 cos φ
The free fall acceleration in the last equation equals g = g 0 R E ω 0 2 cos 2 φ , where RE is the Earth’s radius. The estimation of quantities shows that R E ω 0 2 cos 2 φ g 0 .
The third term in Equation (11) presents a product of the air velocity vertical component and the Earth’s angular velocity horizontal projection. Thus, one can conclude that there are two ways to obtain Equations (1) and (2): either by neglecting the vertical component of the air velocity, or by neglecting the horizontal projection of the Earth’s angular velocity. Although both assumptions lead to identical equations for the geostrophic wind projections, these assumptions are not equivalent. Within the geostrophic model of the atmosphere, on the basis of the scale analysis of variables involved in Equation (3) in the thin-layer approximation, it has usually been concluded that only the angular velocity component normal to the Earth’s surface is essential for the atmospheric dynamics. Thus, the third term in Equations (11) and (13) is usually neglected [5,6]. On the other hand, the importance of the Earth’s angular velocity horizontal projection for the atmospheric dynamics has been demonstrated in [9]. Thus, we can say that the geostrophic state of the atmosphere is a steady state at which the friction forces and the air velocity vertical component can be neglected. In such a situation, Equations (1) and (2) describing the geostrophic state of the atmosphere can be obtained. However, here we have the additional Equation (13) for the description of the geostrophic state. As a matter of fact, the goal of the present study is to reveal the novel aspects of the geostrophic state description resulting due to the accounting of Equation (13).
Thus, contrary to the usual definition of the geostrophic state, in our analysis, we use Equation (13). As a rule, the atmosphere static state equation appears as a third equation in the system describing the geostrophic state (see, e.g., [5,6]). This is due to the fact that the vertical projection of pressure gradient in Equation (13) is many orders of magnitude smaller than the Coriolis acceleration. However, from our point of view, it is incorrect to ignore Equation (13).
Due to the principal importance of this point, here we reconsider the scale analysis of the geostrophic state. Let us introduce the horizontal L and vertical H length scales such as δ = H / L 1 . We have the following estimation of the Equation (11) components:
0 ~ p g x ρ s L + 2 V ω 0 2 W ω 0
The continuity equation leads to
V L ~ W H , W ~ V H L = δ V
i.e., W V . From here, we have
0 ~ p g x ρ s L + 2 V ω 0 ( 1 δ )
As it was noted in [6], to have the same order of magnitude of the horizontal pressure gradient and of the Coriolis acceleration, we must suppose that
p g x ~ 2 ρ s L ω 0 V
Thus, the third component of Equation (11) can be neglected as the component of lower order. Equation (12) gives the same order of magnitude for the pressure disturbance. The estimation of the vertical projection of the pressure gradient was made in [6] in the following way. Assuming that the disturbance of pgz has the same order as the disturbance of pgx (i.e., 2ρs0V), it can be concluded that the ratio of the Coriolis acceleration and the vertical pressure gradient equals δ. Thus, the components of Equation (13) are of lower order than the components of Equations (11) and (12). For that reason, Equation (13) is ignored and the static state equation is used. However, how valid are the estimations presented above?
Suppose that the pressure disturbance has the form:
p g = p m x sin ( k x x ) + p m y sin ( k y y ) + p m z sin ( k z z )
where kxL = π, kyL = π, kzL = π. From Equations (11) and (12), it follows that
p m x ~ p m y ~ 2 ρ s L V ω 0 π
and, from Equation (13), we have
p m z ~ 2 ρ s H V ω 0 π
From here, it follows that p m z p m x ,   p m y . However, if we suppose that the pressure disturbance has the form
p g = p m e k z sin ( k x x + k y y )
where k~1/L, then the vertical pressure gradient will have the same order of magnitude as the horizontal components. In other words, having no expression for the pressure disturbance, we cannot estimate the components of the pressure disturbance gradient.
Keeping in mind the discussion above, the systems (11)–(13) yield the geostrophic wind velocity horizontal projections:
u g = 1 2 ω 0 ρ s sin φ p g y
v g = 1 2 ω 0 ρ s sin φ p g x
u g = 1 2 ω 0 ρ s cos φ p g z α g 2 ω 0 cos φ Δ T
Taking the vector product of Equation (9) and the vector k, we get the expression for the geostrophic wind velocity:
v g = 1 2 ω 0 ρ s ( k , k 0 ) [ k , p g ]
where k is the unit vector directed vertically up in the direction of z-axis perpendicular to the Earth’s geoidal surface; and k0 the unit vector directed in the direction of the Earth rotation angular velocity. As is seen from Equation (25), the geostrophic wind is perpendicular to the pressure gradient and is directed along the isobaric surface.
As it was noted in [6], the isobaric surface disturbance cannot be determined within the geostrophic model, and a more general theory is required. However, as it will be demonstrated, some information about the shape of the disturbed isobaric surface can still be extracted from the analysis of the geostrophic wind.
Consider a special case when p g / x = 0 and p g / y > 0 . The pressure disturbance at steady motion will decrease along the y-axis from the equator to the pole (this is observed in the atmosphere on a global scale), and the geostrophic wind direction will be west to east, i.e., western flow will be observed in this case. It follows from Equations (22) and (24) that, for the initiation of warm air zonal eastward transport, the following condition must be satisfied:
p g / z > α Δ T ρ s g
and the pressure disturbance gradients must obey the relation:
Δ T = 1 ρ s α g ( p g z + ctg φ p g y )
It follows that for the existence of the zonal eastward transport of warm air, the vertical pressure disturbance gradient must be positive and greater than a certain value:
p g z > ctg φ p g y ,   p s y > 0
If this condition is not satisfied, then for the warm air mass, only east wind can be observed and the cold air mass will move in the eastward direction. The same result will take place if p g / z 0 .
In other words, if we assume that the geostrophic state of the atmosphere corresponds to the isobaric surface of certain shape, which still needs to be determined, then at a given isobaric surface shape, the warm air will move westward (on condition that p g / z 0 ) and the cold air will move eastward (on condition that Δ T < 1 α g ρ s p g / z ).

3. Shape of the Disturbed Isobaric Surface

In the general case, it follows from Equations (22)–(24) that the pressure disturbance gradient p s = ( p s / x ,   p s / y ,   p s / z ) has three nonzero components. Thus, if in the static state ( p s = ρ s g ) the pressure gradient is directed along g, then in the general case of the geostrophic state at the arbitrary point of isobaric surface, the total pressure gradient
p = p s + p g
will be deflected from the g direction. In this case, contrary to the static state, the isobaric surface will not be perpendicular to g because it is perpendicular to p . Depending on the signs of components of the pressure disturbance gradient, the total pressure gradient can be deflected from g both towards the rotation axis and away from the rotation axis. In the first case, the isobaric surface will be extended along the Earth’s rotation axis in comparison with the isobaric surface in the static state (Figure 2a). Such a surface we will call the prolate geoid for short. In such a case, the pressure at the poles will be greater than the pressure in the static state. In the second case, the isobaric surface will be flattened at the poles along the Earth’s rotation axis in comparison with the isobaric surface in the static state (Figure 2b). Such a surface we will call the oblate geoid for short. In such a case, the pressure at the poles will be lower than the pressure in the static state.
The above reasoning is evidently inapplicable to the points at the equator and at the poles. Consider the points at the equator. It should be noted that, as sinφ = 0 at the equator, it has usually been concluded from Equations (22) and (23) that the geostrophic approximation does not work at the equator [6]. Actually, such a conclusion is unfounded if we do not know the disturbed isobaric surface. Supposing that the pressure gradient components are zero at the equator, we have the indefinite form (0/0). If we assume, for example, that pg~cos2φ, then p g / y ~ 2   c o s φ   s i n φ , and the indefinite form vanishes.
To address a question of the disturbed isobaric surface shape in the geostrophic state at the equator, consider the rectangular coordinate system shown in Figure 1 with the point of origin located at the equator. The vector Equation (9) is valid at any point. Neglecting the velocity vertical component, write the projections of Equation (9) onto the coordinate axes:
p g / x = 0 ,   p g / y = 0
u g = 1 2 ω 0 y ( 1 ρ s p g z g α Δ T )
Here, at the equator, ω0y = ω0 and the other components of the angular velocity are equal to zero. One can note that Equation (31) coincides with Equation (24); this fact additionally indicates the importance of Equation (24) for the analysis of the geostrophic state of the atmosphere. From Equations (30) and (31), it follows that p is co-directional with p s . Thus, the isobaric surface is perpendicular to g at the equator, i.e., the isobaric surface is parallel to the geoidal isobaric surface of the static state. Equations (30) and (31) should be used at the equator instead of Equations (22)–(24). It follows from Equations (30) and (31) that the geostrophic wind velocity has only one component at the equator. This velocity component is directed along the equator. It follows from Equation (31) that this velocity component is positive if the following condition is satisfied:
p g z > ρ s g α Δ T
If we assume that p g / z = 0 , then Equation (31) dictates that the warm air ( Δ T > 0 ) will move to the negative direction (eastward), and cold air will move to the positive direction.
At any point infinitely near the equator, the meridional component of the velocity will also take place. From the continuity of motion, it can be concluded that the warm air starting from the equator goes along the spiral trajectory in the negative direction towards the pole. Analogously, the cold air will move from the equator in the positive direction along the spiral trajectory towards the pole. It should be noted that within the geostrophic model, the curved trajectory cannot be realized. However, we can consider the local direction of the geostrophic wind at each point.
Let us next analyze the disturbed isobaric surface shape at the poles. To this end, consider the rectangular coordinate system shown in Figure 1 with point of origin located at the pole (R = 0). If Δ T 0 , the atmosphere motion equation at the pole will be
α Δ T g 0 1 ρ e p s | 0 + 2 [ v g ω 0 ] = 0
Taking into account that 2[vgω0]z = 0, the projections of this equation onto the coordinate axes have the form:
1 ρ s p g x | 0 + 2 v g ω 0 = 0
1 ρ s p g y | 0 2 u g ω 0 = 0
α Δ T g 0 1 ρ s p g z = 0
It can be seen that all three components of the pressure disturbance gradient p g are nonzero, and the geostrophic wind velocity is also nonzero. Thus, in the geostrophic state, the total pressure gradient p = p s + p g will be deflected from the g0 direction. In this case, contrary to the static state ( p s = ρ s g 0 ), the isobaric surface will not be perpendicular to the rotation axis (or to g0).
Consider the isobaric surface shape at the poles when overheat is equal to zero ( Δ T = 0 ). As is seen from Equations (34) and (35), two cases are possible. In the first case, the geostrophic wind velocity is equal to zero at the pole. Then, the x- and y-components of the pressure disturbance gradient are equal to zero. In this case, we have the isobaric surface shape in the form of prolate or oblate geoids (Figure 2). In the second case, the geostrophic wind velocity is nonzero at the pole. Then, the x- and y-components of the pressure disturbance gradient are also nonzero. In this case, the pressure disturbance gradient vector p g makes a right angle (π/2) with the static state pressure gradient vector p s . Hence, the resulting pressure gradient vector p will be deflected from vertical (i.e., from p s ). This will lead to the disturbance of the isobaric surface at the pole. The isobaric surface will not be perpendicular to the rotation axis at the pole. The distinct conclusion about the type of the isobaric surface disturbance at the pole cannot be made on the basis of Equations (34)–(36). It is only clear that the pole is the exceptional point (in the sense of smoothness) of the function describing the isobaric surface. Then, the following scenarios are possible. The first scenario is the isobaric surface of an oblate geoid shape with local maxima or minima at the poles. The second scenario is the isobaric surface of a prolate geoid shape with local maxima or minima at the poles. However, if we assume that the baric minimum corresponds to the cyclonic motion, and baric maximum corresponds to the anticyclonic motion, then, as it follows from the continuity of motion, we must exclude the cases of oblate geoid with baric maxima at the poles and prolate geoid with baric minima at the poles. Thus, when the overheat at the equator is positive, the isobaric surface will have the prolate geoidal shape, and the warm air will move along the spiral towards the pole where the local maximum is taking place and the wind velocity has a finite value (Figure 3a). In this case, the polar vortex is of anticyclonic character. When the overheat at the equator is negative, the isobaric surface will have the oblate geoidal shape and the cold air will move along the spiral towards the pole where the local minimum is taking place and the wind velocity has a finite value (Figure 3b). In this case, the polar vortex is of cyclonic character.
Both described scenarios can exist simultaneously. The observations show that near the poles, the atmospheric circulation has the form of a westerly swirling near-surface vortex and an easterly swirling high-altitude vortex [31].
The observations show that the air temperature inside the polar vortex can be both higher and lower than the temperature of the surrounding atmosphere [10,13,16,17]. Consider the influence of the overheat sign on the isobaric surface disturbance at the pole. We start the analysis from the case of zero geostrophic wind velocity. At the positive values of the overheat ( Δ T > 0 ), it follows from Equation (36) that the vertical component of the pressure disturbance gradient is positive p g / z > 0 and directed oppositely to p s . Thus, the resulting pressure gradient decreases. Consider now the case of nonzero geostrophic wind velocity. At the negative values of the overheat, it follows from Equations (34)–(36) that the pressure disturbance gradient p g makes an acute angle with the static state pressure gradient p s (or with g0) as the vertical component p g / z is down-directed. Therefore, the resulting vector p will be directed along the diagonal of the parallelogram made up of these two vectors. The deflection of the resulting vector p from vertical (from p s ) will be smaller than the deflection at zero overheat. It can be concluded that the isobaric surface has more flat minima at the poles. In other words, the negative overheat diminishes the pressure minimum depth at the pole. When the overheat value is positive, the angle between p g and p s will be greater than π/2. Therefore, in this case, the resulting vector p will have a greater deflection from g0 than that in the previous cases. In other words, the positive overheat elevates the pressure maximum at the pole.
These conclusions correlate with the research results described in [32]. According to [32], the polar vortex weakening is accompanied by the cold winters in mid-latitudes. In [32], it is asserted that the polar vortex weakening leads to the fall in temperature, namely, it induces a negative phase of Arctic Oscillation resulting in low temperatures in mid-latitudes. From our point of view, the opposite cause-and-effect relation is taking place. As it was shown, the steady geostrophic state with polar pressure extrema can exist. Any deviation from the geostrophic state leads to the ageostrophic oscillations. Particularly, negative phase of Arctic Oscillation leads to the polar vortex weakening.
The presented analysis does not give any preference to the appearance of either polar pressure maximum or polar pressure minimum; it only demonstrates the possibility of extrema of the pressure field at the poles in the geostrophic state. However, it should be noted that, at the poles, the pressure minima are usually observed and the air motion is of a cyclonic type [10,13,16,17].
Let us find the geostrophic wind divergence:
u g x + v g y = 1 2 ω 0 { y ( 1 ρ g sin φ p g x ) x ( 1 ρ g sin φ p g y ) } = 1 2 ω 0 ρ s sin φ { p g x ( 1 ρ s ρ g y + ctg φ φ y ) p g y ( 1 ρ s ρ g x + ctg φ φ x ) }
As φ / x = 0 , we have
u g x + v g y = 1 2 ω 0 ρ s sin φ [ p g x ( 1 ρ s ρ g y + ctg φ φ y ) 1 ρ s p g y ρ g x ]
From here, it follows that the divergence is nonzero in the geostrophic regime of the atmosphere:
u g x + v g y = 1 2 ω 0 ρ s sin φ [ p g x ( α Δ T y + ctg φ φ y ) + α Δ T x p g y ] 0
Setting the divergence equal to zero, we find the condition at which this takes place:
p g / x p g / y = ( d y d x ) g = ln ρ g x ln ( ρ g sin φ ) y = α Δ T x ln ( ρ g sin φ ) y = 0
It follows that the slope ratio of the tangent to the isobar and the parallel is determined by the horizontal gradients of air density (or by the gradients of the temperature disturbance) along the parallel and meridian. If the temperature does not change along the parallel, then (dy/dx)g = 0. It follows that the geostrophic wind is directed along the parallel ( v g = 0 ) in this case. In the general case Δ T / x 0 , the geostrophic wind divergence is nonzero and both of the velocity components are also nonzero. It follows that the motion along the spiral from the equator towards the pole will take place.
The velocity divergence at the equator is
u g x + v g y = 1 2 ω 0 y x ( 1 ρ g p g z g α Δ T ) = 1 2 ω 0 y [ 1 ρ s 2 p g x z + α Δ T x ( 1 ρ s p g z g ) ] = 1 2 ω 0 y α Δ T x ( 1 ρ s p g z g )
It can be seen that the geostrophic wind divergence equals zero at the equator only if Δ T / x = 0 . In general, the geostrophic wind divergence is nonzero at the equator, i.e., the airflow lines have a source at the equator. The velocity divergence at the pole is
u g x + v g y = 1 2 ω 0 [ x ( 1 ρ g p g y ) y ( 1 ρ g p g x ) ] = 1 2 ω 0 [ p g y x ( 1 ρ g ) p g x y ( 1 ρ g ) ] = α 2 ρ s ω 0 ( p g y Δ T x p g x Δ T y ) = 1 2 ω 0 [ p g y x ( 1 ρ g ) p g x y ( 1 ρ g ) ] = α 2 ρ s ω 0 ( p g y Δ T x p g x Δ T y ) = α 2 ρ s ω 0 [ grad   p g , grad   Δ T ] z
It can be seen that the geostrophic wind divergence equals zero at the pole if either the velocity is zero or the pressure and temperature disturbance gradients are parallel to each other. In general, the geostrophic wind divergence is nonzero at the pole, i.e., the airflow lines have a sink at the pole.
Let us find the geostrophic wind vorticity at the equator:
Ω g z = v g x u g y = 1 2 ω 0 y ( 1 ρ g p g z g α Δ T ) = 1 2 ω 0 [ 1 ρ s 2 p g y z + α Δ T y ( 1 ρ s p g z g ) ] = α 2 ω 0 Δ T y ( 1 ρ s p g z g )
It is seen that the geostrophic wind vorticity is nonzero at the equator in the general case. The geostrophic wind vorticity at the pole is
Ω g z = v g x u g y = 1 2 ω 0 [ x ( 1 ρ g p g x ) + y ( 1 ρ g p g y ) ] = 1 2 ω 0 [ 1 ρ s ( 2 p g x 2 + 2 p g y 2 ) + α ρ s p g x Δ T x + α ρ s p g y Δ T y ] = 1 2 ρ s ω 0 [ 2 p g + α ( p g , Δ T ) ]
where is the horizontal nabla operator. It is seen that the wind vorticity equals zero at the pole only if the velocity is zero. These results confirm the described picture and show that the air motion from the equator pointing towards the pole will take place along the spiral trajectory.
Thus, the presented discussion of the disturbed isobaric surface geometry demonstrates the existence of a pronounced persistent low- or high-pressure area near the poles. As follows from the analysis above, the polar vortices may be an inherent attribute of the geostrophic state of the atmosphere. In other words, the polar vortex corresponds to the basic state of the atmosphere, which can be disturbed by various factors such as the Arctic oscillations, the quasi biannual oscillations, the southern oscillations, El Niño, etc.

4. Conclusions

Thus, in the present work, it has been found that for the wind velocity projection onto the parallel, two equivalent expressions exist: the first demonstrates the velocity dependence on the pressure gradient along the parallel; and the second demonstrates the velocity dependence on the vertical pressure gradient. This result is obtained owing to the use of Equation (13), which has been traditionally neglected. When p s / y < 0 , for the existence of the zonal eastward transport of warm air, the vertical pressure disturbance gradient must be positive and greater than a certain value. Otherwise, for the warm air mass, only an east wind can be detected and the cold air mass will move in the eastward direction.
It has been demonstrated that the following alternative isobaric surface geometries are possible in the geostrophic state. The isobaric surface has a shape of an oblate or prolate geoid and the pressure at the pole is correspondingly lower or higher than the pressure in the static state. The geostrophic wind velocity, divergence and vorticity are equal to zero in these two cases. The prolate geoidal shape of the isobaric surface corresponds to the positive value of the overheating at the equator; and the oblate geoidal shape corresponds to the negative value of the overheating at the equator. The pressure minimum and maximum can occur at the poles in the mentioned cases of the oblate and prolate geoid. In such instances, the geostrophic wind velocity is nonzero at the poles. It follows from here that the polar vortexes can be a special feature of the geostrophic state of the atmosphere. However, the problem of mathematical definition of the exact shape of disturbed isobaric surface in the geostrophic state remains open. There is no doubt that, in addition to the geostrophic nature, various geophysical factors can influence the polar vortex formation process and complicate the resulting picture of the phenomenon.

Acknowledgments

This work was partially supported by the Ministry of Education and Science of the Russian Federation in the framework of the base part of the governmental ordering for scientific research works (order No. 2014/216, project 653).

Author Contributions

Robert Zakinyan planned and supervised the research, co-performed the theoretical analysis, and co-wrote the paper; Arthur Zakinyan co-performed the theoretical analysis and co-wrote the paper; and Roman Ryzhkov co-performed the theoretical analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allen, R.J.; Sherwood, S.C. Warming maximum in the tropical upper troposphere deduced from thermal winds. Nat. Geosci. 2008, 1, 399–403. [Google Scholar] [CrossRef]
  2. Drobinski, P.; Rotunno, R.; Dubos, T. Linear theory of the sea breeze in a thermal wind. Q. J. R. Meteorol. Soc. 2011, 137, 1602–1609. [Google Scholar] [CrossRef]
  3. Gregow, H.; Ruosteenoja, K.; Pimenoff, N.; Jylhä, K. Changes in the mean and extreme geostrophic wind speeds in Northern Europe until 2100 based on nine global climate models. Int. J. Climatol. 2011, 32, 1834–1846. [Google Scholar] [CrossRef]
  4. Kimura, N.; Wakatsuchi, M. Relationship between sea-ice motion and geostrophic wind in the Northern Hemisphere. Geophys. Res. Lett. 2000, 27, 3735–3738. [Google Scholar] [CrossRef]
  5. Gill, A.E. Atmosphere-Ocean Dynamics; Academic Press: San Diego, CA, USA, 1982. [Google Scholar]
  6. Pedlosky, J. Geophysical Fluid Dynamics; Springer: New York, NY, USA, 1987. [Google Scholar]
  7. Zdunkowski, W.; Bott, A. Dynamics of the Atmosphere: A Course in Theoretical Meteorology; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  8. Monin, A.S. Theoretical Geophysical Fluid Dynamics; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1990. [Google Scholar]
  9. White, A.A.; Bromley, R.A. Dynamically consistent, quasi-hydrostatic equations for global models with a complete representation of the Coriolis force. Q. J. R. Meteorol. Soc. 1995, 121, 399–418. [Google Scholar] [CrossRef]
  10. Schoeberl, M.R.; Hartmann, D.L. The dynamics of the stratospheric polar vortex and its relation to springtime ozone depletions. Science 1991, 251, 46–52. [Google Scholar] [CrossRef] [PubMed]
  11. Harvey, V.L.; Pierce, R.B.; Fairlie, T.D.; Hitchman, M.H. A climatology of stratospheric polar vortices and anticyclones. J. Geophys. Res. 2002, 107, 4442. [Google Scholar] [CrossRef]
  12. Kolstad, E.W.; Breiteig, T.; Scaife, A.A. The association between stratospheric weak polar vortex events and cold air outbreaks in the Northern Hemisphere. Q. J. R. Meteorol. Soc. 2010, 136, 886–893. [Google Scholar] [CrossRef]
  13. Mitchell, D.M.; Montabone, L.; Thomson, S.; Read, P.L. Polar vortices on Earth and Mars: A comparative study of the climatology and variability from reanalyses. Q. J. R. Meteorol. Soc. 2015, 141, 550–562. [Google Scholar] [CrossRef] [PubMed]
  14. Orton, G.S.; Yanamandra-Fisher, P.A. Saturn’s temperature field from high-resolution middle-infrared imaging. Science 2005, 307, 696–698. [Google Scholar] [CrossRef] [PubMed]
  15. Luszcz-Cook, S.H.; de Pater, I.; Ádámkovics, M.; Hammel, H.B. Seeing double at Neptune’s south pole. Icarus 2010, 208, 938–944. [Google Scholar] [CrossRef]
  16. Dyudina, U.A.; Ingersoll, A.P.; Ewald, S.P.; Vasavada, A.R.; West, R.A.; Baines, K.H.; Momary, T.W.; Del Genio, A.D.; Barbara, J.M.; Porco, C.C.; et al. Saturn’s south polar vortex compared to other large vortices in the Solar System. Icarus 2009, 202, 240–248. [Google Scholar] [CrossRef]
  17. Teanby, N.A.; de Kok, R.; Irwin, P.G.J.; Osprey, S.; Vinatier, S.; Gierasch, P.J.; Read, P.L.; Flasar, F.M.; Conrath, B.J.; Achterberg, R.K.; et al. Titan’s winter polar vortex structure revealed by chemical tracers. J. Geophys. Res. 2008, 113, E12003. [Google Scholar] [CrossRef]
  18. Luz, D.; Berry, D.L.; Piccioni, G.; Drossart, P.; Politi, R.; Wilson, C.F.; Erard, S.; Nuccilli, F. Venus’s southern polar vortex reveals precessing circulation. Science 2011, 332, 577–580. [Google Scholar] [CrossRef] [PubMed]
  19. Cavallo, S.M.; Hakim, G.J. Potential vorticity diagnosis of a tropopause polar cyclone. Mon. Weather Rev. 2009, 137, 1358–1371. [Google Scholar] [CrossRef]
  20. Cavallo, S.M.; Hakim, G.J. Physical mechanisms of tropopause polar vortex intensity change. J. Atmos. Sci. 2013, 70, 3359–3373. [Google Scholar] [CrossRef]
  21. Wirth, V. Diabatic heating in an axisymmetric cut-off cyclone and related stratosphere-troposphere exchange. Q. J. R. Meteorol. Soc. 1995, 121, 127–147. [Google Scholar] [CrossRef]
  22. Limpasuvan, V.; Hartmann, D.L.; Thompson, D.W.J.; Jeev, K.; Yung, Y.L. Stratosphere-troposphere evolution during polar vortex intensification. J. Geophys. Res. 2005, 110, D24101. [Google Scholar] [CrossRef]
  23. Watson, P.A.G.; Gray, L.J. How does the quasi-biennial oscillation affect the stratospheric polar vortex? J. Atmos. Sci. 2014, 71, 391–409. [Google Scholar] [CrossRef]
  24. Methven, J.; Berrisford, P. The slowly evolving background state of the atmosphere. Q. J. R. Meteorol. Soc. 2015, 141, 2237–2258. [Google Scholar] [CrossRef]
  25. Nakamura, N.; Solomon, A. Finite-amplitude wave activity and mean flow adjustments in the atmospheric general circulation. Part II: Analysis in the isentropic coordinate. J. Atmos. Sci. 2011, 68, 2783–2799. [Google Scholar] [CrossRef]
  26. Drazin, P.G. Introduction to Hydrodynamic Stability; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  27. Van der Toorn, R.; Zimmerman, J.T.F. On the spherical approximation of the geopotential in geophysical fluid dynamics and the use of a spherical coordinate system. Geophys. Astrophys. Fluid Dyn. 2008, 102, 349–371. [Google Scholar] [CrossRef]
  28. White, A.A.; Staniforth, A.; Wood, N. Spheroidal coordinate systems for modeling global atmosphere. Q. J. R. Meteorol. Soc. 2008, 134, 261–270. [Google Scholar] [CrossRef]
  29. De la Torre Juárez, M.; Fisher, B.M.; Orton, G.S. Large scale geostrophic winds with a full representation of the Coriolis force: Application to IR observations of the upper Jovian troposphere. Geophys. Astrophys. Fluid Dyn. 2002, 96, 87–114. [Google Scholar] [CrossRef]
  30. Sidorenkov, N.S. The Interaction between Earth’s Rotation and Geophysical Processes; Wiley-VCH: Weinheim, Germany, 2009. [Google Scholar]
  31. Frauenfeld, O.W.; Davis, R.E. Northern Hemisphere circumpolar vortex trends and climate change implications. J. Geophys. Res. 2003, 108, 4423. [Google Scholar] [CrossRef]
  32. Kim, B.-M.; Son, S.-W.; Min, S.-K.; Jeong, J.-H.; Kim, S.-J.; Zhang, X.; Shim, T.; Yoon, J.-H. Weakening of the stratospheric polar vortex by Arctic sea-ice loss. Nat. Commun. 2014, 5, 4646. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Illustration of the coordinate system.
Figure 1. Illustration of the coordinate system.
Atmosphere 07 00126 g001
Figure 2. Isobaric surfaces in the atmosphere geostrophic model: 1—undisturbed; 2—disturbed. (a) the prolate geoid; and (b) the oblate geoid.
Figure 2. Isobaric surfaces in the atmosphere geostrophic model: 1—undisturbed; 2—disturbed. (a) the prolate geoid; and (b) the oblate geoid.
Atmosphere 07 00126 g002
Figure 3. Isobaric surfaces in the atmosphere geostrophic model: 1—undisturbed; 2—disturbed. (a) the local maxima at the poles when the overheat at the equator is positive; and (b) the local minima at the poles when the overheat at the equator is negative.
Figure 3. Isobaric surfaces in the atmosphere geostrophic model: 1—undisturbed; 2—disturbed. (a) the local maxima at the poles when the overheat at the equator is positive; and (b) the local minima at the poles when the overheat at the equator is negative.
Atmosphere 07 00126 g003

Share and Cite

MDPI and ACS Style

Zakinyan, R.; Zakinyan, A.; Ryzhkov, R. Phases of the Isobaric Surface Shapes in the Geostrophic State of the Atmosphere and Connection to the Polar Vortices. Atmosphere 2016, 7, 126. https://doi.org/10.3390/atmos7100126

AMA Style

Zakinyan R, Zakinyan A, Ryzhkov R. Phases of the Isobaric Surface Shapes in the Geostrophic State of the Atmosphere and Connection to the Polar Vortices. Atmosphere. 2016; 7(10):126. https://doi.org/10.3390/atmos7100126

Chicago/Turabian Style

Zakinyan, Robert, Arthur Zakinyan, and Roman Ryzhkov. 2016. "Phases of the Isobaric Surface Shapes in the Geostrophic State of the Atmosphere and Connection to the Polar Vortices" Atmosphere 7, no. 10: 126. https://doi.org/10.3390/atmos7100126

APA Style

Zakinyan, R., Zakinyan, A., & Ryzhkov, R. (2016). Phases of the Isobaric Surface Shapes in the Geostrophic State of the Atmosphere and Connection to the Polar Vortices. Atmosphere, 7(10), 126. https://doi.org/10.3390/atmos7100126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop