Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The shallow water equation and the vorticity equation for a change in height of the topography

  • ChaoJiu Da,

    Affiliations School of Mathematics and Computer Science Institute, Northwest University for Nationalities, Lanzhou, Gansu, China, College of Atmospheric Sciences, Lanzhou University, Lanzhou, Gansu, China

  • BingLu Shen,

    Affiliation College of Atmospheric Sciences, Chengdu University of Information Technology, Chengdu, Sichuan, China

  • PengCheng Yan,

    Affiliation Institute of Arid Meteorology of CMA, Lanzhou, Gansu, China

  • DeShan Ma,

    Affiliation School of Mathematics and Computer Science Institute, Northwest University for Nationalities, Lanzhou, Gansu, China

  • Jian Song

    dmwsj@163.com

    Affiliation College of Science, Inner Mongolia University of Technology, Hohhot, Inner Mongolia, China

Abstract

We consider the shallow water equation and the vorticity equations for a variable height of topography. On the assumptions that the atmosphere is incompressible and a constant density, we simplify the coupled dynamic equations. The change in topographic height is handled as the sum of the inherent and changing topography using the perturbation method, together with appropriate boundary conditions of the atmosphere, to obtain the relationship between the relative height of the flow, the inherent topography and the changing topography. We generalize the conservation of the function of relative position, and quantify the relationship between the height of the topography and the relative position of a fluid element. If the height of the topography increases (decreases), the relative position of a fluid element descends (ascends). On this basis, we also study the relationship between the vorticity and the topography to find the vorticity decreasing (increasing) for an increasing (decreasing) height of the topography.

Introduction

Atmospheric motion is governed by laws of physics described by partial differential equations, which control changes of motion and the thermal state of the atmosphere. Initial and boundary conditions are required to solve these partial differential equations. Initial conditions are defined based on the initial flow distribution, while the boundary conditions are distinguished by the inner and outer boundary conditions. For the motion of the atmosphere, only the lower boundary conditions are required, which is taken as the surface of the topography. As a result of natural variability or human factors, the height of the topography may change, which can alter the variation in atmospheric circulation, the local climate, or even the global climate. The change of the height of the topography may reflect the effect of human activities on the atmospheric motion on a global scale. We focus here on the change in height of the topography, and its effect on the generalized dynamic processes of the atmosphere.

Liao et al. studied the influence of the topography on rainfall, concluding that the change in the height of the topography not only has an effect on the size and distribution of the horizontal and vertical flow fields of the atmosphere, but also alters atmospheric microphysical process, which affects the local climate [1]. The change in height of the topography can also lead to the occurrence of disastrous weather locally [2]. Researchers have long realized that the dynamic uplift of the topography on warm and humid air causes heavy rain [3,4]. For example, the change in the height of the topography is inversely correlated with the precipitation intensity of a ridge [5]. Because atmospheric flow is a highly nonlinear system, the small-scale variation of the height of the topography may lead to nonlinear instability, Duan et al. did a systematic study on the nonlinear instability of atmospheric motion [6,7]. Nonlinearity is the source of the chaotic behavior of the atmosphere, which mathematical theory attributes to instabilities in the solutions to the governing differential equations, leading to difficulties in numerical weather prediction, for the nonlinearity, Sun et al. did the mathematics and physics fundamental research, and gained a series of conclusion [813], the research result can partly explain nonlinear effect of the differential equations. Feng GuoLin and Huang JianPing et al. used historical observational data to improve the results of numerical simulations, and possible physical mechanism of water vapor transport [1419]. The change in the topography also impacts the forecasts of numerical models [20]. For a slowly changing height of the topography, Da ChaoJiu et al. studied the shallow water equations when the local horizontal divergence changes slowly, modifying the equilibrium between the local horizontal divergence and the local change in the thickness of the atmosphere [2123].

The influence of the topography on the motion of atmosphere is mainly thermal and dynamic effect, the former is mainly the action on the temperature and humidity of the atmosphere, the latter is mainly the action on the velocity and pressure. In this paper, we only do the research on the change of the topography height.

The shallow water equation for a change of the topography

For the barotropic flow, the simplest model is the shallow water model. We assume that the whole atmosphere is one layer of constant density, incompressible, satisfies the static approximation, together with the β-plane approximation. Fig 1 gives the conceptual diagram showing the change in the height of the topography. For a single layer, uniform-density fluid (ρ = constant) above the plane z = 0, the height of the fluid surface is h(x,y,t). In view of the effect of the earth on the atmosphere and ocean, the force of the potential function Φ can be thought of as the vector perpendicular to the plane z = 0. Here, the rotation of the Earth is parallel to the z-coordinate axis. The topographical boundary has two parts, one being the inherent underlying surface of the atmosphere or the ocean hB(x,y,t), the other being the part that slowly changes in time hS(x,y,t), written as (1)

thumbnail
Fig 1. Conceptual diagram showing the change in height of the topography.

https://doi.org/10.1371/journal.pone.0178184.g001

Here, the parameter σ takes a small value, which represents the gradual change in topography. This is based on the requirement of the perturbation method, and also meet the slowly change character of the topography, and can describe the real motion of the atmosphere. The velocity components parallel to the coordinate axes x, y and z are u, v and w, respectively, assuming furthermore the fluid is inviscid [24,25].

From the assumptions of incompressibility and a constant density, the continuity equation can be simplified to (2A) and the horizontal momentum equations are (2B) and (2C)

Here, define the depth of fluid as [24,25] (3)

Integrating the continuity Eq (2A) with respect to z above the interval [HB(x,y,t),h], while noting the normal velocity on the rigid surface z = HB(x,y,t) is zero [26], and the kinematic and the dynamic conditions on the free surface z = h(x,y,t), we obtain (4) which is consistent with the shallow water equations presented formally in reference [25]. The difference is that the function H in reference [25] lacks the temporal variation term σhS(x,y,t) appearing in Eq (4) for consideration of the evolution of the height of the topography. Eqs (2) and (4) are our modified shallow water equations. The unknowns u, v and h are associated with the parameters defining σ, which is relatively small. Supposing that the variable u is the sum of two terms (5A) then the first term u0(x,y,t) is independent of the change of the topography, while the second term u1(x,y,t) is associated with the change of the topography. Using the parameter σ to describe the influence of the topography, the variables v, w and h can similarly be written as (5B) (5C) and (5D)

As the horizontal velocity is bigger than the vertical velocity for two or three orders of magnitude for the motion of the atmosphere, the horizontal velocity expand to order σ, while the vertical velocity to σ2, this will help to ensure high accuracy. Inserting Eq (5) into the shallow water Eqs (2) and (4), and expanding to obtain terms of O(σ0), gives (6A) (6B) and (6C)

Invoking a zero normal-velocity component on the rigid base to obtain the vertical velocity w(x,y,t,σ), where w0 is the term of O(σ0), gives (6D)

Using Eq (6C) to eliminate in Eq (6D) results in (7)

In this way, the function w0 can be written as a function of the topography, where the term has the form (8) which can also be interpreted as the zero-order expansion of the full derivative.

The O(σ1) equations are (9A) (9B) (9C) and (9D)

From Eq (6C), we have (10) which is substituted into Eq (9C) to give (11) where (12)

This equation can also be interpreted as the first-order expansion of the full derivative, where is defined similarly to Eq (8), but is not repeated here. Substituting Eq (10) and Eq (11) into Eq (9D) gives (13) which is w1 expressed in terms of the topography, similarly to Eq (7).

The O(σ2) equations are (14A) (14B) (14C) and (14D)

Substituting Eq (11) into Eq (14D) gives (15)

Adding Eq (7) to Eq (13), multiplying by σ, and then adding to Eq (15) multiplied by σ2 gives (16)

Hence, the relationship between the vertical position z, the topography and the horizontal velocity is now given. If the topography is constant, that is to say σ = 0, Eq (16) can be written as (17) or (18) which illustrates that the relative-position function regulated the motion of each fluid element conserves in the shallow water, which is a well-known condition. The function is the relative height of the fluid element from the bottom boundary, clearly being zero on the bottom boundary (for z = hB) and one on the free surface (for z = h). Meanwhile, the horizontal velocity is independent of z, meaning that the flow moves parallel to z. Moreover, the relative position of a fluid element within a column remains unchanged for a lengthening or shortening column [24,25]. If the height of the topography changes with time, Eq (16) gives the relationship between the position z of the fluid element, the topography and the velocity, and is complex and highly nonlinear. The first line on the right-hand side of Eq (16) is the zero-order approximation, which lacks the changing of the topography, with only the inherent topography being present. The second and third lines on the right-hand side of Eq (16) correspond to the first-order approximation, where a change in the topography is accounted for through hS, including the nonlinear interaction between hS and hB. The fourth and fifth lines on the right-hand side of Eq (16) represent the second-order approximation, which contains the topography term hS and nonlinear interaction between hS and hB. In general, the change of the topography is small, because the parameter characterizing the degree of the change of the topography σ is small. For these reasons, Eq (16) can be simplified by omitting the second-order terms. As σhS in is also small, this term may also be omitted. As the term is the rate of change of slowly varying topography hS with respect to the longitude x, can also be omitted; and the term is the rate of change of hS with respect to the latitude y, can also be omitted. Hence, Eq (16) may be simplified to (19) by which (5A) and (5B) give (20) for simplifying Eq (19) to (21) or (22A) so that Eq (22A) may be written as (22B) where the differential operator in Eq (22) is defined after Eq (8). In comparing Eq (22) with Eq (18), we see that the relative-position function is not conserved. The change of the topography through the relative position function illustrates that if the topography hS changes with the time, this function will adjust to adapt to the change of the topography. Therefore, Eq (22) indicates the motion of the atmosphere more accurately than (18). To help understand Eq (22), we study one special case by supposing u0, v0 and h0 in Eq (6) have the solutions (23) where D is constant, and Eq (23) is a specific solution of Eq (6). This means the atmosphere in the horizontal direction is stationary, and the topography is quasi horizontal, the case is simple, but has the characteristics of the motion of the atmosphere. With this assumption, Eq (9) can be simplified to (24A) (24B) (24C) and (24D)

Eq (22B) also simplified to (25)

As the inherent topography hB is time-independent, Eq (25) has the new form (26) where DhB > 0. If , this indicates hS increases in time, meaning that the height of the topography increases. Thus, the term is negative, implying that for the stretching or contracting of a fluid element, the relative position of the fluid element in the fluid column descends to balance the increasing height of the topography. Conversely, if , and hS decreases with time, the height of the topography decreases for positive . Hence, the relative position of a fluid element within a fluid column must ascend to balance a decline in the height of the topography.

The vorticity equation for a change of the topography

The components of vorticity are written as (27A) (27B) and (27C)

For the shallow water model, u and v are z-independent, while a scale analysis tells us that the horizontal component of relative vorticity is smaller than the vertical component. To study the vertical component, we write (28) and only consider the horizontal momentum Eq (2) to obtain (29) which defines the behavior of the relative vorticity with respect to the divergence of the absolute vortex filament [24, 25]. The deformation of Eq (4) is (30) which, with the help of Eqs (29) and (30), we eliminate to obtain (31)

Eqs (31) is the conservative laws on vorticity, this indicates the vorticity is conserved in the case of the changing topography. Finally, inserting Eqs (3) and (28) into Eq (31) gives (32)

Hence, if the height of the topography increases (decreases), meaning that hS increases (decreases) in time, the relative vorticity ζ must decrease (increase).

Conclusions

Under the influence of global warming and other damaging human activities, the height of the topography may change, which will have an impact on the motion of the atmosphere. A method for quantitatively describing this impact is given here. In targeting the change of the topography in time with the help of the perturbation method, together with the pertinent boundary conditions, we quantitatively study the changing of the relative-position function and the relative vorticity ζ for a changing topography. We obtain the following conclusions: 1. For a changing topography, the relative-position function is not conserved and must therefore adjust to the topography; 2. In cases for which the height of the topography increases, the relative position of a fluid element in a fluid column must descend; conversely, if the height of the topography decreases, the relative position of a fluid element in a fluid column must ascend; 3. If the height of the topography increases (decreases), the vorticity decreases (increases).

Author Contributions

  1. Conceptualization: CJD.
  2. Formal analysis: CJD.
  3. Funding acquisition: CJD.
  4. Methodology: CJD.
  5. Project administration: CJD.
  6. Resources: CJD.
  7. Software: BLS.
  8. Supervision: DSM JS.
  9. Validation: BLS PCY.
  10. Visualization: BLS.
  11. Writing – original draft: BLS PCY.
  12. Writing – review & editing: BLS PCY.

References

  1. 1. Fei Liao, YaMin Hu, YanChao Hong. Numerical Study for the Influences of Orographic Dynamic on Cloud and Precipitation in North China. Plateau Meteorology, 2009, 28(1): 115–126
  2. 2. Hui Xu, RongHua Jin. Analyese of Influence of Terrain on Freezing-Rain Weather in Hunan in the early 2008. Plateau Meteorology, 2010, 29(4): 957–967
  3. 3. JiPing Chao. On the dynamics of orographically produced finite perturbations in a baroclinic westerlies. Acta Meteor Sinica, 1957, 28(4): 59–70
  4. 4. DuZhen Ye, YouXi Gao. Qinghai-Xizang plateau meteorology. Beijing: MingWen Press, 1966, 1–9
  5. 5. BaoGui Bi, YueWei Liu, ZeChun Li. Study on Influence of the Mechanical forcing of mesoscale Topography on the Extremely heavy Rainfall in Southern Shanxi on 8–9 June. Plateau Meteorology, 2006, 25(3): 485–494
  6. 6. Mu Mu. Some Advances In the Study of the Nonlinear Instability of the Atmospheric Motions. Scientia Atmospherica Sinica, 1995, 19(4): 494–509
  7. 7. Mu Mu, BoYu Chen, FeiFan Zhou, YanShan Yu. Methods and Uncertainties Of Meteorological Forecast. Meteorological Monthly, 2011, 37(1): 1–13
  8. 8. Sun Gui-Quan, Wang Cui-Hua, Wu Ze-Yan. Pattern dynamics of a Gierer–Meinhardt model with spatial effects. Nonlinear Dynamics, 2017, 88(2): 1385–1396
  9. 9. Sun Gui-Quan. Mathematical modeling of population dynamics with Allee effect. Nonlinear Dynamics, 2016, 85 (5): 1–12
  10. 10. Sun Gui-Quan, Wu Ze-Yan, Jin Zhen, Wang Zhen. Influence of isolation degree of spatial patterns on persistence of populations. Nonlinear Dynamics, 2016, 83 (1): 811–819
  11. 11. Li Li, Jin Zhen, Li Jing. Periodic solutions in a herbivore-plant system with time delay and spatial diffusion, Applied Mathematical Modelling 40 (2016) 4765–4777
  12. 12. Li Li, Patch invasion in a spatial epidemic model, Applied Mathematics and Computation 258 (2015) 342–349
  13. 13. Sun Gui-Quan, Zhang Juan, Song Li-Peng, Jin Zhen, Li Bai-Lian, Pattern formation of a spatial predator-prey system, Applied Mathematics and Computation 218 (2012) 11151–11162
  14. 14. Huang J P, Yi Y H, Wang S W. The experiment of seasonal prediction using the analogy-dynamical model. Science in China (B), 1992, 35 (2): 207–216
  15. 15. Wu YP, Feng GL, Li BL. Interactions of Multiple Atmospheric Circulation Drive the Drought in Tarim River Basin. Scientific Reports,2016, 6: 26470 pmid:27198665
  16. 16. Wu Y P, Feng GL. A new algorithm for seasonal precipitation forecast based on global atmospheric hydrological water budget, Applied Mathematics and Computation,2015, 268: 478–488
  17. 17. Guo-Lin Feng, and Wu Yong-Ping. Signal of Acceleration and Physical Mechanism of Water Cycle in Xinjiang, China. PloS one, 2016, 11: e0167387 pmid:27907078
  18. 18. Wu Yong-Ping, Shen Yong-ping, Li B. Larry. Possible Physical Mechanism of Water Vapor Transport over Tarim River Basin. Ecological Complexity, 2012, 9: 63–70
  19. 19. Feng G L, Sun S P, Zhao J H,Zheng Z H. Analysis of stable components for extended-range (10–30 days) weather forecast: A case study of continuous overcast-rainy process in early 2009 over the mid-lower reaches of the Yangtze River. Science China Earth Sciences, 2013, 56(9): 1576–1587
  20. 20. Bin Chen, YaoCun Zhang, YuGuo Ding. The Effect of Terrain Fluctuation on the Computed Surface Long-wave Radiation in the Models. Plateau Meteorology, 2006, 25(3): 406–412
  21. 21. Chaojiu Da, ShuPeng Sun, Jian Song, LianGui Yang. Variable CoefficientKdV Equation for the Amplitude of Nonlinear Solitary Rossby Waves In a Sort of Time-Dependent Zonal Flow. Plateau Meteorology, 2011, 30(2): 349–354
  22. 22. ChaoJiu Da, JiFan Chou. the KdV equation with a forcing term for the evolution of the Amplitude of Rossby Waves along a Slowly Changing Topography. Acta Phys. Sin. 2008, (57)4: 2595–2599
  23. 23. ChaoJiu Da, AiXia Feng, ZhiQiang Gong, jian Song. The dynamic modification of the shallow water equation on the slowly changing underlying surface condition. Acta Phys. Sin. 2013, (62)3: 039202
  24. 24. Bing Wang, HenYi Weng. Geophysical Fluid Dynamics. Beijing: China Ocean Press, 1981, 7–141
  25. 25. Pedlosky Joseph. Geophysical Fluid Dynamics. Springer Verlag, 1986, 7–9
  26. 26. BingRong Guo. The Introduction of Linear and nonlinear waves. Beijing: China Meteorological Press, 1990, 30–60