Chemistry, isotope values (δD, δ18O, δ34S(SO4)) and temperatures of the water inflows in two Gotthard tunnels, Swiss Alps

Water inflows in the Gotthard Highway Tunnel and in the Gotthard Exploration Tunnel are meteoric waters infiltrating at different elevations, on both sides of an important orographic divide. Limited interaction of meteoric waters with gneissic rocks produces Ca-HCO3 and Na-Ca-HCO3 waters, whereas prolonged interaction of meteoric waters with the same rocks generates Na-HCO3 to Na-SO4 waters. Waters circulating in Triassic carbonate-evaporite rocks have a Ca-SO4 composition. Calcium-Na-SO4 waters are also present. They can be produced through interaction of either Na-HCO3 waters with anhydrite or Ca-SO4 waters with a local gneissic rock, as suggested by reaction path modeling. An analogous simulation indicates that Na-HCO3 waters are generated through interaction of Ca-HCO3 waters with a local gneissic rock. The two main SO4-sources present in the Alps are leaching of upper Triassic sulfate minerals and oxidative dissolution of sulfide minerals of crystalline rocks. Values of delta S-34(SO4) < <similar to>+ 9 parts per thousand, are due to oxidative dissolution of sulfide minerals, whereas delta S-34(SO4) > similar to+ 9 parts per thousand are controlled either by bacterial SO4 reduction or leaching of upper Triassic sulfate minerals. Most waters have temperatures similar to the expected values for a geothermal gradient of 22 degreesC/km and are close to thermal equilibrium with rocks. However relatively large, descending flows of cold waters and ascending flows of warm waters are present in both tunnels and determine substantial cooling and heating, respectively, of the interacting rocks. The most import upflow zone of warm, Na-rich waters is below Guspisbach, in the Gotthard Highway Tunnel, at 6.2-9.0 km from the southern portal. These warm waters have equilibrium temperatures of 65-75 degreesC and therefore constitute an important low-enthalpy geothermal resource. (C) 2001 Elsevier Science Ltd. All rights reserved.

Thermal energetics of the New-Guinean moss-forest rat (Rattus niobe) in comparison with other tropical murid rodents.

The thermal energetics of rodents from cool, wet tropical highlands are poorly known. Metabolic rate, body temperature and thermal conductance were measured in the moss-forest rat, Rattus niobe (Rodentia), a small murid endemic to the highlands of New Guinea. These data were evaluated in the context of the variation observed in the genus Rattus and among tropical murids. In 7 adult R. niobe, basal metabolic rate (BMR) averaged 53.6±6.6mLO2h(-1), or 103% of the value predicted for a body mass of 42.3±5.8g. Compared to other species of Rattus, R. niobe combines a low body temperature (35.5±0.6°C) and a moderately low minimal wet thermal conductance cmin (5.88±0.7mLO2h(-1)°C(-1), 95% of predicted) with a small size, all of which lead to reduced energy expenditure in a constantly cool environment. The correlations of mean annual rainfall and temperature, altitude and body mass with BMR, body temperature and cmin were analyzed comparatively among tropical Muridae. Neither BMR, nor cmin or body temperature correlated with ambient temperature or altitude. Some of the factors which promote high BMR in higher latitude habitats, such as seasonal exposure to very low temperature and short reproductive season, are lacking in wet montane tropical forests. BMR increased with rainfall, confirming a pattern observed among other assemblages of mammals. This correlation was due to the low BMR of several desert adapted murids, while R. niobe and other species from wet habitats had a moderate BMR.

Water chemistry and isotope composition of the Acquarossa thermal system, Ticino, Switzerland

The thermal springs of Acquarossa and the nearby mineral springs of Soia have outlet temperatures of 12 degrees to 25 degrees C, TDS of 2290 to 3000 mg/kg and Ca-SO4 to Ca-SO4-HCO3 composition. Chemical geothermometers suggest reservoir temperatures close to 60 degrees C. P-CO2 values at depth are estimated to range from 0.3 to 2 bar. delta D and delta(18)O values indicate a meteoric origin and recharge elevations of 1600 +/- 150 m above sea level (a.s.l.) for these thermal and mineral waters. All these waters discharge from the overturned limb of the Simano nappe, probably dose to the contact between basement and underlying cover rocks. They therefore represent rain waters that descend slowly, heat at depth and locally rise relatively quickly to the surface, preserving part of their physical and chemical characteristics. (C) 1999 CNR. Published by Elsevier Science Ltd. All rights reserved.

SIMULIT : un modèle de simulation pour l'analyse et la planification de l'activité hospitalière

Le modèle développé à l'Institut universitaire de médecine sociale et préventive de Lausanne utilise un programme informatique pour simuler les mouvements d'entrées et de sorties des hôpitaux de soins généraux. Cette simulation se fonde sur les données récoltées de routine dans les hôpitaux; elle tient notamment compte de certaines variations journalières et saisonnières, du nombre d'entrées, ainsi que du "Case-Mix" de l'hôpital, c'est-à-dire de la répartition des cas selon les groupes cliniques et l'âge des patients.

Intravenous streptomycin dosing regimen in a patient undergoing hemodialysis: Plasma level monitoring and pharmacokinetic simulation

Introduction: Streptomycin, as other aminoglycosides, exhibits concentration-dependent bacterial killing but has a narrow therapeutic window. It is primarily eliminated unchanged by the kidneys. Data and dosing information to achieve a safe regimen in patients with chronic renal failure undergoing hemodialysis (HD) are scarce. Although main adverse reactions are related to prolonged, elevated serum concentrations, literature recommendation is to administer streptomycin after each HD. Patients (or Materials) and Methods: We report the case of a patient with end-stage renal failure, undergoing HD, who was successfully treated with streptomycin for gentamicin-resistant Enterococcus faecalis bacteremia with prosthetic arteriovenous fistula infection. Streptomycin was administered intravenously 7.5 mg/kg, 3 hours before each dialysis (3 times a week) during 6 weeks in combination with amoxicillin. Streptomycin plasma levels were monitored with repeated blood sampling before, after, and between HD sessions. A 2-compartment model was used to reconstruct the concentration time profile over days on and off HD. Results: Streptomycin trough plasma-concentration was 2.8 mg/L. It peaked to 21.4 mg/L 30 minutes after intravenous administration, decreased to 18.2 mg/L immediately before HD, and dropped to 4.5 mg/L at the end of a 4-hour HD session. Plasma level increased again to 5.7 mg/L 2 hours after the end of HD and was 2.8 mg/L 48 hours later, before the next administration and HD. The pharmacokinetics of streptomycin was best described with a 2-compartment model. The computer simulation fitted fairly well to the observed concentrations during or between HD sessions. Redistribution between the 2 compartments after the end of HD reproduced the rebound of plasma concentrations after HD. No significant toxicity was observed during treatment. The outcome of the infection was favorable, and no sign of relapse was observed after a follow-up of 3 months. Conclusion: Streptomycin administration of 7.5 mg/kg 3 hours before HD sessions in a patient with end-stage renal failure resulted in an effective and safe dosing regimen. Monitoring plasma levels along with pharmacokinetic simulation document the suitability of this dosing scheme, which should replace current dosage recommendations for streptomycin in HD.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Tax simulation in the SHP

This paper describes a simulation package designed to estimate the annual income taxes paid by respondents of the Swiss Household Panel (SHP). In Switzerland, the 26 cantons have their own tax system. Additionally, tax levels vary between the over 2000 municipalities and over time. The simulation package takes account of this complexity by building on existing tables on tax levels which are provided by the Swiss Federal Tax Administration Office. Because these are limited to a few types of households and only 812 municipalities, they have to be extended to cover all households and municipalities. A further drawback of these tables is that they neglect several deductions. The tax simulation package fills this gap by taking additionally account of deductions for children, double-earner couples, third pillar and support for dependent persons according to cantonal legislation. The resulting variable on direct taxes not only serves to calculate household income net of taxes, but can also be a variable for analysis by its own account.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Nuclear well logging in soils. Calibration of slim hole nuclear tools gamma-gamma and neutron-thermal neutron.

The calculation of elasticity parameters by sonic and ultra sonic wave propagation in saturated soils using Biot's theory needs the following variables : forpiation density and porosity (p, ø), compressional and shear wave velocities (Vp, Vs), fluid density, viscosity and compressibility (Pfi Ilfi Ki), matrix density and compressibility (p" K), The first four parameters can be determined in situ using logging probes. Because fluid and matrix characteristics are not modified during core extraction, they can be obtained through laboratory measurements. All parameters necessitate precise calibrations in various environments and for specific range of values encountered in soils. The slim diameter of boreholes in shallow geophysics and the high cost of petroleum equipment demand the use of specific probes, which usually only give qualitative results. The measurement 'of density is done with a gamma-gamma probe and the measurement of hydrogen index, in relation to porosity, by a neutron probe. The first step of this work has been carried out in synthetic formations in the laboratory using homogeneous media of known density and porosity. To establish borehole corrections different casings have been used. Finally a comparison between laboratory and in situ data in cored holes of known geometry and casing has been performed.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.