Determination of hydraulic properties of a tropical soil of Hawaii using column experiments and inverse modeling

A method for determining soil hydraulic properties of a weathered tropical soil (Oxisol) using a medium-sized column with undisturbed soil is presented. The method was used to determine fitting parameters of the water retention curve and hydraulic conductivity functions of a soil column in support of a pesticide leaching study. The soil column was extracted from a continuously-used research plot in Central Oahu (Hawaii, USA) and its internal structure was examined by computed tomography. The experiment was based on tension infiltration into the soil column with free outflow at the lower end. Water flow through the soil core was mathematically modeled using a computer code that numerically solves the one-dimensional Richards equation. Measured soil hydraulic parameters were used for direct simulation, and the retention and soil hydraulic parameters were estimated by inverse modeling. The inverse modeling produced very good agreement between model outputs and measured flux and pressure head data for the relatively homogeneous column. The moisture content at a given pressure from the retention curve measured directly in small soil samples was lower than that obtained through parameter optimization based on experiments using a medium-sized undisturbed soil column.

Hydraulic conductivity and water retention in leptosols-regosols and saprolite derived from sandstone, Brazil

Leptosols and Regosols are soils with a series of restrictions for use, mainly related to the effective depth, which have been poorly studied in Brazil. These soils, when derived from sedimentary rocks should be treated with particular care to avoid environmental damage such as aquifer contamination. The purpose of this study was to verify the behavior of hydraulic conductivity and water retention capacity in profiles of Leptosols and Regosols derived from sandstone of the Caturrita formation in Rio Grande do Sul state. The morphology, particle size distribution, porosity, soil density (Ds), saturated hydraulic conductivity (Ks), basic water infiltration in the field (BI) and water retention were determined in soil and saprolite samples of six soil profiles. High Ds, low macroporosity and high microporosity were observed in the profiles, resulting in a low Ks and BI, even under conditions of sandy texture and a highly fractured saprolite layer. The variation coefficients of data of Ks and BI were high among the studied profiles and between replications of a same profile. Water retention of the studied soils was higher in Cr layers than in the A horizons and the volume of plant-available water greater and variable among A horizons and Cr layers.

Hydraulic conductivity and soil-sewage sludge interactions

One of the main problems faced by humanity is pollution caused by residues resulting from the production and use of goods, e.g, sewage sludge. Among the various alternatives for its disposal, the agricultural use seems promising. The purpose of this study was to evaluate the hydraulic conductivity and interaction of soil with sandy-silty texture, classified as Spodosols, from the Experimental Station Itapirema - IPA, in Goiana, state of Pernambuco, in mixtures with sewage sludge from the Mangueira Sewage Treatment Station, in the city of Recife, Pernambuco at rates of 25, 50 and 75 Mg ha-1. Tests were conducted to let water percolate the natural saturated soil and soil-sludge mixtures to characterize their physical, chemical, and microstructural properties as well as hydraulic conductivity. Statistical data analysis showed that the presence of sewage sludge in soils leads to an increase of the < 0.005 mm fraction, reduction in real specific weight and variation in optimum moisture content from 11.60 to 12.90 % and apparent specific dry weight from 17.10 and 17.50 kN m-3. In the sludge-soil mixture, the quartz grains were covered by sludge and filling of the empty soil macropores between grains. There were changes in the chemical characteristics of soil and effluent due to sewage sludge addition and a small decrease in hydraulic conductivity. The results indicate the possibility that soil acidity influenced the concentrations of the elements found in the leachate, showing higher levels at higher sludge doses. It can be concluded that the leaching degree of potentially toxic elements from the sewage sludge treatments does not harm the environment.

Soil thermal diffusivity estimated from data of soil temperature and single soil component properties

Under field conditions, thermal diffusivity can be estimated from soil temperature data but also from the properties of soil components together with their spatial organization. We aimed to determine soil thermal diffusivity from half-hourly temperature measurements in a Rhodic Kanhapludalf, using three calculation procedures (the amplitude ratio, phase lag and Seemann procedures), as well as from soil component properties, for a comparison of procedures and methods. To determine thermal conductivity for short wave periods (one day), the phase lag method was more reliable than the amplitude ratio or the Seemann method, especially in deeper layers, where temperature variations are small. The phase lag method resulted in coherent values of thermal diffusivity. The method using properties of single soil components with the values of thermal conductivity for sandstone and kaolinite resulted in thermal diffusivity values of the same order. In the observed water content range (0.26-0.34 m³ m-3), the average thermal diffusivity was 0.034 m² d-1 in the top layer (0.05-0.15 m) and 0.027 m² d-1 in the subsurface layer (0.15-0.30 m).

Stabilization of soil hydraulic properties under a long term no-till system

The area under the no-tillage system (NT) has been increasing over the last few years. Some authors indicate that stabilization of soil physical properties is reached after some years under NT while other authors debate this. The objective of this study was to determine the effect of the last crop in the rotation sequence (1st year: maize, 2nd year: soybean, 3rd year: wheat/soybean) on soil pore configuration and hydraulic properties in two different soils (site 1: loam, site 2: sandy loam) from the Argentinean Pampas region under long-term NT treatments in order to determine if stabilization of soil physical properties is reached apart from a specific time in the crop sequence. In addition, we compared two procedures for evaluating water-conducting macroporosities, and evaluated the efficiency of the pedotransfer function ROSETTA in estimating the parameters of the van Genuchten-Mualem (VGM) model in these soils. Soil pore configuration and hydraulic properties were not stable and changed according to the crop sequence and the last crop grown in both sites. For both sites, saturated hydraulic conductivity, K0, water-conducting macroporosity, εma, and flow-weighted mean pore radius, R0ma, increased from the 1st to the 2nd year of the crop sequence, and this was attributed to the creation of water-conducting macropores by the maize roots. The VGM model adequately described the water retention curve (WRC) for these soils, but not the hydraulic conductivity (K) vs tension (h) curve. The ROSETTA function failed in the estimation of these parameters. In summary, mean values of K0 ranged from 0.74 to 3.88 cm h-1. In studies on NT effects on soil physical properties, the crop effect must be considered.

Physical-hydraulic properties of a sandy loam typic paleudalf soil under organic cultivation of 'montenegrina' mandarin (Citrus deliciosa Tenore)¹

Citrus plants are the most important fruit species in the world, with emphasis to oranges, mandarins and lemons. In Rio Grande do Sul, Brazil, most fruit production is found on small properties under organic cultivation. Soil compaction is one of the factors limiting production and due to the fixed row placement of this crop, compaction can arise in various manners in the interrows of the orchard. The aim of this study was to evaluate soil physical properties and water infiltration capacity in response to interrow management in an orchard of mandarin (Citrus deliciosa Tenore 'Montenegrina') under organic cultivation. Interrow management was performed through harrowing, logs in em "V", mowing, and cutting/knocking down plants with a knife roller. Soil physical properties were evaluated in the wheel tracks of the tractor (WT), between the wheel tracks (BWT), and in the area under the line projection of the canopy (CLP), with undisturbed soil samples collected in the 0.00-0.15, 0.15-0.30, 0.30-0.45, and 0.45-0.60 m layers, with four replicates. The soil water infiltration test was performed using the concentric cylinder method, with a maximum time of 90 min for each test. In general, soil analysis showed a variation in the physical-hydraulic properties of the Argissolo Vermelho-Amarelo distrófico arênico (sandy loam Typic Paleudalf) in the three sampling sites in all layers, regardless of the management procedure in the interrows. Machinery traffic leads to heterogeneity in the soil physical-hydraulic properties in the interrows of the orchard. Soil porosity and bulk density are affected especially in the wheel tracks of the tractor (WT), which causes a reduction in the constant rate of infiltration and in the accumulated infiltration of water in this sampling site. The use of the disk harrow and mower leads to greater harmful effects on the soil, which can interfere with mandarin production.

Multiplan Spreadsheets for Highway and Hydraulic Engineering Computations in County Engineering Offices, HR-284, 1987

Multiplan spreadsheet solutions were developed for a set of hydraulic and highway engineering computations of common interest to county engineers. These include earthwork, vertical and horizontal curves, staking superelevated curves and sign inventories for highways. The hydraulic applications were ditch flow, runoff, culvert size and stage discharge.

Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.

Bayesian-MCMC inversion of time-lapse geophysical data to characterize vadose zone hydraulic behaviour

Geophysical methods have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, time-lapse geophysical data, when coupled with a hydrological model and inverted stochastically, may allow for the effective estimation of subsurface hydraulic parameters and their corresponding uncertainties. In this study, we use a Bayesian Markov-chain-Monte-Carlo (MCMC) inversion approach to investigate how much information regarding vadose zone hydraulic properties can be retrieved from time-lapse crosshole GPR data collected at the Arrenaes field site in Denmark during a forced infiltration experiment.

Impact of controlled changes in grain size and pore space characteristics on the hydraulic conductivity and spectral induced polarization response of "proxies" of saturated alluvial sediments

Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial unconsolidated sedimentary deposits, which host the bulk of the world's readily accessible groundwater resources. Here, we present the results of laboratory SIP measurements on industrial-grade, saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel, which can be regarded as proxies for widespread alluvial deposits. We altered the pore space characteristics by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples. In general, the results indicate a clear connection between the SIP response and the granulometric as well as pore space characteristics. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the entire range of considered grain sizes. The results do, however, also indicate that the detailed nature of these relations depends strongly on variations in the pore space characteristics, such as, for example, the degree of compaction. The results of this study underline the complexity of the origin of the SIP signal as well as the difficulty to relate it to a single structural factor of a studied sample, and hence raise some fundamental questions with regard to the practical use of SIP measurements as site- and/or sample-independent predictors of the hydraulic conductivity.

Automation of DEM Cutting for Hydrologic/Hydraulic Modeling, TR-631, 2015

Hydrologic analysis is a critical part of transportation design because it helps ensure that hydraulic structures are able to accommodate the flow regimes they are likely to see. This analysis is currently conducted using computer simulations of water flow patterns, and continuing developments in elevation survey techniques result in higher and higher resolution surveys. Current survey techniques now resolve many natural and anthropogenic features that were not practical to map and, thus, require new methods for dealing with depressions and flow discontinuities. A method for depressional analysis is proposed that uses the fact that most anthropogenically constructed embankments are roughly more symmetrical with greater slopes than natural depressions. An enforcement method for draining depressions is then analyzed on those depressions that should be drained. This procedure has been evaluated on a small watershed in central Iowa, Walnut Creek of the South Skunk River, HUC12 # 070801050901, and was found to accurately identify 88 of 92 drained depressions and place enforcements within two pixels, although the method often tries to drain prairie pothole depressions that are bisected by anthropogenic features.

Scale Effects in Hydraulic Model Tests of Rock Protected Structures, HR-119, 1970

A laboratory investigation was undertaken to determine the limiting model Reynolds number above which the scour behavior of rock protected structures can be reproduced in hydraulic models scaled according to the Froude criterion. A submerged jet was passed over an initially full scour pocket containing uniform glass spheres and the rate of scour was measured as a function of time. The dimensions of the scour pocket and jet and the particle diameters were varied as needed to maintain strict geometric similarity. For each of two different Froude numbers the Reynolds number was varied over a wide range. The normalized scour rate was found to be practically independent of the Reynolds number, R, (based on the jet velocity and particle diameter) at values of R above about 2.5 x 10^3, and to decrease with Rat smaller values. A grid placed in the jet was found to have a very strong effect on the scour rate. In an attempt to explain the effect of R on the scour behavior, turbulent pressure and velocity fluctuations were measured in air flows and water flows, respectively, over rigid scour pockets having the same geometry as those formed in the scour experiments. The normalized spectra of the fluctuations were found to be nearly independent of R, but the flow pattern was found to be very sensitive to the inlet condition, the jet deflecting upward or downward in a not wholly explainable manner. This indicates that scour behavior can be modeled only if the approach flow is also accurately modeled.

A non-linear Bayesian approach for upscaling local-scale hydraulic conductivity measurements based on local- and regional-scale geophysical data

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed an upscaling procedure based on a Bayesian sequential simulation approach. This method is then applied to the stochastic integration of low-resolution, regional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this upscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.