Vegetation impact on stream chemical fluxes: Mule Hole watershed (South India)

The proportion of chemical elements passing through vegetation prior to being exported in a stream was quantified for a forested tropical watershed(Mule Hole, South India) using an extensive hydrological and geochemical monitoring at several scales. First, a solute annual mass balance was established at the scale of the soil-plant profile for assessing the contribution of canopy interaction and litter decay to the solute fluxes of soil inputs (overland flow) and soil outputs (pore water flow as seepages). Second, based on the respective contributions of overland flow and seepages to the stream flow as estimated by a hydrological lumped model, we assigned the proportion of chemical elements in the stream that transited through the vegetation at both flood event (End Member Mixing Analysis) and seasonal scales. At the scale of the 1D soil-plant profile, leaching from the canopy constituted the main source of K above the ground surface. Litter decay was the main source of Si, whereas alkalinity, Ca and Mg originated in the same proportions from both sources. The contribution of vegetation was negligible for Na. Within the soil, all elements but Na were removed from the pore water in proportions varying from 20% for Cl to 95% for K: The soil output fluxes corresponded to a residual fraction of the infiltration fluxes. The behavior of K, Cl, Ca and Mg in the soil-plant profile can be explained by internal cycling, as their soil output fluxes were similar to the atmospheric inputs. Na was released from soils as a result of Na-plagioclase weathering and accompanied by additional release of Si. Concentration of soil pore water by evapotranspiration might limit the chemical weathering in the soil. Overall, the solute K, Ca, Mg, alkalinity and Si fluxes associated with the vegetation turnover within the small experimental watershed represented 10-15 times the solute fluxes exported by the stream, of which 83-97% transited through the vegetation. One important finding is that alkalinity and Si fluxes at the outlet were not linked to the ``current weathering'' of silicates in this watershed. These results highlight the dual effect of the vegetation cover on the solute fluxes exported from the watershed: On one hand the runoff was limited by evapotranspiration and represented only 10% of the annual rainfall, while on the other hand, 80-90% of the overall solute flux exported by the stream transited through the vegetation. The approach combining geochemical monitoring and accurate knowledge of the watershed hydrological budget provided detailed understanding of several effects of vegetation on stream fluxes: (1) evapotranspiration (limiting), (2) vertical transfer through vegetation from vadose zone to ground surface (enhancing) and (3) redistribution by throughfalls and litter decay. It provides a good basis for calibrating geochemical models and more precisely assessing the role of vegetation on soil processes. (C) 2014 Elsevier Ltd. All rights reserved.

Planar reflection of gaseous detonations

Pipes containing flammable gaseous mixtures may be subjected to internal detonation. When the detonation normally impinges on a closed end, a reflected shock wave is created to bring the flow back to rest. This study built on the work of Karnesky (2010) and examined deformation of thin-walled stainless steel tubes subjected to internal reflected gaseous detonations. A ripple pattern was observed in the tube wall for certain fill pressures, and a criterion was developed that predicted when the ripple pattern would form. A two-dimensional finite element analysis was performed using Johnson-Cook material properties; the pressure loading created by reflected gaseous detonations was accounted for with a previously developed pressure model. The residual plastic strain between experiments and computations was in good agreement.

During the examination of detonation-driven deformation, discrepancies were discovered in our understanding of reflected gaseous detonation behavior. Previous models did not accurately describe the nature of the reflected shock wave, which motivated further experiments in a detonation tube with optical access. Pressure sensors and schlieren images were used to examine reflected shock behavior, and it was determined that the discrepancies were related to the reaction zone thickness extant behind the detonation front. During these experiments reflected shock bifurcation did not appear to occur, but the unfocused visualization system made certainty impossible. This prompted construction of a focused schlieren system that investigated possible shock wave-boundary layer interaction, and heat-flux gauges analyzed the boundary layer behind the detonation front. Using these data with an analytical boundary layer solution, it was determined that the strong thermal boundary layer present behind the detonation front inhibits the development of reflected shock wave bifurcation.

Implications of using On-Farm Flood Flow Capture to recharge groundwater and mitigate flood risks along the Kings River, CA

Two large hydrologic issues face the Kings Basin, severe and chronic overdraft of about 0.16M ac-ft annually, and flood risks along the Kings River and the downstream San Joaquin River. Since 1983, these floods have caused over $1B in damage in today’s dollars. Capturing flood flows of sufficient volume could help address these two pressing issues which are relevant to many regions of the Central Valley and will only be exacerbated with climate change. However, the Kings River has high variability associated with flow magnitudes which suggests that standard engineering approaches and acquisition of sufficient acreage through purchase and easements to capture and recharge flood waters would not be cost effective. An alternative approach investigated in this study, termed On-Farm Flood Flow Capture, involved leveraging large areas of private farmland to capture flood flows for both direct and in lieu recharge. This study investigated the technical and logistical feasibility of best management practices (BMPs) associated with On-Farm Flood Flow Capture. The investigation was conducted near Helm, CA, about 20 miles west of Fresno, CA. The experimental design identified a coordinated plan to determine infiltration rates for different soil series and different crops; develop a water budget for water applied throughout the program and estimate direct and in lieu recharge; provide a preliminary assessment of potential water quality impacts; assess logistical issues associated with implementation; and provide an economic summary of the program. At check locations, we measured average infiltration rates of 4.2 in/d for all fields and noted that infiltration rates decreased asymptotically over time to about 2 – 2.5 in/d. Rates did not differ significantly between the different crops and soils tested, but were found to be about an order of magnitude higher in one field. At a 2.5 in/d infiltration rate, 100 acres are required to infiltrate 10 CFS of captured flood flows. Water quality of applied flood flows from the Kings River had concentrations of COC (constituents of concern; i.e. nitrate, electrical conductivity or EC, phosphate, ammonium, total dissolved solids or TDS) one order of magnitude or more lower than for pumped groundwater at Terranova Ranch and similarly for a broader survey of regional groundwater. Applied flood flows flushed the root zone and upper vadose zone of nitrate and salts, leading to much lower EC and nitrate concentrations to a depth of 8 feet when compared to fields in which more limited flood flows were applied or for which drip irrigation with groundwater was the sole water source. In demonstrating this technology on the farm, approximately 3,100 ac-ft was diverted, primarily from April through mid-July, with about 70% towards in lieu and 30% towards direct recharge. Substantial flood flow volumes were applied to alfalfa, wine grapes and pistachio fields. A subset of those fields, primarily wine grapes and pistachios, were used primarily to demonstrate direct recharge. For those fields about 50 – 75% of water applied was calculated going to direct recharge. Data from the check studies suggests more flood flows could have been applied and infiltrated, effectively driving up the amount of water towards direct recharge. Costs to capture flood flows for in lieu and direct recharge for this project were low compared to recharge costs for other nearby systems and in comparison to irrigating with groundwater. Moreover, the potentially high flood capture capacity of this project suggests significant flood avoidance costs savings to downstream communities along the Kings and San Joaquin Rivers. Our analyses for Terranova Ranch suggest that allocating 25% or more flood flow water towards in lieu recharge and the rest toward direct recharge will result in an economically sustainable recharge approach paid through savings from reduced groundwater pumping. Two important issues need further consideration. First, these practices are likely to leach legacy salts and nitrates from the unsaturated zone into groundwater. We develop a conceptual model of EC movement through the unsaturated zone and estimated through mass balance calculations that approximately 10 kilograms per square meter of salts will be flushed into the groundwater through displacing 12 cubic meters per square meter of unsaturated zone pore water. This flux would increase groundwater salinity but an equivalent amount of water added subsequently is predicted as needed to return to current groundwater salinity levels. All subsequent flood flow capture and recharge is expected to further decrease groundwater salinity levels. Second, the project identified important farm-scale logistical issues including irrigator training; developing cropping plans to integrate farming and recharge activities; upgrading conveyance; and quantifying results. Regional logistical issues also exist related to conveyance, integration with agricultural management, economics, required acreage and Operation and Maintenance (O&M).

Characterisation of virus transport and attenuation in epikarst using short pulse and prolonged injection multi-tracer testing

Attenuation processes controlling virus fate and transport in the vadose zone of karstified systems can strongly influence groundwater quality. This research compares the breakthrough of two bacteriophage tracers (H40/1 and T7), with contrasting properties, at subsurface monitoring points following application onto an overlying composite sequence of thin organic soil and weathered limestone (epikarst). Short pulse multi-tracer test results revealed that T7 (Source concentration, Co=1.8x106pfu/mL) and H40/1(Co=5.9x106pfu/mL) could reach sampling points 10m below ground less than 30 minutes after tracer application. Contrasting deposition rates, determined from simulated tracer responses, reflected the potential of the ground to differentially attenuate viruses. Prolonged application of both T7 (Co=2.3x104pfu/mL) and H40/1 (Co=1.3x105pfu/mL) over a five hour period during a subsequent test, in which ionic strength levels observed at monitoring points rose consistently, corresponded to a rapid rise in T7 levels, followed by a gradual decline before the end of tracer injection; this reflected reaction-limited deposition in the system. T7’s response contrasted with that of H40/1, whose concentration remained constant over a three hour period before declining dramatically prior to the end of tracer injection. Subsequent application of lower ionic strength tracer-free flush water generated a rapid rise in H40/1 levels and a more gradual release of T7. Results highlight the benefits of employing prolonged injection multi-tracer tests for identifying processes not apparent from conventional short pulse tests. Study findings demonstrate that despite rapid transport rates, the epikarst is capable of physicochemical filtration of viruses and their remobilization, depending on virus type and hydrochemical conditions.

Comparative tracing experiments to investigate epikarst structural and compositional heterogeneity

Comparative tracer testing may be used to evaluate the vulnerability of groundwater to specific contaminants by comparing reactive tracer response to that of a simultaneously injected non-reactive “conservative” substance. Conversely, knowledge of tracer reaction with specific materials permits information about subsurface heterogeneity to be inferred. A series of tests completed in the vadose zone overlying a limestone aquifer employed a cocktail of particles along with reactive and non-reactive solute tracers to investigate transport rates between the ground surface and monitoring points approximately 10 m below ground. Short pulse tests revealed both solutes and particulate contaminants could travel at rates of over 10 m/h. Comparison of particle (microorganisms) and non-reactive solute tracer breakthrough revealed that particle tracers experience pore exclusion resulting in higher peak relative concentrations which arrive earlier than those of the solute. Prolonged tracer injection during subsequent experiments confirmed the response observed and illustrated that over 40 % of flow paths between injection and monitoring points were inaccessible to particles, but could allow solutes to pass through them. Similarly, the difference in response between various reactive tracers demonstrated tracers reached monitoring points via multiple flow paths and suggests geochemical heterogeneity plays an important role in influencing tracer behaviour. The results of this investigation highlight the complexity of water flow through the epikarst and the vulnerability of groundwater in karst aquifers to contamination when soil cover is thin to absent.

Geophysical and hydrogeological characterisation of the impacts of on-site wastewater treatment discharge to groundwater in a poorly productive bedrock acquifer

Contaminants discharging from on-site wastewater treatment systems (OSWTSs) can impact groundwater quality, threatening human health and surface water ecosystems. Risk of negative impacts becomes elevated in areas of extreme vulnerability with high water tables, where thin unsaturated intervals limit vadose zone attenuation. A combined geophysical/hydrogeological investigation into the effects of an OSWTS, located over a poorly productive aquifer (PPA) with thin subsoil cover, aimed to characterise effluent impacts on groundwater. Groundwater, sampled from piezometers down-gradient of the OSWTS percolation area displayed spatially erratic, yet temporally consistent, contaminant distributions. Electrical resistivity tomography identified an area of gross groundwater contamination close to the percolation area and, when combined with seismic refraction and water quality data, indicated that infiltrating effluent reaching the water table discharged to a deeper more permeable zone of weathered shale resting on more competent bedrock. Subsurface structure, defined by geophysics, indicated that elevated chemical and microbiological contaminant levels encountered in groundwater samples collected from piezometers, down-gradient of sampling points with lower contaminant levels, corresponded to those locations where piezometers were screened close to the weathered shale/competent rock interface; those immediately up-gradient were too shallow to intercept this interval, and thus the more impacted zone of the contaminant plume. Intermittent occurrence of faecal indicator bacteria more than 100 m down gradient of the percolation area suggested relatively short travel times. Study findings highlight the utility of geophysics as part of multidisciplinary investigations for OSWTS contaminant plume characterisation, while also demonstrating the capacity of effluent discharging to PPAs to impact groundwater quality at distance. Comparable geophysical responses observed in similar settings across Ireland suggest the phenomena observed in this study are more widespread than previously suspected.

Using geotechnical and LIDAR spatial monitoring to determine key environmental slope instability thresholds of a Jurassic coastal landslide: Straidkilly Point, Northern Ireland, UK.

The Antrim Coast Road stretching from the seaport of Larne in the East of Northern Ireland to the famous Giant’s Causeway in the North has a well-deserved reputation for being one of the most spectacular roads in Europe (Day, 2006). At various locations along the route, fluid interactions between the problematic geology, Jurassic Lias Clay and Triassic Mudstone overlain by Cretaceous Limestone and Tertiary Basalt, and environmental variables result in frequent instances of slope instability within the vadose zone. During such instances of instability, debris flows and composite mudflows encroach on the carriageway posing a hazard to road users. This paper examines the site investigative, geotechnical and spatial analysis techniques currently being implemented to monitor slope stability for one site at Straidkilly Point, Glenarm, Northern Ireland. An in-depth understanding of the geology was obtained via boreholes, resistivity surveys and laboratory testing. Environmental variables recorded by an on-site weather station were correlated with measured pore water pressure and soil moisture infiltration dynamic data.
Terrestrial LiDAR (TLS) was applied to the slope for the monitoring of failures, with surveys carried out on a bi-monthly basis. TLS monitoring allowed for the generation of Digital Elevation Models (DEMs) of difference, highlighting areas of recent movement, erosion and deposition. Morphology parameters were generated from the DEMs and include slope, curvature and multiple measures of roughness. Changes in the structure of the slope coupled with morphological parameters are characterised and linked to progressive failures from the temporal monitoring. In addition to TLS monitoring, Aerial LiDARi datasets were used for the spatio-morphological characterisation of the slope on a macro scale. Results from the geotechnical and environmental monitoring were compared with spatial data obtained through Terrestrial and Airborne LiDAR, providing a multi-faceted approach to slope stability characterization, which facilitates more informed management of geotechnical risk by the Northern Ireland Roads Service.

Stable hydrogen and oxygen isotope composition of waters from mine tailings in different climatic environments

The stable isotope composition of waters (delta H-2, delta O-18) can be used as a natural tracer of hydrologic processes in systems affected by acid mine drainage. We investigated the delta H-2 and delta O-18 values of pore waters from four oxidizing sulfidic mine tailings impoundments in different climatic regions of Chile (Piuquenes at La Andina with Alpine climate, Cauquenes and Caren at El Teniente with Mediterranean climate, and Talabre at the Chuquicamata deposit with hyperarid climate). No clear relationship was found between altitude and isotopic composition. The observed displacement of the tailings pore waters from the local meteoric water line toward higher delta O-18 values (by similar to +2% delta O-18 relative to delta H-2) is partly due to water-rock interaction processes, including hydration and O-isotope exchange with sulfates and Fe(III) oxyhydroxides produced by pyrite oxidation. In most tailings, from the saturated zone toward the surface, isotopically different zones can be distinguished. Zone I is characterized by an upward depletion of H-2 and O-18 in the pore waters from the saturated zone and the lowermost vadose zone, due to ascending diffused isotopically light water triggered by the constant loss of water vapor by evaporation at the surface. In zone II, the capillary flow of a mix of vapor and liquid water causes an evaporative isotopic enrichment in H-2 and O-18. At the top of the tailings in dry climate a zone III between the capillary zone and the surface contains isotopically light diffused and atmospheric water vapor. In temperate climates, the upper part of the profile is affected by recent rainfall and zone III may not differ isotopically from zone II.

Petrography, Geochemistry and Diagenesis of Coral Deposits of Kavaratti and Minicoy Islands, Lakshadweep, India

In the present thesis the petrographic, geochemical and digenetic variability of the sediments in the islands of Kavaratti and Minicoy has been investigated .The beach profile studies show that in the Kavaratti lagoon beach, the slope is steeper in the southern and south central part than in the northern end’s marginal deposition is taking place in the northern end of the Kavaratti island, whereas a marginal erosion is observed at the southern end. In Minicoy the slope of the lagoon beach is gentle in the south and is slightly steeper at the northern part of the beach. The southern and northern beach sections show a marginal deposition. Based on the mineralogical and geochemical studies it is concluded that the main digenetic changes observed is the transformation of aragonite to LMC.The transformation takes place mainly in the vadose zone and is caused by the abundance of fresh water infiltration.

Analysis of water-level response to rainfall and implications for recharge pathways in the Chalk aquifer, SE England

Water table response to rainfall was investigated at six sites in the Upper, Middle and Lower Chalk of southern England. Daily time series of rainfall and borehole water level were cross-corretated to investigate seasonal variations in groundwater-level response times, based on periods of 3-month duration. The time tags (in days) yielding significant correlations were compared with the average unsaturated zone thickness during each 3-month period. In general, for cases when the unsaturated zone was greater than 18 m thick, the time tag for a significant water-level response increased rapidly once the depth to the water table exceeded a critical value, which varied from site to site. For shallower water tables, a linear relationship between the depth to the water table and the water-level response time was evident. The observed variations in response time can only be partially accounted for using a diffusive model for propagation through the unsaturated matrix, suggesting that some fissure flow was occurring. The majority of rapid responses were observed during the winter/spring recharge period, when the unsaturated zone is thinnest and the unsaturated zone moisture content is highest, and were more likely to occur when the rainfall intensity exceeded 5 mm/day. At some sites, a very rapid response within 24 h of rainfall was observed in addition to the longer term responses even when the unsaturated zone was up to 64 m thick. This response was generally associated with the autumn period. The results of the cross-correlation analysis provide statistical support for the presence of fissure flow and for the contribution of multiple pathways through the unsaturated zone to groundwater recharge. (c) 2006 Elsevier B.V. All rights reserved.

Modelling solute transport in structured soils: performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO−3, Cl−, PO3−4) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed (one-way ANOVA, p < 0.001 indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale.

Quantifying effects of soil heterogeneity on groundwater pollution at four sites in the USA

Four sites located in the north-eastern region of the United States of America have been chosen to investigate the impacts of soil heterogeneity in the transport of solutes (bromide and chloride) through the vadose zone (the zone in the soil that lies below the root zone and above the permanent saturated groundwater). A recently proposed mathematical model based on the cumulative beta distribution has been deployed to compare and contrast the regions' heterogeneity from multiple sample percolation experiments. Significant differences in patterns of solute leaching were observed even over a small spatial scale, indicating that traditional sampling methods for solute transport, for example the gravity pan or suction Iysimeters, or more recent inventions such as the multiple sample percolation systems may not be effective in estimating solute fluxes in soils when a significant degree of soil heterogeneity is present. Consequently, ignoring soil heterogeneity in solute transport studies will likely result in under- or overprediction of leached fluxes and potentially lead to serious pollution of soils and/or groundwater. The cumulative beta distribution technique is found to be a versatile and simple technique of gaining valuable information regarding soil heterogeneity effects on solute transport. It is also an excellent tool for guiding future decisions of experimental designs particularly in regard to the number of samples within one site and the number of sampling locations between sites required to obtain a representative estimate of field solute or drainage flux.

Concepts and dimensionality in modeling unsaturated water flow and solute transport

Many environmental studies require accurate simulation of water and solute fluxes in the unsaturated zone. This paper evaluates one- and multi-dimensional approaches for soil water flow as well as different spreading mechanisms to model solute behavior at different scales. For quantification of soil water fluxes,Richards equation has become the standard. Although current numerical codes show perfect water balances, the calculated soil water fluxes in case of head boundary conditions may depend largely on the method used for spatial averaging of the hydraulic conductivity. Atmospheric boundary conditions, especially in the case of phreatic groundwater levels fluctuating above and below a soil surface, require sophisticated solutions to ensure convergence. Concepts for flow in soils with macro pores and unstable wetting fronts are still in development. One-dimensional flow models are formulated to work with lumped parameters in order to account for the soil heterogeneity and preferential flow. They can be used at temporal and spatial scales that are of interest to water managers and policymakers. Multi-dimensional flow models are hampered by data and computation requirements.Their main strength is detailed analysis of typical multi-dimensional flow problems, including soil heterogeneity and preferential flow. Three physically based solute-transport concepts have been proposed to describe solute spreading during unsaturated flow: The stochastic-convective model (SCM), the convection-dispersion equation (CDE), and the fraction aladvection-dispersion equation (FADE). A less physical concept is the continuous-time random-walk process (CTRW). Of these, the SCM and the CDE are well established, and their strengths and weaknesses are identified. The FADE and the CTRW are more recent,and only a tentative strength weakness opportunity threat (SWOT)analysis can be presented at this time. We discuss the effect of the number of dimensions in a numerical model and the spacing between model nodes on solute spreading and the values of the solute-spreading parameters. In order to meet the increasing complexity of environmental problems, two approaches of model combination are used: Model integration and model coupling. Amain drawback of model integration is the complexity of there sulting code. Model coupling requires a systematic physical domain and model communication analysis. The setup and maintenance of a hydrologic framework for model coupling requires substantial resources, but on the other hand, contributions can be made by many research groups.