Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

The sampling of ignitable liquids on suspects' hands

In arson cases, the collection and detection of traces of ignitable liquids on a suspect's hands can provide information to a forensic investigation. Police forces currently lack a simple, robust, efficient and reliable solution to perform this type of swabbing. In this article, we describe a study undertaken to develop a procedure for the collection of ignitable liquid residues on the hands of arson suspects. Sixteen different collection supports were considered and their applicability for the collection of gasoline traces present on hands and their subsequent analysis in a laboratory was evaluated. Background contamination, consisting of volatiles emanating from the collection supports, and collection efficiencies of the different sampling materials were assessed by passive headspace extraction with an activated charcoal strip (DFLEX device) followed by gas chromatography-mass spectrometry (GC-MS) analysis. After statistical treatment of the results, non-powdered latex gloves were retained as the most suitable method of sampling. On the basis of the obtained results, a prototype sampling kit was designed and tested. This kit is made of a three compartment multilayer bag enclosed in a sealed metal can and containing three pairs of non-powdered latex gloves: one to be worn by the sampler, one consisting of a blank sample and the last one to be worn by the person suspected to have been in contact with ignitable liquids. The design of the kit was developed to be efficient in preventing external and cross-contaminations.

Crop sequences in no-tillage system: effects on soil fertility and soybean, maize and rice yield

Decomposing crop residues in no-tillage system can alter soil chemical properties, which may consequently influence the productivity of succession crops. The objective of this study was to evaluate soil chemical properties and soybean, maize and rice yield, grown in the summer, after winter crops in a no-tillage system. The experiment was carried out in Jaboticabal, SP, Brazil (21 ° 15 ' 22 '' S; 48 ° 18 ' 58 '' W) on a Red Latosol (Oxisol), in a completely randomized block design, in strip plots with three replications. The treatments consisted of four summer crop sequences (maize monocrop, soybean monocrop, soybean/maize rotation and rice/bean/cotton rotation) combined with seven winter crops (maize, sunflower, oilseed radish, pearl millet, pigeon pea, grain sorghum and sunn hemp). The experiment began in September 2002. After the winter crops in the 2005/2006 growing season and before the sowing of summer crops in the 2006/2007 season, soil samples were collected in the layers 0-2.5; 2.5-5.0; 5-10; 10-20; and 20-30 cm. Organic matter, pH, P, K+, Ca2+, Mg2+, and H + Al were determined in each soil sample. In the summer soybean/maize rotation and in maize the organic matter contents and P levels were lower, in the layers 0-10 cm and 0-20 cm, respectively. Summer rice/bean/cotton rotation increased soil K levels at 0-10 cm depth when sunn hemp and oilseed radish had previously been grown in the winter, and in the 0-2.5 cm layer for millet. Sunn hemp, millet, oilseed radish and sorghum grown in the winter increased organic matter contents in the soil down to 30 cm. Higher P levels were found at the depths 0-2.5 cm and 0-5 cm, respectively, when sunn hemp and oilseed radish were grown in the winter. Highest grain yields for soybean in monoculture were obtained in succession to winter oilseed radish and sunn hemp and in rotation with maize, after oilseed radish, sunn hemp and millet. Maize yields were highest in succession to winter oilseed radish, millet and pigeon pea. Rice yields were lowest when grown after sorghum.

A physical frequency-dependent compact model for RF integrated inductors

A frequency-dependent compact model for inductors in high ohmic substrates, which is based on an energy point-of-view, is developed. This approach enables the description of the most important coupling phenomena that take place inside the device. Magnetically induced losses are quite accurately calculated and coupling between electric and magnetic fields is given by means of a delay constant. The later coupling phenomenon provides a modified procedure for the computation of the fringing capacitance value, when the self-resonance frequency of the inductor is used as a fitting parameter. The model takes into account the width of every metal strip and the pitch between strips. This enables the description of optimized layout inductors. Data from experiments and electromagnetic simulators are presented to test the accuracy of the model.

Improvement of the quality factor of RF integrated inductors by layout optimization

A systematic method to improve the quality (Q) factor of RF integrated inductors is presented in this paper. The proposed method is based on the layout optimization to minimize the series resistance of the inductor coil, taking into account both ohmic losses, due to conduction currents, and magnetically induced losses, due to eddy currents. The technique is particularly useful when applied to inductors in which the fabrication process includes integration substrate removal. However, it is also applicable to inductors on low-loss substrates. The method optimizes the width of the metal strip for each turn of the inductor coil, leading to a variable strip-width layout. The optimization procedure has been successfully applied to the design of square spiral inductors in a silicon-based multichip-module technology, complemented with silicon micromachining postprocessing. The obtained experimental results corroborate the validity of the proposed method. A Q factor of about 17 have been obtained for a 35-nH inductor at 1.5 GHz, with Q values higher than 40 predicted for a 20-nH inductor working at 3.5 GHz. The latter is up to a 60% better than the best results for a single strip-width inductor working at the same frequency.

Phosphorus and root distribution and corn growth as related to long-term tillage systems and fertilizer placement

Soil and fertilizer management during cultivation can affect crop productivity and profitability. Long-term experiments are therefore necessary to determine the dynamics of nutrient and root distribution as related to soil profile, as well as the effects on nutrient uptake and crop growth. An 18-year experiment was conducted at the Federal University of Rio Grande do Sul State (UFRGS), in Eldorado do Sul, Brazil, on Rhodic Paleudult soil. Black oat and vetch were planted in the winter and corn in the summer. The soil management methods were conventional, involving no-tillage and strip tillage techniques and broadcast, row-and strip-applied fertilizer placement (triple superphosphate). Available P (Mehlich-1) and root distribution were determined in soil monoliths during the corn grain filling period. Corn shoot dry matter production and P accumulation during the 2006/2007 growing season were determined and the efficiency of P utilization calculated. Regardless of the degree of soil mobilization, P and roots were accumulated in the fertilized zone with time, mainly in the surface layer (0-10 cm). Root distribution followed P distribution for all tillage systems and fertilizer treatments. Under no-tillage, independent of the fertilizer placement, the corn plants developed more roots than in the other tillage systems. Although soil tillage systems and fertilizer treatments affected P and root distribution throughout the soil profile, as well as P absorption and corn growth, the efficiency of P utilization was not affected.

Availability of soil water under tillage systems, mulch management and citrus rootstocks

The increased availability of soil water is important for the management of non-irrigated orange orchards. The objective of this study was to evaluate the availability of soil water in a Haplorthox (Rhodic Ferralsol) under different tillage systems used for orchard plantation, mulch management and rootstocks in a "Pêra" orange orchard in northwest Paraná, Brazil. An experiment in a split-split-plot design was established in 2002, in an area cultivated with Brachiaria brizantha grass in which three tillage systems (no tillage, conventional tillage and strip-tillage) were used for orchard plantation. This grass was mowed twice a year between the rows, representing two mulch managements in the split plots (no mulching and mulching in the plant rows). The split-split-plots were represented by two rootstocks ("Rangpur" lime and "Cleopatra" mandarin). The soil water content in the plant rows was evaluated in the 0-20 cm layer in 2007 and at 0-20 and 20-40 cm in 2008-2009. The effect of soil tillage systems prior to implantation of orange orchards on soil water availability was less pronounced than mulching and the rootstocks. The soil water availability was lower when "Pêra" orange trees were grafted on "Cleopatra" mandarin than on "Rangpur" lime rootstocks. Mulching had a positive influence on soil water availability in the sandy surface layer (0-20 cm) and sandy clay loam subsurface (20-40 cm) of the soil in the spring. The production of B. brizantha between the rows and residue disposal in the plant rows as mulch increased water availability to the "Pêra" orange trees.

Mobility of inorganic and organic phosphorus forms under different levels of phosphate and poultry litter fertilization in soils

The eutrophication of aquifers is strongly linked to the mobility of P in soils. Although P mobility was considered irrelevant in a more distant past, more recent studies have shown that P, both in organic (Po) and inorganic forms (Pi), can be lost by leaching and eluviation through the soil profile, particularly in less weathered and/or sandier soils with low P adsorption capacity. The purpose of this study was to determine losses of P forms by leaching and eluviation from soil columns. Each column consisted of five PVC rings (diameter 5 cm, height 10 cm), filled with two soil types: a clayey Red-Yellow Latosol and a sandy loam Red-Yellow Latosol, which were exposed to water percolation. The soils were previously treated with four P rates (as KH2PO4 ) to reach 0, 12.5, 25.0 and 50 % of the maximum P adsorption capacity (MPAC). The P source was homogenized with the whole soil volume and incubated for 60 days. After this period the soils were placed in the columns; the soil of the top ring was mixed with five poultry litter rates of 0, 20, 40, 80, and 160 t ha-1 (dry weight basis). Treatments consisted of a 4 x 5 x 2 factorial scheme corresponding to four MPAC levels, five poultry litter rates, two soils, with three replications, arranged in a completely randomized block design. Deionized water was percolated through the columns 10 times in 35 days to simulate about 1,200 mm rainfall. In the leachate of each column the inorganic P (reactive P, Pi) and organic P forms (unreactive P, Po) were determined. At the end of the experiment, the columns were disassembled and P was extracted with the extractants Mehlich-1 (HCl 0.05 mol L-1 and H2SO4 0.0125 mol L-1) and Olsen (NaHCO3 0.5 mol L-1; pH 8.5) from the soil of each ring. The Pi and Po fractions were measured by the Olsen extractant. It was found that under higher poultry litter rates the losses of unreactive P (Po) were 6.4 times higher than of reactive P (Pi). Both the previous P fertilization and increasing poultry litter rates caused a vertical movement of P down the soil columns, as verified by P concentrations extracted by Mehlich-1 and NaHCO3 (Olsen). The environmental critical level (ECL), i.e., the P soil concentration above which P leaching increases exponentially, was 100 and 150 mg dm-3 by Mehlich-1 and 40 and 60 mg dm-3 by Olsen, for the sandy loam and clay soils, respectively. In highly weathered soils, where residual P is accumulated by successive crops, P leaching through the profile can be significant, particularly when poultry litter is applied as fertilizer.

Groundwater Availability Modeling of the Lower Dakota Aquifer in Northwest Iowa, October 2008

The 2007 Iowa General Assembly, recognizing the increased demand for water to support the growth of industries and municipalities, approved funding for the first year of a multi-year evaluation and modeling of Iowa’s major aquifers by the Iowa Department of Natural Resources. The task of conducting this evaluation and modeling was assigned to the Iowa Geological and Water Survey (IGWS). The first aquifer to be studied was the Lower Dakota aquifer in a sixteen county area of northwest Iowa.

Soil microbial biomass under different management and tillage systems of permanent intercropped cover species in an orange orchard

To mitigate soil erosion and enhance soil fertility in orange plantations, the permanent protection of the inter-rows by cover species has been suggested. The objective of this study was to evaluate alterations in the microbial biomass, due to different soil tillage systems and intercropped cover species between rows of orange trees. The soil of the experimental area previously used as pasture (Brachiaria humidicola) was an Ultisol (Typic Paleudult) originating from Caiuá sandstone in the northwestern part of the State of Paraná, Brazil. Two soil tillage systems were evaluated: conventional tillage (CT) in the entire area and strip tillage (ST) (strip width 2 m), in combination with different ground cover management systems. The citrus cultivar 'Pera' orange (Citrus sinensis) grafted onto 'Rangpur' lime rootstock was used. Soil samples were collected after five years of treatment from a depth of 0-15 cm, under the tree canopy and in the inter-row, in the following treatments: (1) CT and an annual cover crop with the leguminous species Calopogonium mucunoides; (2) CT and a perennial cover crop with the leguminous peanut Arachis pintoi; (3) CT and an evergreen cover crop with Bahiagrass Paspalum notatum; (4) CT and a cover crop with spontaneous Brachiaria humidicola grass vegetation; and (5) ST and maintenance of the remaining grass (pasture) of Brachiaria humidicola. Soil tillage and the different cover species influenced the microbial biomass, both under the tree canopy and in the inter-row. The cultivation of brachiaria increased C and N in the microbial biomass, while bahiagrass increased P in the microbial biomass. The soil microbial biomass was enriched in N and P by the presence of ground cover species and according to the soil P content. The grass species increased C, N and P in the soil microbial biomass from the inter-row more than leguminous species.

Airborne GPS, 1996

The paper gives the theory of airborne GPs related to Photogrammetry and the results of a self calibration used to validate the theory. Accordingly, no ground control points are required for mapping using a strip or block of photographs provided the site is within 10 Km of the calibration site.

Estimation of water percolation by different methods using TDR

Detailed knowledge on water percolation into the soil in irrigated areas is fundamental for solving problems of drainage, pollution and the recharge of underground aquifers. The aim of this study was to evaluate the percolation estimated by time-domain-reflectometry (TDR) in a drainage lysimeter. We used Darcy's law with K(θ) functions determined by field and laboratory methods and by the change in water storage in the soil profile at 16 points of moisture measurement at different time intervals. A sandy clay soil was saturated and covered with plastic sheet to prevent evaporation and an internal drainage trial in a drainage lysimeter was installed. The relationship between the observed and estimated percolation values was evaluated by linear regression analysis. The results suggest that percolation in the field or laboratory can be estimated based on continuous monitoring with TDR, and at short time intervals, of the variations in soil water storage. The precision and accuracy of this approach are similar to those of the lysimeter and it has advantages over the other evaluated methods, of which the most relevant are the possibility of estimating percolation in short time intervals and exemption from the predetermination of soil hydraulic properties such as water retention and hydraulic conductivity. The estimates obtained by the Darcy-Buckingham equation for percolation levels using function K(θ) predicted by the method of Hillel et al. (1972) provided compatible water percolation estimates with those obtained in the lysimeter at time intervals greater than 1 h. The methods of Libardi et al. (1980), Sisson et al. (1980) and van Genuchten (1980) underestimated water percolation.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Water Balance Components in Covered and Uncovered Soil Growing Irrigated Muskmelon

ABSTRACT Knowledge of the terms (or processes) of the soil water balance equation or simply the components of the soil water balance over the cycle of an agricultural crop is essential for soil and water management. Thus, the aim of this study was to analyze these components in a Cambissolo Háplico (Haplocambids) growing muskmelon (Cucumis melo L.) under drip irrigation, with covered and uncovered soil, in the municipality of Baraúna, State of Rio Grande do Norte, Brazil (05º 04’ 48” S, 37º 37’ 00” W). Muskmelon, variety AF-646, was cultivated in a flat experimental area (20 × 50 m). The crop was spaced at 2.00 m between rows and 0.35 m between plants, in a total of ten 50-m-long plant rows. At points corresponding to ⅓ and ⅔ of each plant row, four tensiometers (at a distance of 0.1 m from each other) were set up at the depths of 0.1, 0.2, 0.3, and 0.4 m, adjacent to the irrigation line (0.1 m from the plant row), between two selected plants. Five random plant rows were mulched using dry leaves of banana (Musa sp.) along the drip line, forming a 0.5-m-wide strip, which covered an area of 25 m2 per of plant row with covered soil. In the other five rows, there was no covering. Thus, the experiment consisted of two treatments, with 10 replicates, in four phenological stages: initial (7-22 DAS - days after sowing), growing (22-40 DAS), fruiting (40-58 DAS) and maturation (58-70 DAS). Rainfall was measured with a rain gauge and water storage was estimated by the trapezoidal method, based on tensiometer readings and soil water retention curves. For soil water flux densities at 0.3 m, the tensiometers at the depths of 0.2, 0.3, and 0.4 m were considered; the tensiometer at 0.3 m was used to estimate soil water content from the soil water retention curve at this depth, and the other two to calculate the total potential gradient. Flux densities were calculated through use of the Darcy-Buckingham equation, with hydraulic conductivity determined by the instantaneous profile method. Crop actual evapotranspiration was calculated as the unknown of the soil water balance equation. The soil water balance method is effective in estimating the actual evapotranspiration of irrigated muskmelon; there was no significant effect of soil coverage on capillary rise, internal drainage, crop actual evapotranspiration, and muskmelon yield compared with the uncovered soil; the transport of water caused by evaporation in the uncovered soil was controlled by the break in capillarity at the soil-atmosphere interface, which caused similar water dynamics for both management practices applied.

Estimate of Groundwater Recharge Based on Water Balance in The Unsaturated Soil Zone

ABSTRACT Groundwater management depends on the knowledge on recharge rates and water fluxes within aquifers. The recharge is one of the water cycle components most difficult to estimate. As a result, despite the chosen method, the estimates are subject to uncertainties that can be identified by means of comparison with other approaches. In this study, groundwater recharge estimates based on the water balance in the unsaturated zone is assessed. Firstly, the approach is evaluated by comparing the results with those of another method. Then, the estimates are used as inputs in a transient groundwater flow model in order to assess how the water table would respond to the obtained recharges rates compared to measured levels. The results suggest a good performance of the adopted approach and, despite some inherent limitations, it has advantages over other methods since the data required are easier to obtain.