The propagation mechanism of a vortex on the β plane

The propagation velocity and propagation mechanism for vortices on a β plane are determined for a reduced-gravity model by integrating the momentum equations over the β plane. Isolated vortices, vortices in a background current, and initial vortex propagation from rest are studied. The propagation mechanism for isolated anticyclones as well as cyclones, which has been lacking up to now, is presented. It is shown that, to first order, the vortex moves to generate a Coriolis force on the mass anomaly of the vortex to compensate for the force on the vortex due to the variation of the Coriolis parameter. Only the mass anomaly of the vortex is of importance, because the Coriolis force due to the motion of the bulk of the layer moving with the vortex is almost fully compensated by the Coriolis force on the motion of the exterior flow. Because the mass anomaly of a cyclone is negative the force and acceleration have opposite sign. The role of dipolar structures in steadily moving vortices is discussed, and it is shown that their overall structure is fixed by the steady westward motion of the mass anomaly. Furthermore, it is shown that reduced-gravity vortices are not advected with a background flow. The reason for this behavior is that the background flow changes the ambient vorticity gradient such that the vortex obtains an extra self-propagation term that exactly cancels the advection by the background flow. Last, it is shown that a vortex initially at rest will accelerate equatorward first, after which a westward motion is generated. This result is independent of the sign of the vortex.

Flow structure and variability in the subtropical Indian Ocean: Instability of the South Indian Ocean Countercurrent

The origin of the eddy variability around the 25°S band in the Indian Ocean is investigated. We have found that the surface circulation east of Madagascar shows an anticyclonic subgyre bounded to the south by eastward flow from southwest Madagascar, and to the north by the westward flowing South Equatorial Current (SEC) between 15° and 20°S. The shallow, eastward flowing South Indian Ocean Countercurrent (SICC) extends above the deep reaching, westward flowing SEC to 95°E around the latitude of the high variability band. Applying a two-layer model reveals that regions of large vertical shear along the SICC-SEC system are baroclinically unstable. Estimates of the frequencies (3.5–6 times/year) and wavelengths (290–470 km) of the unstable modes are close to observations of the mesoscale variability derived from altimetry data. It is likely then that Rossby wave variability locally generated in the subtropical South Indian Ocean by baroclinic instability is the origin of the eddy variability around 25°S as seen, for example, in satellite altimetry.

A link between low-frequency mesoscale eddy variability around Madagascar and the large-scale Indian Ocean variability

A connection is shown to exist between the mesoscale eddy activity around Madagascar and the large-scale interannual variability in the Indian Ocean. We use the combined TOPEX/Poseidon-ERS sea surface height (SSH) data for the period 1993–2003. The SSH-fields in the Mozambique Channel and east of Madagascar exhibit a significant interannual oscillation. This is related to the arrival of large-scale anomalies that propagate westward along 10°–15°S in response to the Indian Ocean dipole (IOD) events. Positive (negative) SSH anomalies associated to a positive (negative) IOD phase induce a shift in the intensity and position of the tropical and subtropical gyres. A weakening (strengthening) results in the intensity of the South Equatorial Current and its branches along east Madagascar. In addition, the flow through the narrows of the Mozambique Channel around 17°S increases (decreases) during periods of a stronger and northward (southward) extension of the subtropical (tropical) gyre. Interaction between the currents in the narrows and southward propagating eddies from the northern Channel leads to interannual variability in the eddy kinetic energy of the central Channel in phase with the one in the SSH-field.

Observations of the inter-ocean exchange around South Africa

There is growing evidence that the interocean exchange south of Africa is an important link in the global overturning circulation of the ocean, the so‐called ocean conveyer belt. At this location, warm and salty Indian Ocean waters enter the South Atlantic and are pulled by currents that eventually reach the North Atlantic, where water cools and sinks. A major contributor to the exchange is the frequent shedding of ring eddies from the termination of the Agulhas Current south of the tip of Africa. This shedding is controlled by developments far upstream in the Indian Ocean, and variations in this ‘Agulhas Leakage’ can lead to changes in the rate and stability of the Atlantic overturning, with possible associated global climate variations [Weijer et al., 1999]. Regional climate variations in the tropical and subtropical Indian Ocean are known to affect the whole system of the Agulhas Current, including the interocean exchanges. This article reports on some of the seminal results of ongoing multinational, multidisciplinary projects that explore these issues.

Data assimilation for marine monitoring and prediction: The MERCATOR operational assimilation systems and the MERSEA developments

During the past 15 years, a number of initiatives have been undertaken at national level to develop ocean forecasting systems operating at regional and/or global scales. The co-ordination between these efforts has been organized internationally through the Global Ocean Data Assimilation Experiment (GODAE). The French MERCATOR project is one of the leading participants in GODAE. The MERCATOR systems routinely assimilate a variety of observations such as multi-satellite altimeter data, sea-surface temperature and in situ temperature and salinity profiles, focusing on high-resolution scales of the ocean dynamics. The assimilation strategy in MERCATOR is based on a hierarchy of methods of increasing sophistication including optimal interpolation, Kalman filtering and variational methods, which are progressively deployed through the Syst`eme d’Assimilation MERCATOR (SAM) series. SAM-1 is based on a reduced-order optimal interpolation which can be operated using ‘altimetry-only’ or ‘multi-data’ set-ups; it relies on the concept of separability, assuming that the correlations can be separated into a product of horizontal and vertical contributions. The second release, SAM-2, is being developed to include new features from the singular evolutive extended Kalman (SEEK) filter, such as three-dimensional, multivariate error modes and adaptivity schemes. The third one, SAM-3, considers variational methods such as the incremental four-dimensional variational algorithm. Most operational forecasting systems evaluated during GODAE are based on least-squares statistical estimation assuming Gaussian errors. In the framework of the EU MERSEA (Marine EnviRonment and Security for the European Area) project, research is being conducted to prepare the next-generation operational ocean monitoring and forecasting systems. The research effort will explore nonlinear assimilation formulations to overcome limitations of the current systems. This paper provides an overview of the developments conducted in MERSEA with the SEEK filter, the Ensemble Kalman filter and the sequential importance re-sampling filter.

Atmospheric composition of West Africa: highlights from the AMMA program

The atmospheric composition of West Africa reflects the interaction of various dynamical and chemical systems (i.e. biogenic, urban, convective and long-range transport) with signatures from local to continental scales. Recent measurements performed during the African Monsoon Multidisciplinary Analyses (AMMA) observational periods in 2005 and 2006 provide new data which has allowed new insight into the processes within these systems that control the distribution of ozone and its precursors. Using these new data and recently published results, we provide an overview of these systems with a particular emphasis on ozone distributions over West Africa during the wet season.

Radiative forcing due to changes in ozone and methane caused by the transport sector

The year 2000 radiative forcing (RF) due to changes in O3 and CH4 (and the CH4-induced stratospheric water vapour) as a result of emissions of short-lived gases (oxides of nitrogen (NOx), carbon monoxide and non-methane hydrocarbons) from three transport sectors (ROAD, maritime SHIPping and AIRcraft) are calculated using results from five global atmospheric chemistry models. Using results from these models plus other published data, we quantify the uncertainties. The RF due to short-term O3 changes (i.e. as an immediate response to the emissions without allowing for the long-term CH4 changes) is positive and highest for ROAD transport (31mWm-2) compared to SHIP (24 mWm-2) and AIR (17 mWm-2) sectors in four of the models. All five models calculate negative RF from the CH4 perturbations, with a larger impact from the SHIP sector than for ROAD and AIR. The net RF of O3 and CH4 combined (i.e. including the impact of CH4 on ozone and stratospheric water vapour) is positive for ROAD (+16(±13)(one standard deviation) mWm-2) and AIR (+6(±5) mWm-2) traffic sectors and is negative for SHIP (-18(±10) mWm-2) sector in all five models. Global Warming Potentials (GWP) and Global Temperature change Potentials (GTP) are presented for AIR NOx emissions; there is a wide spread in the results from the 5 chemistry models, and it is shown that differences in the methane response relative to the O3 response drive much of the spread.

Comparing correlations of continuous observations from two independent populations using a sequential approach

A sequential study design generally makes more efficient use of available information than a fixed sample counterpart of equal power. This feature is gradually being exploited by researchers in genetic and epidemiological investigations that utilize banked biological resources and in studies where time, cost and ethics are prominent considerations. Recent work in this area has focussed on the sequential analysis of matched case-control studies with a dichotomous trait. In this paper, we extend the sequential approach to a comparison of the associations within two independent groups of paired continuous observations. Such a comparison is particularly relevant in familial studies of phenotypic correlation using twins. We develop a sequential twin method based on the intraclass correlation and show that use of sequential methodology can lead to a substantial reduction in the number of observations without compromising the study error rates. Additionally, our approach permits straightforward allowance for other explanatory factors in the analysis. We illustrate our method in a sequential heritability study of dysplasia that allows for the effect of body mass index and compares monozygotes with pairs of singleton sisters. Copyright (c) 2006 John Wiley & Sons, Ltd.

Self-assembled biomimetic [2Fe2S]-hydrogenase-based photocatalyst for molecular hydrogen evolution

The large-scale production of clean energy is one of the major challenges society is currently facing. Molecular hydrogen is envisaged as a key green fuel for the future, but it becomes a sustainable alternative for classical fuels only if it is also produced in a clean fashion. Here, we report a supramolecular biomimetic approach to form a catalyst that produces molecular hydrogen using light as the energy source. It is composed of an assembly of chromophores to a bis(thiolate)-bridged diiron ([2Fe2S]) based hydrogenase catalyst. The supramolecular building block approach introduced in this article enabled the easy formation of a series of complexes, which are all thoroughly characterized, revealing that the photoactivity of the catalyst assembly strongly depends on its nature. The active species, formed from different complexes, appears to be the [Fe-2(mu-pdt)(CO)(4){PPh2(4-py)}(2)] (3) with 2 different types of porphyrins (5a and 5b) coordinated to it. The modular supramolecular approach was important in this study as with a limited number of building blocks several different complexes were generated.

Effects of utterance length on lip kinematics in aphasia

Most existing models of language production and speech motor control do not explicitly address how language requirements affect speech motor functions, as these domains are usually treated as separate and independent from one another. This investigation compared lip movements during bilabial closure between five individuals with mild aphasia and five age and gender-matched control speakers when the linguistic characteristics of the stimuli were varied by increasing the number of syllables. Upper and lower lip movement data were collected for mono-, bi- and tri-syllabic nonword sequences using an AG 100 EMMA system. Each task was performed under both normal and fast rate conditions. Single articulator kinematic parameters (peak velocity, amplitude, duration,and cyclic spatio-temporal index) were measured to characterize lip movements. Results revealed that compared to control speakers, individuals with aphasia showed significantly longer movement duration and lower movement stability for longer items (bi- and tri-syllables). Moreover, utterance length affected the lip kinematics, in that the monosyllables had smaller peak velocities, smaller amplitudes and shorter durations compared to bi- and trisyllables, and movement stability was lowest for the trisyllables. In addition, the rate-induced changes (smaller amplitude and shorter duration with increased rate) were most prominent for the short items (i.e., monosyllables). These findings provide further support for the notion that linguistic changes have an impact on the characteristics of speech movements, and that individuals with aphasia are more affected by such changes than control speakers.

Word frequency and bigram frequency effects on linguistic processing and speech motor performance in individuals with aphasia and normal speakers

Models of normal word production are well specified about the effects of frequency of linguistic stimuli on lexical access, but are less clear regarding the same effects on later stages of word production, particularly word articulation. In aphasia, this lack of specificity of down-stream frequency effects is even more noticeable because there is relatively limited amount of data on the time course of frequency effects for this population. This study begins to fill this gap by comparing the effects of variation of word frequency (lexical, whole word) and bigram frequency (sub-lexical, within word) on word production abilities in ten normal speakers and eight mild–moderate individuals with aphasia. In an immediate repetition paradigm, participants repeated single monosyllabic words in which word frequency (high or low) was crossed with bigram frequency (high or low). Indices for mapping the time course for these effects included reaction time (RT) for linguistic processing and motor preparation, and word duration (WD) for speech motor performance (word articulation time). The results indicated that individuals with aphasia had significantly longer RT and WD compared to normal speakers. RT showed a significant main effect only for word frequency (i.e., high-frequency words had shorter RT). WD showed significant main effects of word and bigram frequency; however, contrary to our expectations, high-frequency items had longer WD. Further investigation of WD revealed that independent of the influence of word and bigram frequency, vowel type (tense or lax) had the expected effect on WD. Moreover, individuals with aphasia differed from control speakers in their ability to implement tense vowel duration, even though they could produce an appropriate distinction between tense and lax vowels. The results highlight the importance of using temporal measures to identify subtle deficits in linguistic and speech motor processing in aphasia, the crucial role of phonetic characteristics of stimuli set in studying speech production and the need for the language production models to account more explicitly for word articulation.

Speech motor control in fluent and dysfluent speech production of an individual with apraxia of speech and Broca’s aphasia

Apraxia of speech (AOS) is typically described as a motor-speech disorder with clinically well-defined symptoms, but without a clear understanding of the underlying problems in motor control. A number of studies have compared the speech of subjects with AOS to the fluent speech of controls, but only a few have included speech movement data and if so, this was primarily restricted to the study of single articulators. If AOS reflects a basic neuromotor dysfunction, this should somehow be evident in the production of both dysfluent and perceptually fluent speech. The current study compared motor control strategies for the production of perceptually fluent speech between a young woman with apraxia of speech (AOS) and Broca’s aphasia and a group of age-matched control speakers using concepts and tools from articulation-based theories. In addition, to examine the potential role of specific movement variables on gestural coordination, a second part of this study involved a comparison of fluent and dysfluent speech samples from the speaker with AOS. Movement data from the lips, jaw and tongue were acquired using the AG-100 EMMA system during the reiterated production of multisyllabic nonwords. The findings indicated that although in general kinematic parameters of fluent speech were similar in the subject with AOS and Broca’s aphasia to those of the age-matched controls, speech task-related differences were observed in upper lip movements and lip coordination. The comparison between fluent and dysfluent speech characteristics suggested that fluent speech was achieved through the use of specific motor control strategies, highlighting the potential association between the stability of coordinative patterns and movement range, as described in Coordination Dynamics theory.

Nonlinear data assimilation in geosciences: an extremely efficient particle filter

Almost all research fields in geosciences use numerical models and observations and combine these using data-assimilation techniques. With ever-increasing resolution and complexity, the numerical models tend to be highly nonlinear and also observations become more complicated and their relation to the models more nonlinear. Standard data-assimilation techniques like (ensemble) Kalman filters and variational methods like 4D-Var rely on linearizations and are likely to fail in one way or another. Nonlinear data-assimilation techniques are available, but are only efficient for small-dimensional problems, hampered by the so-called ‘curse of dimensionality’. Here we present a fully nonlinear particle filter that can be applied to higher dimensional problems by exploiting the freedom of the proposal density inherent in particle filtering. The method is illustrated for the three-dimensional Lorenz model using three particles and the much more complex 40-dimensional Lorenz model using 20 particles. By also applying the method to the 1000-dimensional Lorenz model, again using only 20 particles, we demonstrate the strong scale-invariance of the method, leading to the optimistic conjecture that the method is applicable to realistic geophysical problems. Copyright c 2010 Royal Meteorological Society

Efficient non-linear data assimilation in geophysical fluid dynamics

New ways of combining observations with numerical models are discussed in which the size of the state space can be very large, and the model can be highly nonlinear. Also the observations of the system can be related to the model variables in highly nonlinear ways, making this data-assimilation (or inverse) problem highly nonlinear. First we discuss the connection between data assimilation and inverse problems, including regularization. We explore the choice of proposal density in a Particle Filter and show how the ’curse of dimensionality’ might be beaten. In the standard Particle Filter ensembles of model runs are propagated forward in time until observations are encountered, rendering it a pure Monte-Carlo method. In large-dimensional systems this is very inefficient and very large numbers of model runs are needed to solve the data-assimilation problem realistically. In our approach we steer all model runs towards the observations resulting in a much more efficient method. By further ’ensuring almost equal weight’ we avoid performing model runs that are useless in the end. Results are shown for the 40 and 1000 dimensional Lorenz 1995 model.

Measuring the impact of observations on the predictability of the Kuroshio Extension in a shallow-water model

In this paper sequential importance sampling is used to assess the impact of observations on a ensemble prediction for the decadal path transitions of the Kuroshio Extension (KE). This particle filtering approach gives access to the probability density of the state vector, which allows us to determine the predictive power — an entropy based measure — of the ensemble prediction. The proposed set-up makes use of an ensemble that, at each time, samples the climatological probability distribution. Then, in a post-processing step, the impact of different sets of observations is measured by the increase in predictive power of the ensemble over the climatological signal during one-year. The method is applied in an identical-twin experiment for the Kuroshio Extension using a reduced-gravity shallow water model. We investigate the impact of assimilating velocity observations from different locations during the elongated and the contracted meandering state of the KE. Optimal observations location correspond to regions with strong potential vorticity gradients. For the elongated state the optimal location is in the first meander of the KE. During the contracted state of the KE it is located south of Japan, where the Kuroshio separates from the coast.