| Summary: | O conhecimento da distribuição da umidade no solo para diferentes tipos de solos e diferentes vazões pode ser aplicado na otimização do dimensionamento de um sistema de irrigação, no manejo da água na zona radicular, bem como, auxiliar na aplicação eficiente de fertilizantes na rizosfera. Portanto, objetivou-se com a presente pesquisa desenvolver um modelo numérico bidimensional capaz de simular a distribuição da umidade no perfil de solo para diferentes sistemas de irrigação. Sendo que para tal, o modelo utiliza o método das diferenças finitas, mediante uma discretização da solução da equação de Richards e o método iterativo de Picard, o qual foi utilizado para garantir a conservação da massa. O modelo numérico ainda torna possível considerar a extração de água pelas plantas e a evaporação da água na superfície do solo. E para avaliar o seu desempenho foi feita uma análise de sensibilidade. O modelo matemático foi desenvolvido no Departamento de Engenharia de Biossistemas, pertencente à Escola Superior de Agricultura \"Luiz de Queiroz\" - ESALQ/USP. A linguagem de programação utilizada foi Fortran 90. A estrutura computacional foi feita de forma a permitir ao usuário a entrada de informações tais como: a) dados do perfil do solo, no que se refere às suas propriedades físico-hídricas, b) informações sobre o tipo de irrigação, tempo de simulação, tempo de aplicação de água via irrigação. O programa permite ao usuário a opção entre simular a extração de água pela as plantas ou não. Para a validação do modelo foram utilizados os dados obtidos por Rivera (2004), em um experimento conduzindo no antigo Departamento de Engenharia Rural da Escola Superior de Agricultura \"Luiz de Queiroz\". O desempenho do modelo foi avaliado com base nos parâmetros estatísticos: índice de concordância de Willmott; coeficiente de eficiência, raiz quadrada do erro médio, erro médio e o erro máximo absoluto. A análise de sensibilidade do modelo foi avaliada em função dos parâmetros: densidade de fluxo de água, condutividade hidráulica do solo saturado, alfa e n (parâmetros de ajuste da curva de retenção). Pelos resultados obtidos, o modelo apresentou um desempenho satisfatório na simulação dos perfis de umidade em relação aos dados medidos em condições de experimento em solo não saturado. Em relação à análise de sensibilidade, o modelo apresentou maior sensibilidade aos parâmetros de fluxo de entrada e ao parâmetro n da equação de retenção, indicando que tais parâmetros de entrada, necessitam ser determinados com maior precisão. Houve baixa sensibilidade ao parâmetro \"alfa\" da curva de retenção. === The soil moisture distribution for different soil types under different flow regimes can be used to optimize the design of an irrigation system, water manage water in the root zone, or increase the precision of the application of fertilizers in the rhizosphere. Therefore, the objective of this research was to develop a two-dimensional numerical model capable of simulating the soil moisture distribution in soil profile for different irrigation systems. The model is a finite difference solution of the Richards equation solution, in which the Picard iteration method is used to ensure the conservation of mass. The numerical model also takes water extraction by plants and water evaporation on the soil surface into account. The mathematical model was developed in the Department of Biosystems Engineering, College of Agriculture \"Luiz de Queiroz\", ESALQ / USP. It was coded in Fortran 90, and structured to allow interactive input of a) physical and hydraulic properties soil profile, and b) type of irrigation, irrigation application time, and simulation. The simulation of water extraction by the plants is optional. And to evaluate its performance a sensitivity analysis was done. Data obtained in an experiment carried out in the older Rural Engineering Department, College of Agriculture \"Luiz de Queiroz\" by Rivera (2004), were used to validate the model. Model performance was evaluated using the Willmott index of agreement; efficiency coefficient, root mean square error, average error and the maximum absolute error. Sensitivity analyses were performed on the density of water flow, saturated hydraulic conductivity, and on alpha and n (setting parameters of the soil water retention curve). From the results obtained, the model satisfactory simulated the measured soil moisture distributions in the experimental soil moisture profiles. The model was most sensitive to variations in the parameters of the input stream and to n, a parameter in the soil water retention equation (van Genuchten model, 1980). This results is indicative that these input parameters need to be determined with greater accuracy. The model output was not very sensitive to alpha, a shape parameter in the soil water retention curve.
|
|---|
