Fluid flow processes at basin scale

S u b s u r face fluid flow plays a significant role in many geologic processes and is increasingly being studied in the scale of sedimentary basins and geologic time perspective. Many economic resources such as petroleum and mineral deposits are products of basin scale fluid flow operating over large periods of time. Such ancient flow systems can be studied through analysis of diagenetic alterations and fluid inclusions to constrain physical and chemical conditions of fluids and rocks during their paleohy d r og e o l ogic evolution. Basin simulation models are useful to complement the paleohy d r og e o l ogic record preserved in the rocks and to derive conceptual models on hydraulic basin evolution and generation of economic resources. Different types of fluid flow regimes may evo l ve during basin evolution. The most important with respect to flow rates and capacity for transport of solutes and thermal energy is gr avitational fluid flow driven by the topographic configuration of a basin. Such flow systems require the basin to be elevated above sea level. Consolidational fluid flow is the principal fluid migration process in basins below sea level, caused by loading of compressible rocks. Flow rates of such systems are several orders of magnitude below topogr a p hy driven flow. Howeve r, consolidation may create significant fluid ove rpressure. Episodic dewatering of ove rpressured compart m e n t s m ay cause sudden fluid release with elevated flow velocities and may cause a transient local thermal and chemical disequilibrium betwe e n fluid and rock. This paper gives an ove rv i ew on subsurface fluid flow processes at basin scale and presents examples related to the Pe n e d è s basin in the central Catalan continental margin including the offshore Barcelona half-graben and the compressive South-Pyrenean basin.

Self-expandable metal stents in the treatment of acute esophageal variceal bleeding.

Acute variceal bleeding (AVB) is a life-threatening complication in patients with cirrhosis. Hemostatic therapy of AVB includes early administration of vasoactive drugs that should be combined with endoscopic therapy, preferably banding ligation. However, failure to control bleeding or early rebleed within 5 days still occurs in 15-20% of patients with AVB. In these cases, a second endoscopic therapy may be attempted (mild bleeding in a hemodynamically stable patient) or we can use a balloon tamponade as a bridge to definitive derivative treatment (i.e., a transjugular intrahepatic portosystemic shunt). Esophageal balloon tamponade provides initial control in up to 80% of AVB, but it carries a high risk of major complications, especially in cases of long duration of tamponade (>24 h) and when tubes are inserted by inexperienced staff. Preliminary reports suggest that self-expandable covered esophageal metallic stents effectively control refractory AVB (i.e., ongoing bleeding despite pharmacological and endoscopic therapy or massive bleeding precluding endoscopic therapy) with a low incidence of complications. Thus, covered self-expanding metal stents may represent an alternative to the Sengstaken-Blakemore balloon for the temporary control of bleeding in treatment failures. Further studies are required to determine the role of this new device in AVB.

Error models in hydrogeology applications

Notre consommation en eau souterraine, en particulier comme eau potable ou pour l'irrigation, a considérablement augmenté au cours des années. De nombreux problèmes font alors leur apparition, allant de la prospection de nouvelles ressources à la remédiation des aquifères pollués. Indépendamment du problème hydrogéologique considéré, le principal défi reste la caractérisation des propriétés du sous-sol. Une approche stochastique est alors nécessaire afin de représenter cette incertitude en considérant de multiples scénarios géologiques et en générant un grand nombre de réalisations géostatistiques. Nous rencontrons alors la principale limitation de ces approches qui est le coût de calcul dû à la simulation des processus d'écoulements complexes pour chacune de ces réalisations. Dans la première partie de la thèse, ce problème est investigué dans le contexte de propagation de l'incertitude, oú un ensemble de réalisations est identifié comme représentant les propriétés du sous-sol. Afin de propager cette incertitude à la quantité d'intérêt tout en limitant le coût de calcul, les méthodes actuelles font appel à des modèles d'écoulement approximés. Cela permet l'identification d'un sous-ensemble de réalisations représentant la variabilité de l'ensemble initial. Le modèle complexe d'écoulement est alors évalué uniquement pour ce sousensemble, et, sur la base de ces réponses complexes, l'inférence est faite. Notre objectif est d'améliorer la performance de cette approche en utilisant toute l'information à disposition. Pour cela, le sous-ensemble de réponses approximées et exactes est utilisé afin de construire un modèle d'erreur, qui sert ensuite à corriger le reste des réponses approximées et prédire la réponse du modèle complexe. Cette méthode permet de maximiser l'utilisation de l'information à disposition sans augmentation perceptible du temps de calcul. La propagation de l'incertitude est alors plus précise et plus robuste. La stratégie explorée dans le premier chapitre consiste à apprendre d'un sous-ensemble de réalisations la relation entre les modèles d'écoulement approximé et complexe. Dans la seconde partie de la thèse, cette méthodologie est formalisée mathématiquement en introduisant un modèle de régression entre les réponses fonctionnelles. Comme ce problème est mal posé, il est nécessaire d'en réduire la dimensionnalité. Dans cette optique, l'innovation du travail présenté provient de l'utilisation de l'analyse en composantes principales fonctionnelles (ACPF), qui non seulement effectue la réduction de dimensionnalités tout en maximisant l'information retenue, mais permet aussi de diagnostiquer la qualité du modèle d'erreur dans cet espace fonctionnel. La méthodologie proposée est appliquée à un problème de pollution par une phase liquide nonaqueuse et les résultats obtenus montrent que le modèle d'erreur permet une forte réduction du temps de calcul tout en estimant correctement l'incertitude. De plus, pour chaque réponse approximée, une prédiction de la réponse complexe est fournie par le modèle d'erreur. Le concept de modèle d'erreur fonctionnel est donc pertinent pour la propagation de l'incertitude, mais aussi pour les problèmes d'inférence bayésienne. Les méthodes de Monte Carlo par chaîne de Markov (MCMC) sont les algorithmes les plus communément utilisés afin de générer des réalisations géostatistiques en accord avec les observations. Cependant, ces méthodes souffrent d'un taux d'acceptation très bas pour les problèmes de grande dimensionnalité, résultant en un grand nombre de simulations d'écoulement gaspillées. Une approche en deux temps, le "MCMC en deux étapes", a été introduite afin d'éviter les simulations du modèle complexe inutiles par une évaluation préliminaire de la réalisation. Dans la troisième partie de la thèse, le modèle d'écoulement approximé couplé à un modèle d'erreur sert d'évaluation préliminaire pour le "MCMC en deux étapes". Nous démontrons une augmentation du taux d'acceptation par un facteur de 1.5 à 3 en comparaison avec une implémentation classique de MCMC. Une question reste sans réponse : comment choisir la taille de l'ensemble d'entrainement et comment identifier les réalisations permettant d'optimiser la construction du modèle d'erreur. Cela requiert une stratégie itérative afin que, à chaque nouvelle simulation d'écoulement, le modèle d'erreur soit amélioré en incorporant les nouvelles informations. Ceci est développé dans la quatrième partie de la thèse, oú cette méthodologie est appliquée à un problème d'intrusion saline dans un aquifère côtier. -- Our consumption of groundwater, in particular as drinking water and for irrigation, has considerably increased over the years and groundwater is becoming an increasingly scarce and endangered resource. Nofadays, we are facing many problems ranging from water prospection to sustainable management and remediation of polluted aquifers. Independently of the hydrogeological problem, the main challenge remains dealing with the incomplete knofledge of the underground properties. Stochastic approaches have been developed to represent this uncertainty by considering multiple geological scenarios and generating a large number of realizations. The main limitation of this approach is the computational cost associated with performing complex of simulations in each realization. In the first part of the thesis, we explore this issue in the context of uncertainty propagation, where an ensemble of geostatistical realizations is identified as representative of the subsurface uncertainty. To propagate this lack of knofledge to the quantity of interest (e.g., the concentration of pollutant in extracted water), it is necessary to evaluate the of response of each realization. Due to computational constraints, state-of-the-art methods make use of approximate of simulation, to identify a subset of realizations that represents the variability of the ensemble. The complex and computationally heavy of model is then run for this subset based on which inference is made. Our objective is to increase the performance of this approach by using all of the available information and not solely the subset of exact responses. Two error models are proposed to correct the approximate responses follofing a machine learning approach. For the subset identified by a classical approach (here the distance kernel method) both the approximate and the exact responses are knofn. This information is used to construct an error model and correct the ensemble of approximate responses to predict the "expected" responses of the exact model. The proposed methodology makes use of all the available information without perceptible additional computational costs and leads to an increase in accuracy and robustness of the uncertainty propagation. The strategy explored in the first chapter consists in learning from a subset of realizations the relationship between proxy and exact curves. In the second part of this thesis, the strategy is formalized in a rigorous mathematical framework by defining a regression model between functions. As this problem is ill-posed, it is necessary to reduce its dimensionality. The novelty of the work comes from the use of functional principal component analysis (FPCA), which not only performs the dimensionality reduction while maximizing the retained information, but also allofs a diagnostic of the quality of the error model in the functional space. The proposed methodology is applied to a pollution problem by a non-aqueous phase-liquid. The error model allofs a strong reduction of the computational cost while providing a good estimate of the uncertainty. The individual correction of the proxy response by the error model leads to an excellent prediction of the exact response, opening the door to many applications. The concept of functional error model is useful not only in the context of uncertainty propagation, but also, and maybe even more so, to perform Bayesian inference. Monte Carlo Markov Chain (MCMC) algorithms are the most common choice to ensure that the generated realizations are sampled in accordance with the observations. Hofever, this approach suffers from lof acceptance rate in high dimensional problems, resulting in a large number of wasted of simulations. This led to the introduction of two-stage MCMC, where the computational cost is decreased by avoiding unnecessary simulation of the exact of thanks to a preliminary evaluation of the proposal. In the third part of the thesis, a proxy is coupled to an error model to provide an approximate response for the two-stage MCMC set-up. We demonstrate an increase in acceptance rate by a factor three with respect to one-stage MCMC results. An open question remains: hof do we choose the size of the learning set and identify the realizations to optimize the construction of the error model. This requires devising an iterative strategy to construct the error model, such that, as new of simulations are performed, the error model is iteratively improved by incorporating the new information. This is discussed in the fourth part of the thesis, in which we apply this methodology to a problem of saline intrusion in a coastal aquifer.

Architectural characterization of a delta-front reservoir analogue combining ground penetrating radar and electrical resistivity tomography: Roda Sandstone (Lower Eocene, Graus-Tremp basin, Spain)

Three-dimensional reconstruction of reservoir analogues can be improved combining data from different geophysical methods. Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT) data are valuable tools, since they provide subsurface information from internal architecture and facies distribution of sedimentary rock bodies, enabling the upgrading of depositional models and heterogeneity reconstruction. The Lower Eocene Roda Sandstone is a well-known deltaic complex widely studied as a reservoir analogue that displays a series of sandstone wedges with a general NE to SW progradational trend. To provide a better understanding of internal heterogeneity of a 10m-thick progradational delta-front sandstone unit, 3D GPR data were acquired. In addition, common midpoints (CMP) to measure the sandstone subsoil velocity, test profiles with different frequency antennas (25, 50 and 100MHz) and topographic data for subsequent correction in the geophysical data were also obtained. Three ERT profiles were also acquired to further constrain GPR analysis. These geophysical results illustrate the geometry of reservoir analogue heterogeneities both depositional and diagenetic in nature, improving and complementing previous outcrop-derived data. GPR interpretation using radar stratigraphy principles and attributes analysis provided: 1)tridimensional geometry of major stratigraphic surfaces that define four units in the GPR Prism, 2) image the internal architecture of the units and their statistical study of azimuth and dips, useful for a quick determination of paleocurrent directions. These results were used to define the depositional architecture of the progradational sandbody that shows an arrangement in very-high-frequency sequences characterized by clockwise paleocurrent variations and decrease of the sedimentary flow, similar to those observed at a greater scale in the same system. This high-frequency sequential arrangement has been attributed to the autocyclic dynamics of a supply-dominated delta- front where fluvial and tidal currents are in competition. The resistivity models enhanced the viewing of reservoir quality associated with cement distribution caused by depositional and early diagenetic processes related to the development of transgressive and regressive systems tracts in igh-frequency sequences.

Effect of flushing frequency on emitter clogging in microirrigation with effluents

Flushing is an important maintenance task that removes accumulated particles in microirrigation laterals that can help to reduce clogging problems. The effect of three dripline flushing frequency treatments (no flushing, one flushing at the end of each irrigation period, and a monthly flushing during the irrigation period) was studied in surface and subsurface drip irrigation systems that operated using a wastewater treatment plant effluent for three irrigation periods of 540 h each. The irrigation systems had two different emitters, one pressure compensating and the other not, both molded and welded onto the interior dripline wall, placed in laterals 87 meters long. Dripline flow of the pressure compensating emitter increased 8% over time, while in the nonpressure compensating emitter, dripline flow increased 25% in the surface driplines and decreased 3% in the subsurface driplines by the emitter clogging. Emitter clogging was affected primarily by the interactions between emitter location, emitter type, and flushing frequency treatment. The number of completely clogged emitters was affected by the interaction between irrigation system and emitter type. There was an average of 3.7% less totally clogged emitters in flushed surface driplines with the pressurecompensating emitter as compared to flushed subsurface laterals with the nonpressure compensating emitter

Determinação da composição isotópica de estrôncio em águas naturais: exemplos de sua aplicação em águas subsuperficiais da zona costeira na região Bragantina-PA

Analytical procedures used for determining the concentrations and isotope composition of strontium in subsurface waters, by mass spectrometry, are described. Sampling was performed in coastal plateaus, salt marsh and mangrove environments in the coastal region of Pará. Coastal plateau waters have delta87Sr between 1.51 and 6.26‰ and Sr concentration bellow 58 ppb. Salt marsh waters show delta87Sr between 0.55 and 0.90‰ and Sr concentration between 93 and 114 ppm, while mangrove waters have delta87Sr‰ around zero and Sr concentration above 15 ppm. Differences in the 87Sr/86Sr ratio in these subsurface waters are detected, as well as seasonal variations in the coastal pleteau waters.

Estudi per a la construcció de l'EDAR mitjançant filtres verds de la Nou de Gaià (Tarragonès)

El tractament de les aigües en nuclis menors de 2000 habitants es troba pendent de completar per part de l’Agència Catalana de l’Aigua més concretament al corresponent Pla de Sanejament d’Aigües Residuals Urbanes (PSARU). El nucli de La Nou de Gaià (al Tarragonès) es troba pendent de la construcció de la corresponent instal·lació de sanejament, projectada al 2007. Alternativament a les depuradores tradicionals basades en l’ús de formigó (o materials alternatius) i en la despesa elèctrica per assegurar una aeració i una evacuació dels fangs generats, existeixen tecnologies “toves”. Aquestes tecnologies, també conegudes com a “verdes”, es basen en imitar els sistemes naturals maximitzant el seu potencial d’autodepuració. A grans trets existeixen dos formes de depurar les aigües de forma ecològica”: llacunatges (existeix una capa d’aigua lliure) i filtres verds. El present estudi es basa en l’aplicació de filtres verds de morfologia vertical i flux subsuperficial, plantat amb canyes dels generes Scirpus o Phragmites. El resultat han estat 4 bases de 35*35 per a tractar un cabal de 150 m3/d i una població equivalent de 1272.

Avaliação do desempenho de surfactantes para a solubilização de fases líquidas não aquosas em meio aquoso

The presence of non-aqueous phase liquids (NAPLs) in the subsurface is a threat to public health as well as a serious environmental issue. NAPLs may remain adsorbed or form lenses floating on aquifers causing long-term contaminations. Surfactants may increase NAPLs solubility, enhancing the pump-and-treatment performance. Size, shape, hydration and ionization degree of the micelles define the affinity and the space available for the solubilization of a particular contaminating agent. The tests carried out at laboratory scale, taking into account the NAPL to be removed and the medium characteristics were useful to select surfactants and evaluate their efficiency as NAPLs solubilizers.

Distribution of chemical compartments of soil organic matter and c stocks of a cambisol from south Brazil as affected by Pinus afforestation

Distribution and stocks of soil organic matter (SOM) compartments after Pinus monoculture introduction in a native pasture area of a Cambisol, Santa Catarina, Brazil, were investigated. Pinus introduction increased soil acidity, content of exchangeable Al+3 and diminished soil nutrients. Nevertheless, soil C stock increased in all humic fractions of the 0-5 cm layer after Pinus afforestation. In the subsurface, the vegetation change only promoted SOM redistribution from the NaOH-extractable humic substances to a less hydrophobic humin fraction. Under Pinus, soil organo-mineral interactions were relevant up to a 15 cm depth, while in pasture environment, this mechanism occurred mainly in the surface layer.

Genetic diversity among isolates of stemphylium solani from cotton

The fungus Stemphylium solani causes leaf blight of tomato (Lycopersicon esculentum) in Brazil. In recent years, severe epidemics of a new leaf blight of cotton (Gossipium hyrsutum) caused by S. solani occurred in three major cotton-growing Brazilian states (PR, MT and GO). Molecular analysis was performed to assess the genetic diversity among the S. solani isolates from cotton, and to verify their relationship with representative S. solani isolates from tomato. Random amplified polymorphic DNA (RAPD) markers and internal transcribed spacers of ribosomal DNA (rDNA) were used to compare 33 monosporic isolates of S. solani (28 from cotton and five from tomato). An isolate of Alternaria macrospora from cotton was also used for comparison. RAPD analysis showed the presence of polymorphism between the genera and the species. The A. macrospora and the S. solani isolates from cotton and tomato were distinct from each other, and fell into separate groups. Variation by geographic region was observed for the tomato isolates but not for the cotton isolates. Amplifications of the ITS region using the primer pair ITS4/ITS5 resulted in a single PCR product of approximately 600 bp for all the isolates. Similarly, when amplified fragments were digested with eight restriction enzymes, identical banding patterns were observed for all the isolates. Hence, rDNA analysis revealed no inter-generic or intra-specific variation. The genetic difference observed between the cotton and the tomato isolates provides evidence that S. solani attacking cotton in Brazil belongs to a distinct genotype.

Telomere and microsatellite primers reveal diversity among Sclerotinia sclerotiorum isolates from Brazil

Sclerotinia sclerotiorum, the causal agent of white mold, is a problem of winter bean (Phaseolus vulgaris) production in Brazil under center-pivot irrigation. Isolates of S. sclerotiorum were obtained from a center-pivot-irrigated field near Guaíra-SP, Brazil. Mycelial compatibility group (MCG) studies revealed the presence of only two MCG. PCR/RFLP analysis of the ITS1-5.8S-ITS2 ribosomal subunit regions of these field isolates of S. sclerotiorum failed to show any genetic differences between these two MCGs. DNA amplification with a chromosomal telomere sequence-based primer and one microsatellite primer revealed genetic polymorphisms among isolates within the same MCG. Isolates taken from beans and two other crops from another region of Brazil showed the same two MCG and had identical banding patterns for the telomere and microsatellite primers. These findings support the use of telomere sequence-based primers for revealing genotypic differences among S. sclerotiorum isolates.

Genetic analyses of Rhizoctonia solani isolates from Phaseolus vulgaris grown in the Atlantic Rainforest Region of São Paulo, Brazil

Rhizoctonia solani isolates obtained from common beans (Phaseolus vulgaris) grown in the mountainous Atlantic Rainforest (Mata Atlântica) region of São Paulo, Brazil, were analyzed to determine their genetic diversity using internal transcribed spacer (ITS), microsatellite and telomere sequence-based PCR primers. Restriction digestion of the ITS1/5.8S/ITS2 ribosomal regions yielded unique banding patterns specific for AG4 and its subgroups. The ITS restriction digestion (ITS/RFLP), telomere and microsatellite primers identified five to 11 genotypes within the isolates of R. solani. While all isolates were pathogenic on beans, there was no correlation found between genotypic differences and pathogenicity. The different PCR primers revealed a number of isolates that were genetically similar. Some of these genetic groups were supported by more than one of the primers utilized in this study, thus confirming their relationship.

Detection and coat protein gene characterization of an isolate of Grapevine virus B from corky bark-affected grapevines in Southern Brazil

An isolate of Grapevine virus B (GVB), obtained by indexing Vitis labrusca and V. vinifera grapevines on the indicator LN33, was transmitted mechanically to several Nicotiana species. The virus was partially purified from N. cavicola and the coat protein estimated at 23 kDa by SDS-PAGE. In negatively stained leaf extracts of experimentally inoculated N. cavicola and N. occidentalis, flexuous particles with cross banding were observed, predominantly measuring 750-770 x 12 nm, with a modal length of 760 nm. Decoration indicated a clear, positive reaction against AS-GVB. In DAS-ELISA, GVB was detected in N. cavicola and grapevine extracts, and Western blots showed homologous and cross reaction of GVB and GVA antisera with GVB coat protein. Using specific primers for GVB, a fragment of 594 bp, comprising the coat protein gene coding for 197 amino acids, was amplified by RT-PCR with viral RNA extracted from GVB-infected N. occidentalis. The nucleotide and the deduced amino acid sequences of the coat protein gene showed high identities with Italian and Japanese isolates of GVB.

Development of a scar marker for Pierce's Disease strains of Xylella fastidiosa

The objective of this research was to develop a primer for a polymerase chain reaction specific for Xylella fastidiosa strains that cause Pierce's Disease (PD) in grapes (Vitis vinifera). The DNA amplification of 23 different strains of X. fastidiosa, using a set of primers REP1-R (5'-IIIICGICGIATCCIGGC-3') and REP 2 (5'-ICGICTTATCIGGCCTAC-3') using the following program: 94 ºC/2 min; 35 X (94 ºC/1 min, 45 ºC/1 min and 72 ºC/1 min and 30 s) 72 ºC/5 min, produced a fragment of 630 bp that differentiated the strains that cause disease in grapes from the other strains. However, REP banding patterns could not be considered reliable for detection because the REP1-R and REP 2 primers correspond to repetitive sequences, which are found throughout the bacterial genome. The amplified product of 630 bp was eluted from the agarose gel, purified and sequenced. The nucleotide sequence information was used to identify and synthesize an specific oligonucleotide for X. fastidiosa strains that cause Pierce's Disease denominated Xf-1 (5'-CGGGGGTGTAGGAGGGGTTGT-3') which was used jointly with the REP-2 primer at the following conditions: 94 ºC/2 min; 35 X (94 ºC/1 min, 62 ºC/1 min; 72 ºC/1 min and 30 s) 72 ºC/10 min. The DNAs isolated from strains of X. fastidiosa from other hosts [almond (Prumus amygdalus), citrus (Citrus spp.), coffee (Coffea arabica), elm (Ulmus americana), mulberry (Morus rubra), oak (Quercus rubra), periwinkle wilt (Catharantus roseus), plums (Prunus salicina) and ragweed (Ambrosia artemisiifolia)] and also from other Gram negative and positive bacteria were submitted to amplification with a pair of primers Xf-1/REP 2 to verify its specificity. A fragment, about 350 bp, was amplified only when the DNA from strains of X. fastidiosa isolated from grapes was employed.