Association of sulfur dioxide exposure with circulatory system deaths in a medium-sized city in Brazil

There is a demonstrable association between exposure to air pollutants and deaths due to cardiovascular diseases. The objective of this study was to estimate the effects of exposure to sulfur dioxide on mortality due to circulatory diseases in individuals 50 years of age or older residing in São José dos Campos, SP. This was a time-series ecological study for the years 2003 to 2007 using information on deaths due to circulatory disease obtained from Datasus reports. Data on daily levels of pollutants, particulate matter, sulfur dioxide (SO2), ozone, temperature, and humidity were obtained from the São Paulo State Environmental Agency. Moving average models for 2 to 7 days were calculated by Poisson regression using the R software. Exposure to SO2 was analyzed using a unipollutant, bipollutant or multipollutant model adjusted for mean temperature and humidity. The relative risks with 95%CI were obtained and the percent decrease in risk was calculated. There were 1928 deaths with a daily mean (± SD) of 1.05 ± 1.03 (range: 0-6). Exposure to SO2 was significantly associated with mortality due to circulatory disease: RR = 1.04 (95%CI = 1.01 to 1.06) in the 7-day moving average, after adjusting for ozone. There was an 8.5% decrease in risk in the multipollutant model, proportional to a decrease of SO2 concentrations. The results of this study suggest that residents of medium-sized Brazilian cities with characteristics similar to those of São José dos Campos probably have health problems due to exposure to air pollutants.

Approximation de la distribution a posteriori d'un modèle Gamma-Poisson hiérarchique à effets mixtes

La méthode que nous présentons pour modéliser des données dites de "comptage" ou données de Poisson est basée sur la procédure nommée Modélisation multi-niveau et interactive de la régression de Poisson (PRIMM) développée par Christiansen et Morris (1997). Dans la méthode PRIMM, la régression de Poisson ne comprend que des effets fixes tandis que notre modèle intègre en plus des effets aléatoires. De même que Christiansen et Morris (1997), le modèle étudié consiste à faire de l'inférence basée sur des approximations analytiques des distributions a posteriori des paramètres, évitant ainsi d'utiliser des méthodes computationnelles comme les méthodes de Monte Carlo par chaînes de Markov (MCMC). Les approximations sont basées sur la méthode de Laplace et la théorie asymptotique liée à l'approximation normale pour les lois a posteriori. L'estimation des paramètres de la régression de Poisson est faite par la maximisation de leur densité a posteriori via l'algorithme de Newton-Raphson. Cette étude détermine également les deux premiers moments a posteriori des paramètres de la loi de Poisson dont la distribution a posteriori de chacun d'eux est approximativement une loi gamma. Des applications sur deux exemples de données ont permis de vérifier que ce modèle peut être considéré dans une certaine mesure comme une généralisation de la méthode PRIMM. En effet, le modèle s'applique aussi bien aux données de Poisson non stratifiées qu'aux données stratifiées; et dans ce dernier cas, il comporte non seulement des effets fixes mais aussi des effets aléatoires liés aux strates. Enfin, le modèle est appliqué aux données relatives à plusieurs types d'effets indésirables observés chez les participants d'un essai clinique impliquant un vaccin quadrivalent contre la rougeole, les oreillons, la rub\'eole et la varicelle. La régression de Poisson comprend l'effet fixe correspondant à la variable traitement/contrôle, ainsi que des effets aléatoires liés aux systèmes biologiques du corps humain auxquels sont attribués les effets indésirables considérés.

Caractérisation fonctionnelle chez le poisson zèbre de l'isoforme protéique WNK1/HSN2 mutée dans la neuropathie héréditaire sensitive et autonome de type 2

La neuropathie humaine sensitive et autonome de type 2 (NHSA 2) est une pathologie héréditaire rare caractérisée par une apparition précoce des symptômes et une absence d’affectation motrice. Cette pathologie entraîne la perte de perception de la douleur, de la chaleur et du froid ainsi que de la pression (toucher) dans les membres supérieurs et inférieurs et est due à des mutations autosomales récessives confinées à l’exon HSN2 de la protéine kinase à sérine/thréonine WNK1 (with-no-lysine protein kinase 1). Cet exon spécifique permettrait de conférer une spécificité au système nerveux à l’isoforme protéique WNK1/HSN2. La kinase WNK1 est étudiée en détails, en particulier au niveau du rein, mais son rôle au sein du système nerveux demeure inconnu. Considérant le début précoce de la neuropathie et le manque d’innervation sensorielle révélé par des biopsies chez les patients NHSA2, notre hypothèse de recherche est que les mutations tronquantes menant à la NHSA de type 2 causent une perte de fonction de l’isoforme WNK1/HSN2 spécifique au système nerveux entraînant un défaut dans le développement du système nerveux sensoriel périphérique. Chez l’embryon du poisson zèbre, WNK1/HSN2 est exprimé au niveau des neuromastes de la ligne latérale postérieure, un système mécanosensoriel périphérique. Nous avons obtenu des embryons knockdown pour WNK1/HSN2 par usage d’oligonucléotides morpholino antisens (AMO). Nos trois approches AMO ont révélé des embryons présentant des défauts d’établissement au niveau de la ligne latérale postérieure. Afin de déterminer la voie pathogène impliquant l’isoforme WNK1/HSN2, nous nous sommes intéressés à l’interaction rapportée entre la kinase WNK1 et le co-transporteur neuronal KCC2. Ce dernier est une cible de phosphorylation de WNK1 et son rôle dans la promotion de la neurogenèse est bien connu. Nous avons détecté l’expression de KCC2 au niveau de neuromastes de la ligne latérale postérieure et observé une expression accrue de KCC2 chez les embryons knockdown pour WNK1/HSN2 à l’aide de RT-PCR semi-quantitative. De plus, une sur-expression d’ARN humain de KCC2 chez des embryons a produit des défauts dans la ligne latérale postérieure, phénocopiant le knockdown de WNK1/HSN2. Ces résultats furent validés par un double knockdown, produisant des embryons n’exprimant ni KCC2, ni WNK1/HSN2, dont le phénotype fut atténué. Ces résultats nous mènent à suggérer une voie de signalisation où WNK1/HSN2 est en amont de KCC2, régulant son activation, et possiblement son expression. Nous proposons donc que la perte de fonction de l’isoforme spécifique cause un débalancement dans les niveaux de KCC2 activée, menant à une prolifération et une différenciation réduites des progéniteurs neuronaux du système nerveux périphérique. Les défauts associés à la NHSA de type 2 seraient donc de nature développementale et non neurodégénérative.

Approximate approximations for the Poisson and the Stokes equations

The method of approximate approximations is based on generating functions representing an approximate partition of the unity, only. In the present paper this method is used for the numerical solution of the Poisson equation and the Stokes system in R^n (n = 2, 3). The corresponding approximate volume potentials will be computed explicitly in these cases, containing a one-dimensional integral, only. Numerical simulations show the efficiency of the method and confirm the expected convergence of essentially second order, depending on the smoothness of the data.

On exact time averages of a massive Poisson particle

In this work we study, under the Stratonovich definition, the problem of the damped oscillatory massive particle subject to a heterogeneous Poisson noise characterized by a rate of events, lambda(t), and a magnitude, Phi, following an exponential distribution. We tackle the problem by performing exact time averages over the noise in a similar way to previous works analysing the problem of the Brownian particle. From this procedure we obtain the long-term equilibrium distributions of position and velocity as well as analytical asymptotic expressions for the injection and dissipation of energy terms. Considerations on the emergence of stochastic resonance in this type of system are also set forth.

Petrogenesis of Upper Paleozoic post-collisional peraluminous leucogranites, Sierra de Ancasti, northwest Argentina

The Santa Rosa and Sauce Guacho plutons are two post-collisional peraluminous Late Devonian to Early Carboniferous leucogranites that intruded the banded schists of the Ancasti Formation. The leucogranites are composed of microcline phenocrysts along with quartz, plagioclase, muscovite, biotite, ilmenite, tourmaline, apatite, monazite and zircon. Their geochemical composition is consistent with S-type granites and mineralogically they belong to MPG granites (muscovite-peraluminous granites). It is proposed that granite magma generation was related to shear zones that concentrated fluids in the metasedimentary crust during a collision or transcurrent tectonics. U-Pb analyses on monazite gave an age of 369.8 +/- 5.3 Ma, while Sm/Nd isotopic data yield epsilon(Nd(t)) values of -5.3 for Sauce Guacho and -5.7 for Santa Rosa indicating crustal provenance. Nd model ages between 1,544 and 1,571 Ma are within the range of magmatic rocks from the Lower Ordovician Famatinian Arc in the Central Sierras Pampeanas.

pH at the micellar interface: Synthesis of pH probes derived from salicylic acid, acid-base dissociation in sodium dodecyl sulfate micelles, and Poisson-Boltzmann simulation

The study of the H+ concentration at the micellar interface is a convenient system for modeling the distribution of H+ at interfaces. We have synthesized salicylic acid derivatives to analyze the proton dissociation of both the carboxylic and phenol groups of' the probes, determining spectrophotometrically the apparent pK(a)'s (pK(ap)) in sodium dodecyl Sulfate, SDS, micelles with and without added salt. The synthesized probes were 2-hydroxy-5-(2-trimethylammoniumacetyl)benzoate; 2-hydroxy-5-(2-dimethylhexadecylammoniumacetyl)benzoate- 2-hydroxy-5-(2-dimethylhexadecylammoniumhexanoyl)benzoate-, 2-hydroxy-5-(2-diniethylhexadecylammoniumundecanoyl)betizoate; 2-hydroxy-5-acetylbenzoic acids and 2-hydroxy-5-dodecanoylbenzoic acid. Upon incorporation into SDS micelles the pK(ap)'s of both carboxylic and phenol groups increased by ca. 3 pH units and NaCl addition caused a decrease in the probe-incorporated pKap. The experimental results were fitted with a cell model Poisson-Boltzmann (P-B) equation taking in consideration the effect of salt on the aggregation number of SDS and using the distance of' the dissociating group as a parameter. The conformations of the probes were analyzed theoretically using two dielectric constants, e.g., 2 and 78. Both the P-B analysis and conformation calculations can be interpreted by assuming that the acid groups dissociate very close to, or at, the interface. Our results are consistent with the assumption that the intrinsic pK(a)'s of both carboxylic and phenol groups of the salicylic acid probes used here can be taken as those in water. Using this assumption the micellar and salt effects on the pKap's of the (trialkylammonium)benzoate probes were described accurately using a cell model P-B analysis. (c) 2005 Elsevier B.V. All rights reserved.

Dynamical capture in the Pluto-Charon system

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

Aligning the CMS muon chambers with the muon alignment system during an extended cosmic ray run

The alignment system for the muon spectrometer of the CMS detector comprises three independent subsystems of optical and analog position sensors. It aligns muon chambers with respect to each other and to the central silicon tracker. System commissioning at full magnetic field began in 2008 during an extended cosmic ray run. The system succeeded in tracking muon detector movements of up to 18 mm and rotations of several milliradians under magnetic forces. Depending on coordinate and subsystem, the system achieved chamber alignment precisions of 140-350 μm and 30-200 μrad, close to the precision requirements of the experiment. Systematic errors on absolute positions are estimated to be 340-590 μm based on comparisons with independent photogrammetry measurements. © 2010 IOP Publishing Ltd and SISSA.

Measurement of masses in the tt̄ system by kinematic endpoints in pp collisions at √s = 7 TeV

A simultaneous measurement of the top-quark, W-boson, and neutrino masses is reported for tt̄ events selected in the dilepton final state from a data sample corresponding to an integrated luminosity of 5.0 fb-1 collected by the CMS experiment in pp collisions at √s = 7 TeV. The analysis is based on endpoint determinations in kinematic distributions. When the neutrino and W-boson masses are constrained to their world-average values, a top-quark mass value of Mt = 173.9 ± 0.9 (stat)+1.7 -2.1(syst.) GeV is obtained. When such constraints are not used, the three particle masses are obtained in a simultaneous fit. In this unconstrained mode the study serves as a test of mass determination methods that may be used in beyond standard model physics scenarios where several masses in a decay chain may be unknown and undetected particles lead to underconstrained kinematics. © 2013 CERN for the benefit of the CMS collaboration.

A Bayesian destructive weighted Poisson cure rate model and an application to a cutaneous melanoma data

In this article, we propose a new Bayesian flexible cure rate survival model, which generalises the stochastic model of Klebanov et al. [Klebanov LB, Rachev ST and Yakovlev AY. A stochastic-model of radiation carcinogenesis - latent time distributions and their properties. Math Biosci 1993; 113: 51-75], and has much in common with the destructive model formulated by Rodrigues et al. [Rodrigues J, de Castro M, Balakrishnan N and Cancho VG. Destructive weighted Poisson cure rate models. Technical Report, Universidade Federal de Sao Carlos, Sao Carlos-SP. Brazil, 2009 (accepted in Lifetime Data Analysis)]. In our approach, the accumulated number of lesions or altered cells follows a compound weighted Poisson distribution. This model is more flexible than the promotion time cure model in terms of dispersion. Moreover, it possesses an interesting and realistic interpretation of the biological mechanism of the occurrence of the event of interest as it includes a destructive process of tumour cells after an initial treatment or the capacity of an individual exposed to irradiation to repair altered cells that results in cancer induction. In other words, what is recorded is only the damaged portion of the original number of altered cells not eliminated by the treatment or repaired by the repair system of an individual. Markov Chain Monte Carlo (MCMC) methods are then used to develop Bayesian inference for the proposed model. Also, some discussions on the model selection and an illustration with a cutaneous melanoma data set analysed by Rodrigues et al. [Rodrigues J, de Castro M, Balakrishnan N and Cancho VG. Destructive weighted Poisson cure rate models. Technical Report, Universidade Federal de Sao Carlos, Sao Carlos-SP. Brazil, 2009 (accepted in Lifetime Data Analysis)] are presented.

Ausbreitung einer Infektion durch ein System Brownscher Bewegungen

In dieser Arbeit werden wir ein Modell untersuchen, welches die Ausbreitung einer Infektion beschreibt. Bei diesem Modell werden zunächst Partikel gemäß eines Poissonschen Punktprozesses auf der reellen Achse verteilt. Bis zu einem gewissen Punkt auf der reellen Achse sind alle Partikel von einer Infektion befallen. Während sich nicht infizierte Partikel nicht bewegen, folgen die infizierten Partikel den Pfaden von voneinander unabhängigen Brownschen Bewegungen und verbreitet die Infektion dabei an den Orten, welche sie betreten. Wenn sie dabei auf ein nicht infiziertes Partikel treffen, ist dieses von diesem Moment an auch infiziert und beginnt ebenfalls, dem Pfad einer Brownschen Bewegung zu folgen und die Infektion auszubreiten. Auf diese Art verschiebt sich nun der am weitesten rechts liegende Ort R_t, an dem die Infektion bereits verbreitet wurde. Wir werden mit Hilfe des subadditiven Ergodensatzes zeigen, dass sich dieser Ort mit linearer Geschwindigkeit fortbewegt. Ferner werden wir eine obere und eine untere Schranke für die Ausbreitungsgeschwindkeit angeben. Danach werden wir zeigen, dass der Prozess Regenerationszeiten hat, nämlich solche zufällige Zeiten, zu denen er eine Art Neustart unter speziellen Startbedingungen durchführt. Wir werden diese für eine weitere Charakterisierung der Ausbreitungsgeschwingkeit nutzen. Ferner erhalten wir durch die Regenerationszeiten auch einen Zentralen Grenzwertsatz für R_t und können zeigen, dass die Verteilung der infizierten Partikel aus Sicht des am weitesten rechts liegenden infizierten Ortes gegen eine invariante Verteilung konvergiert.

Low Work-Function Thermionic Emission and Orbital-Motion-Limited Ion Collection at Bare-Tether Cathodic Contact

Una amarra electrodinámica (electrodynamic tether) opera sobre principios electromagnéticos intercambiando momento con la magnetosfera planetaria e interactuando con su ionosfera. Es un subsistema pasivo fiable para desorbitar etapas de cohetes agotadas y satélites al final de su misión, mitigando el crecimiento de la basura espacial. Una amarra sin aislamiento captura electrones del plasma ambiente a lo largo de su segmento polarizado positivamente, el cual puede alcanzar varios kilómetros de longitud, mientras que emite electrones de vuelta al plasma mediante un contactor de plasma activo de baja impedancia en su extremo catódico, tal como un cátodo hueco (hollow cathode). En ausencia de un contactor catódico activo, la corriente que circula por una amarra desnuda en órbita es nula en ambos extremos de la amarra y se dice que ésta está flotando eléctricamente. Para emisión termoiónica despreciable y captura de corriente en condiciones limitadas por movimiento orbital (orbital-motion-limited, OML), el cociente entre las longitudes de los segmentos anódico y catódico es muy pequeño debido a la disparidad de masas entre iones y electrones. Tal modo de operación resulta en una corriente media y fuerza de Lorentz bajas en la amarra, la cual es poco eficiente como dispositivo para desorbitar. El electride C12A7 : e−, que podría presentar una función de trabajo (work function) tan baja como W = 0.6 eV y un comportamiento estable a temperaturas relativamente altas, ha sido propuesto como recubrimiento para amarras desnudas. La emisión termoiónica a lo largo de un segmento así recubierto y bajo el calentamiento de la operación espacial, puede ser más eficiente que la captura iónica. En el modo más simple de fuerza de frenado, podría eliminar la necesidad de un contactor catódico activo y su correspondientes requisitos de alimentación de gas y subsistema de potencia, lo que resultaría en un sistema real de amarra “sin combustible”. Con este recubrimiento de bajo W, cada segmento elemental del segmento catódico de una amarra desnuda de kilómetros de longitud emitiría corriente como si fuese parte de una sonda cilíndrica, caliente y uniformemente polarizada al potencial local de la amarra. La operación es similar a la de una sonda de Langmuir 2D tanto en los segmentos catódico como anódico. Sin embargo, en presencia de emisión, los electrones emitidos resultan en carga espacial (space charge) negativa, la cual reduce el campo eléctrico que los acelera hacia fuera, o incluso puede desacelerarlos y hacerlos volver a la sonda. Se forma una doble vainas (double sheath) estable con electrones emitidos desde la sonda e iones provenientes del plasma ambiente. La densidad de corriente termoiónica, variando a lo largo del segmento catódico, podría seguir dos leyes distintas bajo diferentes condiciones: (i) la ley de corriente limitada por la carga espacial (space-charge-limited, SCL) o (ii) la ley de Richardson-Dushman (RDS). Se presenta un estudio preliminar sobre la corriente SCL frente a una sonda emisora usando la teoría de vainas (sheath) formada por la captura iónica en condiciones OML, y la corriente electrónica SCL entre los electrodos cilíndricos según Langmuir. El modelo, que incluye efectos óhmicos y el efecto de transición de emisión SCL a emisión RDS, proporciona los perfiles de corriente y potencial a lo largo de la longitud completa de la amarra. El análisis muestra que en el modo más simple de fuerza de frenado, bajo condiciones orbitales y de amarras típicas, la emisión termoiónica proporciona un contacto catódico eficiente y resulta en una sección catódica pequeña. En el análisis anterior, tanto la transición de emisión SCL a RD como la propia ley de emisión SCL consiste en un modelo muy simplificado. Por ello, a continuación se ha estudiado con detalle la solución de vaina estacionaria de una sonda con emisión termoiónica polarizada negativamente respecto a un plasma isotrópico, no colisional y sin campo magnético. La existencia de posibles partículas atrapadas ha sido ignorada y el estudio incluye tanto un estudio semi-analítico mediante técnica asintóticas como soluciones numéricas completas del problema. Bajo las tres condiciones (i) alto potencial, (ii) R = Rmax para la validez de la captura iónica OML, y (iii) potencial monotónico, se desarrolla un análisis asintótico auto-consistente para la estructura de plasma compleja que contiene las tres especies de cargas (electrones e iones del plasma, electrones emitidos), y cuatro regiones espaciales distintas, utilizando teorías de movimiento orbital y modelos cinéticos de las especies. Aunque los electrones emitidos presentan carga espacial despreciable muy lejos de la sonda, su efecto no se puede despreciar en el análisis global de la estructura de la vaina y de dos capas finas entre la vaina y la región cuasi-neutra. El análisis proporciona las condiciones paramétricas para que la corriente sea SCL. También muestra que la emisión termoiónica aumenta el radio máximo de la sonda para operar dentro del régimen OML y que la emisión de electrones es mucho más eficiente que la captura iónica para el segmento catódico de la amarra. En el código numérico, los movimientos orbitales de las tres especies son modelados para potenciales tanto monotónico como no-monotónico, y sonda de radio R arbitrario (dentro o más allá del régimen de OML para la captura iónica). Aprovechando la existencia de dos invariante, el sistema de ecuaciones Poisson-Vlasov se escribe como una ecuación integro-diferencial, la cual se discretiza mediante un método de diferencias finitas. El sistema de ecuaciones algebraicas no lineal resultante se ha resuelto de con un método Newton-Raphson paralelizado. Los resultados, comparados satisfactoriamente con el análisis analítico, proporcionan la emisión de corriente y la estructura del plasma y del potencial electrostático. ABSTRACT An electrodynamic tether operates on electromagnetic principles and exchanges momentum through the planetary magnetosphere, by continuously interacting with the ionosphere. It is a reliable passive subsystem to deorbit spent rocket stages and satellites at its end of mission, mitigating the growth of orbital debris. A tether left bare of insulation collects electrons by its own uninsulated and positively biased segment with kilometer range, while electrons are emitted by a low-impedance active device at the cathodic end, such as a hollow cathode, to emit the full electron current. In the absence of an active cathodic device, the current flowing along an orbiting bare tether vanishes at both ends and the tether is said to be electrically floating. For negligible thermionic emission and orbital-motion-limited (OML) collection throughout the entire tether (electron/ion collection at anodic/cathodic segment, respectively), the anodic-to-cathodic length ratio is very small due to ions being much heavier, which results in low average current and Lorentz drag. The electride C12A7 : e−, which might present a possible work function as low as W = 0.6 eV and moderately high temperature stability, has been proposed as coating for floating bare tethers. Thermionic emission along a thus coated cathodic segment, under heating in space operation, can be more efficient than ion collection and, in the simplest drag mode, may eliminate the need for an active cathodic device and its corresponding gas-feed requirements and power subsystem, which would result in a truly “propellant-less” tether system. With this low-W coating, each elemental segment on the cathodic segment of a kilometers-long floating bare-tether would emit current as if it were part of a hot cylindrical probe uniformly polarized at the local tether bias, under 2D probe conditions that are also applied to the anodic-segment analysis. In the presence of emission, emitted electrons result in negative space charge, which decreases the electric field that accelerates them outwards, or even reverses it, decelerating electrons near the emitting probe. A double sheath would be established with electrons being emitted from the probe and ions coming from the ambient plasma. The thermionic current density, varying along the cathodic segment, might follow two distinct laws under different con ditions: i) space-charge-limited (SCL) emission or ii) full Richardson-Dushman (RDS) emission. A preliminary study on the SCL current in front of an emissive probe is presented using the orbital-motion-limited (OML) ion-collection sheath and Langmuir’s SCL electron current between cylindrical electrodes. A detailed calculation of current and bias profiles along the entire tether length is carried out with ohmic effects considered and the transition from SCL to full RDS emission is included. Analysis shows that in the simplest drag mode, under typical orbital and tether conditions, thermionic emission provides efficient cathodic contact and leads to a short cathodic section. In the previous analysis, both the transition between SCL and RDS emission and the current law for SCL condition have used a very simple model. To continue, considering an isotropic, unmagnetized, colissionless plasma and a stationary sheath, the probe-plasma contact is studied in detail for a negatively biased probe with thermionic emission. The possible trapped particles are ignored and this study includes both semianalytical solutions using asymptotic analysis and complete numerical solutions. Under conditions of i) high bias, ii) R = Rmax for ion OML collection validity, and iii) monotonic potential, a self-consistent asymptotic analysis is carried out for the complex plasma structure involving all three charge species (plasma electrons and ions, and emitted electrons) and four distinct spatial regions using orbital motion theories and kinetic modeling of the species. Although emitted electrons present negligible space charge far away from the probe, their effect cannot be neglected in the global analysis for the sheath structure and two thin layers in between the sheath and the quasineutral region. The parametric conditions for the current to be space-chargelimited are obtained. It is found that thermionic emission increases the range of probe radius for OML validity and is greatly more effective than ion collection for cathodic contact of tethers. In the numerical code, the orbital motions of all three species are modeled for both monotonic and non-monotonic potential, and for any probe radius R (within or beyond OML regime for ion collection). Taking advantage of two constants of motion (energy and angular momentum), the Poisson-Vlasov equation is described by an integro differential equation, which is discretized using finite difference method. The non-linear algebraic equations are solved using a parallel implementation of the Newton-Raphson method. The results, which show good agreement with the analytical results, provide the results for thermionic current, the sheath structure, and the electrostatic potential.

Operational short term health impact assessment of air pollution modelling system over Europe

The last decade, scientific studies have indicated an association between air pollution to which people are exposed and wide range of adverse health outcomes. We have developed a tool which is based on a model (MM5-CMAQ) running over Europe with 50 km spatial resolution, based on EMEP annual emissions, to produce a short-term forecast of the impact on health. In order to estimate the mortality change (forecasted for the next 24 hours) we have chosen a log-linear (Poisson) regression form to estimate the concentration-response function. The parameters involved in the C-R function have been estimated based on epidemiological studies, which have been published. Finally, we have derived the relationship between concentration change and mortality change from the C-R function which is the final health impact function.

Single-Server Queueing System with Markov-Modulated Arrivals and Service Times

Key words: Markov-modulated queues, waiting time, heavy traffic.