Introduction. Stochastic physics and climate modelling

Finite computing resources limit the spatial resolution of state-of-the-art global climate simulations to hundreds of kilometres. In neither the atmosphere nor the ocean are small-scale processes such as convection, clouds and ocean eddies properly represented. Climate simulations are known to depend, sometimes quite strongly, on the resulting bulk-formula representation of unresolved processes. Stochastic physics schemes within weather and climate models have the potential to represent the dynamical effects of unresolved scales in ways which conventional bulk-formula representations are incapable of so doing. The application of stochastic physics to climate modelling is a rapidly advancing, important and innovative topic. The latest research findings are gathered together in the Theme Issue for which this paper serves as the introduction.

Comparison of stochastic parameterisation approaches in a single-column model

We discuss and test the potential usefulness of single-column models (SCMs) for the testing of stochastic physics schemes that have been proposed for use in general circulation models (GCMs). We argue that although single column tests cannot be definitive in exposing the full behaviour of a stochastic method in the full GCM, and although there are differences between SCM testing of deterministic and stochastic methods, SCM testing remains a useful tool. It is necessary to consider an ensemble of SCM runs produced by the stochastic method. These can be usefully compared to deterministic ensembles describing initial condition uncertainty and also to combinations of these (with structural model changes) into poor man's ensembles. The proposed methodology is demonstrated using an SCM experiment recently developed by the GCSS (GEWEX Cloud System Study) community, simulating transitions between active and suppressed periods of tropical convection.

Stochastic resonance in a nonlinear model of a rotating, stratified shear flow, with a simple stochastic inertia-gravity wave parameterization

We report on a numerical study of the impact of short, fast inertia-gravity waves on the large-scale, slowly-evolving flow with which they co-exist. A nonlinear quasi-geostrophic numerical model of a stratified shear flow is used to simulate, at reasonably high resolution, the evolution of a large-scale mode which grows due to baroclinic instability and equilibrates at finite amplitude. Ageostrophic inertia-gravity modes are filtered out of the model by construction, but their effects on the balanced flow are incorporated using a simple stochastic parameterization of the potential vorticity anomalies which they induce. The model simulates a rotating, two-layer annulus laboratory experiment, in which we recently observed systematic inertia-gravity wave generation by an evolving, large-scale flow. We find that the impact of the small-amplitude stochastic contribution to the potential vorticity tendency, on the model balanced flow, is generally small, as expected. In certain circumstances, however, the parameterized fast waves can exert a dominant influence. In a flow which is baroclinically-unstable to a range of zonal wavenumbers, and in which there is a close match between the growth rates of the multiple modes, the stochastic waves can strongly affect wavenumber selection. This is illustrated by a flow in which the parameterized fast modes dramatically re-partition the probability-density function for equilibrated large-scale zonal wavenumber. In a second case study, the stochastic perturbations are shown to force spontaneous wavenumber transitions in the large-scale flow, which do not occur in their absence. These phenomena are due to a stochastic resonance effect. They add to the evidence that deterministic parameterizations in general circulation models, of subgrid-scale processes such as gravity wave drag, cannot always adequately capture the full details of the nonlinear interaction.

Supply, demand and a failure of understanding: addressing the culture clash between archaeologists' expectations for training and employment in 'academia' versus 'practice

A university degree is effectively a prerequisite for entering the archaeological workforce in the UK. Archaeological employers consider that new entrants to the profession are insufficiently skilled, and hold university training to blame. But university archaeology departments do not consider it their responsibility to deliver fully formed archaeological professionals, but rather to provide an education that can then be applied in different workplaces, within and outside archaeology. The number of individuals studying archaeology at university exceeds the total number working in professional practice, with many more new graduates emerging than archaeological jobs advertised annually. Over-supply of practitioners is also a contributing factor to low pay in archaeology. Steps are being made to provide opportunities for vocational training, both within and outside the university system, but archaeological training and education within the universities and subsequently the archaeological labour market may be adversely impacted upon by the introduction of variable top-up student fees.

Individual differences in some (but not all) medial prefrontal regions reflect cognitive demand while regulating unpleasant emotion.

The present study investigated the premise that individual differences in autonomic physiology could be used to specify the nature and consequences of information processing taking place in medial prefrontal regions during cognitive reappraisal of unpleasant pictures. Neural (blood oxygenation level-dependent functional magnetic resonance imaging) and autonomic (electrodermal [EDA], pupil diameter, cardiac acceleration) signals were recorded simultaneously as twenty-six older people (ages 64–66 years) used reappraisal to increase, maintain, or decrease their responses to unpleasant pictures. EDA was higher when increasing and lower when decreasing compared to maintaining. This suggested modulation of emotional arousal by reappraisal. By contrast, pupil diameter and cardiac acceleration were higher when increasing and decreasing compared to maintaining. This suggested modulation of cognitive demand. Importantly, reappraisal-related activation (increase, decrease > maintain) in two medial prefrontal regions (dorsal medial frontal gyrus and dorsal cingulate gyrus) was correlated with greater cardiac acceleration (increase, decrease > maintain) and monotonic changes in EDA (increase > maintain > decrease). These data indicate that these two medial prefrontal regions are involved in the allocation of cognitive resources to regulate unpleasant emotion, and that they modulate emotional arousal in accordance with the regulatory goal. The emotional arousal effects were mediated by the right amygdala. Reappraisal-related activation in a third medial prefrontal region (subgenual anterior cingulate cortex) was not associated with similar patterns of change in any of the autonomic measures, thus highlighting regional specificity in the degree to which cognitive demand is reflected in medial prefrontal activation during reappraisal.

A comparative method to evaluate and validate stochastic parametrizations

There is a growing interest in using stochastic parametrizations in numerical weather and climate prediction models. Previously, Palmer (2001) outlined the issues that give rise to the need for a stochastic parametrization and the forms such a parametrization could take. In this article a method is presented that uses a comparison between a standard-resolution version and a high-resolution version of the same model to gain information relevant for a stochastic parametrization in that model. A correction term that could be used in a stochastic parametrization is derived from the thermodynamic equations of both models. The origin of the components of this term is discussed. It is found that the component related to unresolved wave-wave interactions is important and can act to compensate for large parametrized tendencies. The correction term is not proportional to the parametrized tendency. Finally, it is explained how the correction term could be used to give information about the shape of the random distribution to be used in a stochastic parametrization. Copyright © 2009 Royal Meteorological Society

A 3D stochastic cloud model for investigating the radiative properties of inhomogeneous cirrus clouds

(From author). Comments: First 3D stochastic/fractal model of cirrus; first detailed analysis & explanation of power spectra of ice water content, including first observations of 50-km scale break and mixing-induced steepening of spectrum; first demonstration of the potential effect of wind shear on radiative fluxes by changing fall-streak orientation. Has spawned work on the effect of 3D photon transport on the radiative effects of cirrus clouds.