3D geoelectrical monitoring to follow hydric soil behaviour during rainfall

RESUME Durant les dernières années, les méthodes électriques ont souvent été utilisées pour l'investigation des structures de subsurface. L'imagerie électrique (Electrical Resistivity Tomography, ERT) est une technique de prospection non-invasive et spatialement intégrée. La méthode ERT a subi des améliorations significatives avec le développement de nouveaux algorithmes d'inversion et le perfectionnement des techniques d'acquisition. La technologie multicanale et les ordinateurs de dernière génération permettent la collecte et le traitement de données en quelques heures. Les domaines d'application sont nombreux et divers: géologie et hydrogéologie, génie civil et géotechnique, archéologie et études environnementales. En particulier, les méthodes électriques sont souvent employées dans l'étude hydrologique de la zone vadose. Le but de ce travail est le développement d'un système de monitorage 3D automatique, non- invasif, fiable, peu coûteux, basé sur une technique multicanale et approprié pour suivre les variations de résistivité électrique dans le sous-sol lors d'événements pluvieux. En raison des limitations techniques et afin d'éviter toute perturbation physique dans la subsurface, ce dispositif de mesure emploie une installation non-conventionnelle, où toutes les électrodes de courant sont placées au bord de la zone d'étude. Le dispositif le plus approprié pour suivre les variations verticales et latérales de la résistivité électrique à partir d'une installation permanente a été choisi à l'aide de modélisations numériques. Les résultats démontrent que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle et plus apte à détecter les variations latérales et verticales de la résistivité électrique, et cela malgré la configuration non-conventionnelle des électrodes. Pour tester l'efficacité du système proposé, des données de terrain ont été collectées sur un site d'étude expérimental. La technique de monitorage utilisée permet de suivre le processus d'infiltration 3D pendant des événements pluvieux. Une bonne corrélation est observée entre les résultats de modélisation numérique et les données de terrain, confirmant par ailleurs que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle. La nouvelle technique de monitorage 3D de résistivité électrique permet de caractériser les zones d'écoulement préférentiel et de caractériser le rôle de la lithologie et de la pédologie de manière quantitative dans les processus hydrologiques responsables d'écoulement de crue. ABSTRACT During the last years, electrical methods were often used for the investigation of subsurface structures. Electrical resistivity tomography (ERT) has been reported to be a useful non-invasive and spatially integrative prospecting technique. The ERT method provides significant improvements, with the developments of new inversion algorithms, and the increasing efficiency of data collection techniques. Multichannel technology and powerful computers allow collecting and processing resistivity data within few hours. Application domains are numerous and varied: geology and hydrogeology, civil engineering and geotechnics, archaeology and environmental studies. In particular, electrical methods are commonly used in hydrological studies of the vadose zone. The aim of this study was to develop a multichannel, automatic, non-invasive, reliable and inexpensive 3D monitoring system designed to follow electrical resistivity variations in soil during rainfall. Because of technical limitations and in order to not disturb the subsurface, the proposed measurement device uses a non-conventional electrode set-up, where all the current electrodes are located near the edges of the survey grid. Using numerical modelling, the most appropriate arrays were selected to detect vertical and lateral variations of the electrical resistivity in the framework of a permanent surveying installation system. The results show that a pole-dipole array has a better resolution than a pole-pole array and can successfully follow vertical and lateral resistivity variations despite the non-conventional electrode configuration used. Field data are then collected at a test site to assess the efficiency of the proposed monitoring technique. The system allows following the 3D infiltration processes during a rainfall event. A good correlation between the results of numerical modelling and field data results can be observed since the field pole-dipole data give a better resolution image than the pole-pole data. The new device and technique makes it possible to better characterize the zones of preferential flow and to quantify the role of lithology and pedology in flood- generating hydrological processes.

Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers

Recent technological advances in remote sensing have enabled investigation of the morphodynamics and hydrodynamics of large rivers. However, measuring topography and flow in these very large rivers is time consuming and thus often constrains the spatial resolution and reach-length scales that can be monitored. Similar constraints exist for computational fluid dynamics (CFD) studies of large rivers, requiring maximization of mesh-or grid-cell dimensions and implying a reduction in the representation of bedform-roughness elements that are of the order of a model grid cell or less, even if they are represented in available topographic data. These ``subgrid'' elements must be parameterized, and this paper applies and considers the impact of roughness-length treatments that include the effect of bed roughness due to ``unmeasured'' topography. CFD predictions were found to be sensitive to the roughness-length specification. Model optimization was based on acoustic Doppler current profiler measurements and estimates of the water surface slope for a variety of roughness lengths. This proved difficult as the metrics used to assess optimal model performance diverged due to the effects of large bedforms that are not well parameterized in roughness-length treatments. However, the general spatial flow patterns are effectively predicted by the model. Changes in roughness length were shown to have a major impact upon flow routing at the channel scale. The results also indicate an absence of secondary flow circulation cells in the reached studied, and suggest simpler two-dimensional models may have great utility in the investigation of flow within large rivers. Citation: Sandbach, S. D. et al. (2012), Application of a roughness-length representation to parameterize energy loss in 3-D numerical simulations of large rivers, Water Resour. Res., 48, W12501, doi: 10.1029/2011WR011284.

Occupational exposure of truck drivers to dust and polynuclear aromatic hydrocarbons: a pilot study in Geneva, Switzerland.

The exposure to dust and polynuclear aromatic hydrocarbons (PAH) of 15 truck drivers from Geneva, Switzerland, was measured. The drivers were divided between "long-distance" drivers and "local" drivers and between smokers and nonsmokers and were compared with a control group of 6 office workers who were also divided into smokers and nonsmokers. Dust was measured on 1 workday both by a direct-reading instrument and by sampling. The local drivers showed higher exposure to dust (0.3 mg/m3) and PAH than the long-distance drivers (0.1 mg/m3), who showed no difference with the control group. This observation may be due to the fact that the local drivers spend more time in more polluted areas, such as streets with heavy traffic and construction sites, than do the long-distance drivers. Smoking does not influence exposure to dust and PAH of professional truck drivers, as measured in this study, probably because the ventilation rate of the truck cabins is relatively high even during cold days (11-15 r/h). The distribution of dust concentrations was shown in some cases to be quite different from the expected log-normal distribution. The contribution of diesel exhaust to these exposures could not be estimated since no specific tracer was used. However, the relatively low level of dust exposure dose not support the hypothesis that present day levels of diesel exhaust particulates play a significant role in the excess occurrence of lung cancer observed in professional truck drivers.

Enhanced Proinflammatory Cytokine Response to Bacterial Lipopolysaccharide in the Adult Male Rat after either Neonatal or Prepubertal Ablation of Biological Testosterone Activity.

A sex steroid-dependent modulation of the immune function in mammals is accepted, and evidence suggests that while estrogens enhance, androgens inhibit the immune response. The aim of this study was to explore in the adult male rat the effect of either neonatal flutamide (FTM) treatment or prepubertal orchidectomy (ODX) on endocrine markers in the basal condition and peripheral tumor necrosis factor alpha (TNFα) levels during inflammatory stress. For these purposes, (1) 5-day-old male rats were subcutaneously injected with either sterile vehicle alone or containing 1.75 mg FTM, and (2) 25-day-old male rats were sham operated or had ODX. Rats were sacrificed (at 100 days of age) in the basal condition for determination of peripheral metabolite levels. Additional rats were intravenously injected with bacterial lipopolysaccharide (LPS; 25 μg/kg body weight, i.v.) and bled for up to 4 h. Data indicate that (1) ODX increased peripheral glucocorticoid levels and reduced those of testosterone, whereas FTM-treated rats displayed low circulating leptin concentrations, and (2) LPS-induced TNFα secretion in plasma was significantly enhanced in the FTM and ODX groups. Our study supports that neonatal FTM treatment affected adiposity function, and adds data maintaining that androgens have a suppressive role in proinflammatory cytokine release in plasma during inflammation.

Combining probabilistic land-use change and tree population dynamics modelling to simulate responses in mountain forests

Altitudinal tree lines are mainly constrained by temperature, but can also be influenced by factors such as human activity, particularly in the European Alps, where centuries of agricultural use have affected the tree-line. Over the last decades this trend has been reversed due to changing agricultural practices and land-abandonment. We aimed to combine a statistical land-abandonment model with a forest dynamics model, to take into account the combined effects of climate and human land-use on the Alpine tree-line in Switzerland. Land-abandonment probability was expressed by a logistic regression function of degree-day sum, distance from forest edge, soil stoniness, slope, proportion of employees in the secondary and tertiary sectors, proportion of commuters and proportion of full-time farms. This was implemented in the TreeMig spatio-temporal forest model. Distance from forest edge and degree-day sum vary through feed-back from the dynamics part of TreeMig and climate change scenarios, while the other variables remain constant for each grid cell over time. The new model, TreeMig-LAb, was tested on theoretical landscapes, where the variables in the land-abandonment model were varied one by one. This confirmed the strong influence of distance from forest and slope on the abandonment probability. Degree-day sum has a more complex role, with opposite influences on land-abandonment and forest growth. TreeMig-LAb was also applied to a case study area in the Upper Engadine (Swiss Alps), along with a model where abandonment probability was a constant. Two scenarios were used: natural succession only (100% probability) and a probability of abandonment based on past transition proportions in that area (2.1% per decade). The former showed new forest growing in all but the highest-altitude locations. The latter was more realistic as to numbers of newly forested cells, but their location was random and the resulting landscape heterogeneous. Using the logistic regression model gave results consistent with observed patterns of land-abandonment: existing forests expanded and gaps closed, leading to an increasingly homogeneous landscape.

Plasma angiotensin-(1-8) octapeptide measurement to assess acute angiotensin-converting enzyme inhibition with captopril administered parenterally to normal subjects.

The blood pressure (BP), heart rate (HR), and humoral effects of single intravenous (i.v.) doses of the angiotensin-converting enzyme (ACE) inhibitor captopril was investigated in five normotensive healthy volunteers. Each subject received at 1-week intervals a bolus dose of either captopril (1, 5, and 25 mg) or its vehicle. The study was conducted in a single-blind fashion, and the order of treatment phases was randomized. The different doses of captopril had no acute effect on BP and HR. They induced a dose-dependent decrease in plasma ACE activity and plasma angiotensin II levels. The angiotensin-(1-8) octapeptide was isolated by solid-phase extraction and high-performance liquid chromatography (HPLC) prior to radioimmunoassay (RIA). All three doses of captopril reduced circulating angiotensin II levels within 15 min of drug administration. Only with the 25-mg dose was the angiotensin II concentration below the detection limit at 15 min and still significantly reduced 90 min after drug administration. Simultaneous and progressive decreases in plasma aldosterone levels were observed both with ACE inhibition and during vehicle injection, but the relative fall was more pronounced after captopril administration. No adverse reaction was noticed. These results demonstrate that captopril given parenterally blocks the renin-angiotensin system in a dose-dependent manner. Only with the dose of 25 mg was the inhibition of plasma-converting enzyme activity and the reduction of plasma angiotensin II sustained for at least 1 1/2 h.

Use of high-resolution geophysical data to characterize heterogeneous aquifers: influence of data integration method on hydrological predictions

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. The investigation of exactly how much benefit can be brought by geophysical data in terms of its effect on hydrological predictions, however, has received considerably less attention in the literature. Here, we examine the potential hydrological benefits brought by a recently introduced simulated annealing (SA) conditional stochastic simulation method designed for the assimilation of diverse hydrogeophysical data sets. We consider the specific case of integrating crosshole ground-penetrating radar (GPR) and borehole porosity log data to characterize the porosity distribution in saturated heterogeneous aquifers. In many cases, porosity is linked to hydraulic conductivity and thus to flow and transport behavior. To perform our evaluation, we first generate a number of synthetic porosity fields exhibiting varying degrees of spatial continuity and structural complexity. Next, we simulate the collection of crosshole GPR data between several boreholes in these fields, and the collection of porosity log data at the borehole locations. The inverted GPR data, together with the porosity logs, are then used to reconstruct the porosity field using the SA-based method, along with a number of other more elementary approaches. Assuming that the grid-cell-scale relationship between porosity and hydraulic conductivity is unique and known, the porosity realizations are then used in groundwater flow and contaminant transport simulations to assess the benefits and limitations of the different approaches.

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

Cooperation in centralised CATI Household Panel Surveys : A Contact-based Multilevel Analysis to Examine Interviewer, Respondent, and Fieldwork Process Effects

In this research, we analyse the contact-specific mean of the final cooperation probability, distinguishing on the one hand between contacts with household reference persons and with other eligible household members, and on the other hand between first and later contacts. Data comes from two Swiss Household Panel surveys. The interviewer-specific variance is higher for first contacts, especially in the case of the reference person. For later contacts with the reference person, the contact-specific variance dominates. This means that interaction effects and situational factors are decisive. The contact number has negative effects on the performance of contacts with the reference person, positive in the case of other persons. Also time elapsed since the previous contact has negative effects in the case of reference persons. The result of the previous contact has strong effects, especially in the case of the reference person. These findings call for a quick completion of the household grid questionnaire, assigning the best interviewers to conducting the first contact. While obtaining refusals has negative effects, obtaining other contact results has only weak effects on the interviewer's subsequent contact outcome. Using the same interviewer for contacts has no positive effects.

Response to comment on "the geometric structure of the brain fiber pathways".

In response to Catani et al., we show that corticospinal pathways adhere via sharp turns to two local grid orientations; that our studies have three times the diffusion resolution of those compared; and that the noted technical concerns, including crossing angles, do not challenge the evidence of mathematically specific geometric structure. Thus, the geometric thesis gives the best account of the available evidence.

A system to test the toxicity of brake wear particles

Fine particulate matter from traffic increases mortality and morbidity. An important source of traffic particles is brake wear. American studies reported cars to emit break wear particles at a rate of about 11mg/km to 20mg/km of driven distance. A German study estimated that break wear contributes about 12.5% to 21% of the total traffic particle emissions. The goal of this study was to build a system that allows the study of brake wear particle emissions during different braking behaviours of different car and brake types. The particles should be characterize in terms of size, number, metal, and elemental and organic carbon composition. In addition, the influence of different deceleration schemes on the particle composition and size distribution should be studied. Finally, this system should allow exposing human cell cultures to these particles. An exposure-box (0.25 cubic-m volume) was built that can be mounted around a car's braking system. This allows exposing cells to fresh brake wear particles. Concentrations of particle numbers, mass and surface, metals, and carbon compounds were quantified. Tests were conducted with A549 lung epithelial cells. Five different cars and two typical braking behaviours (full stop and normal deceleration) were tested. Particle number and size distribution was analysed for the first six minutes. In this time, two braking events occurred. Full stop produced significantly higher particle concentrations than normal deceleration (average of 23'000 vs. 10'400 #/cm3, p= 0.016). The particle number distribution was bi-modal with one peak at 60 to 100 nm (depending on the tested car and braking behaviour) and a second peak at 200 to 400 nm. Metal concentrations varied depending on the tested car type. Iron (range of 163 to 15'600 μg/m3) and Manganese (range of 0.9 to 135 μg/m3) were present in all samples, while Copper was absent in some samples (<6 to 1220 μg/m3). The overall "fleet" metal ratio was Fe:Cu:Mn = 128:14:1. Temperature and humidity varied little. A549-cells were successfully exposed in the various experimental settings and retained their viability. Culture supernatant was stored and cell culture samples were fixated to test for inflammatory response. Analysis of these samples is ongoing. The established system allowed testing brake wear particle emissions from real-world cars. The large variability of chemical composition and emitted amounts of brake wear particles between car models seems to be related to differences between brake pad compositions of different producers. Initial results suggest that the conditions inside the exposure box allow exposing human lung epithelial cells to freshly produced brake wear particles.

Segmentation and grid generation for numerical simulations of vad connections with patient-specific data

Objectives: We are interested in the numerical simulation of the anastomotic region comprised between outflow canula of LVAD and the aorta. Segmenta¬tion, geometry reconstruction and grid generation from patient-specific data remain an issue because of the variable quality of DICOM images, in particular CT-scan (e.g. metallic noise of the device, non-aortic contrast phase). We pro¬pose a general framework to overcome this problem and create suitable grids for numerical simulations.Methods: Preliminary treatment of images is performed by reducing the level window and enhancing the contrast of the greyscale image using contrast-limited adaptive histogram equalization. A gradient anisotropic diffusion filter is applied to reduce the noise. Then, watershed segmentation algorithms and mathematical morphology filters allow reconstructing the patient geometry. This is done using the InsightToolKit library (www.itk.org). Finally the Vascular Model¬ing ToolKit (www.vmtk.org) and gmsh (www.geuz.org/gmsh) are used to create the meshes for the fluid (blood) and structure (arterial wall, outflow canula) and to a priori identify the boundary layers. The method is tested on five different patients with left ventricular assistance and who underwent a CT-scan exam.Results: This method produced good results in four patients. The anastomosis area is recovered and the generated grids are suitable for numerical simulations. In one patient the method failed to produce a good segmentation because of the small dimension of the aortic arch with respect to the image resolution.Conclusions: The described framework allows the use of data that could not be otherwise segmented by standard automatic segmentation tools. In particular the computational grids that have been generated are suitable for simulations that take into account fluid-structure interactions. Finally the presented method features a good reproducibility and fast application.