Adaptive MCMC methods with applications in environmental and geophysical models

This work presents new, efficient Markov chain Monte Carlo (MCMC) simulation methods for statistical analysis in various modelling applications. When using MCMC methods, the model is simulated repeatedly to explore the probability distribution describing the uncertainties in model parameters and predictions. In adaptive MCMC methods based on the Metropolis-Hastings algorithm, the proposal distribution needed by the algorithm learns from the target distribution as the simulation proceeds. Adaptive MCMC methods have been subject of intensive research lately, as they open a way for essentially easier use of the methodology. The lack of user-friendly computer programs has been a main obstacle for wider acceptance of the methods. This work provides two new adaptive MCMC methods: DRAM and AARJ. The DRAM method has been built especially to work in high dimensional and non-linear problems. The AARJ method is an extension to DRAM for model selection problems, where the mathematical formulation of the model is uncertain and we want simultaneously to fit several different models to the same observations. The methods were developed while keeping in mind the needs of modelling applications typical in environmental sciences. The development work has been pursued while working with several application projects. The applications presented in this work are: a winter time oxygen concentration model for Lake Tuusulanjärvi and adaptive control of the aerator; a nutrition model for Lake Pyhäjärvi and lake management planning; validation of the algorithms of the GOMOS ozone remote sensing instrument on board the Envisat satellite of European Space Agency and the study of the effects of aerosol model selection on the GOMOS algorithm.

A spatial modelling framework for assessing climate change impacts on freshwater ecosystems: Response of brown trout (Salmo trutta L.) biomass to warming water temperature

Mountain regions worldwide are particularly sensitive to on-going climate change. Specifically in the Alps in Switzerland, the temperature has increased twice as fast than in the rest of the Northern hemisphere. Water temperature closely follows the annual air temperature cycle, severely impacting streams and freshwater ecosystems. In the last 20 years, brown trout (Salmo trutta L) catch has declined by approximately 40-50% in many rivers in Switzerland. Increasing water temperature has been suggested as one of the most likely cause of this decline. Temperature has a direct effect on trout population dynamics through developmental and disease control but can also indirectly impact dynamics via food-web interactions such as resource availability. We developed a spatially explicit modelling framework that allows spatial and temporal projections of trout biomass using the Aare river catchment as a model system, in order to assess the spatial and seasonal patterns of trout biomass variation. Given that biomass has a seasonal variation depending on trout life history stage, we developed seasonal biomass variation models for three periods of the year (Autumn-Winter, Spring and Summer). Because stream water temperature is a critical parameter for brown trout development, we first calibrated a model to predict water temperature as a function of air temperature to be able to further apply climate change scenarios. We then built a model of trout biomass variation by linking water temperature to trout biomass measurements collected by electro-fishing in 21 stations from 2009 to 2011. The different modelling components of our framework had overall a good predictive ability and we could show a seasonal effect of water temperature affecting trout biomass variation. Our statistical framework uses a minimum set of input variables that make it easily transferable to other study areas or fish species but could be improved by including effects of the biotic environment and the evolution of demographical parameters over time. However, our framework still remains informative to spatially highlight where potential changes of water temperature could affect trout biomass. (C) 2015 Elsevier B.V. All rights reserved.-

Accessibility, equity and the sharing of water resources

While the supply of water to dry or arid mountain regions has long been a major challenge, the on-going processes of climatic and socio-economic change currently affecting the hydrosystems of the Alps raise the spectre of renewed pressure on water resources and possible local shortages. In such a context, questions relating to fair distribution of water are all the more sensitive given the tendency to neglect the social dimension of sustainability. The present paper makes both a conceptual and empirical contribution to this debate by analysing a system of distribution that has a long experience of water scarcity management: the community governance models traditionally linked to the irrigation channels, or bisses, typical of the Swiss Alpine canton of Valais. More specifically, we evaluate these models in terms of accessibility and equity, characteristics that we use to operationalize the notion of 'fair distribution'. We examine these dimensions in three case studies with a view to highlighting the limitations of the aforementioned models. Indeed, despite their cooperative and endogenous nature, they tend to not only exclude certain members of the population, but also to reproduce rather than reduce social inequalities within the community. In general, these results challenge the rosy picture generally found in the literature relating to these community governance models.

Non-wood massojen vedenpoiston arvioiminen märkäpuristuksessa

Yksivuotisten kasvien (non-wood) kuitua verrataan usein lehtipuukuituihin. Käytetyimpiä non-wood kasveja ovat vehnän olki, bambu, järviruoko ja bagassi. Non-wood massan erottaa puumassasta kuitenkin korkea silikaattipitoisuus sekä parenkyymisolupitoisuus, joka antaa massalle korkean hienoainepitoisuuden. Tämä yhdessä korkean hemiselluloosapitoisuuden ja lyhyen kuidun pituuden kanssa heikentävät voimakkaasti non-wood massan vedenpoisto-ominaisuuksia. Non-wood kuidulla voidaan korvata lehtipuumassaa hienopapereissa. Non-wood kuitu antaa paperille hyvän opasiteetin, korkean valonsirontakertoimen sekä sileän painopinnan. Massaan lisättävä pitkäkuituinen havupuumassa parantaa ajetta¬vuutta paperikoneella ja helpottaa massan vedenpoistoa. Non-wood massan vedenpoistoa voidaan tehostaa esimerkiksi poistamalla osa hienoaineesta, käyttämällä non-wood massalle sopivaa keittotapaa sekä käyttämällä märkä¬puristuksessa pitkänippityyppistä puristinratkaisua. Myös non-wood kuidun kuivaaminen parantaa vedenpoistoa. Tässä tutkimuksessa kirjallisuusosassa keskityttiin yleisimpiin paperin valmistuksessa käytettäviin non-wood kuidun lähteisiin, märkäpuristuksen teoriaan ja tapoihin tehostaa vedenpoistoa. Kokeellisessa osassa tutkittiin vehnänolkimassan käyttäytymistä märkäpuristuksessa erilaisten ominaisuuksien pohjalta. Tutkimuksen kohteena oli non-wood massan keittotapa (hapan/alkali), hienoainepitoisuus, silikaattipitoisuus sekä kuivattu/kuivaamaton kuitu. Vertailun vuoksi tutkimuksessa oli mukana myös yksi järviruokomassa. Tuloksista huomataan, että non-wood massan vedenpoistoon vaikuttaa hienoainepitoisuus, kuidun kuivaus sekä massan valmistustapa. Järviruokomassan veden¬poisto on tehokkaampaa kuin vehnänolkimassan paremman kuitu¬koostumuksensa takia. Jos hienopaperimassassa korvataan lehtipuumassaa non-wood kuidulla maksimissaan 40 %, massan vedenpoistoa voidaan hyvin arvioida erilaisten suotautuvuusmittojen, kuten freeneksen, vedenpidätyskyvyn ja suotautumisajan, avulla.

Nutrient enrichment and water depth modify consumer control in rocky shore macroalgal communities

The objective of my thesis is to assess mechanisms of ecological community control in macroalgal communities in the Baltic Sea. In the top-down model, predatory fish feed on invertebrate mesograzers, releasing algae partly from grazing pressure. Such a reciprocal relationship is called trophic cascade. In the bottom-up model, nutrients increase biomass in the food chain. The nutrients are first assimilated by algae and, via food chain, increase also abundance of grazers and predators. Previous studies on oceanic shores have described these two regulative mechanisms in the grazer - alga link, but how they interact in the trophic cascades from fish to algae is still inadequately known. Because the top-down and bottom-up mechanisms are predicted to depend on environmental disturbances, such as wave stress and light, I have studied these models at two distinct water depths. There are five factorial field experiments behind the thesis, which were all conducted in the Finnish Archipelago Sea. In all the experiments, I studied macroalgal colonization - either density, filament length or biomass - on submerged colonization substrates. By excluding predatory fish and mesograzers from the algal communities, the studies compared the strength of the top-down control to natural algal communities. A part of the experimental units were, in addition, exposed to enriched nitrogen and phosphorus concentrations, which enabled testing of bottom-up control. These two models of community control were further investigated in shallow (<1 m) and deep (ca. 3 m) water. Moreover, the control mechanisms were also expected to depend on grazer species. Therefore different grazer species were enclosed into experimental units and their impacts on macroalgal communities were followed specifically. The community control in the Baltic rocky shores was found to follow theoretical predictions, which have not been confirmed by field studies before. Predatory fish limited grazing impact, which was seen as denser algal communities and longer algal filaments. Nutrient enrichment increased density and filament length of annual algae and, thus, changed the species composition of the algal community. The perennial alga Fucus vesiculosusA and the red alga Ceramium tenuicorne suffered from the increased nutrient availabilities. The enriched nutrient conditions led to denser grazer fauna, thereby causing strong top-down control over both the annual and perennial macroalgae. The strength of the top-down control seemed to depend on the density and diversity of grazers and predators as well as on the species composition of macroalgal assemblages. The nutrient enrichment led to, however, weaker limiting impact of predatory fish on grazer fauna, because fish stocks did not respond as quickly to enhanced resources in the environment as the invertebrate fauna. According to environmental stress model, environmental disturbances weaken the top-down control. For example, on a wave-exposed shore, wave stress causes more stress to animals close to the surface than deeper on the shore. Mesograzers were efficient consumers at both the depths, while predation by fish was weaker in shallow water. Thus, the results supported the environmental stress model, which predicts that environmental disturbance affects stronger the higher a species is in the food chain. This thesis assessed the mechanisms of community control in three-level food chains and did not take into account higher predators. Such predators in the Baltic Sea are, for example, cormorant, seals, white-tailed sea eagle, cod and salmon. All these predatory species were recently or are currently under intensive fishing, hunting and persecution, and their stocks have only recently increased in the region. Therefore, it is possible that future densities of top predators may yet alter the strengths of the controlling mechanisms in the Baltic littoral zone.

Adaptation of a Gaussia princeps Luciferase reporter system in Candida albicans for in vivo detection in the Galleria mellonella infection model.

For the past 10 years, mini-host models and in particular the greater wax moth Galleria mellonella have tended to become a surrogate for murine models of fungal infection mainly due to cost, ethical constraints and ease of use. Thus, methods to better assess the fungal pathogenesis in G. mellonella need to be developed. In this study, we implemented the detection of Candida albicans cells expressing the Gaussia princeps luciferase in its cell wall in infected larvae of G. mellonella. We demonstrated that detection and quantification of luminescence in the pulp of infected larvae is a reliable method to perform drug efficacy and C. albicans virulence assays as compared to fungal burden assay. Since the linearity of the bioluminescent signal, as compared to the CFU counts, has a correlation of R(2) = 0.62 and that this method is twice faster and less labor intensive than classical fungal burden assays, it could be applied to large scale studies. We next visualized and followed C. albicans infection in living G. mellonella larvae using a non-toxic and water-soluble coelenterazine formulation and a CCD camera that is commonly used for chemoluminescence signal detection. This work allowed us to follow for the first time C. albicans course of infection in G. mellonella during 4 days.

In Vivo Longitudinal (1)H MRS Study of Transgenic Mouse Models of Prion Disease in the Hippocampus and Cerebellum at 14.1 T.

In vivo (1)H MR spectroscopy allows the non invasive characterization of brain metabolites and it has been used for studying brain metabolic changes in a wide range of neurodegenerative diseases. The prion diseases form a group of fatal neurodegenerative diseases, also described as transmissible spongiform encephalopathies. The mechanism by which prions elicit brain damage remains unclear and therefore different transgenic mouse models of prion disease were created. We performed an in vivo longitudinal (1)H MR spectroscopy study at 14.1 T with the aim to measure the neurochemical profile of Prnp -/- and PrPΔ32-121 mice in the hippocampus and cerebellum. Using high-field MR spectroscopy we were able to analyze in details the in vivo brain metabolites in Prnp -/- and PrPΔ32-121 mice. An increase of myo-inositol, glutamate and lactate concentrations with a decrease of N-acetylaspartate concentrations were observed providing additional information to the previous measurements.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.

Non-Linear Journal Bearing Model for Analysis of Superharmonic Vibrations of Rotor Systems

A rotating machine usually consists of a rotor and bearings that supports it. The nonidealities in these components may excite vibration of the rotating system. The uncontrolled vibrations may lead to excessive wearing of the components of the rotating machine or reduce the process quality. Vibrations may be harmful even when amplitudes are seemingly low, as is usually the case in superharmonic vibration that takes place below the first critical speed of the rotating machine. Superharmonic vibration is excited when the rotational velocity of the machine is a fraction of the natural frequency of the system. In such a situation, a part of the machine’s rotational energy is transformed into vibration energy. The amount of vibration energy should be minimised in the design of rotating machines. The superharmonic vibration phenomena can be studied by analysing the coupled rotor-bearing system employing a multibody simulation approach. This research is focused on the modelling of hydrodynamic journal bearings and rotorbearing systems supported by journal bearings. In particular, the non-idealities affecting the rotor-bearing system and their effect on the superharmonic vibration of the rotating system are analysed. A comparison of computationally efficient journal bearing models is carried out in order to validate one model for further development. The selected bearing model is improved in order to take the waviness of the shaft journal into account. The improved model is implemented and analyzed in a multibody simulation code. A rotor-bearing system that consists of a flexible tube roll, two journal bearings and a supporting structure is analysed employing the multibody simulation technique. The modelled non-idealities are the shell thickness variation in the tube roll and the waviness of the shaft journal in the bearing assembly. Both modelled non-idealities may cause subharmonic resonance in the system. In multibody simulation, the coupled effect of the non-idealities can be captured in the analysis. Additionally one non-ideality is presented that does not excite the vibrations itself but affects the response of the rotorbearing system, namely the waviness of the bearing bushing which is the non-rotating part of the bearing system. The modelled system is verified with measurements performed on a test rig. In the measurements the waviness of bearing bushing was not measured and therefore it’s affect on the response was not verified. In conclusion, the selected modelling approach is an appropriate method when analysing the response of the rotor-bearing system. When comparing the simulated results to the measured ones, the overall agreement between the results is concluded to be good.

Modeling reaction kinetics and mass transfer in ozonation in water solutions

This dissertation is based on four articles dealing with modeling of ozonation. The literature part of this considers some models for hydrodynamics in bubble column simulation. A literature review of methods for obtaining mass transfer coefficients is presented. The methods presented to obtain mass transfer are general models and can be applied to any gas-liquid system. Ozonation reaction models and methods for obtaining stoichiometric coefficients and reaction rate coefficients for ozonation reactions are discussed in the final section of the literature part. In the first article, ozone gas-liquid mass transfer into water in a bubble column was investigated for different pH values. A more general method for estimation of mass transfer and Henry’s coefficient was developed from the Beltrán method. The ozone volumetric mass transfer coefficient and the Henry’s coefficient were determined simultaneously by parameter estimation using a nonlinear optimization method. A minor dependence of the Henry’s law constant on pH was detected at the pH range 4 - 9. In the second article, a new method using the axial dispersion model for estimation of ozone self-decomposition kinetics in a semi-batch bubble column reactor was developed. The reaction rate coefficients for literature equations of ozone decomposition and the gas phase dispersion coefficient were estimated and compared with the literature data. The reaction order in the pH range 7-10 with respect to ozone 1.12 and 0.51 the hydroxyl ion were obtained, which is in good agreement with literature. The model parameters were determined by parameter estimation using a nonlinear optimization method. Sensitivity analysis was conducted using object function method to obtain information about the reliability and identifiability of the estimated parameters. In the third article, the reaction rate coefficients and the stoichiometric coefficients in the reaction of ozone with the model component p-nitrophenol were estimated at low pH of water using nonlinear optimization. A novel method for estimation of multireaction model parameters in ozonation was developed. In this method the concentration of unknown intermediate compounds is presented as a residual COD (chemical oxygen demand) calculated from the measured COD and the theoretical COD for the known species. The decomposition rate of p-nitrophenol on the pathway producing hydroquinone was found to be about two times faster than the p-nitrophenol decomposition rate on the pathway producing 4- nitrocatechol. In the fourth article, the reaction kinetics of p-nitrophenol ozonation was studied in a bubble column at pH 2. Using the new reaction kinetic model presented in the previous article, the reaction kinetic parameters, rate coefficients, and stoichiometric coefficients as well as the mass transfer coefficient were estimated with nonlinear estimation. The decomposition rate of pnitrophenol was found to be equal both on the pathway producing hydroquinone and on the path way producing 4-nitrocathecol. Comparison of the rate coefficients with the case at initial pH 5 indicates that the p-nitrophenol degradation producing 4- nitrocathecol is more selective towards molecular ozone than the reaction producing hydroquinone. The identifiability and reliability of the estimated parameters were analyzed with the Marcov chain Monte Carlo (MCMC) method. @All rights reserved. No part of the publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior permission of the author.

Catalytic gasification of glycerol in supercritical water

The conversion of glycerol in supercritical water (SCW) was studied at 510-550 °C and a pressure of 350 bars using both a bed of inert and non-porous ZrO2 particles (hydrothermal experiments), and a bed of a 1% Ru/ZrO2 catalyst. Experiments were conducted with a glycerol concentration of 5 wt% in a continuous isothermal fixed-bed reactor at a residence time between 2 and 10 s. Hydrothermolysis of glycerol formed water-soluble products such as acetaldehyde, acetic acid, hydroxyacetone and acrolein, and gases like H2, CO and CO2. The catalyst enhanced the formation of acetic acid, inhibited the formation of acrolein, and promoted gasification of the glycerol decomposition products. Hydrogen and carbon oxides were the main gases produced in the catalytic experiments, with minor amounts of methane and ethylene. Complete glycerol conversion was achieved at a residence time of 8.5 s at 510 °C, and at around 5 s at 550 °C with the 1 wt% Ru/ZrO2 catalyst. The catalyst was not active enough to achieve complete gasification since high yields of primary products like acetic acid and acetaldehyde were still present. Carbon balances were between 80 and 60% in the catalytic experiments, decreasing continuously as the residence time was increased. This was attributed partially to the formation of methanol and acetaldehyde, which were not recovered and analyzed efficiently in our set-up, but also to the formation of carbon deposits. Carbon deposition was not observed on the catalyst particles but on the surface of the inert zirconia particles, especially at high residence time. This was related to the higher concentration of acetic acid and other acidic species in the catalytic experiments, which may polymerize to form tar-like carbon precursors. Because of carbon deposition, hydrogen yields were significantly lower than expected; for instance at 550 °C the hydrogen yield potential was only 50% of the stoichiometric value.

Effects of cobalt-chromium everolimus eluting stents or bare metal stent on fatal and non-fatal cardiovascular events: patient level meta-analysis

Objectives: To examine the safety and effectiveness of cobalt-chromium everolimus eluting stents compared with bare metal stents. Design: Individual patient data meta-analysis of randomised controlled trials. Cox proportional regression models stratified by trial, containing random effects, were used to assess the impact of stent type on outcomes. Hazard ratios with 95% confidence interval for outcomes were reported. Data sources and study selection: Medline, Embase, the Cochrane Central Register of Controlled Trials. Randomised controlled trials that compared cobalt-chromium everolimus eluting stents with bare metal stents were selected. The principal investigators whose trials met the inclusion criteria provided data for individual patients. Primary outcomes: The primary outcome was cardiac mortality. Secondary endpoints were myocardial infarction, definite stent thrombosis, definite or probable stent thrombosis, target vessel revascularisation, and all cause death. Results: The search yielded five randomised controlled trials, comprising 4896 participants. Compared with patients receiving bare metal stents, participants receiving cobalt-chromium everolimus eluting stents had a significant reduction of cardiac mortality (hazard ratio 0.67, 95% confidence interval 0.49 to 0.91; P=0.01), myocardial infarction (0.71, 0.55 to 0.92; P=0.01), definite stent thrombosis (0.41, 0.22 to 0.76; P=0.005), definite or probable stent thrombosis (0.48, 0.31 to 0.73; P<0.001), and target vessel revascularisation (0.29, 0.20 to 0.41; P<0.001) at a median follow-up of 720 days. There was no significant difference in all cause death between groups (0.83, 0.65 to 1.06; P=0.14). Findings remained unchanged at multivariable regression after adjustment for the acuity of clinical syndrome (for instance, acute coronary syndrome v stable coronary artery disease), diabetes mellitus, female sex, use of glycoprotein IIb/IIIa inhibitors, and up to one year v longer duration treatment with dual antiplatelets. Conclusions: This meta-analysis offers evidence that compared with bare metal stents the use of cobalt-chromium everolimus eluting stents improves global cardiovascular outcomes including cardiac survival, myocardial infarction, and overall stent thrombosis.

Effects of cobalt-chromium everolimus eluting stents or bare metal stent on fatal and non-fatal cardiovascular events: patient level meta-analysis

Objectives: To examine the safety and effectiveness of cobalt-chromium everolimus eluting stents compared with bare metal stents. Design: Individual patient data meta-analysis of randomised controlled trials. Cox proportional regression models stratified by trial, containing random effects, were used to assess the impact of stent type on outcomes. Hazard ratios with 95% confidence interval for outcomes were reported. Data sources and study selection: Medline, Embase, the Cochrane Central Register of Controlled Trials. Randomised controlled trials that compared cobalt-chromium everolimus eluting stents with bare metal stents were selected. The principal investigators whose trials met the inclusion criteria provided data for individual patients. Primary outcomes: The primary outcome was cardiac mortality. Secondary endpoints were myocardial infarction, definite stent thrombosis, definite or probable stent thrombosis, target vessel revascularisation, and all cause death. Results: The search yielded five randomised controlled trials, comprising 4896 participants. Compared with patients receiving bare metal stents, participants receiving cobalt-chromium everolimus eluting stents had a significant reduction of cardiac mortality (hazard ratio 0.67, 95% confidence interval 0.49 to 0.91; P=0.01), myocardial infarction (0.71, 0.55 to 0.92; P=0.01), definite stent thrombosis (0.41, 0.22 to 0.76; P=0.005), definite or probable stent thrombosis (0.48, 0.31 to 0.73; P<0.001), and target vessel revascularisation (0.29, 0.20 to 0.41; P<0.001) at a median follow-up of 720 days. There was no significant difference in all cause death between groups (0.83, 0.65 to 1.06; P=0.14). Findings remained unchanged at multivariable regression after adjustment for the acuity of clinical syndrome (for instance, acute coronary syndrome v stable coronary artery disease), diabetes mellitus, female sex, use of glycoprotein IIb/IIIa inhibitors, and up to one year v longer duration treatment with dual antiplatelets. Conclusions: This meta-analysis offers evidence that compared with bare metal stents the use of cobalt-chromium everolimus eluting stents improves global cardiovascular outcomes including cardiac survival, myocardial infarction, and overall stent thrombosis.

Thermodynamic study of the solubility of sodium naproxen in some ethanol + water mixtures

By using the van't Hoff and Gibbs equations the apparent thermodynamic functions Gibbs energy, enthalpy, and entropy of solution for sodium naproxen in ethanol + water cosolvent mixtures, were evaluated from solubility data determined at temperatures from (278.15 to 308.15) K. The drug solubility was greatest in neat water and lowest in neat ethanol at all the temperatures studied. By means of non-linear enthalpy-entropy compensation analysis, it follows that the dissolution process of this drug in ethanol-rich mixtures is entropy-driven, whereas, in water-rich mixtures the process is enthalpy-driven.

Solution thermodynamics of methocarbamol in some ethanol + water mixtures

Apparent thermodynamic functions, Gibbs energy, enthalpy and entropy of solution and mixing, for methocarbamol in ethanol + water mixtures, were evaluated from solubility data determined at temperatures from 293.15 K to 313.15 K and from calorimetric values of drug fusion. The drug solubility was greatest in the mixtures with 0.70 or 0.80 mass fraction of ethanol and lowest in neat water across all temperatures studied. Non-linear enthalpy-entropy compensation was found for the dissolution processes. Accordingly, solution enthalpy drives the respective processes in almost all the solvent systems analyzed.