Predictability of the rotating annulus

During my DPhil I worked with the thermally-driven rotating annulus (right), a laboratory analogue for a generic planetary atmosphere. We can use the annulus to study methods for forecasting or data assimilation in current use or in development using a real fluid with a non-idealised model under laboratory conditions. The laboratory is a bridge between analytical systems, where new methods are first tested, and large atmospheric models, where they are eventually applied.

Baroclinic instability is an important mechanism for the development of large-scale waves and eddies in the atmospheres of the Earth (left) and other planets such as Mars. These waves and eddies are associated with the large-scale transport of heat and momentum; their chaotic evolution in the Earth's atmosphere determines the fundamental limit on how far ahead it is possible to forecast weather and (to a lesser degree) the climate. On Mars, however, baroclinic weather systems are significantly more coherent, and at times its atmosphere appears to be more predictable than the Earth's.

This research had three main aims:

  • Measure the predictability of fluid flow in the rotating annulus, particularly under experimental conditions broadly equivalent to the atmospheres of Earth and Mars. This was done by applying common weather forecasting techniques to the laboratory experiment.
  • Identify what physical mechanisms cause the predictions to fail over time and how this depends on the type and complexity of the flow.
  • Improve weather forecasting techniques by using the laboratory experiment under controlled and reproducible conditions to test methods used by weather and climate forecasters.
  • Of the two kinds of predictability identified by Lorenz (1975), predictability of the "second kind" - i.e. the "climate" problem: determining the type of flow observed for a particular experimental setup - had already been well characterised for the rotating annulus over several decades. Predictability of the "first kind", however - i.e. the "weather" problem: predicting the evolution of a particular flow over time from a given initial state - had not been well explored until this work.

    The intrinsic predictability of the most important quasi-periodic and chaotic flow regimes observed in the laboratory was investigated using numerical simulations with MORALS, the Met Office/Oxford Rotating Annulus Laboratory Simulation (Young & Read, 2008a). To determine the most unstable flow patterns we used the breeding vector method, which is a well-established technique for evaluating the intrinsic predictability of the state of the Earth's atmosphere and has been used by several operational weather forecasting centres worldwide in their everyday weather forecasts. The mechanisms causing predictability to break down were then studied using ensemble prediction techniques (Young & Read, 2008b). We used the analysis correction method to assimilate data into the model (Young & Read, 2013) in order to create initial states for our prediction experiments (Young, 2010, Young & Read, 2016).

    The laboratory analogue

    The rotating annulus is a classic laboratory experiment that has been used since the 1950s to study fluid flows whose dynamics are similar to airflow in the midlatitudes of planetary atmospheres.

    The experiment rotates about a vertical axis and is heated at the vertical boundaries of the fluid, generating a horizontal temperature gradient across the fluid. This is similar to the atmosphere, which is heated by the Sun more at low latitudes than at high latitudes (below), and also rotates about an axis. By changing these two parameters in the experiment a wide range of behaviour can be seen, from simple waves (right) to turbulence and chaotic flow. In this sense it is similar to the Earth's atmosphere.

    Left: The atmospheric analogue to the laboratory experiment [NASA]. Right: Composite image of the South Pole showing planetary-scale waves similar to those that can be produced in the annulus [NASA/Galileo/PIA00729].

    Using the annulus avoids much of the complicated behaviour seen in the atmosphere, such as the difference between land and sea, the chemistry of different molecules in the air, and the biosphere. These effects are secondary to rotation, gravity, and heating, although some later experiments have included them. This simplicity makes it easier to study the physics behind weather and climate. The laboratory context also allows interpretation of the system behaviour without the difficulties associated with imperfect and incomplete observations of atmospheric circulation, and avoids some of the complexities resulting from a wide range of interacting scales and sub-systems.

    Data assimilation

    In Young & Read (2013) we used the analysis correction method to assimilate laboratory measurements of fluid velocities into MORALS. The dataset containing these measurements is available via ORA-Data. We analyzed the 2S (steady flow, azimuthal wavenumber 2) and 3AV (amplitude vacillating) regular flow regimes, and the 3SV (structural vacillation) weakly chaotic flow regime at higher rotation rates.

    The assimilated horizontal velocities were combined with the model to produce model states that were both consistent with our physical understanding (as represented by the numerical model) and consistent with the observations. As well as showing that data assimilation worked in this context, we examined a number of specific scenarios:

  • A wavenumber transition between two flow regimes (below)
  • Using assimilation to propagate information from data-rich to data-poor regions
  • The response of the assimilation to a strong disturbance in the flow
  • A vortex-shedding instability phenomenon at the highest high rotation rate (below)
  • By using unobserved variables that are only available via the assimilation (such as temperature and the vertical structure of the velocity field), we found that the mean azimuthal flow is weakened by the action of eddies, which partly explained why vortices were shed at the very highest rotation rates but not at lower rotation.

    Assimilation across a wavenumber transition from 3AV to 2S flow. Top: Horizontal velocity analysis (blue) and observations (red). Bottom: Azimuthal sections at mid-height/mid-radius; analysis is blue, observations are black, and a free-running simulation with no assimilation is red (Young & Read, 2013).

    Sequences of analyses showing vortex shedding at high rotation rate in 3SV flow. Colours show stream function and arrows are observations. Frames are separated by 10s. (Young & Read, 2013).

    Predictability

    In Young & Read (2016) we used our ensemble prediction system with breeding vectors and data assimlation to forecast the rotating annulus in steady 2S, amplitude vacillating 3AV, and structurally vacillating 3SV flow regimes, verifying the forecasts against laboratory data.

    Block diagram showing stages of the annulus ensemble prediction scheme (Young & Read, 2016).

    The steady wave flow regime was predictable until the end of the available data. Forecasts in the amplitude and structural vacillation flow regimes lost quality and skill by a combination of wave drift and wavenumber transition. Amplitude vacillation was predictable up to several hundred seconds ahead, and structural vacillation was predictable for a few hundred seconds. The wavenumber transitions were partly explained by hysteresis in the rotating annulus experiment and model.

    Sequence of spaghetti plots showing a forecast in the 3AV flow regime. Red vectors show observations, and contours show individual ensemble members. Frames are separated by 100s. The forecast ensemble splits into two parts, one remaining wavenumber 3 while the other transitions to wavenumber 2. (Young & Read, 2016).

    Example of breeding vectors used to start an ensemble forecast in the 3AV flow regime. These represent the uncertainty in the flow at each point in the flow at the start of the forecast (Young & Read, 2016).


    Examples of ensemble forecasts showing radial velocity at a single point in the flow for 2S (top) and 3AV (bottom) flow regimes. Black lines show forecast ensemble members, dot-dashed line is ensemble mean, and observations (interpolated to the model grid point shown) are dotted with grey bars indicating observational error. The 2S forecast has very small spread but tracks the observations well. The 3AV forecast splits into two after about 1000s as some of the ensemble members transition from wavenumber 3 to wavenumber 2 flow (Young & Read, 2016).