Analysis of flow and solute transport problem in porous media are affected by uncertainty inbuilt both in boundary conditions and spatial variability in system parameters. The experimental investigation reveals that the parameters may vary in various scales by several orders. These affect the solute plume characteristics in field-scale problem and cause uncertainty in the prediction of concentration. The main focus of the present thesis is to analyze the probabilistic behavior of solute concentration in three dimensional(3-D) heterogeneous porous media. The framework for the probabilistic analysis has been developed using perturbation approach for both spectral based analytical and finite element based numerical method. The results of the probabilistic analysis are presented either in terms of solute plume characteristics or prediction uncertainty of the concentration. After providing a brief introduction on the role of stochastic analysis in subsurface hydrology in chapter 1, a detailed review of the literature is presented to establish the existing state-of-art in the research on the probabilistic analysis of flow and transport in simple and complex heterogeneous porous media in chapter 2. The literature review is mainly focused on the methods of solution of the stochastic differential equation. Perturbation based spectral method is often used for probabilistic analysis of flow and solute transport problem. Using this analytical method a nonlocal equation is solved to derive the expression of the spatial plume moments. The spatial plume moments represent the solute movement, spreading in an average sense. In chapter 3 of the present thesis, local dispersivity if also assumed to be random space function along with hydraulic conductivity. For various correlation coefficients of the random parameters, the results in terms of the field scale effective dispersivity are presented to demonstrate the effect of local dispersivity variation in space. The randomness of local dispersivity is found to reduce the effective fields scale dispersivity. The transverse effective macrodispersivity is affected more than the longitudinal effective macrodispersivity due to random spatial variation of local dispersivity. The reduction in effective field scale longitudinal dispersivity is more for positive correlation coefficient. The applicability of the analytical method, which is discussed in earlier chapter, is limited to the simple boundary conditions. The solution by spectral method in terms of statistical moments of concentration as a function of space and time, require higher dimensional integration. Perturbation based stochastic finite element method(SFEM) is an alternative method for performing probabilistic analysis of concentration. The use of this numerical method for performing probabilistic analysis of concentration. The use of this numerical method is non common in the literature of stochastic subsurface hydrology. The perturbation based SFEM which uses FEM for spatial discretization of the steady state flow and Laplace transform for the solute transport equation, is developed in chapter 4. The SFEM is formulated using Taylor series of the dependent variable upto second-order term. This results in second-order accurate mean and first-order accurate standard deviation of concentration. In this study the governing medium properties viz. hydraulic Conductivity, dispersivity, molecular diffusion, porosity, sorption coefficient and decay coefficient are considered to vary randomly in space. The accuracy of results and computational efficiency of the SFEM are compared with Monte Carle Simulation method(MCSM) for both I-D and 3-D problems. The comparison of results obtained hby SFEM and MCSM indicates that SFEM is capable in providing reasonably accurate mean and standard deviation of concentration. The Laplace transform based SFEM is simpler and advantageous since it does not require any stability criteria for choosing the time step. However it is not applicable for nonlinear transport problems as well as unsteady flow conditions. In this situation, finite difference method is adopted for the time discretization. The first part of the Chapter 5, deals with the formulation of time domain SFEM for the linear solute transport problem. Later the SFEM is extended for a problem which involve uncertainty of both system parameters and boundary/source conditions. For the flow problem, the randomness in the boundary condition is attributed by the random spatial variation of recharge at the top of the domain. The random recharge is modeled using mean, standard deviation and 2-D spatial correlation function. It is observed that even for the deterministic recharge case, the behavior of prediction uncertainty of concentration in the space is affected significantly due to the variation of flow field. When the effect of randomness of recharge condition is included, the standard deviation of concentration increases further. For solute transport, the concentration input at the source is modeled as a time varying random process. Two types of random source at the source is modeled as a time varying random process. Two types of random source condition are considered, firstly the amount of solute mass released at uniform time interval is random and secondly the source is treated as a Poission process. For the case of multiple random mass releases, the stochastic response function due to stochastic system is obtained by using SFEM. Comparing the results for the two type of random sources, it sis found that the prediction uncertainty is more when it is modeled as a Poisson process. The probabilistic analysis of nonlinear solute transport problem using MCSM is often requires large computational cost. The formulation of the alternative efficient method, SFEM, for nonlinear solute transport problem is presented in chapter 6. A general Langmuir-Freundlich isotherm is considered to model the equilibrium mass transfer between aqueous and sorbed phase. In the SFEM formulation, which uses the Taylor Series expansion, the zeroth-order derivatives of concentration are obtained by solving nonlinear algebraic equation. The higher order derivatives are obtained by solving linear equation. During transport, the nonlinear sorbing solutes is characterized by sharp solute fronts with a traveling wave behavior. Due to this the prediction uncertainty is significantly higher. The comparison of accuracy and computational efficiency of SFEM with MCSM for I-D and 3-D problems, reveals that the performance of SFEM for nonlinear problem is good and similar to the linear problem. In Chapter 7, the nonlinear SFEM is extended for probabilistic analysis of biodegrading solute, which is modeled by a set of PDEs coupled with nonlinear Monod type source/sink terms. In this study the biodegradation problem involves a single solute by a single class of microorganisms coupled with dynamic microbial growth is attempted using this methods. The temporal behavior of mean and standard deviation of substrate concentration are not monotonic, they show peaks before reaching lower steady state value. A comparison between the SFEM and MCSM for the mean and standard deviation of concentration is made for various stochastic cases of the I-D problem. In most of the cases the results compare reasonably well. The analysis of probabilistic behavior of substrate concentration for different correlation coefficient between the physical parameters(hydraulic conductivity, porosity, dispersivity and diffusion coefficient) and the biological parameters(maximum substrate utilization rate and the coefficient of cell decay) is performed. It is observed that the positive correlation between the two sets of parameters results in a lower mean and significantly higher standard deviation of substrate concentration. In the previous chapters, the stochastic analysis pertaining to the prediction uncertainty of concentration has been presented for simple problem where the system parameters are modeled as statistically homogeneous random. The experimental investigations in a small watershed, point towards a complex in geological substratum. It has been observed through the 2-D electrical resistivity imaging that the interface between the layers of high conductive weathered zone and low conductive clay is very irregular and complex in nature. In chapter 8 a theoretical model based on stochastic approach is developed to stimulate the complex geological structure of the weathered zone, using the 2-D electrical image. The statistical parameters of hydraulic conductivity field are estimated using the data obtained from the Magnetic Resonance Sounding(MRS) method. Due to the large complexity in the distribution of weathered zone, the stochastic analysis of seepage flux has been carried out by using MCSM. A batter characterization of the domain based on sufficient experimental data and suitable model of the random conductivity field may help to use the efficient SFEM. The flow domain is modeled as (i) an unstructured random field consisting of a single material with spatial heterogeneity, and (ii) a structured random field using 2-D electrical imaging which is composed of two layers of different heterogeneous random hydraulic properties. The simulations show that the prediction uncertainty of seepage flux is comparatively less when structured modeling framework is used rather than the unstructured modeling. At the end, in chapter 9 the important conclusions drawn from various chapters are summarized.

<p>In the recent decades, therehave been frequent conflicts between groundwater waterresources and environmentally hazardous activities. Newmethodologies for aiding decision-making in groundwater impactassessment and protection areneeded and in which issues ofincreased awareness, better understanding of the groundwaterresources processes, and validation of predictive mathematicalmodels are addressed.</p><p>A framework fordecision–aid, based on predictive simulations that a)predicts the environmental impacts b) provides the totaleconomical value c) visualises the impacts and the groundwaterproperties and d) describes the uncertainties in the results isproposed herein. The framework can be applied in environmentalimpact assessments, strategic environmental assessments andprotection and management of water resources. The results ofthe model are used as feedback for determining new scenarios,depending on the required uncertainties, and if the plannedactivity is sustainable, and/or fulfils the legislative andpolicy measures. This framework is applied to a particular casestudy, Nybroåsen, in the south-eastern part of Sweden,where the highway E22 is constructed through the importantglaciofluvial esker aquifer, passing the protection zone of thewater supply for the Kalmar municipality.</p><p>The impacts from the new highwayand the existing road have been predicted by two-dimensionalphysically based time-variant flow and solute groundwatermodelling. The results, breakthrough curves of contaminantconcentration in wells and maps of concentration distributions,as well as travel times, flow paths, and capture zones forwells determined by particle tracking have been presented.</p><p>The constructed model of theNybroåsen study area was calibrated by comparing observedand simulated groundwater levels for 15 observation wells forten years of measurements. The model has been evaluated bothgraphically and numerically and the calibration target wasfulfilled for 13 of the 15 observation wells. The model workincludes investigations of the catchment information, a waterbalance study, simulation of the groundwater recharge,consideration of the unsaturated zone by a numerical columnsimulation, and sensitivity analysis.</p><p>From the sensitivity analysis ofthe flow and transport parameters, it has been shown that theuncertainties are mainly due to the hydraulic conductivity.Comparison of the derived conductivity from the steady-stateautomatic calibration and the time-variant calibration showedthat there are major differences in the derived parameters,which illustrates the importance of a time dependentcalibration over both wet and dry periods and in more than onepoint in the area of interest of the model predictions.</p><p>In addition, a multi-criteriadecision analysis has been carried out for four roadalternatives (including the new highway E22) and the existingroad in the case study concerned. The multi-criteria decisionaid is applied as an illustration of how it can be used in thestudy area to identify a) interest groups of actors and theirconcerns b) ranking of alternative road scenarios according toactors’preferences and c) coalition groups of actors<i>i.e.</i>groups that have similar views with regard to theroad alternatives.</p><p><b>Keywords:</b>Physically-based groundwater modelling,contamination, flow and solute transport, glaciofluvialdeposits, Nybroåsen, Sweden, and multi-criteriadecision-aid.</p>

Transit times through hydrologic systems vary in time, but the nature of that variability is not well understood. Transit times variability was investigated in a 1 m(3) sloping lysimeter, representing a simplified model of a hillslope receiving periodic rainfall events for 28 days. Tracer tests were conducted using an experimental protocol that allows time-variable transit time distributions (TTDs) to be calculated from data. Observed TTDs varied with the storage state of the system, and the history of inflows and outflows. We propose that the observed time variability of the TTDs can be decomposed into two parts: internal variability associated with changes in the arrangement of, and partitioning between, flow pathways; and external variability driven by fluctuations in the flow rate along all flow pathways. These concepts can be defined quantitatively in terms of rank StorAge Selection (rSAS) functions, which is a theory describing lumped transport dynamics. Internal variability is associated with temporal variability in the rSAS function, while external is not. The rSAS function variability was characterized by an inverse storage effect, whereby younger water is released in greater proportion under wetter conditions than drier. We hypothesize that this effect is caused by the rapid mobilization of water in the unsaturated zone by the rising water table. Common approximations used to model transport dynamics that neglect internal variability were unable to reproduce the observed breakthrough curves accurately. This suggests that internal variability can play an important role in hydrologic transport dynamics, with implications for field data interpretation and modeling.

Fluoroquinolones (FQ) are broad-spectrum charged antimicrobials exhibiting excellent tissue/fluid permeation. Thus, FQ disposition depends essentially on active transport and facilitative diffusion. Although most early transporter studies investigating renal elimination of FQs have focused on apical efflux of FQs from renal proximal tubule cell (RPTC) into urine, their basolateral uptake mechanism(s) from blood into RPTC (i.e., first step to tubular secretion) has not yet been explored in detail. Renally expressed SLC22 members: organic anion (OATs) and cation (OCTs) transporters are known to transport such small organic ionic substrates (molecular weight ~400 Da). Hence it is of interest to explore the role of these basolateral transporters in renal elimination of FQs, and to further quantitatively assess their impact in clinically observed FQ drug-drug interactions (DDI). An initial systematic review of clinical literature for FQs (n=18) demonstrated substantial differences among their renal clearance (CLren~46-fold) and unbound renal clearance (CLrenu~20-fold), and suggested that tubular secretion and reabsorption could be major determinants of FQ half-life, efficacy, and DDIs. FQs (n=13) identified from the above review were investigated by in-vitro transport studies using stably transfected cell lines, for potential interactions with organic cation [human (h) OCT1, hOCT2 and hOCT3] and anion [mouse (m) and hOAT3, hOAT1; and hOAT4] transporters. Further, kinetic inhibition studies were conducted to determine inhibition potency (Ki/IC50 values) for those FQs exhibiting significant OCT/OAT inhibition in preliminary interaction experiments. Gatifloxacin, moxifloxacin, prulifloxacin, and sparfloxacin were determined to be competitive inhibitors of hOCT1 with Ki = 250±18, 161±19, 136±33, and 94±8 μM, respectively. Moxifloxacin competitively inhibited hOCT3-mediated uptake, Ki = 1,598±146 μM. Enoxacin, fleroxacin, levofloxacin, lomefloxacin, moxifloxacin, prulifloxacin, and sparfloxacin exhibited competitive inhibition for mOat3 with Ki = 396±15, 817±31, 515±22, 539±27, 1356±114, 299±35, 205±12 μM, respectively. Fleroxacin and pefloxacin were found to inhibit hOAT1 with IC50 = 2228±84 and 1819±144 respectively. Despite expression in enterocytes, hepatocytes, and RPTC, hOCT3 does not appear to contribute significantly to FQ disposition. However, due to hepatic and potential RPTC expression, hOCT1 could play an important role in elimination of these antimicrobials. Among renally expressed OATs in humans, hOAT1 and hOAT3 are likely to be involved in FQ elimination.

This project was undertaken to gain a better understanding of the hydration behaviors of gas phase ions from solutions containing amino acids and peptides. In order to characterize their hydration behavior, the molecules of interest in solutions were first converted into gas phase ions by electrospray ionization (ESI). The completely desolvated ions were then deliberately dispersed into an inert bath gas, usually nitrogen, containing accurately known concentrations of solvent vapor. The resulting mixtures of ions and bath gas were subsequently passed into a vacuum chamber by way of an adiabatic supersonic free jet expansion. The cooling during that expansion caused solvation of the ions, the extent of which was determined by a quadrupole mass analyzer. Mass analysis of the solute ions in the absence of vapor showed peaks with the mass to charge ratios corresponding to the desolvated ions. On the other hand, mass spectrometric analyses of ions in the presence of solvent vapor showed sequences of peaks corresponding to the solvated ions with varying numbers of water molecules. The extent of the ion solvation was controlled by varying the concentration of solvent vapor in the bath gas. Two different scales were proposed for the evaluation of the relative affinities of amino acids for water molecules. One was based primarily on the assumption that the affinities of amino acids for water molecules are directly proportional to their gas phase solvation rate constants (k). An alternative approach produced an affinity scale based on the extent of ion hydration occurred during the free jet expansion. It was found that the addition of a polar solvent vapor to the bath gas at low concentrations substantially enhanced the production of the bare solute ions from the evaporating charged droplets. This remarkable result not only provided a means to increase the ion production and thus detection sensitivity of mass spectrometric analyses, but also yielded important information regarding the ion formation mechanism of ESI. Additional studies revealed that the extent of the increase in ion yield was directly related to the charge state and molecular weight of the solute ions. In sum, this evidence strongly indicated that gas phase ions produced from charged droplets, as in electrospray ionization, must proceed by the sequence of events assumed in the Ion Evaporation Model proposed by Iribarne and Thomson rather than in the Charged Residue Model originally proposed by Malcolm Dole and coworkers. The hydration behaviors of electrospray ions from peptides with similar primary amino acid sequences and capable of forming ions with more than one charge state were also investigated. In a study with dipeptides, the extent of hydration was found to vary widely and to depend not only on the chemical composition of the ions but also on their configurations and charge states. The results obtained with lysine oligomers clearly indicated that the number of charges on an ion played an important role in the solvation process. An exception to this generalization was found in an experiment with multiply protonated pentalysine ions. For example, the quadruply protonated monomers of that species were found to undergo charge reduction via proton exchange with the surrounding water molecules in such a way as to maximize the distance between charges on the molecule, thereby reducing the internal repulsive forces.The hydration study of angiotensin II and III showed that while the former has an additional hydrophilic amino acid on the N-terminus, the latter peptide was more hydrophilic. This result suggests that the hydrophilicities of peptides are not a simple sum of the hydrophilicities of the individual amino acid components. As further evidence of interaction complexity, the Magic Number Clusters containing 21 water molecules were obtained with the doubly protonated angiotensin III, but not with the doubly protonated angiotensin II. Taken together, these observations seem to indicate that the multiply charged ions of angiotensin II and III had different structural conformations.

Experiments using electrical resistivity tomography (ERT) have shown promising results in reducing the uncertainty of solute plume characteristics related to estimates based on the analysis of local point measurements only. To explore the similarities and differences between two cross-borehole ERT inversion approaches for characterizing salt tracer plumes, namely the classical smoothness-constrained inversion and a geostatistically based approach, we performed two-dimensional synthetic experiments. Simplifying assumptions about the solute transport model and the electrical forward and inverse model allowed us to study the sensitivity of the ERT inversion approaches towards a variety of basic conditions, including the number of boreholes, measurement schemes, contrast between the plume and background electrical conductivity, use of a priori knowledge, and point conditioning. The results show that geostatistically based and smoothness-constrained inversions of electrical resistance data yield plume characteristics of similar quality, which can be further improved when point measurements are incorporated and advantageous measurement schemes are chosen. As expected, an increased number of boreholes included in the ERT measurement layout can highly improve the quality of inferred plume characteristics, while in this case the benefits of point conditioning and advantageous measurement schemes diminish. Both ERT inversion approaches are similarly sensitive to the noise level of the data and the contrast between the solute plume and background electrical conductivity, and robust with regard to biased input parameters, such as mean concentration, variance, and correlation length of the plume. Although sophisticated inversion schemes have recently become available, in which flow and transport as well as electrical forward models are coupled, these schemes effectively rely on a relatively simple geometrical parameterization of the hydrogeological model. Therefore, we believe that standard uncoupled ERT inverse approaches, like the ones discussed and assessed in this paper, will continue to be important to the imaging and characterization of solute plumes in many real-world applications. (C) 2016 Hohai University. Production and hosting by Elsevier B.V.

This study focuses on numerical modeling of fluid flow and solute transport in rough-walled rock fractures and fracture-matrix systems, with the main aim to investigate the impacts of fracture surface roughness on flow and transport processes in rock fractures. Both 2D and 3D fracture models were built from laser-scanned surface tomography of a real granite rock sample, to consider realistic features of surface tomography and potential asperity contacts. The flow was simulated by directly solving the Navier-Stokes equations (NSE) and the transport was modeled by solving the advection-dispersion equation (ADE) in the entire domain of fracture-matrix system, including matrix diffusion process. Such direct simulations provided detailed flow and concentration fields for quantitatively analysis of flow and transport behavior. The detailed analysis of surface roughness decomposition, complex flow patterns (i.e., channeling, transverse and eddy flows), effective advective flow apertures, effective transmissivity, effective dispersivity, residence time, transport resistance and specific surface area demonstrated significant impacts of realistic fracture surface roughness on fluid flow and solute transport processes in rock fractures. The results show that the surface roughness and shear displacement caused asperity contacts significantly enhance nonlinearity and complexity of flow and transport processes in rough-walled fractures and fracture-matrix systems. The surface roughness also causes invasion flows in intersected fractures which enhance solute mixing at fracture intersections. Therefore, the fracture surface roughness is an important source of uncertainty in application of such simplified models like cubic law (CL) for fluid flow and analytical solutions for solute transport in rock fractures. The research conducted advances our understanding of realistic flow and transport processes in natural fractured rocks. The results are useful for model validation/extension, uncertainty analysis/quantification and laboratory experiments design in the context of various applications related to fracture flow and transport. / Denna studie fokuserar på numerisk modellering av vätskeflöde och transport av lösta ämnen i frakturer med ojämna väggar samt fraktur-matrissystem, med det huvudsakliga syftet att undersöka effekterna av frakturernas ytjämnhet på flödes- och transportprocesser i bergsfrakturer. Både 2D och 3D modeller skapades utifrån laser skannad tomografi av ett verkligt bergartsprov av granit, för att överväga de realistiska egenskaperna hos ytan och potentiell skrovlighet. Flödet simulerades genom att lösa Navier-Stokes ekvationer (NSE) och transporten modellerades genom att lösa advektion-dispersion ekvation (ADE) i hela domänen av fraktur-matrissystemet, inklusive diffusions process i matrisen. Sådana direkta simuleringar resulterade i detaljerade flödes- och koncentrationsfält för att kvantitativt kunna analysera flödet och transportbeteendet. En detaljerad analys av upplösningen av ytjämnhet, komplexa flödesmönster (dvs kanalisering, tvärgående och virvelströmmar), effektiv advektiv flödesöppning, effektiv transmissivitet, effektiv dispersivitet, uppehållstid, transport motstånd och specifik yta visade signifikanta effekter av realistiska ojämna frakturväggar på vätskeflöde och lösta transportprocesser i bergssprickor. Resultaten visar att ytjämnhet och skjuvningssystemsorsakade asperitetskontakter avsevärt förbättrar olinjäritet och komplexitet av flödes- och transportprocesser i frakturer med ojämna väggar samt fraktur-matrissystem. Ytråheten orsakar också intrång av flöde i tvärgående frakturer vilket ökar blandingen av lösta ämnen i korsningarna. Därför är ytjämnhet av frakturerna en viktig källa till osäkerhet i tillämpningen av sådana förenklade modeller som kubisk lag (CL) för vätskeflöde och analytiska lösningar för transport av lösta ämnen i bergsfrakturer. Studien har ökat förståelsen för realistiska flödes- och transportprocesser i naturligt sprucket berg. Resultaten är användbara för modellvalidering/förlängning, osäkerhetsanalys/kvantifiering och design av laboratorieexperiment i samband med olika tillämpningar av flöde och transport i bergsfrakturer. / <p>QC 20161010</p>

O presente trabalho teve por objetivo estudar o tempo necessário para que o fluxo da água no solo sob condições de saturação atinja o equilíbrio dinâmico, medindose a condutividade hidráulica saturada (K0) do solo. A condição de equilíbrio dinâmico é a etapa que antecede a aplicação do íon no solo em estudos de eluição de solutos. Utilizaram-se amostras deformadas de três solos com diferentes granulometrias: muito argilosa, média e arenosa, coletadas em área experimental da ESALQ/USP. As amostras foram acondicionadas em colunas de PVC, e utilizou-se água destilada e deaerada para saturação e posterior escoamento nas colunas, conforme a montagem experimental do permeâmetro de carga constante. Os resultados demonstraram que para as amostras dos solos muito argiloso e médio, o equilíbrio foi atingido após 15 dias de drenagem e para o arenoso essa condição só foi verificada após 27 dias, sendo a variabilidade de K0 maior nos primeiros seis dias de ensaio. Concluindo-se, portanto pela inadequação da utilização de um tempo fixo para a adição dos solutos. / The purpose of this work is to study the time for the water flow in soil under saturation to reach the steady-state, by measuring the saturated hydraulic conductivity (K0) of soil. The steady-state is the stage that precedes the application of the ion solution in the soil sample in solute elution studies. Disturbed soil samples with different granulometries were used: very clayey, clay-sandy and sandy soil, collected in experimental area of the ESALQ/USP. The samples were conditioned in PVC columns, and distilled and deaerated water were used for saturation and subsequent drainage in the columns, the test is in accordance with the experimental assembly of constant head permeameter. The data demonstrated that for very clayey and clay-sandy samples, the steady-state was reached after 15 days and for sandy soil this condition was verified after 27 days, being the variability of K0 bigger in the first six days of the test. Therefore it was concluded that is inadequate to use a fixed time for the addition of solutes in soil.

Saline water evaporation from porous media with the associated salt precipitation patterns is frequently observed in a number of industrial and environmental applications and it is important in a variety of topics including, but not limited to, water balance and land-atmosphere interaction, terrestrial ecosystem functioning, geological carbon storage, and preservation of historical monuments. The excess accumulation of salt in soil is a global problem and is one of the most widespread soil degradation processes. Thus, it is important to understand the dominant mechanisms controlling saline water evaporation from porous media. This process is controlled by the transport properties of the porous medium, the external conditions, and the properties of the evaporating fluid. During saline water evaporation from porous media, the capillary induced liquid flow transports the solute towards the evaporation surface while diffusive transport tends to spread the salt homogeneously thorough the porous medium. Therefore, the solute distribution is influenced by the competition between the diffusive and convective transport. As water evaporates, salt concentration in the pore space increases continually until it precipitates. The formation of precipitated salt adds to the complexity of the description of saline water evaporation from porous media. In this dissertation, the effects of salt concentration, type of salt, and the presence of precipitated salt, on the evaporation dynamics have been investigated. The obtained results show that the precipitated salt has a porous structure and it evolves as the drying progresses. The presence of porous precipitated salt at the surface causes top-supplied creeping of the evaporating solution, feeding the growth of subsequent crystals. This could be visualized by thermal imaging in the form of appearance and disappearance of cold-spots on the surface of the porous medium, brought about by preferential water evaporation through the salt crust. My results show that such a phenomenon influences the dynamics of saline water evaporation from porous media. Moreover, a simple but effective tool was developed in this dissertation capable of describing the effects of ambient temperature, relative humidity, type of salt and its concentration, on the evaporative fluxes. Additionally, pore-scale data obtained by synchrotron x-ray tomography was used to study ion transport during saline water evaporation from porous media in 4D (3D space + time). Using iodine K-edge dual energy imaging, the ion concentration at pore scale with a high temporal and spatial resolution could be quantified. This enabled us to reveal the mechanisms controlling solute transport during saline water evaporation from porous media and extend the corresponding physical understanding of this process. Within this context, the effects of particle size distribution on the dispersion coefficient were investigated together with the evolution of the dispersion coefficient as the evaporation process progresses. The results reported in this dissertation shed new insight on the physics of saline water evaporation from porous media and its complex dynamics. The results of this dissertation have been published in 3 peer-reviewed journal papers together with one additional manuscript which is currently under review.

O presente trabalho teve por objetivo estudar o tempo necessário para que o fluxo da água no solo sob condições de saturação atinja o equilíbrio dinâmico, medindose a condutividade hidráulica saturada (K0) do solo. A condição de equilíbrio dinâmico é a etapa que antecede a aplicação do íon no solo em estudos de eluição de solutos. Utilizaram-se amostras deformadas de três solos com diferentes granulometrias: muito argilosa, média e arenosa, coletadas em área experimental da ESALQ/USP. As amostras foram acondicionadas em colunas de PVC, e utilizou-se água destilada e deaerada para saturação e posterior escoamento nas colunas, conforme a montagem experimental do permeâmetro de carga constante. Os resultados demonstraram que para as amostras dos solos muito argiloso e médio, o equilíbrio foi atingido após 15 dias de drenagem e para o arenoso essa condição só foi verificada após 27 dias, sendo a variabilidade de K0 maior nos primeiros seis dias de ensaio. Concluindo-se, portanto pela inadequação da utilização de um tempo fixo para a adição dos solutos. / The purpose of this work is to study the time for the water flow in soil under saturation to reach the steady-state, by measuring the saturated hydraulic conductivity (K0) of soil. The steady-state is the stage that precedes the application of the ion solution in the soil sample in solute elution studies. Disturbed soil samples with different granulometries were used: very clayey, clay-sandy and sandy soil, collected in experimental area of the ESALQ/USP. The samples were conditioned in PVC columns, and distilled and deaerated water were used for saturation and subsequent drainage in the columns, the test is in accordance with the experimental assembly of constant head permeameter. The data demonstrated that for very clayey and clay-sandy samples, the steady-state was reached after 15 days and for sandy soil this condition was verified after 27 days, being the variability of K0 bigger in the first six days of the test. Therefore it was concluded that is inadequate to use a fixed time for the addition of solutes in soil.