The analysis of ranked observations of soil structure using indicator geostatistics

Structure is an important physical feature of the soil that is associated with water movement, the soil atmosphere, microorganism activity and nutrient uptake. A soil without any obvious organisation of its components is known as apedal and this state can have marked effects on several soil processes. Accurate maps of topsoil and subsoil structure are desirable for a wide range of models that aim to predict erosion, solute transport, or flow of water through the soil. Also such maps would be useful to precision farmers when deciding how to apply nutrients and pesticides in a site-specific way, and to target subsoiling and soil structure stabilization procedures. Typically, soil structure is inferred from bulk density or penetrometer resistance measurements and more recently from soil resistivity and conductivity surveys. To measure the former is both time-consuming and costly, whereas observations made by the latter methods can be made automatically and swiftly using a vehicle-mounted penetrometer or resistivity and conductivity sensors. The results of each of these methods, however, are affected by other soil properties, in particular moisture content at the time of sampling, texture, and the presence of stones. Traditional methods of observing soil structure identify the type of ped and its degree of development. Methods of ranking such observations from good to poor for different soil textures have been developed. Indicator variograms can be computed for each category or rank of structure and these can be summed to give the sum of indicator variograms (SIV). Observations of the topsoil and subsoil structure were made at four field sites where the soil had developed on different parent materials. The observations were ranked by four methods and indicator and the sum of indicator variograms were computed and modelled for each method of ranking. The individual indicators were then kriged with the parameters of the appropriate indicator variogram model to map the probability of encountering soil with the structure represented by that indicator. The model parameters of the SIVs for each ranking system were used with the data to krige the soil structure classes, and the results are compared with those for the individual indicators. The relations between maps of soil structure and selected wavebands from aerial photographs are examined as basis for planning surveys of soil structure. (C) 2007 Elsevier B.V. All rights reserved.

Mineral and geochemical characterization of a leptic aluandic soil and a thapto aluandic-ferralsol developed on trachytes in Mount Bambouto (Cameroon volcanic line)

Mineral and geochemical investigations were carried out on soil samples and fresh rock (trachytes) from two selected soil profiles (TM profile on leptic aluandic soils and TL profile on thapto aluandic-ferralsols) from Mount Bambouto to better understand geochemical processes and mineral paragenesis involved in the development of soils in this environment. In TM profile, the hydrated halloysites and goethite occur in the weathered saprolite boulders of BC horizon while dehydrated halloysite, gibbsite and goethite dominate the soils matrices of BC and A horizons. In TL profile, the dehydrated halloysites and goethite are the most abundant secondary minerals in the weathered saprolites of C and BC horizons while gibbsite, hematite and kaolinite occur in the soil matrices of BC, B and A horizons. The highest gibbsite content is in the platy nodules of B horizon. In both soil profiles, organo-metal complexes (most likely of AI and Fe) are present in the surface A horizon. Geochemically, between the fresh rock and the weathered saprolites in both soils, SiO2, K2O, CaO, Na2O and MgO contents decrease strongly while Fe2O3 and Al2O3 tend to accumulate. The molar ratio of SiO2/Al2O3 (Ki) and the sum of Ca, Mg, K and Na ions (TRB) also decreases abruptly between fresh rocks and the weathered saprolites, but increases significantly at the soil surface. The TM profile shows intense Al enrichment whereas the TL profile highlights enrichment in both AI and Fe as the weathering progresses upwards. Both soil profiles are enriched in Ni, Cu, Ba and Co and depleted in U, Th, Ta, Hf, Y, Sr, Pb, Zr and Zn relative to fresh rock. They also show a relatively low fractionation of the rare earth elements (REE: La, Nd, Sm, Eu, Tb, Yb and Lu), except for Ce which tends to be enriched in soils compared to CI chondrite. All these results give evidence of intense hydrolysis at soil deep in Mount Bambouto resulting in the formation of halloysite which progressively transforms into gibbsite and/or dehydrated halloysite. At the soil surface, the prominent pedogenetic process refers to andosolization with formation of organo-metal complexes. In TL profile, the presence of kaolinite in soil matrices BC and B horizons is consistent with ferralitization at soil deep. In conclusion, soil forming processes in Mount Bambouto are strongly influenced by local climate: (i) in the upper mountain (>2000 m), the fresh, misty and humid climate favors andosolization; whereas (ii) in the middle lands (1700-2000 m) with a relatively dry climate, both andosolization at the soil surface and ferralitization at soil deep act together. (C) 2009 Elsevier B.V. All rights reserved.

The potential for the use of agent-based models in ecotoxicology

This chapter introduces ABMs, their construction, and the pros and cons of their use. Although relatively new, agent-basedmodels (ABMs) have great potential for use in ecotoxicological research – their primary advantage being the realistic simulations that can be constructed and particularly their explicit handling of space and time in simulations. Examples are provided of their use in ecotoxicology primarily exemplified by different implementations of the ALMaSS system. These examples presented demonstrate how multiple stressors, landscape structure, details regarding toxicology, animal behavior, and socioeconomic effects can and should be taken into account when constructing simulations for risk assessment. Like ecological systems, in ABMs the behavior at the system level is not simply the mean of the component responses, but the sum of the often nonlinear interactions between components in the system; hence this modeling approach opens the door to implementing and testing much more realistic and holistic ecotoxicological models than are currently used.

The use of image analysis to estimate population growth rate in Daphnia magna

1. Population growth rate (PGR) is central to the theory of population ecology and is crucial for projecting population trends in conservation biology, pest management and wildlife harvesting. Furthermore, PGR is increasingly used to assess the effects of stressors. Image analysis that can automatically count and measure photographed individuals offers a potential methodology for estimating PGR. 2. This study evaluated two ways in which the PGR of Daphnia magna, exposed to different stressors, can be estimated using an image analysis system. The first method estimated PGR as the ratio of counts of individuals obtained at two different times, while the second method estimated PGR as the ratio of population sizes at two different times, where size is measured by the sum of the individuals' surface areas, i.e. total population surface area. This method is attractive if surface area is correlated with reproductive value (RV), as it is for D. magna, because of the theoretical result that PGR is the rate at which the population RV increases. 3. The image analysis system proved reliable and reproducible in counting populations of up to 440 individuals in 5 L of water. Image counts correlated well with manual counts but with a systematic underestimate of about 30%. This does not affect accuracy when estimating PGR as the ratio of two counts. Area estimates of PGR correlated well with count estimates, but were systematically higher, possibly reflecting their greater accuracy in the study situation. 4. Analysis of relevant scenarios suggested the correlation between RV and body size will generally be good for organisms in which fecundity correlates with body size. In these circumstances, area estimation of PGR is theoretically better than count estimation. 5. Synthesis and applications. There are both theoretical and practical advantages to area estimation of population growth rate when individuals' reproductive values are consistently well correlated with their surface areas. Because stressors may affect both the number and quality of individuals, area estimation of population growth rate should improve the accuracy of predicting stress impacts at the population level.

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.

Infection of bryozoans by Tetracapsuloides bryosalmonae at sites endemic for salmonid proliferative kidney disease

Laboratory-reared colonies of the bryozoans Fredericella sultana and Plumatella fungosa were placed upstream of 2 fish farms endemic for salmonid proliferative kidney disease (PKD) to assess rates of infection of bryozoans by Tetra caps uloides bryosalmonae, the causative agent of PKD. Colonies were deployed in the field for 8 trial periods of 2 wk each throughout the summer of 2001. Following each trial, bryozoan colonies were maintained in laboratory culture for 28 d and were regularly monitored for infection by searching for sac stages of T bryosalmonae. Infections were never identified by observations of sac stages, however positive PCR results and sequencing of cultured material confirmed that cryptic infections were present in colonies of both species deployed at one site. The possibility that PCR results reflected contamination of surfaces of bryozoans can be excluded, given the short period of spore viability of T bryosalmonae. Highest rates of infection occurred when 4 of 23 colonies of F sultana and 1 of 12 colonies of P. fungosa were infected during the period 10 to 24 July. No infections were detected from mid-August to late October at this site. None of the colonies at the other site became infected throughout the period of study. Our data provide the first estimates of infection rates of bryozoans by T bryosalmonae. Additionally, they provide evidence that a cryptic stage can be maintained within bryozoan hosts for a period of 4 to 6 wk.

Hydroponic isotope labelling of entire plants (HILEP) for quantitative plant proteomics; an oxidative stress case study

Hydroponic isotope labelling of entire plants (HILEP) is a cost-effective method enabling metabolic labelling of whole and mature plants with a stable isotope such as N-15. By utilising hydroponic media that contain N-15 inorganic salts as the sole nitrogen source, near to 100% N-15-labelling of proteins can be achieved. In this study, it is shown that HILEP, in combination with mass spectrometry, is suitable for relative protein quantitation of seven week-old Arabidopsis plants submitted to oxidative stress. Protein extracts from pooled N-14- and N-15-hydroponically grown plants were fractionated by SDS-PAGE, digested and analysed by liquid chromatography electrospray ionisation tandem mass spectrometry (LC-ESI-MS/MS). Proteins were identified and the spectra of N-14/N-15 peptide pairs were extracted using their m/z chromatographic retention time, isotopic distributions, and the m/z difference between the N-14 and N-15 peptides. Relative amounts were calculated as the ratio of the sum of the peak areas of the two distinct N-14 and N-15 peptide isotope envelopes. Using Mascot and the open source trans-proteomic pipeline (TPP), the data processing was automated for global proteome quantitation down to the isoform level by extracting isoform specific peptides. With this combination of metabolic labelling and mass spectrometry it was possible to show differential protein expression in the apoplast of plants submitted to oxidative stress. Moreover, it was possible to discriminate between differentially expressed isoforms belonging to the same protein family, such as isoforms of xylanases and pathogen-related glucanases (PR 2). (C) 2008 Elsevier Ltd. All rights reserved.

LIPGENE food-exchange model for alteration of dietary fat quantity and quality in free-living participants from eight European countries

Controlled human intervention trials are required to confirm the hypothesis that dietary fat quality may influence insulin action. The aim was to develop a food-exchange model, suitable for use in free-living volunteers, to investigate the effects of four experimental diets distinct in fat quantity and quality: high SFA (HSFA); high MUFA (HMUFA) and two low-fat (LF) diets, one supplemented with 1.24g EPA and DHA/d (LFn-3). A theoretical food-exchange model was developed. The average quantity of exchangeable fat was calculated as the sum of fat provided by added fats (spreads and oils), milk, cheese, biscuits, cakes, buns and pastries using data from the National Diet and Nutrition Survey of UK adults. Most of the exchangeable fat was replaced by specifically designed study foods. Also critical to the model was the use of carbohydrate exchanges to ensure the diets were isoenergetic. Volunteers from eight centres across Europe completed the dietary intervention. Results indicated that compositional targets were largely achieved with significant differences in fat quantity between the high-fat diets (39.9 (SEM 0.6) and 38.9 (SEM 0.51) percentage energy (%E) from fat for the HSFA and HMUFA diets respectively) and the low-fat diets (29.6 (SEM 0.6) and 29.1 (SEM 0.5) %E from fat for the LF and LFn-3 diets respectively) and fat quality (17.5 (SEM 0.3) and 10.4 (SEM 0.2) %E front SFA and 12.7 (SEM 0.3) and 18.7 (SEM 0.4) %E MUFA for the HSFA and HMUFA diets respectively). In conclusion, a robust, flexible food-exchange model was developed and implemented successfully in the LIPGENE dietary intervention trial.

Profiling of composition and metabolic activities of the colonic microflora of growing pigs fed diets supplemented with prebiotic oligosaccharides

It is evident that quantitative information on different microbial groups and their contribution in terms of activity in the gastrointestinal (GI) tract of humans and animals is required in order to formulate functional diets targeting improved gut function and host health. In this work, quantitative information on levels and spatial distributions of Bacteroides spp, Eubacterium spp, Clostridium spp, Escherichia coli, Bifidobacterium spp and Lactobacillus/Enterococcus spp. along the porcine large intestine was investigated using 16S rRNA targeted probes and fluorescent in situ hybridisation (FISH). Caecum, ascending colon (AC) and rectum luminal digesta from three groups of individually housed growing pigs fed either a corn-soybean basal diet (CON diet) or a prebiotic diet containing 10 g/kg oligofructose (FOS diet) or trans-galactooligosaccharides (TOS diet) at the expense of cornstarch were analysed. DAPI staining was used to enumerate total number of cells in the samples. Populations of total cells, Bacteroides, Eubacterium, Clostridium and Bifidobacterium, declined significantly (P < 0.05) from caecum to rectum, and were not affected by dietary treatments. Populations of Lactobacillus/ Enterococcus and E coli did not differ throughout the large intestine. The relative percent (%) contribution of each bacterial group to the total cell count did not differ between caecum and rectum, with the exception of Eubacterium that was higher in the AC digesta. FISH analysis showed that the sum of all bacterial groups made up a small percentage of the total cells, which was 12.4%, 21.8% and 10.3% in caecum, AC and rectum, respectively. This supports the view that in swine, the diversity of GI microflora might be higher compared to other species. In terms of microflora metabolic activity, the substantially higher numerical trends seen in FOS and TOS treatments regarding total volatile fatty acid, acetate concentrations and glycolytic activities, it could be postulated that FOS and TOS promoted saccharolytic activities in the porcine colon. (c) 2006 Elsevier Ltd. All rights reserved.

Lack of association between central adiposity and lipaemia in UK Sikh men

OBJECTIVE: To determine whether the positive statistical associations between measures of total and regional adiposity and measures of glucose, insulin and triacylglycerol ( TAG) metabolism reported in Caucasian men, are also observed in UK Sikhs. DESIGN: A matched cross-sectional study in which each volunteer provided a blood sample after a 12-h overnight fast and had anthropometric measurements taken. SUBJECTS: A total of 55 healthy Caucasian and 55 healthy UK Sikh men were recruited. The Caucasian and Sikh men were matched for age ( 48.7 +/- 10.9 and 48.3 +/- 10.0 y, respectively) and body mass index (BMI) ( 26.1 +/- 2.8 and 26.3 +/- 3.2 kg/m(2), respectively). MEASUREMENTS: Anthropometric measurements were performed to assess total and regional fat depots. The concentrations of plasma total cholesterol, high-density cholesterol (HDL- C), low-density cholesterol (LDL-C) and small dense LDL (LDL3), TAG, glucose, fasting insulin (ins) and nonesterified fatty acids (NEFA) were analysed in fasted plasma. Surrogate measures of insulin resistance (HOMA-IR) and insulin sensitivity (RQUICKI) were calculated from insulin and glucose (HOMA-IR) and insulin, glucose and NEFA ( RQUICKI) measurements. RESULTS: The Sikh men had significantly higher body fat, with the sum of the four skinfold measurements (Ssk) ( P = 0.0001) and subscapular skinfold value (P = 0.009) higher compared with the Caucasian men. The Sikh volunteers also had characteristics of the metabolic syndrome: lower HDL-C (P = 0.07), higher TAG (P = 0.004), higher % LDL3 (P = 0.0001) and insulin resistance (P = 0.05). Both ethnic groups demonstrated positive correlations between insulin and waist circumference (Caucasian: r = 0.661, P = 0.0001; Sikh: r = 0.477, P = 0.0001). The Caucasian men also demonstrated significant positive correlations between central adiposity (r = 0.275, P = 0.04), other measures of adiposity (BMI and suprailiac skinfold) and plasma TAG, whereas the Sikh men showed no correlation for central adiposity (r = 0.019, ns) and TAG with a trend to a negative relationship between other measures ( Ssk and suprailiac) which reached near significance for subscapular skinfold and TAG (r = - 0.246, P = 0.007). The expected positive association between insulin and TAG was observed in the Caucasian men (r = 0.318, P = 0.04) but not in the Sikh men (r = 0.011, ns). CONCLUSIONS: In the Caucasian men, the expected positive association between plasma TAG and centralized body fat was observed. However, a lack of association between centralized, or any other measure of adiposity, and plasma TAG was observed in the matched Sikh men, although both ethnic groups showed the positive association between centralized body fat and insulin resistance, which was less strong for Sikhs. These findings in the Sikh men were not consistent with the hypothesis that there is a clear causal relationship between body fat and its distribution, insulin resistance, and lipid abnormalities associated with the metabolic syndrome, in this ethnic group.

The automatic identification and prioritisation of criminal networks from police crime data

The identification of criminal networks is not a routine exploratory process within the current practice of the law enforcement authorities; rather it is triggered by specific evidence of criminal activity being investigated. A network is identified when a criminal comes to notice and any associates who could also be potentially implicated would need to be identified if only to be eliminated from the enquiries as suspects or witnesses as well as to prevent and/or detect crime. However, an identified network may not be the one causing most harm in a given area.. This paper identifies a methodology to identify all of the criminal networks that are present within a Law Enforcement Area, and, prioritises those that are causing most harm to the community. Each crime is allocated a score based on its crime type and how recently the crime was committed; the network score, which can be used as decision support to help prioritise it for law enforcement purposes, is the sum of the individual crime scores.

Integrating control systems defined on the frame bundles of the space forms

This paper considers left-invariant control systems defined on the orthonormal frame bundles of simply connected manifolds of constant sectional curvature, namely the space forms Euclidean space E-3, the sphere S-3 and Hyperboloid H-3 with the corresponding frame bundles equal to the Euclidean group of motions SE(3), the rotation group SO(4) and the Lorentz group SO(1, 3). Orthonormal frame bundles of space forms coincide with their isometry groups and therefore the focus shifts to left-invariant control systems defined on Lie groups. In this paper a method for integrating these systems is given where the controls are time-independent. In the Euclidean case the elements of the Lie algebra se(3) are often referred to as twists. For constant twist motions, the corresponding curves g(t) is an element of SE(3) are known as screw motions, given in closed form by using the well known Rodrigues' formula. However, this formula is only applicable to the Euclidean case. This paper gives a method for computing the non-Euclidean screw motions in closed form. This involves decoupling the system into two lower dimensional systems using the double cover properties of Lie groups, then the lower dimensional systems are solved explicitly in closed form.

Parallel Monte Carlo approach for integration of the rendering equation

This paper is addressed to the numerical solving of the rendering equation in realistic image creation. The rendering equation is integral equation describing the light propagation in a scene accordingly to a given illumination model. The used illumination model determines the kernel of the equation under consideration. Nowadays, widely used are the Monte Carlo methods for solving the rendering equation in order to create photorealistic images. In this work we consider the Monte Carlo solving of the rendering equation in the context of the parallel sampling scheme for hemisphere. Our aim is to apply this sampling scheme to stratified Monte Carlo integration method for parallel solving of the rendering equation. The domain for integration of the rendering equation is a hemisphere. We divide the hemispherical domain into a number of equal sub-domains of orthogonal spherical triangles. This domain partitioning allows to solve the rendering equation in parallel. It is known that the Neumann series represent the solution of the integral equation as a infinity sum of integrals. We approximate this sum with a desired truncation error (systematic error) receiving the fixed number of iteration. Then the rendering equation is solved iteratively using Monte Carlo approach. At each iteration we solve multi-dimensional integrals using uniform hemisphere partitioning scheme. An estimate of the rate of convergence is obtained using the stratified Monte Carlo method. This domain partitioning allows easy parallel realization and leads to convergence improvement of the Monte Carlo method. The high performance and Grid computing of the corresponding Monte Carlo scheme are discussed.

On the spectra of certain integro-differential-delay problems with applications in neurodynamics

We investigate the spectrum of certain integro-differential-delay equations (IDDEs) which arise naturally within spatially distributed, nonlocal, pattern formation problems. Our approach is based on the reformulation of the relevant dispersion relations with the use of the Lambert function. As a particular application of this approach, we consider the case of the Amari delay neural field equation which describes the local activity of a population of neurons taking into consideration the finite propagation speed of the electric signal. We show that if the kernel appearing in this equation is symmetric around some point a= 0 or consists of a sum of such terms, then the relevant dispersion relation yields spectra with an infinite number of branches, as opposed to finite sets of eigenvalues considered in previous works. Also, in earlier works the focus has been on the most rightward part of the spectrum and the possibility of an instability driven pattern formation. Here, we numerically survey the structure of the entire spectra and argue that a detailed knowledge of this structure is important within neurodynamical applications. Indeed, the Amari IDDE acts as a filter with the ability to recognise and respond whenever it is excited in such a way so as to resonate with one of its rightward modes, thereby amplifying such inputs and dampening others. Finally, we discuss how these results can be generalised to the case of systems of IDDEs.

The full-spectrum correlated k method for longwave atmospheric radiative transfer: treatment of gaseous overlap

The correlated k-distribution (CKD) method is widely used in the radiative transfer schemes of atmospheric models, and involves dividing the spectrum into a number of bands and then reordering the gaseous absorption coefficients within each one. The fluxes and heating rates for each band may then be computed by discretizing the reordered spectrum into of order 10 quadrature points per major gas, and performing a pseudo-monochromatic radiation calculation for each point. In this paper it is first argued that for clear-sky longwave calculations, sufficient accuracy for most applications can be achieved without the need for bands: reordering may be performed on the entire longwave spectrum. The resulting full-spectrum correlated k (FSCK) method requires significantly fewer pseudo-monochromatic calculations than standard CKD to achieve a given accuracy. The concept is first demonstrated by comparing with line-by-line calculations for an atmosphere containing only water vapor, in which it is shown that the accuracy of heating-rate calculations improves approximately in proportion to the square of the number of quadrature points. For more than around 20 points, the root-mean-squared error flattens out at around 0.015 K d−1 due to the imperfect rank correlation of absorption spectra at different pressures in the profile. The spectral overlap of m different gases is treated by considering an m-dimensional hypercube where each axis corresponds to the reordered spectrum of one of the gases. This hypercube is then divided up into a number of volumes, each approximated by a single quadrature point, such that the total number of quadrature points is slightly fewer than the sum of the number that would be required to treat each of the gases separately. The gaseous absorptions for each quadrature point are optimized such they minimize a cost function expressing the deviation of the heating rates and fluxes calculated by the FSCK method from line-by-line calculations for a number of training profiles. This approach is validated for atmospheres containing water vapor, carbon dioxide and ozone, in which it is found that in the troposphere and most of the stratosphere, heating-rate errors of less than 0.2 K d−1 can be achieved using a total of 23 quadrature points, decreasing to less than 0.1 K d−1 for 32 quadrature points. It would be relatively straightforward to extend the method to include other gases.