Establishment of new murine embryonic stem cell lines for the generation of mouse models of human genetic diseases

Embryonic stem cells are totipotent cells derived from the inner cell mass of blastocysts. Recently, the development of appropriate culture conditions for the differentiation of these cells into specific cell types has permitted their use as potential therapeutic agents for several diseases. In addition, manipulation of their genome in vitro allows the creation of animal models of human genetic diseases and for the study of gene function in vivo. We report the establishment of new lines of murine embryonic stem cells from preimplantation stage embryos of 129/Sv mice. Most of these cells had a normal karyotype and an XY sex chromosome composition. The pluripotent properties of the cell lines obtained were analyzed on the basis of their alkaline phosphatase activity and their capacity to form complex embryoid bodies with rhythmically contracting cardiomyocytes. Two lines, USP-1 and USP-3, with the best in vitro characteristics of pluripotency were used in chimera-generating experiments. The capacity to contribute to the germ line was demonstrated by the USP-1 cell line. This cell line is currently being used to generate mouse models of human diseases.

Serum metabolic profiling of human gastric cancer based on gas chromatography/mass spectrometry

Research on molecular mechanisms of carcinogenesis plays an important role in diagnosing and treating gastric cancer. Metabolic profiling may offer the opportunity to understand the molecular mechanism of carcinogenesis and help to non-invasively identify the potential biomarkers for the early diagnosis of human gastric cancer. The aims of this study were to explore the underlying metabolic mechanisms of gastric cancer and to identify biomarkers associated with morbidity. Gas chromatography/mass spectrometry (GC/MS) was used to analyze the serum metabolites of 30 Chinese gastric cancer patients and 30 healthy controls. Diagnostic models for gastric cancer were constructed using orthogonal partial least squares discriminant analysis (OPLS-DA). Acquired metabolomic data were analyzed by the nonparametric Wilcoxon test to find serum metabolic biomarkers for gastric cancer. The OPLS-DA model showed adequate discrimination between cancer and non-cancer cohorts while the model failed to discriminate different pathological stages (I-IV) of gastric cancer patients. A total of 44 endogenous metabolites such as amino acids, organic acids, carbohydrates, fatty acids, and steroids were detected, of which 18 differential metabolites were identified with significant differences. A total of 13 variables were obtained for their greatest contribution in the discriminating OPLS-DA model [variable importance in the projection (VIP) value >1.0], among which 11 metabolites were identified using both VIP values (VIP >1) and the Wilcoxon test. These metabolites potentially revealed perturbations of glycolysis and of amino acid, fatty acid, cholesterol, and nucleotide metabolism of gastric cancer patients. These results suggest that gastric cancer serum metabolic profiling has great potential in detecting this disease and helping to understand its metabolic mechanisms.

Percent body fat is a better predictor of cardiovascular risk factors than body mass index

The objective of the present study was to evaluate the predictive values of percent body fat (PBF) and body mass index (BMI) for cardiovascular risk factors, especially when PBF and BMI are conflicting. BMI was calculated by the standard formula and PBF was determined by bioelectrical impedance analysis. A total of 3859 ambulatory adult Han Chinese subjects (2173 males and 1686 females, age range: 18-85 years) without a history of cardiovascular diseases were recruited from February to September 2009. Based on BMI and PBF, they were classified into group 1 (normal BMI and PBF, N = 1961), group 2 (normal BMI, but abnormal PBF, N = 381), group 3 (abnormal BMI, but normal PBF, N = 681), and group 4 (abnormal BMI and PBF, N = 836). When age, gender, lifestyle, and family history of obesity were adjusted, PBF, but not BMI, was correlated with blood glucose and lipid levels. The odds ratio (OR) and 95% confidence interval (CI) for cardiovascular risk factors in groups 2 and 4 were 1.88 (1.45-2.45) and 2.06 (1.26-3.35) times those in group 1, respectively, but remained unchanged in group 3 (OR = 1.32, 95%CI = 0.92-1.89). Logistic regression models also demonstrated that PBF, rather than BMI, was independently associated with cardiovascular risk factors. In conclusion, PBF, and not BMI, is independently associated with cardiovascular risk factors, indicating that PBF is a better predictor.

Computational models for and from biology : simple gene assembly and reaction systems

The advancement of science and technology makes it clear that no single perspective is any longer sufficient to describe the true nature of any phenomenon. That is why the interdisciplinary research is gaining more attention overtime. An excellent example of this type of research is natural computing which stands on the borderline between biology and computer science. The contribution of research done in natural computing is twofold: on one hand, it sheds light into how nature works and how it processes information and, on the other hand, it provides some guidelines on how to design bio-inspired technologies. The first direction in this thesis focuses on a nature-inspired process called gene assembly in ciliates. The second one studies reaction systems, as a modeling framework with its rationale built upon the biochemical interactions happening within a cell. The process of gene assembly in ciliates has attracted a lot of attention as a research topic in the past 15 years. Two main modelling frameworks have been initially proposed in the end of 1990s to capture ciliates’ gene assembly process, namely the intermolecular model and the intramolecular model. They were followed by other model proposals such as templatebased assembly and DNA rearrangement pathways recombination models. In this thesis we are interested in a variation of the intramolecular model called simple gene assembly model, which focuses on the simplest possible folds in the assembly process. We propose a new framework called directed overlap-inclusion (DOI) graphs to overcome the limitations that previously introduced models faced in capturing all the combinatorial details of the simple gene assembly process. We investigate a number of combinatorial properties of these graphs, including a necessary property in terms of forbidden induced subgraphs. We also introduce DOI graph-based rewriting rules that capture all the operations of the simple gene assembly model and prove that they are equivalent to the string-based formalization of the model. Reaction systems (RS) is another nature-inspired modeling framework that is studied in this thesis. Reaction systems’ rationale is based upon two main regulation mechanisms, facilitation and inhibition, which control the interactions between biochemical reactions. Reaction systems is a complementary modeling framework to traditional quantitative frameworks, focusing on explicit cause-effect relationships between reactions. The explicit formulation of facilitation and inhibition mechanisms behind reactions, as well as the focus on interactions between reactions (rather than dynamics of concentrations) makes their applicability potentially wide and useful beyond biological case studies. In this thesis, we construct a reaction system model corresponding to the heat shock response mechanism based on a novel concept of dominance graph that captures the competition on resources in the ODE model. We also introduce for RS various concepts inspired by biology, e.g., mass conservation, steady state, periodicity, etc., to do model checking of the reaction systems based models. We prove that the complexity of the decision problems related to these properties varies from P to NP- and coNP-complete to PSPACE-complete. We further focus on the mass conservation relation in an RS and introduce the conservation dependency graph to capture the relation between the species and also propose an algorithm to list the conserved sets of a given reaction system.

Models to study intermediate filament dynamics and functions

Intermediate filaments are part of the cytoskeleton and nucleoskeleton; they provide cells with structure and have important roles in cell signalling. The IFs are a large protein family with more than 70 members; each tightly regulated and expressed in a cell type-specific manner. Although the IFs have been known and studied for decades, our knowledge about their specific functions is still limited, despite the fact that mutations in IF genes cause numerous severe human diseases. In this work, three IF proteins are examined more closely; the nuclear lamin A/C and the cytoplasmic nestin and vimentin. In particular the regulation of lamin A/C dynamics, the role of nestin in muscle and body homeostasis as well as the functions and evolutionary aspects of vimentin are investigated. Together this data highlights some less well understood functions of these IFs. We used mass-spectrometry to identify inter-phase specific phosphorylation sites on lamin A. With the use of genetically engineered lamin A protein in combination with high resolution microscopy and biochemical methods we discovered novel roles for this phosphorylation in regulation of lamin dynamics. More specifically, our data suggests that the phosphorylation of certain amino acids in lamin A determines the localization and dynamics of the protein. In addition, we present results demonstrating that lamin A regulates Cdk5-activity. In the second study we use mice lacking nestin to gain more knowledge of this seldom studied protein. Our results show that nestin is essential for muscle regeneration; mice lacking nestin recover more slowly from muscle injury and show signs of spontaneous muscle regeneration, indicating that their muscles are more sensitive to stresses and injury. The absence of nestin also leads to decreased over-all muscle mass and slower body growth. Furthermore, nestin has a role in controlling testicle homeostasis as nestin-/- male mice show a greater variation in testicle size. The common fruit fly Drosophila melanogaster lacks cytoplasmic IFs as most insects do. By creating a fly that expresses human vimentin we establish a new research platform for vimentin studies, as well as provide a new tool for the studies of IF evolution.

Interactions of diploid embryonal carcinoma cells and normal cell mass cells with 2.5 day mouse embryos following aggregation

Spontaneous teratocarcinomas are ovarian or testicular tumors which have their origins in germ cells. The tumors contain a disorganized array of benign differentiated cells as well as an undifferentiated population of malignant stem cells, the embryonal carcinoma or EC cells. These pluripotent stem cells in tissue culture share many properties with the transient pluripotent cells of the early embryo, and might therefore serve as models for the investigation of developmental events ill vitro. The property of EC cells of prime interest in this study is an in vivo phenomenon. Certain EC cell lines are known to be regulated ill vivo and to differentiate normally in association with normal embryonic cells, resulting in chimeric mice. These mice have two genetically distinct cell populations, one of which is derived from the originally malignant EC cells. This has usually been accomplished by injection of the EC cells into the Day 3 blastocyst. In this study, the interactions between earlier stage embryos and EC cells have been tested by aggregating clumps of EC cells with Day 2 embryos. The few previous aggregation studies produced a high degree of abnormality in chimeric embryos, but the EC cells employed had known chromosomal abnormalities. In this study, two diploid EC cell lines (P19 and Pi0) were aggregated with 2.5 day mouse embryos, and were found to behave quite differently in the embryonic environment. P19 containing aggregates generally resorbed early, and the few embryos recovered at midgestation were normal and non-chimeric. Pi0 containing aggregates survived in high numbers to midgestation, and the Pi0 cells were very successful in colonizing the embryo. All these embryos were chimeric, and the contribution by the EC cells to each chimera was very high. However, these heavily chimeric embryos were all abnormal. Blastocyst injection had previously produced some abnormal embryos with high Pl0 contributions in addition to the live born mice, which had lower EC contributions. This study now adds more support to the hypothesis that high EC contributions may be incompatible with normal development. The possibility that the abnormalities were due to the mixing of temporally asynchronous embryonic cell types in the aggregates was tested by aggregating normal pluripotent cells taken from 3.5 day embryos with 2.5 day embryos. Early embryo loss was very high, and histological studies showed that the majority of these embryos died by 6.5 days development. Some embryos escaped this early death such that some healthy chimeras were recovered, in contrast to recovery of abnormal chimeric embryos following Pl0-morula aggregations, and non-chimeric embryos following P19-morula aggregations. This somewhat surprising adverse effect on development following aggregation of normal cell types suggests that there are developmental difficulties associated with the mixing of asynchronous cell types in aggregates. However, the greater magnitude of the adverse effects when the aggregates contained tumor derived cells suggests that EC cells should not be considered the complete equivalent of the pluripotent cells of the early embryo.

The relationship between body mass index and breast cancer recurrence/progression and breast cancer-specific death

The aim of this study was to describe the nonlinear association between body mass index (BMI) and breast cancer outcomes and to determine whether BMI improves prediction of outcomes. A cohort of906 breast cancer patients diagnosed at Henry Ford Health System, Detroit (1985-1990) were studied. The median follow-up was 10 years. Multivariate logistic regression was used to model breast cancer recurrence/progression and breast cancer-specific death. Restricted cubic splines were used to model nonlinear effects. Receiver operator characteristic areas under the curves (ROC AUC) were used to evaluate prediction. BMI was nonlinearly associated with recurrence/progression and death (p= 0.0230 and 0.0101). Probability of outcomes increased with increase or decrease ofBMI away from 25. BMI splines were suggestive of improved prediction of death. The ROC AUCs for nested models with and without BMI were 0.8424 and 0.8331 (p= 0.08). I f causally associated, modifying patients BMI towards 25 may improve outcomes.

Some concepts and models usefulninthe analysis of discrete life time data

This thesis entitled Reliability Modelling and Analysis in Discrete time Some Concepts and Models Useful in the Analysis of discrete life time data.The present study consists of five chapters. In Chapter II we take up the derivation of some general results useful in reliability modelling that involves two component mixtures. Expression for the failure rate, mean residual life and second moment of residual life of the mixture distributions in terms of the corresponding quantities in the component distributions are investigated. Some applications of these results are also pointed out. The role of the geometric,Waring and negative hypergeometric distributions as models of life lengths in the discrete time domain has been discussed already. While describing various reliability characteristics, it was found that they can be often considered as a class. The applicability of these models in single populations naturally extends to the case of populations composed of sub-populations making mixtures of these distributions worth investigating. Accordingly the general properties, various reliability characteristics and characterizations of these models are discussed in chapter III. Inference of parameters in mixture distribution is usually a difficult problem because the mass function of the mixture is a linear function of the component masses that makes manipulation of the likelihood equations, leastsquare function etc and the resulting computations.very difficult. We show that one of our characterizations help in inferring the parameters of the geometric mixture without involving computational hazards. As mentioned in the review of results in the previous sections, partial moments were not studied extensively in literature especially in the case of discrete distributions. Chapters IV and V deal with descending and ascending partial factorial moments. Apart from studying their properties, we prove characterizations of distributions by functional forms of partial moments and establish recurrence relations between successive moments for some well known families. It is further demonstrated that partial moments are equally efficient and convenient compared to many of the conventional tools to resolve practical problems in reliability modelling and analysis. The study concludes by indicating some new problems that surfaced during the course of the present investigation which could be the subject for a future work in this area.

Impact of sward type, cutting date and conditioning temperature on mass and energy flows in the integrated generation of solid fuel and biogas from semi-natural grassland

Energy production from biomass and the conservation of ecologically valuable grassland habitats are two important issues of agriculture today. The combination of a bioenergy production, which minimises environmental impacts and competition with food production for land with a conversion of semi-natural grasslands through new utilization alternatives for the biomass, led to the development of the IFBB process. Its basic principle is the separation of biomass into a liquid fraction (press fluid, PF) for the production of electric and thermal energy after anaerobic digestion to biogas and a solid fraction (press cake, PC) for the production of thermal energy through combustion. This study was undertaken to explore mass and energy flows as well as quality aspects of energy carriers within the IFBB process and determine their dependency on biomass-related and technical parameters. Two experiments were conducted, in which biomass from semi-natural grassland was conserved as silage and subjected to a hydrothermal conditioning and a subsequent mechanical dehydration with a screw press. Methane yield of the PF and the untreated silage was determined in anaerobic digestion experiments in batch fermenters at 37°C with a fermentation time of 13-15 and 27-35 days for the PF and the silage, respectively. Concentrations of dry matter (DM), ash, crude protein (CP), crude fibre (CF), ether extract (EE), neutral detergent fibre (NDF), acid detergent fibre (ADF), acid detergent ligning (ADL) and elements (K, Mg, Ca, Cl, N, S, P, C, H, N) were determined in the untreated biomass and the PC. Higher heating value (HHV) and ash softening temperature (AST) were calculated based on elemental concentration. Chemical composition of the PF and mass flows of all plant compounds into the PF were calculated. In the first experiment, biomass from five different semi-natural grassland swards (Arrhenaterion I and II, Caricion fuscae, Filipendulion ulmariae, Polygono-Trisetion) was harvested at one late sampling (19 July or 31 August) and ensiled. Each silage was subjected to three different temperature treatments (5°C, 60°C, 80°C) during hydrothermal conditioning. Based on observed methane yields and HHV as energy output parameters as well as literature-based and observed energy input parameters, energy and green house gas (GHG) balances were calculated for IFBB and two reference conversion processes, whole-crop digestion of untreated silage (WCD) and combustion of hay (CH). In the second experiment, biomass from one single semi-natural grassland sward (Arrhenaterion) was harvested at eight consecutive dates (27/04, 02/05, 09/05, 16/05, 24/05, 31/05, 11/06, 21/06) and ensiled. Each silage was subjected to six different treatments (no hydrothermal conditioning and hydrothermal conditioning at 10°C, 30°C, 50°C, 70°C, 90°C). Energy balance was calculated for IFBB and WCD. Multiple regression models were developed to predict mass flows, concentrations of elements in the PC, concentration of organic compounds in the PF and energy conversion efficiency of the IFBB process from temperature of hydrothermal conditioning as well as NDF and DM concentration in the silage. Results showed a relative reduction of ash and all elements detrimental for combustion in the PC compared to the untreated biomass of 20-90%. Reduction was highest for K and Cl and lowest for N. HHV of PC and untreated biomass were in a comparable range (17.8-19.5 MJ kg-1 DM), but AST of PC was higher (1156-1254°C). Methane yields of PF were higher compared to those of WCD when the biomass was harvested late (end of May and later) and in a comparable range when the biomass was harvested early and ranged from 332 to 458 LN kg-1 VS. Regarding energy and GHG balances, IFBB, with a net energy yield of 11.9-14.1 MWh ha-1, a conversion efficiency of 0.43-0.51, and GHG mitigation of 3.6-4.4 t CO2eq ha-1, performed better than WCD, but worse than CH. WCD produces thermal and electric energy with low efficiency, CH produces only thermal energy with a low quality solid fuel with high efficiency, IFBB produces thermal and electric energy with a solid fuel of high quality with medium efficiency. Regression models were able to predict target parameters with high accuracy (R2=0.70-0.99). The influence of increasing temperature of hydrothermal conditioning was an increase of mass flows, a decrease of element concentrations in the PC and a differing effect on energy conversion efficiency. The influence of increasing NDF concentration of the silage was a differing effect on mass flows, a decrease of element concentrations in the PC and an increase of energy conversion efficiency. The influence of increasing DM concentration of the silage was a decrease of mass flows, an increase of element concentrations in the PC and an increase of energy conversion efficiency. Based on the models an optimised IFBB process would be obtained with a medium temperature of hydrothermal conditioning (50°C), high NDF concentrations in the silage and medium DM concentrations of the silage.

Compositional evolution with mass transfer in closed systems

Evolution of compositions in time, space, temperature or other covariates is frequent in practice. For instance, the radioactive decomposition of a sample changes its composition with time. Some of the involved isotopes decompose into other isotopes of the sample, thus producing a transfer of mass from some components to other ones, but preserving the total mass present in the system. This evolution is traditionally modelled as a system of ordinary di erential equations of the mass of each component. However, this kind of evolution can be decomposed into a compositional change, expressed in terms of simplicial derivatives, and a mass evolution (constant in this example). A rst result is that the simplicial system of di erential equations is non-linear, despite of some subcompositions behaving linearly. The goal is to study the characteristics of such simplicial systems of di erential equa- tions such as linearity and stability. This is performed extracting the compositional dif ferential equations from the mass equations. Then, simplicial derivatives are expressed in coordinates of the simplex, thus reducing the problem to the standard theory of systems of di erential equations, including stability. The characterisation of stability of these non-linear systems relays on the linearisation of the system of di erential equations at the stationary point, if any. The eigenvelues of the linearised matrix and the associated behaviour of the orbits are the main tools. For a three component system, these orbits can be plotted both in coordinates of the simplex or in a ternary diagram. A characterisation of processes with transfer of mass in closed systems in terms of stability is thus concluded. Two examples are presented for illustration, one of them is a radioactive decay

L'evolució dels models de turisme litoral: el Regne Unit, la Costa brava i Cuba

La tesi es centra en una anàlisi prèvia de les implicacions de la teoria del cicle de vida aplicada a les destinacions turístiques pel geògraf canadenc Richard W. Butler (1982) i en un estudi del significat de l'aparició del postmodernisme i les seves repercussions en el turisme, a la vegada que es resegueix la història que ha fet possible l'aparició del turisme massificat de sol i platja. Tot això és la base teòrica imprescindible per poder estudiar empíricament les destinacions de Saltburn (Regne Unit), l'Estartit (Costa Brava) i Cayo Coco (Cuba), utilitzant el cicle de vida per escriure la seva història i veure les polítiques de desenvolupament i regeneració turística que s'han seguit i si aquestes es poden qualificar de postfordistes. La principal conclusió és que el concepte del cicle de vida es limita a ser una eina de diagnòstic a posteriori. El cicle de vida, com a instrument prescriptiu, no serveix. Es específic per a cada destinació, amb unes etapes i punts d'inflexió que només es poden establir en restrospectiva. El cicle de vida és una eina descriptiva molt útil per entendre l'evolució de les destinacions turístiques i els seus mercats, però vigilant de no caure en falses exploracions o perillosos determinismes.

A review of the theoretical basis for bulk mass flux convective parameterization

Most parameterizations for precipitating convection in use today are bulk schemes, in which an ensemble of cumulus elements with different properties is modelled as a single, representative entraining-detraining plume. We review the underpinning mathematical model for such parameterizations, in particular by comparing it with spectral models in which elements are not combined into the representative plume. The chief merit of a bulk model is that the representative plume can be described by an equation set with the same structure as that which describes each element in a spectral model. The equivalence relies on an ansatz for detrained condensate introduced by Yanai et al. (1973) and on a simplified microphysics. There are also conceptual differences in the closure of bulk and spectral parameterizations. In particular, we show that the convective quasi-equilibrium closure of Arakawa and Schubert (1974) for spectral parameterizations cannot be carried over to a bulk parameterization in a straightforward way. Quasi-equilibrium of the cloud work function assumes a timescale separation between a slow forcing process and a rapid convective response. But, for the natural bulk analogue to the cloud-work function (the dilute CAPE), the relevant forcing is characterised by a different timescale, and so its quasi-equilibrium entails a different physical constraint. Closures of bulk parameterization that use the non-entraining parcel value of CAPE do not suffer from this timescale issue. However, the Yanai et al. (1973) ansatz must be invoked as a necessary ingredient of those closures.

Over-parameterised, uncertain 'mathematical marionettes' - How can we best use catchment water quality models? An example of an 80-year catchment-scale nutrient balance

The Integrated Catchment Model of Nitrogen (INCA-N) was applied to the River Lambourn, a Chalk river-system in southern England. The model's abilities to simulate the long-term trend and seasonal patterns in observed stream water nitrate concentrations from 1920 to 2003 were tested. This is the first time a semi-distributed, daily time-step model has been applied to simulate such a long time period and then used to calculate detailed catchment nutrient budgets which span the conversion of pasture to arable during the late 1930s and 1940s. Thus, this work goes beyond source apportionment and looks to demonstrate how such simulations can be used to assess the state of the catchment and develop an understanding of system behaviour. The mass-balance results from 1921, 1922, 1991, 2001 and 2002 are presented and those for 1991 are compared to other modelled and literature values of loads associated with nitrogen soil processes and export. The variations highlighted the problem of comparing modelled fluxes with point measurements but proved useful for identifying the most poorly understood inputs and processes thereby providing an assessment of input data and model structural uncertainty. The modelled terrestrial and instream mass-balances also highlight the importance of the hydrological conditions in pollutant transport. Between 1922 and 2002, increased inputs of nitrogen from fertiliser, livestock and deposition have altered the nitrogen balance with a shift from possible reduction in soil fertility but little environmental impact in 1922, to a situation of nitrogen accumulation in the soil, groundwater and instream biota in 2002. In 1922 and 2002 it was estimated that approximately 2 and 18 kg N ha(-1) yr(-1) respectively were exported from the land to the stream. The utility of the approach and further considerations for the best use of models are discussed. (C) 2008 Elsevier B.V. All rights reserved.

Aqueous geochemistry and oxygen isotope compositions of acid mine drainage from the Río Tinto, SW Spain, highlight inconsistencies in current models

The Rio Tinto river in SW Spain is a classic example of acid mine drainage and the focus of an increasing amount of research including environmental geochemistry, extremophile microbiology and Mars-analogue studies. Its 5000-year mining legacy has resulted in a wide range of point inputs including spoil heaps and tunnels draining underground workings. The variety of inputs and importance of the river as a research site make it an ideal location for investigating sulphide oxidation mechanisms at the field scale. Mass balance calculations showed that pyrite oxidation accounts for over 93% of the dissolved sulphate derived from sulphide oxidation in the Rio Tinto point inputs. Oxygen isotopes in water and sulphate were analysed from a variety of drainage sources and displayed delta O-18((SO4-H2O)) values from 3.9 to 13.6 parts per thousand, indicating that different oxidation pathways occurred at different sites within the catchment. The most commonly used approach to interpreting field oxygen isotope data applies water and oxygen fractionation factors derived from laboratory experiments. We demonstrate that this approach cannot explain high delta O-18((SO4-H2O)) values in a manner that is consistent with recent models of pyrite and sulphoxyanion oxidation. In the Rio Tinto, high delta O-18((SO4-H2O)) values (11.2-13.6 parts per thousand) occur in concentrated (Fe = 172-829 mM), low pH (0.88-1.4), ferrous iron (68-91% of total Fe) waters and are most simply explained by a mechanism involving a dissolved sulphite intermediate, sulphite-water oxygen equilibrium exchange and finally sulphite oxidation to sulphate with O-2. In contrast, drainage from large waste blocks of acid volcanic tuff with pyritiferous veins also had low pH (1.7). but had a low delta O-18((SO4-H2O)) value of 4.0 parts per thousand and high concentrations of ferric iron (Fe(III) = 185 mM, total Fe = 186 mM), suggesting a pathway where ferric iron is the primary oxidant, water is the primary source of oxygen in the sulphate and where sulphate is released directly from the pyrite surface. However, problems remain with the sulphite-water oxygen exchange model and recommendations are therefore made for future experiments to refine our understanding of oxygen isotopes in pyrite oxidation. (C) 2009 Elsevier B.V. All rights reserved.

Over-parameterised, uncertain 'mathematical marionettes' - How can we best use catchment water quality models? An example of an 80-year catchment-scale nutrient balance

The Integrated Catchment Model of Nitrogen (INCA-N) was applied to the River Lambourn, a Chalk river-system in southern England. The model's abilities to simulate the long-term trend and seasonal patterns in observed stream water nitrate concentrations from 1920 to 2003 were tested. This is the first time a semi-distributed, daily time-step model has been applied to simulate such a long time period and then used to calculate detailed catchment nutrient budgets which span the conversion of pasture to arable during the late 1930s and 1940s. Thus, this work goes beyond source apportionment and looks to demonstrate how such simulations can be used to assess the state of the catchment and develop an understanding of system behaviour. The mass-balance results from 1921, 1922, 1991, 2001 and 2002 are presented and those for 1991 are compared to other modelled and literature values of loads associated with nitrogen soil processes and export. The variations highlighted the problem of comparing modelled fluxes with point measurements but proved useful for identifying the most poorly understood inputs and processes thereby providing an assessment of input data and model structural uncertainty. The modelled terrestrial and instream mass-balances also highlight the importance of the hydrological conditions in pollutant transport. Between 1922 and 2002, increased inputs of nitrogen from fertiliser, livestock and deposition have altered the nitrogen balance with a shift from possible reduction in soil fertility but little environmental impact in 1922, to a situation of nitrogen accumulation in the soil, groundwater and instream biota in 2002. In 1922 and 2002 it was estimated that approximately 2 and 18 kg N ha(-1) yr(-1) respectively were exported from the land to the stream. The utility of the approach and further considerations for the best use of models are discussed. (C) 2008 Elsevier B.V. All rights reserved.