Dynamics II: lecture 9

The daily observed maximum and minimum temperatures is often compared to the “normal” temperatures based upon the 30-year average. Climate averages provide a context for something like “this winter will be wetter (or drier, or colder, or warmer, etc.) than normal.

The daily observed maximum and minimum temperatures is often compared to the “normal” temperatures based upon the 30-year average. Climate averages provide a context for something like “this winter will be wetter (or drier, or colder, or warmer, etc.) than normal.

It has been said "Climate is what you expect. Weather is what you get.”

Here: general approach to the mean and fluctuations in the climate system. Indeed, the Brownian motion approach is a helpful analogue for weather and climate.

The stuff the world is made of

World views: Elementary particles, quantum mechanics, theory of relativity

Limit of divisibility (Democritus, Aristotle: Matter is not a continuous whole: “The world cannot be composed of infinitely small particles”)

The Roman Lucretius's scientific poem

“You will see a multitude of tiny particles mingling in a multitude of ways… their dancing is an actual indication of underlying movements of matter that are hidden from our sight… It originates with the atoms which move of themselves”

Although the mingling motion of dust particles is caused largely by air currents, the glittering, tumbling motion of small dust particles is, indeed, caused chiefly by true Brownian dynamics.

Jan Ingenhousz had described the irregular motion of coal dust particles on the surface of alcohol in 1785.

Botanist Robert Brown in 1827: studying pollen particles floating in water under the microscope. He then observed minute particles within the vacuoles of the pollen grains executing a jittery motion. By repeating the experiment with particles of dust, he was able to rule out that the motion was due to pollen particles being 'alive', although the origin of the motion was yet to be explained.

Film Brown

Thorvald N. Thiele in 1880 in a paper on the method of least squares. This was followed independently by Louis Bachelier in 1900 “The theory of speculation”, stock and option markets.

However, it was Albert Einstein's (in his 1905 paper) and Marian Smoluchowski's (1906) independent research of the problem that brought the solution to the attention of physicists, and presented it as a way to indirectly confirm the existence of atoms and molecules.

The confirmation of Einstein's theory constituted empirical progress for the kinetic theory of heat. In essence, Einstein showed that the motion can be predicted directly from the kinetic model of thermal equilibrium. The importance of the theory lay in the fact that it confirmed the kinetic theory's account of the second law of thermodynamics as being an essentially statistical law.

$$\rho(x, t+\tau) = \rho(x, t) \cdot 1 + 0 + \frac{\partial^2 \rho}{\partial x^2} \cdot \int_{-\infty}^{+\infty} \frac{\Delta^2}{2} \cdot \phi(\Delta) \, \mathrm{d} \Delta \quad + \quad ...$$

The integral in the first term is equal to one by the definition of probability, and the second and other even terms (i.e. first and other odd moments) vanish because of space symmetry.

$$\frac{\partial\rho}{\partial t} = \frac{\partial^2 \rho}{\partial x^2} \cdot \int_{-\infty}^{+\infty} \frac{\Delta^2}{2\, \tau} \cdot \phi(\Delta) \, \mathrm{d} \Delta + \mathrm{higher\;order\; even\; moments}$$

Where the second moment of probability of displacement $\Delta,$ is interpreted as diffusivity $$D = \int \frac{\Delta^2}{2\, \tau} \cdot \phi(\Delta) \, \mathrm{d} \Delta$$

Then the density of Brownian particles $\rho$ at point x at time t satisfies the diffusion equation: $$\frac{\partial\rho}{\partial t}=D\cdot \frac{\partial^2\rho}{\partial x^2}$$

Density of Brownian particles $\rho$ satisfies the diffusion equation: $$\frac{\partial\rho}{\partial t}=D\cdot \frac{\partial^2\rho}{\partial x^2}$$

Assuming that N particles start from the origin at the initial time $$\rho(x,t)=\frac{N}{\sqrt{4\pi Dt}}e^{-\frac{x^2}{4Dt}}.$$

This expression allowed Einstein to calculate the moments directly. The first moment is seen to vanish, meaning that the Brownian particle is equally likely to move to the left as it is to move to the right. The second moment is, however, non-vanishing, being given by

$$\overline{x^2}=2\,D\,t.$$

This expresses the mean squared displacement in terms of the time elapsed and the diffusivity. From this expression Einstein argued that the displacement of a Brownian particle is not proportional to the elapsed time, but rather to its square root.

In a stochastic framework of climate theory $$\frac{d}{dt} x(t) = f(x) + g(x) \xi ,$$ where $\xi$ is a stationary stochastic process and the functions f, g describe the climate dynamics. The properties of the random force are described through its distribution and its correlation properties at different times.

The process $\xi$ is assumed to have a Gaussian distribution of zero average, $$< \xi (t) > = 0$$ and to be $\delta$-correlated in time, $$< \xi(t) \xi(t +\tau) > = \delta (\tau) \label{deltacorr}$$ where $\delta$ is the delta function defined by $$\int f (x) \, \delta (x-x_0) \, dx = f(x_0) \quad .$$

The brackets indicate an average over realizations of the random force.

Using the ergodic hypothesis, the ensemble average
$$\left\langle \qquad \right\rangle$$

can be expressed as the time average $$\quad \lim_{T \to \infty} \frac{1}{T} \int_{-T/2}^{T/2} dt$$

For a Gaussian process only the average and second moment need to be specified since all higher moments can be expressed in terms of the first two. Note that the dependence of the correlation function on the time difference $\tau$ assumes that $\xi$ is a stationary process.

Imagine a so-called red-noise process $$\frac{dx}{dt} = -\lambda x + \xi \,$$