92 research outputs found

    Estimation of water content profiles by inverse analysis of TDR waveforms: application to infer soil hydraulic properties

    Get PDF
    La medida de las propiedades hidráulicas del suelo, curva de retención (θ(ψ)) y conductividad hidráulica (K) tiene una importancia fundamental para la simulación de procesos hidrológicos. La técnica de Reflectometría de Dominio Temporal (TDR) es una herramienta ampliamente utilizada para la medida no destructiva de contenido volumétrico de agua en el suelo (θ) y de conductividad eléctrica (σ). El objetivo de este trabajo es desarrollar una nueva metodología basada en el uso de la técnica TDR para estimar las propiedades hidráulicas del suelo (α, n and K) por análisis inverso de la dinámica de los perfiles de humedad (WCPs ) durante un proceso de infiltración de agua. Los WCPs se estiman a partir del análisis inverso de ondas TDR empleando un modelo físico de propagación electromagnética. Posteriormente, los parámetros α, n and K se calculan empleando una interfaz HYDRUS-1D-Matlab por medio del análisis inverso de los TDR-WCPs. Esta interfaz calcula los parámetros hidráulicos a partir del mejor ajuste entre los TDR-WCP registrados y los simulados por HYDRUS-1D. Para este fin, se emplea el método de optimización de fuerza bruta, el cual permite barrer un rango amplio de parámetros hidráulicos. El método fue probado en tres medios porosos distintos (tierra franca tamizada a 2 mm, arena y microesferas de vidrio) durante un proceso de infiltración. Las propiedades hidráulicas estimadas con este método se compararon con aquellas medidas en los mismos medios porosos empleando técnicas convencionales de laboratorio: cámaras de presión-TDR y mini-infiltrómetro de disco. Aunque se obtuvieron resultados satisfactorios para la mediad de K, los resultados obtenidos para la medida de n y α fueron imprecisos. Estas discrepancias se pueden atribuir a las siguientes causas (i) el uso de una función unimodal en lugar de una bimodal en HYDRUS-1D; (ii) el fenómeno de histéresis del suelo; (iii) incertidumbres en la función “puente” empleada para estimar WCP a partir del modelo de simulación de ondas TDR. Este método necesita nuevos esfuerzos para mejorar su precisión y así poder probarse en muestras de suelo inalterado

    Using a sprinkler infiltrometer and the GAML model to predict moving sprinkler performance in the field

    Get PDF
    The Green-Ampt (1911) equation has been regarded as one of the foremost infiltration models. Mein and Larson’s (1971) work extended its use to modelling infiltration under rainfall conditions, known as the GAML model, and Chu (1987) further extended its use to time-varying application rates such as occurs under moving sprinkler systems. However, Chu only demonstrated the efficacy of his work using simple, idealised application patterns that are not seen in the real world. This project, then, sought to extend Chu’s work by testing it in the field using real sprinkler performance data. Sprinkler performance data, using Nelson brand centre-pivot S3000 sprinkler heads, was collected for the project in the hydraulics laboratory at USQ, Toowoomba. A sprinkler infiltrometer was used in the field to determine modified GAML model parameters, per Chu (1986). A computer program written in Matlab, based on the graphical methods of Chu (1987), used the laboratory sprinkler data and the modified GAML parameters to make a prediction of the runoff that would be generated from a specified time-varying application rate. A mobile sprinkler rig was constructed to deliver the time-varying application rate of water in the field. A new concept for a sprinkler infiltrometer was tested in the course of the project. Initial work in the laboratory appeared promising but the concept failed to meet expectations in the field, principally due to wind interference. Consequently a small droplet-forming sprinkler infiltrometer was constructed and used for all of the field testing. The process of determining the modified GAML model parameters was reasonably successful. However, predictions of runoff by the computer model were consistently far larger than that measured. This was believed to be for two key reasons. Firstly, whilst Chu’s (1987) method worked well for simple application patterns, it appeared to struggle with real data and so the predicted runoff by the computer program could only be regarded with suspicion. Secondly, there were significant difficulties collecting all of the runoff from the soil plots in the field. These two reasons for the disparity between predicted and measured runoff meant that this project was not able to conclusively affirm or reject Chu’s (1987) method for applying the GAML model as being suitable for use with real moving sprinkler systems in the field

    Rainfall-Runoff Processes in Tropical Urban Environments

    Get PDF
    Ph.DDOCTOR OF PHILOSOPH

    Upscaling of water flow and mass transport in a tropical soil: numerical, laboratory and field studies

    Full text link
    Los modelos numéricos son herramientas fundamentales para realizar predicciones de muchos problemas enfrentados por ingenieros geotécnicos y geoambientales. Sin embargo, para que estos modelos puedan realizar predicciones confiables, los parámetros de entrada del modelo deben ser estimados considerando el efecto escala. En este contexto, esta tesis se concentra en las reglas del cambio de escala de los parámetros de flujo y transporte de masa en un suelo tropical a través de estudios numéricos, de laboratorio y de campo. Esta está organizada en cuatro partes. Primero, la heterogeneidad, correlación y correlación cruzada entre los parámetros de transporte de solutos (dispersividad, ¿, y coeficiente de partición, Kd) y las propiedades del suelo fueron estudiadas en detalle. En esta parte fue verificado que la conductividad hidráulica (K) y los parámetros de transporte de solutos son altamente heterogéneos, mientras que las propiedades del suelo no lo son. La correlación espacial de ¿ y K con variables estadísticamente significativas fue estudiada. Este resultado probablemente podrá mejorar la estimación en casos de estudios de pequeña escala debido a que solo fue observada correlaciones de hasta 2,5 m. Este estudio fue un primer intento de evaluar la variación espacial en el coeficiente de correlación de los parámetros de transporte de un soluto reactivo y de un no reactivo, indicando las variables más relevantes y aquella que debería ser incluida en estudios futuros. En la segunda parte, el efecto escala en K, dispersividad y coeficiente de partición de potasio y clorito fue estudiado experimentalmente a través de experimentos de laboratorio y de campo. El objetivo de esta parte fue contribuir a la discusión sobre el efecto escala en K, ¿ y Kd, y entender como estos parámetros se comportan con el cambio de escala de medición. La dispersividad tiende a aumentar con la altura de la muestra de manera exponencial. El coeficiente de partición tiende a aumentar con la altura, el diámetro y el volumen de la muestra. Estas diferencias encontradas en los parámetros de acuerdo con la escala de medición deben ser considerados cuando estos valores sean usados posteriormente como datos de entrada de modelos numéricos; de otra manera, las respuestas pueden ser malinterpretadas. Tercero, análisis estocásticos tridimensionales de cambio de escala de la conductividad hidráulica fueron realizados usando los métodos de promedios simples y de Laplace con piel para una variedad de tamaños de bloques usando mediciones reales de K. En esta parte son demostrados los errores que pueden ser introducidos al usar métodos determinísticos de cambio de escala usando promedios simples de las mediciones de K sin llevar en consideración la correlación espacial. La aplicación muestra que la heterogeneidad de K puede ser incorporada en la práctica diaria del modelador geotécnico. Los aspectos que considerar durante un proceso de cambio de escala también son discutidos. Finalmente, la dependencia del exponente de la norma-p como función del tamaño del bloque fue analizada. En la última parte, una aplicación de cambio de escala estocástico del coeficiente de dispersión hidrodinámica D y del factor de retardo R fue realizada usando datos reales con el objetivo de reducir la falta de casos de investigación experimental de cambio de escala de parámetros de transporte de solutos reactivos. El cambio de escala de D fue realizado usando el método de macrodispersión. El método de promedio simple de norma-p fue usado para realizar el cambio de escala de R. Una buena propagación de incertidumbres fue alcanzada. Métodos simples de cambio de escala pueden ser introducidos en la práctica del modelaje usando programas comerciales de transporte y conseguir reproducir el transporte en escala gruesa, pero puede requerir correcciones con el objetivo de reducir el efecto de suavizado de la heterogeneidad causado por elNumerical models are becoming fundamental tools to predict a range of complex problems faced by geotechnical and geo-environmental engineers. However, to render the model reliable for future predictions, the model input parameters must be determined with consideration of the scale effects. In this context, this thesis focuses on upscaling of water flow and mass transport in a tropical soil by means of numerical, laboratory and field studies. This thesis is organized in four parts. First, the heterogeneity, correlation and cross-correlation between solute transport parameters (dispersivity, ¿, and partition coefficient, Kd) and soil properties were studied in detail. In this part, it was verified that the hydraulic conductivity (K) and solute transport parameters are highly heterogeneous, while soil properties not. Spatial correlation of ¿, K, and statistically significant variables were studied, and it would probably improve the estimation only in a small-scale study, since the spatial correlation were only observed up to 2.5 m. This study was a first attempt to evaluate the spatial variation in the correlation coefficient of transport parameters of a reactive and a nonreactive solute, indicating the more relevant variables and the one that should be included in future studies. In the second part, scale effect on K, dispersivity and partition coefficient of potassium and chloride is studied experimentally by means of laboratory and field experiments. The purpose of was to contribute to the discussion about scale effects on K, ¿ and Kd and understanding how these parameters behave with the change in the scale of measurement. Results shows that K increases with scale, regardless of the method of measurement. Dispersivity trends to increases exponentially with the sample height. Partition coefficient, tend to increase with sample length, diameter and volume. These differences in the parameters according to the scale of measurement must be considered when these observations are later used as input to numerical models, otherwise the responses can be misrepresented. Third, stochastic analysis of three-dimensional hydraulic conductivity upscaling was performed using a simple average and the Laplacian-with-skin methods for a variety of block sizes using real K measurements. In this part it was demonstrated the errors that can be introduced by using a deterministic upscaling using simple averages of the measured K without accounting for the spatial correlation. The application shows that K heterogeneity can be incorporated in the daily practice of the geotechnical modeler. The aspects to consider when performing the upscaling were also discussed. Finally, the dependence of the exponent of the p-norm as a function of the block size was analyzed. In the last part, an application of stochastic upscaling of hydrodynamic dispersion coefficient (D) and retardation factor (R) was performed using real data aiming to reduce the lack in experimental upscaling of reactive solute transport research. Upscaling of D was done using macrodispersion method. Simple average method based on p-norm was used to perform R upscaling. A good propagation of the uncertainties was achieved. Simple upscaling methods can be incorporated to the modeling practice using commercial transport codes and properly reproduce de transport at coarse scale but may require corrections to reduce smoothing of the heterogeneity caused by the upscaling procedure.Els models numèrics s'estan constituint en eines fonamentals per a realitzar prediccions d'una àmplia gamma de problemes enfrontats per enginyers geotècnics i geoambientales. No obstant açò, perquè aquests models puguen realitzar prediccions fiables, els paràmetres d'entrada del model han de considerar l'efecte escala. En aquest context, aquesta tesi es concentra en les regles del canvi d'escala dels paràmetres de flux i transport de massa en un sòl tropical a través d'estudis numèrics, de laboratori i de camp. Aquesta tesi està organitzada en quatre parts. Primer, l'heterogeneïtat, correlació i correlació creuada entre els paràmetres de transport de soluts (dispersivitat, ¿, i coeficient de partició, Kd) i les propietats del sòl van ser estudiades detalladament. En aquesta part va ser verificat que la conductivitat hidràulica (K) i els paràmetres de transport de soluts són altament heterogenis, mentre que les propietats del sòl no ho són. La correlació espacial de ¿ i K amb variables estadísticament significatives va ser estudiada. Aquest resultat probablement podrà millorar l'estimació en casos d'estudis de xicoteta escala a causa que solament va ser observada correlacions de fins a 2,5 m. Aquest estudi va ser un primer intent d'avaluar la variació espacial en el coeficient de correlació dels paràmetres de transport d'un solut reactiu i d'un no reactiu, indicant les variables més rellevants i aquelles que haurien de ser inclosas en estudis futurs. En la segona part, l'efecte escala en K, dispersivitat i coeficient de partició de potassi i clorito va ser estudiat experimentalment a través d'experiments de laboratori i de camp. L'objectiu d'aquesta part va ser contribuir a la discussió sobre l'efecte escala en K, ¿ i Kd, i entendre com aquests paràmetres es comporten amb el canvi d'escala de mesurament. La dispersivitat tendeix a augmentar amb l'altura de la mostra, és a dir, amb la longitud del transport, de manera exponencial. El coeficient de partició tendeix a augmentar amb l'altura, el diàmetre i el volum de la mostra. Aquestes diferències en els paràmetres d'acord amb l'escala de mesurament han de ser considerats quan aquests valors siguen usats posteriorment com a dades d'entrada de models numèrics; d'una altra manera, les respostes poden ser malament interpretades. Tercer, anàlisis estocàstiques tridimensionals de canvi d'escala de la conductivitat hidràulica van ser realitzats usant els mètodes de mitjanes simples i de Laplace amb pell per a una varietat de grandàries de blocs usant mesuraments reals de K. En aquesta part són demostrats els errors que poden ser introduïts en usar mètodes determinístics de canvi d'escala usant mitjanes simples dels mesuraments de K sense tindre en consideració la correlació espacial. L'aplicació mostra que l'heterogeneïtat de K pot ser incorporada en la pràctica diària del modelador geotècnic. Els aspectes a considerar durant un procés de canvi d'escala també són discutits. Finalment, la dependència de l'exponent de la norma-p com a funció de la grandària del bloc va ser analitzada. En l'última part, una aplicació de canvi d'escala estocàstic del coeficient de dispersió hidrodinámica D i del factor de retard R va ser realitzada usant dades reals amb l'objectiu de reduir la falta de casos de recerca experimental de canvi d'escala de paràmetres de transport de soluts reactius. El canvi d'escala de D va ser realitzat usant el mètode de macrodispersió. El mètode de mitjana simple de norma-p va ser usat per a realitzar el canvi d'escala de R. Una bona propagació d'incerteses va ser aconseguida. Mètodes simples de canvi d'escala poden ser introduïts en la pràctica de la modelació usant programes comercials de transport i aconseguir reproduir el transport en escala gruixuda, però pot requerir correccions amb l'objectiu de reduir l'efecte de suavitzat de l'heterogeneïtat causat pel procés de canvi d'escala.Almeida De Godoy, V. (2018). Upscaling of water flow and mass transport in a tropical soil: numerical, laboratory and field studies [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/102405TESI

    Comparing the Hydraulic Properties of Forested and Grassed Soils on an Experimental Hillslope in a Mediterranean Environment

    Get PDF
    AbstractThis experimental research compares the physical and hydraulic properties of two adjacent soils, one covered with a native forest of Mediterranean maquis, and the other with spontaneous grass. The latter replaced the previous natural forest. The aim is to quantify the significant differences in the soil properties caused by the removal of the natural vegetation. Although the soil texture was similar in the different land uses, the soil under the forest had a higher organic matter content, a lower apparent density and a higher water content at saturation than the grassed soil. The analysis of the water retention characteristics indicated that the retained water content of the forest soil exceeded that of the grassed soil in the range from saturation to -50cm of water tension. This suggests that changing the land use altered the soil pore structure within this range. The hydraulic conductivity of the forest soil exceeded that of the grassed soil at water tensions of -10, -5 and -3cm. Conversely the hydraulic conductivity of the grassed soil was similar to that of the forest soil at -1cm of water tension and at saturation. This result was probably due to the hydraulic activation of the desiccation cracks in the grassed soil. This increased the amount of infiltrated water in saturated and near-saturated soil conditions.This work shows that changes in land use have an unfavorable impact on the physical and hydraulic properties of the soil. Soil covered with grass is more vulnerability to water erosion than that under forest, and there is likely to be general worsening of flow regimes

    Development and analysis of the Soil Water Infiltration Global database

    Get PDF
    In this paper, we present and analyze a novel global database of soil infiltration measurements, the Soil Water Infiltration Global (SWIG) database. In total, 5023 infiltration curves were collected across all continents in the SWIG database. These data were either provided and quality checked by the scientists who performed the experiments or they were digitized from published articles. Data from 54 different countries were included in the database with major contributions from Iran, China, and the USA. In addition to its extensive geographical coverage, the collected infiltration curves cover research from 1976 to late 2017. Basic information on measurement location and method, soil properties, and land use was gathered along with the infiltration data, making the database valuable for the development of pedotransfer functions (PTFs) for estimating soil hydraulic properties, for the evaluation of infiltration measurement methods, and for developing and validating infiltration models. Soil textural information (clay, silt, and sand content) is available for 3842 out of 5023 infiltration measurements ( ∼ 76%) covering nearly all soil USDA textural classes except for the sandy clay and silt classes. Information on land use is available for 76% of the experimental sites with agricultural land use as the dominant type ( ∼ 40%). We are convinced that the SWIG database will allow for a better parameterization of the infiltration process in land surface models and for testing infiltration models. All collected data and related soil characteristics are provided online in *.xlsx and *.csv formats for reference, and we add a disclaimer that the database is for public domain use only and can be copied freely by referencing it. Supplementary data are available at https://doi.org/10.1594/PANGAEA.885492 (Rahmati et al., 2018). Data quality assessment is strongly advised prior to any use of this database. Finally, we would like to encourage scientists to extend and update the SWIG database by uploading new data to it

    Development and analysis of the Soil Water Infiltration Global database.

    Get PDF
    In this paper, we present and analyze a novel global database of soil infiltration measurements, the Soil Water Infiltration Global (SWIG) database. In total, 5023 infiltration curves were collected across all continents in the SWIG database. These data were either provided and quality checked by the scientists who performed the experiments or they were digitized from published articles. Data from 54 different countries were included in the database with major contributions from Iran, China, and the USA. In addition to its extensive geographical coverage, the collected infiltration curves cover research from 1976 to late 2017. Basic information on measurement location and method, soil properties, and land use was gathered along with the infiltration data, making the database valuable for the development of pedotransfer functions (PTFs) for estimating soil hydraulic properties, for the evaluation of infiltration measurement methods, and for developing and validating infiltration models. Soil textural information (clay, silt, and sand content) is available for 3842 out of 5023 infiltration measurements (~76%) covering nearly all soil USDA textural classes except for the sandy clay and silt classes. Information on land use is available for 76% of the experimental sites with agricultural land use as the dominant type (~40%). We are convinced that the SWIG database will allow for a better parameterization of the infiltration process in land surface models and for testing infiltration models. All collected data and related soil characteristics are provided online in *.xlsx and *.csv formats for reference, and we add a disclaimer that the database is for public domain use only and can be copied freely by referencing it. Supplementary data are available at https://doi.org/10.1594/PANGAEA.885492 (Rahmati et al., 2018). Data quality assessment is strongly advised prior to any use of this database. Finally, we would like to encourage scientists to extend and update the SWIG database by uploading new data to it

    Prediction of water retention of soils from the humid tropics by the nonparametric k-nearest neighbor approach

    Get PDF
    Nonparametric approaches such as the k-nearest neighbor (k-NN) approach are considered attractive for pedotransfer modeling in hydrology; however, they have not been applied to predict water retention of highly weathered soils in the humid tropics. Therefore, the objectives of this study were: to apply the k-NN approach to predict soil water retention in a humid tropical region; to test its ability to predict soil water content at eight different matric potentials; to test the benefit of using more input attributes than most previous studies and their combinations; to discuss the importance of particular input attributes in the prediction of soil water retention at low, intermediate, and high matric potentials; and to compare this approach with two published tropical pedotransfer functions (PTFs) based on multiple linear regression (MLR). The overall estimation error ranges generated by the k-NN approach were statistically different but comparable to the two examined MLR PTFs. When the best combination of input variables (sand + silt + clay + bulk density + cation exchange capacity) was used, the overall error was remarkably low: 0.0360 to 0.0390 m(3) m(-3) in the dry and very wet ranges and 0.0490 to 0.0510 m(3) m(-3) in the intermediate range (i.e., -3 to -50 kPa) of the soil water retention curve. This k-NN variant can be considered as a competitive alternative to more classical, equation-based PTFs due to the accuracy of the water retention estimation and, as an added benefit, its flexibility to incorporate new data without the need to redevelop new equations. This is highly beneficial in developing countries where soil databases for agricultural planning are at present sparse, though slowly developing
    corecore