Predicting mesh density for adaptive modelling of the global atmosphere

The shallow water equations are solved using a mesh of polygons on the sphere, which adapts infrequently to the predicted future solution. Infrequent mesh adaptation reduces the cost of adaptation and load-balancing and will thus allow for more accurate mapping on adaptation. We simulate the growth of a barotropically unstable jet adapting the mesh every 12 h. Using an adaptation criterion based largely on the gradient of the vorticity leads to a mesh with around 20 per cent of the cells of a uniform mesh that gives equivalent results. This is a similar proportion to previous studies of the same test case with mesh adaptation every 1–20 min. The prediction of the mesh density involves solving the shallow water equations on a coarse mesh in advance of the locally refined mesh in order to estimate where features requiring higher resolution will grow, decay or move to. The adaptation criterion consists of two parts: that resolved on the coarse mesh, and that which is not resolved and so is passively advected on the coarse mesh. This combination leads to a balance between resolving features controlled by the large-scale dynamics and maintaining fine-scale features.

Adjoint goal-based error norms for adaptive mesh ocean modelling

Flow in the world's oceans occurs at a wide range of spatial scales, from a fraction of a metre up to many thousands of kilometers. In particular, regions of intense flow are often highly localised, for example, western boundary currents, equatorial jets, overflows and convective plumes. Conventional numerical ocean models generally use static meshes. The use of dynamically-adaptive meshes has many potential advantages but needs to be guided by an error measure reflecting the underlying physics. A method of defining an error measure to guide an adaptive meshing algorithm for unstructured tetrahedral finite elements, utilizing an adjoint or goal-based method, is described here. This method is based upon a functional, encompassing important features of the flow structure. The sensitivity of this functional, with respect to the solution variables, is used as the basis from which an error measure is derived. This error measure acts to predict those areas of the domain where resolution should be changed. A barotropic wind driven gyre problem is used to demonstrate the capabilities of the method. The overall objective of this work is to develop robust error measures for use in an oceanographic context which will ensure areas of fine mesh resolution are used only where and when they are required. (c) 2006 Elsevier Ltd. All rights reserved.

Why high seed densities within buried mesh bags may overestimate depletion rates of soil seed banks

1. Estimates of seed bank depletion rates are essential for modelling and management of plant populations. The seed bag burial method is often used to measure seed mortality in the soil. However, the density of seeds within seed bags is higher than densities in natural seed banks, which may elevate levels of pathogens and influence seed mortality. The aim of this study was to quantify the effects of fungi and seed density within buried mesh bags on the mortality of seeds. Striga hermonthica was chosen as the study species because it has been widely studied but different methods for measuring seed mortality in the soil have yielded contradictory estimates. 2. Seed bags were buried in soil and exhumed at regular time intervals to monitor mortality of the seeds in three field experiments during two rainy seasons. The effect of fungal activity on seed mortality was evaluated in a fungi exclusion experiment. Differences in seed-to-seed interaction were obtained by using two and four densities within the seed bags in consecutive years. Densities were created by mixing 1000 seeds with 0, 10, 100 or 1000 g of coarse sand. 3. The mortality rate was significantly lower when fungi were excluded, indicating the possible role of pathogenic fungi. 4. Decreasing the density of seeds in bags significantly reduced seed mortality, most probably because of decreased seed-to-seed contamination by pathogenic fungi. 5. Synthesis and applications. Models of plant populations in general and annual weeds in particular often use values from the literature for seed bank depletion rates. These depletion rates have often been estimated by the seed bag burial method, yet seed density within seed bags may be unrealistically high. Consequently, estimates of seed mortality rates may be too high because of an overestimation of the effects of soil or seed-borne pathogens. Species that have been classified from such studies as having short-lived seed banks may need to be re-assessed using realistic densities either within seed bags or otherwise. Similarly, models of seed bank dynamics based on such overestimated depletion rates may lead to incorrect conclusions regarding the seed banks and, perhaps, the management of weeds and rare species.

A non-uniform mesh scheme for compressible flow

A non-uniform mesh scheme is presented for the computation of compressible flows governed by the Euler equations of gas dynamics. The scheme is based on flux-difference splitting and represents an extension of a similar scheme designed for uniform meshes. The numerical results demonstrate that little, if any, spurious oscillation occurs as a result of the non-uniformity of the mesh; and importantly, shock speeds are computed correctly.

Progressive erosion of genetic and epigenetic variation in callus-derived cocoa (Theobroma cacao) plants

Relatively little is known about the timing of genetic and epigenetic forms of somaclonal variation arising from callus growth. We surveyed for both types of change in cocoa (Theobroma cacao) plants regenerated from calli of various ages, and also between tissues from the source trees. For genetic change, we used 15 single sequence repeat (SSR) markers from four source trees and from 233 regenerated plants. For epigenetic change, we used 386 methylation-sensitive amplified polymorphism (MSAP) markers on leaf and explant (staminode) DNA from two source trees and on leaf DNA from 114 regenerants. Genetic variation within source trees was limited to one slippage mutation in one leaf. Regenerants were far more variable, with 35% exhibiting at least one mutation. Genetic variation initially accumulated with culture age but subsequently declined. MSAP (epigenetic) profiles diverged between leaf and staminode samples from source trees. Multivariate analysis revealed that leaves from regenerants occupied intermediate eigenspace between leaves and staminodes of source plants but became progressively more similar to source tree leaves with culture age. Statistical analysis confirmed this rather counterintuitive finding that leaves of ‘late regenerants’ exhibited significantly less genetic and epigenetic divergence from source leaves than those exposed to short periods of callus growth.

A moving-mesh finite element method and its application to the numerical solution of phase-change problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase. The algorithm uses a distributed conservation principle to determine nodal mesh velocities, which are then used to move the nodes. The nodal values are obtained from an ALE (Arbitrary Lagrangian-Eulerian) equation, which represents a generalization of the original algorithm presented in Applied Numerical Mathematics, 54:450--469 (2005). Having described the details of the generalized algorithm it is validated on two test cases from the original paper and is then applied to one-phase and, for the first time, two-phase Stefan problems in one and two space dimensions, paying particular attention to the implementation of the interface boundary conditions. Results are presented to demonstrate the accuracy and the effectiveness of the method, including comparisons against analytical solutions where available.