Analysis of convectively-generated gravity waves in mesoscale model simulations and wind profiler observations

The characteristics of convectively-generated gravity waves during an episode of deep convection near the coast of Wales are examined in both high resolution mesoscale simulations [with the (UK) Met Oce Unified Model] and in observations from a Mesosphere-Stratosphere-Troposphere (MST) wind profiling Doppler radar. Deep convection reached the tropopause and generated vertically propagating, high frequency waves in the lower stratosphere that produced vertical velocity perturbations O(1 m/s). Wavelet analysis is applied in order to determine the characteristic periods and wavelengths of the waves. In both the simulations and observations, the wavelet spectra contain several distinct preferred scales indicated by multiple spectral peaks. The peaks are most pronounced in the horizontal spectra at several wavelengths less than 50 km. Although these peaks are most clear and of largest amplitude in the highest resolution simulations (with 1 km horizontal grid length), they are also evident in coarser simulations (with 4 km horizontal grid length). Peaks also exist in the vertical and temporal spectra (between approximately 2.5 and 4.5 km, and 10 to 30 minutes, respectively) with good agreement between simulation and observation. Two-dimensional (wavenumber-frequency) spectra demonstrate that each of the selected horizontal scales contains peaks at each of preferred temporal scales revealed by the one- dimensional spectra alone.

An objective climatology of the dynamical forcing of polar lows in the Nordic seas

A new objective climatology of polar lows in the Nordic (Norwegian and Barents) seas has been derived from a database of diagnostics of objectively identified cyclones spanning the period January 2000 to April 2004. There are two distinct parts to this study: the development of the objective climatology and a characterization of the dynamical forcing of the polar lows identified. Polar lows are an intense subset of polar mesocyclones. Polar mesocyclones are distinguished from other cyclones in the database as those that occur in cold air outbreaks over the open ocean. The difference between the wet-bulb potential temperature at 700 hPa and the sea surface temperature (SST) is found to be an effective discriminator between the atmospheric conditions associated with polar lows and other cyclones in the Nordic seas. A verification study shows that the objective identification method is reliable in the Nordic seas region. After demonstrating success at identifying polar lows using the above method, the dynamical forcing of the polar lows in the Nordic seas is characterized. Diagnostics of the ratio of mid-level vertical motion attributable to quasi-geostrophic forcing from upper and lower levels (U/L ratio) are used to determine the prevalence of a recently proposed category of extratropical cyclogenesis, type C, for which latent heat release is crucial to development. Thirty-one percent of the objectively identified polar low events (36 from 115) exceeded the U/L ratio of 4.0, previously identified as a threshold for type C cyclones. There is a contrast between polar lows to the north and south of the Nordic seas. In the southern Norwegian Sea, the population of polar low events is dominated by type C cyclones. These possess strong convection and weak low-level baroclinicity. Over the Barents and northern Norwegian seas, the well-known cyclogenesis types A and B dominate. These possess stronger low-level baroclinicity and weaker convection.

Temperature and moist-stability effects on midlatitude orographic precipitation

Idealized, convection-resolving simulations of moist orographic flows are conducted to investigate the influence of temperature and moist stability on the drying ratio (DR), defined as the fraction of the impinging water mass removed as orographic precipitation. In flow past a long ridge, where most of the air rises over the barrier rather than detouring around it, DR decreases as the surface temperature (Ts) increases, even as the orographic cap cloud becomes statically unstable at higher Ts and develops embedded convection. This behaviour is explained by a few physical principles: (1) the Clausius–Clapeyron equation dictates that the normalized condensation rate decreases as the flow gets warmer, (2) the replacement of ice-phase precipitation growth with warm-rain processes decreases the efficiency by which condensate is converted to precipitation, thereby lowering precipitation efficiency, and (3) embedded convection acts more to vertically redistribute moisture than to enhance precipitation. Over an isolated mountain, the effects of (1) and (2) are counteracted by moisture deflection around the barrier, which is stronger in the colder, more stable flows.

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.

Influence of south Atlantic sea surface temperatures on rainfall variability and extremes over southern Africa

It is generally agreed that changing climate variability, and the associated change in climate extremes, may have a greater impact on environmentally vulnerable regions than a changing mean. This research investigates rainfall variability, rainfall extremes, and their associations with atmospheric and oceanic circulations over southern Africa, a region that is considered particularly vulnerable to extreme events because of numerous environmental, social, and economic pressures. Because rainfall variability is a function of scale, high-resolution data are needed to identify extreme events. Thus, this research uses remotely sensed rainfall data and climate model experiments at high spatial and temporal resolution, with the overall aim being to investigate the ways in which sea surface temperature (SST) anomalies influence rainfall extremes over southern Africa. Extreme rainfall identification is achieved by the high-resolution microwave/infrared rainfall algorithm dataset. This comprises satellite-derived daily rainfall from 1993 to 2002 and covers southern Africa at a spatial resolution of 0.1° latitude–longitude. Extremes are extracted and used with reanalysis data to study possible circulation anomalies associated with extreme rainfall. Anomalously cold SSTs in the central South Atlantic and warm SSTs off the coast of southwestern Africa seem to be statistically related to rainfall extremes. Further, through a number of idealized climate model experiments, it would appear that both decreasing SSTs in the central South Atlantic and increasing SSTs off the coast of southwestern Africa lead to a demonstrable increase in daily rainfall and rainfall extremes over southern Africa, via local effects such as increased convection and remote effects such as an adjustment of the Walker-type circulation.

Flight periodicity and the vertical distribution of high-altitude moth migration over southern Britain

The continuous operation of insect-monitoring radars in the UK has permitted, for the first time, the characterization of various phenomena associated with high-altitude migration of large insects over this part of northern Europe. Previous studies have taken a case-study approach, concentrating on a small number of nights of particular interest. Here, combining data from two radars, and from an extensive suction- and light-trapping network, we have undertaken a more systematic, longer-term study of diel flight periodicity and vertical distribution of macro-insects in the atmosphere. Firstly, we identify general features of insect abundance and stratification, occurring during the 24-hour cycle, which emerge from four years’ aggregated radar data for the summer months in southern Britain. These features include mass emigrations at dusk and to a lesser extent at dawn, and daytime concentrations associated with thermal convection. We then focus our attention on the well-defined layers of large nocturnal migrants that form in the early evening, usually at heights of 200–500 m above ground. We present evidence from both radar and trap data that these nocturnal layers are composed mainly of noctuid moths, with species such as Noctua pronuba, Autographa gamma, Agrotis exclamationis, A. segetum, Xestia c-nigrum and Phlogophora meticulosa predominating.

Uncertainties in future projections of extreme precipitation in the Indian monsoon region

Uncertainties in changes to the spatial distribution and magnitude of the heaviest extremes of daily monsoon rainfall over India are assessed in the doubled CO2 climate change scenarios in the IPCC Fourth Assessment Report. Results show diverse changes to the spatial pattern of the 95th and 99th subseasonal percentiles, which are strongly tied to the mean precipitation change during boreal summer. In some models, the projected increase in heaviest rainfall over India at CO2 doubling is entirely predictable based upon the surface warming and the Clausius–Clapeyron relation, a result which may depend upon the choice of convection scheme. Copyright © 2009 Royal Meteorological Society and Crown Copyright

Boundary integral methods in high frequency scattering

In this article we review recent progress on the design, analysis and implementation of numerical-asymptotic boundary integral methods for the computation of frequency-domain acoustic scattering in a homogeneous unbounded medium by a bounded obstacle. The main aim of the methods is to allow computation of scattering at arbitrarily high frequency with finite computational resources.

Mean climate and transience in the tropics of the UGAMP GCM - sensitivity to convective parametrization

The sensitivity of the UK Universities Global Atmospheric Modelling Programme (UGAMP) General Circulation Model (UGCM) to two very different approaches to convective parametrization is described. Comparison is made between a Kuo scheme, which is constrained by large-scale moisture convergence, and a convective-adjustment scheme, which relaxes to observed thermodynamic states. Results from 360-day integrations with perpetual January conditions are used to describe the model's tropical time-mean climate and its variability. Both convection schemes give reasonable simulations of the time-mean climate, but the representation of the main modes of tropical variability is markedly different. The Kuo scheme has much weaker variance, confined to synoptic frequencies near 4 days, and a poor simulation of intraseasonal variability. In contrast, the convective-adjustment scheme has much more transient activity at all time-scales. The various aspects of the two schemes which might explain this difference are discussed. The particular closure on moisture convergence used in this version of the Kuo scheme is identified as being inappropriate.

Condition number estimates for combined potential boundary integral operators in acoustic scattering

We study the classical combined field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle, namely the indirect formulation due to Brakhage-Werner/Leis/Panic, and the direct formulation associated with the names of Burton and Miller. We obtain lower and upper bounds on the condition numbers for these formulations, emphasising dependence on the frequency, the geometry of the scatterer, and the coupling parameter. Of independent interest we also obtain upper and lower bounds on the norms of two oscillatory integral operators, namely the classical acoustic single- and double-layer potential operators.

A review of the theoretical basis for bulk mass flux convective parameterization

Most parameterizations for precipitating convection in use today are bulk schemes, in which an ensemble of cumulus elements with different properties is modelled as a single, representative entraining-detraining plume. We review the underpinning mathematical model for such parameterizations, in particular by comparing it with spectral models in which elements are not combined into the representative plume. The chief merit of a bulk model is that the representative plume can be described by an equation set with the same structure as that which describes each element in a spectral model. The equivalence relies on an ansatz for detrained condensate introduced by Yanai et al. (1973) and on a simplified microphysics. There are also conceptual differences in the closure of bulk and spectral parameterizations. In particular, we show that the convective quasi-equilibrium closure of Arakawa and Schubert (1974) for spectral parameterizations cannot be carried over to a bulk parameterization in a straightforward way. Quasi-equilibrium of the cloud work function assumes a timescale separation between a slow forcing process and a rapid convective response. But, for the natural bulk analogue to the cloud-work function (the dilute CAPE), the relevant forcing is characterised by a different timescale, and so its quasi-equilibrium entails a different physical constraint. Closures of bulk parameterization that use the non-entraining parcel value of CAPE do not suffer from this timescale issue. However, the Yanai et al. (1973) ansatz must be invoked as a necessary ingredient of those closures.

Further characterization of muscarinic agonist-induced epileptiform bursting activity in immature rat piriform cortex, in vitro

The characteristics of muscarinic acetylcholine receptor agonist-induced epileptiform bursting seen in immature rat piriform cortex slices in vitro were further investigated using intracellular recording, with particular focus on its postnatal age-dependence (P+14-P+30), pharmacology, site(s) of origin and the likely contribution of the muscarinic acetylcholine receptor agonist-induced post-stimulus slow afterdepolarization and gap junction functionality toward its generation. The muscarinic agonist, oxotremorine-M (10 microM), induced rhythmic bursting only in immature piriform cortex slices; however, paroxysmal depolarizing shift amplitude, burst duration and burst incidence were inversely related to postnatal age. No significant age-dependent changes in neuronal membrane properties or postsynaptic muscarinic responsiveness accounted for this decline. Burst incidence was higher when recorded in anterior and posterior regions of the immature piriform cortex. In adult and immature neurones, oxotremorine-M effects were abolished by M1-, but not M2-muscarinic acetylcholine receptor-selective antagonists. Rostrocaudal lesions, between piriform cortex layers I and II, or layer III and endopiriform nucleus in adult or immature slices did not influence oxotremorine-M effects; however, the slow afterdepolarization in adult (but not immature) lesioned slices was abolished. Gap junction blockers (carbenoxolone or octanol) disrupted muscarinic bursting and diminished the slow afterdepolarization in immature slices, suggesting that gap junction connectivity was important for bursting. Our data show that neural networks within layers II-III function as primary oscillatory circuits for burst initiation in immature rat piriform cortex during persistent muscarinic receptor activation. Furthermore, we propose that muscarinic slow afterdepolarization induction and gap junction communication could contribute towards the increased epileptiform susceptibility of this brain area.

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.

Jupiter's and Saturn's convectively driven banded jets in the laboratory

The banded patterns of cloud and wind are among the most striking features of the atmospheres of Jupiter and Saturn, but their dynamical origin remains poorly understood. Most approaches towards understanding zonation so far (also in the terrestrial oceans) have used highly idealized models to show that it might originate from dynamical anisotropy in a shallow turbulent fluid layer due to the planetary β-effect. Here we report the results of laboratory experiments, conducted on a 14-m diameter turntable, which quantitatively confirm that multiple zonal jets may indeed be generated and maintained by this mechanism in the presence of deep convection and a topographic β-effect. At the very small values of Ekman number (≤2 × 10−5) and large local Reynolds numbers (≥2000, based on jet scales) achieved, the kinetic energy spectra suggest the presence of both energy-cascading and enstrophy-cascading inertial ranges in addition to the zonation near twice the Rhines wave number.

On the dependence of boundary-layer ventilation on frontal type

Case study simulations with idealized tracers have been used to determine the relationship between the dynamics and conceptual representations of different midlatitude frontal systems and the amount, distribution, and time scale of boundary layer ventilation by these systems. The key features of ventilation by a kata– and ana–cold frontal system are found to be quantitatively and also often qualitatively similar to the main ventilation pathways, which are the conveyor belts, cloud head, and other convection. The conveyor belts and cloud head occur within cloud, implying that they can be identified using satellite imagery. Differences in the transport by the two systems can be related to their conceptual representations and include a sensitive dependence on the diurnal cycle for the kata- but not the ana-cold frontal case.