Mathematical Modeling of the Influence of the Wind Field Structure in the Atmosphere on the Cloud Formation Processes

The state of the physics of convective clouds and cloud seeding is discussed briefly. It is noted that at the present time there is a transition from the stage of investigation of “elementary” processes in the clouds to the stage of studying the formation of macroand microstructural characteristics of clouds as a whole, taking into account their system properties. The main directions of the development of cloud physics at the upcoming stage of its development are discussed. The paper points out that one of these areas is the determination of the structure-forming factors for the clouds and the study of their influence on their formation and evolution. It is noted that one of such factors is the interaction of clouds with their surrounding atmosphere, and the main method of studying its role in the processes of cloud formation is mathematical modeling. A three-dimensional nonstationary model of convective clouds is presented with a detailed account of the processes of thermohydrodynamics and microphysics, which is used for research. The results of modeling the influence of the wind field structure in the atmosphere on the formation and evolution of clouds are presented. It is shown that the dynamic characteristics of the atmosphere have a significant effect on the formation of macroand microstructural characteristics of convective clouds: the more complex the structure of the wind field in the atmosphere (i.e., the more intense the interaction of the atmosphere and the cloud), the less powerful the clouds are formed.


Introduction
The present period of time is a transition period for cloud physics and cloud seeding: a gradual transition from the stage of studying the "elementary" processes in the clouds to the stage of studying the formation and evolution of clouds as a whole, taking into account their system properties [1] [2].Obviously, such a transition is natural, because cloud physics is not limited to the study of individual processes in the clouds, and there are still many factors affecting the formation and evolution of clouds and requiring study.It should also be noted that the transition of cloud physics to the next stage of evolution will take a certain period of time.This is due to the fact that, on the one hand, such "elementary" processes, whose role in the formation and evolution of clouds are great, remained unexplored or insufficiently studied, and, on the other hand, it is necessary to formulate research directions at the next stage of the evolution of cloud physics, the tasks of these directions, develop methodologies and methods for solving these problems.Taking into account that convective clouds belong to complex physical systems, one of the directions of research of a new stage of its evolution is the determination of the main structure-forming factors for clouds and the study of their influence on the formation of their macro-and microstructural characteristics.Such factors, as noted in [1] [2], are the interaction of clouds with their surrounding atmosphere and the interaction of processes in the clouds (the properties of hierarchy and emergence of systems) [3] [4].
In this paper, we present some results of studies on the effect of the structure of the wind field in the atmosphere on the formation of macro-and microstructural characteristics of convective clouds.It is one of the most important mechanisms of interaction of clouds with their surrounding atmosphere.
Note that this problem has long been in the field of view of researchers.As an example, we can mention papers [5]- [11], which are devoted to the investigation of the possible influence on the processes of cloud formation of various mechanisms of interaction of the atmosphere with clouds.The list of these works can be continued, but without dwelling on this, we note that the authors of these works did not distinguish the belonging of this problem to a new stage in the evolution of cloud physics.In addition, the possibilities of the methods that they used were significantly limited.Therefore, these studies should be carried out from a new perspective and using more effective methods and methodologies.
In this regard, it is important to note that the study of the influence of this factor on the processes of cloud-and precipitation formation processes is possible only on the basis of numerical modeling using full three-dimensional cloud models.In this paper we use for this purpose a three-dimensional numerical model of mixed convective clouds with a detailed account of the processes [1] [2] [12] [13].The model is effectively used to study various issues related to the Atmospheric and Climate Sciences processes of cloud formation.In particular, on the basis of this model, the formation of the thermodynamic, microstructural and electrical characteristics of convective clouds was studied [14] [15] [16], and the formation of their microstructural characteristics was studied in [17] [18].Of the works devoted to the study of the role of the interaction of processes in the clouds on the formation of their microstructural characteristics, it can be noted [19].It presents the results of modeling the influence of deformations in the clouds (formation of the interaction of processes in the clouds) on the formation of microstructural characteristics of hail clouds.We also note that a new method based on the use of the "digital atmosphere" is used to form the input data of the model (initial conditions) [20].

Research Method
The present work uses a three-dimensional model with detailed account of processes, including electrical ones, to study the effect of the wind field structure in the atmosphere on the formation of macro-and microstructural characteristics of convective clouds.
We present a system of equations describing processes in the clouds to get an idea of the model, and for a more detailed acquaintance with the model, we can recommend the papers [4] [5] [6].
Border conditions: The system of equations is applied for the space-time domain: The following notation is used: ( ) w r -velocity vector components; l-inertial force parameter; ( ) ( ) , V m V m -steady rate of fall of liquid and solid particles; -changes in the distribution function of droplets due to microphysical condensation processes, coagulation of droplets, accretion of droplets and crystals, crushing and freezing, respectively; -changes in the distribution function of crystals due to sublimation, accretion and freezing of droplets; The boundaries of the spatial domain are denoted by 0, , 0, x y L L и 0, z L The remaining notations are given in the works [4] [5] [6].
Powerful convective clouds observed in the North Caucasus on 02.09.2010 and 07.06.2012 and accompanied by the fall of hail were selected for the study.Calculations were made of the formation and evolution of clouds for various wind field structures in the atmosphere, which included the actual structure constructed from the aerological sounding data (the Mineralnye Vody Airport, North Caucasus, Russia) of the atmosphere on the days indicated, and the model wind structures constructed using the same data but modified as planned authors.The remaining parameters of the atmosphere remained unchanged during the calculations.In Figure 1, the real (left) and model (right) wind field structures in the atmosphere are given as an example (02.06.2010).The difference between the real and model field structures, as you can see, is that in the second case, at all levels, the wind velocity vector in the atmosphere is directed along the axis OX.
Comparative analysis of the characteristics of the model clouds (i.e. based on the calculation results) to study the effect of the wind field structure on the processes of cloud-and precipitation formation processes at the same time instants was carried out: the characteristics of the velocity fields of vertical air movements in the cloud and near cloudiness, radar reflectivity, as well as the maximum values of cloud parameters.

Results of Calculations
Let us dwell on the results of calculations of the cloud corresponding to the real structure of the wind field in the atmosphere.In the Figure 2     At the same time, the shape of the zone of upflows in the cloud and its inclination to the OX axis change more slowly with time.As for the vertical movements of air around the cloud, the figure shows that the structure of the velocity field of vertical air currents has undergone significant changes.Firstly, in the course of the evolution of the cloud, the area covered by the ascending and descending movements of the air expanded significantly, and secondly, the velocity field structure of these movements became noticeably more complex.As for the maximum values of the velocity of vertical airflows in the cloud, they remain relatively small-17.9m/s for 20th min.and 12.8 m/s for 30th min of the cloud evolution, i.e. there is a slight decrease in the value of this parameter.Figure 4 and Figure 5 show the isolines of the rate of ascending and descending air movements in and around the cloud at the same time in the same vertical plane passing through the center of the cloud along the axis OX.But calculations were made for the model structure of the wind field in the atmosphere, i.e. for the case when the direction of the wind velocity vector at all levels coincides with the direction of the axis OX.In dark colors, the isosurfaces corresponding to the velocities of the ascending currents W = 35.0 and 34.1 m/s are identified in the figures.Comparison of Figure 4 and Figure 5 shows that the structures of the vertical air movement fields at the points in time under consideration differ markedly, and there is also a noticeable increase in the rate of descending air movements.As for the structure of the field of vertical air movements, one can  notice a significant change in it: the descending air currents in the considered plane under the influence of wind in the atmosphere are localized for 30th min in the upper part of the cloud and on its leeward side.
Comparison of Figure 2 and Figure 3 with Figure 4 and Figure 5 shows that the velocity fields of vertical air movements around the cloud corresponding to different structures of the wind field in the atmosphere differ qualitatively: the structures of vertical air movements in the cloud and near-cloud space corresponding to the real structure of the wind in the atmosphere, are noticeably more complex.As for the maximum rate of ascending airflow in the cloud, attention is drawn to the fact that the values of this parameter at the points of time corresponding to the model structure of the wind field in the atmosphere are significantly higher than its values corresponding to the actual structure of the wind field in the atmosphere.
As the results of calculations show, at all stages of the evolution of the cloud, the influence of the wind field structure in the atmosphere on the zone of ascending currents in the cloud is insignificant, which, in our opinion, is associated with large values of the air velocity in this zone.This zone behaves as an obstacle, which is flowed by horizontal movements of air.Based on the results of calculations, it can be noted that the complication of the structure of the wind field in the atmosphere, in particular, the change in wind speed and direction with altitude is a factor hampering the evolution of convection.This may be due to the fact that in this case the interaction of the cloud and its surrounding atmosphere becomes more intense: the exchange of energy and mass between them takes a more intense character.And the flow from the surrounding atmosphere into a cloud of colder air with less water vapor will hamper its evolution.The values and fields of other parameters of the cloud, in particular, water content and ice content, depend in a significant way on the structure of the wind field in the atmosphere.In this case, the water content in the cloud is localized in the region of ascending currents, which agrees with the theory of formation of the liquid-drop fraction in the cloud, and the level of localization of the increased water content may change with time depending on the rate of ascending currents.The region of localization of glaciation is located at higher levels in comparison with the region of localization of droplets.It is important to note that these areas may overlap depending on the stage of the evolution of the cloud.As an example, Figure 6    The obtained results were supplemented with the results of calculations of the maximum values of the cloud parameters at different instants corresponding to different wind field structures in the atmosphere.These parameter values also give useful information about the effect of the wind field structure in the atmosphere on the formation and evolution of clouds.Table 1 gives the maximum values of the main cloud parameters obtained as a result of calculations.The following designations are used in the table: W max , Q max and Z max are the maximum values of the rate of air upflows and water content in the cloud, and also the reflectivity, HW max , HQ max and HZ max are the levels at which these characteristics are located.It can be seen from the table that the maximum values of the velocity of air up flows W max at all stages of the cloud evolution are much larger in the case of the model structure of the wind field in the atmosphere, i.e. in the case when the interaction of the cloud with the surrounding atmosphere is less intense.Accordingly, the heights in which these characteristics are located in the cloud are noticeably higher in this case than in the case of the actual structure of the wind field in the atmosphere.As it is known, the importance of water content plays an important role in the formation and evolution of clouds.It can be seen from the table that the maximum values of this parameter at all times are much larger and located higher in the cloud evaluating in the case of the model structure of the wind field in the atmosphere.In the same way, the maximum glacier values behave in relation to the structure of the wind field in the atmosphere.
Taking into account the maximum water content and altitudes on which they are located in the cloud, it can be noted that in the case when the direction of the wind does not vary in height, the conditions for crystal growth and the formation of hailstones are more favorable.This can be evidenced by the values of maximum reflectivity at different times, corresponding to different structures of the wind field in the atmosphere.In the first case, the maximum value of this pressure; θ -average potential temperature; R-gas constant; ( ) s r -specific air humidity; ( ) S Q r -the total ratio of the mixture of liquid and solid phases in the cloud; ( ) z σ -a parameter that takes into account the change in air density with altitude; ( ) P z и ( ) T r -respectively, pressure and temperature; p c -heat capacity of air at constant pressure; , , K C З L L L -respectively, the specific heat of condensation, subli- mation and freezing;( ) ( ) ( ) of dimensionless pressure, potential temperature and specific humidity from their background values in the Atmospheric and Climate Sciences ambient atmosphere the diffusion of steam into droplets and crystals; dropping water that freezes per unit time in a unit volume of air; formation of fragments during spontaneous freezing of supercooled cloud droplets and their accretion with crystals; 1 2 and I I -sources of droplets and crystals; АВ I -source of artificial crystals under cloud seeding; 0 ε -dielectric constant of va- cuum.
and Figure3

Figure 1 .
Figure 1.Real (left) and model (right) structure of the wind field in the atmosphere, used for calculations.

Figure 3 .
Figure 3. Isolines of different values of the upflow velocity in the vertical plane (30th min).

Figure 4 .
Figure 4. Isolines of different values of the upflow velocity and isosurface W = 35.0m /s in the vertical plane (20th min).

Figure 5 .
Figure 5. Isolines of vertical stream velocity and isosurface W = 34.1 m/s in the vertical plane (30th min).
and Figure 7 show isolines of water content at the 20th min of the cloud evolution, corresponding to the real and model structures of the wind field in the atmosphere.The figures show that the water-localization zones corresponding to different wind structures in the atmosphere differ qualitatively.Significantly, the values (including the maximum) of this parameter are also different.Comparison of the figures shows that in the case of the model structure of the wind field in the atmosphere, the size of the region of localization of water content in the cloud, as well as its values, is much larger than in the cloud evaluating in the atmosphere with real wind field structure.In Figure 7, it can be seen that the isolines of water content reach the surface of the earth, i.e. in the case of the model structure of the wind field in the atmosphere, precipitation is observed in the cloud.This is explained by the fact that in connection with the more intense airflow in the cloud in the case of the model structure of the wind field in the atmosphere, the formation of precipitation particles also becomes more intense.W = 35.0m/s against the background of the isolines of water content in the vertical plane (20th min).

Figure 6 .
Figure 6.Isolines of water content in the vertical plane (20th min).

Figure 7 .
Figure 7.The isosurface of the vertical component of the velocity.W = 35.0m/s against the background of the isolines of water content in the vertical plane (20th min).