Abstract

The evolution of nonlinear gravity solitary waves in the atmosphere is related to the formation of severe weather. The nonlinear concentration of gravity solitary wave leads to energy accumulation, which further forms the disastrous weather phenomenon such as squall line. This paper theoretically proves that the formation of squall line in baroclinic nonstatic equilibrium atmosphere can be reduced to the fission process of algebraic gravity solitary waves described by the (2 + 1)-dimensional generalized Boussinesq-BO (B-BO) equation. Compared with previous models describing isolated waves, the Boussinesq-BO model can describe the propagation process of waves in two media, which is more suitable for actual atmospheric conditions. In order to explore more structural features of this solitary wave, the derived integer order model is transformed into the more practical time fractional-order model by using the variational method. By obtaining the exact solution and the conservation laws, the fission properties of algebraic gravity solitary waves are discussed. When the disturbance with limited width appears along the low-level jet stream, these solitary waves can be excited. When the disturbance intensity and width reach a certain value, solitary wave formation takes place, which is exactly the squall line or thunderstorm formation observed in the atmosphere.

1. Introduction

On June 4, 2016, a squall line weather occurred in northern Sichuan, China, lasting for 12 hours. On that day, from 09:00 to 20:00 Beijing time, a large area of catastrophic winds hit Aba, Guangyuan, Mianyang, Nanchong, and Bazhong in Sichuan province and the northern part of Chongqing province, with the maximum speed reaching 33.5 , causing serious loss of life and property. Among them, the squall line’s strong convection cell generated local thunderstorm strong wind which capsized Guangyuan city Bailong Lake cruise ship, killing 15 people. So, the research of squall lines has important theoretical meaning and real value. But, the research on the formation and prediction of squall line phenomenon is very difficult, which is a problem.

The phenomenon of gravitational waves occurring in the upper atmosphere is the most common and important dynamic process. They propagate in the middle atmosphere and generate energy and momentum transfer between atmospheres at different altitudes, which leads to the process of energy coupling between atmospheres [1]. The pressure jump, strong wind, strong convergence, and upward movement of gravity waves are caused by nonlinear action in the dispersion process. Mesoscale strong convection and rainstorm systems often have the characteristics of gravity wave [2]. Therefore, the gravity wave is an important factor in the formation of many mesoscale weather phenomena.

Many previous studies have shown that the squall line phenomenon is closely related to gravity solitary waves. In the 1960s, Long [3] proposed that the amplitude of atmospheric fluctuation satisfied the Korteweg–de Vries (KdV) equation in the first time. In 1978, Li [4] obtained the periodic and elliptic cosine solutions of the nonlinear atmospheric gravity waves, nonlinear Rossby waves, and elliptic cosine waves from the nonlinear atmospheric motion equations. Liu and Liu [5] found that the amplitude of gravity solitary waves was proportional to the propagation speed in 1983. Then, Bu et al. [6] and Wang [7] studied gravity waves by using the nonlinear Korteweg–de Vries (KdV) equation from the nonlinear atmospheric dynamics equation. In 2002, Lott and Plougonven [8] calculated the EP flux using the WKB approximation method and parameterized the gravity wave source. In recent years, with the rapid development of numerical models, many experts at home and abroad began to use weather models, cloud models, and climate models to research the development of gravity waves. Jewett et al. [9] and Plougonven and Snyder [10] successfully simulated storm, frontal, and gravity waves with mesoscale models MM5 and WRF, respectively. In 2011, Lane and Zhang [11] used a two-dimensional convective cloud model to study the mechanism of convection-induced gravity waves. Regarding the formation of squall line, Li and Xue [12] and Li [13] pointed out that the nonlinear concentration of gravity solitary waves was the formation mechanism of squall line. Yao et al. [14] also analyzed the formation and maintenance of strong squall line in detail in 2005. Then, in 2014, Srinivasan et al. [15] studied the gravitational wave characteristics of squall line during its propagation. Stephan and Alexander [16] discussed the relationship between the formation of summer squall line and gravity solitary waves. Recently, in 2000, Scinocca and Ford [17] found that large gravity solitary waves forced unstable layered shear flows.

It can be seen from previous studies that the research on gravity waves and squall line mainly focuses on numerical simulation, while the research on theoretical analysis is very few. We know that the basic dynamic equations describing the motion of baroclinic atmosphere are more and more complex [18, 19], and the gravity waves of baroclinic atmosphere [20, 21] are seldom studied. However, because of the baroclinic characteristics of the atmosphere, the baroclinic problem in the actual atmosphere is inevitable. In addition, these studies describe classical gravity solitary waves [22, 23]. In fact, some classical gravity solitary wave models, such as the Korteweg–de Vries (KdV) equation, modified Korteweg–de Vries (mKdV) equation, and Boussinesq equation, describe the propagation of gravity solitary waves in a certain direction [2427]. But, true gravity solitary waves travel in both directions [28]. Therefore, we need to establish a new model to replace the classical gravity solitary waves and explore the formation mechanism of squall line so as to better adapt to the actual atmospheric conditions.

In recent decades, fractional partial differential equations have been increasingly used to describe problems in optical and thermal systems, rheological and material and mechanical systems, signal processing and system identification, control and robotics, and other applications [2931]. Fractional-order equations are divided into time fractional equations and space fractional equations [3234] and have various definitions. Because fractional derivatives are historically dependent and nonlocal, they can more accurately describe complex physical and dynamic system processes in nature. The solution of fractional-order equation has attracted more and more attention of scholars, who have studied many effective solutions, such as the (G’/G)-expansion method [35], Khater method [36], sun-equation method [37], functional variable method [38], modified extended Tanh method [39], Godunov-type method [40], Lie group analysis method [41], and so on [4245].

2. Perturbation Expansion of Atmospheric Dynamics Equation Set

In order to obtain a new integer-order model to describe the evolution of algebraic gravity solitary waves, we start from the classical dynamic equation set and carry out perturbation expansion and perturbation analysis. A general gas is not under extreme conditions, and the properties of a real gas are very close to that of an ideal gas. In addition, just to simplify our calculation, the gases that we are dealing with here are all in ideal gas states, so they all satisfy the ideal gas equation, the Clapeyron–Mendeleev equation. The basic dynamic equation set of baroclinic nonstatic equilibrium atmosphere is transformed to the dimensionless form (Appendix A) as follows:

In the vicinity of the low-level jet, the basic flow tends to have strong horizontal shear, while when it is away from the low-level jet, the shear of the basic flow is small. Therefore, we can assume that the basic flow of the low-level jet is almost constant and the region is divided into two parts and ; the corresponding shear flow is written aswhere is a constant, and the boundary conditions can be written as

This is equivalent to assuming that the boundary of the atmosphere at is a rigid boundary, and the disturbance in the region is only generated by the disturbance in region , so the disturbance disappears at infinity.

First of all, the situation of domain is studied. Introducing the multiscale slow time and space transformations [46, 47],and supposing that , we have the following small parameter expansions:where α is the order of magnitude parameter of 1.

Substituting equations (4) and (5) into equation (1) yields

Taking the zero-order approximation of ε, we have

The above equation indicates that the basic flow is balanced. First of all, we take the first-order approximations of ε as

Suppose the equation has the following solutions in separate variable form:

Then, we plug these solutions in equation (8); an equation that only depends on can be obtained:where the specific form of is shown in Appendix B. In addition to that, the second-order approximations of ε are given as

In the same way, we assume that the separable solutions of equation (11) have the following form:

Substituting the solutions into equation (11) yields

We find that the governing model of gravity solitary waves cannot be derived, so we need to proceed to the higher-order equation of ε:

Eliminating of equation (14), we can obtain an equation about :where the specific forms of are shown in Appendix B.

Substituting equations (9) and (12) into equation (15) yields

Then, we have

Because the homogeneous part of equation (10) is the same as that of equation (16), we can obtainwhere the specific form of is shown in Appendix B.

Integrating equation (19) with respect to from 0 to yields

Secondly, we introduce the following transformation in the domain :and are rewritten as

Substituting (21) and (22) into (1), the following can be obtained:

When, from equation (23), we obtain a Laplace equation for :

Also, the boundary conditions are

Therefore, the solution of equation (24) under boundary conditions can be obtained:where is the principal value of Cauchy integral equation (26). Taking the derivative with respect to y for both sides of equation (26) yields

At the point , matching the value of zone and zone , we obtain

Applying equation (27), we obtain the following:where . Furthermore, we haveand the boundary conditions

In conclusion, substituting (31) and (32) into (20), a new equation is obtained as follows:where is the Hilbert transformation and the coefficients are seen in Appendix B.

The new equation is the (2 + 1)-dimensional generalized Boussinesq-BO equation. Previous researchers have studied the (1 + 1)-dimensional Boussinesq equation, which can only describe the propagation of algebraic gravity solitary waves on a line, and the results described are quite different from the actual environment. Our new model can describe the propagation of algebraic gravity solitary waves on a surface that is more in line with the actual natural atmospheric environment.

In equation (33), denotes the dispersion effect and and denote the nonlinear effects, which shows that the generalized Boussinesq-BO equation contains two kinds of processes of gravity solitary waves, namely, dispersion process and nonlinear process. The dispersion and nonlinear processes in barotropic atmosphere are the main causes of squall line formation. Therefore, the specific formation mechanism of squall line can be obtained by researching this new model.

3. The (2 + 1)-Dimensional Time-Fractional Generalized Boussinesq-BO Equation

In this section, we obtain the (2 + 1)-dimensional time-fractional generalized Boussinesq-BO equation by using the Agrawal method, semi-inverse method, and fractional variational principle. At the beginning, we introduce some definitions as follows.

Definition 1. (see [48]). The Riemann–Liouville fractional derivation of a function is defined by

Definition 2. (see [48]). The fractional integration by parts rule is defined as follows:where is the Riemann–Liouville fractional derivation.
According to the (2 + 1)-dimensional generalized Boussinesq-BO equation,It is worth noting that we cannot obtain the ideal form of variation by integral transformation of the (2 + 1)-dimensional time-fractional generalized Boussinesq-BO equation. Suppose , where refers to a potential function. The potential equation of the generalized Boussinesq-BO equation iswhere and are considered as fixed functions because of their specificity. Then, we give the function of the potential equation as follows:where are Lagrangian multipliers and can be determined later, and R and T are the boundary of space and the limit of time, respectively.
Using integration by parts for equation (38) and taking , equation (38) is rewritten aswhereThen, taking the variation of this function, we haveApplying the variation optimum condition, i.e., , we can obtain the following:Because equation (42) is equivalent to equation (37), the unknown constants becomeTherefore, the Lagrangian form of the integer-order generalized Boussinesq-BO equation is as follows:Similarly, the Lagrangian form of the (2 + 1)-dimensional time fractional generalized Boussinesq-BO equation is expressed aswhere is Riemann–Liouville fractional partial derivative of α with respect to t. Also, the function of the time-fractional generalized Boussinesq-BO equation can be given asUsing Agrawal’s method [49, 50], the variation of functional equation (46) yieldsApplying the fractional integration by parts (Definition 2) and considering , we can obtainObviously, optimizing the variation of the functional , i.e., , the Euler–Lagrange equation for the (2 + 1)-dimensional time-fractional generalized Boussinesq-BO equation can be expressed asSubstituting equation (45) into the Euler–Lagrange equation, we havewhere . Again, substituting the potential function leads toThis is the (2 + 1)-dimensional time-fractional generalized Boussinesq-BO equation, which is an extension of the (2 + 1)-dimensional integer-order Boussinesq-BO equation. The time-fractional model is nonlocal, which makes it more suitable than the integer-order model to describe the propagation process of algebraic gravity solitary waves in baroclinic atmosphere.

4. A Solution of the Fractional Generalized Boussinesq-BO Equation

In this section, the exact solution of the (2 + 1)-dimensional time fractional generalized Boussinesq-BO equation is obtained by applying the trial function method. We can explore some propagation characteristics of gravity solitary waves by calculation.

First, we transform the partial differential equation into an ordinary differential equation by assuming that u is a function of , equation (51) can be rewritten as

Integrating twice with respect to η and assuming the integral constant to be zero, we have

Then, the solution of equation (53) has the following form:

Substituting the solution into equation (53) yieldswhere the relation used is

And then, finally, by calculating, we can figure out

It is well known that the real gravity solitary waves are spread along the two directions. Thus, when c is positive, the gravity solitary waves propagate towards the right; when c is negative, the gravity solitary waves propagate in the opposite direction. This is a special feature of the Boussinesq-Bo equation compared with the BO equation. Compared with the solution of Camassa–Holm equation [28], the solution form of the Boussinesq-Bo equation is very different not only because of the difference in the nonlinear terms of the two equations but also because of the steps involved in solving them. The solution of the Boussinesq-Bo equation is simpler than the solution of the Camassa–Holme equation.

5. Fission of Waves and Formation of Squall Lines

In Section 4, a solution of the fractional generalized Boussinesq-BO equation is obtained. In addition to solving equations, the conservation law of equation is something that a dynamic equation must be explored because it plays an important role [51]. In this section, conservation laws of the (2 + 1)-dimensional time fractional generalized Boussinesq-BO equation are derived from the variational principle. Then, most importantly, we will combine the conservation laws and the exact solution to study the fission process of gravitational solitary waves and the formation mechanism of squall line.

The first step is to obtain the conservation law of the equation. Because the Lagrangian density of equation (51) does not exist, we introduce a new variable to construct an equal formula as follows:

Letting be a section, the first prolongation of ψ be denoted as

The Lagrange density of equation (58) is as follows:where

Combining the Lagrangian density, the action functional is given aswhere is an open set of X, represents the space of independent variables, and represents the space of dependent variables. Introducing a vector J yields

Transform a section to a family section that depends on the parameter κ [52]. Then, we calculate the variation for the action functional as follows:where

If , and are supported by compactness, then .

Assuming yields the first conservation law:where means the mass of the algebraic gravity solitary waves, and equation (66) shows that the mass of the algebraic gravity solitary waves is conserved.

Then, from , we can see from equation (64) that the variation τ gives the local energy conservation law:wherewhere the equations indicate that the energy of the algebraic gravity solitary waves is conserved.

Also, the variation ξ gives the local momentum conservation law:wherewhere the equations indicate that the momentum of the algebraic gravity solitary waves is conserved.

Second, combining the conservation laws and the exact solution that we have, we forecast the amplitude of gravity solitary waves emerging from a special type of initial value. The initial value is given as

Without loss of generality, assuming that the initial gravitational solitary waves gradually change to n solitary waves with different amplitudes and phases as follows:where

According to equations (52) and (71), we obtain

Then, at , applying the conservation laws, the conserved quantities , and can be calculated from equations (71) and (74) as

In a given period of time, the conserved quantities can be thought of as the sum of the conserved quantities related to each emerging gravity solitary wave. According to equations (58)–(62) and (72) and (74), we have

Equating corresponding quantities in equations (75) and (76), as well as taking the three gravity solitary waves as the object of study and assuming that , we obtain

Finally, according to the above calculation results, we can see that the two algebraic gravity solitary waves which are generated due to fission can be described clearly. We will demonstrate the fission process of algebraic gravity solitary waves by figures with (Figures 13).

As can be seen from Figure 1, a solitary wave described by the equation is propagating forward at a fixed wave speed. In this way, the solitary wave is characterized by a pulsing nature and stable propagation in the atmosphere. In fact, many of these waves have been observed on the low jet stream in recent years.

In the atmosphere, the velocity of algebraic gravity solitary waves depends on the amplitude of the waves. Under certain conditions, the larger the amplitude of the wave, the higher the propagation speed of the gravity isolated wave, which reflects the characteristics of the fast propagation of the gravity solitary waves with large amplitude. According to the conservation laws, the energy of a gravity solitary wave is conserved during its propagation. When the energy dispersion velocity of the gravity solitary waves is slower than the propagation velocity of the waves, the propagation velocity of the gravity solitary waves is faster for the long waves. For short waves, gravity waves travel more slowly or even backwards. When the nonlinear energy and wave energy dispersion reach a balance, the stable algebraic gravity solitary wave will be formed in the atmosphere.

As can be seen from Figure 2, with dispersion and nonlinear changes, when a disturbance occurs within a certain range, a single gravitational solitary wave will be excited and split into a number of stable solitary waves and oscillatory wave train. Because of the amplitude of the solitary wave moving block, after the interaction of the initial time, at some time T, the solitary wave of amplitude a will appear at the position of . Therefore, the gravity solitary waves are arranged into an array according to their intensity, that is, the strong waves are in front and the weak waves are behind, as shown in Figure 3.

The above discussion shows that the disturbance in a limited range in the atmosphere can form a solitary wave array composed of several solitary waves, and the increase of both the width and intensity of the initial disturbance can increase the number, intensity, and velocity of solitary waves in the solitary wave array. This is consistent with the squall line observed in the low-altitude air and the formation of thunderstorm group caused by the squall line.

All in all, the fission of algebraic gravity solitary waves will produce the squall line phenomenon, that is to say, the newly derived (2 + 1)-dimensional fractional generalized Boussinesq-BO equation can reasonably describe the generation and evolution of squall line.

6. Conclusions

In barotropic nonstatic equilibrium atmosphere, squall line evolution is essentially a kind of algebraic gravity solitary wave under the interaction of nonlinear process and dispersion process, which can be described by the Boussinesq-BO equation. For the first time, we derived the (2 + 1)-dimensional generalized Boussinesq-BO equation from the basic equation of atmospheric dynamics by using the multiscale analysis and perturbation method. Compared with the previous equations, the generalized Boussinesq-BO equation described the properties of solitary waves similar to those observed in the atmosphere.

In addition, we extended the generalized Boussinesq-BO equation to the fractional order. By solving the fractional generalized Boussinesq-BO equation and studying the conservation laws, we find that the squall line and thunderstorm formation observed in the atmosphere can be explained as a solitary wave formation occurred by the fission of gravity solitary waves excited by a disturbance source.

Appendix

A. Dimensionless Process of the Basic Dynamic Equation Set of Baroclinic Nonstatic Equilibrium Atmosphere

We expressed the pressure gradient force and gravity in the vertical direction as the sum of the disturbance pressure gradient and buoyancy and used the Boussinesq approximation. The dimensionless process of the basic dynamic equation set in the baroclinic nonstatic equilibrium atmosphere iswhere t is the time variable and u, , and are the three fluid velocity components in the x-, y-, and z-directions, respectively. denotes the atmospheric pressure, f denotes the Coriolis parameter, denotes the gravitational acceleration, and are the potential and density of environmental flow field, respectively, and .

The characteristic scales and dimensionless quantities related to variables are introduced as follows:where is the characteristic pressure on the ground, H is the height of homogeneous atmosphere, and and represent changes in pressure in the horizontal and vertical directions, respectively. The ones with a prime at the top right are dimensionless variables.

By substituting equation (A.2) into equation (A.1), dimensionless equations can be obtained as follows:

We assumed that , , and introduced a small parameter , where N is the Brunt–Vaisala frequency and . If the tip sign is omitted, equation (A.3) is changed into

B. The Concrete Forms of Complex Parameters

The specific form of is

The specific forms of are

The specific form of is

The coefficients of equation (33) are

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was supported by the National Natural Science Foundation of China (Grant no. 11762011) and the Natural Science Foundation of Inner Mongolia (Grant no. 2017MS0108).