1. Introduction
There have been two approaches in flow dynamics: one is the deterministic Newtonian approach and the other is probabilistic (variational), stated in this article as variational field Lagrangian formalism, or in short, variational formalism. The governing equations, conventionally used in meteorology and oceanography, are deterministic, which can be derived from the variational formalism (sections 2, 3, and 4 and the appendix).
The fundamental methodology of the classical variational formalism, the Euler–Lagrange (E–L) procedure with Hamilton’s least-action principle coincides with the methodology that has been broadly used, not only in classical mechanics (Lamb 1932; Bateman 1932; Oden and Reddy 1976; Sasaki 1955), but also in modern advanced physics, with Lagrangian density designed in a physically appropriate form for each problem (Feynman and Hibbs 1965; Kaku 1993; Weinberg 1995; Nair 2005; Mandl and Shaw 2010). The research cited here represents a great number of available references with respect to field Lagrangian formalism used to describe essential mechanisms in both classical and modern physics.
Feynman and Hibbs (1965) and section 5 of this manuscript clearly show the classical limit of the quantum theory to describe the classical systems. The order of magnitudes of the dimensions, mass, times, and so on of the classical system is so large that the action L (of the notation used in this article) is enormous, which is estimated 0.5 × 10−2 mksa or ρ = 10−3 g−1 and |v| = 10 m s−1. Since the Plank’s constant ħ = h/(2π) = 1.054 × 10−34 mksa rationalized system (J, s), the magnitude of L/ħ is on the order of 1032, so that the quantum field formalism has not been considered normally in the classical flow formalism. We will reconsider it, showing a possible need for both systems, especially their interactions between classical and quantum systems.
The entropic balance theory was originally developed on the same basis as the variational field Lagrangian formalism to derive the same equation of Clebsch transformation (Lamb 1932). The Clebsch transformation equation was developed in the variational formalism but has not been used directly for actual phenomena because of the unobservable, nonmeteorological, but mathematical Lagrange multiplier. The Lagrange multipliers for the constraints of mass and entropy that appeared in the Clebsch transformation are analogous to the mathematical vector potential (and gauge field) of the theoretical process found known as the “Aharonov–Bohm effect” (or “A-B effect” or “A effect”) (Feynman et al. 1964; Feynman and Hibbs 1965; Landsman 1998; section 5 herein). Its experimental verification has been difficult and has not been made until two decades later (Tonomura 1997).
The entropic balance theory, changing the unobservable nonmeteorological Lagrange multiplier to observable meteorological rotational flow velocity with entropy, is now applicable to tornadic phenomena and is made applicable to solve mysteries of tornadogenesis. Furthermore, the balance is a newly found one, because it is different from the other known balance conditions, such as hydrostatic, (quasi-) geostrophic, cyclostrophic, Boussinesq, and anelastic balance.
The variational formalism is similar to the one observed by Euler in 1736 and referred to as the Euler equation (Oden and Reddy 1976; Lanczos 1970) and later as the Gateaux derivative (Gateaux 1913) in linear approximation with respect to the hypothetical variation that has two parameters, one infinitesimal and the other arbitrary, and the Frechet differential in normed space that deals with a variety of broader problems. A similar mathematical manipulation leads to a full set of dynamical and thermodynamical nonlinear equations of the ideal flow (Lamb 1932; Bateman 1932; Sasaki 1955). The E–L equations are all prognostic except for one that is diagnostic, so-called by Lamb as the Clebsch’s transformation (Clebsch 1857) of flow velocity, although the terminology was not commonly used by authors who derived the same or a similar equation.
In the action Lagrangian of the variational formalism by Bateman (1932) and Lamb (1932), the thermodynamic property is assumed to be barotropic, that is, the flow pressure is a function of density only, in order to derive the equation of motion. Baroclinicity is included by Sasaki (1955, 1999, 2009, 2010) and Dutton (1976). Salmon (1988, 1998)’s Hamiltonian formalism leads to the Clebsch transformation equation. Note that Hamiltonian formalism with a Poisson bracket is a powerful tool for quantization. Dutton compared two types of variations, Lagrangian (material) variation and Eulerian variation, and obtained similar results to those of variational formalism. Sasaki further modified the Clebsch transformation and developed the entropic balance equation that is made applicable to real meteorological phenomena (section 4). The entropic balance theory was developed by introducing the two hypotheses in the Lagrangian density: hypothesis 1 of discontinuous sudden cloud physical phase changes, compared with longer time scales of a tornado, a supercell, and daily synoptic weather, and hypothesis 2, concerned with the ensemble property of such weather systems. Of course, it also includes the continuous change of entropy flux through the boundary of the domain.
In sections 2, 3, and 4 and also in the appendix, we discuss Lagrangian density and the action Lagrangian physically appropriate for the tornado mechanism and derive the Clebsch transformation and entropic balance equation.
2. Action Lagrangian of tornadic storm

Note that the notation S instead of L is commonly used in the variational formalism, but S is used for entropy and L is used for the action Lagrangian in this study.
3. Variational formalism
In the application of (3.4) and (3.7) to the Lagrangian [(2.3a)] and the action Lagrangian [(2.3b)], we find that there exists a sole diagnostic E–L equation, that is, the Clebsch transformation equation and the entropic balance equation, among all other prognostic E–L equations. The Clebsch transformation has not been applied to real meteorological data in the past. The entropic balance equation is developed from the Clebsch transformation and is found to explain several mysteries of the tornadogenesis mechanism. We will discuss both the Clebsch transformation and the entropic balance equation and how the entropic balance equation plays roles to explain the tornadogenesis mechanism in the next section.
4. Diagnostic E–L equation: Clebsch transformation, entropic balance equation, and tornadogenesis
Solution space. The domain of full solution of the E–L equations is schematically shown in the solution space by the thick solid line. It includes nonstationary state (NSS), stationary state (SS), and the solution of the diagnostic E–L equation (DS). The solution in the domain covered by DS and SS has long-lasting property mathematically similar to the attractor.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
Equations (4.5)–(4.7) are the basic equations of the entropic balance theory, which reveal various important properties of supercells and tornadoes and furthermore the wraparound mechanism of tornadogenesis. In practice, S is estimated from radar reflectivity analysis, and vR is obtained from the variational data assimilation of Doppler velocity radar observation, as planned to be reported separately, but estimated approximately from Doppler velocity observation. The approximate estimate was determined to have worked well to estimate vorticity in testing. Note that the wind under entropic balance is much more general (there is less assumption) than the thermal wind balance that is under the geostrophic or quasigeostrophic balance.
Figure 2 illustrates schematically the entropic right-hand rule. The diagnostic velocity equation [(4.5)] is universal for the ideal flow. The vorticity equation [(4.6)] derived from (4.5) is conveniently demonstrated by the mutually orthogonal vector relation, similar to the so-called Fleming’s right-hand law of electromagnetic fields (called the “entropic right-hand rule” by the author), among the orthogonal variables of the spatial three dimensions: the vorticity ω, the entropy gradient (1/S)∇S, the rotational flow velocity component −S∇β, denoted by vβ or vR, and the divergent component −∇α, denoted vα or vD. These notations are used in the figure illustrations in this article. Based on the entropic balance theory, the wraparound mechanism is introduced by Sasaki (1999, 2009, 2010) to explain explicitly the nonlinear process of tornadogenesis. The results are consistent with tornadic storm observations and successful tornado simulations of phenomena in tornadic storms, such as overshooting hydrometeors against the upper-level westerlies. The development of a mesocyclone, a hook echo, and a wall cloud is easily explained by the entropic balance theory. The almost discontinuous transition from supercell to tornado suggested from visual observation and data analysis is also easily explained by the entropic balance theory as a transition from baroclinic to barotropic stages. Furthermore, an increase of the relative helicity to 1 (its maximum value) as a result of computer simulation is explained by the entropic balance theory. From the entropic balance theory, we can easily explain the observation of the tornado touching the ground in the perpendicular direction.
Entropic right-hand rule. The rule shows the mutually orthogonal vector relation, similar to the so-called Fleming’s right-hand rule of electromagnetic fields, which we call the entropic right-hand rule, among the orthogonal variables of three spatial dimensions: entropy gradient (1/S)∇S, rotational component of flow velocity vR or vβ := −S∇β, and vorticity ω.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
For simplicity, the vorticity in balance with the entropic gradient and rotational component of flow velocity as given by (4.6) derived from the entropic balance theory is schematically illustrated as the entropic right-hand rule (Fig. 2). Note that when the gradient of entropy is of larger magnitude, it is baroclinic (Fig. 3), while when it is of smaller magnitude, it is barotropic (Fig. 4). A schematic diagram of a supercell with a tornado is shown in Fig. 5, in which various known features of the supercell and tornado, together with a newly found mechanism, such as the overshooting of hydrometeors against the upper headwind westerlies, are explained by the entropic balance theory in the caption of the figure. Figure 6 shows a horizontal cross section of entropic source and sink that produces a baroclinic gradient in balance with the mesocyclone vortex at middle levels according to the entropic balance theory. The existence of a single vortex (Fig. 7) and that of multiple vortices (Fig. 8, for four vortices) are easily understood to exist by applying the entropic right-hand rule (Fig. 2) to each vortex.
Baroclinic case of the entropic balance expressed as the right-hand rule.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
Barotropic case of the entropic balance expressed as the right-hand rule.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
Schematic diagram of tornadogenesis based on the entropic balance theory. Meandering westerlies transport water vapor evaporated from Gulf Stream deeper inland to the central U.S. Great Plains area and meet with the dry air to onset a supercell. The moisture of the southerly flow condenses and releases the latent heat to the surrounding air resulting in an entropic source that is shown by the circled plus sign, called the entropic source here. The hydrometeors, such as raindrops and ice particles, created by the condensation in the entropic source area are lifted by the updraft of the thermal convection of the storm reaching near the cloud top and are blown away toward the downstream (eastern) side of the storm. However, some of the lifted hydrometeors are overshot toward the upstream direction against the upper-air westerlies. This is due to the upper-air horizontal vortex, where the rotational flow direction is shown by the arrow with double solid lines. The overshot hydrometeors will fall down, evaporating because of the dry surrounding air, thus cooling the air. The descending hydrometeors with cooled air meet with a dry midlevel southwesterly jet and are cooled further to produce the rear-frank downdraft. Thus, an entropic sink forms nearly at the same level as the entropic source. The horizontal spatial gradient entropy is generated by the source–sink pair. A vortex (mesocyclone) is formed, and the wraparound mechanism is organized. The wraparound mechanism becomes activated by the mesocyclone that existed between the entropic source and sink. The wraparound mechanism is a nonlinear process (see Fig. 9), similar to the attractor, that generates a hook echo, a low-level mesocyclone, a wall cloud, and a tornado. These processes are explained by the diagnostic E–L equation, the entropic right-hand rule, and the wraparound mechanism that are derived by the entropic balance theory.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
An entropic vortex exists in balance with the gradient of entropy in the baroclinic field. The horizontal vortex is formed due to the slanted vertical entropy gradient at the upper levels above the entropic source and overshoots the hydrometeors upstream against the headwind westerlies. The vertical vortex is formed at the middle levels due to the horizontal gradient between the source and the sink of entropy. The entropic vortex formation is explained by the entropic right-hand rule derived from the entropic balance theory.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
Schematic illustration of single-vortex tornadogenesis explained by the wraparound mechanism. The wraparound mechanism is modeled by the entropic right-hand rule that is derived from the entropic balance theory.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
Schematic illustration of multiple-vortex tornadogenesis, for four vortices as an example, explained by the wraparound mechanism.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
Transition from the mesocyclone, in baroclinic state, shown in Figs. 5 and 6 to the single vortex and/or multiple vortices, in barotropic state, as shown in Figs. 7 and 8, occurs through the wraparound process. The newly found wraparound mechanism for tornadogenesis is analogous to a nonlinear process, the so-called Baker’s transformation, and the transition is discontinuous from baroclinic to barotropic stages by trapping the entropic sink core inside the vortex, like a nonlinear attractor (Figs. 1 and 9). Note also that the wraparound mechanism is two-dimensional, while Baker’s transformation is one-dimensional.
Wraparound mechanism (two dimensional) analogous to the nonlinear Baker’s transformation (one dimensional). The wraparound mechanism is modeled by the entropic right-hand rule that is derived from the entropic balance theory.
Citation: Journal of the Atmospheric Sciences 71, 6; 10.1175/JAS-D-13-0211.1
The axisymmetric single vortex is theoretically well supported by the Noether’s rotational invariant theorem (Kaku 1993) with the Lagrangian density [(2.3a)] and the application of the variational formalism. The existence of multiple vortices is also proven by the Noether’s rotational invariant theorem extending it to application for sectional rotation.
Many visual observations of cases of single-vortex and multiple-vortex tornadoes are recorded. For instance, an excellent example is shown by the case of the tornado in central Oklahoma area on 3 May 1999. Multiple-vortex funnels were observed near Chickasha, Oklahoma, and a single EF5 tornado was observed near Moore, Oklahoma, likely from the same supercell during a few hours.
5. Conservation constraint of the classical system: Encompassing classical flow and quantum field mechanics

Accordingly, the inclusion of the quantum flow formalism based on Feynman and Hibbs (1965) may serve as a comprehensive way for the variational field Lagrangian approach, a new area of research. The variations in calculus of variations and the classical variational formalism are hypothetical. However, this study suggests that the variations are physical, relating to quantum variations and their interaction with classical systems. It would cast some light on the still mysterious lightning and microphysical processes of storm clouds (Petersen et al. 2008). It would justify the use of the variational formalism better than a deterministic Newtonian approach for the classical system.
Acknowledgments
The author expresses his appreciation for the encouragement and support given to this research by Dr. Kelvin Droegemeier, Vice President for Research, Dr. Robert Palmer, Director, Advanced Radar Research Center, University of Oklahoma, and Dr. Seon Ki Park of Ewha Women’s University, and also to Weathernews, Inc., for generous support. Also, the author expresses his appreciation to faculty and scientists especially of NWC, University of Oklahoma, University of Tokyo, Kyoto University, Japan, for their useful discussion and to Vivek Mahale and Mark Laufersweiler of the University of Oklahoma for his assistance with computer graphics. I would like to acknowledge that Springer publisher gave permission to the author for the use of nine figures for this manuscript from the published book Data Assimilation for Atmospheric, Oceanic, and Hydrostatic Applications (Vol. 2), Chapter 18: Entropic Balance Theory and Radar Observation for Prospective Tornado Data Assimilation, Yoshi K. Sasaki, Mathew R. Kumjian, and Bradley M. Isom.
APPENDIX
Prognostic E–L Equations
REFERENCES
Bateman, H., 1932: Partial Differential Equations of Mathematical Physics. Cambridge University Press, 522 pp.
Clebsch, A., 1857: Über eine allgemeine Transformation der hydrodynamischen Gleichungen. Z. Reine Angew. Math., 54, 293–312.
Dutton, J. A., 1976: The Ceaseless Wind: An Introduction to the Theory of Atmospheric Motion.McGraw-Hill, 579 pp.
Feynman, R. P., and A. R. Hibbs, 1965: Quantum Mechanics and Path Integrals.McGraw-Hill, 365 pp.
Feynman, R. P., R. B. Leighton, and M. Sands, 1964: The Feynman Lectures on Physics. Vol. 2. Addison-Wesley Publishing, 592 pp.
Gateaux, R., 1913: Sur les fonctionnelles continues et les fonctionnelles analytiques. C. R. Acad. Sci., 157, 325–327.
Kaku, M., 1993: Quantum Field Theory: A Modern Introduction.Oxford University Press, 785 pp.
Kashiwa, T., 2006: New Edition Exercise Quantum Field Theory (in Japanese). Publishing Company Science, 223 pp.
Lamb, H., 1932: Hydrodynamics. 2nd ed. Cambridge University Press, 738 pp.
Lanczos, C., 1970: Variational Principles of Mechanics.4th ed. University of Toronto Press, 418 pp.
Landsman, N. P., 1998: Mathematical Topics between Classical and Quantum Mechanics.Springer-Verlag, 529 pp.
Mandl, F., and G. Shaw, 2010: Quantum Field Theory.2nd ed. John Wiley and Sons, 478 pp.
Nair, V. P., 2005: Quantum Field Theory: A Modern Perspective. Springer, 557 pp.
Oden, J. T., and J. N. Reddy, 1976: Variational Methods in Theoretical Mechanics.Springer-Verlag, 302 pp.
Petersen, D. M., M. Bailey, W. H. Beasley, and J. Hallett, 2008: A brief review of the problem of lightning initiation and a hypothesis of initial lightning leader formation. J. Geophys. Res., 113, D17205, doi:10.1029/2007JD009036.
Salmon, R., 1988: Hamiltonian fluid mechanics. Annu. Rev. Fluid Mech., 20, 225–256, doi:10.1146/annurev.fl.20.010188.001301.
Salmon, R., 1998: Lectures on Geophysical Fluid Dynamics.Oxford University Press, 378 pp.
Sasaki, Y. K., 1955: A fundamental study of numerical prediction based on the variational principle. J. Meteor. Soc. Japan, 33, 262–275.
Sasaki, Y. K., 1999: Tornado and hurricane: Needs of accurate prediction and effective dissemination of information (in Japanese). J. Visualization Soc. Japan,19, 187–192.
Sasaki, Y. K., 2009: Real challenge of data assimilation for tornadogenesis. Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications, S.-K. Park and L. Xu, Eds., Springer-Verlag, 97–126.
Sasaki, Y. K., 2010: Entropic Balance Theory and Tornadogenesis.NOVA Science Publishers, 39 pp.
Sasaki, Y. K., M. R. Kumjian, and B. M. Isom, 2013: Entropic balance theory and radar observation for prospective tornado data assimilation. Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications, S.-K. Park and L. Xu, Eds., Springer-Verlag, 429–464.
Tonomura, A., 1997: To See Gauge Field: To Reveal the Micron World by Electronic Wave (in Japanese). Koudansha, 260 pp.
Weinberg, S., 1995: Foundations. Vol. 1, The Quantum Theory of Fields, Cambridge University Press, 609 pp.
Yoshida, T., 2000: Maxwell: Dancing of Field and Particle (in Japanese). Kyoritsu Publishing, 136 pp.