Dissertations / Theses on the topic 'HYDRUS-1D'

Tools for predicting rain garden performance are limited, particularly for predicting seasonal performance. Changes in temperature cause changes in the viscosity of water, infiltration rates, and evapotranspiration rates. A variably-saturated soil model, HYDRUS-1D, was calibrated and validated using observed ponding depth and soil moisture data for the Philadelphia Zoo Rain Garden, owned and operated by the Philadelphia Water Department (PWD). Warm and cold seasons were simulated with typical meteorological data and temperature-adjusted saturated hydraulic conductivity values. Design-storm simulations confirmed that the rain garden is over-performing. Maximum capacity of the system was simulated by increasing the loading ratio until overtopping occurred or ponding remained longer than 24 hours. This study will demonstrate how modeling of an operational urban rain garden offers a realistic picture of performance and could be used as a tool for informing regulations and design.

Modelos matemáticos são comumente utilizados no estudo da dinâmica da água no solo não-saturado. A principal equação para se quantificar esta dinâmica é a equação diferencial de Richards. Sua solução direta é impossível na maioria dos casos, necessitando de métodos numéricos, dos quais se destaca a utilização das funções de condutividade e de retenção de água de Van Genuchten - Mualem para obtenção da solução numérica. Diante disto, o objetivo deste trabalho foi apresentar um novo método para a obtenção dos parâmetros de Van Genuchten - Mualem, que utiliza a modelagem inversa de dados de teores de água medidos periodicamente e obtidos por experimento de evaporação assistido por atenuação de radiação gama, para a modelagem inversa foi utilizado o software Hydrus-1D. O método foi testado em amostras com diferentes texturas, colhidas em 11 localidades na região de Piracicaba-SP. Dos conjuntos de exemplares analisados, apenas um apresentou resultado insatisfatório, concluindo que o método é válido
Mathematical models are commonly used in studies of water dynamics in unsaturated soil. The main equation to quantify water the dynamics is the differential Richards equation. Its analytical solution is impossible in almost all cases, requiring numerical methods, among which the Van Genuchten - Mualem water conductivity and water retention functions are frequently used to obtain the numerical solution. The objective of this work is to present a new method to obtain the Van Genuchten - Mualem parameters, using the inverse modeling of water content data measured periodically by gamma radiation attenuation in evaporating samples. Hydrus-1D software was used for the inverse modeling. The method was tested in samples with different textures, collected in 11 locations in the region of Piracicaba, state of São Paulo, Brazil. Of the sets of samples analyzed, only one presented an unsatisfactory result, concluding that the method is valid

The percolation of water through waste rock dumps at mine sites can lead to the production of a leachate with high concentrations of dissolved metals, sulfate and nitrogen compounds. It is important to understand how water flows in waste rock dumps in order to predict the environmental impact of this leachate on recipients. The dynamics of percolation and leachate discharge are controlled by climatological conditions at the site, where relatively large flows in northern Sweden correspond to snowmelt during late Spring. Rock dumps are often tens of meters in height, resulting in an unsaturated water flow system through heterogeneous material. Hence, the simulation of leachate generation requires an accurate representation of the subsurface materials as well as the flow processes, where water flow in waste rock dumps is dominated by matrix flow with macropore flow being of secondary importance. Matrix flow is rather slow and may thus potentially yield relatively high concentrations of contaminants in the leachate, in response to precipitation and snow melt. This study uses Hydrus-1D to predict leachate generation from a small-scale waste rock dump in Kiruna in terms of discharge magnitude and timing. The 3-dimensional geometry of the waste rock dump is approximated by summing simulations from 1225 one-dimensional columns of different length, with a surface area of 1 m2 each. There are four output parameters that are compared between the model results and measured data: snow accumulation, water content, temperature, and discharge. There are some discrepancies between the model results and field measurements, most likely due to uncertainties in the input parameters (especially waste rock properties), limitations in the Hydrus-1D model (i.e. freeze-thaw dynamics), and assumptions that are used in constructing the conceptual model. For better agreement between model results and measured data, a new modelling approach is recommended, potentially using a different program than Hydrus-1D.

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, objetivando descrever os processos físicos envolvidos no transporte desses produtos no perfil do solo. Nesse sentido, a presente pesquisa teve como objetivo a aplicação do modelo Hydrus-1D para simulação do movimento da água e dos íons nitrato e potássio, em condições de laboratório, utilizando-se para tal colunas segmentadas preenchidas com dois tipos de solo não-saturado (Latossolo Vermelho Amarelo e Nitossolo Vermelho), bem como avaliar o desempenho e sensibilidade do referido modelo. Para a obtenção dos parâmetros de transporte de cada soluto, foram elaboradas as curvas de distribuição de efluentes (Breakthrough Curves BTC). Os parâmetros de transporte obtidos e as condições de contorno de cada ensaio foram inseridos no modelo Hydrus- 1D para realização das simulações. O desempenho do modelo foi avaliado com base nos 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. A sensibilidade do modelo foi avaliada conforme o método proposto por McCuen e Snyder (1986) e a análise foi aplicada aos parâmetros: fluxo de entrada, condutividade hidráulica do solo saturado, teor de água na saturação, alfa e n (parâmetros de ajuste da curva de retenção), coeficiente de distribuição e dispersividade. Os resultados experimentais mostraram que o deslocamento do íon nitrato acompanhou a frente de molhamento e em relação ao potássio, observou-se uma maior retenção nas camadas superficiais da coluna de solo, para ambos os solos. Pôde-se concluir que o modelo Hydrus-1D foi eficiente para simulações de deslocamento de potássio e água para ambos os solos estudados e mediante a utilização do modelo de equilíbrio para o transporte de solutos no solo, o modelo Hydrus-1D não foi eficiente para simular o deslocamento de nitrato, para ambos os solos. Os procedimentos experimentais para estimativa de parâmetros de transporte, como as curvas de distribuição de efluentes, foram suficientes para descrever a movimentação de potássio no solo, gerando informações de entrada precisas para os modelos de simulação. Para a simulação do deslocamento do íon nitrato, concluiu-se que o modelo Hydrus-1D requer um processo experimental mais detalhado, sendo necessária a estimativa de um maior número de parâmetros de transporte. Em relação à análise de sensibilidade, ao simular os teores de potássio e água no solo, o modelo apresentou maior sensibilidade aos parâmetros teor de água na saturação e fluxo de entrada da solução. Tais parâmetros, portanto, precisam ser determinados com maior precisão. Houve baixa sensibilidade aos parâmetros condutividade hidráulica do solo saturado e dispersividade, para ambos os solos estudados.
The concern about the fate of chemical products and water, applied to the soil, has been motivating several researchers to develop and apply theoretical models, aiming to describe the physical processes involved in the transport of those products in soil profile. The present research had as objective the application of the model Hydrus- 1D for water and solute (nitrate and potassium) simulation profile, in laboratory conditions, using soil columns filled with two types of unsaturated soil (Oxisol (Haplustox) and Hapludox), as well to evaluate the acting and sensitivity of the model. The transport parameters were obtained for each solute by breakthrough curves (BTCs). The transport parameters obtained and the boundary conditions were inserted in the Hydrus-1D model to realize the simulations. The acting of the Hydrus-1D model was evaluated using the statistical indicators: maximum error, mean absolute error, normalized root mean-square error, coefficient of residual mass, determination coefficient, efficiency and Willmott concordance index. The model sensitivity was evaluated by the method proposed by McCuen and Snyder (1986) and applied to the parameters: input flow, soil saturated hydraulic conductivity, water content (saturation point), alpha and n (soil water retention curve parameters), distribution coefficient and dispersivity. The obtained results, experimentally, showed the nitrate displacement following the wetting front and, in relation to the potassium ion, a larger retention was observed at the superficial soil columns layers. Therefore, was possible to conclude that Hydrus-1D model was efficient for both water and potassium displacement simulations, for both studied soils and by the balance model for solute transport in soil, the Hydrus- 1D model was not efficient to simulate the nitrate displacement in both soils. The experimental procedures to estimate transport parameters, by the BTCs, were enough to describe the potassium movement in soil, generating input information necessary to the model simulation. To the simulation of the nitrate displacement, was possible to conclude that the model Hydrus-1D requests a more detailed experimental process, being necessary the estimate of a larger number of transport parameters. In relation to the sensitivity analysis of the Hydrus-1D model, when simulated the content of potassium and water at the soil columns, was observed that the model shows more sensitivity about to the parameters: water content (saturation point) and input flow. This parameters need to be estimated with more precision. There was a low sensitivity to the parameters soil saturated hydraulic conductivity and dispersivity for both studied soils.

Groundwater is a crucial part of the hydrological cycle and is an important source for drinking water, irrigation and industry, particularly during droughts. With climate change, the hydrological variability is predicted to increase, making predictions for recharge and groundwater storage even more important to implement and to maintain sustainable water use. This study examines the importance of model spin-up in simulating recharge using the Hydrus-1D computer model. The focus is on two previously made Hydrus-1D models that represent end members in climate and hydrology; one which is a natural grassland in a semi-arid climate, while the other is a low impact development (LID) bioswale site in a Mediterranean climate. The main goal of this study is to characterize the range and causes of spin-up behavior as well as to analyze the extent of the effects that the spin-up process has on the recharge simulations. Although there has been some research on spin-up behavior for surface-water models, there is still a knowledge gap regarding the effects of model spin-up on vadose zone models simulating recharge. The initial conditions varied using three parameters for each of the two models: time (3, 15 and 30 years), initial moisture (θ = 0.1, 0.2 and 0.3) and precipitation (25% drier than historical data, historical 30-year data, 25% wetter than historical data). The output from these spin-ups were then used as initial conditions in simulating recharge using the 15-year models. The study found that the impact of spin-up is significant in the natural grassland site where there is a slow response between atmospheric forcings and recharge and where there is a relatively thick vadose zone. Especially spin-up time showed great variability and there is an inverse relationship between spin-up time and magnitude of recharge, where the longer spin-ups had lower recharge rates. Initial water content and precipitation did not result in different recharge amounts for the LID model. Length of spin-up only had very small differences in recharge for the LID models, indicating they are less sensitive to changes in initial spin-up parameters.

A demanda alimentar vem aumentando de acordo com o crescimento populacional e com isso há uma necessidade de que as práticas agrícolas tornem-se cada vez mais intensas e tecnificadas, mediante o incremento, entre outros, da utilização de insumos químicos. Porém, uma vez aplicados de maneira desorganizada, sem a devida preocupação com os possíveis impactos aos recursos naturais, tais produtos podem vir a provocar danos ao solo, contaminação de águas subterrâneas (mediante sua lixiviação) e em alguns casos, podem ser responsáveis pelo aumento da salinidade em alguns solos da região do semiárido do nordeste brasileiro. Nesse sentido, devido ao fato de envolverem processos físicos na natureza, o estudo da dinâmica desses produtos químicos têm motivado diversos pesquisadores a aplicarem ferramentas matemáticas (modelos matemáticos) na área de engenharia de água e solo, buscando entender a correlação entre a dinâmica da água e o movimento de solutos no perfil do solo. Portanto, a presente pesquisa tem como objetivo a aplicação do modelo matemático HYDRUS-1D para a simulação do movimento da água e do íon potássio, em condições de laboratório, utilizando-se colunas preenchidas com dois tipos de material de solos salinos e não saturados: Argissolo Amarelo (S1) e Latossolo Vermelho Amarelo (S2). Além disso, buscou-se também avaliar tanto o desempenho do referido modelo, em condições salinas, quanto proceder a uma análise de sensibilidade. Para alcançar tais objetivos, foram conduzidas as seguintes etapas: 1) elaboração de curvas de distribuição de efluentes (Breakthrough Curves - BTC) para obtenção dos parâmetros de transporte do íon potássio, com a aplicação de soluções de 1000 ppm, 2000 ppm e 3000 ppm de concentração de potássio, nos dois tipos de solos salinos, 2) Aplicação dessas soluções de potássio na coluna preenchida com o mesmo material de solo salino não saturado, 3) Simulação dos perfis de água e potássio para os diferentes materiais de solo salino pelo modelo HYDRUS-1D e 4) Avaliação de desempenho e análise de sensibilidade do modelo HYDRUS-1D. A avaliação de desempenho foi feita levando-se em conta os seguintes índices 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 e a análise de sensibilidade foi feita mediante a determinação do erro padrão, por meio de variações positivas (+ 10%, + 20%, + 30%, + 40%, + 50%) e negativas (- 10%, - 20%, - 10%, - 40%, - 50%), dos valores de: umidade volumétrica do solo saturado, condutividade hidráulica e os parâmetros n e ? do modelo de van Genuchten (1980). Diante dos resultados obtidos, pôde-se perceber que os procedimentos experimentais (BTC\'s) para a estimativa dos parâmetros de transporte do íon potássio em condições de solo salino, foram suficientes para descrever a movimentação do íon potássio no solo, gerando informações de entrada precisas para os modelos de simulação. Em termos do modelo HYDRUS-1D, observou-se que o deslocamento do íon potássio acompanhou a frente de molhamento e que o modelo foi eficiente nas simulações de deslocamento do íon potássio e da água para ambos os solos, em condições salinas. A análise de sensibilidade evidenciou que o modelo apresentou-se sensível, às variações negativas dos dados de entrada: umidade volumétrica do solo saturado, condutividade hidráulica do solo saturado e parâmetros \"n\" e \"?\" do modelo de van Genuchten (1980).
The increased demand for food due to population growth requires that agricultural practices become increasingly intensive and very technical, including the increased use of agricultural chemicals (fertilizers). If improperly applied without considering possible impacts on natural resources, agricultural chemicals may lead to soil and groundwater contamination through their leaching from the soil root zone. They may be responsible also for increased salinity in some soils of semiarid regions in northeastern Brazil. Since many transient physical and chemical processes affect their transport in the subsurface, mathematical models have become popular tools in soil and water engineering and management in order to understand the correlation between water dynamics and solute movement in soils. Thus, this research aimed at using the HYDRUS-1D software package to simulate water and potassium movement, under laboratory conditions, in unsaturated saline soil columns filled with two soil types: an Ultisol (S1) and an Oxisol (S2). Comparisons were made with experimental data while also a sensitivity analysis was carried out to evaluate the effect of various parameters on solute transport under saline conditions. For this purpose the following studies were performed: 1) Measurement of solute breakthrough curves (BTCs) to estimate the transport parameters of the potassium ion by applying potassium solutions of 1,000 ppm, 2,000 ppm and 3,000 ppm to both soil types, 2) Application of similar potassium solutions to columns containing the same saline unsaturated soils, 3) Simulation of water and potassium distributions for different saline soil materials using HYDRUS-1D, and 4) Performance evaluation and sensitivity analyses of the HYDRUS-1D numerical model. The performance evaluation was conducted using the following statistical indices: maximum error, mean absolute error, normalized root mean square error, coefficient of residual mass, coefficient of determination, efficiency and Willmott\'s concordance index. The sensitivity analyses considered standard deviations resulting from positive and negative changes (+ 10% + 20% + 30% + 40% + 50%) (-10%, -20%, -10%, -40%, -50%) in the values of the saturated volumetric soil moisture content, the saturated hydraulic conductivity, and the parameters n and ? of van Genuchten\'s (1980) model for the unsaturated soil hydraulic functions. Results indicate that the experimental procedures (BTCs) for estimating the transport parameters of potassium for saline soil conditions were sufficient to describe potassium ion transport in the soils by generating the required input information for the simulation models. Relative to HYDRUS-1D model, the displacement of potassium was found to follow closely the wetting front, with the model providing a very efficient means for simulating the movement of both water and potassium in the two soils during saline conditions. The sensitivity analysis showed that the model was relatively sensitive to negative variations of the input data, notably the saturated water content, the saturated hydraulic conductivity, and the n and ? soil hydraulic parameters.

Made available in DSpace on 2014-06-12T17:42:16Z (GMT). No. of bitstreams: 2 arquivo7111_1.pdf: 1388033 bytes, checksum: 645508c7ef101ae31fed17e866992644 (MD5) license.txt: 1748 bytes, checksum: 8a4605be74aa9ea9d79846c1fba20a33 (MD5) Previous issue date: 2011
Faculdade de Amparo à Ciência e Tecnologia do Estado de Pernambuco
Na atualidade, metade da população mundial reside nos grandes centros urbanos. Grande parcela desse crescimento tem ocorrido em países em desenvolvimento. No Brasil, já se verifica um contingente de cerca de 80% da população residindo nas áreas urbanas. Neste crescente cenário de urbanização, impactos ambientais e socioeconômicos decorrentes de eventos hidrológicos têm sido recorrentes, afetando grande parte da população. O aumento da impermeabilização reduz as taxas de infiltração, que por sua vez leva à diminuição das taxas de recarga para os aquíferos e à diminuição do escoamento de base. O escoamento superficial é intensificado, aumentando em velocidade e, a frequência e magnitude dos picos de cheia, levando ocasionalmente às enchentes. Nesse contexto, têm sido empregados os telhados verdes em várias partes do mundo principalmente com finalidades estéticas de valorização do espaço urbano e para melhoria do conforto ambiental. Essas áreas verdes podem servir também para detenção do escoamento superficial, minimizando as enchentes urbanas. A simulação da dinâmica da água no solo do telhado verde realizada no programa Hydrus 1-D, a partir das características do sistema do teto verde implantado e dados obtidos em campo, proporcionou a caracterização da dinâmica da água em seu perfil de solo, fornecendo subsídios quanto ao desempenho deste dispositivo no amortecimento do escoamento superficial oriundo dos telhados

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior
Soils subjected to successive pig slurry (PS) amendments present zinc (Zn) and copper (Cu) increases, mostly on soil surface layer. The continuous application of PS causes accumulation of these heavy metals (HM) in soil, increasing risks for environmental contamination. Therefore, development of methods to estimate Zn and Cu accumulation in surface layer and their movement in soil is very important for the management planning of PS and of the agricultural soils amended with these residues. The main objective of this work was to model the vertical transport of Zn and Cu in profile of agricultural soils receiving successive application of PS at high, moderate and low doses, evaluating the effects of soil tillage and water absorption by roots on the movement of these HM. We also aimed to predict the long-term accumulation of HM in soil surface and their transport into groundwater, comparing estimated results with the threshold values established for agricultural soils and groundwater by the Conselho Nacional do Meio Ambiente (CONAMA). Two studies were performed on field experiments composed by different PS doses, the first located on an Alfisol in Santa Maria RS, and the second on an Oxisol in Campos Novos SC. Hydro-physical and chemical characteristics of these soils were measured and, accompanied by atmospheric and PS application data sets, were introduced into HYDRUS-1D for parameterization of the two-site model used. In the first study, simulations for next 100 years were performed under scenarios combining two PS doses and four time intervals between soil tillage of the arable layer. The solute transport model used on these simulations was already validated for this local soil in another study. In the second study, the validation of a solute transport model was carried out primarily, following the same approaches used for the Alfisol. Afterwards, root water uptake and root growth modules were introduced to the model. The future scenarios were simulated using applications of different PS doses during the next 50 years. The solute transport model previously validated for the Alfisol was also validated for the Oxisol, and the introduction of the root modules in HYDRUS-1D produced even better results. This increased the usefulness of the model for its use in simulations related to Zn and Cu transport on other contaminated soils. Future scenarios exhibit a great accumulation of these two HM in soil surface layer during the simulated times. Moreover, when soils are submitted to high PS annual doses, its Cu concentrations reached the maximum values established by CONAMA for agricultural soils in approximately 94 years on the Alfisol and in 29 years on the Oxisol. Reduction in PS doses and sporadic soil tillage are strategies that slow the increase of these HM concentrations in soil surface layer, decreasing so the environmental risks and making PS additions on agricultural soils viable for longer time.
Solos submetidos a sucessivas aplicações de dejeto líquido de suínos (DLS) apresentam incrementos nos seus teores de zinco (Zn) e cobre (Cu), principalmente na camada superficial. A continuidade dessa prática faz com que esses metais pesados (MP) se acumulem no solo, aumentando o potencial de contaminação do ambiente. Dessa forma, o desenvolvimento de métodos que estimem o acúmulo superficial e o movimento de Zn e Cu no solo é de extrema importância no planejamento do manejo dos DLS e dos solos agrícolas que recebem esses resíduos. O trabalho objetivou modelar o transporte vertical de Zn e Cu no perfil de solos agrícolas contaminados por aplicações sucessivas de altas, moderadas e baixas doses de DLS, avaliando o efeito do revolvimento do solo e da absorção de água pelas raízes no movimento desses MP, prevendo, no longo prazo, seus acúmulos na camada superficial do solo e suas transferências rumo ao lençol freático, confrontando os resultados estimados com os respectivos limites de concentração para solos agrícolas e águas subterrâneas estabelecidos pelo Conselho Nacional do Meio Ambiente (CONAMA). Foram realizados dois estudos em experimentos de campo com diferentes doses de DLS, o primeiro localizado sobre um Argissolo, em Santa Maria RS, e o segundo sobre um Latossolo, em Campos Novos SC. Diversas características físico-hídricas e químicas foram determinadas nesses solos para, juntamente com os dados atmosféricos e do histórico de aplicação de DLS nessas áreas, serem introduzidos no HYDRUS-1D para parametrização do modelo químico de dois sítios utilizado. No primeiro estudo foram realizadas simulações para os próximos 100 anos, considerando cenários combinando duas doses de DLS e quatro intervalos de tempo entre revolvimentos periódicos da camada arável do solo. O modelo de transporte de solutos utilizado nessas simulações já havia sido validado para o solo deste local em outro trabalho. No segundo estudo foi feita a validação de um modelo de transporte de solutos, seguindo as mesmas premissas do modelo utilizado no Argissolo. Posteriormente, módulos de absorção de água e crescimento de raízes foram incorporados ao modelo. Realizaram-se então simulações futuras de cenários considerando aplicações de diferentes doses de DLS durante os próximos 50 anos. O modelo de transporte de solutos validado previamente no Argissolo foi validado também no Latossolo, apresentando resultados ainda melhores quando os módulos de raízes foram introduzidos no HYDRUS-1D. Isso aumenta a abrangência desse modelo para simular o transporte de Zn e Cu em outros solos contaminados. Os cenários futuros mostram que ocorre grande acúmulo desses MP na camada superficial do solo durante os períodos simulados. Além disso, quando os solos são submetidos a altas doses anuais de DLS, as suas concentrações de Cu atingem os valores máximos estipulados pelo CONAMA, para solos agrícolas, em aproximadamente 94 anos no Argissolo e 29 anos no Latossolo. Reduzir as doses de DLS e revolver o solo esporadicamente são alternativas que retardam o aumento da concentração desses MP na camada superficial do solo, diminuindo o potencial de impacto ao ambiente e tornando a adição de DLS em áreas agrícolas viável por mais tempo.

This research investigated the water use of different plant types in urban landscapes, nitrogen (N) and water transport in turf, and potential N leaching from urban landscapes to ground water. In the first study, three landscape treatments integrating different types of plants—woody, herbaceous perennial, turf—and putative water use classifications—Mesic, Mixed, Xeric—were grown in large drainage lysimeters. Each landscape plot was divided into woody, turf, and herbaceous perennial plant hydrozones and irrigated for optimum water status over two years, with water use measured using a water balance approach. For woody plants and herbaceous perennials, canopy cover, rather than plant type or water use classification, was the key determinant of water use relative to reference evapotranspiration (ETo) under well-watered conditions. For turf, monthly evapotranspiration (ETa) followed a trend linearly related to ETo. In the second study, water transport parameters were calibrated using an inverse simulation with Kentucky bluegrass (KBG). Subsequently, those parameters were applied to simulate water use by tall fescue (TF) and buffalograss (BG) turfgrasses using numerical modeling (Hydrus-1D). By using the calibrated soil hydraulic parameters obtained from the water transport simulation, N transport and transformation was modeled with Hydrus- 1D under different irrigation rates and different fertilization rates. Different soil texture scenarios were also simulated to demonstrate the influence of soil texture on N leaching. In the third study, the simulated N-leaching from different soil textures was integrated into a Geographic Information System (GIS) approach to estimate NO3-N leaching mass from urban turf areas. Nitrate-N leaching risks to ground water under overirrigation and overfertilization scenarios and efficient irrigation and fertilization scenarios were estimated. The results showed improvement of turf irrigation and fertilization management may decrease N-leaching significantly and greatly decrease the risk of ground water being contaminated by NO3-N leaching in the Salt Lake Valley.

Water is undeniably among the most important natural resources and the most critical in semi-arid regions like the Intermountain West of the United States. Such regions are characterized by low precipitation, the majority of which is transferred to the atmosphere from the soil and vegetation as evapotranspiration (ET). Quantification of ET is thus crucial for understanding the balance of water within the region, which is important for efficiently planning the available water resources. This study was motivated towards advancing the estimation of actual ET (ETA) in mountain ecosystems, where the variation in different types of vegetation and non-uniformity of soil including considerable stone content creates challenges for estimating water use as ET. With the aim of addressing the effect of stone content in controlling soil moisture and ET, this study examined the influence of stone content on bulk soil hydraulic properties. An averaging model referred to as a binary mixing model was used to describe the way in which water is held and released in stony soil. This approach was based on the individual hydraulic behavior of the background soil and of the stones within the soil. The effect of soil stone content on ETA was evaluated by accounting for the water retention properties of stones in the soil using a numerical simulation model (HYDRUS-1D). The results revealed overestimation of simulated ETA when effects of stone content were not accounted for in comparison to ETA measured by the state-of-the-art “eddy covariance” measurement method for ETA. An even larger-scale model was evaluated, named the Noah-Multiphysics (Noah-MP) land surface model. The land surface model was run using different arrangements of complexity to determine the importance of stone content information on simulation results. The version of the model with information about stone content along with detailed soil properties was able to provide the best Noah-MP prediction of ET. The study suggests that improvement in representation of soil properties including stone content information, can substantially advance the ability of numerical and land surface models to more accurately simulate soil water flow and ETA.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES
This research was conducted with the objective of studying the dynamics of soil temperature and water content as a function of the different amounts of black oat residue (Avena sativa Schreb.) in the soil cover, through existing relationships with meteorological variables, as well as the simulation of these dynamics using Hydrus-1D. The experiment was carried out in an Argissolo Vermelho Distrófico arênico (Rhodic Paleudalf) in Santa Maria - RS. The treatments consisted of the use of 3, 6 and 9 Mg ha-1 of oat straw in soil cover and an uncovered soil treatment (0 Mg ha-1). Soil temperature and water content at different depths, in addition to meteorological conditions, were measured during the period from December 1, 2014 to December 6, 2015. The effect of straw on soil temperature was analyzed at different depths, correlating covered soil temperature with uncovered soil temperature. The equations used for estimating soil surface temperature were adjusted by the least square method, where the surface temperature is estimated as a function of air temperature and incident solar radiation, and its performance compared to the sinusoidal equation. The hydraulic and thermal parameters required for the simulation of soil temperature and water content were optimized by the reverse solution and validated on uncovered soils. The optimized parameters were used to simulate the temperature and water content in soils with different amounts of straw with two contour conditions for the soil surface temperature, one measure and the other estimated. The straw affected the soil temperature down to 50 cm deep. The uncovered soil temperature at which inversion of the straw effect occurred was 19.6 ° C at the surface and 15.8 ° C at 50 cm depth. The presence of straw keeps the soil warmer than the soil without straw at the beginning of the day, but prevents it from warming similar to the ground without straw for the rest of the day, causing in turn that, bare soil gets warmer during the day and cools more at night compared to the covered soil. This temperature was similar in all amounts of straw, but decreased in depth. The temperature of the soil surface can be estimated as a function of air temperature and solar radiation. In uncovered soil the proposed equation presents better performance compared to the sinusoidal equation. In covered soils the sinusoidal equation had better performance. The optimal hydraulic and thermal parameters of the soil improved the simulations of temperature and soil water content. The soil temperature with different amounts of oat straw can be simulated using Hydrus-1D, and the use of the contour condition with the measured surface temperature performed better than the estimated surface temperature.
Esta pesquisa foi conduzida com o objetivo de estudar a dinâmica da temperatura e do conteúdo de água do solo em função de diferentes quantidades de resíduos de aveia preta disposta em cobertura por meio das relações existente com variáveis meteorológicas, bem como a representação destes comportamentos por meio de simulação utilizando o Hydrus-1D. O experimento foi conduzido em um Argissolo Vermelho Distrófico arênico em Santa Maria - RS. Os tratamentos consistiram da utilização de três diferentes níveis de cobertura do solo (3, 6 e 9 Mg ha-1) com palha de aveia (Avena sativa) e um tratamento com solo descoberto (0 Mg ha-1). A temperatura do solo e o conteúdo volumétrico de água, em diferentes profundidades, além das condições meteorológicas, foram mensurados durante o período de 01/12/2014 a 06/12/2015. O efeito da cobertura do solo com resíduos vegetais na temperatura do solo foi analisado nas diferentes profundidades correlacionando a temperatura do solo com palha e a temperatura do solo descoberto. O método dos mínimos quadrados foi usado para gerar equações de estimativa da temperatura da superfície do solo em função da temperatura do ar e da radiação solar e o seu desempenho foi comparado com a equação senoidal. Os parâmetros hidráulicos e térmicos, necessários para a simulação da temperatura e do conteúdo de água do solo, com o Hydrus-1D, foram otimizados por meio da solução inversa e validados em solos descobertos. Os parâmetros otimizados foram usados para a simulação da temperatura e do conteúdo de água em solos com diferentes quantidades de palha com duas condições de contorno para a temperatura da superfície do solo, uma medida e outra estimada. A palha exerce influência sobre a temperatura do solo até 50 cm de profundidade. A temperatura do solo descoberto em que ocorre a inversão do efeito da palha foi de 19,6 °C na superfície e 15,8 °C a 50 cm de profundidade. A presença de cobertura de palha mantém o solo mais aquecido que o solo sem palha no início do dia, mas impede que ele se aqueça igual ao solo sem palha no restante do dia. Assim resultando em uma percepção de que a palha aquece e resfria o solo em relação à temperatura do solo desnudo. Essa temperatura foi semelhante em todas as quantidades de palha, porém decresceu em profundidade. A estimativa da temperatura da superfície do solo pode ser estimada em função da temperatura do ar e da radiação solar. Em solo descoberto as equações propostas tem melhor desempenho que a equação senoidal. Em solos cobertos a equação senoidal teve melhor desempenho. Os parâmetros hidráulicos e térmicos do solo otimizados melhoraram as simulações da temperatura e conteúdo de água do solo. A simulação da temperatura do solo com diferentes quantidades de cobertura vegetal pode ser simulada usando o Hydrus-1D, sendo que o uso da condição de contorno com a temperatura da superfície medido teve melhor desempenho que o utilizando a temperatura da superfície estimada.

Land application of wastes has become increasingly popular, to promote nutrient recycling and environmental protection, with soil functioning as a partial barrier between wastes and groundwater. Dairy shed effluent (DSE), may contain a wide variety of pathogenic micro-organisms, including bacteria (e.g. Salmonella paratyphyi, Escherichia coli. and Campylobacter), protozoa and viruses. Groundwater pathogen contamination resulting from land-applied DSE is drawing more attention with the intensified development of the dairy farm industry in New Zealand. The purpose of this research was to investigate the fate and transport of bacterial indicator-faecal coliform (FC) from land-applied DSE under different irrigation practices via field lysimeter studies, using two water irrigation methods (flood and sprinkler) with contrasting application rates, through the 2005-2006 irrigation season. It was aimed at better understanding, quantifying and modelling of the processes that govern the removal of microbes in intact soil columns, bridging the gap between previous theoretical research and general farm practices, specifically for Templeton soil. This study involved different approaches (leaching experiments, infiltrometer measurements and a dye infiltration study) to understand the processes of transient water flow and bacterial transport; and to extrapolate the relationships between bacterial transport and soil properties (like soil structure, texture), and soil physical status (soil water potential ψ and volumetric water content θ). Factors controlling FC transport are discussed. A contaminant transport model, HYDRUS-1D, was applied to simulate microbial transport through soil on the basis of measured datasets. This study was carried out at Lincoln University’s Centre for Soil and Environmental Quality (CSEQ) lysimeter site. Six lysimeters were employed in two trials. Each trial involved application of DSE, followed by a water irrigation sequence applied in a flux-controlled method. The soil columns were taken from the site of the new Lincoln University Dairy Farm, Lincoln, Canterbury. The soil type is Templeton fine sandy loam (Udic-Ustochrept, coarse loamy, mixed, mesic). Vertical profiles (at four depths) of θ and ψ were measured during leaching experiments. The leaching experiments directly measured concentrations of chemical tracer (Br⁻ or Cl⁻) and FC in drainage. Results showed that bacteria could readily penetrate through 700 mm deep soil columns, when facilitated by water flow. In the first (summer) trial, FC in leachate as high as 1.4×10⁶ cfu 100 mL⁻¹ (similar to the DSE concentration), was detected in one lysimeter that had a higher clay content in the topsoil, immediately after DSE application, and before any water irrigation. This indicates that DSE flowed through preferential flow paths without significant treatment or reduction in concentrations. The highest post-irrigation concentration was 3.4×10³ cfu 100 mL⁻¹ under flood irrigation. Flood irrigation resulted in more bacteria and Br⁻ leaching than spray irrigation. In both trials (summer and autumn) results showed significant differences between irrigation treatments in lysimeters sharing similar drainage class (moderate or moderately rapid). Leaching bacterial concentration was positively correlated with both θ and ψ, and sometimes drainage rate. Greater bacterial leaching was found in the one lysimeter with rapid whole-column effective hydraulic conductivity, Keff, for both flood and spray treatments. Occasionally, the effect of Keff on water movement and bacterial transport overrode the effect of irrigation. The ‘seasonal condition’ of the soil (including variation in initial water content) also influenced bacterial leaching, with less risk of leaching in autumn than in summer. A tension infiltrometer experiment measured hydraulic conductivity of the lysimeters at zero and 40 mm suction. The results showed in most cases a significant correlation between the proportion of bacteria leached and the flow contribution of the macropores. The higher the Ksat, the greater the amount of drainage and bacterial leaching obtained. This research also found that this technique may exclude the activity of some continuous macropores (e.g., cracks) due to the difference of initial wetness which could substantially change the conductivity and result in more serious bacterial leaching in this Templeton soil. A dye infiltration study showed there was great variability in water flow patterns, and most of the flow reaching deeper than 50 cm resulted from macropores, mainly visible cracks. The transient water flow and transport of tracer (Br⁻) and FC were modelled using the HYDRUS-1D software package. The uniform flow van Genuchten model, and the dual-porosity model were used for water flow and the mobile-immobile (MIM) model was used for tracer and FC transport. The hydraulic and solute parameters were optimized during simulation, on the basis of measured datasets from the leaching experiments. There was evidence supporting the presence of macropores, based on the water flow in the post-DSE application stage. The optimised saturated water content (θs) decreased during the post-application process, which could be explained in terms of macropore flow enhanced by irrigation. Moreover, bacterial simulation showed discrepancies in all cases of uniform flow simulations at the very initial stage, indicating that non-equilibrium processes were dominant during those short periods, and suggesting that there were strong dynamic processes involving structure change and subsequently flow paths. It is recommended that management strategies to reduce FC contamination following application of DSE in these soils must aim to decrease preferential flow by adjusting irrigation schemes. Attention needs to be given to a) decreasing irrigation rates at the beginning of each irrigation; b) increasing the number of irrigations, by reducing at the same time the amount of water applied and the irrigation rate at each irrigation; c) applying spray irrigation rather than flood irrigation.

The distributions of calcite rich horizons within dryland soils are commonly used as paleoclimate proxies. Comprehensive conceptual and mathematical models of calcite accumulation in soils are required to accurately interpret and calibrate these proxies. A conceptual model for calcite accumulation is already well established: As water percolates through a soil, it dissolves minerals, such as calcite, transporting the soluble minerals downward. As soil water is removed by evaporation and transpiration, the water solution becomes supersaturated resulting in precipitation of calcite at depth. The impacts of dynamic plant growth and microbial respiration have not yet been simulated in numerical models for calcite accumulation but are likely important because of their influence on variables governing calcite solubility. The soil hydraulic modeling software, HYDRUS-1D, simulates water and solute transfer through a soil column, accounting for variations in all previously studied variables (temperature, water content, soil pCO₂) while additionally simulating vegetation-soil interactions. Five separate sensitivity studies were conducted to determine the importance for calcite accumulation of 1) soil texture, 2) plant growth, 3) plant phenology, 4) atmospheric CO₂ concentrations, and 5) the proximal variables that control calcite dissolution and precipitation: soil CO₂, soil water content, and soil temperature. In each modeling simulation, calcite was leached from the top several cm and redistributed deeper in the soil after 20 years. Soils with courser texture yield deeper (+20cm), more diffuse calcite horizons, as did simulations with bare soil compared to vegetated soil. The phenology of plant communities (late spring versus late summer growth) resulted in soil calcite accumulation at temperatures differing by at least 10°C. Changes in atmospheric CO₂ concentrations do not affect the soil calcite distribution. Variations in concentration of soil CO₂, rather than soil water content, have the greatest direct effect on calcite solubility. The most significant time periods of annual accumulation also corresponded with positive water fluxes resulting from high matric potential at the surface. Transpiration and evaporation moisture sinks caused solution to travel upward from higher to lower soil CO₂ concentrations, causing CO₂ de-gassing and calcite accumulation. This pathway describes a new qualitative mechanism for soil calcite formation and should be included in the conceptual model.
text

碩士
國立中央大學
應用地質研究所
98
Infiltration test is an important technique to investigate the hydrologic properties in soils, while the unsaturated flow numerical models can be an efficient tool to understand flow mechanism when the physical properties and conditions of soils are changed. This study employs HYDRUS-1D model and self developed Levenberg-Marquardt algorithm to inversely estimate parameters in the van Genuchten formula. The developed LM model was tested based on the synthetically generated data from HYDRUS-1D model. Additionally, the initial guesses of soil parameters in LM model were systematically analyzed to obtain general insight into the laboratory and field scale applications. The tested model was then used to estimate soil parameters of laboratory sand box under different infiltration scenarios. During the infiltration tests, the suction or pressure and water content were measured with tensiometers TDR(Time Domain Reflectometry). Such measured data can thus be used in the LM model to estimate the soil parameters. The estimated parameters were compared with those based on trial and error approach. Field scale double-ring infiltration was tested on the river bank of Zhuoshui river near His-Luo. The saturated infiltration rate(fc), initial infiltration rate (f0), and infiltration parameter by using the in Horton formula were estimated by the developed LM model. The objective of the field scale test is to obtain the saturated infiltration rate(f0) to compare with the results of laboratory experiment. Estimate Van Genuchten parameters by combining the laboratory infiltration experiment of field soil sample and saturated infiltration parameters. The results of the study reveal that the developed LM model associated with HYDRUS-1D model can accurately estimate van Genuchten soil parameters. The results of model test shows that the lower initial guess values would lead to better estimations of residual water content (θr) and α values. The saturated hydraulic conductivity (Ks) is the most sensitive parameter while the initial guess n values are relatively insignificant. The LM inversion of tests show that LM model associated with HYDRUS-1D model can estimate accurately van Genuchten for laboratory scale tests soil parameters. The field scale multiple depth double-ring infiltration tests reveals that the saturated infiltration rates(f0) of Zhuoshui river near His-Luo have one order of magnitude difference for different depths such depth difference of hydraulic conductivity (Ks) (0.0014 and 0.00015 cm/s) were also shown in laboratory scale tests. Based on the soil sample from Zhuoshui river the sand box experiment data were then used in LM model to estimate soil parameters. This simulation results showed that the value of α is from 0.01688 to 0.09289(cm-1), the value of n from 2.528 to 5.886(-), the material of the soil can be classified as sand and silt.

碩士
國立臺灣大學
生物環境系統工程學研究所
107
Due to the impact of climate change and the rising frequency of extreme weather have caused the problem of food security and water resources. Precision agriculture in greenhouses is gradually becoming a trend because it is less affected by external environmental factors and easier to manage. In order to effectively improve the utilization of water resources and reduce wasting water, this study planted cherry tomatoes in the greenhouse and installed monitoring instruments to collect data. The soil temperature and volumetric water content during tomato growth can be predicted effectively by combining the weather forecast. The Hydrus-1D, random forest method and ICON were used to simulate the change of soil temperature and volumetric water content in the greenhouse. The results of the simulation show good performances among the three models. Further combined with weather forecast data, three models were applied to predict soil temperatures and volumetric moisture contents. Both of predicted soil temperature and volumetric water content by Hydrus-1D followed the trend of real measured data, but Hydrus-1D can’t predict the watering situation. The predicted performance of soil temperature by ICON was low. Nevertheless, both predictions by random forest method were better than others. This study simulated and predicted the soil temperatures and volumetric water contents by applying Hydrus-1D, random forest method and ICON to improve greenhouse management and increase water-use efficiency.

Indiana University-Purdue University Indianapolis (IUPUI)
Non-rainfall waters such as fog and dew are considered as important source of water in drylands, and the knowledge of possible sources of its formation is very important to make future predictions. Prior studies have suggested the presence of radiation fog in drylands; however, its formation mechanism still remains unclear. There have been earlier studies on the effects of fog on soil moisture dynamics and groundwater recharge. On the contrary, no research has yet been conducted to understand the contribution of soil moisture and groundwater to fog formation. This study, therefore, for the first time intends to examine such possibility in a fog-dominated dryland ecosystem, the Namib Desert. The study was conducted at three sites representing two different land forms (sand dunes and gravel plains) in the Namib Desert. This thesis is divided into two parts: the first part examines evidences of fog formation through water vapor movement using field observations, and the second part simulates water vapor transport using HYDRUS-1D model. In the first part of the study, soil moisture, soil temperature and air temperature data were analyzed, and the relationships between these variables were taken as one of the key indicators for the linkage between soil water and fog formation. The analysis showed that increase in soil moisture generally corresponds to similar increase in air or soil temperature near the soil surface, which implied that variation in soil moisture might be the result of water vapor movement (evaporated soil moisture or groundwater) from lower depths to the soil surface. In the second part of the study, surface fluxes of water vapor were simulated using the HYDRUS-1D model to explore whether the available surface flux was sufficient to support fog formation. The actual surface flux and cumulative evaporation obtained from the model showed positive surface fluxes of water vapor. Based on the field observations and the HYDRUS-1D model results, it can be concluded that water vapor from soil layers and groundwater is transported through the vadose zone to the surface and this water vapor likely contributes to the formation of non-rainfall waters in fog-dominated drylands, like the Namib Desert.