Nonlinear saturation of baroclinic instability: part I: the two-layer model

A rigorous bound is derived which limits the finite-amplitude growth of arbitrary nonzonal disturbances to an unstable baroclinic zonal flow within the context of the two-layer model. The bound is valid for conservative (unforced) flow, as well as for forced-dissipative flow that when the dissipation is proportional to the potential vorticity. The method used to derive the bound relies on the existence of a nonlinear Liapunov (normed) stability theorem for subcritical flows, which is a finite-amplitude generalization of the Charney-Stern theorem. For the special case of the Philips model of baroclinic instability, and in the limit of infinitesimal initial nonzonal disturbance amplitude, an improved form of the bound is possible which states that the potential enstrophy of the nonzonal flow cannot exceed ϵβ2, where ϵ = (U − Ucrit)/Ucrit is the (relative) supereriticality. This upper bound turns out to be extremely similar to the maximum predicted by the weakly nonlinear theory. For unforced flow with ϵ < 1, the bound demonstrates that the nonzonal flow cannot contain all of the potential enstrophy in the system; hence in this range of initial supercriticality the total flow must remain, in a certain sense, “close” to a zonal state.

Unpicking antecedents of CRM adoption: a two-stage model

Most CRM work focuses on consumer applications. This paper addresses the operational adoption issues facing the organisation deploying CRM practices. There are a plethora of challenges facing organisations when adopting CRM. Previous research is limited to either examining the CRM adoption process at an individual/employees level or an organisational level. Hence, in this paper the myriad of organisational, marketing and technical antecedents that seem to impinge upon employee perceptions and organisational implementation of CRM are structured in a two-stage model. Using a stratified sample of ten organisations across four sectors, seven hypotheses are tested on data collected from 301 practitioners. A two-stage model is analysed using structural equation modelling. Findings reveal that CRM implementation relates to employee perceptions of CRM. This paper deepens our understanding of organisational practices to adopt CRM, so as an organisation properly profits from the expected benefits of CRM.

QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments.

Nonlinear saturation of baroclinic instability: part II: continuously stratified fluid

Rigorous upper bounds are derived that limit the finite-amplitude growth of arbitrary nonzonal disturbances to an unstable baroclinic zonal flow in a continuously stratified, quasi-geostrophic, semi-infinite fluid. Bounds are obtained bath on the depth-integrated eddy potential enstrophy and on the eddy available potential energy (APE) at the ground. The method used to derive the bounds is essentially analogous to that used in Part I of this study for the two-layer model: it relies on the existence of a nonlinear Liapunov (normed) stability theorem, which is a finite-amplitude generalization of the Charney-Stern theorem. As in Part I, the bounds are valid both for conservative (unforced, inviscid) flow, as well as for forced-dissipative flow when the dissipation is proportional to the potential vorticity in the interior, and to the potential temperature at the ground. The character of the results depends on the dimensionless external parameter γ = f02ξ/β0N2H, where ξ is the maximum vertical shear of the zonal wind, H is the density scale height, and the other symbols have their usual meaning. When γ ≫ 1, corresponding to “deep” unstable modes (vertical scale ≈H), the bound on the eddy potential enstrophy is just the total potential enstrophy in the system; but when γ≪1, corresponding to ‘shallow’ unstable modes (vertical scale ≈γH), the eddy potential enstrophy can be bounded well below the total amount available in the system. In neither case can the bound on the eddy APE prevent a complete neutralization of the surface temperature gradient which is in accord with numerical experience. For the special case of the Charney model of baroclinic instability, and in the limit of infinitesimal initial eddy disturbance amplitude, the bound states that the dimensionless eddy potential enstrophy cannot exceed (γ + 1)2/24&gamma2h when γ ≥ 1, or 1/6;&gammah when γ ≤ 1; here h = HN/f0L is the dimensionless scale height and L is the width of the channel. These bounds are very similar to (though of course generally larger than) ad hoc estimates based on baroclinic-adjustment arguments. The possibility of using these kinds of bounds for eddy-amplitude closure in a transient-eddy parameterization scheme is also discussed.

Evaluation of satellite-based and model re-analysis rainfall estimates for Uganda

The dependence of much of Africa on rain fed agriculture leads to a high vulnerability to fluctuations in rainfall amount. Hence, accurate monitoring of near-real time rainfall is particularly useful, for example in forewarning possible crop shortfalls in drought-prone areas. Unfortunately, ground based observations are often inadequate. Rainfall estimates from satellite-based algorithms and numerical model outputs can fill this data gap, however rigorous assessment of such estimates is required. In this case, three satellite based products (NOAA-RFE 2.0, GPCP-1DD and TAMSAT) and two numerical model outputs (ERA-40 and ERA-Interim) have been evaluated for Uganda in East Africa using a network of 27 rain gauges. The study focuses on the years 2001 to 2005 and considers the main rainy season (February to June). All data sets were converted to the same temporal and spatial scales. Kriging was used for the spatial interpolation of the gauge data. All three satellite products showed similar characteristics and had a high level of skill that exceeded both model outputs. ERA-Interim had a tendency to overestimate whilst ERA-40 consistently underestimated the Ugandan rainfall.

Some numerical experiments on the effect of the variation of static stability in two-layer quasi-geostrophic models

A method to solve a quasi-geostrophic two-layer model including the variation of static stability is presented. The divergent part of the wind is incorporated by means of an iterative procedure. The procedure is rather fast and the time of computation is only 60–70% longer than for the usual two-layer model. The method of solution is justified by the conservation of the difference between the gross static stability and the kinetic energy. To eliminate the side-boundary conditions the experiments have been performed on a zonal channel model. The investigation falls mainly into three parts: The first part (section 5) contains a discussion of the significance of some physically inconsistent approximations. It is shown that physical inconsistencies are rather serious and for these inconsistent models which were studied the total kinetic energy increased faster than the gross static stability. In the next part (section 6) we are studying the effect of a Jacobian difference operator which conserves the total kinetic energy. The use of this operator in two-layer models will give a slight improvement but probably does not have any practical use in short periodic forecasts. It is also shown that the energy-conservative operator will change the wave-speed in an erroneous way if the wave-number or the grid-length is large in the meridional direction. In the final part (section 7) we investigate the behaviour of baroclinic waves for some different initial states and for two energy-consistent models, one with constant and one with variable static stability. According to the linear theory the waves adjust rather rapidly in such a way that the temperature wave will lag behind the pressure wave independent of the initial configuration. Thus, both models give rise to a baroclinic development even if the initial state is quasi-barotropic. The effect of the variation of static stability is very small, qualitative differences in the development are only observed during the first 12 hours. For an amplifying wave we will get a stabilization over the troughs and an instabilization over the ridges.

Mean motions induced by baroclinic instability in a jet

A study is made of the zonal-mean motions induced by a growing baroclinic wave in several contexts, under the framework of three different analysis schemes: the conventional Eulerian mean (EM), the transformed Eulerian mean (TEM), and the generalized Lagrangian mean (GLM). The effect of meridional shear in the initial jet on these induced mean motions is considered by treating the instability problem in the context of the two-layer model. The conceptual simplicity of the TEM formulation is shown to be useful in diagnosing the dynamics of instability, much as it has been found helpful in many problems of wave, mean-flow interaction. In addition, it is found that the TEM vertical velocity is a very good indicator of the GLM vertical velocity. However, the GLM meridional velocity is always convergent towards the centre of instability activity, and is not at all well represented by the nondivergent TEM meridional velocity. In comparing the results with Uryu's (1979) calculation of the GLM circulation induced by a growing Eady wave, it is found that the inclusion of meridional jet shear in the present work leads to some strikingly different effects in the GLM zonal wind acceleration. In the case of pure baroclinic instability treated by Uryu, the Eulerian and Stokes accelerations nearly cancel each other in the centre of the channel, leaving a weak Lagrangian acceleration opposed to the Eulerian one. In the more general case of mixed baroclinic-barotropic instability, however, the Eulerian and Stokes accelerations can reinforce one another, leading to a very strong Lagrangian zonal wind

The coevolution of two phytoplankton species on a single resource: allelopathy as a pseudo-mixotrophy

Without the top-down effects and the external/physical forcing, a stable coexistence of two phytoplankton species under a single resource is impossible — a result well known from the principle of competitive exclusion. Here I demonstrate by analysis of a mathematical model that such a stable coexistence in a homogeneous media without any external factor would be possible, at least theoretically, provided (i) one of the two species is toxin producing thereby has an allelopathic effect on the other, and (ii) the allelopathic effect exceeds a critical level. The threshold level of allelopathy required for the coexistence has been derived analytically in terms of the parameters associated with the resource competition and the nutrient recycling. That the extra mortality of a competitor driven by allelopathy of a toxic species gives a positive feed back to the algal growth process through the recycling is explained. And that this positive feed back plays a pivotal role in reducing competition pressures and helping species succession in the two-species model is demonstrated. Based on these specific coexistence results, I introduce and explain theoretically the allelopathic effect of a toxic species as a ‘pseudo-mixotrophy’—a mechanism of ‘if you cannot beat them or eat them, just kill them by chemical weapons’. The impact of this mechanism of species succession by pseudo-mixotrophy in the form of alleopathy is discussed in the context of current understanding on straight mixotrophy and resource-species relationship among phytoplankton species.

A statistical model for sea surface diurnal warming driven by numercial weather predictions fluxes and winds

A statistical model is derived relating the diurnal variation of sea surface temperature (SST) to the net surface heat flux and surface wind speed from a numerical weather prediction (NWP) model. The model is derived using fluxes and winds from the European Centre for Medium-Range Weather Forecasting (ECMWF) NWP model and SSTs from the Spinning Enhanced Visible and Infrared Imager (SEVIRI). In the model, diurnal warming has a linear dependence on the net surface heat flux integrated since (approximately) dawn and an inverse quadratic dependence on the maximum of the surface wind speed in the same period. The model coefficients are found by matching, for a given integrated heat flux, the frequency distributions of the maximum wind speed and the observed warming. Diurnal cooling, where it occurs, is modelled as proportional to the integrated heat flux divided by the heat capacity of the seasonal mixed layer. The model reproduces the statistics (mean, standard deviation, and 95-percentile) of the diurnal variation of SST seen by SEVIRI and reproduces the geographical pattern of mean warming seen by the Advanced Microwave Scanning Radiometer (AMSR-E). We use the functional dependencies in the statistical model to test the behaviour of two physical model of diurnal warming that display contrasting systematic errors.

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.

Ranking efficiency units in DEA using bootstrapping an applied analysis for Slovenian farm data

This article explores how data envelopment analysis (DEA), along with a smoothed bootstrap method, can be used in applied analysis to obtain more reliable efficiency rankings for farms. The main focus is the smoothed homogeneous bootstrap procedure introduced by Simar and Wilson (1998) to implement statistical inference for the original efficiency point estimates. Two main model specifications, constant and variable returns to scale, are investigated along with various choices regarding data aggregation. The coefficient of separation (CoS), a statistic that indicates the degree of statistical differentiation within the sample, is used to demonstrate the findings. The CoS suggests a substantive dependency of the results on the methodology and assumptions employed. Accordingly, some observations are made on how to conduct DEA in order to get more reliable efficiency rankings, depending on the purpose for which they are to be used. In addition, attention is drawn to the ability of the SLICE MODEL, implemented in GAMS, to enable researchers to overcome the computational burdens of conducting DEA (with bootstrapping).

A DFT study of free radicals formed from artemisinin and related compounds

Ab initio calculations using density functional theory have shown that the reactions that occur between artemisinin, 1, a cyclic trioxane active against malaria, and some metal ions and complexes lead to a series of radicals which are probably responsible for its therapeutic activity. In particular it has been shown that the interaction of Fe(H) with artemisinin causes the O-O bond to be broken as indeed does Fe(III) and Cu(I), while Zn(II) does not. Calculations were carried out with Fe(II) in several different forms including the bare ion, [Fe(H2O)(5)](2+) and [FeP(Im)] (P, porphyrin; Im, imadazole) and similar results were obtained. The resulting oxygen-based radicals are readily converted to more stable carbon-based radicals and/or. stable products. Similar radicals and products are also formed from two simple model trioxanes 2 and 3 that show little or no therapeutic action against malaria although some subtle differences were obtained. This suggests that the scaffold surrounding the pharmacophore may be involved in molecular recognition events allowing efficient uptake of this trioxane warhead into the parasite. (C) 2004 Elsevier B.V. All rights reserved.

Structural analysis of PEO-PBO copolymer monolayers at the air-water interface

X-ray reflectivity (XR) and grazing incidence X-ray diffraction (GIXD) have been used to examine an oxyethylene-b-oxybutylene (E23B8) copolymer film at the air-water interface. The XR data were fitted using both a one- and a two-layer model that outputted the film thickness, roughness, and electron density. The best fit to the experimental data was obtained using a two-layer model (representing the oxyethylene and oxybutylene blocks, respectively), which showed a rapid thickening of the copolymer film at pressures above 7 mN/m. The large roughness values found indicate a significant degree of intermixing between the blocks and back up the GIXD data, which showed no long range lateral ordering within the layer. It was found from the electron density model results that there is a large film densification at 7 mN/m, possibly suggesting conformational changes within the film, even though no such change occurs on the pressure-area isotherm at the same surface pressure.

Comparative characterisation by atomic force microscopy and ellipsometry of soft and solid thin films

Ellipsometry and atomic force microscopy (AFM) were used to study the film thickness and the surface roughness of both 'soft' and solid thin films. 'Soft' polymer thin films of polystyrene and poly(styrene-ethylene/butylene-styrene) block copolymer were prepared by spin-coating onto planar silicon wafers. Ellipsometric parameters were fitted by the Cauchy approach using a two-layer model with planar boundaries between the layers. The smooth surfaces of the prepared polymer films were confirmed by AFM. There is good agreement between AFM and ellipsometry in the 80-130 nm thickness range. Semiconductor surfaces (Si) obtained by anisotropic chemical etching were investigated as an example of a randomly rough surface. To define roughness parameters by ellipsometry, the top rough layers were treated as thin films according to the Bruggeman effective medium approximation (BEMA). Surface roughness values measured by AFM and ellipsometry show the same tendency of increasing roughness with increased etching time, although AFM results depend on the used window size. The combined use of both methods appears to offer the most comprehensive route to quantitative surface roughness characterisation of solid films. Copyright (c) 2007 John Wiley & Sons, Ltd.

Selfishness and altruism can coexist when help is subject to diminishing returns

Altruism and selfishness are 30–50% heritable in man in both Western and non-Western populations. This genetically based variation in altruism and selfishness requires explanation. In non-human animals, altruism is generally directed towards relatives, and satisfies the condition known as Hamilton's rule. This nepotistic altruism evolves under natural selection only if the ratio of the benefit of receiving help to the cost of giving it exceeds a value that depends on the relatedness of the individuals involved. Standard analyses assume that the benefit provided by each individual is the same but it is plausible in some cases that as more individuals contribute, help is subject to diminishing returns. We analyse this situation using a single-locus two-allele model of selection in a diploid population with the altruistic allele dominant to the selfish allele. The analysis requires calculation of the relationship between the fitnesses of the genotypes and the frequencies of the genes. The fitnesses vary not only with the genotype of the individual but also with the distribution of phenotypes amongst the sibs of the individual and this depends on the genotypes of his parents. These calculations are not possible by direct fitness or ESS methods but are possible using population genetics. Our analysis shows that diminishing returns change the operation of natural selection and the outcome can now be a stable equilibrium between altruistic and selfish alleles rather than the elimination of one allele or the other. We thus provide a plausible genetic model of kin selection that leads to the stable coexistence in the same population of both altruistic and selfish individuals. This may explain reported genetic variation in altruism in man.