Research Article Volume 11 Issue 3
Designing the similarity solution associated with nonlinear zonal atmospheric waves in a thin circular layer in the form of the Fibonacci spirals
Ranis Ibragimov, Eleanor Blake
Wenatchee Valley College, Wenatchee, Washington, USA
Correspondence: Ranis Ibragimov, Wenatchee Valley College, Wenatchee, Washington, USA
Received: November 15, 2022 | Published: November 23, 2022
Citation: Ibragimov R, Blake E. Designing the similarity solution associated with nonlinear zonal atmospheric waves in a thin circular layer in the form of the Fibonacci spirals. J Aquac Mar Biol. 2022;11(3):125‒133. DOI: 10.15406/jamb.2022.11.00346
Download PDF
Abstract
It is shown that a non-uniform scaling transformation group represents a particular class of the exact solution of the nonlinear shallow water equations corresponding to the two-dimensional atmospheric waves within a circular geometry. The obtained solution has a Fibonacci spiral-like form. We have used the parameters of the model to associate to obtained solutions with the Fibonacci spiral. It is shown that the physically relevant part of the solution matches exactly the Fibonaccy spiral.
Keywords: atmospheric waves, Fibonacci sequence, shallow water theory, exact solutions
Introduction
Probably, the most known relationship that occurs in the smallest, to the largest objects in nature is the Fibonacci sequence of numbers
is defined by using the recursive relation
with the assumed values
. It can be observed in shells of snails, in the pattern found on the wings of dragonflies or galaxies, as shown in Figure 1. The symmetry and appearance of Fibonacci spiral patterns in nature attracted scientists, but a mathematical or physical explanation for their common occurrence in nature is still known.1
Figure 1 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
In this paper we observe some correlation between two-dimensional atmospheric equatorial waves and the Fibonacci spiral. It is known that equatorial waves play an important role in the climate system as well as in the weather and climate changes. Equatorial waves are being associated with large-scale perturbations of the atmospheric motion extending around an equator.2–6 Also, it is believed that the equatorial waves can be an important component of the long-term mean upwelling at the tropical tropopause, widespread changes in the climate system and in absorbing harmful solar ultraviolet radiation.3,7–13 Also, equatorial waves can be an important component in the controlling the stratospheric temperature via atmospheric radiative heating14 and in improving seasonal weather forecasting, as pointed out in.8,12,15–21
Various methods of method of visualizing the atmospheric equatorial waves are discussed in.5,6
In this paper, we aim to visualize atmospheric equatorial waves that are being described by the particular class of the solution that is deduced as the similarity solution of the nonlinear shallow water equations representing the free boundary of two-dimensional model describing the nonstationary motion of an incompressible perfect fluid propagating around a solid circle of a sufficiently large radius with the gravity directed to it’s center. The solution of the interest was obtained by means of Lie symmetries first reported in.22 Acquaintance with group analysis is important for constructing and investigating nonlinear mathematical models of natural and engineering problems. Numerous physical phenomena can be investigated using Lie symmetries to unearth various group invariant solutions and conservation laws that provide significant physical insight into the problem. Moreover, models of natural phenomena can often be described directly in group theoretic terms. Differential equations, conservation laws, solutions to boundary value problems, and so forth can be derived from the group invariance principle.
We consider a two dimensional motion of an incompressible perfect fluid which has a free boundary η and a solid bottom represented by a circle of radius . The fluid is propagating within a curricular domain bounded below by a sufficiently large solid circle and above by a free boundary. It is assumed that the motion is irrotational and the pressure on a free boundary is constant. It is also assumed that the distance from the origin of the circle to the unperturbed free boundary is small compared to the radius of the circle, as shown schematically in Figure 2.
Figure 2 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
The model is described by means of polar coordinates
,
, in which
is a polar angle,
is the distance from the origin of the circle and
is the radius of the circle representing the radius of a planet,
, where
is time,
is a constant undisturbed level of the equatorial atmospheric layer from the center of the planet, and
is the unknown level of disturbance of a free boundary. Hence, irrotational motion of a perfect fluid is confined in the domain
, which is bounded by a solid circular boundary
and a perturbed free boundary
. As mentioned above, it s assumed that the unperturbed level of atmospheric level
is small compared to the radius and the homogeneous gravity field is given by the vector
which is assumed to be a constant and directed to the center of the circle.
Shallow water approximation: In what follows, it is assumed that the fluid motion is potential in the domain of the motion which allows to introduce the stream function
via
(1)
Notice that
is a harmonic function in
, since we assumed that the flow is irrotational. Next, we define the average velocity
in terms of the stream function
as
(2)
In order to reduce the number of parameters, we introduce dimensionless variables:
,(3)
We next introduce the parameter
(4)
The hypothesis that water is "shallow" is based on the assumption that the parameter is "small".
In particular the relation between the boundaries of the dimensional and nondimensional variables can be summarized in the Table 1 below Equ_(5)
Dimensional variables
|
Relation
|
Non dimensional variables
|
Solid boundary:
|
|
|
Unperturbed free surface:
|
|
|
Perturbed free surface:
|
|
|
Hereafter the symbol "hat" is omitted.
So, in the nondimensional variables, the free boundary is given by
and thus the functions
and
are two unknown functions whereas is a given parameter. Although shallow water theory is usually related to the case when the water depth is small relative to the wavelengths of the waves, we find it more appropriate to choose the radius as a natural physical scale since, in the frame of the present model, we consider waves with wavelengths of the order of the radius of a planet.
Representing the stream function
by the series expansion:
and application of the Lagrange’s method allows one to write the boundary problem as the following nondimensional system of nonlinear shallow water equations, which are a higher-order analogue of the Su - Gardner equations:23
,(6)
,(7)
where the subscripts denote partial derivatives in which independent variables
and
denote the time and the polar angle, respectively, the dependent variables are the average velocity
and the level
of the atmosphere perturbed from
, whereas
is a small parameter.
Due to the fact that
is a small parameter, the terms of order
can be considered as small perturbation to the zeroth order terms (unperturbed model). Our main concern is a simplified version of the model, in which the perturbations (nonlinear terms at
, are ignored and we consider the following unperturbed system. As has been shown in,24,25 elimination of
and
from the terms of the equations (6)-(7) and ignore the terms with
. This leads to the following unperturbed system:
,(8)
.(9)
It can be checked by direct substitution that one particular exact solution of the unperturbed shallow water system (8)-(9) is
(10)
where
is the intensity of the vortex (source) localized at the center of the planet and is related with the the rotation rate of the planet (for example, the angular velocity of the earth is
by the equation
Since the solution (10) corresponds to the constant flow with an undisturbed circular free surface, we call it a trivial solution. At certain extent, the above ansatz (10) can be associated with the polar vortex, which represents a very powerful whirlpool swirling steadily around the planet’s poles at all times.
Understanding and predicting break downs and overall dynamic structure of the polar vortex is important for improving seasonal forecasting.21 If the undulating path of the west to east atmospheric flow generated by the polar vortex could be predicted, then weather could be predicted too – not just for a week or two, but for an entire season.26 However, as pointed out in,27 until now understanding of a such correlation has been based on observations and statistical modeling only rather than on the knowledge of its physical foundation. Because of the periodic seasonal breakdowns of the polar vortex and the lack of continuous data for the input of such interactions, statistical approach generally cannot provide the weather patterns accurately enough.28,29
Up to the present days, our knowledge of gravity wave sources and properties in the polar region is very much limited because collecting the observations is generally difficult because of harsh natural environments.30
Nontrivial similarity solution
Detailed presentations of the theory of symmetries and invariant solutions of differential equations can be found elsewhere.31–38 For convenience, we summarize the basic notation from calculus of Lie group analysis in Appendix, which represents a simplified version of the overview of basic concepts of Lie symmetry groups.
A simple inspection shows that the system (8)-(9) admits the infinite-dimensional Lie algebra composed by the operators
(11)
In particular, the operator
is admitted due to the invariance of the system (8) - (9) under the Galilean transformation group
(12)
With an arbitrary parameter
. The operator
indicates that the system (8) - (9) is invariant under the non-uniform scaling transformation group.
(13)
where
is another arbitrary parameter and the operator is responsible for a specific non-scaling symmetry of the one-dimensional shallow water model (8) - (9).
The term similarity solution refers to invariant solutions based on scaling transformations (12). It can be checked by direct substitution that the shallow water model (8)-(9) is invariant under the transformation (12). One calculates the invariant
of the group
by solving the first-order linear partial differential equation
(14)
The latter equation (14) gives three functionally independent invariants corresponding to the transformation (12):
Accordingly, we look for the invariant solution of the (8) - (9) in the form
(15)
Direct substitution of the presentation (15) in the system (8) - (9) yields the following nonlinear ordinary differential equations:
(16)
(17)
In case when
the equation (17) can be written as the following coupled equations
(18)
where
satisfies the nonlinear differential equation of the second order,
(19)
Integration of the equation (19) yields
(20)
where
is a constant and so
is given implicitly by equation
(21)
where
is an arbitrary constant.
In particular, if the function
is known, the function
can be expressed in terms of
as
(22)
The equation (21) provides an implicit representation of the function
and hence the functions
and
is (15) due to the equations (18). We calculate the integral and distinguish the following three cases:
Case k=0
When k=0 the equation (21) has the form
(23)
In this case, the equations (15) provide us with the solution
(24)
and thus the exact solution of the system (8) - (9) is given by
(25)
Case k>0
When
, the functions
and
are given by
(26)
where
and
is given implicitly by equation
(27)
Case k<0
When
the functions
and
are given by
(28)
where
and
is given implicitly by equation
(29)
The function
is defined subject to constraint
(30)
For example, the function
is not defined for
We will approximate the functions
for large and small values of
and also for large and small values of
Asymptotic analysis
We will analyze the asymptotic behavior of the similarity solution (15) for large and small
as well as for large and small values
Approximation of the similarity solution for large and small W
We next approximate the function
given by equation (27) in the limit
We start with
where
(31)
First, we consider the limiting behavior of
(32)
and for
since
we look for
Using
we find
(33)
After applying the L’Hopital’s rule 3 times, we get
(34)
Next, with the limiting behavior of
(35)
For
since
we get
(36)
Using this, we can find two approximations for
(37)
since
(38)
(39)
since
(40)
This analysis shows that there are two approximations for
for very large
and for very small
In particular,
but
does not exist since
is subject to the constraint
(Figure 3).
Figure 3 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
Approximation of the similarity solution for large and small
Approximation of
by
and
for larger and small values of
We observe that
as illustrated in panel (a) and
as illustrated in panel (b). Panel (c) is used to compare the curves
and
with
for the value of
We note that the function
represents a better approximation of
for larger values of
whereas the function
represents a better approximation of
for smaller values values of
which is demonstrated in Figure ??. In particular, the panel (a) shows the curve
and
for
and
It can be shown that
The panel (b) compares the curve
and v for
and
It can be shown that
subject to the constraint (30). Finally, the panel (c) is used to compare
with
and
for
For this reason we call the exact solution (15) an attractor solution in the sense that
is defined as the smallest unit which cannot be itself decomposed into two or more attractors with distinct basins of attraction. This restriction is necessary since, in general, a dynamical system may have multiple attractors, each with its own basin of attraction.
Nonlinear Analysis
We start with analyzing the exact attractor solution
representing the deviation of the free boundary (see equation (15)) form the unperturbed boundary
for the case
at different values of time
As shown in Figure 4, the attractor solution is decreasing function of
for a fixed values of
but it is increasing function of
for a fixed values of time
We will also apply the numerical technique to compare the attractor solution
Figure 4 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
(41)
with the exact solution
(42)
the both cases when
and
The results are shown in Figures 5,6. Attractor solution
for the case
and differentt valaues of time
and
Comparison of the solution
for
with the attractor solution
at
and different values of
Panel (a) shows the attractor solution (red line) and the solution
evaluated at
(blue line) and
evaluated at
(green line) and panel (a) compares these solutions for the
and 1 as well, but evaluated at
In particular, the Figure 5 is used to compare the solution
given by (42) for
with the attractor solution given by (41) at
and different values of
Panel (a) shows the attractor solution (red line) and the solution
given by (41) evaluated at
(blue line) and
(green line). We can observe that the attractor solution represents a betterb approximation of the exact solution given by (41) at
then at
as expected from the asymptotic analysis showing that
A similar conclusion holds for larger values of time t, as illustrated in Panel (b) of this figure, showing the solutions evaluated at
Figure 5 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
Similarly, the Figure 6 is used to compre the nonlinear solution
given by (42) for
with the attractor solution
given by (41) at
and different values of
In this plot, the attractor solution is plotted by the red line and the solution
evaluated at
is shown by the blue line where as
evaluated at
is shown by the green line. We note that the similarity solution
is defined for
and
, which covers the whole range of
Figure 6 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
Fibonacci spirals
We associate the parameter k with the Fibonacci sequence
defined by
with
and
Here each value
corresponds to the circular arc
where
so that
corresponds to the circular arc
with
Also
corresponds to the circular arc
with
and
corresponds to the circular arc
with
and so on. Next, the parameter
is chosen for each value of
such that the endpoints of the arc segments
connect smoothly with the next circular segment
for any given value of time
so that
where we call
model hypermarameter and it is tuned for the given predictive model represented by the Fibonacci spiral. We are tuning the hypermarameter in order to discover how the parameters
and
of the model result in the given prediction.
For example, for
the values of
corresponding to the first seven terms of Fibonacci sequence
are approximated by the values shown in Table 2 Equ_(43).
The values of
are also shown in Figure 7 and the corresponding few circular arcs are shown in Figure 8. We observe that
(44)
which is also seen from the Figure 7.
Figure 7 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
In particular, based on the results for the hyperparameter
as shown in Table 2, we visualize our exact solution (15) for
for
in Figure 8 as the sequence of circular arcs plotted at the values of time given by
. For example, when
we evaluate our exact solution as a function of
at the value of time
Similarly, when
we evaluate our exact solution as a function of
at the value of time
and so forth.
Figure 8 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
Figure 9 is used to compare the attractor solution
with the numerical approximation of the perturbation
for
and
and
as the function of time
We note that
represents the deviation of the free boundary from the unperturbed state
so the feasible domain of
is given by
(45)
Figure 9 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
where
is a small parameter since it is assumed that the unperturbed level of atmospheric “depth” is small compared to the radius of the Earth. As Figure ?? shows,
is decreasing function of time
but it is increasing function of a polar angle
which means that our nontrivial solution is valid for only very small values of
Comparison of the free boundary
and the Fibonacci spiral for the first three values of
Figure 10 is used to compare the free boundary given by
and the Fibonacci spiral for the first four values of
namely, for
According to the values chosen as shown in Table 2, panel (a) shows the segment of the Fibonacci spiral (green circles) and the free boundary (blue solid line) evaluated at the values of
and
Panel (b) compares the segment of the Fibonacci spiral and the free boundary evaluated at the values of
and
Panel (c) compares the segment of the Fibonacci spiral and the free boundary evaluated at the values of
and
Finally, panel (d) compares the segment of the Fibonacci spiral and the free boundary evaluated at the values of
and
Figure 10 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
As it was discussed earlier, the solution for the free boundary perturbation
is only valid in the domain
defined by (45), i.e. in a very narrow band
For example, if
we can see from Figure (??) that our solution is valid for a very small values of
bounded by
Thus, for example, if
the solution is bounded by
So, in this domain, the Fibonacci spiral corresponding to
is a very good approximation of the nontrivial solution
as shown in panel (a) of the Figure (??). For larger values of
we observe that the solution
diverges from the Fibonacci spiral.
Discussion
In this paper, we have analyzed and visualized the exact invariant solution of the nonlinear simplified version of the shallow water equations (8) - (9)
(46)
which are being used to simulate equatorial atmospheric waves of planetary scales. Our model is represented by the Cauchy–Poisson free boundary problem on the nonstationary motion of a perfect incompressible fluid circulating around a vortex field approximated by a circle of a large radius and the gravity is directed to the center of the circle. The solution for the free boundary perturbation
is only valid in a very narrow band
as also shown schematically in Figure 12. We have shown that within this band, Fibonacci spiral corresponding to
is a very good approximation of the nontrivial solution
as shown in Figure 11. However, we observe that the solution
diverges from the Fibonacci spiral for increasing values of
. In other words, the physically relevant part of the solution matches exactly the Fibonaccy spiral.
Figure 11 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
Figure 12 N=57; Epidemiological distribution of the pathological fractures, traumatic fractures, and nonunion.
We also remark that the higher order of the shallow water equations can be derived for our free boundary model and the higher-order approximation has the following form:
(47)
(48)
It is one of our goals for further studies to study the higher-order approximation (47) - (48) and see if the Fibonacci spiral can be a better approximation of the free boundary model for larger values of
We also note that the model (8) - (9) is invariant under the group of transformations
(49)
where the functions
solve the system of first-order equations
(50)
If we write the operator (49) as an infinitesimal transformation
(51)
The equations (50) show that the infinitesimal transformation (51) changes the variables
by adding to them an arbitrary solution of the system of linear differential equations
(52)
and that the system (52) admits the infinitesimal transformation (51).
Hence the operator (49) is admitted both by the nonlinear system (8) - (9) and by the linear system (52). This fact predicts a possibility to map the nonlinear equations (8) - (9) to the linear system (52) by an appropriate change of variables. It is important to note that the linear system is homogeneous, i.e. invariant under the uniform dilation
produced by the operator
from (11).
We conclude that the operator (49) and the operator
from (11) are responsible for the possibility of mapping of the nonlinear system (8) - (9) to the linear homogeneous system (52). Mapping the nonlinear model (8) - (9) to a linear system will be studies in the forthcoming paper.
Acknowledgments
We would like to thank the Dean of STEM & Social Sciences at Wenatchee Valley College, Holly Bringman and Electronic Media Producer Lead, Brett McGinnis for providing us with all the necessary facilities for the research, great support, and encouragement.
Appendix: Outline of methods from Lie group analysis
Basic concepts from Lie group analysis of differential equations that are used in the present paper are assembled here.For further information regarding Lie groups and their applications to the theory of differential equations, the reader should consult the various classical and modern texts in the field, such as.31–38
Definition of one-parameter groups: Let
(A1)
be a one-parameter family of invertible transformations of points
into points
Here
is a real parameter from a neighborhood of
and we impose the condition that Transformation (A1) is an identity if and only if
i.e.,
(A2)
The set
of transformations (A1) satisfying Condition (A2) is called a (local) one-parameter group of transformations in
if the successive action of two transformations is identical to the action of a third transformation from
i.e., if the function
satisfies the following group property:
(A3)
where
b (A4)
With a smooth function
defined for sufficiently small
and
The group parameter
in the transformation (A1) can be changed so that the function (A4) becomes
In other words, the group property (A3) can be written, upon choosing an appropriate parameter
(called a canonical parameter) in the form
(A5)
Group Generator. Let
be a group of transformations (A1) satisfying the condition (A2) and the group property (A5). Expanding the functions
into Taylor series near
and keeping only the linear terms in
one obtains the infinitesimal transformation of the group
(A6)
where
(A7)
The first-order linear differential operator
(A8)
is known as the generator of the group
Invariants: A function
is said to be an invariant of the group
if for each point
is is constant along the trajectory determined by the totality of transformed points
The function
is an invariant of the group
with Generator (A8) if and only if
(A9)
Hence any one-parameter group has exactly
functionally independent invariants (basis of invariants). One can take them to be the left-hand sides of
first integrals
of the characteristic equations for linear partial differential equation (A9). Then any other invariant is a function of
Invariant equations: We say that a system of equations
(A10)
Is invariant with respect to the group
(or admits the group
) if the transformations (A1) of the group
map any solution of Eqs. (A10) into a solution of the same equations, i.e.,
(A11)
whenever
solves Eqs. (A10). The group
with the generator (A8) is admitted by Eqs. (A10) if and only if
(A12)
where the symbol
means evaluated on the solutions of Eqs. (A10).
If
is a collection of independent variables
dependent variables
and partial derivatives
of
with respect to
up to certain order, where
then (A10) is a system of partial differential equations
(A13)
Furthermore, if the transformations (A1) are obtained by the transformations of the independent and dependent variables
(A14)
and the extension of (A14) to all derivatives
etc. involved in the differential equations (A13), then Eqs. (A11) define a group
of transformations (A14) admitted by the differential equations (A13). In other words, an admitted group does not change the form of the system of differential equations (A13). The generator of the admitted group
is termed an infinitesimal symmetry (or simply symmetry) of the differential equations (A13). Equations (A12) serve for obtaining the infinitesimal symmetries and are known as the determining equations. These equations are linear and homogeneous and therefore the set
of its solutions is a vector space. integration of determining equations often provides several linearly independent infinitesimal symmetries. Moreover, the determining equations have a specific property that guarantees that the set
is closed with respect to the commutator
Due to this property
is called a Lie algebra. If the dimension of the vector space
is equal to
the space is denoted by
and is called an
dimensional Lie algebra. An
dimensional Lie algebra
generates a group depending on
parameters which is called an v parameter group.
Invariant solutions. Let the differential equations (A13) admit a multi-parameter group
and let
be a subgroup of
A solution
(A15)
of Eqs. (A13) is called an
invariant solution (termed for brevity an invariant solution) if Eqs. (A15) are invariant with respect to the subgroup
If
is a one-parameter group and has the generator
then the
invariant solutions are constructed by calculating a basis of invariants
39,40
References
- Li C, Ji A, Cao Z. Stressed Fibonacci spiral patterns of definite chirality. Applied Physics Letters. 2007;90:164102.
- Ibragimov R, Mohazzabi P, Roembke R, et al. Asymptotic stability of the polar vortex perturbed by harmonic waves describing atmospheric gravity waves circulating in an equatorial plane of a spherical planet. Math Model Nat Phenom. 2018;13(4).
- Holton J, Haynes P, McIntyre M, et al. Stratosphere–troposphere exchange. Rev Geophys. 1995;33(4):403–439.
- Newman P. Chemistry and dynamics of the Antarctic ozone hole, in the stratosphere: Dynamics, fransport, and Ghemistry. Geophys Monogr. 2010.
- Ibragimov R. Nonlinear viscous fluid patterns in a thin rotating spherical domain and applications. Phys Fluids. 2011:23(12):123102.
- Ibragimov R, Gunag L. Splitting phenomenon of a higher-order shallow water theory associated with a longitudinal planetary wave. Dynam Atmos Ocean. 2015;69:1–11.
- Yulaeva E, Wallace J. The signature of ENSO in global temperatures and pre- cipitation fields derived from the Microwave Sounding Unit. J Glimate. 1994;7(11):1719–1736.
- Zhang Y, Li J, Zhou J. The Relationship between Polar Vortex and Ozone Depletion in the Antarctic Stratosphere during the Period 1979–2016. Advances in Meteorology. 2017:1–12.
- Ibragimov R, Lauren D, Stimmel M, et al. Visualization of Exact Invariant Solutions Associated with Atmospheric Waves in a Thin Circular Layer. Journal of Applied Mathematics and Physics. 2003;9(5):1321–1328.
- Newman P. Chemistry and dynamics of the Antarctic ozone hole, in the stratosphere: Dynamics, fransport, and Ghemistry. Geophys Monogr. 2010.
- Shindell DT, Schmidt GA. Southern hemisphere climate response to ozone changes and greenhouse gas increases. Geophys Res Gett. 2004.
- Lubin D, Jensen E. Effects of clouds and stratospheric ozone depletion on ultraviolet radiation trends. Nature. 2002;377:710–713.
- Chipperfield M, Dhomse S, Feng W, et al. Quantifying the ozone and ultraviolet benefits already achieved by the Montreal Protocol. Nat Commun. 2015;6:7233.
- Baldwin M, Dunkerton T. Stratospheric Harbingers of Anomalous Weather Regimes. Science. 2001;294(5542):581–584.
- El Niño, La Niña. New Zealand’s National Institute of Water and Atmospheric Research. 2016.
- Becker J. How Much Do El Niño and La Niña Affect Our Weather? Scientific American. 2016.
- Angell J, Korshover J. Variation in size and location of the 300 mb north circumpolar vortex between 1963 and 1975. Monthly Weather Review. 1977;105(1)19–25.
- Sun L, Wu H, Li X. Visualization of Exact Invariant Solutions Associated with Atmospheric Waves in a Thin Circular Layer. Journal of Applied Mathematics and Physics. 2006;9(5):52–62.
- McCreary J. Eastern tropical ocean response to changing wind systems: with Applications to El Niño. J Phys Oceanogr. 1976;6(5)632–645
- Murray F. Dynamic stability in the stratosphere. J Geophys Res. 1960;65(10):3273–3305.
- Alexander M, Tsuda T, Vincent R. Latitudinal Variations Observed in Gravity Waves with Short Vertical Wavelengths. J Atmos Sci. 2002;59:1394–1404.
- Dickson LE. Differential equations from the group standpoint. Ann Math. 1924;25(4):287–378.
- Vassada A, Horst S, Kennedy M, et al. Cassini imaging of Saturn: Southern hemisphere winds and vorticities. J Geophys Res. 2006;111:5004–5017.
- Su C, Gardner C. Derivation of the Korteweg-de Vries equations and Burgers equations. J Math Phys. 1968;10(3):536–539.
- Schoeberl M, Hartmann D. The dynamics of the stratospheric polar vortex and its relation to springtime ozone depletions. Science. 1991;251:46–52.
- Yang X. Variability of the northern circumpolar vortex and its association with climate anomaly in China. Acta Meteorologica Sinica. 2016;26(2):135–142.
- Sigmond M, Scinocca J, Kharin V, et al. Enhanced seasonal forecast skill following stratospheric sudden warmings. Nat Geosci. 2013;6:98–102.
- Shi C, Xu T, Guo D, et al. Modulating effects of planetary wave 3 on a stratospheric sudden warming event in 2005. J Atmos Sci. 2017;74(5):1549–1559.
- Shi C, Xu T, Li H, et al. 2017: The Role of Rossby-Wave Propagation in a North American Extreme Cold Event. Advances in Meteorology. 2017: 4635849.
- Cohen A. An Introduction to the Lie Theory of One-Parameter Groups With Applications to the Solution of Differential Equations. DC Heath. New York; 1911.
- Dickson LE. Differential equations from the group standpoint. Ann Math. 1924;25(4):287–378.
- Ovsyannikov LV. Group Properties of Differential Equations, Siberian Branch, USSR Academy of Sciences. Novosibirsk. 1962.
- Olver PJ. Applications of Lie Groups to Differential Equations. Springer-Verlag: New York; 1986.
- Bluman G, Kumei S. On invariance properties of the wave equations. J Math Phys. 1987;28(2):307.
- Stephani H. Differential Equations: Their Solution Using Symmetries. Cambridge Univ Press: Cambridge; 1989.
- Cantwell BJ. Introduction to Symmetry Analysis. Cambridge University Press: Cambridge; 2002.
- Bluman GW, Anco SC. Symmetry and Integration Methods for Differential Equations. Springer-Verlag: New York; 2002.
- Kovalev V, Pustovalov V, Senashov S. Lie-Bäcklund group for the nonlinear geometrical optics equations. Differential equations. 1993;29(10):1751–1764.
- Yulaeva E, Wallace J. The signature of ENSO in global temperatures and pre- cipitation fields derived from the Microwave Sounding Unit. J Glimate. 1994;7(11):1719–1736.
- Zhang Y, Li J, Zhou J. The Relationship between Polar Vortex and Ozone Depletion in the Antarctic Stratosphere during the Period 1979–2016. Advances in Meteorology. 2017:1–12.
©2022 Ibragimov, et al. This is an open access article distributed under the terms of the,
which
permits unrestricted use, distribution, and build upon your work non-commercially.