Symmetric stability of compressible zonal flows on a generalized equatorial β plane

Sufficient conditions are derived for the linear stability with respect to zonally symmetric perturbations of a steady zonal solution to the nonhydrostatic compressible Euler equations on an equatorial � plane, including a leading order representation of the Coriolis force terms due to the poleward component of the planetary rotation vector. A version of the energy–Casimir method of stability proof is applied: an invariant functional of the Euler equations linearized about the equilibrium zonal flow is found, and positive definiteness of the functional is shown to imply linear stability of the equilibrium. It is shown that an equilibrium is stable if the potential vorticity has the same sign as latitude and the Rayleigh centrifugal stability condition that absolute angular momentum increase toward the equator on surfaces of constant pressure is satisfied. The result generalizes earlier results for hydrostatic and incompressible systems and for systems that do not account for the nontraditional Coriolis force terms. The stability of particular equilibrium zonal velocity, entropy, and density fields is assessed. A notable case in which the effect of the nontraditional Coriolis force is decisive is the instability of an angular momentum profile that decreases away from the equator but is flatter than quadratic in latitude, despite its satisfying both the centrifugal and convective stability conditions.

Neuromantic : from semi-manual to semi-automatic reconstruction of neuron morphology

The ability to create accurate geometric models of neuronal morphology is important for understanding the role of shape in information processing. Despite a significant amount of research on automating neuron reconstructions from image stacks obtained via microscopy, in practice most data are still collected manually. This paper describes Neuromantic, an open source system for three dimensional digital tracing of neurites. Neuromantic reconstructions are comparable in quality to those of existing commercial and freeware systems while balancing speed and accuracy of manual reconstruction. The combination of semi-automatic tracing, intuitive editing, and ability of visualizing large image stacks on standard computing platforms provides a versatile tool that can help address the reconstructions availability bottleneck. Practical considerations for reducing the computational time and space requirements of the extended algorithm are also discussed.

On the generalized eigenvalue problem for the Rossby wave vertical velocity in the presence of mean flow and topography

In a series of papers, Killworth and Blundell have proposed to study the effects of a background mean flow and topography on Rossby wave propagation by means of a generalized eigenvalue problem formulated in terms of the vertical velocity, obtained from a linearization of the primitive equations of motion. However, it has been known for a number of years that this eigenvalue problem contains an error, which Killworth was prevented from correcting himself by his unfortunate passing and whose correction is therefore taken up in this note. Here, the author shows in the context of quasigeostrophic (QG) theory that the error can ulti- mately be traced to the fact that the eigenvalue problem for the vertical velocity is fundamentally a non- linear one (the eigenvalue appears both in the numerator and denominator), unlike that for the pressure. The reason that this nonlinear term is lacking in the Killworth and Blundell theory comes from neglecting the depth dependence of a depth-dependent term. This nonlinear term is shown on idealized examples to alter significantly the Rossby wave dispersion relation in the high-wavenumber regime but is otherwise irrelevant in the long-wave limit, in which case the eigenvalue problems for the vertical velocity and pressure are both linear. In the general dispersive case, however, one should first solve the generalized eigenvalue problem for the pressure vertical structure and, if needed, diagnose the vertical velocity vertical structure from the latter.

The full infinite dimensional moment problem on semi-algebraic sets of generalized functions

We consider a generic basic semi-algebraic subset S of the space of generalized functions, that is a set given by (not necessarily countably many) polynomial constraints. We derive necessary and sufficient conditions for an infinite sequence of generalized functions to be realizable on S, namely to be the moment sequence of a finite measure concentrated on S. Our approach combines the classical results about the moment problem on nuclear spaces with the techniques recently developed to treat the moment problem on basic semi-algebraic sets of Rd. In this way, we determine realizability conditions that can be more easily verified than the well-known Haviland type conditions. Our result completely characterizes the support of the realizing measure in terms of its moments. As concrete examples of semi-algebraic sets of generalized functions, we consider the set of all Radon measures and the set of all the measures having bounded Radon–Nikodym density w.r.t. the Lebesgue measure.

Interactive effect of floral abundance and semi-natural habitats on pollinators in field beans (Vicia faba)

Pollination services are economically important component of agricultural biodiversity which enhance the yield and quality of many crops. An understanding of the suitability of extant habitats for pollinating species is crucial for planning management actions to protect and manage these service providers. In a highly modified agricultural ecosystem, we tested the effect of different pollination treatments (open, autonomous self- and wind-pollination) on pod set, seed set, and seed weight in field beans (Vicia faba). We also investigated the effect of semi-natural habitats and flower abundance on pollinators of field beans. Pollinator sampling was undertaken in ten field bean fields along a gradient of habitat complexity; CORINE land cover classification was used to analyse the land use patterns between 500–3000 m around the sites. Total yield from open-pollination increased by 185% compared to autonomous self-pollination. There was positive interactive effect of local flower abundance and cover of semi-natural habitats on overall abundance of pollinators at 1500 and 2000 m, and abundance of bumblebees (Bombus spp.) at 1000–2000 m. In contrast, species richness of pollinators was only correlated with flower abundance and not with semi-natural habitats. We did not find a link between pod set from open-pollination and pollinator abundance, possibly due to variations in the growing conditions and pollinator communities between sites. We conclude that insect pollination is essential for optimal bean yields and therefore the maintenance of semi-natural habitats in agriculture-dominated landscapes should ensure stable and more efficient pollination services in field beans.

The compositional and configurational heterogeneity of matrix habitats shape woodland carabid communities in wooded-agricultural landscapes

Landscape heterogeneity (the composition and configuration of matrix habitats) plays a major role in shaping species communities in wooded-agricultural landscapes. However, few studies consider the influence of different types of semi-natural and linear habitats in the matrix, despite their known ecological value for biodiversity. Objective To investigate the importance of the composition and configuration of matrix habitats for woodland carabid communities and identify whether specific landscape features can help to maintain long-term populations in wooded-agricultural environments. Methods Carabids were sampled from woodlands in 36 tetrads of 4 km2 across southern Britain. Landscape heterogeneity including an innovative representation of linear habitats was quantified for each tetrad. Carabid community response was analysed using ordination methods combined with variation partitioning and additional response trait analyses. Results Woodland carabid community response was trait-specific and better explained by simultaneously considering the composition and configuration of matrix habitats. Semi-natural and linear features provided significant refuge habitat and functional connectivity. Mature hedgerows were essential for slow-dispersing carabids in fragmented landscapes. Species commonly associated with heathland were correlated with inland water and woodland patches despite widespread heathland conversion to agricultural land, suggesting that species may persist for some decades when elements representative of the original habitat are retained following landscape modification. Conclusions Semi-natural and linear habitats have high biodiversity value. Landowners should identify features that can provide additional resources or functional connectivity for species relative to other habitat types in the landscape matrix. Agri-environment options should consider landscape heterogeneity to identify the most efficacious changes for biodiversity.

Influence of molecular composition on the development of microstructure from sheared polyethylene melts: Molecular and lamellar templating

A combination of in situ and ex situ X-ray scattering techniques and transmission electron microscopy has been used to study the crystallization behaviour of polyethylene, following the imposition of melt shear. In the case of a branched material, the imposition of shear flow up to a rate of 30 s(-1) was found to induce no anisotropy. Although shearing the linear material only ever induced a very small degree of anisotropy in the melt, for shear rates > 0.15 s(-1), subsequent crystallization resulted in increasing anisotropy. Blends of the above two polyethylenes were produced, in which the linear material constituted the minority fraction (similar to 10%). Isothermal crystallization at temperatures where extensive crystallization of the branched material does not occur demonstrated that the behaviour of the linear component of the sheared blend mirrored that of the linear polyethylene alone. However, in addition, it was found that when crystallized in the presence of an oriented morphology, the branched polymer also formed anisotropic structures. We have termed the process templating, in which the crystallization behaviour of the bulk of the system (similar to 90% branched material) is completely altered (spherulitic to oriented lamellar) by mapping it onto a pre-existing minority structure (similar to 10% linear polymer). (c) 2006 Elsevier Ltd. All rights reserved.

Boundary integral equations on unbounded rough surfaces: Fredholmness and the finite section method

We consider a class of boundary integral equations that arise in the study of strongly elliptic BVPs in unbounded domains of the form $D = \{(x, z)\in \mathbb{R}^{n+1} : x\in \mathbb{R}^n, z > f(x)\}$ where $f : \mathbb{R}^n \to\mathbb{R}$ is a sufficiently smooth bounded and continuous function. A number of specific problems of this type, for example acoustic scattering problems, problems involving elastic waves, and problems in potential theory, have been reformulated as second kind integral equations $u+Ku = v$ in the space $BC$ of bounded, continuous functions. Having recourse to the so-called limit operator method, we address two questions for the operator $A = I + K$ under consideration, with an emphasis on the function space setting $BC$. Firstly, under which conditions is $A$ a Fredholm operator, and, secondly, when is the finite section method applicable to $A$?

A record of Holocene climate change from lake geochemical analyses in southeastern Arabia

Lacustrine sediments from southeastern Arabia reveal variations in lake level corresponding to changes in the strength and duration of Indian Ocean Monsoon (IOM) summer rainfall and winter cyclonic rainfall. The late glacial/Holocene transition of the region was characterised by the development of mega-linear dunes. These dunes became stabilised and vegetated during the early Holocene and interdunal lakes formed in response to the incursion of the IOM at approximately 8500 cal yr BP with the development of C3 dominated savanna grasslands. The IOM weakened ca. 6000 cal yr BP with the onset of regional aridity, aeolian sedimentation and dune reactivation and accretion. Despite this reduction in precipitation, the take was maintained by winter dominated rainfall. There was a shift to drier adapted C4 grasslands across the dune field. Lake sediment geochemical analyses record precipitation minima at 8200, 5000 and 4200 cal yr BP that coincide with Bond events in the North Atlantic. A number of these events correspond with changes in cultural periods, suggesting that climate was a key mechanism affecting human occupation and exploitation of this region. (c) 2006 University of Washington. All rights reserved.

Effects of dietary supplementation with selenium enriched yeast or sodium selenite on selenium tissue distribution and meat quality in lambs

The objective was to determine the concentration of total selenium (Se) and the proportion of total Se comprised as selenomethionine (SeMet) and selenocysteine (SeCys), as well as meat quality in terms of oxidative stability in post mortem tissues of lambs offered diets with an increasing dose rate of selenized enriched yeast (SY), or sodium selenite (SS). Fifty lambs were offered, for a period of 112 d, a total mixed ration which had either been supplemented with SY (0, 0.11, 0.21 or 0.31 mg/kg DM to give total Se contents of 0.19, 0.3, 0.4 and 0.5 mg Se/kg DM for treatments T1, T2, T3 and T4, respectively) or SS (0.11 mg/kg DM to give 0.3 mg Se/kg DM total Se [T5]). At enrolment and at 28, 56, 84 and 112 d following enrolment, blood samples were taken for Se and Se species determination, as well as glutathione peroxidase (GSH-Px) activity. At the end of the study lambs were euthanased and samples of heart, liver, kidney, and skeletal muscle were retained for Se and Se species determination. Tissue GSH-Px activity and thiobarbituric acid reactive substances (TBARS) were determined in Longissimus Thoracis. The incorporation into the diet of ascending concentrations of Se as SY increased whole blood total Se and the proportion of total Se comprised as SeMet, and erythrocyte GSH-Px activity. Comparable doses of SS supplementation did not result in significant differences between these parameters. With the exception of kidney tissue, all other tissues showed a dose dependant response to increasing concentrations of dietary SY, such that total Se and SeMet increased. Selenium content of Psoas Major was higher in animals fed SY when compared to a similar dose of SS, indicating improvements in Se availability and retention. There were no significant treatment effects on meat quality assessments GHS-Px and TBARS, reflecting the lack of difference in the proportion of total Se that was comprised as SeCys. However, oxidative stability improved marginally with ascending tissue Se content, providing an indication of a linear dose response whereby TBARS improved with ascending SY inclusion.

Multi-mode approximations to wave scattering by an uneven bed

Approximations to the scattering of linear surface gravity waves on water of varying quiescent depth are Investigated by means of a variational approach. Previous authors have used wave modes associated with the constant depth case to approximate the velocity potential, leading to a system of coupled differential equations. Here it is shown that a transformation of the dependent variables results in a much simplified differential equation system which in turn leads to a new multi-mode 'mild-slope' approximation. Further, the effect of adding a bed mode is examined and clarified. A systematic analytic method is presented for evaluating inner products that arise and numerical experiments for two-dimensional scattering are used to examine the performance of the new approximations.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

Diagnosis of variability and trends in a global precipitation dataset using a physically motivated statistical model

A physically motivated statistical model is used to diagnose variability and trends in wintertime ( October - March) Global Precipitation Climatology Project (GPCP) pentad (5-day mean) precipitation. Quasi-geostrophic theory suggests that extratropical precipitation amounts should depend multiplicatively on the pressure gradient, saturation specific humidity, and the meridional temperature gradient. This physical insight has been used to guide the development of a suitable statistical model for precipitation using a mixture of generalized linear models: a logistic model for the binary occurrence of precipitation and a Gamma distribution model for the wet day precipitation amount. The statistical model allows for the investigation of the role of each factor in determining variations and long-term trends. Saturation specific humidity q(s) has a generally negative effect on global precipitation occurrence and with the tropical wet pentad precipitation amount, but has a positive relationship with the pentad precipitation amount at mid- and high latitudes. The North Atlantic Oscillation, a proxy for the meridional temperature gradient, is also found to have a statistically significant positive effect on precipitation over much of the Atlantic region. Residual time trends in wet pentad precipitation are extremely sensitive to the choice of the wet pentad threshold because of increasing trends in low-amplitude precipitation pentads; too low a choice of threshold can lead to a spurious decreasing trend in wet pentad precipitation amounts. However, for not too small thresholds, it is found that the meridional temperature gradient is an important factor for explaining part of the long-term trend in Atlantic precipitation.

Perturbed Rankine vortices in surface quasi-geostrophic dynamics

An analytical dispersion relation is derived for linear perturbations to a Rankine vortex governed by surface quasi-geostrophic dynamics. Such a Rankine vortex is a circular region of uniform anomalous surface temperature evolving under quasi-geostrophic dynamics with uniform interior potential vorticity. The dispersion relation is analysed in detail and compared to the more familiar dispersion relation for a perturbed Rankine vortex governed by the Euler equations. The results are successfully verified against numerical simulations of the full equations. The dispersion relation is relevant to problems including wave propagation on surface temperature fronts and the stability of vortices in quasi-geostrophic turbulence.

Initial boundary value problems for the nonlinear Schrödinger equation

A new spectral method for solving initial boundary value problems for linear and integrable nonlinear partial differential equations in two independent variables is applied to the nonlinear Schrödinger equation and to its linearized version in the domain {x≥l(t), t≥0}. We show that there exist two cases: (a) if l″(t)<0, then the solution of the linear or nonlinear equations can be obtained by solving the respective scalar or matrix Riemann-Hilbert problem, which is defined on a time-dependent contour; (b) if l″(t)>0, then the Riemann-Hilbert problem is replaced by a respective scalar or matrix problem on a time-independent domain. In both cases, the solution is expressed in a spectrally decomposed form.