Avaliação do modelo HYDRUS-1D na simulação do transporte de água e potássio em colunas preenchidas com solos tropicais

A preocupação com o destino de produtos químicos e água, aplicados ao solo, tem motivado vários pesquisadores a desenvolverem e aplicarem modelos teóricos, a fim de descrever os processos físicos envolvidos no transporte desses produtos no perfil do solo. Nesse sentido, esta pesquisa teve como objetivo aplicar o modelo HYDRUS-1D, bem como avaliar sua performance, em simulações do deslocamento do potássio e água (umidade do solo), em colunas segmentadas, preenchidas com dois tipos de solos tropicais, em condições não saturadas (Latossolo Vermelho- - Amarelo e Nitossolo Vermelho). Os parâmetros de transporte do potássio foram obtidos por meio de curvas de distribuição de efluentes (Breakthrough Curves (BTC)). O desempenho do modelo foi avaliado com base nos seguintes parâmetros estatísticos: erro máximo, erro absoluto médio, raiz quadrada do erro médio normalizado, coeficiente de massa residual, coeficiente de determinação, eficiência e índice de concordância de Willmott. Diante dos resultados obtidos, pôde-se concluir que o modelo HYDRUS-1D foi eficiente nas simulações de deslocamento do potássio e da água, em relação aos dois materiais de solo estudados.

Transport of iodide in structured clay–loam soil under maize during irrigation experiments analyzed using HYDRUS model

Transport of radioactive iodide 131I− in a structured clay loam soil under maize in a final growing phase was monitored during five consecutive irrigation experiments under ponding. Each time, 27 mm of water were applied. The water of the second experiment was spiked with 200 MBq of 131I− tracer. Its activity was monitored as functions of depth and time with Geiger-Müller (G-M) detectors in 11 vertically installed access tubes. The aim of the study was to widen our current knowledge of water and solute transport in unsaturated soil under different agriculturally cultivated settings. It was supposed that the change in 131I− activity (or counting rate) is proportional to the change in soil water content. Rapid increase followed by a gradual decrease in 131I− activity occurred at all depths and was attributed to preferential flow. The iodide transport through structured soil profile was simulated by the HYDRUS 1D model. The model predicted relatively deep percolation of iodide within a short time, in a good agreement with the observed vertical iodide distribution in soil. We found that the top 30 cm of the soil profile is the most vulnerable layer in terms of water and solute movement, which is the same depth where the root structure of maize can extend.

Optimization of the irrigation time and irrigation frequency by using Hydrus-2D and a capacitance FDR sensor

The evolution of water content on a sandy soil during the sprinkler irrigation campaign, in the summer of 2010, of a field of sugar beet crop located at Valladolid (Spain) is assessed by a capacitive FDR (Frequency Domain Reflectometry) EnviroScan. This field is one of the experimental sites of the Spanish research center for the sugar beet development (AIMCRA). The objective of the work focus on monitoring the soil water content evolution of consecutive irrigations during the second two weeks of July (from the 12th to the 28th). These measurements will be used to simulate water movement by means of Hydrus-2D. The water probe logged water content readings (m3/m3) at 10, 20, 40 and 60 cm depth every 30 minutes. The probe was placed between two rows in one of the typical 12 x 15 m sprinkler irrigation framework. Furthermore, a texture analysis at the soil profile was also conducted. The irrigation frequency in this farm was set by the own personal farmer 0 s criteria that aiming to minimizing electricity pumping costs, used to irrigate at night and during the weekend i.e. longer irrigation frequency than expected. However, the high evapotranspiration rates and the weekly sugar beet water consumption—up to 50mm/week—clearly determined the need for lower this frequency. Moreover, farmer used to irrigate for six or five hours whilst results from the EnviroScan probe showed the soil profile reaching saturation point after the first three hours. It must be noted that AIMCRA provides to his members with a SMS service regarding weekly sugar beet water requirement; from the use of different meteorological stations and evapotranspiration pans, farmers have an idea of the weekly irrigation needs. Nevertheless, it is the farmer 0 s decision to decide how to irrigate. Thus, in order to minimize water stress and pumping costs, a suitable irrigation time and irrigation frequency was modeled with Hydrus-2D. Results for the period above mentioned showed values of water content ranging from 35 and 30 (m3/m3) for the first 10 and 20cm profile depth (two hours after irrigation) to the minimum 14 and 13 (m3/m3) ( two hours before irrigation). For the 40 and 60 cm profile depth, water content moves steadily across the dates: The greater the root activity the greater the water content variation. According to the results in the EnviroScan probe and the modeling in Hydrus-2D, shorter frequencies and irrigation times are suggested.

Modeling water movement in horizontal columns using fractal theory

Fractal mathematics has been used to characterize water and solute transport in porous media and also to characterize and simulate porous media properties. The objective of this study was to evaluate the correlation between the soil infiltration parameters sorptivity (S) and time exponent (n) and the parameters dimension (D) and the Hurst exponent (H). For this purpose, ten horizontal columns with pure (either clay or loam) and heterogeneous porous media (clay and loam distributed in layers in the column) were simulated following the distribution of a deterministic Cantor Bar with fractal dimension H" 0.63. Horizontal water infiltration experiments were then simulated using Hydrus 2D software. The sorptivity (S) and time exponent (n) parameters of the Philip equation were estimated for each simulation, using the nonlinear regression procedure of the statistical software package SAS®. Sorptivity increased in the columns with the loam content, which was attributed to the relation of S with the capillary radius. The time exponent estimated by nonlinear regression was found to be less than the traditional value of 0.5. The fractal dimension estimated from the Hurst exponent was 17.5 % lower than the fractal dimension of the Cantor Bar used to generate the columns.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Avaliação de lavagem parcial de sais em neossolo flúvico irrigado, utilizando modelagem computacional

Este trabalho foi realizado em um vale aluvial irrigado do semiárido do Brasil (Estado de Pernambuco). O objetivo foi investigar o efeito de lâminas de lixiviação e de precipitação no controle da salinidade do solo em lisímetros de drenagem. Com base em simulações numéricas, utilizando-se do modelo HYDRUS-1D, avaliaram-se o fluxo e o transporte de sais na zona não saturada. Mediante medidas de campo de evapotranspiração real e precipitação, foi possível inferir que o modelo se apresentou como ferramenta adequada nas simulações realizadas, para diferentes cenários de lâminas de lixiviação. Verificou-se, experimental e numericamente, que a lixiviação parcial pode ser utilizada como alternativa de manejo para a redução da salinidade na zona radicular, desde que precipitações efetivas possam complementar o processo de lavagem.

Monitoring and modelling the three-dimensional flow of water under drip irrigation

The interpretation of soil water dynamics under drip irrigation systems is relevant for crop production as well as on water use and management. In this study a three-dimensional representation of the flow of water under drip irrigation is presented. The work includes analysis of the water balance at point scale as well as area-average, exploring uncertainties in water balance estimations depending on the number of locations sampled. The water flow was monitored by detailed profile water content measurements before irrigation, after irrigation and 24 h later with a dense array of soil moisture access tubes radially distributed around selected drippers. The objective was to develop a methodology that could be used on selected occasions to obtain 'snap shots' of the detailed three-dimensional patterns of soil moisture. Such patterns are likely to be very complex, as spatial variability will be induced for a number of reasons, such as strong horizontal gradients in soil moisture, variations between individual sources in the amount of water applied and spatial variability is soil hydraulic properties. Results are compared with a widely used numerical model, Hydrus-2D. The observed dynamic of the water content distribution is in good agreement with model simulations, although some discrepancies concerning the horizontal distribution of the irrigation bulb are noted due to soil heterogeneity. (c) 2006 Elsevier B.V. All rights reserved.

Effects of litter manipulation in a tropical Eucalyptus plantation on leaching of mineral nutrients, dissolved organic nitrogen and dissolved organic carbon

Although many studies have shown that soil solution chemistry can be a reliable indicator of biogeochemical cycling in forest ecosystems, the effects of litter manipulations on the fluxes of dissolved elements in gravitational soil solutions have rarely been investigated. We estimated the fluxes of NH4-N, NO3-N, K, Ca, Mg, Na, Cl, dissolved organic nitrogen (DON) and dissolved organic carbon (DOC) over the first two years after re-planting Eucalyptus trees in the coastal area of Congo. Two treatments were replicated in two blocks after clear-cutting 7-year-old stands: in treatment R, all the litter above the mineral soil was removed before planting, and in a double slash (DS) treatment, the amount of harvest residues was doubled. The soil solutions were sampled down to a depth of 4 m and the water fluxes were estimated using the Hydrus 1D model parameterized from soil moisture measurements in 4 plots. Isotopic and spectroscopic analytical techniques were used to assess the changes in dissolved organic matter (DOM) properties throughout the transfer in the soil. The first year after planting, the fluxes of NH4-N, K, Ca, Mg, Na, Cl and DOC in the topsoil of the DS treatment were 2-5 times higher than in R, which showed that litter was a major source of dissolved nutrients. Nutrient fluxes in gravitational solutions decreased sharply in the second year after planting, irrespective of the soil depth, as a result of intense nutrient uptake by Eucalyptus trees. Losses of dissolved nutrients were noticeably low in these Eucalyptus plantations despite a low cation exchange capacity, a coarse soil texture and large amounts of harvest residues left on-site at the clear cut in the DS treatment. All together, these results clarified the strong effect of litter manipulation observed on eucalypt growth in Congolese sandy soils. DOM fluxes, as well as changes in delta C-13, C:N and aromaticity of DOM throughout the soil profile showed that the organic compounds produced in the litter layer were mainly consumed by microorganisms or retained in the topsoil. Below a depth of 15 cm, most of the DOC and the DON originated from the first 2 cm of the soil and the exchanges between soil solutions and soil organic matter were low. (C) 2014 Elsevier B.V. All rights reserved.

Methoden zur Analyse, Interpretation und Modellierung des elektrischen Widerstandes von archäologischen Steinstrukturen

In der Archäologie werden elektrische Widerstandsmessungen routinemäßig zur Prospektion von Fundstellen eingesetzt. Die Methode ist kostengünstig, leicht anwendbar und liefert in den meisten Fällen zuverlässige und leicht zu interpretierende Ergebnisse. Dennoch kann die Methode die archäologischen Strukturen in manchen Fällen nur teilweise oder gar nicht abbilden, wenn die bodenphysikalischen und bodenchemischen Eigenschaften des Bodens und der archäologischen Strukturen dies nicht zulassen. Der spezifische elektrische Widerstand wird durch Parameter wie Wassergehalt, Bodenstruktur, Bodenskelett, Bodentextur, Salinität und Bodentemperatur beeinflusst. Manche dieser Parameter, wie z.B. der Wassergehalt und die Bodentemperatur, unterliegen einer saisonalen Veränderung. Die vorliegende Arbeit untersucht den spezifischen elektrischen Widerstand von archäologischen Steinstrukturen und evaluiert die Möglichkeit, auf Grundlage von Geländemessungen und Laboranalysen archäologische Strukturen und Böden als numerische Modelle darzustellen. Dazu wurde eine Kombination von verschiedenen bodenkundlichen, geoarchäologischen und geophysikalischen Methoden verwendet. Um archäologische Strukturen und Bodenprofile als numerische Widerstandsmodelle darstellen zu können, werden Informationen zur Geometrie der Strukturen und ihren elektrischen Widerstandswerten benötigt. Dabei ist die Qualität der Hintergrundinformationen entscheidend für die Genauigkeit des Widerstandsmodells. Die Geometrie der Widerstandsmodelle basiert auf den Ergebnissen von Rammkernsondierungen und archäologische Ausgrabungen. Die an der Ausbildung des elektrischen Widerstands beteiligten Parameter wurden durch die Analyse von Bodenproben gemessen und ermöglichen durch Pedotransfer-Funktion, wie die Rhoades-Formel, die Abschätzung des spezifischen elektrischen Widerstandes des Feinbodens. Um den Einfluss des Bodenskeletts auf den spezifischen elektrischen Widerstand von Bodenprofilen und archäologischen Strukturen zu berechnen, kamen die Perkolationstheorie und die Effective Medium Theory zum Einsatz. Die Genauigkeit und eventuelle Limitierungen der Methoden wurden im Labor durch experimentelle Widerstandsmessungen an ungestörten Bodenproben und synthetischen Materialien überprüft. Die saisonale Veränderung des Wassergehalts im Boden wurde durch numerische Modelle mit der Software HYDRUS simuliert. Die hydraulischen Modelle wurden auf Grundlage der ermittelten bodenkundlichen und archäologischen Stratigraphie erstellt und verwenden die Daten von lokalen Wetterstationen als Eingangsparameter. Durch die Kombination der HYDRUS-Ergebnisse mit den Pedotransfer-Funktionen konnte der Einfluss dieser saisonalen Veränderung auf die Prospektionsergebnisse von elektrischen Widerstandsmethoden berechnet werden. Die Ergebnisse der Modellierungsprozesse wurden mit den Geländemessungen verglichen. Die beste Übereinstimmung zwischen Modellergebnissen und den Prospektionsergebnissen konnte für die Fallstudie bei Katzenbach festgestellt werden. Bei dieser wurden die Modelle auf Grundlage von archäologischen Grabungsergebnissen und detaillierten bodenkundlichen Analysen erstellt. Weitere Fallstudien zeigen, dass elektrische Widerstandsmodelle eingesetzt werden können, um den Einfluss von ungünstigen Prospektionsbedingungen auf die Ergebnisse der elektrischen Widerstandsmessungen abzuschätzen. Diese Informationen unterstützen die Planung und Anwendung der Methoden im Gelände und ermöglichen eine effektivere Interpretation der Prospektionsergebnisse. Die präsentierten Modellierungsansätze benötigen eine weitere Verifizierung durch den Vergleich der Modellierungsergebnisse mit detailliertem geophysikalischem Gelände-Monitoring von archäologischen Fundstellen. Zusätzlich könnten elektrische Widerstandsmessungen an künstlichen Mauerstrukturen unter kontrollierten Bedingungen zur Überprüfung der Modellierungsprozesse genutzt werden.

Tree species and diversity effects on soil water seepage in a tropical plantation

Plant diversity has been shown to influence the water cycle of forest ecosystems by differences in water consumption and the associated effects on groundwater recharge. However, the effects of biodiversity on soil water fluxes remain poorly understood for native tree species plantations in the tropics. Therefore, we estimated soil water fluxes and assessed the effects of tree species and diversity on these fluxes in an experimental native tree species plantation in Sardinilla (Panama). The study was conducted during the wet season 2008 on plots of monocultures and mixtures of three or six tree species. Rainfall and soil water content were measured and evapotranspiration was estimated with the Penman-Monteith equation. Soil water fluxes were estimated using a simple soil water budget model considering water input, output, and soil water and groundwater storage changes and in addition, were simulated using the physically based one-dimensional water flow model Hydrus-1D. In general, the Hydrus simulation did not reflect the observed pressure heads, in that modeled pressure heads were higher compared to measured ones. On the other hand, the results of the water balance equation (WBE) reproduced observed water use patterns well. In monocultures, the downward fluxes through the 200 cm-depth plane were highest below Hura crepitans (6.13 mm day−1) and lowest below Luehea seemannii (5.18 mm day−1). The average seepage rate in monocultures (±SE) was 5.66 ± 0.18 mm day−1, and therefore, significantly higher than below six-species mixtures (5.49 ± 0.04 mm day−1) according to overyielding analyses. The three-species mixtures had an average seepage rate of 5.63 ± 0.12 mm day−1 and their values did not differ significantly from the average values of the corresponding species in monocultures. Seepage rates were driven by the transpiration of the varying biomass among the plots (r = 0.61, p = 0.017). Thus, a mixture of trees with different growth rates resulted in moderate seepage rates compared to monocultures of either fast growing or slow growing tree species. Our results demonstrate that tree-species specific biomass production and tree diversity are important controls of seepage rates in the Sardinilla plantation during the wet season.

Field-scale water transport in unsaturated crystalline rock

Safe disposal of toxic wastes in geologic formations requires minimal water and gas movement in the vicinity of storage areas, Ventilation of repository tunnels or caverns built in solid rock can desaturate the near field up to a distance of meters from the rock surface, even when the surrounding geological formation is saturated and under hydrostatic pressures. A tunnel segment at the Grimsel test site located in the Aare granite of the Bernese Alps (central Switzerland) has been subjected to a resaturation and, subsequently, to a controlled desaturation, Using thermocouple psychrometers (TP) and time domain reflectometry (TDR), the water potentials psi and water contents theta were measured within the unsaturated granodiorite matrix near the tunnel wall at depths between 0 and 160 cm. During the resaturation the water potentials in the first 30 cm from the rock surface changed within weeks from values of less than -1.5 MPa to near saturation. They returned to the negative initial values during desaturation, The dynamics of this saturation-desaturation regime could be monitored very sensitively using the thermocouple psychrometers, The TDR measurements indicated that water contents changed dose to the surface, but at deeper installation depths the observed changes were within the experimental noise. The field-measured data of the desaturation cycle were used to test the predictive capabilities of the hydraulic parameter functions that were derived from the water retention characteristics psi(theta) determined in the laboratory. A depth-invariant saturated hydraulic conductivity k(s) = 3.0 x 10(-11) m s(-1) was estimated from the psi(t) data at all measurement depths, using the one-dimensional, unsaturated water flow and transport model HYDRUS Vogel er al., 1996, For individual measurement depths, the estimated k(s) varied between 9.8 x 10(-12) and 6.1 x 10(-11) The fitted k(s) values fell within the range of previously estimated k(s) for this location and led to a satisfactory description of the data, even though the model did not include transport of water vapor.

Implementation of the Heated Pulsed Theory using actively heated fiber optics: measurements of soil volumetric water content and soil volumetric heat capacity

Existe una creciente necesidad de hacer el mejor uso del agua para regadío. Una alternativa eficiente consiste en la monitorización del contenido volumétrico de agua (θ), utilizando sensores de humedad. A pesar de existir una gran diversidad de sensores y tecnologías disponibles, actualmente ninguna de ellas permite obtener medidas distribuidas en perfiles verticales de un metro y en escalas laterales de 0.1-1,000 m. En este sentido, es necesario buscar tecnologías alternativas que sirvan de puente entre las medidas puntuales y las escalas intermedias. Esta tesis doctoral se basa en el uso de Fibra Óptica (FO) con sistema de medida de temperatura distribuida (DTS), una tecnología alternativa de reciente creación que ha levantado gran expectación en las últimas dos décadas. Específicamente utilizamos el método de fibra calentada, en inglés Actively Heated Fiber Optic (AHFO), en la cual los cables de Fibra Óptica se utilizan como sondas de calor mediante la aplicación de corriente eléctrica a través de la camisa de acero inoxidable, o de un conductor eléctrico simétricamente posicionado, envuelto, alrededor del haz de fibra óptica. El uso de fibra calentada se basa en la utilización de la teoría de los pulsos de calor, en inglés Heated Pulsed Theory (HPP), por la cual el conductor se aproxima a una fuente de calor lineal e infinitesimal que introduce calor en el suelo. Mediante el análisis del tiempo de ocurrencia y magnitud de la respuesta térmica ante un pulso de calor, es posible estimar algunas propiedades específicas del suelo, tales como el contenido de humedad, calor específico (C) y conductividad térmica. Estos parámetros pueden ser estimados utilizando un sensor de temperatura adyacente a la sonda de calor [método simple, en inglés single heated pulsed probes (SHPP)], ó a una distancia radial r [método doble, en inglés dual heated pulsed probes (DHPP)]. Esta tesis doctoral pretende probar la idoneidad de los sistemas de fibra óptica calentada para la aplicación de la teoría clásica de sondas calentadas. Para ello, se desarrollarán dos sistemas FO-DTS. El primero se sitúa en un campo agrícola de La Nava de Arévalo (Ávila, España), en el cual se aplica la teoría SHPP para estimar θ. El segundo sistema se desarrolla en laboratorio y emplea la teoría DHPP para medir tanto θ como C. La teoría SHPP puede ser implementada con fibra óptica calentada para obtener medidas distribuidas de θ, mediante la utilización de sistemas FO-DTS y el uso de curvas de calibración específicas para cada suelo. Sin embargo, la mayoría de aplicaciones AHFO se han desarrollado exclusivamente en laboratorio utilizando medios porosos homogéneos. En esta tesis se utiliza el programa Hydrus 2D/3D para definir tales curvas de calibración. El modelo propuesto es validado en un segmento de cable enterrado en una instalación de fibra óptica y es capaz de predecir la respuesta térmica del suelo en puntos concretos de la instalación una vez que las propiedades físicas y térmicas de éste son definidas. La exactitud de la metodología para predecir θ frente a medidas puntuales tomadas con sensores de humedad comerciales fue de 0.001 a 0.022 m3 m-3 La implementación de la teoría DHPP con AHFO para medir C y θ suponen una oportunidad sin precedentes para aplicaciones medioambientales. En esta tesis se emplean diferentes combinaciones de cables y fuentes emisoras de calor, que se colocan en paralelo y utilizan un rango variado de espaciamientos, todo ello en el laboratorio. La amplitud de la señal y el tiempo de llegada se han observado como funciones del calor específico del suelo. Medidas de C, utilizando esta metodología y ante un rango variado de contenidos de humedad, sugirieron la idoneidad del método, aunque también se observaron importantes errores en contenidos bajos de humedad de hasta un 22%. La mejora del método requerirá otros modelos más precisos que tengan en cuenta el diámetro del cable, así como la posible influencia térmica del mismo. ABSTRACT There is an increasing need to make the most efficient use of water for irrigation. A good approach to make irrigation as efficient as possible is to monitor soil water content (θ) using soil moisture sensors. Although, there is a broad range of different sensors and technologies, currently, none of them can practically and accurately provide vertical and lateral moisture profiles spanning 0-1 m depth and 0.1-1,000 m lateral scales. In this regard, further research to fulfill the intermediate scale and to bridge single-point measurement with the broaden scales is still needed. This dissertation is based on the use of Fiber Optics with Distributed Temperature Sensing (FO-DTS), a novel approach which has been receiving growing interest in the last two decades. Specifically, we employ the so called Actively Heated Fiber Optic (AHFO) method, in which FO cables are employed as heat probe conductors by applying electricity to the stainless steel armoring jacket or an added conductor symmetrically positioned (wrapped) about the FO cable. AHFO is based on the classic Heated Pulsed Theory (HPP) which usually employs a heat probe conductor that approximates to an infinite line heat source which injects heat into the soil. Observation of the timing and magnitude of the thermal response to the energy input provide enough information to derive certain specific soil thermal characteristics such as the soil heat capacity, soil thermal conductivity or soil water content. These parameters can be estimated by capturing the soil thermal response (using a thermal sensor) adjacent to the heat source (the heating and the thermal sources are mounted together in the so called single heated pulsed probe (SHPP)), or separated at a certain distance, r (dual heated pulsed method (DHPP) This dissertation aims to test the feasibility of heated fiber optics to implement the HPP theory. Specifically, we focus on measuring soil water content (θ) and soil heat capacity (C) by employing two types of FO-DTS systems. The first one is located in an agricultural field in La Nava de Arévalo (Ávila, Spain) and employ the SHPP theory to estimate θ. The second one is developed in the laboratory using the procedures described in the DHPP theory, and focuses on estimating both C and θ. The SHPP theory can be implemented with actively heated fiber optics (AHFO) to obtain distributed measurements of soil water content (θ) by using reported soil thermal responses in Distributed Temperature Sensing (DTS) and with a soil-specific calibration relationship. However, most reported AHFO applications have been calibrated under laboratory homogeneous soil conditions, while inexpensive efficient calibration procedures useful in heterogeneous soils are lacking. In this PhD thesis, we employ the Hydrus 2D/3D code to define these soil-specific calibration curves. The model is then validated at a selected FO transect of the DTS installation. The model was able to predict the soil thermal response at specific locations of the fiber optic cable once the surrounding soil hydraulic and thermal properties were known. Results using electromagnetic moisture sensors at the same specific locations demonstrate the feasibility of the model to detect θ within an accuracy of 0.001 to 0.022 m3 m-3. Implementation of the Dual Heated Pulsed Probe (DPHP) theory for measurement of volumetric heat capacity (C) and water content (θ) with Distributed Temperature Sensing (DTS) heated fiber optic (FO) systems presents an unprecedented opportunity for environmental monitoring. We test the method using different combinations of FO cables and heat sources at a range of spacings in a laboratory setting. The amplitude and phase-shift in the heat signal with distance was found to be a function of the soil volumetric heat capacity (referred, here, to as Cs). Estimations of Cs at a range of θ suggest feasibility via responsiveness to the changes in θ (we observed a linear relationship in all FO combinations), though observed bias with decreasing soil water contents (up to 22%) was also reported. Optimization will require further models to account for the finite radius and thermal influence of the FO cables, employed here as “needle probes”. Also, consideration of the range of soil conditions and cable spacing and jacket configurations, suggested here to be valuable subjects of further study and development.

Comparison among monitoring strategies to assess water flow dynamic and soil hydraulic properties in agricultural soils

Irrigated agriculture is usually performed in semi-arid regions despite scarcity of water resources. Therefore, optimal irrigation management by monitoring the soil is essential, and assessing soil hydraulic properties and water flow dynamics is presented as a first measure. For this purpose, the control of volumetric water content, θ, and pressure head, h, is required. This study adopted two types of monitoring strategies in the same experimental plot to control θ and h in the vadose zone: i) non-automatic and more time-consuming; ii) automatic connected to a datalogger. Water flux was modelled with Hydrus-1D using the data collected from both acquisition strategies independently (3820 daily values for the automatic; less than 1000 for the non-automatic). Goodness-of-fit results reported a better adjustment in case of automatic sensors. Both model outputs adequately predicted the general trend of θ and h, but with slight differences in computed annual drainage (711 mm and 774 mm). Soil hydraulic properties were inversely estimated from both data acquisition systems. Major differences were obtained in the saturated volumetric water content, θs, and the n and α van Genuchten model shape parameters. Saturated hydraulic conductivity, Ks, shown lower variability with a coefficient of variation range from 0.13 to 0.24 for the soil layers defined. Soil hydraulic properties were better assessed through automatic data acquisition as data variability was lower and accuracy was higher.

Fate and Transport of Naproxen in a Sandy Aquifer Material: Saturated Column Studies and Model Evaluation

Naproxen-C14H14O3 is a nonsteroidal anti-inflammatory drug which has been found at detectable concentrations in wastewater, surface water, and groundwater. Naproxen is relatively hydrophilic and is in anionic form at pH between 6 and 8. In this study, column experiments were performed using an unconsolidated aquifer material from an area near Barcelona (Spain) to assess transport and reaction mechanisms of Naproxen in the aquifer matrix under different pore water fluxes. Results were evaluated using HYDRUS-1D, which was used to estimate transport parameters. Batch sorption isotherms for Naproxen conformed with the linear model with a sorption coefficient of 0.42 (cm3 g−1), suggesting a low sorption affinity. Naproxen breakthrough curves (BTCs) measured in soil columns under steady-state, saturated water flow conditions displayed similar behavior, with no apparent hysteresis in sorption or dependence of retardation (R, 3.85-4.24) on pore water velocities. Soil sorption did not show any significant decrease for increasing flow rates, as observed from Naproxen recovery in the effluent. Sorption parameters estimated by the model suggest that Naproxen has a low sorption affinity to aquifer matrix. Most sorption of Naproxen occurred on the instantaneous sorption sites, with the kinetic sorption sites representing only about 10 to 40% of total sorption.

Artificial pumping errors in the Kool-Parker scaling model of soil moisture hysteresis

Hysteresis models that eliminate the artificial pumping errors associated with the Kool-Parker (KP) soil moisture hysteresis model, such as the Parker-Lenhard (PL) method, can be computationally demanding in unsaturated transport models since they need to retain the wetting-drying history of the system. The pumping errors in these models need to be eliminated for correct simulation of cyclical systems (e.g. transport above a tidally forced watertable, infiltration and redistribution under periodic irrigation) if the soils exhibit significant hysteresis. A modification is made here to the PL method that allows it to be more readily applied to numerical models by eliminating the need to store a large number of soil moisture reversal points. The modified-PL method largely eliminates any artificial pumping error and so essentially retains the accuracy of the original PL approach. The modified-PL method is implemented in HYDRUS-1D (version 2.0), which is then used to simulate cyclic capillary fringe dynamics to show the influence of removing artificial pumping errors and to demonstrate the ease of implementation. Artificial pumping errors are shown to be significant for the soils and system characteristics used here in numerical experiments of transport above a fluctuating watertable. (c) 2005 Elsevier B.V. All rights reserved.