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.

On the initiation of surface waves by turbulent shear flow

An analytical model is developed for the initial stage of surface wave generation at an air-water interface by a turbulent shear flow in either the air or in the water. The model treats the problem of wave growth departing from a flat interface and is relevant for small waves whose forcing is dominated by turbulent pressure fluctuations. The wave growth is predicted using the linearised and inviscid equations of motion, essentially following Phillips [Phillips, O.M., 1957. On the generation of waves by turbulent wind. J. Fluid Mech. 2, 417-445], but the pressure fluctuations that generate the waves are treated as unsteady and related to the turbulent velocity field using the rapid-distortion treatment of Durbin [Durbin, P.A., 1978. Rapid distortion theory of turbulent flows. PhD thesis, University of Cambridge]. This model, which assumes a constant mean shear rate F, can be viewed as the simplest representation of an oceanic or atmospheric boundary layer. For turbulent flows in the air and in the water producing pressure fluctuations of similar magnitude, the waves generated by turbulence in the water are found to be considerably steeper than those generated by turbulence in the air. For resonant waves, this is shown to be due to the shorter decorrelation time of turbulent pressure in the air (estimated as proportional to 1/Gamma), because of the higher shear rate existing in the air flow, and due to the smaller length scale of the turbulence in the water. Non-resonant waves generated by turbulence in the water, although being somewhat gentler, are still steeper than resonant waves generated by turbulence in the air. Hence, it is suggested that turbulence in the water may have a more important role than previously thought in the initiation of the surface waves that are subsequently amplified by feedback instability mechanisms.

Solar modulation in surface atmospheric electricity

The solar wind modulates the flux of galactic cosmic rays impinging on Earth inversely with solar activity. Cosmic ray ionisation is the major source of air’s electrical conductivity over the oceans and well above the continents. Differential solar modulation of the cosmic ray energy spectrum modifies the cosmic ray ionisation at different latitudes,varying the total atmospheric columnar conductance. This redistributes current flow in the global atmospheric electrical circuit, including the local vertical current density and the related surface potential gradient. Surface vertical current density and potential gradient measurements made independently at Lerwick Observatory,Shetland,from 1978 to 1985 are compared with modelled changes in cosmic ray ionisation arising from solar activity changes. Both the lower troposphere atmospheric electricity quantities are significantly increased at cosmic ray maximum(solar minimum),with a proportional change greater than that of the cosmic ray change.

Taxonomical vs. functional responses of bee communities to fire in two contrasting climatic regions

Valuable insights into mechanisms of community responses to environmental change can be gained by analysing in tandem the variation in functional and taxonomic composition along environmental gradients. We assess the changes in species and functional trait composition (i.e. dominant traits and functional diversity) of diverse bee communities in contrasting fire-driven systems in two climatic regions: Mediterranean (scrub habitats in Israel) and temperate (chestnut forests in southern Switzerland). In both climatic regions, there were shifts in species diversity and composition related to post-fire age. In the temperate region, functional composition responded markedly to fire; however, in the Mediterranean, the taxonomic response to fire was not matched by functional replacement. These results suggest that greater functional stability to fire in the Mediterranean is achieved by replacement of functionally similar species (i.e. functional redundancy) which dominate under different environmental conditions in the heterogeneous landscapes of the region. In contrast, the greater functional response in the temperate region was attributed to a more rapid post-fire vegetation recovery and shorter time-window when favourable habitat was available relative to the Mediterranean. Bee traits can be used to predict the functional responses of bee communities to environmental changes in habitats of conservation importance in different regions with distinct disturbance regimes. However, predictions cannot be generalized from one climatic region to another where distinct habitat configurations occur.

Rapeseed cytoplasm gives advantage in wild relatives and complicates genetically modified crop biocontainment

Biocontainment methods for genetically modified crops closest to commercial reality (chloroplast transformation, male sterility) would be compromised (in absolute terms) by seed-mediated gene flow leading to chloroplast capture. Even in these circumstances, however, it can be argued that biocontainment still represses transgene movement, with the efficacy depending on the relative frequency of seed-and pollen-mediated gene flow. In this study, we screened for crop-specific chloroplast markers from rapeseed (Brassica napus) amongst sympatric and allopatric populations of wild B. oleracea in natural cliff-top populations and B. rapa in riverside and weedy populations. We found only modest crop chloroplast presence in wild B. oleracea and in weedy B. rapa, but a surprisingly high incidence in sympatric (but not in allopatric) riverside B. rapa populations. Chloroplast inheritance models indicate that elevated crop chloroplast acquisition is best explained if crop cytoplasm confers selective advantage in riverside B. rapa populations. Our results therefore imply that chloroplast transformation may slow transgene recruitment in two settings, but actually accelerate transgene spread in a third. This finding suggests that the appropriateness of chloroplast transformation for biocontainment policy depends on both context and geographical location.

Gene flow, risk assessment and the environmental release of transgenic plants

The release of genetically modified plants is governed by regulations that aim to provide an assessment of potential impact on the environment. One of the most important components of this risk assessment is an evaluation of the probability of gene flow. In this review, we provide an overview of the current literature on gene flow from transgenic plants, providing a framework of issues for those considering the release of a transgenic plant into the environment. For some plants gene flow from transgenic crops is well documented, and this information is discussed in detail in this review. Mechanisms of gene flow vary from plant species to plant species and range from the possibility of asexual propagation, short- or long-distance pollen dispersal mediated by insects or wind and seed dispersal. Volunteer populations of transgenic plants may occur where seed is inadvertently spread during harvest or commercial distribution. If there are wild populations related to the transgenic crop then hybridization and eventually introgression in the wild may occur, as it has for herbicide resistant transgenic oilseed rape (Brassica napus). Tools to measure the amount of gene flow, experimental data measuring the distance of pollen dispersal, and experiments measuring hybridization and seed survivability are discussed in this review. The various methods that have been proposed to prevent gene flow from genetically modified plants are also described. The current "transgenic traits'! in the major crops confer resistance to herbicides and certain insects. Such traits could confer a selective advantage (an increase in fitness) in wild plant populations in some circumstances, were gene flow to occur. However, there is ample evidence that gene flow from crops to related wild species occurred before the development of transgenic crops and this should be taken into account in the risk assessment process.

Modelling oil absorption during post-frying cooling - I: model development

Several studies have highlighted the importance of the cooling period in oil absorption in deep-fat fried products. Specifically, it has been established that the largest proportion of oil which ends up into the food, is sucked into the porous crust region after the fried product is removed from the oil bath, stressing the importance of this time interval. The main objective of this paper was to develop a predictive mechanistic model that can be used to understand the principles behind post-frying cooling oil absorption kinetics, which can also help identifying the key parameters that affect the final oil intake by the fried product. The model was developed for two different geometries, an infinite slab and an infinite cylinder, and was divided into two main sub-models, one describing the immersion frying period itself and the other describing the post-frying cooling period. The immersion frying period was described by a transient moving-front model that considered the movement of the crust/core interface, whereas post-frying cooling oil absorption was considered to be a pressure driven flow mediated by capillary forces. A key element in the model was the hypothesis that oil suction would only begin once a positive pressure driving force had developed. The mechanistic model was based on measurable physical and thermal properties, and process parameters with no need of empirical data fitting, and can be used to study oil absorption in any deep-fat fried product that satisfies the assumptions made.

Determination of digesta flow entering the omasal canal of dairy cows using different marker systems

Four studies were conducted to compare the effect of four indigestible markers (LiCoEDTA, Yb-acetate, Cr-mordanted straw and indigestible neutral-detergent fibre (INDF)) and three marker systems on the flow of digesta entering the omasal canal of lactating dairy cows. Samples of digesta aspirated from the omasal canal were pooled and separated using filtration and high-speed centrifugation into three fractions defined as the liquid phase, small particulate and large particulate matter. Co was primarily associated with the liquid phase, Yb was concentrated in small particulate matter, whilst Cr and INDF were associated with large particles. Digesta flow was calculated based on single markers or using the reconstitution system based on combinations of two (Co + Yb, Co + Cr and Co + INDF) or three markers (Co + Yb + Cr and Co + Yb + INDF). Use of single markers resulted in large differences between estimates of organic matter (OM) flow entering the omasal canal suggesting that samples were not representative of true digesta. Digesta appeared to consist of at least three phases that tended to separate during sampling. OM was concentrated in particulate matter, whilst the liquid phase consisted mainly of volatile fatty acids and inorganic matter. Yb was intimately associated with nitrogenous compounds, whereas Cr and INDF were concentrated in fibrous material. Current data indicated that marker systems based on Yb in combination with Cr or INDF are required for the accurate determination of OM, N and neutral-detergent fibre flow. In cases where the flow of water-soluble nutrients entering the omasal canal is also required, the marker system should also include Co.

Results of applying sensor fusion to a control system using optic flow

The results from applying a sensor fusion process to an adaptive controller used to balance all inverted pendulum axe presented. The goal of the sensor fusion process was to replace some of the four mechanical measurements, which are known to be sufficient inputs for a linear state feedback controller to balance the system, with optic flow variables. Results from research into the psychology of the sense of balance in humans were the motivation for the investigation of this new type of controller input. The simulated model of the inverted pendulum and the virtual reality environments used to provide the optical input are described. The successful introduction of optical information is found to require the preservation of at least two of the traditional input types and entail increased training till-le for the adaptive controller and reduced performance (measured as the time the pendulum remains upright)

A policy model for secure information flow

When a computer program requires legitimate access to confidential data, the question arises whether such a program may illegally reveal sensitive information. This paper proposes a policy model to specify what information flow is permitted in a computational system. The security definition, which is based on a general notion of information lattices, allows various representations of information to be used in the enforcement of secure information flow in deterministic or nondeterministic systems. A flexible semantics-based analysis technique is presented, which uses the input-output relational model induced by an attacker's observational power, to compute the information released by the computational system. An illustrative attacker model demonstrates the use of the technique to develop a termination-sensitive analysis. The technique allows the development of various information flow analyses, parametrised by the attacker's observational power, which can be used to enforce what declassification policies.

Link between DOC in near surface peat and stream water in an upland catchment

Hydrologic transport of dissolved organic carbon (DOC) from peat soils may differ to organo-mineral soils in how they responded to changes in flow, because of differences in soil profile and hydrology. In well-drained organo-mineral soils, low flow is through the lower mineral layer where DOC is absorbed and high flow is through the upper organic layer where DOC is produced. DOC concentrations in streams draining organo-mineral soils typically increase with flow. In saturated peat soils, both high and low flows are through an organic layer where DOC is produced. Therefore, DOC in stream water draining peat may not increase in response to changes in flow as there is no switch in flow path between a mineral and organic layer. To verify this, we conducted a high-resolution monitoring study of soil and stream water at an upland peat catchment in northern England. Our data showed a strong positive correlation between DOC concentrations at − 1 and − 5 cm depth and stream water, and weaker correlations between concentrations at − 20 to − 50 cm depth and stream water. Although near surface organic material appears to be the key source of stream water DOC in both peat and organo-mineral soils, we observed a negative correlation between stream flow and DOC concentrations instead of a positive correlation as DOC released from organic layers during low and high flow was diluted by rainfall. The differences in DOC transport processes between peat and organo-mineral soils have different implications for our understanding of long-term changes in DOC exports. While increased rainfall may cause an increase in DOC flux from peat due to an increase in water volume, it may cause a decrease in concentrations. This response is contrary to expected changes in DOC exports from organo-mineral soils, where increase rainfall is likely to result in an increase in flux and concentration.

Insights into redox cycling of sulfur and iron in peatlands using high-resolution diffusive equilibrium thin film (DET) gel probe sampling

High spatial resolution vertical profiles of pore-water chemistry have been obtained for a peatland using diffusive equilibrium in thin films (DET) gel probes. Comparison of DET pore-water data with more traditional depth-specific sampling shows good agreement and the DET profiling method is less invasive and less likely to induce mixing of pore-waters. Chloride mass balances as water tables fell in the early summer indicate that evaporative concentration dominates and there is negligible lateral flow in the peat. Lack of lateral flow allows element budgets for the same site at different times to be compared. The high spatial resolution of sampling also enables gradients to be observed that permit calculations of vertical fluxes. Sulfate concentrations fall at two sites with net rates of 1.5 and 5.0nmol cm− 3 day− 1, likely due to a dominance of bacterial sulfate reduction, while a third site showed a net gain in sulfate due to oxidation of sulfur over the study period at an average rate of 3.4nmol cm− 3 day− 1. Behaviour of iron is closely coupled to that of sulfur; there is net removal of iron at the two sites where sulfate reduction dominates and addition of iron where oxidation dominates. The profiles demonstrate that, in addition to strong vertical redox related chemical changes, there is significant spatial heterogeneity. Whilst overall there is evidence for net reduction of sulfate within the peatland pore-waters, this can be reversed, at least temporarily, during periods of drought when sulfide oxidation with resulting acid production predominates.

An efficient shock capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.