Balanced Ocean-Data Assimilation near the Equator

The question is addressed whether using unbalanced updates in ocean-data assimilation schemes for seasonal forecasting systems can result in a relatively poor simulation of zonal currents. An assimilation scheme, where temperature observations are used for updating only the density field, is compared to a scheme where updates of density field and zonal velocities are related by geostrophic balance. This is done for an equatorial linear shallow-water model. It is found that equatorial zonal velocities can be detoriated if velocity is not updated in the assimilation procedure. Adding balanced updates to the zonal velocity is shown to be a simple remedy for the shallow-water model. Next, optimal interpolation (OI) schemes with balanced updates of the zonal velocity are implemented in two ocean general circulation models. First tests indicate a beneficial impact on equatorial upper-ocean zonal currents.

A regime-dependent balanced control variable based on potential vorticity

In this paper it is argued that rotational wind is not the best choice of leading control variable for variational data assimilation, and an alternative is suggested and tested. A rotational wind parameter is used in most global variational assimilation systems as a pragmatic way of approximately representing the balanced component of the assimilation increments. In effect, rotational wind is treated as a proxy for potential vorticity, but one that it is potentially not a good choice in flow regimes characterised by small Burger number. This paper reports on an alternative set of control variables which are based around potential vorticity. This gives rise to a new formulation of the background error covariances for the Met Office's variational assimilation system, which leads to flow dependency. It uses similar balance relationships to traditional schemes, but recognises the existence of unbalanced rotational wind which is used with a new anti-balance relationship. The new scheme is described and its performance is evaluated and compared to a traditional scheme using a sample of diagnostics.

Comment on ``Piecewise potential vorticity inversion: Elementary tests''

Egger (2008) constructs some idealised experiments to test the usefulness of piecewise potential vorticity inversion (PPVI) in the diagnosis of Rossby wave dynamics and baroclinic development. He concludes that, ``PPVI does not help us to understand the dynamics of linear Rossby waves. It provides local tendencies of the streamfunction which are unrelated to the true ones. The same way, the motion of baroclinic waves in shear flow cannot be understood by using PPVI. Moreover, the effect of boundary temperatures as determined by PPVI is unrelated to the flow evolution.'' He goes further in arguing that we should not consider velocities as ``induced'' by PV anomalies defined by carving up the global domain. However, these conclusions partly reflect the limitations of his idealised experiments and the manner in which the PV components were partitioned from one another.

A high-order arbitrarily unstructured finite-volume model of the global atmosphere: tests solving the shallow-water equations

Simulations of the global atmosphere for weather and climate forecasting require fast and accurate solutions and so operational models use high-order finite differences on regular structured grids. This precludes the use of local refinement; techniques allowing local refinement are either expensive (eg. high-order finite element techniques) or have reduced accuracy at changes in resolution (eg. unstructured finite-volume with linear differencing). We present solutions of the shallow-water equations for westerly flow over a mid-latitude mountain from a finite-volume model written using OpenFOAM. A second/third-order accurate differencing scheme is applied on arbitrarily unstructured meshes made up of various shapes and refinement patterns. The results are as accurate as equivalent resolution spectral methods. Using lower order differencing reduces accuracy at a refinement pattern which allows errors from refinement of the mountain to accumulate and reduces the global accuracy over a 15 day simulation. We have therefore introduced a scheme which fits a 2D cubic polynomial approximately on a stencil around each cell. Using this scheme means that refinement of the mountain improves the accuracy after a 15 day simulation. This is a more severe test of local mesh refinement for global simulations than has been presented but a realistic test if these techniques are to be used operationally. These efficient, high-order schemes may make it possible for local mesh refinement to be used by weather and climate forecast models.

Inertia-Gravity Waves Emitted from Balanced Flow: Observations, Properties, and Consequences

This paper describes laboratory observations of inertia–gravity waves emitted from balanced fluid flow. In a rotating two-layer annulus experiment, the wavelength of the inertia–gravity waves is very close to the deformation radius. Their amplitude varies linearly with Rossby number in the range 0.05–0.14, at constant Burger number (or rotational Froude number). This linear scaling challenges the notion, suggested by several dynamical theories, that inertia–gravity waves generated by balanced motion will be exponentially small. It is estimated that the balanced flow leaks roughly 1% of its energy each rotation period into the inertia–gravity waves at the peak of their generation. The findings of this study imply an inevitable emission of inertia–gravity waves at Rossby numbers similar to those of the large-scale atmospheric and oceanic flow. Extrapolation of the results suggests that inertia–gravity waves might make a significant contribution to the energy budgets of the atmosphere and ocean. In particular, emission of inertia–gravity waves from mesoscale eddies may be an important source of energy for deep interior mixing in the ocean.

The influence of different artificial soil types on the acute toxicity of carbendazim to the earthworm Eisenia fetida in laboratory toxicity tests

Field populations of earthworms have shown a varied response in mortality to the fungicide carbendazim, the toxic reference substance used in agrochemical field trials. The aim of this study was to determine the influence of soil conditions as a potential cause of this variation. Laboratory acute toxicity tests were conducted using a range of artificial soils with varying soil components (organic matter, clay, pH and moisture). Batch adsorption/desorption studies were run to determine the influence of the soil properties on carbendazim behaviour. Adsorption was shown to be correlated with organic matter content and pH and this in turn could be linked to Eisenia fetida mortality, with lower mortality occurring with increased adsorption. Overall while E.fetida mortality did vary significantly between several of the soils the calculated LC50 values in the different soils did not cover a wide range (6.04-16.00 mg kg(-1)), showing that under these laboratory conditions soil components did not greatly influence carbendazim toxicity to E.fietida. (c) 2007 Elsevier Masson SAS. All rights reserved.

A capture–recapture approach for screening using two diagnostic tests with availability of disease status for the test positives only

The article considers screening human populations with two screening tests. If any of the two tests is positive, then full evaluation of the disease status is undertaken; however, if both diagnostic tests are negative, then disease status remains unknown. This procedure leads to a data constellation in which, for each disease status, the 2 × 2 table associated with the two diagnostic tests used in screening has exactly one empty, unknown cell. To estimate the unobserved cell counts, previous approaches assume independence of the two diagnostic tests and use specific models, including the special mixture model of Walter or unconstrained capture–recapture estimates. Often, as is also demonstrated in this article by means of a simple test, the independence of the two screening tests is not supported by the data. Two new estimators are suggested that allow associations of the screening test, although the form of association must be assumed to be homogeneous over disease status. These estimators are modifications of the simple capture–recapture estimator and easy to construct. The estimators are investigated for several screening studies with fully evaluated disease status in which the superior behavior of the new estimators compared to the previous conventional ones can be shown. Finally, the performance of the new estimators is compared with maximum likelihood estimators, which are more difficult to obtain in these models. The results indicate the loss of efficiency as minor.

The effect of consumption volume on profile and liking of oral nutritional supplements of varied sweetness: sequential profiling and boredom tests

Oral nutrition supplements (ONS) are routinely prescribed to those with, or at risk of, malnutrition. Previous research identified poor compliance due to taste and sweetness. This paper investigates taste and hedonic liking of ONS, of varying sweetness and metallic levels, over consumption volume; an important consideration as patients are prescribed large volumes of ONS daily. A sequential descriptive profile was developed to determine the perception of sensory attributes over repeat consumption of ONS. Changes in liking of ONS following repeat consumption were characterised by a boredom test. Certain flavour (metallic taste, soya milk flavour) and mouthfeel (mouthdrying, mouthcoating) attributes built up over increased consumption volume (p 0.002). Hedonic liking data from two cohorts, healthy older volunteers (n = 32, median age 73) and patients (n = 28, median age 85), suggested such build-up was disliked. Efforts made to improve the palatability of ONS must take account of the build up of taste and mouthfeel characteristics over increased consumption volume.

Phyto-activity of biocides used to manipulate herbivory: tests of three pesticides on fourteen plant species

1. We tested three pesticides used for field manipulations of herbivory for direct phytoactive effects on the germination and growth of 14 herbaceous plant species selected to provide a range of life-history strategies and functional groups. 2. We report three companion experiments: (A) Two insecticides, chlorpyrifos (granular soil insecticide) and dimethoate (foliar spray), were applied in fully-factorial combination to pot-germinated individuals of 12 species. (B) The same fully-factorial design was used to test for direct effects on the germination of four herbaceous legumes. (C) The molluscicide, metaldehyde, was tested for direct effects on the germination and growth of six plant species. 3. The insecticides had few significant effects on growth and germination. Dimethoate acted only on growth stimulating Anisantha sterilis, Sonchus asper and Stellaria graminea. In contrast, chlorpyrifos acted on germination increasing the germination of Trifolium dubium and Trifolium pratense. There was also a significant interactive effect of chlorpyrifos and dimethoate on the germination of T pratense. However, all. effects were relatively small in magnitude and explanatory power. The molluscicide had no significant effect on plant germination or growth. 4. The small number and size of direct effects of the pesticides on plant performance is encouraging for the use of these pesticides in manipulative experiments on herbivory, especially for the molluscicide. However, a smatt number of direct (positive) effects of the insecticides on some plant species need to be taken into account when interpreting field manipulations of herbivory with these compounds, and emphasises the importance of conducting tests for direct phyto-active effects. (C) 2004 Elsevier GmbH. All rights reserved.

A capture-recapture approach for screening using two diagnostic tests with availability of disease status for the test positives only

The article considers screening human populations with two screening tests. If any of the two tests is positive, then full evaluation of the disease status is undertaken; however, if both diagnostic tests are negative, then disease status remains unknown. This procedure leads to a data constellation in which, for each disease status, the 2 x 2 table associated with the two diagnostic tests used in screening has exactly one empty, unknown cell. To estimate the unobserved cell counts, previous approaches assume independence of the two diagnostic tests and use specific models, including the special mixture model of Walter or unconstrained capture-recapture estimates. Often, as is also demonstrated in this article by means of a simple test, the independence of the two screening tests is not supported by the data. Two new estimators are suggested that allow associations of the screening test, although the form of association must be assumed to be homogeneous over disease status. These estimators are modifications of the simple capture-recapture estimator and easy to construct. The estimators are investigated for several screening studies with fully evaluated disease status in which the superior behavior of the new estimators compared to the previous conventional ones can be shown. Finally, the performance of the new estimators is compared with maximum likelihood estimators, which are more difficult to obtain in these models. The results indicate the loss of efficiency as minor.

Detection of frailty in Weibull lifetime data using Outlier tests

Heterogeneity in lifetime data may be modelled by multiplying an individual's hazard by an unobserved frailty. We test for the presence of frailty of this kind in univariate and bivariate data with Weibull distributed lifetimes, using statistics based on the ordered Cox-Snell residuals from the null model of no frailty. The form of the statistics is suggested by outlier testing in the gamma distribution. We find through simulation that the sum of the k largest or k smallest order statistics, for suitably chosen k , provides a powerful test when the frailty distribution is assumed to be gamma or positive stable, respectively. We provide recommended values of k for sample sizes up to 100 and simple formulae for estimated critical values for tests at the 5% level.

Tests of MULTIMODE calculations of rovibrational energies of CH4

We report variational calculations of rovibrational energies of CH4 using the code MULTIMODE and an ab initio force field of Schwenke and Partridge. The systematic convergence of the energies with respect to the level of mode coupling is presented. Converged vibrational energies calculated using the five-mode representation of the potential for zero total angular momentum are compared with previous, benchmark calculations based on Radau coordinates using this force field for zero total angular momentum and for J = 1. Very good agreement with the previous benchmark calculations is found. (c) 2006 Elsevier B.V. All rights reserved.