# Weather and Climate Manipulation as an Optimal Control for Adaptive Dynamical Systems

Authors: Sergei A. Soldatenko
Journal: Complexity
 (2017)

---

## Abstract

The weather and climate manipulation is examined as an optimal control problem for the earth climate system, which is considered as a complex adaptive dynamical system. Weather and climate manipulations are actually amorphous operations. Since their objectives are usually formulated vaguely, the expected results are fairly unpredictable and uncertain. However, weather and climate modification is a purposeful process and, therefore, we can formulate operations to manipulate weather and climate as the optimization problem within the framework of the optimal control theory. The complexity of the earth’s climate system is discussed and illustrated using the simplified low-order coupled chaotic dynamical system. The necessary conditions of optimality are derived for the large-scale atmospheric dynamics. This confirms that even a relatively simplified control problem for the atmospheric dynamics requires significant efforts to obtain the solution.


---

## Body

## 1. Introduction

Weather and climate manipulation is defined as a deliberate intervention in the natural evolution of physical processes in the atmosphere and other components of the earth’s climate system (ECS) in order to achieve some desired results [1, 2]. Here, as usual, the ECS (S) is understood as a complex large-scale physical system that consists of five basic and interacting constituent complex subsystems, namely [3],(i)
atmosphere (A), the gaseous and aerosol envelope of the earth that propagates from the land, water bodies, and ice-covered surface outward to space,
(ii)
hydrosphere (H), the oceans and other water bodies on the surface of our planet and water that is underground and in the atmosphere,
(iii)
cryosphere (C), the sea ice, freshwater ice, snow cover, glaciers, ice caps and ice sheets, and permafrost,
(iv)
lithosphere (L), the solid, external part of our planet,
(v)
biosphere (B), the part of our planet where life exists; that is,(1)
S
=
A
∪
H
∪
C
∪
L
∪
B
.Each component of the ECS is characterized by a finite set of variables, usually called state variables, whose values at a given time determine the instant state of the ECS. All components of the ECS are extremely complex physical systems that in turn also are composed of constituent elements (subsystems). It should be noted that the current status of climate system studies is characterized by dramatic progress achieved over the last few decades in understanding basic physical processes and fundamental feedback mechanisms that control the climate formation and change (thanks to the advances in computing hardware and algorithms and due to a substantial increase in the volume of climatological data [4–8]). The detailed retrospective analysis of climate science can be found, for example, in [9]. However, our knowledge for climate remains incomplete due to the enormous complexity of the ECS. A deeper understanding of current climate, its variability, and change is among the greatest challenges facing modern science that can help not only to project the future climate but also to develop methods for controlling the ECS.The objectives and methods of weather modification are distinctly different than those of climate manipulation. This distinction is due to the fact that the terms “weather” and “climate” have different meanings. Weather is defined as the daily conditions of the atmosphere in terms of such atmospheric variables as temperature, humidity, wind direction and velocity, surface pressure, clouds, and precipitation. In contrast, climate represents the ensemble of states traversed by the climate system over a sufficiently long enough temporal interval. Here, the ensemble means not only the set of system states, but also a probability measure defined on this set. Climate, roughly speaking, can be considered as the “average” weather, in terms of mean and variance, in a certain geographical location over many years.Weather modification is aimed at altering weather conditions and weather phenomena locally or regionally. The successful examples of weather modification operations are cloud seeding to increase precipitations, the removal of fogs, stratiform and convective clouds, the suppression of hail, and lightning (e.g., [10–14]). The large-scale weather modification experiments such as the reduction of the intensity of hurricanes were less successful [15–17]. Thus weather modification aims to intentionally affect processes and phenomena only in the atmosphere, the most volatile component of the ECS. Climate manipulation known also as a climate engineering or geoengineering was suggested as a response measure on the global climate change that currently occurs on our planet [3]. Geoengineering aiming at stabilizing the global climate and preventing further warming is a planetary-scale intentional intervention on the ECS and represents a multidisciplinary problem of extreme complexity, which is currently explored mostly theoretically because the consequences and unintended side effects of climate manipulation operations are very uncertain. Geoengineering requires further consideration of not only scientific and technical issues, but also the legal and ethical aspects and limitations. It is expected that the interest in climate and weather manipulation will continue to grow, so that it makes it necessary to develop the appropriate theoretical framework. To date a number of methods and technologies were proposed by scientists and engineers to stabilize the ECS and to mitigate the global warming (e.g., [18–24]).The realization of both weather modification and climate engineering operations is a goal-directed process, having specific objectives that should be somehow formulated and certainly achieved. Mathematically, weather and climate manipulation can be viewed as an optimization problem [25, 26]. Consequently, instead of examining weather and climate manipulation as amorphous actions without clearly defined objectives, boundaries, and methods, we can consider operations to modify weather and climate as an optimal control problem for the ECS and its components.It is important to make the following comments concerning controllability of the ECS and its components and, in particular, the atmosphere. Controllability can be thought of as a property of a control system that indicates the ability to transfer a system from one state to another, using admissible operations (e.g., [27]). From the perspective of physics, the ECS is a controllable system, and its key driver is the incoming shortwave solar radiation. On the geological time scale, the earth’s climate has been highly variable throughout history as a result of both external forcing and internal perturbations [28] (e.g., variations in the chemical and aerosol composition of the atmosphere and aperiodic El Niño-Southern Oscillation). Since the incoming solar radiation can be viewed as a primary control variable, geoengineering seeks to reduce global warming by using physically feasible ability to manipulate the climate via changes in the amount of incoming solar radiation (by varying reflection or absorption) as well as natural feedback mechanisms that exist in the ECS. However, geoengineering is not the same as control, and thus the term “engineering” rather than “control” is currently used in the literature to emphasize this difference. Considering geoengineering within the formal framework of the optimal control theory, we are faced with a number of very hard mathematical problems, such as the proof of the controllability property for any particular climate model.Similar problems arise when we consider weather modification. For example, we can create the artificial precipitation by cloud seeding; however the formal proof of the existence of the optimal solution and the proof of the controllability property for an applied mathematical model are problems of extreme difficulty. For this reason, in weather modification, the term “intentional (deliberate) action” is usually used instead of “control.” Thus, we have to make the distinction between control as an engineering problem and control as a mathematical problem.In this paper, the weather and climate manipulation is examined as an optimal control problem for the climate system, which is considered as a complex adaptive dynamical system. The complexity of the climate system is discussed and illustrated using the simplified low-order coupled chaotic dynamical system. Finally, as an example, the necessary conditions of optimality are derived for the large-scale atmospheric dynamics. This shows that even a relatively simplified control problem for the atmospheric dynamics requires significant efforts to obtain the solution.


## 2. Complexity in the Earth’s Climate System

Let us make some preliminary comments. The term “system” generally refers to a goal-oriented set of interconnected and interdependent elements that operate together to achieve some objectives [29]. The system is called complex if it possesses such characteristics as emergent behavior, nonlinearity and high sensitivity to initial conditions and/or to perturbations, self-organization, chaotic behavior, feedback loops, spontaneous order, robustness, hierarchical structure, and some others [30]. Complexity in systems appears in particular from spatiotemporal interactions between their components. These nonlinear interactions lead to the appearance of new dynamical properties (e.g., synchronous oscillations and other structural changes) that cannot be observed by exploring constituent elements individually. Let us go one step forward. Complex systems include a special class of systems that have the capacity to adopt to system’s environment. These systems are known as complex adaptive systems (e.g., [31, 32]). In complex adaptive system its parts are linked together in such a way that the entire system as a whole has the capacity to transform fundamentally the interrelations and interdependences between its components, the collective behavior of a system, and also the behavior of individual components due to the external forcing. Complex adaptive systems are dynamical systems since they evolve and change over time. These systems have a number of properties that include the following [32, 33]: coevolution, connectivity, suboptimality, requisite variety, iteration, edge of chaos and, certainly, emergence, and self-organization.The ECS is a large-scale and peculiar dynamical system that possesses numerous unique physical, dynamical, and chemical properties (e.g., [7, 8, 33]). The study and a fortiori control of the ECS are a highly difficult problem. There is no doubt that the ECS can be viewed as a complex adaptive system. This is because of the following.(i)
The ECS is a complex system combining the atmosphere, hydrosphere, cryosphere, land, and biota together with global biochemical cycles (first of all, cycles of CO2, N2O, and CH4) and aerosols. The climate system’s components are heterogeneous thermodynamical systems characterized by specific variables that determine their states. Subsystems of the ECS have strong differences in their structure, dynamics, physics, and chemistry. They embrace processes with different temporal and spatial scales and link together via numerous coupling physical mechanisms, which can be weak and strong. Each constituent component of the ECS is also a complex and possibly also adaptive system. For example, the following weather systems could be viewed as subsystems of an overarching atmospheric system: tropical cyclones, hurricanes, mid-latitude cyclones and anticyclones, clouds, precipitations, and so forth.
(ii)
The ECS has a hierarchical structure. Each component of the ECS can be characterized by a specific response time. The belonging of a certain component to the ECS is determined by the ratio between the temporal scale of processes under consideration and the response time of the ECS components. For example, the atmosphere, which has a response time of about one month in the troposphere, can be considered as a sole component of the ECS for processes with temporal scales of days to weeks. In this case, oceans, land surface, and ice cover are considered as the boundary conditions and/or external forcing. If we study processes which have temporal scales of months to years, the atmosphere and ocean must be included in the ECS together with sea ice.
(iii)
The ECS has a large number of positive and negative feedback mechanisms, for example, ice-albedo feedback (positive feedback), water vapor feedback (positive feedback), cloud feedback (both positive and negative feedbacks), carbon cycle feedback (negative feedback), and feedback due to Arctic methane release (positive feedback).
(iv)
Physical and dynamical processes in the ECS cover a broad spectrum of temporal and spatial scales. Time scales are varied from seconds to decades, and spatial spectrum of dynamical processes covers molecular to planetary scales. Dynamical processes in the ECS and its components are nonlinear. Subsystems of the ECS interact with one another nonlinearly producing under certain conditions a chaotic behavior of subsystems and the overall climate system.
(v)
The ECS and its components inherently have emergent properties. Examples of atmospheric emergent phenomena include but not limited to clouds, large-scale eddies (cyclones and anticyclones), and small-scale vortices such as tornados. Examples of climate emergent phenomena are the El Niño-Southern Oscillation, which is quasiperiodic irregular variation in ocean surface temperature over the Pacific in tropics that strongly influences the global climate and ocean circulation patterns. Natural emergent phenomena appear spontaneously under certain favorable conditions.
(vi)
The ECS is thermodynamically open and nonisolated system because it exchanges energy with its surroundings. However, the ECS is a closed system for the exchange of matter with outer space. The energy that drives the ECS is a solar energy. The ECS is affected by changes in external driving forces, which imply natural causes such as solar activity variations and volcanic activities and human-made changes in chemical composition of the atmosphere. However, the impact of the ECS on the outer space is insignificant. Currently, changes in climate are mostly affected by variations in the atmospheric composition of particles and gases. In the Arctic the role of changes in albedo (reflection coefficient) is also tangible. In spite of the fact that the most significant heat-trapping gas is atmospheric water vapor, the most influential gas component to affect the climate change is CO2, which comprises about 70 percentage of points of the global warming potential [3].
(vii)
The components of the ECS are also nonisolated systems. They act as cascading systems and interact with each other in various ways including the transfer of momentum, sensible and latent heat, gases, and particles. All together they compose the climate system, which is a unique large-scale natural system.
(viii)
Dynamical processes in the ECS periodically and irregularly oscillate due to both internal factors (natural oscillations) and external forcing (forced oscillations). Natural fluctuations are due to internal instability (e.g., convective, barotropic, and baroclinic instabilities) with respect to stochastic perturbations. Human impacts, both intentional and unintentional, are considered as external forcing.Certainly, there are other specific properties of the ECS that should be taken into account while studying climate as a complex adaptive system.With regard to weather and climate manipulation, it is necessary to underline that the response of nonlinear systems to external perturbations is completely different, both quantitatively and qualitatively, than the response of linear systems [34]. Recall that operations to manipulate weather and climate represent external forcing with respect to the ECS. This is very important in terms of assessing the effects of weather and climate manipulation operations.To study and simulate the ECS we should assign some mathematical object that is an abstract representation of the real climate system. This object is usually known as an “ideal” model of the system of interest. The climate system model represents a set of interacting and interdepending subsystems (agents). The number of these subsystems is determined by the objectives of problem under consideration. For example, to study the large-scale climate variability the model can include the following major agents: tropical, mid-latitude and polar troposphere, stratosphere, land ice, oceans and sea ice, surface and boundary layers, hydrological cycle, clouds of various types, precipitations, aerosols, CO2 and CH4 cycles, solar radiation, and terrestrial emission. Other subsystems of the ECS (e.g., vegetation, land surface, and biota) can be considered as the boundary conditions and external forcing.The major components of the ECS are physical continuum and their evolution can mathematically be described by a set of multidimensional nonlinear differential equations in partial derivatives:(2)
∂
t
φ
r
,
t
=
M
φ
r
,
t
,
α
r
,
t
,
φ
0
,
t
=
φ
0
r
,where φ is a state vector of a system, α is a parameter vector, r
∈
D
⊂
R
3 is a vector of spatial variables, t
∈
0
,
T is the time, M is a nonlinear differential operator that describes the evolution of a system, and φ
0 is initial conditions. Since (2) describes a continuous medium, the state vector φ is infinite-dimensional. The solution to such infinite-dimensional system cannot be found analytically, and we need to employ numerical methods. To obtain numerical solution, the original set of (2) is replaced with discrete spatiotemporal approximations using some appropriate technique (e.g., finite-difference method and Galerkin approach). It is very important to underline that a large number of physical processes and cycles (model agents), which also may be viewed as a system, cannot be explicitly represented in the climate model due to its discrete spatiotemporal structure. State-of-the-art models of the ECS are generally unable to realistically simulate processes on spatial scales of the order of twice the model grid length [35]. Such thermodynamical, physical, and chemical processes and cycles are parameterized, that is, expressed parametrically using simplified description. Some of these newly introduced parameters can be considered as controls (see below). By varying control parameters, we can formally manipulate the ECS and its components.Models of the ECS due to their extreme complexity require extensive computational resources to study the behavior of climate system and its components by computer simulations. However, in order to mimic some features of the ECS behavior and to study its certain properties as a complex system we can use simplified low-dimensional models. For illustrative purposes, we consider a coupled nonlinear model, which is composed of fast (the “atmosphere”) and slow (“the ocean”) components and obtained by coupling of two versions of the original Lorenz system [36] with distinct temporal scales (e.g., [37–39]):(3)
d
x
d
t
=
σ
y
-
x
-
c
a
X
+
k
,
d
y
d
t
=
r
x
-
y
-
x
z
+
c
a
Y
+
k
,
d
z
d
t
=
x
y
-
b
z
+
c
z
Z
,
d
X
d
t
=
ε
σ
Y
-
X
-
c
x
+
k
,
d
Y
d
t
=
ε
r
X
-
Y
-
a
X
Z
+
c
y
+
k
,
d
Z
d
t
=
ε
a
X
Y
-
b
Z
-
c
z
z
.Here lower case and capital letters represent, respectively, the “atmosphere” (fast model) and the “ocean” (slow model); σ, r, and b are the parameters of Lorenz system; c and c
z are parameters that describe the coupling strength between fast and slow models, a is the amplitude scale factor, k is an “uncentering” index, and ε is the time-scale separation parameter. For simplicity but without loss of generality we can assume that a
=
1, k
=
0, and c
=
c
z. The combination of two Lorenz systems allows one to imitate the interaction between fast-oscillating atmosphere and slow-fluctuating ocean [40, 41].The Lorenz system is one of the most well-known nonlinear dynamical systems, which derived from the equations of atmospheric Rayleigh-Bénard convection [42]. Three variables of the Lorenz model, x, y, and z, are, respectively, the intensity of convective motion and horizontal and vertical temperature gradients. The model parameters have the following meanings: σ is the Prandtl number, r is a normalized Rayleigh number, and b is a nondimensional wavenumber. Properties of the Lorenz model are well known and described in the literature (e.g., [43]). We only note that the original Lorenz system can be studied using different values of parameters. However, for the standard parameter values σ
=
10, b
=
8
/
3, and r
=
28, the Lorenz model exhibits chaotic behavior and possesses a two-wing strange attractor, which has a fractal dimension of about 2.06 [36, 43]. For σ
=
10 and b
=
8
/
3, the chaotic behavior is detected if r
>
r
c, where r
c
≈
24.74. Varying the parameter r, we can observe structural changes in the system dynamics. In our computations the standard values of parameters were used. We also assume that ε
=
0.1.The dynamics of coupled system (3) strongly depends on the parameter c. This parameter, referred to as the coupling parameter, describes the interaction degree between the fast and slow systems [44]. In the existing climate models, atmosphere and ocean are commonly coupled by adjusted air-sea fluxes of heat, momentum, and fresh water using various parameterization techniques (e.g., [45]). In this context, the parameter c can be viewed as a simple parameterization of the interaction between two models used in this study. Figures 1–4 show the temporal changes of fast (x and z) and slow (X and Z) variables calculated by the numerical integration of the coupled model equations. Even a cursory examination of these figures reveals that when a coupling strength parameter is varied, the qualitative behavior of the system changes. Qualitative shifts in the dynamical properties of coupled system can be traced by the spectrum of conditional Lyapunov exponents. These exponents, on the one hand, characterize the average rate of exponential divergence of nearby orbits in the phase space, and, on the other hand, they are used to analyse the synchronization with coupled systems [44]. When the coupling strength parameter tends to zero, coupled system (3) has six distinct Lyapunov exponents: two positive, two negative, and two zero. The dependence of two largest Lyapunov exponents on the parameter c is shown in Figure 5 [39]. These exponents are monotonically decreasing functions. Being initially positive, they approach the x-axis at about c
≈
0.8 and at about c
≈
0.95 become negative. Thus, the phase-synchronous regime is observed when c
>
0.95. When c
>
1.0, the system approaches a limit cycle since all six Lyapunov exponents become negative.Figure 1
Fast (x and z) and slow (X and Z) variables for the case of c
=
0.15.
(a)
(b)Figure 2
Fast (x and z) and slow (X and Z) variables for the case of c
=
0.5.
(a)
(b)Figure 3
Fast (x and z) and slow (X and Z) variables for the case of c
=
0.8.
(a)
(b)Figure 4
Fast (x and z) and slow (X and Z) variables for the case of c
=
1.0.
(a)
(b)Figure 5
Two largest Lyapunov exponents.However, the real climate system is a forced system. To explore the influence of external forcing (constant, periodic, aperiodic, stochastic, etc.) on the coupled system dynamics one can add additional terms in the right-hand sides of model equation (3) (e.g., [46]). This allows the testing of various strategies to manipulate the system dynamics.


## 3. Control versus Manipulation of the Earth’s Climate System

Operations to manipulate the weather and climate are mainly considered as intentional actions to influence natural processes within the ECS. Objectives of these actions are commonly formulated vaguely, in general terms and, consequently, expected results are fairly unpredictable and uncertain. However, since weather and climate manipulation is a goal-directed process, we find it useful to formulate operations to manipulate weather and climate as the optimization problem within the framework of the optimal control theory (e.g., [25, 26]). Then, instead of mathematically amorphous problem, we can obtain an optimal control problem for the ECS, which is formulated mathematically accurately. In this context the climate system and/or its components can be viewed as a control system and specially trained human operators will perform the role of controller.The fundamental difference between weather and climate manipulation and control can be illustrated by the combination of block diagrams and transfer functions (e.g., [47, 48]). The transfer function of a system is simply the ratio of the Laplace transform of the system output to the Laplace transform of the system input assuming initial conditions are equal to zero. This approach used, for example, in signal processing, telecommunication, and control theory allows deriving the transfer function of complex system by pure algebraic operations. Transfer functions are applied for the study of linear, time invariant systems. Of cause, models used in weather and climate simulations are extremely complex and nonlinear. However, a number of properties of nonlinear systems can be explored using the assumption that for small perturbations the system is approximately linear. Figure 6 shows the block diagram for weather and climate manipulation operations. The signal Z
s fed to the system is a sum of state variables X
s and man-made perturbations W
s, which serve as intentional actions to manipulate the weather and climate. The outcome of these intentional actions strongly depends on spatiotemporal distribution and magnitude of man-made perturbations. The determination of these perturbation parameters is performed heuristically, with the use of expert estimations. Therefore, the results of weather and climate manipulation are very uncertain and unpredictable. The transfer function for weather and climate manipulation operations is expressed as follows: (4)
T
s
≡
Y
s
X
s
=
G
s
1
+
W
s
X
s
.Note that for climate manipulation W
s
/
X
s
≪
1, while for weather modification the ratio between X
s and W
s can vary.Figure 6
Block diagram for weather and climate manipulation.The block diagram for weather and climate control is shown in Figure7, which reflects both open-loop and closed-loop control concepts. A controller is the essential integral part of control system. Since the atmosphere and other climate system components have a large number of positive and negative feedback mechanisms, we present the transfer function for weather and climate control operations taking into account closed-loop concept: (5)
T
c
s
≡
Y
s
X
s
=
G
s
1
+
C
s
1
+
G
s
C
s
H
s
.Figure 7
Blog diagram for feedback climate and weather control.The transfer function technique is a very helpful tool to analyse feedback mechanisms in the ECS, providing important information relevant to the elaboration of climate control systems (e.g., [49–51]).The ECS refers to a class of the distributed parameter systems, which have great difficulty not only in controlling but even in studying, simulating, and analysing. The reason is that mathematical models of the atmospheric and climate processes in spite of their complexities are not sufficiently accurate and adequate; therefore natural climate processes are not really well identified as control systems. To date, the control problem for the ECS processes remains poorly studied due to its immense complexity and originality. To develop a general theoretical framework for optimal control of the ECS we should to take the following into consideration.(i)
The ECS is a spatially distributed system; consequently the control actions for manipulating the ECS should also be distributed in space; however, the implementation of such controls is faintly developed.
(ii)
Processes in the ECS possess enormous energy potential. It is unlikely that we are able to provide control actions whose energy is comparable to the energy of natural climate processes. However, many processes in the ECS are unstable (or metastable) under almost regular conditions. Therefore, infinitesimal perturbations caused by small external forcing may lead to significant changes in the system dynamics. If the system is unstable, the energy barrier is minimal to change qualitatively the system behavior. The identification of sensitive points in which components of the ECS are unstable (metastable) is a critical issue for the development of methods and technique to control and manipulate the ECS processes.
(iii)
Large spatiotemporal scales and huge energy of climate-driven processes impose rigorous requirements for the accuracy and reliability of control systems, since even insignificant errors in control actions can lead to disastrous consequences.
(iv)
Processes in the ECS are interconnected and interdependent; therefore man-made changes in the dynamics of some processes can result in unpredictable and uncontrollable changes in the entire climate system.
(v)
Control actions to perform geoengineering and weather modification operations must be physically feasible and practically executable.Since we are going to explore the weather and climate manipulation as an optimal control problem, let us consider an abstract dynamical system of the form(6)
x
˙
=
f
x
t
,
u
t
x
∈
X
⊆
R
n
,
u
∈
U
⊂
R
m
,
t
∈
0
,
T
⊂
R
+with the initial conditions(7)
x
0
=
x
0
.Suppose that the state vector x belongs to the class of continuously differentiable functions C
1
0
,
T and the control vector u belongs to the class of piecewise continuous functions C
^
0
,
T. We can assume that control vector depends on the state of a system u
t
=
g
t
,
x
t, which means that system (6) is a closed-loop control system, representing the ECS. Suppose also that system (6) is controllable and control vector belongs to the set of admissible controls u
∈
U
⊂
U. It is necessary to emphasize that the set of admissible controls U is determined on the foundation of physical feasibility and technical implementability of methods of weather and climate control. The optimal control aim is to synthesize the control law that ensures the accomplishment of the control objective. Formally an optimal control problem statement can be written as follows.Find the control vectoru
∗
:
0
,
T
↦
U
⊂
U and the orbit x
∗
:
0
,
T
↦
X
⊂
X, which is generated by u
∗, such that the given performance index (cost function) J
x
,
u reaches a minimum (maximum) value; that is,(8)
J
x
,
u
=
f
t
0
,
x
t
0
,
T
,
x
T
+
∫
0
T
F
x
t
,
u
t
,
t
d
t
⟶
m
i
n
.An integrand of the performance indexF
:
X
×
U
→
R. We denote by X a set of state variables on which the cost function is defined. Equation (6) represents the dynamic constraints. Performance index (8) corresponds to the Bolza problem, while the problems of Mayer and Lagrange are its special cases. The formulation of cost function is not trivial and depends on the problem under consideration. Unfortunately, there are no universal approaches regarding how the performance index can be specified. Optimal control problem (8) is nonlinear and can be solved by Pontryagin’s maximum principle, dynamical programming, or classical approaches of the variational calculus.As an example, consider an optimal control problem for the large-scale atmospheric dynamics that is characterized by the slowly moving Rossby waves. These planetary-scale waves strongly influence weather conditions over large geographical regions. Under some assumptions the Rossby wave dynamics at the so-called equivalent-barotropic level of the atmosphere can be described by the following vorticity equation [52, 53]:(9)
∂
t
ζ
+
u
∂
x
ζ
+
f
+
v
∂
y
ζ
+
f
=
f
A
P
ω
0
,where u and v are the components of wind velocity in the directions x and y, respectively, f is the Coriolis parameter, ζ
=
∂
x
v
-
∂
y
u is the vertical component of the relative vorticity, P
=
10
5 Pa is a “standard” surface pressure, A
ξ is an empirical function used to describe the vertical variation of the horizontal wind velocity, ξ
=
p
/
P is a normalized pressure that plays the role of vertical coordinate, and ω
0 is a vertical velocity at the level ξ
=
1, which corresponds to the earth’s surface. Equation (9) is a forced barotropic vorticity equation since its right-hand side is nonzero. The forcing term is determined by the vertical motion at lower boundary ω
0. The natural forcing of the atmospheric Rossby waves is of thermal and orographic origin. These two mechanisms are described by the right-hand side term in (9), that is, by the vertical velocity ω
0. Therefore, this velocity can be considered as a control variable. Equation (9) can be transformed into the inhomogeneous Helmholtz equation [52]: (10)
∇
2
q
-
1
L
0
2
q
=
F
.Here q
=
∂
ψ
/
∂
t, where ψ
=
Φ
/
f is a geostrophic stream function and Φ is the geopotential; 1
/
L
0
2
=
A
R
c
f
/
P is a coefficient that implicitly depends on ω
0; R
c is a correlation coefficient that describes statistical relationship between the surface pressure and the mid-tropospheric geopotential; F
=
-
J
ψ
,
∇
2
ψ
+
f, where J
a
,
b
=
∂
a
/
∂
x
∂
b
/
∂
y
-
∂
a
/
∂
y
∂
b
/
∂
x is Jacobian.Equation (10) is considered in a closed domain D of the x
y-plane with a piecewise continuous boundary ∂
D. The boundary curve in parametric form is represented as x
=
γ
1
σ and y
=
γ
2
σ, where γ
1 and γ
2 are continuous functions of the parameter σ. Suppose that the coordinate origin coincides with the North Pole; then x
=
r
cos
⁡
ϑ, and y
=
r
sin
⁡
ϑ, where r
,
ϑ is a polar coordinate of a generic point on the circle. The solution to (10) is defined on the temporal interval 0
,
T. Specifying ψ
=
const along the boundary ∂
D and introducing the new dependent variable z
=
∂
ψ
/
∂
t
δ
t, where δ
t is the time step of numerical integration, one can obtain the Dirichlet problem:(11)
∇
2
z
-
1
L
0
2
z
=
G
in
D
,
z
=
0
on
∂
D
,where G
=
F
δ
t. Problem (11) is solved numerically to obtain the stream function(12)
ψ
δ
t
=
ψ
0
+
z
t
=
0
,where ψ
0 is the initial condition. Then one can calculate the forecast of geopotential field Φ
=
f
·
ψ
δ
t, which characterizes the dynamics of Rossby waves. Note that the forward difference is applied on the first time step and then a central difference can be employed.Consider the control system characterized by the state variablez
=
z
x
,
y and control variable u
=
u
x
,
y, taking into account (11): (13)
∇
2
z
-
1
L
0
2
z
+
u
=
G
in
D
,
z
=
0
on
∂
D
.Note that the control u represents the additional vertical velocity at the earth’s surface. Let us define the performance index as follows:(14)
J
x
,
u
=
1
2
∬
D
z
-
z
∗
2
+
η
u
2
d
D
,where the pair z
,
u satisfies (13) and z
∗ is the desired spatial distribution of z in the domain D. The optimal control problem is formulated as follows: Find the control
u
∗
∈
U
and the corresponding system state
z
∗
∈
Z
such that cost function (14) is minimized.HereU is a set of admissible controls and Z is a set of some additional constraints on the state variable z. From the standpoint of atmospheric physics, the control problem aims to manipulate the phase speed of Rossby waves by changing the vertical velocity at the earth’s surface. Indeed the linearization of the vorticity equation around the unperturbed westerly zonal flow u
0 and then representation of the solution to the linearized equation in the form of plane waves(15)
ψ
′
x
,
y
,
t
=
R
e
Ψ
0
e
i
k
x
-
c
t
cos
⁡
l
ygive the following expression for the phase speed of Rossby waves:(16)
c
=
u
0
k
2
+
l
2
-
β
-
u
0
′
′
k
2
+
l
2
+
1
/
L
0
2
,where Ψ
0 is the amplitude of perturbation of stream function ψ
′, k and l are the zonal and meridional wavenumbers, β
=
∂
f
/
∂
y is the Rossby parameter, and u
0
′
′
=
d
2
u
0
/
d
y
2. The latest formula shows that the vertical structure of the atmosphere via the term 1
/
L
0
2 contributes to the dispersion properties of Rossby waves. Meanwhile, the term 1
/
L
0
2 depends on the correlation coefficient R
c, which relates the vertical velocity ω
0 with the stream function at the equivalent-barotropic level of the atmosphere. Consequently, the phase velocity of Rossby waves is in a function of ω
0.By introducing new dependent and parametric variables (13) can be written in the normal form [54]:(17)
∂
z
1
∂
x
-
z
2
=
0
,
∂
z
1
∂
y
-
z
3
=
0
,
∂
z
2
∂
x
-
χ
1
=
0
,
∂
z
2
∂
y
-
χ
2
=
0
,
∂
z
3
∂
x
-
χ
2
=
0
,
∂
z
3
∂
y
+
χ
1
+
u
-
1
L
0
2
z
1
-
G
=
0
.Adjoin (17) to the cost function (14) with spatial-varying Lagrange multipliers λ
1, λ
2, λ
3, μ
1, μ
2, and μ
3:(18)
J
A
=
∬
D
1
2
z
1
-
z
∗
2
+
η
u
2
+
λ
1
z
2
-
∂
z
1
∂
x
+
μ
1
z
3
-
∂
z
1
∂
y
+
λ
2
χ
1
-
∂
z
2
∂
x
+
μ
2
χ
2
-
∂
z
2
∂
y
+
λ
3
χ
2
-
∂
z
3
∂
x
+
μ
3
-
χ
1
-
u
+
1
L
0
2
+
G
-
∂
z
3
∂
y
d
D
.The Hamiltonian H associated with the optimal control problem is of the form (19)
H
=
1
2
z
1
-
z
∗
2
+
η
u
2
+
λ
1
z
2
+
λ
2
χ
1
+
λ
3
χ
2
+
μ
1
z
3
+
μ
2
χ
2
+
μ
3
-
χ
1
-
u
+
1
L
0
2
z
1
+
G
.The augmented cost function J
A is then rewritten as(20)
J
A
=
∬
D
H
-
λ
1
∂
z
1
∂
x
-
λ
2
∂
z
2
∂
x
-
λ
3
∂
z
3
∂
x
-
μ
1
∂
z
1
∂
y
-
μ
2
∂
z
2
∂
y
-
μ
3
∂
z
3
∂
y
d
Dor in the vector form(21)
J
A
=
∬
D
H
-
λ
∂
z
∂
x
-
μ
∂
z
∂
y
d
D
,where z
=
z
1
,
z
2
,
z
3, λ
=
λ
1
,
λ
2
,
λ
3, and μ
=
μ
1
,
μ
2
,
μ
3. By adding and subtracting the term z
∂
λ
/
∂
x
+
∂
μ
/
∂
y to the integrand of the latest equation, the augmented cost function J
A takes the form (22)
J
A
=
∬
D
H
+
∂
λ
∂
x
+
∂
μ
∂
y
z
d
D
-
∬
D
∂
∂
x
λ
z
+
∂
∂
y
μ
z
d
D
.Using Green’s theorem, the second integral in this expression becomes a line integral, and then J
A can be rewritten as (23)
J
A
=
∬
D
H
+
∂
λ
∂
x
+
∂
μ
∂
y
z
d
D
+
∫
∂
D
λ
γ
2
′
σ
-
μ
γ
1
′
σ
z
d
σ
.The first variation in J
A caused by variations in u and z is(24)
δ
J
A
=
∫
D
∂
H
∂
z
+
∂
λ
∂
x
+
∂
μ
∂
y
δ
z
+
∂
H
∂
u
1
d
D
+
∫
∂
D
λ
γ
2
′
σ
-
μ
γ
1
′
σ
δ
z
d
σ
,where u
1
=
u
,
χ
1
,
χ
2. We can select both λ and μ to cause the first term in parentheses in the double integral to vanish: (25)
∂
λ
∂
x
+
∂
μ
∂
y
=
-
∂
H
∂
z
.The function z
1 does not change along the boundary ∂
D (z
1
=
0); therefore its variation δ
z
1
=
0 on ∂
D. Functions z
2 and z
3 are not defined on the boundary ∂
D; therefore we can require the relationship λ
γ
2
′
σ
-
μ
γ
1
′
σ
=
0 to hold along ∂
D. Then in polar coordinate we can obtain the boundary conditions for (17):(26)
λ
2
cos
⁡
ϑ
+
μ
2
sin
⁡
ϑ
=
0
,
λ
3
cos
⁡
ϑ
+
μ
3
sin
⁡
ϑ
=
0
.Then the first variation in J
A becomes(27)
δ
J
A
=
∫
D
∂
H
∂
u
1
δ
u
1
d
D
.At the stationary point δ
J
A
=
0, which can be achieved if(28)
∂
H
∂
u
1
=
0
.Equations (25) and (28) with boundary conditions (13) and (26) are the necessary conditions of optimality that can be represented in the expanded form as(29)
∂
λ
1
∂
x
+
∂
μ
1
∂
y
=
-
z
1
-
z
∗
-
μ
3
L
0
2
,
∂
λ
2
∂
x
+
∂
μ
2
∂
y
=
-
λ
1
,
∂
λ
3
∂
x
+
∂
μ
3
∂
y
=
-
μ
1
,
η
u
-
μ
3
=
0
,
λ
2
-
μ
3
=
0
,
λ
3
+
μ
2
=
0
.Thus, to find the control u that minimizes cost function (14) one must solve the set of 12 partial differential equations (17) and (29) with 12 unknowns and with specified boundary conditions:(30)
z
1
∂
D
=
0
,
λ
2
∂
D
=
μ
2
∂
D
=
0
.For particular case of the optimal control problem (η
=
0) we obtain μ
3
=
λ
2
=
0 and the necessary conditions for optimality take the form(31)
∂
λ
1
∂
x
+
∂
μ
1
∂
y
=
-
z
1
-
z
∗
,
∂
μ
2
∂
y
=
-
λ
1
,
∂
λ
3
∂
x
=
-
μ
1
,
λ
3
+
μ
2
=
0
.Since λ
3
=
-
μ
2, the second and third equations in (31) are similar to the first two equations of system (17); therefore, the following assumptions can be made [55]:(32)
λ
1
=
z
3
α
,
μ
1
=
z
2
α
,
λ
3
=
-
μ
2
=
z
1
α
,where α is an arbitrary constant. These assumptions are compatible with boundary conditions (30). The substitution of (32) into the first equation of system (32) gives(33)
∂
2
z
1
∂
x
∂
y
-
∂
2
z
1
∂
x
∂
y
=
-
α
z
1
-
z
∗
.Therefore, we have that α
z
1
-
z
∗
=
0 in the domain D. It follows that the cost function attains its absolute minimum under constraint (13). In order to find the optimal control u that minimizes the cost function, one should solve the following problem in the domain D:(34)
u
=
-
∇
2
z
∗
-
1
L
0
2
z
∗
+
G
0
,
z
∗
∂
D
=
0
.The Rossby waves propagate westward relative to the mean quasi-zonal flow. If the optimal control objective is to reach the stationarity of these waves with respect to the surface by manipulating the lower boundary vertical velocityω
0, we should solve the following equation to obtain the optimal control: u
=
G
0, since in this case z
∗
≡
∂
ψ
/
∂
t
=
0.Thus, we mathematically accurately proved the existence of the absolute minimum of cost function (14) under constraint (13). In the problem considered above, the vertical velocity at the earth’s surface was chosen as control variable. Generally, however, the determination of feasible control variables is not a trivial problem. The “instability-sensitivity” approach suggested in [56] can serve as one of the possible ways to substantiate the choice of control variables in weather and climate control. This method is theoretically based on the fact that some processes in the ECS are inherently unstable under almost regular conditions. Since in mathematical models used in climate studies control actions that manipulate the weather and climate can be expressed via variations in model parameters that act as controls, sensitivity analysis of atmospheric and climate models allows the determination of parameters that can be considered as controls. Using mathematical models that describe various types of natural instability phenomena, one can find the necessary conditions for instability and estimate the influence of various model parameters on the development of instability. The obtained results provide very important information that can be taken into account in designing control systems for atmospheric and climate processes.


## 4. Concluding Remarks

Since the ECS possesses a full spectrum of properties of complex adaptive systems, the exploration of the atmosphere and other components of the ECS as complex adaptive systems allows the consideration of weather and climate manipulation as an optimization problem within the scope of the optimal control theory. However, the complexity of the ECS makes the control problem extremely difficult to solve. In this paper, the necessary conditions of optimality are derived for control of the large-scale atmospheric dynamics. This demonstrates that even a relatively simplified control problem for the atmospheric Rossby waves requires significant efforts to obtain the solution. The study of the ECS as a complex adaptive system, the exploration of its feedback mechanisms, and interactions between climate system components represent the most appropriate approach to estimate the effectiveness and feasibility of potential methods to control the weather and climate. The use of “instability-sensitivity” approach allows the determination of feasible control variables.In our subsequent papers we intend to show how the optimal control theory makes it possible to obtain the mathematically accurate solution to the problem of weather and climate manipulation.



---

*Source: 4615072-2017-01-12.xml*