Submitted by (lucia.rodrigues@ufrpe.br) on 2016-10-05T13:49:43Z No. of bitstreams: 1 Jose Graciliano da Silva Junior.pdf: 1829762 bytes, checksum: 76597753b818f8772c987849bed6f86e (MD5) / Made available in DSpace on 2016-10-05T13:49:43Z (GMT). No. of bitstreams: 1 Jose Graciliano da Silva Junior.pdf: 1829762 bytes, checksum: 76597753b818f8772c987849bed6f86e (MD5) Previous issue date: 2009-03-09 / Conselho Nacional de Pesquisa e Desenvolvimento Científico e Tecnológico - CNPq / In the Brazilian semi-arid the family agriculture comprises a financial support to small poor communities, contributing for reduction of the rural exodus. However, due to the adverse climatic conditions, as the concentration of the rainfall in some months during the year, the irrigation becomes indispensable for the practising agriculture. The objective of the work was of gaging the model SALTMED during the station of growth of the cabbage so that the same can be used to analyze alternatives of integrated handling of irrigation in systems of the semiarid considering aspects of the natural resources (soil and water), climate, culture and overhead irrigation in a frequency of watering of 24 hours and to validate him/it frequently for the treatments of 12 and 48 hours in areas with soil of the alluvial type of medium texture with salinity potentials. The experiment was driven in irrigated lot of Our Finance Mrs. of Rosario, municipal district of Pesqueira-PE, Rural of the State of Pernambuco. The irrigation sheets, for the frequencies of 12, 24 and 48 hours, they were calculated with base in the culture coefficient (Kc) and in the reference evapotranspiração of the (Eto) from pan. During the cycle of the culture the humidity contents were monitored in the depths of 20 and 40 cm that were compared with values made calculations by SALTMED for calibration and validation. Also the salinity of the soil was measured, for the treatment frequently of 24 hours. The productivity of the culture of the cabbage was evaluated for the different frequencies. It was verified experimentally that the three frequencies used in the study can be used as handling alternative for reduction of the salinity of the soil, with prominence for them frequency of 12 and 48 hours. The results also showed that the model SALTMED answered satisfactorily in the prediction of the humidity data and of salinity of the soil, for the tested conditions. Therefore the model can be used for forecast of handling of similar irrigation, in way to provide larger productivities of cultures sensitive olerícolas the salinity, varying sheet of water, quality of the water and types of overhead irrigations. / No semi-árido brasileiro, a pequena agricultura familiar constitui suporte financeiro a pequenas comunidades pobres, vindo a contribuir para redução do êxodo rural. Entretanto, devido às condições climáticas adversas pela concentração da chuva em alguns meses durante o ano, a irrigação se torna indispensável para o desenvolvimento da agricultura. O objetivo desse trabalho foi de calibrar o modelo SALTMED durante a estação de crescimento do repolho para que o mesmo possa ser utilizado para analisar alternativas de manejo integrado de irrigação em sistemas do semi-árido considerando aspectos dos recursos naturais (solo e água), clima, cultura e sistema de irrigação numa freqüência de rega de 24 horas e validá-lo para os tratamentos com freqüência de 12 e 48 horas em áreas com solo do tipo aluvial de textura média com susceptibilidade à salinização. O experimento foi conduzido em lote irrigado da Fazenda Nossa Senhora do Rosário, município de Pesqueira, Agreste do Estado de Pernambuco. As lâminas de irrigação, para as frequencias de 12, 24 e 48 horas, foram calculadas com base no coeficiente de cultura (Kc) e na evapotranspiração de referência do Tanque Classe A (Eto). Durante o ciclo da cultura foram monitorados os conteúdos de umidade nas profundidades de 20 e 40 cm que foram comparados com valores calculados pelo SALTMED para calibração e validação. Também foi medida a salinidade do solo, para o tratamento com frequência de 24 horas. A produtividade da cultura do repolho foi avaliada para as diferentes frequências.Verificou-se experimentalmente que as três freqüências utilizadas no estudo podem ser utilizadas como alternativa de manejo para redução da salinidade do solo, com destaque para as frequência de 12 e 48 horas. Os resultados também mostraram que o modelo SALTMED respondeu satisfatoriamente na predição dos dados de umidade e de salinidade do solo, para as condições testadas. Portanto o modelo pode ser utilizado para previsão de manejo de irrigação semelhante, de maneira a proporcionar maiores produtividades de culturas olerícolas sensíveis a salinidade, variando lâmina de água, qualidade da água e tipos de sistemas de irrigação.

Este trabalho estudou o impacto do nitrato, proveniente de uma fossa negra desativada nas zonas não-saturada e saturada no Aqüífero Adamantina (Urânia, SP), através do acompanhamento hidrogeoquímico dos íons maiores (incluindo a série nitrogenada), menores, isótopos estáveis (15NNO3, 15NN2O, 18ON2O, 18ONO3, 18OO2, 13CCO2) e gases (O2, CO2, N2O e CH4). Os três anos de monitoramento da fossa negra, desativada desde 2002, se deu pela construção de uma estação experimental, constituída de um poço cacimba com 11,20 m de profundidade. Foram instalados, ao longo de sua parede até a zona saturada, 12 tensiômetros e 12 lisímetros de sucção (0,5 a 9,0 m), bem como um poço de monitoramento instalado no fundo da estação, para o monitoramento dos processos hidráulicos e da qualidade da água infiltrada. Os resultados obtidos sobre a hidráulica de infiltração mostraram que a zona nãosaturada apresenta condutividade hidráulica média vertical entre 3,1x10-6 a 1,4x10-5 m/s. As tensões capilares medidas indicaram a ocorrência de variações laterais e em profundidade, cuja origem está relacionada às heterogeneidades do material geológico. Tais constatações possibilitaram a identificação de três diferentes zonas hidráulicas: a primeira, entre 0,5 e 2,0 m; a segunda, de 2,0 a 7,0 m e a terceira, a partir de 7,0 m até o nível d´água. A taxa de recarga na área é de 171 mm, o que corresponde a 13% da precipitação anual. As velocidades das frentes de molhamento indicaram um tempo de trânsito de 90 dias para a água se deslocar de 0,5 até 9,0 m de profundidade. As análises químicas da água ao longo do perfil da zona não-saturada e do poço de monitoramento, apontaram concentrações altas de nitrato ao longo de todo o perfil (até 2028,94 mg/L). Os parâmetros físico-químicos, isotópicos e de gases analisados apresentaram variações, sazonais e em profundidade, de suas concentrações. Os perfis de concentração de gases, isótopos e série nitrogenada indicam que estes parâmetros apresentam um comportamento ?zigue-zague? sendo os maiores valores a 2,0; 4,0 e 6,0 m de profundidade. Este comportamento seria condicionado pelas características texturais do material geológico, o que leva a considerar a existência de áreas de menor permeabilidade nestes pontos, restringindo o acesso de águas ricas em oxigênio e de gases ao longo da zona não-saturada e co-existência de zonas redutoras em meio a zonas oxidantes (microcosmos). Este monitoramento detalhado possibilitou o estabelecimento dos principais processos hidrogeoquímicos atuantes na zona não-saturada: i) oxidação da matéria orgânica, ii) amonificação, iii) nitrificação, iv) metanogênese e v) desnitrificação. Um fator muito importante a ser levantado é a simultaneidade entre os mecanismos de nitrificação e desnitrificação na zona não-saturada. Várias foram as constatações que reforçam esta hipótese, exemplificadas pelas seguintes relações: i) 15NN2O e 18ON2O; ii) 15NN2O e 15NNO3 e iii) 18ONO3 e15NNO3. Os valores de ?15NNO3 obtidos foram de +6,80 a +30,09? (estação úmida) e +9,54 a +23,25? (estação seca). Valores mais enriquecidos de ?15NNO3, comparados à assinatura isotópica reportada para dejetos humanos (+8,1? a +13,1?), em conjunto com os resultados de ?18ONO3 (-3,60 a 4,50 ?), razão de fracionamento destes isótopos e os coeficientes de regressão, obtidos em outros trabalhos, reforçaram a coexistência de desnitrificação e nitrificação na área estudada. Os resultados de ?15NN2O e ?18OO2 também são indicativos de desnitrificação. As assinaturas destes dois isótopos (-16,16 a -11,94? e 28,05 a 30,69?, respectivamente), confrontados aos valores reportados em outros trabalhos (-37 a -11? e -21 a +57?), indicam que o N2O, detectado na área, seria produzido por esse processo. Contudo, cabe salientar que o uso de isótopos de 18OO2 e 15NN2O constitui-se ainda em técnica muito recente em estudos hidrogeológicos, carecendo de trabalhos que stabeleçam as assinaturas isotópicas para os processos de nitrificação e desnitrificação considerando as diferentes fontes de contaminação de nitrato. / The present work studied the impact of nitrate from an inactive cesspit to the unsaturated and saturated zones of the Adamantina Aquifer (Urânia, São Paulo, Brazil), through the hydrogeochemistry of major ions (including the nitrogen series), minor ions, stable isotopes (15NNO3, 15NN2O, 18ON2O, 18ONO3, 18OO2, 13CCO2) and gases (O2, CO2, N2O and CH4). The cesspit, which has been inactive since 2002, was monitored for three years after the construction of a monitoring station. The station is comprised by a 11.20m deep dug well, in which 12 tensiometers and 12 suction lysimeters were installed from the surface down to the saturated zone (0.5-9.0 m). A monitoring well was installed at the bottom of the monitoring station. The station was used to examine the hydraulic processes and the quality of the infiltrating water. The results from the hydraulics of the water infiltration show that the unsaturated zone presents average hydraulic conductivity between 3.1x10-6 and 1.4x10-5 m/s. The capillary tension varies laterally and vertically, which is attributed to the heterogeneities of the geological material. These results allowed the identification of three different hydraulic zones: the first one, between the depths of 0.5 and 2.0 m; the second one, between 2.0 and 7.0 m and the third one, from the depth of 7.0 m down to the water level. The recharge rate of the area is 171 mm, which corresponds to 13% of the annual precipitation. The velocities of the wetting fronts indicated a travel time of 90 days for the water to move between the depths of 0.5 m and 9.0m. The chemical analyses of the water of the unsaturated and saturated zones indicated high concentration of nitrate, along the entire unsaturated zone (up to 2028,94 mg/L). The physico-chemical, isotopic and gases data vary seasonally and spatially. The concentration of gases, isotopes and nitrogen species show a ?zig-zag? behavior, with peaks of concentration at the depths 2.0; 4.0 and 6.0 m. This behavior is a consequence of the textural features of the geological material, with less permeable areas occuring at these depths, which can restrict the access of oxygen- and gases-rich waters along the unsaturated zone and the coexistence of reducing zones among oxidizing ones (microcosms). The detailed monitoring provided data to establish the main hydrogeochemical processes occuring in the unsaturated zone: i) oxidation of organic matter, ii) ammonification, iii) nitrification, iv) methanogenesis and v) denitrification. A very important fact is the simultaneous occurrence of nitrification and denitrification processes in the unsaturated zone. This fact is verified by a series of observations, which can be demonstrated by the following relationships: i) 15NN2O and 18ON2O; ii) 15NN2O and 15NNO3 and iii) 18ONO3 and15NNO3. The values obtained for ?15NNO3 are between +6.80 and +30.09? (wet season) and between +9.54 and +23.25? (dry season). The values of ?15NNO3, which are more enriched than the isotopic fingerprint for sewage (+8.1? to +13.1?), along with the results for ?18ONO3 (-3.60 to 4.50 ?), the fractionation ratio of this isotopes and the linear regression coefficient, obtained by other authors, reinforced the coexistence of denitrification and nitrification in the studied area. The results of ?15NN2O and ?18OO2 also indicated the presence of denitrification. The signatures of these two isotopes (-16.16 to -11.94? and 28.05 to 30.69?, respectivelly), compared to values reported by other authors (-37 to -11? and -21 to +57?), indicate that the N2O detected in the area is produced by this process. However, it is worth emphasizing that the use of 18OO2 and 15NN2O isotopes is still a very new technique in hydrogeological studies, and more studies are needed to establish the isotopic fingerprints of the nitrification and denitrification processes, considering the different sources of nitrate contamination.

Este trabalho estudou o impacto do nitrato, proveniente de uma fossa negra desativada nas zonas não-saturada e saturada no Aqüífero Adamantina (Urânia, SP), através do acompanhamento hidrogeoquímico dos íons maiores (incluindo a série nitrogenada), menores, isótopos estáveis (15NNO3, 15NN2O, 18ON2O, 18ONO3, 18OO2, 13CCO2) e gases (O2, CO2, N2O e CH4). Os três anos de monitoramento da fossa negra, desativada desde 2002, se deu pela construção de uma estação experimental, constituída de um poço cacimba com 11,20 m de profundidade. Foram instalados, ao longo de sua parede até a zona saturada, 12 tensiômetros e 12 lisímetros de sucção (0,5 a 9,0 m), bem como um poço de monitoramento instalado no fundo da estação, para o monitoramento dos processos hidráulicos e da qualidade da água infiltrada. Os resultados obtidos sobre a hidráulica de infiltração mostraram que a zona nãosaturada apresenta condutividade hidráulica média vertical entre 3,1x10-6 a 1,4x10-5 m/s. As tensões capilares medidas indicaram a ocorrência de variações laterais e em profundidade, cuja origem está relacionada às heterogeneidades do material geológico. Tais constatações possibilitaram a identificação de três diferentes zonas hidráulicas: a primeira, entre 0,5 e 2,0 m; a segunda, de 2,0 a 7,0 m e a terceira, a partir de 7,0 m até o nível d´água. A taxa de recarga na área é de 171 mm, o que corresponde a 13% da precipitação anual. As velocidades das frentes de molhamento indicaram um tempo de trânsito de 90 dias para a água se deslocar de 0,5 até 9,0 m de profundidade. As análises químicas da água ao longo do perfil da zona não-saturada e do poço de monitoramento, apontaram concentrações altas de nitrato ao longo de todo o perfil (até 2028,94 mg/L). Os parâmetros físico-químicos, isotópicos e de gases analisados apresentaram variações, sazonais e em profundidade, de suas concentrações. Os perfis de concentração de gases, isótopos e série nitrogenada indicam que estes parâmetros apresentam um comportamento ?zigue-zague? sendo os maiores valores a 2,0; 4,0 e 6,0 m de profundidade. Este comportamento seria condicionado pelas características texturais do material geológico, o que leva a considerar a existência de áreas de menor permeabilidade nestes pontos, restringindo o acesso de águas ricas em oxigênio e de gases ao longo da zona não-saturada e co-existência de zonas redutoras em meio a zonas oxidantes (microcosmos). Este monitoramento detalhado possibilitou o estabelecimento dos principais processos hidrogeoquímicos atuantes na zona não-saturada: i) oxidação da matéria orgânica, ii) amonificação, iii) nitrificação, iv) metanogênese e v) desnitrificação. Um fator muito importante a ser levantado é a simultaneidade entre os mecanismos de nitrificação e desnitrificação na zona não-saturada. Várias foram as constatações que reforçam esta hipótese, exemplificadas pelas seguintes relações: i) 15NN2O e 18ON2O; ii) 15NN2O e 15NNO3 e iii) 18ONO3 e15NNO3. Os valores de ?15NNO3 obtidos foram de +6,80 a +30,09? (estação úmida) e +9,54 a +23,25? (estação seca). Valores mais enriquecidos de ?15NNO3, comparados à assinatura isotópica reportada para dejetos humanos (+8,1? a +13,1?), em conjunto com os resultados de ?18ONO3 (-3,60 a 4,50 ?), razão de fracionamento destes isótopos e os coeficientes de regressão, obtidos em outros trabalhos, reforçaram a coexistência de desnitrificação e nitrificação na área estudada. Os resultados de ?15NN2O e ?18OO2 também são indicativos de desnitrificação. As assinaturas destes dois isótopos (-16,16 a -11,94? e 28,05 a 30,69?, respectivamente), confrontados aos valores reportados em outros trabalhos (-37 a -11? e -21 a +57?), indicam que o N2O, detectado na área, seria produzido por esse processo. Contudo, cabe salientar que o uso de isótopos de 18OO2 e 15NN2O constitui-se ainda em técnica muito recente em estudos hidrogeológicos, carecendo de trabalhos que stabeleçam as assinaturas isotópicas para os processos de nitrificação e desnitrificação considerando as diferentes fontes de contaminação de nitrato. / The present work studied the impact of nitrate from an inactive cesspit to the unsaturated and saturated zones of the Adamantina Aquifer (Urânia, São Paulo, Brazil), through the hydrogeochemistry of major ions (including the nitrogen series), minor ions, stable isotopes (15NNO3, 15NN2O, 18ON2O, 18ONO3, 18OO2, 13CCO2) and gases (O2, CO2, N2O and CH4). The cesspit, which has been inactive since 2002, was monitored for three years after the construction of a monitoring station. The station is comprised by a 11.20m deep dug well, in which 12 tensiometers and 12 suction lysimeters were installed from the surface down to the saturated zone (0.5-9.0 m). A monitoring well was installed at the bottom of the monitoring station. The station was used to examine the hydraulic processes and the quality of the infiltrating water. The results from the hydraulics of the water infiltration show that the unsaturated zone presents average hydraulic conductivity between 3.1x10-6 and 1.4x10-5 m/s. The capillary tension varies laterally and vertically, which is attributed to the heterogeneities of the geological material. These results allowed the identification of three different hydraulic zones: the first one, between the depths of 0.5 and 2.0 m; the second one, between 2.0 and 7.0 m and the third one, from the depth of 7.0 m down to the water level. The recharge rate of the area is 171 mm, which corresponds to 13% of the annual precipitation. The velocities of the wetting fronts indicated a travel time of 90 days for the water to move between the depths of 0.5 m and 9.0m. The chemical analyses of the water of the unsaturated and saturated zones indicated high concentration of nitrate, along the entire unsaturated zone (up to 2028,94 mg/L). The physico-chemical, isotopic and gases data vary seasonally and spatially. The concentration of gases, isotopes and nitrogen species show a ?zig-zag? behavior, with peaks of concentration at the depths 2.0; 4.0 and 6.0 m. This behavior is a consequence of the textural features of the geological material, with less permeable areas occuring at these depths, which can restrict the access of oxygen- and gases-rich waters along the unsaturated zone and the coexistence of reducing zones among oxidizing ones (microcosms). The detailed monitoring provided data to establish the main hydrogeochemical processes occuring in the unsaturated zone: i) oxidation of organic matter, ii) ammonification, iii) nitrification, iv) methanogenesis and v) denitrification. A very important fact is the simultaneous occurrence of nitrification and denitrification processes in the unsaturated zone. This fact is verified by a series of observations, which can be demonstrated by the following relationships: i) 15NN2O and 18ON2O; ii) 15NN2O and 15NNO3 and iii) 18ONO3 and15NNO3. The values obtained for ?15NNO3 are between +6.80 and +30.09? (wet season) and between +9.54 and +23.25? (dry season). The values of ?15NNO3, which are more enriched than the isotopic fingerprint for sewage (+8.1? to +13.1?), along with the results for ?18ONO3 (-3.60 to 4.50 ?), the fractionation ratio of this isotopes and the linear regression coefficient, obtained by other authors, reinforced the coexistence of denitrification and nitrification in the studied area. The results of ?15NN2O and ?18OO2 also indicated the presence of denitrification. The signatures of these two isotopes (-16.16 to -11.94? and 28.05 to 30.69?, respectivelly), compared to values reported by other authors (-37 to -11? and -21 to +57?), indicate that the N2O detected in the area is produced by this process. However, it is worth emphasizing that the use of 18OO2 and 15NN2O isotopes is still a very new technique in hydrogeological studies, and more studies are needed to establish the isotopic fingerprints of the nitrification and denitrification processes, considering the different sources of nitrate contamination.

Devido aos riscos de contaminação dos recursos naturais solo e água, ao alto custo, ao tempo e ao esforço humano nas investigações de campo, os modelos matemáticos, aliados às técnicas numéricas e aos avanços computacionais, constituem uma ferramenta importante na previsão do deslocamento de solutos, contribuindo assim, para o controle de alterações ambientais. No Brasil, a modelação de fluxo e transporte de solutos na zona não-saturada é voltada, quase que exclusivamente, aos problemas relacionados às atividades agrícolas. Entretanto, tão importante quanto a problemática dos produtos químicos nas atividades agrícolas é a questão de poluição e contaminação do solo e da água por chorume, gerado pelos resíduos sólidos domiciliares. Neste trabalho, é desenvolvido e validado um modelo computacional unidimensional para simulação de fluxo e transporte de solutos na zona não-saturada do solo. O modelo matemático é dado pela equação diferencial parcial não-linear de Richards, que rege o movimento de água no solo, e a equação diferencial parcial linear de advecção-dispersão, do transporte de solutos, acompanhadas das condições iniciais e de contorno. A equação de Richards é dada em função do potencial matricial da água e a equação de transporte de solutos estima a evolução temporal da concentração de solutos no perfil do solo. Devido à dificuldade de se obter soluções analíticas destas equações, são resolvidas numericamente pelo método de elementos finitos. As referidas equações são resolvidas utilizando-se malhas uniformes inicialmente. Com a finalidade de obter simulações mais eficientes, a um custo computacional reduzido, é empregada a adaptatividade com refinamento h na malha de elementos finitos. A função interpolação polinomial utilizada é de grau 2 ou maior que garante a conservação de massa. Na equação de Richards, a derivada temporal é aproximada por um quociente de diferença finita e é aplicado o esquema de Euler explícito e na equação de advecção-dispersão, é aproximada por um quociente de diferença finita, aplicando-se o esquema de Euler implícito, devido à linearidade da equação. O sistema operacional é o Linux Ubuntu 32 bits, o ambiente de programação é o PZ, escrito em linguagem de programação C++. Na validação do modelo, utilizam-se dados disponíveis na literatura. Os resultados são comparados, utilizando-se malhas uniformes e malhas adaptativas com refinamento h. Usando-se as malhas uniformes para o problema de Richards e de transporte de potássio, o tempo de execução é de 22 minutos e a memória utilizada de 6164 Kb. Com as malhas adaptadas, o tempo de execução é de 3 minutos e 27 segundos, consumindo 5876 Kb de memória. Houve, portanto, uma redução de 84,32% no tempo de execução, usando-se malhas adaptativas. A utilização da função interpolação polinomial de grau 2 ou maior e o refinamento h, permitem uma boa concordância do modelo na comparação com soluções disponíveis na literatura. / Due to the risks of contamination of soil and water resources, the high cost, time and human effort in the field investigations, the mathematical models, combined with numerical techniques and computational advances, are important tools in forecasting the movement of solutes thereby contributing to the control of environmental alteration. In Brazil, modeling of flow and solute transport in the unsaturated zone is focused, almost exclusively, on problems related to agricultural activities. However, as important as the problematical of chemicals products in agricultural activities is the issue of pollution and contamination of soil and water by leachate, generated by municipal solid wastes. In this work, an one-dimensional computational model for simulation of flow and solute transport in the unsaturated soil has been developed and validated. The mathematical model is given by the Richards\'s non-linear partial differential equation, which determines the movement of water in the soil, and the advection-dispersion linear partial differential equation, of the solute transport, together with initial and boundary conditions. The Richards equation is a function of the water pressure head and the solute transport equation estimate the temporal evolution of the solutes concentration in the soil profile. Due to the difficulty of obtaining analytical solutions of these equations, they are solved numerically using the finite element method. The governing equations are solved using initially a uniform mesh. In order to obtain more efficient simulations with low computational cost, adaptativity with h refinement on the finite element mesh is implemented. The interpolation function is of degree two or higher, assuring mass conservation. In Richards\' equation, the temporal derivative is approximated by Euler explicit finite difference. For the advection-dispersion equation, due to the linearity of the equation, an implicit finite difference scheme is used. The code is written in the programming language C++ based on the programming environment PZ using operating system Linux Ubuntu 32 bit. Model results are validated in comparison with data available in the literature. The results are evaluated using uniform meshes and with h refinement adaptive mesh. Using the uniform meshes for the problem of Richards and transport of potassium, the running time is 22 minutes and 6164 Kb of memory is used. With the adapted meshes, the execution time is 3 minutes and 27 seconds, consuming 5,876 Kb of memory. Therefore there was a reduction of 84.32% in execution time, using adaptive meshes. The interpolation function with degree two or higher and the h refinement, with reduction of the computation time, showed a good agreement in comparison with the literature.

Devido aos riscos de contaminação dos recursos naturais solo e água, ao alto custo, ao tempo e ao esforço humano nas investigações de campo, os modelos matemáticos, aliados às técnicas numéricas e aos avanços computacionais, constituem uma ferramenta importante na previsão do deslocamento de solutos, contribuindo assim, para o controle de alterações ambientais. No Brasil, a modelação de fluxo e transporte de solutos na zona não-saturada é voltada, quase que exclusivamente, aos problemas relacionados às atividades agrícolas. Entretanto, tão importante quanto a problemática dos produtos químicos nas atividades agrícolas é a questão de poluição e contaminação do solo e da água por chorume, gerado pelos resíduos sólidos domiciliares. Neste trabalho, é desenvolvido e validado um modelo computacional unidimensional para simulação de fluxo e transporte de solutos na zona não-saturada do solo. O modelo matemático é dado pela equação diferencial parcial não-linear de Richards, que rege o movimento de água no solo, e a equação diferencial parcial linear de advecção-dispersão, do transporte de solutos, acompanhadas das condições iniciais e de contorno. A equação de Richards é dada em função do potencial matricial da água e a equação de transporte de solutos estima a evolução temporal da concentração de solutos no perfil do solo. Devido à dificuldade de se obter soluções analíticas destas equações, são resolvidas numericamente pelo método de elementos finitos. As referidas equações são resolvidas utilizando-se malhas uniformes inicialmente. Com a finalidade de obter simulações mais eficientes, a um custo computacional reduzido, é empregada a adaptatividade com refinamento h na malha de elementos finitos. A função interpolação polinomial utilizada é de grau 2 ou maior que garante a conservação de massa. Na equação de Richards, a derivada temporal é aproximada por um quociente de diferença finita e é aplicado o esquema de Euler explícito e na equação de advecção-dispersão, é aproximada por um quociente de diferença finita, aplicando-se o esquema de Euler implícito, devido à linearidade da equação. O sistema operacional é o Linux Ubuntu 32 bits, o ambiente de programação é o PZ, escrito em linguagem de programação C++. Na validação do modelo, utilizam-se dados disponíveis na literatura. Os resultados são comparados, utilizando-se malhas uniformes e malhas adaptativas com refinamento h. Usando-se as malhas uniformes para o problema de Richards e de transporte de potássio, o tempo de execução é de 22 minutos e a memória utilizada de 6164 Kb. Com as malhas adaptadas, o tempo de execução é de 3 minutos e 27 segundos, consumindo 5876 Kb de memória. Houve, portanto, uma redução de 84,32% no tempo de execução, usando-se malhas adaptativas. A utilização da função interpolação polinomial de grau 2 ou maior e o refinamento h, permitem uma boa concordância do modelo na comparação com soluções disponíveis na literatura. / Due to the risks of contamination of soil and water resources, the high cost, time and human effort in the field investigations, the mathematical models, combined with numerical techniques and computational advances, are important tools in forecasting the movement of solutes thereby contributing to the control of environmental alteration. In Brazil, modeling of flow and solute transport in the unsaturated zone is focused, almost exclusively, on problems related to agricultural activities. However, as important as the problematical of chemicals products in agricultural activities is the issue of pollution and contamination of soil and water by leachate, generated by municipal solid wastes. In this work, an one-dimensional computational model for simulation of flow and solute transport in the unsaturated soil has been developed and validated. The mathematical model is given by the Richards\'s non-linear partial differential equation, which determines the movement of water in the soil, and the advection-dispersion linear partial differential equation, of the solute transport, together with initial and boundary conditions. The Richards equation is a function of the water pressure head and the solute transport equation estimate the temporal evolution of the solutes concentration in the soil profile. Due to the difficulty of obtaining analytical solutions of these equations, they are solved numerically using the finite element method. The governing equations are solved using initially a uniform mesh. In order to obtain more efficient simulations with low computational cost, adaptativity with h refinement on the finite element mesh is implemented. The interpolation function is of degree two or higher, assuring mass conservation. In Richards\' equation, the temporal derivative is approximated by Euler explicit finite difference. For the advection-dispersion equation, due to the linearity of the equation, an implicit finite difference scheme is used. The code is written in the programming language C++ based on the programming environment PZ using operating system Linux Ubuntu 32 bit. Model results are validated in comparison with data available in the literature. The results are evaluated using uniform meshes and with h refinement adaptive mesh. Using the uniform meshes for the problem of Richards and transport of potassium, the running time is 22 minutes and 6164 Kb of memory is used. With the adapted meshes, the execution time is 3 minutes and 27 seconds, consuming 5,876 Kb of memory. Therefore there was a reduction of 84.32% in execution time, using adaptive meshes. The interpolation function with degree two or higher and the h refinement, with reduction of the computation time, showed a good agreement in comparison with the literature.