Earthworm distribution and abundance predicted by a process-based model

Earthworms are significant ecosystem engineers and are an important component of the diet of many vertebrates and invertebrates, so the ability to predict their distribution and abundance would have wide application in ecology, conservation and land management. Earthworm viability is known to be affected by the availability and quality of food resources, soil water conditions and temperature, but has not yet been modelled mechanistically to link effects on individuals to field population responses. Here we present a novel model capable of predicting the effects of land management and environmental conditions on the distribution and abundance of Aporrectodea caliginosa, the dominant earthworm species in agroecosystems. Our process-based approach uses individual based modelling (IBM), in which each individual has its own energy budget. Individual earthworm energy budgets follow established principles of physiological ecology and are parameterised for A. caliginosa from experimental measurements under optimal conditions. Under suboptimal conditions (e.g. food limitation, low soil temperatures and water contents) reproduction is prioritised over growth. Good model agreement to independent laboratory data on individual cocoon production and growth of body mass, under variable feeding and temperature conditions support our representation of A. caliginosa physiology through energy budgets. Our mechanistic model is able to accurately predict A. caliginosa distribution and abundance in spatially heterogeneous soil profiles representative of field study conditions. Essential here is the explicit modelling of earthworm behaviour in the soil profile. Local earthworm movement responds to a trade-off between food availability and soil water conditions, and this determines the spatiotemporal distribution of the population in the soil profile. Importantly, multiple environmental variables can be manipulated simultaneously in the model to explore earthworm population exposure and effects to combinations of stressors. Potential applications include prediction of the population-level effects of pesticides and changes in soil management e.g. conservation tillage and climate change.

Validation of the Deardorff model for estimating energy balance components for a sugarcane crop

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Enzyme Optimization for Lignocellulose Hydrolysis using Mechanistic Modeling

A novel mechanistic model for the saccharification of cellulose and hemicellulose is utilized to predict the products of hydrolysis over a range of enzyme loadings and times. The mechanistic model considers the morphology of the substrate and the kinetics of enzymes to optimize enzyme concentrations for the enzymatic hydrolysis of cellulose and hemicellulose simultaneously. Substrates are modeled based on their fraction of accessible sites, glucan content, xylan content, and degree of polymerizations. This enzyme optimization model takes into account the kinetics of six core enzymes for lignocellulose hydrolysis: endoglucanase I (EG1), cellobiohydrolase I (CBH1), cellobiohydrolase II (CBH2), and endo-xylanase (EX) from Trichoderma reesei; β-glucosidase (BG), and β-xylosidase (BX) from Aspergillus niger. The model employs the synergistic action of these enzymes to predict optimum enzyme concentrations for hydrolysis of Avicel and ammonia fiber explosion (AFEX) pretreated corn stover. Glucan, glucan + xylan, glucose and glucose + xylose conversion predictions are given over a range of mass fractions of enzymes, and a range of enzyme loadings. Simulation results are compared with optimizations using statistically designed experiments. BG and BX are modeled in solution at later time points to predict the effect on glucose conversion and xylose conversion.

Slaving and release in co-infection control

Background Animal and human infection with multiple parasite species is the norm rather than the exception, and empirical studies and animal models have provided evidence for a diverse range of interactions among parasites. We demonstrate how an optimal control strategy should be tailored to the pathogen community and tempered by species-level knowledge of drug sensitivity with use of a simple epidemiological model of gastro-intestinal nematodes. Methods We construct a fully mechanistic model of macroparasite co-infection and use it to explore a range of control scenarios involving chemotherapy as well as improvements to sanitation. Results Scenarios are presented whereby control not only releases a more resistant parasite from antagonistic interactions, but risks increasing co-infection rates, exacerbating the burden of disease. In contrast, synergisms between species result in their becoming epidemiologically slaved within hosts, presenting a novel opportunity for controlling drug resistant parasites by targeting co-circulating species. Conclusions Understanding the effects on control of multi-parasite species interactions, and vice versa, is of increasing urgency in the advent of integrated mass intervention programmes.

Modelling predicts positive and negative interactions between three Australian tropical tree species in monoculture and binary mixture

Computer modelling promises to be an important tool for analysing and predicting interactions between trees within mixed species forest plantations. This study explored the use of an individual-based mechanistic model as a predictive tool for designing mixed species plantations of Australian tropical trees. The `spatially explicit individually based-forest simulator' (SeXI-FS) modelling system was used to describe the spatial interaction of individual tree crowns within a binary mixed-species experiment. The three-dimensional model was developed and verified with field data from three forest tree species grown in tropical Australia. The model predicted the interactions within monocultures and binary mixtures of Flindersia brayleyana, Eucalyptus pellita and Elaeocarpus grandis, accounting for an average of 42% of the growth variation exhibited by species in different treatments. The model requires only structural dimensions and shade tolerance as species parameters. By modelling interactions in existing tree mixtures, the model predicted both increases and reductions in the growth of mixtures (up to +/-50% of stem volume at 7 years) compared to monocultures. This modelling approach may be useful for designing mixed tree plantations.

Development of novel grain morphology during hot extrusion of magnesium AZ21 alloy

Microstructure and microtexture evolution during static annealing of a hot-extruded AZ21 magnesium alloy was studied. Apart from fine recrystallized equiaxed grains and large elongated deformed grains, a new third kind of abnormal grains that are stacked one after the other in a row parallel to the extrusion direction were observed. The crystallographic misorientation inside these grains was similar to that of the fine recrystallized grains. The large elongated grains exhibited significant in-grain misorientation. A self-consistent mechanistic model was developed to describe the formation of these grain morphologies during dynamic recrystallization (DRX). The texture of pre-extruded material, although lost in DRX, leaves a unique signature which manifests itself in the form of these grain morphologies. The origin of abnormal stacked grains was associated with slow nucleation in pre-extruded grains of a certain orientation. Further annealing resulted in large secondary recrystallized grains with occasional extension twins. (c) 2009 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Dye sensitized visible light degradation of phenolic compounds

The present research work reports the eosin Y (EY) and fluorescein (FL) sensitized visible light degradation of phenol, 4-chlorophenol (CP), 2,4-dichlorophenol (DCP) and 2,4,6-trichlorophenol (TCP) using combustion synthesized nano-TiO2 (CS TiO2). The rate of degradation of the phenolic compounds was higher in the presence of EY/CS TiO2 compared to FL/CS TiO2 system. A detailed mechanism of sensitized degradation was proposed and a mechanistic model for the rate of degradation of the phenolic compound was derived using the pyramidal network reduction technique. It was found that at low initial dye concentrations, the rate of degradation of the phenolic compound is first order in the concentration of the dye, while at high initial dye concentrations, the rate is first order in the concentration of the phenolic compound. The order of degradation of the different phenolic compounds follows: CP > TCP > DCP > phenol. The different phenolic and dye intermediates that were formed during the degradation were identified by liquid chromatography-mass spectrometry (LC-MS) and the most probable pathway of degradation is proposed. (C) 2010 Elsevier B.V. All rights reserved.

Field emission properties of carbon nanotube arrays with defects and impurities

Field emission from carbon nanotubes (CNTs) in the form of arrays or thin films give rise to several strongly correlated process of electromechanical interaction and degradation. Such processes are mainly due to (1) electron-phonon interaction (2) electromechanical force field leading to stretching of CNTs (3) ballistic transport induced thermal spikes, coupled with high dynamic stress, leading to degradation of emission performance at the device scale. Fairly detailed physics based models of CNTs considering the aspects (1) and (2) above have already been developed by these authors, and numerical results indicate good agreement with experimental results. What is missing in such a system level modeling approach is the incorporation of structural defects and vacancies or charge impurities. This is a practical and important problem due to the fact that degradation of field emission performance is indeed observed in experimental I-V curves. What is not clear from these experiments is whether such degradation in the I-V response is due to dynamic reorientation of the CNTs or due to the defects or due to both of these effects combined. Non-equilibrium Green’s function based simulations using a tight-binding Hamiltonian for single CNT segment show up the localization of carrier density at various locations of the CNTs. About 11% decrease in the drive current with steady difference in the drain current in the range of 0.2-0.4V of the gate voltage was reported in literature when negative charge impurity was introduced at various locations of the CNT over a length of ~20nm. In the context of field emission from CNT tips, a simplistic estimate of defects have been introduced by a correction factor in the Fowler-Nordheim formulae. However, a more detailed physics based treatment is required, while at the same time the device-scale simulation is necessary. The novelty of our present approach is the following. We employ a concept of effective stiffness degradation for segments of CNTs, which is due to structural defects, and subsequently, we incorporate the vacancy defects and charge impurity effects in the Green’s function based approach. Field emission induced current-voltage characteristics of a vertically aligned CNT array on a Cu-Cr substrate is then simulated using a detailed nonlinear mechanistic model of CNTs coupled with quantum hydrodynamics. An array of 10 vertically aligned and each 12 m long CNTs is considered for the device scale analysis. Defect regions are introduced randomly over the CNT length. The result shows the decrease in the longitudinal strain due to defects. Contrary to the expected influence of purely mechanical degradation, this result indicates that the charge impurity and hence weaker transport can lead to a different electromechanical force field, which ultimately can reduce the strain. However, there could be significant fluctuation in such strain field due to electron-phonon coupling. The effect of such fluctuations (with defects) is clearly evident in the field emission current history. The average current also decreases significantly due to such defects.

CONFORMATIONAL FEATURES OF SIGMA FACTOR/ANTI-SIGMA COMPLEXES

The association of a factors with the RNA polymerase dictates the expression profile of a bacterial cell. Major changes to the transcription profile are achieved by the use of multiple sigma factors that confer distinct promoter selectivity to the holoenzyme. The cellular concentration of a sigma factor is regulated by diverse mechanisms involving transcription, translation and post-translational events. The number of sigma factors varies substantially across bacteria. The diversity in the interactions between sigma factors also vary-ranging from collaboration, competition or partial redundancy in some cellular or environmental contexts. These interactions can be rationalized by a mechanistic model referred to as the partitioning of a space model of bacterial transcription. The structural similarity between different sigma/anti-sigma complexes despite poor sequence conservation and cellular localization reveals an elegant route to incorporate diverse regulatory mechanisms within a structurally conserved scaffold. These features are described here with a focus on sigma/anti-sigma complexes from Mycobacterium tuberculosis. In particular, we discuss recent data on the conditional regulation of sigma/anti-sigma factor interactions. Specific stages of M. tuberculosis infection, such as the latent phase, as well as the remarkable adaptability of this pathogen to diverse environmental conditions can be rationalized by the synchronized action of different a factors.

Influence of the AgrC-AgrA Complex on the Response Time of Staphylococcus aureus Quorum Sensing

The Staphylococcus aureus agr quorum-sensing system plays a major role in the transition from the persistent to the virulent phenotype. S. aureus agr type I to IV strains are characterized by mutations in the sensor domain of the histidine kinase AgrC and differences in the sequences of the secreted autoinducing peptides (AIP). Here we demonstrate that interactions between the cytosolic domain of AgrC (AgrC(Cyto)) and the response regulator domain of AgrA (AgrA(RR)) dictate the spontaneity of the cellular response to AIP stimuli. The crystal structure of AgrC(Cyto) provided a basis for a mechanistic model of AgrC-AgrA interactions. This model enabled an analysis of the biochemical and biophysical parameters of AgrC-AgrA interactions in the context of the conformational features of the AgrC-AgrA complex. This analysis revealed distinct sequence and conformational features that determine the affinity, specificity, and kinetics of the phosphotransfer reaction. This step, which governs the response time for transcriptional reengineering triggered by an AIP stimulus, is independent of the agr type and similar for agonist and antagonist stimuli. These experimental data could serve as a basis on which to validate simulations of the quorum-sensing response and for strategies that employ the agr quorum-sensing system to combat biofilm formation in S. aureus infections.

Fingerprinting environmental conditions and related stress using stable isotopic composition of rice (Oryza sativa L.) grain organic matter

Sea level rise (SLR) is a primary factor responsible for inundation of low-lying coastal regions across the world, which in turn governs the agricultural productivity. In this study, rice (Oryza sativa L.) cultivated seasonally in the Kuttanad Wetland, a SLR prone region on the southwest coast of India, were analysed for oxygen, hydrogen and carbon isotopic ratios (delta O-18, delta H-2 and delta C-13) to distinguish the seasonal environmental conditions prevalent during rice cultivation. The region receives high rainfall during the wet season which promotes large supply of fresh water to the local water bodies via the rivers. In contrast, during the dry season reduced river discharge favours sea water incursion which adversely affects the rice cultivation. The water for rice cultivation is derived from regional water bodies that are characterised by seasonal salinity variation which co-varies with the delta O-18 and delta H-2 values. Rice cultivated during the wet and the dry season bears the isotopic imprints of this water. We explored the utility of a mechanistic model to quantify the contribution of two prominent factors, namely relative humidity and source water composition in governing the seasonal variation in oxygen isotopic composition of rice grain OM. delta C-13 values of rice grain OM were used to deduce the stress level by estimating the intrinsic water use efficiency (WUEi) of the crop during the two seasons. 1.3 times higher WUE, was exhibited by the same genotype during the dry season. The approach can be extended to other low lying coastal agro-ecosystems to infer the growth conditions of cultivated crops and can further be utilised for retrieving paleo-environmental information from well preserved archaeological plant remains. (c) 2015 Elsevier Ltd. All rights reserved.

Kinetic investigations of heterogeneously catalyzed reactions on the Ru(001) and Pt(110)-(1x2) surfaces

The initial probabilities of activated, dissociative chemisorption of methane and ethane on Pt(110)-(1 x 2) have been measured. The surface temperature was varied from 450 to 900 K with the reactant gas temperature constant at 300 K. Under these conditions, we probe the kinetics of dissociation via trapping-mediated (as opposed to 'direct') mechanism. It was found that the probabilities of dissociation of both methane and ethane were strong functions of the surface temperature with an apparent activation energies of 14.4 kcal/mol for methane and 2.8 kcal/mol for ethane, which implys that the methane and ethane molecules have fully accommodated to the surface temperature. Kinetic isotope effects were observed for both reactions, indicating that the C-H bond cleavage was involved in the rate-limiting step. A mechanistic model based on the trapping-mediated mechanism is used to explain the observed kinetic behavior. The activation energies for C-H bond dissociation of the thermally accommodated methane and ethane on the surface extracted from the model are 18.4 and 10.3 kcal/mol, respectively.

The studies of the catalytic decomposition of formic acid on the Ru(001) surface with thermal desorption mass spectrometry following the adsorption of DCOOH and HCOOH on the surface at 130 and 310 K are described. Formic acid (DCOOH) chemisorbs dissociatively on the surface via both the cleavage of its O-H bond to form a formate and a hydrogen adatom, and the cleavage of its C-O bond to form a carbon monoxide, a deuterium adatom and an hydroxyl (OH). The former is the predominant reaction. The rate of desorption of carbon dioxide is a direct measure of the kinetics of decomposition of the surface formate. It is characterized by a kinetic isotope effect, an increasingly narrow FWHM, and an upward shift in peak temperature with Ɵ_T, the coverage of the dissociatively adsorbed formic acid. The FWHM and the peak temperature change from 18 K and 326 K at Ɵ_T = 0.04 to 8 K and 395 K at Ɵ_T = 0.89. The increase in the apparent activation energy of the C-D bond cleavage is largely a result of self-poisoning by the formate, the presence of which on the surface alters the electronic properties of the surface such that the activation energy of the decomposition of formate is increased. The variation of the activation energy for carbon dioxide formation with Ɵ_T accounts for the observed sharp carbon dioxide peak. The coverage of surface formate can be adjusted over a relatively wide range so that the activation energy for C-D bond cleavage in the case of DCOOH can be adjusted to be below, approximately equal to, or well above the activation energy for the recombinative desorption of the deuterium adatoms. Accordingly, the desorption of deuterium was observed to be governed completely by the desorption kinetics of the deuterium adatoms at low Ɵ_T, jointly by the kinetics of deuterium desorption and C-D bond cleavage at intermediate Ɵ_T, and solely by the kinetics of C-D bond cleavage at high Ɵ_T. The overall branching ratio of the formate to carbon dioxide and carbon monoxide is approximately unity, regardless the initial coverage Ɵ_T, even though the activation energy for the production of carbon dioxide varies with Ɵ_T. The desorption of water, which implies C-O bond cleavage of the formate, appears at approximately the same temperature as that of carbon dioxide. These observations suggest that the cleavage of the C-D bond and that of the C-O bond of two surface formates are coupled, possibly via the formation of a short-lived surface complex that is the precursor to to the decomposition.

The measurement of steady-state rate is demonstrated here to be valuable in determining kinetics associated with short-lived, molecularly adsorbed precursor to further reactions on the surface, by determining the kinetic parameters of the molecular precursor of formaldehyde to its dissociation on the Pt(110)-(1 x 2) surface.

Overlayers of nitrogen adatoms on Ru(001) have been characterized both by thermal desorption mass spectrometry and low-energy electron diffraction, as well as chemically via the postadsorption and desorption of ammonia and carbon monoxide.

The nitrogen-adatom overlayer was prepared by decomposing ammonia thermally on the surface at a pressure of 2.8 x 10^(-6) Torr and a temperature of 480 K. The saturated overlayer prepared under these conditions has associated with it a (√247/10 x √247/10)R22.7° LEED pattern, has two peaks in its thermal desorption spectrum, and has a fractional surface coverage of 0.40. Annealing the overlayer to approximately 535 K results in a rather sharp (√3 x √3)R30° LEED pattern with an associated fractional surface coverage of one-third. Annealing the overlayer further to 620 K results in the disappearance of the low-temperature thermal desorption peak and the appearance of a rather fuzzy p(2x2) LEED pattern with an associated fractional surface coverage of approximately one-fourth. In the low coverage limit, the presence of the (√3 x √3)R30° N overlayer alters the surface in such a way that the binding energy of ammonia is increased by 20% relative to the clean surface, whereas that of carbon monoxide is reduced by 15%.

A general methodology for the indirect relative determination of the absolute fractional surface coverages has been developed and was utilized to determine the saturation fractional coverage of hydrogen on Ru(001). Formaldehyde was employed as a bridge to lead us from the known reference point of the saturation fractional coverage of carbon monoxide to unknown reference point of the fractional coverage of hydrogen on Ru(001), which is then used to determine accurately the saturation fractional coverage of hydrogen. We find that ƟSAT/H = 1.02 (±0.05), i.e., the surface stoichiometry is Ru : H = 1 : 1. The relative nature of the method, which cancels systematic errors, together with the utilization of a glass envelope around the mass spectrometer, which reduces spurious contributions in the thermal desorption spectra, results in high accuracy in the determination of absolute fractional coverages.

Mechanics of deformable glacier beds

My focus in this thesis is to contribute to a more thorough understanding of the mechanics of ice and deformable glacier beds. Glaciers flow under their own weight through a combination of deformation within the ice column and basal slip, which involves both sliding along and deformation within the bed. Deformable beds, which are made up of unfrozen sediment, are prevalent in nature and are often the primary contributors to ice flow wherever they are found. Their granular nature imbues them with unique mechanical properties that depend on the granular structure and hydrological properties of the bed. Despite their importance for understanding glacier flow and the response of glaciers to changing climate, the mechanics of deformable glacier beds are not well understood.

Our general approach to understanding the mechanics of bed deformation and their effect on glacier flow is to acquire synoptic observations of ice surface velocities and their changes over time and to use those observations to infer the mechanical properties of the bed. We focus on areas where changes in ice flow over time are due to known environmental forcings and where the processes of interest are largely isolated from other effects. To make this approach viable, we further develop observational methods that involve the use of mapping radar systems. Chapters 2 and 5 focus largely on the development of these methods and analysis of results from ice caps in central Iceland and an ice stream in West Antarctica. In Chapter 3, we use these observations to constrain numerical ice flow models in order to study the mechanics of the bed and the ice itself. We show that the bed in an Iceland ice cap deforms plastically and we derive an original mechanistic model of ice flow over plastically deforming beds that incorporates changes in bed strength caused by meltwater flux from the surface. Expanding on this work in Chapter 4, we develop a more detailed mechanistic model for till-covered beds that helps explain the mechanisms that cause some glaciers to surge quasi-periodically. In Antarctica, we observe and analyze the mechanisms that allow ocean tidal variations to modulate ice stream flow tens of kilometers inland. We find that the ice stream margins are significantly weakened immediately upstream of the area where ice begins to float and that this weakening likely allows changes in stress over the floating ice to propagate through the ice column.

Analysis of high wheat yields in Northwest China

Grain yields of over 14 Mg ha(-1) were reported in 1978 for spring wheat (Triticum aestivum L.) grown in Northwest China. Understanding the circumstances under which this record yield was achieved may be useful in defining the key factors that lead to high grain yields and in determining the limits to wheat yield. A relatively simple, mechanistic model was used in an effort to simulate the record yield. The model was used as a framework in which various crop traits could be adjusted to match the observed crop growth. The weather that was characterized by cool temperatures and high levels of solar radiation, proved to be especially important in allowing a full-season crop to achieve record yields. Variables defining plant development in the model also had to be set to describe the high yielding cultivar grown in China. Leaf development was defined by the length of a phyllochron, which was set equal to 78 TU (thermal units, base temperature equal to 0 degrees C) based on independent data. The description of grain fill had to be defined to match simulation results with the observations. Two variables, length of the grain-fill period and the grain growth rate, were set in response to the unique traits of this cultivar and the low temperatures during grain development. These simulations led to important suggestions for examining the interaction between cool temperature regimes and developmental traits of wheat cultivars. (C) 1997 Published by Elsevier Science Ltd.

Modelling the Stoichiometric Regulation of C-Rich Toxins in Marine Dinoflagellates

Toxin production in marine microalgae was previously shown to be tightly coupled with cellular stoichiometry. The highest values of cellular toxin are in fact mainly associated with a high carbon to nutrient cellular ratio. In particular, the cellular accumulation of C-rich toxins (i.e., with C:N > 6.6) can be stimulated by both N and P deficiency. Dinoflagellates are the main producers of C-rich toxins and may represent a serious threat for human health and the marine ecosystem. As such, the development of a numerical model able to predict how toxin production is stimulated by nutrient supply/deficiency is of primary utility for both scientific and management purposes. In this work we have developed a mechanistic model describing the stoichiometric regulation of C-rich toxins in marine dinoflagellates. To this purpose, a new formulation describing toxin production and fate was embedded in the European Regional Seas Ecosystem Model (ERSEM), here simplified to describe a monospecific batch culture. Toxin production was assumed to be composed by two distinct additive terms; the first is a constant fraction of algal production and is assumed to take place at any physiological conditions. The second term is assumed to be dependent on algal biomass and to be stimulated by internal nutrient deficiency. By using these assumptions, the model reproduced the concentrations and temporal evolution of toxins observed in cultures of Ostreopsis cf. ovata, a benthic/epiphytic dinoflagellate producing C-rich toxins named ovatoxins. The analysis of simulations and their comparison with experimental data provided a conceptual model linking toxin production and nutritional status in this species. The model was also qualitatively validated by using independent literature data, and the results indicate that our formulation can be also used to simulate toxin dynamics in other dinoflagellates. Our model represents an important step towards the simulation and prediction of marine algal toxicity.