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.

Will Extratropical Storms Intensify in a Warmer Climate?

Extratropical cyclones and how they may change in a warmer climate have been investigated in detail with a high-resolution version of the ECHAM5 global climate model. A spectral resolution of T213 (63 km) is used for two 32-yr periods at the end of the twentieth and twenty-first centuries and integrated for the Intergovernmental Panel on Climate Change (IPCC) A1B scenario. Extremes of pressure, vorticity, wind, and precipitation associated with the cyclones are investigated and compared with a lower-resolution simulation. Comparison with observations of extreme wind speeds indicates that the model reproduces realistic values. This study also investigates the ability of the model to simulate extratropical cyclones by computing composites of intense storms and contrasting them with the same composites from the 40-yr ECMWF Re-Analysis (ERA-40). Composites of the time evolution of intense cyclones are reproduced with great fidelity; in particular the evolution of central surface pressure is almost exactly replicated, but vorticity, maximum wind speed, and precipitation are higher in the model. Spatial composites also show that the distributions of pressure, winds, and precipitation at different stages of the cyclone life cycle compare well with those from ERA-40, as does the vertical structure. For the twenty-first century, changes in the distribution of storms are very similar to those of previous study. There is a small reduction in the number of cyclones but no significant changes in the extremes of wind and vorticity in both hemispheres. There are larger regional changes in agreement with previous studies. The largest changes are in the total precipitation, where a significant increase is seen. Cumulative precipitation along the tracks of the cyclones increases by some 11% per track, or about twice the increase in global precipitation, while the extreme precipitation is close to the globally averaged increase in column water vapor (some 27%). Regionally, changes in extreme precipitation are even higher because of changes in the storm tracks.

Cue combination for 3D location judgements

Cue combination rules have often been applied to the perception of surface shape but not to judgements of object location. Here, we used immersive virtual reality to explore the relationship between different cues to distance. Participants viewed a virtual scene and judged the change in distance of an object presented in two intervals, where the scene changed in size between intervals (by a factor of between 0.25 and 4). We measured thresholds for detecting a change in object distance when there were only 'physical' (stereo and motion parallax) or 'texture-based' cues (independent of the scale of the scene) and used these to predict biases in a distance matching task. Under a range of conditions, in which the viewing distance and position of the tarte relative to other objects was varied, the ration of 'physical' to 'texture-based' thresholds was a good predictor of biases in the distance matching task. The cue combination approach, which successfully accounts for our data, relies on quite different principles from those underlying geometric reconstruction.

Latitude and longitude vertical disparities

The literature on vertical disparity is complicated by the fact that several different definitions of the term “vertical disparity” are in common use, often without a clear statement about which is intended or a widespread appreciation of the properties of the different definitions. Here, we examine two definitions of retinal vertical disparity: elevation-latitude and elevation-longitude disparities. Near the fixation point, these definitions become equivalent, but in general, they have quite different dependences on object distance and binocular eye posture, which have not previously been spelt out. We present analytical approximations for each type of vertical disparity, valid for more general conditions than previous derivations in the literature: we do not restrict ourselves to objects near the fixation point or near the plane of regard, and we allow for non-zero torsion, cyclovergence, and vertical misalignments of the eyes. We use these expressions to derive estimates of the latitude and longitude vertical disparities expected at each point in the visual field, averaged over all natural viewing. Finally, we present analytical expressions showing how binocular eye position—gaze direction, convergence, torsion, cyclovergence, and vertical misalignment—can be derived from the vertical disparity field and its derivatives at the fovea.

Stereo and motion parallax cues in human 3D vision: can they vanish without a trace?

In an immersive virtual reality environment, subjects fail to notice when a scene expands or contracts around them, despite correct and consistent information from binocular stereopsis and motion parallax, resulting in gross failures of size constancy (A. Glennerster, L. Tcheang, S. J. Gilson, A. W. Fitzgibbon, & A. J. Parker, 2006). We determined whether the integration of stereopsis/motion parallax cues with texture-based cues could be modified through feedback. Subjects compared the size of two objects, each visible when the room was of a different size. As the subject walked, the room expanded or contracted, although subjects failed to notice any change. Subjects were given feedback about the accuracy of their size judgments, where the “correct” size setting was defined either by texture-based cues or (in a separate experiment) by stereo/motion parallax cues. Because of feedback, observers were able to adjust responses such that fewer errors were made. For texture-based feedback, the pattern of responses was consistent with observers weighting texture cues more heavily. However, for stereo/motion parallax feedback, performance in many conditions became worse such that, paradoxically, biases moved away from the point reinforced by the feedback. This can be explained by assuming that subjects remap the relationship between stereo/motion parallax cues and perceived size or that they develop strategies to change their criterion for a size match on different trials. In either case, subjects appear not to have direct access to stereo/motion parallax cues.

In vitro degradability of mature and immature leaves of tropical forage legumes differing in condensed tannin and non-starch polysaccharide content and composition

A study was designed to examine the relationships between protein, condensed tannin and cell wall carbohydrate content and composition and the nutritional quality of seven tropical legumes (Desmodium ovalifolium, Flemingia macrophylla, Leucaena leucocephala, L pallida, L macrophylla, Calliandra calothyrsus and Clitotia fairchildiana). Among the legume species studied, D ovalifolium showed the lowest concentration of nitrogen, while L leucocephala showed the highest. Fibre (NDF) content was lowest in C calothyrsus, L Leucocephala and L pallida and highest in L macrophylla, which had no measurable condensed tannins. The highest tannin concentration was found in C calothyrsus. Total non-structural polysaccharides (NSP) varied among legumes species (lowest in C calothyrsus and highest in D ovalifolium), and glucose and uronic acids were the most abundant carbohydrate constituents in all legumes. Total NSP losses were lowest in F macrophylla and highest in L leucocephala and L pallida. Gas accumulation and acetate and propionate levels were 50% less with F macrophylla and D ovalifolium as compared with L leucocephala. The highest levels of branched-chain fatty acids were observed with non-tanniniferous legumes, and negative concentrations were observed with some of the legumes with high tannin content (D ovalifolium and F macrophylla). Linear regression analysis showed that the presence of condensed tannins was more related to a reduction of the initial rate of gas production (0-48 h) than to the final amount of gas produced or the extent (144h) of dry matter degradation, which could be due to differences in tannin chemistry. Consequently, more attention should be given in the future to elucidating the impact of tannin structure on the nutritional quality of tropical forage legumes. (C) 2003 Society of Chemical Industry.

In vitro biological activity of tannins from Acacia and other tree fruits: Correlations with colorimetric and gravimetric phenolic assays

This Study was designed to investigate impact of tannins on in vitro ruminal fermentation parameters as well as relationships between concentration and in vitro biological activity of tannins present in tree fruits. Dry and mature fruits of known phenolic content harvested from Acacia nilotica, A. erubescens, A. erioloba, A. sieberiana, Piliostigima thonningii and Dichrostachys cinerea tree species were fermented with rumen fluid in vitro with or without polyethylene glycol (PEG). Correlation between in vitro biological activity and phenolic concentration was determined. Polyethylene glycol inclusion increased Cumulative gas production from all fruit substrates. The largest Increase (225%) after 48 h incubation was observed in D. cinerea fruits while the least (12.7%) increase was observed in A. erubescens fruits. Organic matter degradability (48 h) was increased by PEG inclusion for all tree species except A. erubescens and P. thonningii. For D. cinerea fruits, colorimetric assays were poorly correlated to Increases In gas production due to PEG treatment. Ytterbium precipitable phenolics (YbPh) were also poorly correlated with response to PEG for A. erioloba and P. thonningii fruits. However, YbPh were strongly and positively correlated to the increase In Cumulative gas production due to PEG for A. erubescens and A. nilotica. Folin-Ciocalteau assayed phenolics (SPh) were not correlated to response to PEG in P. thonningii and A. sieberiana. It was Concluded that the PEG effect oil in vitro fermentation was closely related to some measures of phenolic concentration but the relationships varied with tree species.

The nutritive value of forage sorghum genotypes developed for the dry tropical highlands of Kenya as feed source for ruminants

in vitro studies were conducted on five sorghum genotypes developed for the dry tropical highland climate of Kenya and which can be fed to ruminants fresh or as silage. The five sorghum genotypes consisted of two normal white mid-rib (WMR) genotypes, coded E1291 and E65181, and three brown-midrib (BMR) genotypes, coded Lan-5, Lan-6 and Lan-12. Whole mature plants (herbage plus grain) and silage made from E 1291 were used in the study. An in vitro manual gas production technique was used to compare the nutritive characteristics of these genotypes for ruminants. These sorghums differed significantly in true organic matter degraded (OMDeg), which ranged from 520 to 678 g/kg after 24 h incubation and 706 to 805 g/kg after 72 h incubation. All the BMR sorghums had a higher degradability than the WMR genotype, E6518, and the silage, with Lan-5 having the highest degradability. Methane produced per g OMDeg ranged from 40.6 to 46.4 mL/g after 24 h incubation and from 53.1 to 62.6 mL/g after 72 h incubation. It was similar for all genotypes after 24 h incubation but Lan-12 had the highest methane production after 72 h incubation. After 24 h and 72 h incubation all the genotypes produced a similar total amount of gas per OMDeg (293 to 309 and 357 to 385 mL/g, respectively) with similar total short chain fatty acid concentrations in the liquid digesta (7.8 to 10.4 and 9.5 to 10.3 mmol, respectively) and acetate to propionate ratios of 2.16 to 2.49 and 2.35 to 2.87, respectively. The sorghums showed great potential as ruminant feed sources in the region.

An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genotnes. Several statistical methods have been proposed to study the genetic structure using AFLP's but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F-IS. A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F-IS and F-ST values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.uif-grenoble.fi-/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.

Determination of mitochondrial genetic diversity in mammals

Mitochondrial DNA (mtDNA) is one of the most Popular population genetic markers. Its relevance as an indicator Of Population size and history has recently been questioned by several large-scale studies in animals reporting evidence for recurrent adaptive evolution, at least in invertebrates. Here we focus on mammals, a more restricted taxonomic group for which the issue of mtDNA near neutrality is crucial. By analyzing the distribution of mtDNA diversity across species and relating 4 to allozyme diversity, life-history traits, and taxonomy, we show that (i) mtDNA in mammals (toes not reject the nearly neutral model; (ii) mtDNA diversity, however, is unrelated to any of the 14 life-history and ecological variables that we analyzed, including body mass, geographic range, and The World Conservation Union (IUCN) categorization; (iii) mtDNA diversity is highly variable between mammalian orders and families; (iv) this taxonomic effect is most likely explained by variations of mutation rate between lineages. These results are indicative of a strong stochasticity of effective population size in mammalian species. They Suggest that, even in the absence of selection, mtDNA genetic diversity is essentially unpredictable, knowing species biology, and probably uncorrelated to species abundance.

The population genetic structure of clonal organisms generated by exponentially bounded and fat-tailed dispersal

Long distance dispersal (LDD) plays an important role in many population processes like colonization, range expansion, and epidemics. LDD of small particles like fungal spores is often a result of turbulent wind dispersal and is best described by functions with power-law behavior in the tails ("fat tailed"). The influence of fat-tailed LDD on population genetic structure is reported in this article. In computer simulations, the population structure generated by power-law dispersal with exponents in the range of -2 to -1, in distinct contrast to that generated by exponential dispersal, has a fractal structure. As the power-law exponent becomes smaller, the distribution of individual genotypes becomes more self-similar at different scales. Common statistics like G(ST) are not well suited to summarizing differences between the population genetic structures. Instead, fractal and self-similarity statistics demonstrated differences in structure arising from fat-tailed and exponential dispersal. When dispersal is fat tailed, a log-log plot of the Simpson index against distance between subpopulations has an approximately constant gradient over a large range of spatial scales. The fractal dimension D-2 is linearly inversely related to the power-law exponent, with a slope of similar to -2. In a large simulation arena, fat-tailed LDD allows colonization of the entire space by all genotypes whereas exponentially bounded dispersal eventually confines all descendants of a single clonal lineage to a relatively small area.