Spatial variability of satured soil hydraulic conductivity in the region of Araguaia River - Brazil.

This study evaluates the spatial variability of saturated hydraulic conductivity in the soil in an area of 51,850 ha at the headwaters of the Araguaia River MT/GO. This area is highly vulnerable because it is a location of recharging through natural water infiltration of the Guarani Aquifer System and an area of intense increases in agriculture since its adoption by growers in the last 30 years. Soil samples were collected at 383 points, geographically located by GPS. The samples were collected from depths of 0 - 20 cm and 60 - 80 cm. Exploratory statistics and box-plot were used in the descriptive analysis and semivariogram were constructed to determine the spatial model. The exploratory analysis showed that the mean hydraulic conductivity in the superficial layer was less than at the level of 60-80 cm; however, the greatest variability evaluated with a coefficient of variation also was from this layer. Data tended towards a normal distribution. These results can be explained by the greater soil compaction in the superficial layer. The semivariogram models, adjusted for the two layers, were exponential and demonstrated moderate and strong dependence, with ranges of 5000 and 3000 utm respectively. It was concluded that soil use is influencing the spatial distribution model of the hydraulic conductivity in the region.

Spatial variability of satured soil hydraulic conductivity in the region of Araguaia River - Brazil.

2008

Hydrostratigraphic analysis of the MADE site with full-resolution GPR and direct-push hydraulic profiling

Full-resolution 3D Ground-Penetrating Radar (GPR) data were combined with high-resolution hydraulic conductivity (K) data from vertical Direct-Push (DP) profiles to characterize a portion of the highly heterogeneous MAcro Dispersion Experiment (MADE) site. This is an important first step to better understand the influence of aquifer heterogeneities on observed anomalous transport. Statistical evaluation of DP data indicates non-normal distributions that have much higher similarity within each GPR facies than between facies. The analysis of GPR and DP data provides high-resolution estimates of the 3D geometry of hydrostratigraphic zones, which can then be populated with stochastic K fields. The lack of such estimates has been a significant limitation for testing and parameterizing a range of novel transport theories at sites where the traditional advection-dispersion model has proven inadequate.

Radium and barium removal through blending hydraulic fracturing fluids with acid mine drainage.

Wastewaters generated during hydraulic fracturing of the Marcellus Shale typically contain high concentrations of salts, naturally occurring radioactive material (NORM), and metals, such as barium, that pose environmental and public health risks upon inadequate treatment and disposal. In addition, fresh water scarcity in dry regions or during periods of drought could limit shale gas development. This paper explores the possibility of using alternative water sources and their impact on NORM levels through blending acid mine drainage (AMD) effluent with recycled hydraulic fracturing flowback fluids (HFFFs). We conducted a series of laboratory experiments in which the chemistry and NORM of different mix proportions of AMD and HFFF were examined after reacting for 48 h. The experimental data combined with geochemical modeling and X-ray diffraction analysis suggest that several ions, including sulfate, iron, barium, strontium, and a large portion of radium (60-100%), precipitated into newly formed solids composed mainly of Sr barite within the first ∼ 10 h of mixing. The results imply that blending AMD and HFFF could be an effective management practice for both remediation of the high NORM in the Marcellus HFFF wastewater and beneficial utilization of AMD that is currently contaminating waterways in northeastern U.S.A.

Modeling plant transpiration under limited soil water: comparison of different plant and soil hydraulic parameterizations and preliminary implications for their use in land surface models

Accurate estimates of how soil water stress affects plant transpiration are crucial for reliable land surface model (LSM) predictions. Current LSMs generally use a water stress factor, β, dependent on soil moisture content, θ, that ranges linearly between β = 1 for unstressed vegetation and β = 0 when wilting point is reached. This paper explores the feasibility of replacing the current approach with equations that use soil water potential as their independent variable, or with a set of equations that involve hydraulic and chemical signaling, thereby ensuring feedbacks between the entire soil–root–xylem–leaf system. A comparison with the original linear θ-based water stress parameterization, and with its improved curvi-linear version, was conducted. Assessment of model suitability was focused on their ability to simulate the correct (as derived from experimental data) curve shape of relative transpiration versus fraction of transpirable soil water. We used model sensitivity analyses under progressive soil drying conditions, employing two commonly used approaches to calculate water retention and hydraulic conductivity curves. Furthermore, for each of these hydraulic parameterizations we used two different parameter sets, for 3 soil texture types; a total of 12 soil hydraulic permutations. Results showed that the resulting transpiration reduction functions (TRFs) varied considerably among the models. The fact that soil hydraulic conductivity played a major role in the model that involved hydraulic and chemical signaling led to unrealistic values of β, and hence TRF, for many soil hydraulic parameter sets. However, this model is much better equipped to simulate the behavior of different plant species. Based on these findings, we only recommend implementation of this approach into LSMs if great care with choice of soil hydraulic parameters is taken

Bridging hydraulic diffusivity from aquifer to particle-size scale: a study on loess sediments from southwest Hungary

In situ megascale hydraulic diffusivities (D) of a confined loess aquifer were estimated at various scales (10 <= L <= 1500 m) by a finite difference model, and laboratory microscale diffusivities of a loess sample by empirical formulas. A scatter plot reveals that D fits to a single power function of L, providing that microscale diffusivities are assigned to L = 1 m and that differences in diffusivity observed between micro- and megascales are assigned to medium heterogeneity appraised by variations in the curvature and slope of natural hydraulic head waves propagating through the aquifer. Subsequently, a general power relationship between D and L is defined where the base and exponent terms stand for the aquifer storage capability under a confined regime of flow, for the microscale hydraulic conductivity and specific yield of loess, and for the changes in curvature and slope of hydraulic head waves relative to values defined at unit scale.[GRAPHICS]Editor Z.W. Kundzewicz

Comparison among monitoring strategies to assess water flow dynamic and soil hydraulic properties in agricultural soils

Irrigated agriculture is usually performed in semi-arid regions despite scarcity of water resources. Therefore, optimal irrigation management by monitoring the soil is essential, and assessing soil hydraulic properties and water flow dynamics is presented as a first measure. For this purpose, the control of volumetric water content, θ, and pressure head, h, is required. This study adopted two types of monitoring strategies in the same experimental plot to control θ and h in the vadose zone: i) non-automatic and more time-consuming; ii) automatic connected to a datalogger. Water flux was modelled with Hydrus-1D using the data collected from both acquisition strategies independently (3820 daily values for the automatic; less than 1000 for the non-automatic). Goodness-of-fit results reported a better adjustment in case of automatic sensors. Both model outputs adequately predicted the general trend of θ and h, but with slight differences in computed annual drainage (711 mm and 774 mm). Soil hydraulic properties were inversely estimated from both data acquisition systems. Major differences were obtained in the saturated volumetric water content, θs, and the n and α van Genuchten model shape parameters. Saturated hydraulic conductivity, Ks, shown lower variability with a coefficient of variation range from 0.13 to 0.24 for the soil layers defined. Soil hydraulic properties were better assessed through automatic data acquisition as data variability was lower and accuracy was higher.

(Table 1) Sediment temperature and thermal conductivity at DSDP Leg 86 Holes

Using a new temperature recording instrument recently developed at the Woods Hole Oceanographic Institution, downhole temperature measurements were made at five sites during Deep Sea Drilling Project Leg 86. The instrument, which can be installed in the shoe of the hydraulic piston corer, allows measurements of sediment temperature to be made simultaneously with the collection of sediment cores. A numerical procedure was applied to correct the temperature disturbance caused by the corer's friction with the sediment. Detailed temperature profiles constructed from the data were combined with the measurement of thermal conductivity to calculate heat flow. Heat flow values were generally low at all sites of Leg 86, consistent with the age of the lithosphere (>100 m.y.) in the Northwestern Pacific Basin.

Hydraulic properties of layered soils influence survival of Rhodes grass (Chloris gayana Kunth.) during water stress

Survival of vegetation on soil-capped mining wastes is often impaired during dry seasons due to the limited amount of water stored in the shallow soil capping. Growth and survival of Rhodes grass (Chloris gayana) during soil drying on various layered capping sequences constructed of combinations of topsoil, subsoil, seawater-neutralised residue sand and low grade bauxite was determined in a glasshouse. The aim was to describe the survival of Rhodes grass in terms of plant and soil water relationships. The soil water characteristic curve and soil texture analysis was a good predictor of plant survival. The combination of soil with a high water holding capacity and low soil water diffusivity (e.g. subsoil with high clay contents) with soil having a high water holding capacity and high diffusivity (e.g. residue sand) gave best survival during drying down (up to 88 days without water), whereas topsoil and low grade bauxite were unsuitable (plants died within 18-39 days). Clayey soil improved plant survival by triggering a water stress response during peak evaporative water demand once residue sand dried down and its diffusivity fell below a critical range. Thus, for revegetation in seasonally dry climates, soil capping should combine one soil with low diffusivity and one or more soils with high total water holding capacity and high diffusivity.

The role of recharge and evapotranspiration as hydraulic drivers of ion concentrations in shallow groundwater on Everglades tree islands, Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year-round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration-to-recharge rates were elevated, while low evapotranspiration-to-recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands.

The Role of Recharge and Evapotranspiration as Hydraulic Drivers of Ion Concentrations in Shallow Groundwater on Everglades Tree Islands, Florida (USA)

Recently, evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year-round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007 to 2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration-to-recharge rates were elevated, while low evapotranspiration-to-recharge rates resulted in a substantial increase of the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion, the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands.

Evaluating temporal changes in hydraulic conductivities near karst-terrain dams: Dokan Dam (Kurdistan-Iraq)

Dam sites provide an outstanding opportunity to explore dynamic changes in the groundwater flow regime because of the high hydraulic gradient rapidly induced in their surroundings. This paper investigates the temporal changes of the hydraulic conductivities of the rocks and engineered structures via a thorough analysis of hydrological data collected at the Dokam Dam, Iraq, and a numerical model that simulates the Darcian component of the seepage. Analysis of the data indicates increased seepage with time and suggests that the hydraulic conductivity of the rocks increased as the conductivity of the grout curtain decreased. Conductivity changes on the order of 10-8m/s, in a 20-yr period were quantified using the numerical analysis. It is postulated that the changes in hydraulic properties in the vicinity of Dokan Dam are due to suspension of fine materials, interbedded in small fissures in the rocks, and re-settlement of these materials along the curtain. Consequently, the importance of the grout curtain to minimize the downstream seepage, not only as a result of the conductivity contrast with the rocks, but also as a barrier to suspended clay sediments, is demonstrated. The numerical analysis also helped us to estimate the proportion of the disconnected karstic conduit flow to the overall flow.

Potential hydraulic barrier performance of cyclic organic carbonate modified bentonite complexes against hyper salinity

The effect of adding glycerol carbonate (GC) or propylene carbonate (PC) to sodium (Na)-bentonite on the hydraulic performance of geosynthetic clay liners (GCLs) under hypersaline conditions is examined. Fluid loss (FL), swell index (SI) and solution retention capacity (SRC) measurements were carried out to compare the potential hydraulic performance of these two cyclic organic carbonates (COCs) as bentonite modifiers. A modified FL test enabled quantitative measurement of both the water retention characteristics of untreated and COC modified bentonites as well as calculation of hydraulic conductivity values. Tests under aggressively saline conditions (ionic strength, I ≥ 1 M of NaCl and ≥3 M of CaCl2) showed that at a mass ratio of 1:1 (GC to bentonite), the FL of a GC-Na-bentonite was ≈40–104 mL in NaCl and ≈61–91 mL in CaCl2. This was about 10–20 mL and 70–200 mL, respectively, lower than that of a comparable PC-Na-bentonite (1:1 PC to bentonite) and untreated Na-bentonite. Greater swelling (SI) and greater solution retention capacity (SRC) was observed for the GC treated Na-bentonite compared to untreated Na-bentonite in all salt solutions, and for PC-Na-bentonite at high ionic strength of both NaCl and CaCl2 solutions, demonstrating the superior hydraulic barrier performance of COC-bentonites under severely saline conditions. Experiments conducted in flexible-wall permeameters with I = 3 M CaCl2 showed approximately one order of magnitude lower (∼10−11 m/s vs ∼1.9 × 10−10 m/s) hydraulic conductivity of GC treated bentonite cake compared to the k value of the untreated Na-bentonite cake. Calculated hydraulic conductivity from fluid loss tests estimated the measured values in a conservative way (overestimation).

Hydraulic performance of geosynthetic clay liners to sulfuric acid solutions

The ability of geosynthetic clay liners (GCLs) to contain acidic mining leachates is examined. The results of saturated hydraulic conductivity (k) of two GCLs permeated with sulfuric acid solutions (H2SO4) at 0.015M, 0.125M and 0.5M concentrations are reported. Also, the saturated k values of consolidated (35kPa) bentonite cakes made from sodium bentonite extracted from both GCLs were compared to a commonly used magnesium-sodium form bentonite. Chemical compatibility and effects of pre-hydration and effective stress were assessed as part of this study. Results indicated that an increased acid concentration (ionic strength) increased the k of all tested specimens. The ratio of the k0.5 values for non-prehydrated specimens permeated with 0.5M H2SO4 to the kw values for specimens permeated with deionized (DI) water (k0.5/kw) ranged from 10 to 110. Pre-hydration (50-140% water content) and increased effective stress (35-200kPa) improved the performance of GCLs (lower k). Strong correlations were observed between k and liquid limit and swell index parameters independent of pre-hydration and effective stress in this study. However, care should still be taken when using these correlations to evaluate hydraulic performance because the intrinsic micro-structure properties of bentonite, such as porosity, should also be considered. This work showed that, for example, high SI of bentonite does not translate necessarily to a better hydraulic performance of GCLs.