Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Effects of small surface tension in Hele-Shaw multifinger dynamics: an analytical and numerical study

We study the singular effects of vanishingly small surface tension on the dynamics of finger competition in the Saffman-Taylor problem, using the asymptotic techniques described by Tanveer [Philos. Trans. R. Soc. London, Ser. A 343, 155 (1993)] and Siegel and Tanveer [Phys. Rev. Lett. 76, 419 (1996)], as well as direct numerical computation, following the numerical scheme of Hou, Lowengrub, and Shelley [J. Comput. Phys. 114, 312 (1994)]. We demonstrate the dramatic effects of small surface tension on the late time evolution of two-finger configurations with respect to exact (nonsingular) zero-surface-tension solutions. The effect is present even when the relevant zero-surface-tension solution has asymptotic behavior consistent with selection theory. Such singular effects, therefore, cannot be traced back to steady state selection theory, and imply a drastic global change in the structure of phase-space flow. They can be interpreted in the framework of a recently introduced dynamical solvability scenario according to which surface tension unfolds the structurally unstable flow, restoring the hyperbolicity of multifinger fixed points.

Topographic forcing of flow partition and flow structures at river bifurcations

This paper presents the predicted flow dynamics from the application of a Reynolds-averaged NavierStokes model to a series of bifurcation geometries with morphologies measured during previous flume experiments. The topography of the bifurcations consists of either plane or bedform-dominated beds which may or may not possess discordance between the two bifurcation distributaries. Numerical predictions are compared with experimental results to assess the ability of the numerical model to reproduce the division of flow into the bifurcation distributaries. The hydrodynamic model predicts: (1) diverting fluxes in the upstream channel which direct water into the distributaries; (2) super-elevation of the free surface induced at the bifurcation edge by pressure differences; and (3) counter-rotating secondary circulation cells which develop upstream of the apex of the bifurcation and move into the downstream channels, with water converging at the surface and diverging at the bed. When bedforms are not present, weak transversal fluxes characterize the upstream channel for almost its entire length, associated with clearly distinguishable secondary circulation cells, although these may be under-estimated by the turbulence model used in the solution. In the bedform dominated case, the same hydrodynamic conditions were not observed, with the bifurcation influence restricted and depth scale secondary circulation cells not forming. The results also demonstrate the dominant effect bed discordance has upon flow division between the two distributaries. Finally, results indicate that in bedform dominated rivers. Consequently, we suggest that sand-bed river bifurcations are more likely to have an influence that extends much further upstream and have a greater impact upon water distribution. This may contribute to observed morphological differences between sand-bedded and gravel-bedded braided river networks. Copyright (C) 2012 John Wiley & Sons, Ltd.

Properties of debris flow deposits and source materials compared: implications for debris flow characterization

In the Alps, debris flow deposits generally contain < 5% clay-size particles, and the role of the surface-charged < 2 mu m particles is often neglected, although these particles may have a significant impact on the rheological properties of the interstitial fluid. The objective of this study was to compare debris flow deposits and parent materials from two neighbouring catchments of the Swiss Alps, with special emphasis on the colloidal constituents. The catchments are small in area (4 km(2)), 2.5 km long, similar in morphology, but different in geology. The average slopes are 35-40%. The catchments were monitored for debris flow events and mapped for surface aspect and erosion activity. Debris flow deposits and parent materials were sampled, the clay and silt fractions extracted and the bulk density, < 2 mm fraction bulk density, particle size distribution, chemical composition, cation exchange capacity (CEC) and mineralogy analysed. The results show that the deposits are similar to the parent screes in terms of chemical composition, but differ in terms of: (i) particle size distribution; and (ii) mineralogy, reactivity and density of the < 2 mm fraction. In this fraction, compared with the parent materials the deposits show dense materials enriched in coarse monocrystalline particles, of which the smallest and more reactive particles were leached. The results suggest that deposit samples should not be considered as representative of source or flow materials, particularly with respect to their physical properties.

Quality of surface water related to land use: a case study in a catchment with small farms and intensive vegetable crop production in southern Brazil

Water degradation is strongly related to agricultural activity. The aim of this study was to evaluate the influence of land use and some environmental components on surface water quality in the Campestre catchment, located in Colombo, state of Parana, Brazil. Physical and chemical attributes were analyzed (total nitrogen, ammonium, nitrate, total phosphorus, electrical conductivity, pH, temperature, turbidity, total solids, biological oxygen demand, chemical oxygen demand and dissolved oxygen). Monthly samples of the river water were taken over one year at eight monitoring sites, distributed over three sub-basins. Overall, water quality was worse in the sub-basin with a higher percentage of agriculture, and was also affected by a lower percentage of native forest and permanent preservation area, and a larger drainage area. Water quality was also negatively affected by the presence of agriculture in the riparian zone. In the summer season, probably due to higher rainfall and intensive soil use, a higher concentration of total nitrogen and particulate nitrogen was observed, as well as higher electrical conductivity, pH and turbidity. All attributes, except for total phosphorus, were in compliance with Brazilian Conama Resolution Nº 357/2005 for freshwater class 1. However, it should be noted that these results referred to the base flow and did not represent a discharge condition since most of the water samples were not collected at or near the rainfall event.

Flow structures at an idealized bifurcation : a numerical experiment

River bifurcations are key nodes within braided river systems controlling the flow and sediment partitioning and therefore the dynamics of the river braiding process. Recent research has shown that certain geometrical configurations induce instabilities that lead to downstream mid-channel bar formation and the formation of bifurcations. However, we currently have a poor understanding of the flow division process within bifurcations and the flow dynamics in the downstream bifurcates, both of which are needed to understand bifurcation stability. This paper presents results of a numerical sensitivity experiment undertaken using computational fluid dynamics (CFD) with the purpose of understanding the flow dynamics of a series of idealized bifurcations. A geometric sensitivity analysis is undertaken for a range of channel slopes (0.005 to 0.03), bifurcation angles (22 degrees to 42 degrees) and a restricted set of inflow conditions based upon simulating flow through meander bends with different curvature on the flow field dynamics through the bifurcation. The results demonstrate that the overall slope of the bifurcation affects the velocity of flow through the bifurcation and when slope asymmetry is introduced, the flow structures in the bifurcation are modified. In terms of bifurcation evolution the most important observation appears to be that once slope asymmetry is greater than 0.2 the flow within the steep bifurcate shows potential instability and the potential for alternate channel bar formation. Bifurcation angle also defines the flow structures within the bifurcation with an increase in bifurcation angle increasing the flow velocity down both bifurcates. However, redistributive effects of secondary circulation caused by upstream curvature can very easily counter the effects of local bifurcation characteristics. Copyright (C) 2011 John Wiley & Sons, Ltd.

Detection and Characterization of Preferential Flow Paths in the Downstream Area of a Hazardous Landfill

We have used surface-based electrical resistivity tomography to detect and characterize preferential hydraulic pathways in the immediate downstream area of an abandoned, hazardous landfill. The landfill occupies the void left by a former gravel pit and its base is close to the groundwater table and lacking an engineered barrier. As such, this site is remarkably typical of many small- to medium-sized waste deposits throughout the densely populated and heavily industrialized foreland on both sides of the Alpine arc. Outflows of pollutants lastingly contaminated local drinking water supplies and necessitated a partial remediation in the form of a synthetic cover barrier, which is meant to prevent meteoric water from percolating through the waste before reaching the groundwater table. Any future additional isolation of the landfill in the form of lateral barriers thus requires adequate knowledge of potential preferential hydraulic pathways for outflowing contaminants. Our results, inferred from a suite of tomographically inverted surfaced-based electrical resistivity profiles oriented roughly perpendicular to the local hydraulic gradient, indicate that potential contaminant outflows would predominantly occur along an unexploited lateral extension of the original gravel deposit. This finds its expression as a distinct and laterally continuous high-resistivity anomaly in the resistivity tomograms. This interpretation is ground-truthed through a litholog from a nearby well. Since the probed glacio-fluvial deposits are largely devoid of mineralogical clay, the geometry of hydraulic and electrical pathways across the pore space of a given lithological unit can be assumed to be identical, which allows for an order-of-magnitude estimation of the overall permeability structure. These estimates indicate that the permeability of the imaged extension of the gravel body is at least two to three orders-of-magnitude higher than that of its finer-grained embedding matrix. This corroborates the preeminent role of the high-resistivity anomaly as a potential preferential flow path.

Characterization of surface functional groups present on laboratory-generated and ambient aerosol particles by means of heterogeneous titration reactions

A Knudsen flow reactor has been used to quantify surface functional groups on aerosols collected in the field. This technique is based on a heterogeneous titration reaction between a probe gas and a specific functional group on the particle surface. In the first part of this work, the reactivity of different probe gases on laboratory-generated aerosols (limonene SOA, Pb(NO3)2, Cd(NO3)2) and diesel reference soot (SRM 2975) has been studied. Five probe gases have been selected for the quantitative determination of important functional groups: N(CH3)3 (for the titration of acidic sites), NH2OH (for carbonyl functions), CF3COOH and HCl (for basic sites of different strength), and O3 (for oxidizable groups). The second part describes a field campaign that has been undertaken in several bus depots in Switzerland, where ambient fine and ultrafine particles were collected on suitable filters and quantitatively investigated using the Knudsen flow reactor. Results point to important differences in the surface reactivity of ambient particles, depending on the sampling site and season. The particle surface appears to be multi-functional, with the simultaneous presence of antagonistic functional groups which do not undergo internal chemical reactions, such as acid-base neutralization. Results also indicate that the surface of ambient particles was characterized by a high density of carbonyl functions (reactivity towards NH2OH probe in the range 0.26-6 formal molecular monolayers) and a low density of acidic sites (reactivity towards N(CH3)3 probe in the range 0.01-0.20 formal molecular monolayer). Kinetic parameters point to fast redox reactions (uptake coefficient ?0>10-3 for O3 probe) and slow acid-base reactions (?0<10-4 for N(CH3)3 probe) on the particle surface. [Authors]

Foot skin blood flow following infrainguinal revascularization for critical lower limb ischemia.

INTRODUCTION: The aim of this study was to assess the blood flow in the feet before and after lower limb revascularization using laser Doppler imaging (LDI). METHODS: Ten patients with critical lower limb ischemia were prospectively enrolled from June to October 2004. All patients underwent successful unilateral surgical interventions including above-knee bypass, distal bypass and endarterectomy. Skin blood flow (SBF) over the plantar surface of both forefeet and heels was measured by LDI 24h before and 10 days after revascularization, expressed in perfusion units (PU), and reported as mean+/-SD. RESULTS: Measurements in the forefoot and heel were similar. Before revascularization mean SBF was significantly lower in the ischemic foot (130+/-71 PU) compared to the contralateral foot (212+/-68 PU), p<0.05. After revascularization a significant increase of the SBF in the forefoot (from 135+/-67 to 202+/-86 PU, p=0.001) and hindfoot (from 148+/-58 to 203+/-83, p=0.001) was observed on the treatment side. However, a large decrease of the SBF was seen in forefoot and hindfoot on the untreated side (from 250+/-123 PU to 176+/-83 and from 208+/-116 to 133+/-40, p=0.001, respectively). CONCLUSION: This study confirms the benefits of revascularization in patients with nonhealing foot lesions due to critical limb ischemia. A significant increase of the SBF was observed on the treatment side. However, an unexpected decrease was observed on the untreated side.

Development of Elispot and CFSE flow cytometric assays for detecting beryllium sensitivity

Background: Beryllium sensitization (BeS) is caused by exposure to beryllium in the workplace and may progress to chronic beryllium disease (CBD). This granulomatous lung disorder mimicks sarcoidosis clinically, but is characterized by beryllium specific CD4+ T-cells immune response. BeS is classically detected by beryllium lymphocyte proliferation test (BeLPT), but this assay requires radioactivity and is not very sensitive. In the context of a study aiming to evaluate if CBD patients are misdiagnosed as sarcoidosis patients in Switzerland, we developed EliSpot and CFSE beryllium flow cytometric test. Methods: 23 patients considered as having sarcoidosis (n = 21), CBD (n = 1) and possible CBD (n = 1) were enrolled. Elispot was performed using plate covered with gamma-IFN mAb. Cells were added to wells and incubated overnight at 37 °C with medium (neg ctrl), SEB (pos ctrl) or BeSO4 at 1, 10 and 100 microM. Anti-IFN-gamma biotinylated mAb were added and spots were visualized using streptavidinhorseradish peroxidase and AEC substrate reagent. Results were reported as spot forming unit (SFU). For Beryllium specific CFSE flow cytometry analysis, CFSE labelled cells were cultured in the presence of SEB and 1, 10 or 100 microM BeSO4. Unstimulated CFSE labeled cells were defined as controls. The cells were incubated for 6 days at 37 °C and 5% CO2. Surface labelling of T-lymphocytes and vivid as control of cells viability was performed at the time of harvest. Results: Using EliSpot technology, we were able to detect a BeS in 1/23 enrolled patients with a mean of 780 SFU (cut off value at 50 SFU). This positive result was confirmed using different concentration of BeSO4. Among the 23 patients tested, 22 showed negative results with EliSpot. Using CFSE flow cytometry, 1/7 tested patients showed a positive result with a beryllium specific CD4+ count around 30% versus 45% for SEB stimulation as positif control and 0.6 % for negative control. This patient was the one with a positive EliSpot assay. Conclusions: The preliminary data demonstrated the feasibility of Elispot and CFSE flow cytometry to detect BeS. The patient with a beryllium specific positive EliSpot and CFSE flow cytometry result had been exposed to beryllium at her workplace 20 years ago and is still regularly controlled for her pulmonary status. A positive BeLPT had already been described in 2001 in France for this patient. Further validation of these techniques are in progress.

LiveCount Assay: concomitant measurement of cytolytic activity and phenotypic characterisation of CD8(+) T-cells by flow cytometry.

Tumour immunologists strive to develop efficient tumour vaccination and adoptive transfer therapies that enlarge the pool of tumour-specific and -reactive effector T-cells in vivo. To assess the efficiency of the various strategies, ex vivo assays are needed for the longitudinal monitoring of the patient's specific immune responses providing both quantitative and qualitative data. In particular, since tumour cell cytolysis is the end goal of tumour immunotherapy, routine immune monitoring protocols need to include a read-out for the cytolytic efficiency of Ag-specific cells. We propose to combine current immune monitoring techniques in a highly sensitive and reproducible multi-parametric flow cytometry based cytotoxicity assay that has been optimised to require low numbers of Ag-specific T-cells. The possibility of re-analysing those T-cells that have undergone lytic activity is illustrated by the concomitant detection of CD107a upregulation on the surface of degranulated T-cells. To date, the LiveCount Assay provides the only possibility of assessing the ex vivo cytolytic activity of low-frequency Ag-specific cytotoxic T-lymphocytes from patient material.

Is it appropriate to index glomerular filtration rate for body surface area

Purpose: Obesity is an established independent risk factor for chronic kidney disease. Thus, measurement of glomerular filtration rate (GFR) is important in this population. Traditionally, GFR has been indexed for body surface area (BSA), but this indexation may not be appropriate in obese individuals. Therefore, the objective of the study was to compare absolute GFR with GFR indexed for BSA and with GFR indexed for height. Methods and materials: The study was conducted in 66 families from the Seychelles islands that included several members with hypertension. GFR and effective renal plasma flow (ERPF) were measured using inulin and PAH clearances, respectively. Antihypertensive treatment, if used, was withheld 2 weeks before conducting the clearances. Participants with diabetes mellitus were excluded from the analysis. BSA was calculated using the Dubois formula. We assessed trend across BMI categories using a non parametric test. Results: Participants included 174 women and 127 men. The prevalence of hypertension was 61%, of which 68% were treated. The table shows that absolute GFR, GFR indexed for height, ERPF, filtration fraction were significantly higher across BMI categories. When GFR was indexed for BSA, the association between GFR and BMI categories was lost. Conclusion: Indexing GFR for BSA in overweight and obese individuals leads to a substantial underestimation of GFR. Filtration fraction, which does not depend on BSA, is higher in obese individuals, which suggests glomerular hyperfiltration. Indexing GFR for BSA therefore would mask the underlying glomerular hyperfiltration. As the number of nephrons does not increase with weight gain, absolute GFR represents a better marker of single nephron GFR and is more appropriate.

A framework for hybrid simulations of two-phase flow in porous media

Rough a global coarse problem. Although these techniques are usually employed for problems in which the fine-scale processes are described by Darcy's law, they can also be applied to pore-scale simulations and used as a mathematical framework for hybrid methods that couples a Darcy and pore scales. In this work, we consider a pore-scale description of fine-scale processes. The Navier-Stokes equations are numerically solved in the pore geometry to compute the velocity field and obtain generalized permeabilities. In the case of two-phase flow, the dynamics of the phase interface is described by the volume of fluid method with the continuum surface force model. The MsFV method is employed to construct an algorithm that couples a Darcy macro-scale description with a pore-scale description at the fine scale. The hybrid simulations results presented are in good agreement with the fine-scale reference solutions. As the reconstruction of the fine-scale details can be done adaptively, the presented method offers a flexible framework for hybrid modeling.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work 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 model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Fate of the herbicide glyphosate and its metabolite AMPA in soils and their transfer to surface waters: A multi-scale approach in the Lavaux vineyards, western Switzerland

The use of herbicides in agriculture may lead to environmental problems, such as surface water pollution, with a potential risk for aquatic organisms. The herbicide glyphosate is the most used active ingredient in the world and in Switzerland. In the Lavaux vineyards it is nearly the only molecule applied. This work aimed at studying its fate in soils and its transfer to surface waters, using a multi-scale approach: from molecular (10-9 m) and microscopic scales (10-6 m), to macroscopic (m) and landscape ones (103 m). First of all, an analytical method was developed for the trace level quantification of this widely used herbicide and its main by-product, aminomethylphosphonic acid (AMPA). Due to their polar nature, their derivatization with 9-fluorenylmethyl chloroformate (FMOC-Cl) was done prior to their concentration and purification by solid phase extraction. They were then analyzed by ultra performance liquid chromatography coupled with tandem mass spectrometry (UPLC-MS/MS). The method was tested in different aqueous matrices with spiking tests and validated for the matrix effect correction in relevant environmental samples. Calibration curves established between 10 and 1000ng/l showed r2 values above 0.989, mean recoveries varied between 86 and 133% and limits of detection and quantification of the method were as low as 5 and 10ng/l respectively. At the parcel scale, two parcels of the Lavaux vineyard area, located near the Lutrive River at 6km to the east of Lausanne, were monitored to assess to which extent glyphosate and AMPA were retained in the soil or exported to surface waters. They were equipped at their bottom with porous ceramic cups and runoff collectors, which allowed retrieving water samples for the growing seasons 2010 and 2011. Results revealed that the mobility of glyphosate and AMPA in the unsaturated zone was likely driven by the precipitation regime and the soil characteristics, such as slope, porosity structure and layer permeability discrepancy. Elevated glyphosate and AMPA concentrations were measured at 60 and 80 cm depth at parcel bottoms, suggesting their infiltration in the upper parts of the parcels and the presence of preferential flow in the studied parcels. Indeed, the succession of rainy days induced the gradual saturation of the soil porosity, leading to rapid infiltration through macropores, as well as surface runoff formation. Furthermore, the presence of more impervious weathered marls at 100 cm depth induced throughflows, the importance of which for the lateral transport of the herbicide molecules was determined by the slope steepness. Important rainfall events (>10 mm/day) were clearly exporting molecules from the soil top layer, as indicated by important concentrations in runoff samples. A mass balance showed that total loss (10-20%) mainly occurred through surface runoff (96%) and, to a minor extent, by throughflows in soils (4%), with subsequent exfiltration to surface waters. Observations made in the Lutrive River revealed interesting details of glyphosate and AMPA dynamics in urbanized landscapes, such as the Lavaux vineyards. Indeed, besides their physical and chemical properties, herbicide dynamics at the catchment level strongly depend on application rates, precipitation regime, land use and also on the presence of drains or constructed channels. Elevated concentrations, up to 4970 ng/l, observed just after the application, confirmed the diffuse export of these compounds from the vineyard area by surface runoff during main rain events. From April to September 2011, a total load of 7.1 kg was calculated, with 85% coming from vineyards and minor urban sources and 15% from arable crops. Small vineyard surfaces could generate high concentrations of herbicides and contribute considerably to the total load calculated at the outlet, due to their steep slopes (~10%). The extrapolated total amount transferred yearly from the Lavaux vineyards to the Lake of Geneva was of 190kg. At the molecular scale, the possible involvement of dissolved organic matter (DOM) in glyphosate and copper transport was studied using UV/Vis fluorescence spectroscopy. Combined with parallel factor (PARAFAC) analysis, this technique allowed characterizing DOM of soil and surface water samples from the studied vineyard area. Glyphosate concentrations were linked to the fulvic-like spectroscopic signature of DOM in soil water samples, as well as to copper, suggesting the formation of ternary complexes. In surface water samples, its concentrations were also correlated to copper ones, but not in a significant way to the fulvic-like signature. Quenching experiments with standards confirmed field tendencies in the laboratory, with a stronger decrease in fluorescence intensity for fulvic-like fluorophore than for more aromatic ones. Lastly, based on maximum concentrations measured in the river, an environmental risk for these compounds was assessed, using laboratory tests and ecotoxicity data from the literature. In our case and with the methodology applied, the risk towards aquatic species was found negligible (RF<1).