Introductory remarks

The very first numerical models which were developed more than 20 years ago were drastic simplifications of the real atmosphere and they were mostly restricted to describe adiabatic processes. For prediction of a day or two of the mid tropospheric flow these models often gave reasonable results but the result deteriorated quickly when the prediction was extended further in time. The prediction of the surface flow was unsatisfactory even for short predictions. It was evident that both the energy generating processes as well as the dissipative processes have to be included in numerical models in order to predict the weather patterns in the lower part of the atmosphere and to predict the atmosphere in general beyond a day or two. Present-day computers make it possible to attack the weather forecasting problem in a more comprehensive and complete way and substantial efforts have been made during the last decade in particular to incorporate the non-adiabatic processes in numerical prediction models. The physics of radiational transfer, condensation of moisture, turbulent transfer of heat, momentum and moisture and the dissipation of kinetic energy are the most important processes associated with the formation of energy sources and sinks in the atmosphere and these have to be incorporated in numerical prediction models extended over more than a few days. The mechanisms of these processes are mainly related to small scale disturbances in space and time or even molecular processes. It is therefore one of the basic characteristics of numerical models that these small scale disturbances cannot be included in an explicit way. The reason for this is the discretization of the model's atmosphere by a finite difference grid or the use of a Galerkin or spectral function representation. The second reason why we cannot explicitly introduce these processes into a numerical model is due to the fact that some physical processes necessary to describe them (such as the local buoyance) are a priori eliminated by the constraints of hydrostatic adjustment. Even if this physical constraint can be relaxed by making the models non-hydrostatic the scale problem is virtually impossible to solve and for the foreseeable future we have to try to incorporate the ensemble or gross effect of these physical processes on the large scale synoptic flow. The formulation of the ensemble effect in terms of grid-scale variables (the parameters of the large-scale flow) is called 'parameterization'. For short range prediction of the synoptic flow at middle and high latitudes, very simple parameterization has proven to be rather successful.

The angular momentum constraint on climate sensitivity and downward influence in the middle atmosphere

It is shown that under reasonable assumptions, conservation of angular momentum provides a strong constraint on gravity wave drag feedbacks to radiative perturbations in the middle atmosphere. In the time mean, radiatively induced temperature perturbations above a given altitude z cannot induce changes in zonal mean wind and temperature below z through feedbacks in gravity wave drag alone (assuming an unchanged gravity wave source spectrum). Thus, despite the many uncertainties in the parameterization of gravity wave drag, the role of gravity wave drag in middle-atmosphere climate perturbations may be much more limited than its role in climate itself. This constraint limits the possibilities for downward influence from the mesosphere. In order for a gravity wave drag parameterization to respect the momentum constraint and avoid spurious downward influence, any nonzero parameterized momentum flux at a model lid must be deposited within the model domain, and there must be no zonal mean sponge layer. Examples are provided of how violation of these conditions leads to spurious downward influence. For planetary waves, the momentum constraint does not prohibit downward influence, but it limits the mechanisms by which it can occur: in the time mean, downward influence from a radiative perturbation can only arise through changes in reflection and meridional propagation properties of planetary waves.

Energetics of a symmetric circulation including momentum constraints

A theory of available potential energy (APE) for symmetric circulations, which includes momentum constraints, is presented. The theory is a generalization of the classical theory of APE, which includes only thermal constraints on the circulation. Physically, centrifugal potential energy is included along with gravitational potential energy. The generalization relies on the Hamiltonian structure of the conservative dynamics, although (as with classical APE) it still defines the energetics in a nonconservative framework. It follows that the theory is exact at finite amplitude, has a local form, and can be applied to a variety of fluid models. It is applied here to the f -plane Boussinesq equations. It is shown that, by including momentum constraints, the APE of a symmetrically stable flow is zero, while the energetics of a mechanically driven symmetric circulation properly reflect its causality.

When long-range zero-lag synchronization is feasible in cortical networks

Many studies have reported long-range synchronization of neuronal activity between brain areas, in particular in the beta and gamma bands with frequencies in the range of 14–30 and 40–80 Hz, respectively. Several studies have reported synchrony with zero phase lag, which is remarkable considering the synaptic and conduction delays inherent in the connections between distant brain areas. This result has led to many speculations about the possible functional role of zero-lag synchrony, such as for neuronal communication, attention, memory, and feature binding. However, recent studies using recordings of single-unit activity and local field potentials report that neuronal synchronization may occur with non-zero phase lags. This raises the questions whether zero-lag synchrony can occur in the brain and, if so, under which conditions. We used analytical methods and computer simulations to investigate which connectivity between neuronal populations allows or prohibits zero-lag synchrony. We did so for a model where two oscillators interact via a relay oscillator. Analytical results and computer simulations were obtained for both type I Mirollo–Strogatz neurons and type II Hodgkin–Huxley neurons. We have investigated the dynamics of the model for various types of synaptic coupling and importantly considered the potential impact of Spike-Timing Dependent Plasticity (STDP) and its learning window. We confirm previous results that zero-lag synchrony can be achieved in this configuration. This is much easier to achieve with Hodgkin–Huxley neurons, which have a biphasic phase response curve, than for type I neurons. STDP facilitates zero-lag synchrony as it adjusts the synaptic strengths such that zero-lag synchrony is feasible for a much larger range of parameters than without STDP.

Quantification of chemical and physical processes influencing ozone during long-range transport using a trajectory ensemble

During long-range transport, many distinct processes – including photochemistry, deposition, emissions and mixing – contribute to the transformation of air mass composition. Partitioning the effects of different processes can be useful when considering the sensitivity of chemical transformation to, for example, a changing environment or anthropogenic influence. However, transformation is not observed directly, since mixing ratios are measured, and models must be used to relate changes to processes. Here, four cases from the ITCT-Lagrangian 2004 experiment are studied. In each case, aircraft intercepted a distinct air mass several times during transport over the North Atlantic, providing a unique dataset and quantifying the net changes in composition from all processes. A new framework is presented to deconstruct the change in O3 mixing ratio (Δ O3) into its component processes, which were not measured directly, taking into account the uncertainty in measurements, initial air mass variability and its time evolution. The results show that the net chemical processing (Δ O3chem) over the whole simulation is greater than net physical processing (Δ O3phys) in all cases. This is in part explained by cancellation effects associated with mixing. In contrast, each case is in a regime of either net photochemical destruction (lower tropospheric transport) or production (an upper tropospheric biomass burning case). However, physical processes influence O3 indirectly through addition or removal of precursor gases, so that changes to physical parameters in a model can have a larger effect on Δ O3chem than Δ O3phys. Despite its smaller magnitude, the physical processing distinguishes the lower tropospheric export cases, since the net photochemical O3 change is −5 ppbv per day in all three cases. Processing is quantified using a Lagrangian photochemical model with a novel method for simulating mixing through an ensemble of trajectories and a background profile that evolves with them. The model is able to simulate the magnitude and variability of the observations (of O3, CO, NOy and some hydrocarbons) and is consistent with the time-average OH following air-masses inferred from hydrocarbon measurements alone (by Arnold et al., 2007). Therefore, it is a useful new method to simulate air mass evolution and variability, and its sensitivity to process parameters.

Lightning-induced extensive charge sheets provide long range electrostatic thunderstorm detection

By combining electrostatic measurements of lightning-induced electrostatic field changes with radio frequency lightning location, some field changes from exceptionally distant lightning events are apparent which are inconsistent with the usual inverse cube of distance. Furthermore, by using two measurement sites, a transition zone can be identified beyond which the electric field response reverses polarity. For these severe lightning events, we infer a horizontally extensive charge sheet above a thunderstorm, consistent with a mesospheric halo of several hundred kilometers’ extent.

Tracking the momentum flux of a CME and quantifying its influence on geomagnetically induced currents at Earth

We investigate a coronal mass ejection (CME) propagating toward Earth on 29 March 2011. This event is specifically chosen for its predominately northward directed magnetic field, so that the influence from the momentum flux onto Earth can be isolated. We focus our study on understanding how a small Earth-directed segment propagates. Mass images are created from the white-light cameras onboard STEREO which are also converted into mass height-time maps (mass J-maps). The mass tracks on these J-maps correspond to the sheath region between the CME and its associated shock front as detected by in situ measurements at L1. A time series of mass measurements from the STEREO COR-2A instrument is made along the Earth propagation direction. Qualitatively, this mass time series shows a remarkable resemblance to the L1 in situ density series. The in situ measurements are used as inputs into a three-dimensional (3-D) magnetospheric space weather simulation from the Community Coordinated Modeling Center. These simulations display a sudden compression of the magnetosphere from the large momentum flux at the leading edge of the CME, and predictions are made for the time derivative of the magnetic field (dB/dt) on the ground. The predicted dB/dt values were then compared with the observations from specific equatorially located ground stations and showed notable similarity. This study of the momentum of a CME from the Sun down to its influence on magnetic ground stations on Earth is presented as a preliminary proof of concept, such that future attempts may try to use remote sensing to create density and velocity time series as inputs to magnetospheric simulations.

Simulations of angular momentum evolution in wind-fed CV precursors

It is thought that the secondary stars in cataclysmic variables (CVs) may undergo a period of mass loss in the form of a wind during the evolution of the system (Mullan et al. 1992). This wind is thought to magnetically brake the secondary star with a time-scale ~ 10^8 yr (e.g. van Paradijs 1986). When the secondary’s spin has been brought close to synchronism with the orbit it is possible for tidal torques to lock the secondary in synchronous rotation.