Predictive analysis and mapping of indoor radon concentrations in Switzerland

It is estimated that around 230 people die each year due to radon (222Rn) exposure in Switzerland. 222Rn occurs mainly in closed environments like buildings and originates primarily from the subjacent ground. Therefore it depends strongly on geology and shows substantial regional variations. Correct identification of these regional variations would lead to substantial reduction of 222Rn exposure of the population based on appropriate construction of new and mitigation of already existing buildings. Prediction of indoor 222Rn concentrations (IRC) and identification of 222Rn prone areas is however difficult since IRC depend on a variety of different variables like building characteristics, meteorology, geology and anthropogenic factors. The present work aims at the development of predictive models and the understanding of IRC in Switzerland, taking into account a maximum of information in order to minimize the prediction uncertainty. The predictive maps will be used as a decision-support tool for 222Rn risk management. The construction of these models is based on different data-driven statistical methods, in combination with geographical information systems (GIS). In a first phase we performed univariate analysis of IRC for different variables, namely the detector type, building category, foundation, year of construction, the average outdoor temperature during measurement, altitude and lithology. All variables showed significant associations to IRC. Buildings constructed after 1900 showed significantly lower IRC compared to earlier constructions. We observed a further drop of IRC after 1970. In addition to that, we found an association of IRC with altitude. With regard to lithology, we observed the lowest IRC in sedimentary rocks (excluding carbonates) and sediments and the highest IRC in the Jura carbonates and igneous rock. The IRC data was systematically analyzed for potential bias due to spatially unbalanced sampling of measurements. In order to facilitate the modeling and the interpretation of the influence of geology on IRC, we developed an algorithm based on k-medoids clustering which permits to define coherent geological classes in terms of IRC. We performed a soil gas 222Rn concentration (SRC) measurement campaign in order to determine the predictive power of SRC with respect to IRC. We found that the use of SRC is limited for IRC prediction. The second part of the project was dedicated to predictive mapping of IRC using models which take into account the multidimensionality of the process of 222Rn entry into buildings. We used kernel regression and ensemble regression tree for this purpose. We could explain up to 33% of the variance of the log transformed IRC all over Switzerland. This is a good performance compared to former attempts of IRC modeling in Switzerland. As predictor variables we considered geographical coordinates, altitude, outdoor temperature, building type, foundation, year of construction and detector type. Ensemble regression trees like random forests allow to determine the role of each IRC predictor in a multidimensional setting. We found spatial information like geology, altitude and coordinates to have stronger influences on IRC than building related variables like foundation type, building type and year of construction. Based on kernel estimation we developed an approach to determine the local probability of IRC to exceed 300 Bq/m3. In addition to that we developed a confidence index in order to provide an estimate of uncertainty of the map. All methods allow an easy creation of tailor-made maps for different building characteristics. Our work is an essential step towards a 222Rn risk assessment which accounts at the same time for different architectural situations as well as geological and geographical conditions. For the communication of 222Rn hazard to the population we recommend to make use of the probability map based on kernel estimation. The communication of 222Rn hazard could for example be implemented via a web interface where the users specify the characteristics and coordinates of their home in order to obtain the probability to be above a given IRC with a corresponding index of confidence. Taking into account the health effects of 222Rn, our results have the potential to substantially improve the estimation of the effective dose from 222Rn delivered to the Swiss population.

Improved predictive mapping of indoor radon concentrations using ensemble regression trees based on automatic clustering of geological units.

PURPOSE: According to estimations around 230 people die as a result of radon exposure in Switzerland. This public health concern makes reliable indoor radon prediction and mapping methods necessary in order to improve risk communication to the public. The aim of this study was to develop an automated method to classify lithological units according to their radon characteristics and to develop mapping and predictive tools in order to improve local radon prediction. METHOD: About 240 000 indoor radon concentration (IRC) measurements in about 150 000 buildings were available for our analysis. The automated classification of lithological units was based on k-medoids clustering via pair-wise Kolmogorov distances between IRC distributions of lithological units. For IRC mapping and prediction we used random forests and Bayesian additive regression trees (BART). RESULTS: The automated classification groups lithological units well in terms of their IRC characteristics. Especially the IRC differences in metamorphic rocks like gneiss are well revealed by this method. The maps produced by random forests soundly represent the regional difference of IRCs in Switzerland and improve the spatial detail compared to existing approaches. We could explain 33% of the variations in IRC data with random forests. Additionally, the influence of a variable evaluated by random forests shows that building characteristics are less important predictors for IRCs than spatial/geological influences. BART could explain 29% of IRC variability and produced maps that indicate the prediction uncertainty. CONCLUSION: Ensemble regression trees are a powerful tool to model and understand the multidimensional influences on IRCs. Automatic clustering of lithological units complements this method by facilitating the interpretation of radon properties of rock types. This study provides an important element for radon risk communication. Future approaches should consider taking into account further variables like soil gas radon measurements as well as more detailed geological information.

Soil diffuse degassing and thermal energy fluxes from the southern Lakki plain, Nisyros (Greece)

Two diffuse soil CO2 flux surveys from the southern Lakki plain show that CO2 is mainly released from the hydrothermal explosion craters. The correspondence between high CO2 fluxes and elevated soil temperatures suggests that a flux of hot hydrothermal fluids ascends towards the surface. Steam mostly condenses near the surface and the heat given off is conductively transferred to the atmosphere through the soil, accompanied by a large CO2 flux. Tt was calculated, that 68 t d(-1) of hydrothermal CO2 are released through the total surveyed area of similar to1.3 km(2) Admitting that a steam flux of 2200 t d(-1) accompanies this CO2 flux, the thermal energy released through steam condensation amounts to 58 MW.

Derivation of soil-specific streaming potential electrical parameters from hydrodynamic characteristics of partially saturated soils

Water movement in unsaturated soils gives rise to measurable electrical potential differences that are related to the flow direction and volumetric fluxes, as well as to the soil properties themselves. Laboratory and field data suggest that these so-called streaming potentials may be several orders of magnitudes larger than theoretical predictions that only consider the influence of the relative permeability and electrical conductivity on the self potential (SP) data. Recent work has improved predictions somewhat by considering how the volumetric excess charge in the pore space scales with the inverse of water saturation. We present a new theoretical approach that uses the flux-averaged excess charge, not the volumetric excess charge, to predict streaming potentials. We present relationships for how this effective excess charge varies with water saturation for typical soil properties using either the water retention or the relative permeability function. We find large differences between soil types and the predictions based on the relative permeability function display the best agreement with field data. The new relationships better explain laboratory data than previous work and allow us to predict the recorded magnitudes of the streaming potentials following a rainfall event in sandy loam, whereas previous models predict values that are three orders of magnitude too small. We suggest that the strong signals in unsaturated media can be used to gain information about fluxes (including very small ones related to film flow), but also to constrain the relative permeability function, the water retention curve, and the relative electrical conductivity function.

Radiocarbon dating of modern groundwater: the role of the unsaturated zone

Biological and physical processes occurring in soils may lead to significant isotopic changes between the isotopic compositions of atmospheric CO2 and of soil CO2. Also, during water and gas transport from the soil surface to the water table, isotopic changes likely occur due to numerous physical processes such as gas production and diffusion, water advection, and gas-water-rock interactions. In most cases, these changes are not included in the correction models developed for groundwater dating, whereas they can significantly impact the calculation of the 14C age. We explore the role of these processes using: i) experimental data from two aquifer sites (Fontainebleau sands and Astian sands, France), ii) a distributed model to simulate the 14C activities of soil CO2, and iii) numerical simulations in order to highlight the role of the physical processes.¦The 13C content in soil CO2 showed seasonal variations and highlighted the competition between CO2 production and CO2 diffusion. Their respective contributions played a significant role in defining the isotopic composition of CO2 at the water table. On both study sites, variations of the 14C activity in soil CO2 reflect the competition between the fluxes of root derived-CO2 and organic matter derived-CO2. Since the nuclear weapon tests in the fifties and sixties, soil CO2 became significantly depleted in 14C compared to modern atmospheric CO2. Models that take into account this 14C depletion in soil CO2 for dating modern groundwater would lead to apparent younger 14C ages than models that only consider the 14C activity in atmospheric CO2. Moreover, since 2000-2005, the inverse effect is observed as soil CO2 is enriched in 14C compared to atmospheric CO2.¦Therefore, we conclude that the isotopic composition of CO2 at the water table have to be taken into account for the dating of modern groundwater. This requires a systematic sampling of soil CO2 and the measurement of its 13C and 14C contents. We used this information in a numerical simulation to calculate the evolution of isotopic composition of CO2 from the soil surface to the water table. This simulation integrated physical processes in the unsaturated zone (e.g. CO2 production and diffusion, water advection, etc.) and gas-water-rock interactions.

The role of soil in vegetated gravelly river braid plains: more than just a passive response?

This paper reviews the role of alluvial soils in vegetated gravelly river braid plains. When considering decadal time scales of river evolution, we argue that it becomes vital to consider soil development as an emergent property of the developing ecosystem. Soil processes have been relatively overlooked in accounts of the interactions between braided river processes and vegetation, although soils have been observed on vegetated fluvial landforms. We hypothesise that soil development plays a major role in the transition (speed and pathway) from a fresh sediment deposit to a vegetated soil-covered landform. Disturbance (erosion and/or deposition), vertical sediment structure (process history), vegetation succession, biological activity and water table fluctuation are seen as the main controls on early alluvial soil evolution. Erosion and deposition processes may not only act as soil disturbing agents, but also as suppliers of ecosystem resources, because of their role in delivering and changing access (e.g. through avulsion) to fluxes of water, fine sediments and organic matter. In turn, the associated initial ecosystem may influence further fluvial landform development, such as through the trapping of fine-grained sediments (e.g. sand) by the engineering action of vegetation and the deposit stabilisation by the developing above and belowground biomass. This may create a strong feedback between geomorphological processes, vegetation succession and soil evolution which we summarise in a conceptual model. We illustrate this model by an example from the Allondon River (CH) and identify the research questions that follow.

Open to closed system transition traced through the TDIC isotopic signature at the aquifer recharge stage, implications for groundwater 14C dating

14C dating models are limited when considering recent groundwater for which the carbon isotopic signature of the total dissolved inorganic carbon (TDIC) is mainly acquired in the unsaturated zone. Reducing the uncertainties of dating thus implies a better identification of the processes controlling the carbon isotopic composition of the TDIC during groundwater recharge. Geochemical interactions between gas, water and carbonates in the unsaturated zone were investigated for two aquifers (the carbonate-free Fontainebleau sands and carbonate-bearing Astian sands, France) in order to identify the respective roles of CO2 and carbonates on the carbon isotopic signatures of the TDIC; this analysis is usually approached using open or closed system terms. Under fully open system conditions, the seasonality of the 13C values in the soil CO2 can lead to important uncertainties regarding the so-called "initial 14C activity" used in 14C correction models. In a carbonate-bearing unsaturated zone such as in the Astian aquifer, we show that an approach based on fully open or closed system conditions is not appropriate. Although the chemical saturation between water and calcite occurs rapidly within the first metre of the unsaturated zone, the carbon isotopic contents (δ13C) of the CO2 and the TDIC evolve downward, impacted by the dissolution-precipitation of the carbonates. In this study, we propose a numerical approach to describe this evolution. The δ13C and the A 14C (radiocarbon activity) of the TDIC at the base of the carbonate-hearing unsaturated zone depends on (i) the δ13C and the A 14C of the TDIC in the soil determined by the soil CO2, (ii) the water's residence time in the unsaturated zone and (iii) the carbonate precipitation-dissolution fluxes. In this type of situation, the carbonate δ13C-A 14C evolutions indicate the presence of secondary calcite and permit the calculation of its accretion flux, equal to ~ 4.5 ± 0.5 x 10-9 mol grock-1 yr-1. More generally, for other sites under temperate climate and with similar properties to the Astian sands site, this approach allows for a reliable determination of the carbon isotopic composition at the base of the unsaturated zone as the indispensable "input function" data of the carbon cycle into the aquifer.

A new green fluorescent protein-based bacterial biosensor for analysing phenanthrene fluxes.

The polycyclic aromatic hydrocarbon (PAH)-degrading strain Burkholderia sp. RP007 served as host strain for the design of a bacterial biosensor for the detection of phenanthrene. RP007 was transformed with a reporter plasmid containing a transcriptional fusion between the phnS putative promoter/operator region and the gene encoding the enhanced green fluorescent protein (GFP). The resulting bacterial biosensor--Burkholderia sp. strain RP037--produced significant amounts of GFP after batch incubation in the presence of phenanthrene crystals. Co-incubation with acetate did not disturb the phenanthrene-specific response but resulted in a homogenously responding population of cells. Active metabolism was required for induction with phenanthrene. The magnitude of GFP induction was influenced by physical parameters affecting the phenanthrene flux to the cells, such as the contact surface area between solid phenanthrene and the aqueous phase, addition of surfactant, and slow phenanthrene release from Model Polymer Release System beads or from a water-immiscible oil. These results strongly suggest that the bacterial biosensor can sense different phenanthrene fluxes while maintaining phenanthrene metabolism, thus acting as a genuine sensor for phenanthrene bioavailability. A relationship between GFP production and phenanthrene mass transfer is proposed.

On the diurnal soil water content dynamics during evaporation using dielectric methods

The water content dynamics in the upper soil surface during evaporation is a key element in land-atmosphere exchanges. Previous experimental studies have suggested that the soil water content increases at the depth of 5 to 15 cm below the soil surface during evapo- ration, while the layer in the immediate vicinity of the soil surface is drying. In this study, the dynamics of water content profiles exposed to solar radiative forcing was monitored at a high temporal resolution using dielectric methods both in the presence and absence of evaporation. A 4-d comparison of reported moisture content in coarse sand in covered and uncovered buckets using a commercial dielectric-based probe (70 MHz ECH2O-5TE, Decagon Devices, Pullman, WA) and the standard 1-GHz time domain reflectometry method. Both sensors reported a positive correlation between temperature and water content in the 5- to 10-cm depth, most pronounced in the morning during heating and in the afternoon during cooling. Such positive correlation might have a physical origin induced by evaporation at the surface and redistribution due to liquid water fluxes resulting from the temperature- gradient dynamics within the sand profile at those depths. Our experimental data suggest that the combined effect of surface evaporation and temperature-gradient dynamics should be considered to analyze experimental soil water profiles. Additional effects related to the frequency of operation and to protocols for temperature compensation of the dielectric sensors may also affect the probes' response during large temperature changes.

Controls on 13C and 14C variability in soil CO2

14C dating of groundwater depends on the isotopic composition of both the solid carbonate and the soil CO2 and requires the use of 14C age correction models. To better assess the variability of the 14C activity of soil CO2 (A14Csoil-CO2) and the δ13C of soil CO2 (δ13Csoil-CO2), which are two parameters used in 14C age correction models, we studied the different processes involving carbon isotopes in the soil. The approach used experimental data from two sites in France (Fontainebleau sands and Astian sands) and a steady-state transport model. In most cases, the 14C activity (A14C) of atmospheric CO2 is directly used in the 14C age correction models as the A14Csoil-CO2. However, we demonstrate that since 1950, the evolution of the A14Csoil-CO2 reflects the competition between the fluxes of root-derived CO2 and organic matter-derived CO2. Therefore, the A14Csoil-CO2 must be used to date groundwater that is younger than 60 years old. Moreover, the δ13C of soil CO2 (δ13Csoil-CO2) showed large seasonal variations that must be taken into account in selecting the δ13Csoil-CO2 for 14C age correction models.

Getting the most sulfate from soil: Regulation of sulfate uptake transporters in Arabidopsis.

Sulfur (S) is an essential macronutrient for all living organisms. Plants require large amounts of sulfate for growth and development, and this serves as a major entry point of sulfate into the food web. Plants acquire S in its ionic form from the soil; they have evolved tightly controlled mechanisms for the regulation of sulfate uptake in response to its external and internal availability. In the model plant Arabidopsis thaliana, the first key step in sulfate uptake is presumed to be carried out exclusively by only two high-affinity sulfate transporters: SULTR1;1 and SULTR1;2. A better understanding of the mode of regulation for these two transporters is crucial because they constitute the first determinative step in balancing sulfate in respect to its supply and demand. Here, we review the recent progress achieved in our comprehension of (i) mechanisms that regulate these two high-affinity sulfate transporters at the transcriptional and post-transcriptional levels, and (ii) their structure-function relationship. Such progress is important to enable biotechnological and agronomic strategies aimed at enhancing sulfate uptake and improving crop yield in S-deficient soils.

Fungi, bacteria and soil pH: the oxalate-carbonate pathway as a model for metabolic interaction

The oxalatecarbonate pathway involves the oxidation of calcium oxalate to low-magnesium calcite and represents a potential long-term terrestrial sink for atmospheric CO2. In this pathway, bacterial oxalate degradation is associated with a strong local alkalinization and subsequent carbonate precipitation. In order to test whether this process occurs in soil, the role of bacteria, fungi and calcium oxalate amendments was studied using microcosms. In a model system with sterile soil amended with laboratory cultures of oxalotrophic bacteria and fungi, the addition of calcium oxalate induced a distinct pH shift and led to the final precipitation of calcite. However, the simultaneous presence of bacteria and fungi was essential to drive this pH shift. Growth of both oxalotrophic bacteria and fungi was confirmed by qPCR on the frc (oxalotrophic bacteria) and 16S rRNA genes, and the quantification of ergosterol (active fungal biomass) respectively. The experiment was replicated in microcosms with non-sterilized soil. In this case, the bacterial and fungal contribution to oxalate degradation was evaluated by treatments with specific biocides (cycloheximide and bronopol). Results showed that the autochthonous microflora oxidized calcium oxalate and induced a significant soil alkalinization. Moreover, data confirmed the results from the model soil showing that bacteria are essentially responsible for the pH shift, but require the presence of fungi for their oxalotrophic activity. The combined results highlight that the interaction between bacteria and fungi is essential to drive metabolic processes in complex environments such as soil.

Indirect hydrogen analysis by gas chromatography coupled to mass spectrometry (GC-MS).

Gas chromatography (GC) is an analytical tool very useful to investigate the composition of gaseous mixtures. The different gases are separated by specific columns but, if hydrogen (H2 ) is present in the sample, its detection can be performed by a thermal conductivity detector or a helium ionization detector. Indeed, coupled to GC, no other detector can perform this detection except the expensive atomic emission detector. Based on the detection and analysis of H2 isotopes by low-pressure chemical ionization mass spectrometry (MS), a new method for H2 detection by GC coupled to MS with an electron ionization ion source and a quadrupole analyser is presented. The presence of H2 in a gaseous mixture could easily be put in evidence by the monitoring of the molecular ion of the protonated carrier gas. Copyright © 2013 John Wiley & Sons, Ltd.

Validation of an analytical method for nitrous oxide (N2O) laughing gas by headspace gas chromatography coupled to mass spectrometry (HS-GC-MS): Forensic application to a lethal intoxication.

Drug abuse is a widespread problem affecting both teenagers and adults. Nitrous oxide is becoming increasingly popular as an inhalation drug, causing harmful neurological and hematological effects. Some gas chromatography-mass spectrometry (GC-MS) methods for nitrous oxide measurement have been previously described. The main drawbacks of these methods include a lack of sensitivity for forensic applications; including an inability to quantitatively determine the concentration of gas present. The following study provides a validated method using HS-GC-MS which incorporates hydrogen sulfide as a suitable internal standard allowing the quantification of nitrous oxide. Upon analysis, sample and internal standard have similar retention times and are eluted quickly from the molecular sieve 5Å PLOT capillary column and the Porabond Q column therefore providing rapid data collection whilst preserving well defined peaks. After validation, the method has been applied to a real case of N2O intoxication indicating concentrations in a mono-intoxication.

Rapid high-performance liquid chromatographic determination with fluorescence detection of furosemide in human body fluids and its confirmation by gas chromatography-mass spectrometry.

Furosemide (FD: Lasix) is a loop diuretic which strongly increases both urine flow and electrolyte urinary excretion. Healthy volunteers were administered 40 mg orally (dissolved in water) and concentrations of FD were determined in serum and urine for up to 6 h for eight subjects, who absorbed water at a rate of 400 ml/h. Quantification was performed by HPLC with fluorescence detection (excitation at 233 nm, emission at 389 nm) with a limit of detection of 5 ng/ml for a 300-microliters sample. The elution of FD was completed within 4 min using a gradient of acetonitrile concentration rising from 30 to 50% in 0.08 M phosphoric acid. The delay to the peak serum concentration ranged from 60 to 120 min. FD was still easily measurable in the sera from all subjects 6 h after administration. In urine, the excretion rates reached their maximum between 1 and 3 h. The total amount of FD excreted in the urine averaged 11.2 mg (range 7.6-14.0 mg), with a mean urine volume of 3024 ml (range 2620-3596 ml). Moreover, the urine density was lower than 1.010 (recommended as an upper limit in doping analysis to screen diuretics) only for 2 h. An additional volunteer was administered 40 mg of FD and his urine was collected over a longer period. FD was still detectable 48 h after intake. Gas chromatography-mass spectrometry with different types of ionization was used to confirm the occurrence of FD after permethylation of the extract. Negative-ion chemical ionization, with ammonia as reactant gas, was found to be the most sensitive method of detection.