Simulating Nitrate Leaching Profiles in a Highly Permeable Vadose Zone

An approach is developed to simulate leaching of a dissolved chemical constituent in the vadose zone of an aquifer. Specifically, nitrate loading at the water table for different water table depths, for a range of aquifer permeability values, and for different cases of heterogeneity of the aquifer, are considered. Models from the literature are first used to derive soil–water characteristic curves (water retention and hydraulic conductivity) from a grain size distribution curve for unsaturated conditions. Given infiltration from the surface, the initial conditions for the chemical concentration, and the water content profile, leaching of the chemical in the vadose zone is simulated as a function of both time and depth. The methodology is illustrated for a permeable aquifer. Simulations are undertaken using a finite element code for saturated and unsaturated flow. Different scenarios are simulated depending on the heterogeneity of the aquifer and the depth of the water table. Modeling results show that in the example case studied, nitrate concentration loading at the water table does not depend strongly on the position of the water table, but rather on the material properties of the aquifer. The contribution of this endeavor resides in the methodology which allows a prediction of nitrate leaching using only the grain size property of the aquifer. It allows practitioners to obtain a first assessment of leaching with limited data.     

Numerical Modeling on the Effect of Dissolved Oxygen on Nitrogen Transformation and Transport in Unsaturated Porous System

Nitrogen pollution in groundwater resulting from wastewater application to land is a common problem, and it causes a major threat to groundwater-based drinking water supplies. In this study, a numerical model is developed to study the nitrogen species transport and transformation in unsaturated porous media. Further, a new mass transfer module for dissolved oxygen (DO) is incorporated in the one-dimensional numerical model for nitrogen species transport to describe the fate and transport of nitrogen species, dissolved oxygen, dissolved organic carbon (DOC), and biomass. The spatial and temporal variation of dissolved oxygen is incorporated in the model through the mass transfer from gaseous phase to water phase in an unsaturated porous system. The numerical results of the water flow model and single species and multispecies transport model in an unsaturated zone developed for this purpose have been validated with the available analytical/numerical solution. The developed model is applied in clay loam, silt, and sand soils to analyze the transport behavior of nitrogen species under unsaturated condition. The numerical results suggest that the high rate of oxygen mass transfer from the air phase to the water phase positively increases the dissolved oxygen in the applied wastewater and enhances the nitrification process. Because of this high oxygen mass transfer, the nitrate nitrogen concentration significantly increases in the unsaturated zone and the same is transported to a larger depth at higher simulation period. On the other hand, the low rate of oxygen mass transfer implicitly enhances the denitrification process and finally reduces the nitrate nitrogen concentration in the unsaturated zone. The numerical results also show that the nitrate nitrogen transport is rapid in sandy soil when compared with clay loam and silty soils under high oxygen mass transfer rate. In essence, the high oxygen mass transfer rate significantly increases the nitrate nitrogen in the unsaturated zone, especially at a greater depth at larger time levels and eventually affects the groundwater quality.     

Centrifugal study of zone of influence during air-sparging

Air sparging (AS) is one of the groundwater remediation techniques for remediating volatile organic compounds (VOCs) in saturated soil. However, in spite of the success of air sparging as a remediation technique for the cleanup of contaminated soils, to date, the fundamental mechanisms or the physics of air flow through porous media is not well understood. In this study, centrifugal modeling tests were performed to investigate air flow rates and the evolution of the zone of influence during the air sparging under various g-levels. The test results show that with the increase in sparging pressure the mass flow rate of the air sparging volume increases. The air mass flow rate increases linearly with the effective sparging pressure ratio, which is the difference between sparging pressure and hydrostatic pressure normalized with respect to the effective overburden pressure at the sparging point. Also the slope of mass flow rate with effective sparging pressure ratio increases with higher g-levels. This variation of the slope of mass flow rate of air sparging volume versus effective sparging pressure ratio, M, is linear with g-level confirming that the air flow through soil for a given effective sparging pressure ratio only depends on the g-level. The test results also show that with increasing sparging pressure, the zone of influence (ZOI), which consists of the width at the tip of the cone or lateral intrusion and the cone angle, will lead to an increase in both lateral intrusion and the cone angle. With a further increase in air injection pressure, the cone angle reaches a constant value while the lateral intrusion becomes the main contributor to the enlargement of the ZOI. However, beyond a certain value of effective sparging pressure ratio, there is no further enlargement of the ZOI.     

Two-dimensional finite elements model for selenium transport in saturated and unsaturated zones

A two-dimensional finite element model was developed to simulate species of selenium transport in two dimensions in both saturated and unsaturated soil zones. The model considers water, selenate, selenite, and selenomethionine uptake by plants. It also considers adsorption and desorption, oxidation and reduction, volatilization, and chemical and biological transformations of selenate, selenite, and selenomethionine. In addition to simulating water flow, selenate, selenite, and selenomethionine transport, the model also simulates organic and gaseous selenium transport. The developed model was applied to simulate two different observed field data. The simulation of the observed data was satisfactory, with mean absolute error of 48.5 μg/l and mean relative error of 8.9%.     

庞庄煤矿煤矸石淋溶液对周边水体影响分析

为研究庞庄煤矿张小楼煤矸石堆场淋溶液对其场地和周边地表水、地下水水质影响情况,围绕煤矸石堆场及周边采集13组水样做分析比对。结果表明,煤矸石堆场淋溶液对周边地表水体形成普遍污染,且污染物浓度随着污染距离增加逐步降低,一般化学污染指标主要包括溶解性总固体、硫酸盐、铁、锰、氨氮等,毒理学污染指标包括亚硝酸盐、氟化物、砷,超标倍率约为0.02~5.36。受距离及成井深度等因素影响,地下水体未遭受淋溶液明显的污染和影响。     

Phosphorus losses from arable fields in Sweden—effects of field-specific factors and long-term trends

Accelerated eutrophication of surface water is often caused by high phosphorus (P) losses from agricultural fields. Long-term measurements of P concentrations from arable fields are therefore important for understanding of processes and key factors behind losses. Unfortunately, long-time series are difficult to compare due to high variablity, non-normal distribution, precipitation variability and flow relatedness. The objective was to use locally weighted scatterplot smoothing (LOWESS) and Seasonal Kendall Test to account for discharge effects on concentrations of P forms and suspended soil material, and to identify trends for seven agricultural fields included in a Swedish water quality monitoring program. Soil P content and sorption capacity governed concentrations and losses of dissolved reactive phosphorus (DRP). Other factors, such as preferential flow may influence the importance of these factors. Suspended soil material and unreactive phosphorus (UP) concentrations in water discharge depended on soil texture, slope of the field and total soil P content. Statistically significant trends in concentrations of different P forms were found for two fields. These trends are probably a consequence of changes in P additions and soils P saturation levels.     

海南昌江石碌钴铜矿尾矿库重金属污染环境现状调查

为了解海南省昌江县石碌镇钴铜矿尾矿库环境特点,分别采集了矿砂、背景土壤、废水、纳污水体和植物进行分析。结果表明,钴铜矿尾矿库的土壤受As和Cu污染严重,周边土壤受到不同程度的重金属污染;废水经过尾矿库沉淀处理后,SS、As、Cu、Co、Fe等监测指标均有所降低,尾矿库起到沉淀吸附处理作用;纳污河流石碌河水质良好;尾矿库内自然生长的菖蒲和水竹均受到重金属不同程度污染,尾矿库内的甘蔗中Cr、Pb含量均超过食品中污染物限量标准,表明尾矿库内不宜种植可食用作物。     

Linear alkylbenzene sulfonate (LAS) in water of Lake Dianchi—spatial and seasonal variation, and kinetics of biodegradation

Agroforesty systems, which are recommended as a management option to lower the shallow groundwater level and to reuse saline subsurface drainage waters from the tile-drained croplands in the drainage-impacted areas of Jan Joaquin Valley of California, have resulted in excessive boron buildup in the soil root zone. To assess the efficacy of the long-term impacts of soil boron buildup in agroforesty systems, a mathematical model was developed to simulate non-conservative boron transport. The developed dynamic two-dimensional finite element model simulates water flow and boron transport in saturated–unsaturated soil system, including boron sorption and boron uptake by root-water extraction processes. The simulation of two different observed field data sets by the developed model is satisfactory, with mean absolute error of 1.5 mg/L and relative error of 6.5%. Application of the model to three different soils shows that boron adsorption is higher in silt loam soil than that in sandy loam and clay loam soils. This result agrees with the laboratory experimental observations. The results of the sensitivity analysis indicate that boron uptake by root-water extraction process influences the boron concentration distribution along the root zone. Also, absorption coefficient and maximum adsorptive capacity of a soil for boron are found to be sensitive parameters.     

Agricultural Chemicals in the Alluvial Aquifer of a Typical County of the Arkansas Delta

Statistical methods and a Geographic Information System (GIS) were used to investigate potential indicators of ground water vulnerability to agricultural chemical contamination in a representative area of the Mississippi River alluvial aquifer. A total of 47 wells were sampled for analysis of nitrate, phosphorus, potassium, and 13 pesticides commonly-used in the area. Ten soil and hydrogeologic variables and five ground water vulnerability indices were examined to explain the variations of chemical concentrations. The results showed that no individual soil or hydrogeologic variables or their linear combinations could explain more than 25% of the variation of the chemical concentrations. A quadratic response surface model with the values of confining unit thickness, slope, soil permeability, depth to ground water, and recharge rate accounted for 62% of the variation of nitrate, 43% of P, and 83% of K, suggesting that the interactions among soil and hydrogeologic variables were significant. Observed trends of decreasing nitrate and P concentrations with increasing well depth and/or depth to ground water seemed to correlate with carbonate equilibrium in the aquifer and more reduced environment with depth. In view of uncertainties involved, it was recognized that the limitations associated with input data resolution used in GIS and the formulation of leaching indices limited their use for predicting ground water vulnerability. Misuse of pesticides could be another factor that would complicate the relationships between pesticide concentrations and the vulnerability indices.     

The leachate composition of an old waste dump connected to groundwater: Influence of the reclamation works

This study has been carried out in an uncontrolled dump in the suburbs of Buenos Aires, Argentina. This dump still contains solid wastes of different types in direct contact with the human consumption aquifer. After 18 years of overturning the dump was closed in 1992, and the reclamation works started.Three leachate samplings were carried out in three different zones of the dump, selected according to the topography, the age of residues and the degree of water contact. These samplings were done before the reclamation works, at the end, and four months after they were finished. Several chemical parameters, including heavy metals in solid and dissolved fractions, were analyzed.Multivariate analysis (PCA) was performed to help to understand the effects of the labors on leachate composition. The works could be responsible for a temporal increase of dissolved heavy metals and several chemical variables in the leachate. At the sites of the dump where the modifications were carried out, the levels of some elements analyzed recovered their former values a few months after the reclamations works. At the farthest site, the impact and recovery effects of the works on leachate composition, were delayed.     

Translocation of Sb and Ti in an undisturbed floodplain soil after application of Sb2O3 and TiO2 nanoparticles to the surface

The pore water transport of antimony and titanium, applied as nanoparticles (NPs), was studied by spiking stable suspensions of two different nanomaterials on the surface of an undisturbed floodplain soil. For preparation of stable dispersions, two different strategies were followed. (i) Comparable to those used in industrial applications: titanium dioxide nanoparticles, with an average diameter of 99 nm, were prepared by high-energy ball milling in water, whereas for (ii) antimony trioxide (Sb(2)O(3); average diameter 121 nm) a dispersing agent (sodium salt of poly[(naphthaleneformaldehyde)sulfonate] (pNFS) in water) was used. The upper 17 cm of a floodplain soil (river Rhine, Germany) was sampled using the minimally invasive sediment or fauna incubation experiment (SOFIE? two compartment cell; 3 l volume each), which preserved the pore system of the soil. The cells were equipped with 450 and 100 nm filter probes at different depths providing a non-invasive sampling of the pore water. The pore water was sampled at different times (T = 0, 24, 48, 96 and 196 h) and analysed by inductively coupled plasma quadrupole mass spectrometry (ICP-QMS). Sb and Ti were transported via the pore water of the floodplain soil to a depth of 14 cm, corresponding to the maximum cell depth. The highest Sb concentration in the pore water was detected after 24 h at a depth of 5.5-8 cm. Although the spiked concentration was higher for Ti than for Sb, the total Ti concentration in the pore water of the spiked cell was lower. This indicates a stronger agglomeration of TiO(2) NPs or a more intensive interaction of Ti with the solid matrix and a faster transport of Sb towards deeper soil layers. The results show that metal(loid)s from metal oxide NPs are transported in the soil pore water and, hence, have the potential to act as the source of contamination of deeper soil layers after soil surface contamination.     

Soil Water Content Estimated by Support Vector Machine for the Assessment of Shallow Landslides Triggering: the Role of Antecedent Meteorological Conditions

Soil water content is a key parameter for representing water dynamics in soils. Its prediction is fundamental for different practical applications, such as identifying shallow landslides triggering. Support vector machine (SVM) is a machine learning technique, which can be used to predict the temporal trend of a quantity since training from past data. SVM was applied to a test slope of Oltrepò Pavese (northern Italy), where meteorological parameters coupled with soil water content at different depths (0.2, 0.4, 0.6, 1.0, 1.2, 1.4 m) were measured. Two SVM models were developed for water content assessment: (i) model 1, considering rainfall amount, air temperature, air humidity, net solar radiation, and wind speed; (ii) model 2, considering the same predictors of model 1 together with antecedent condition parameters (cumulated rainfall of 7, 30, and 60 days; mean air temperature of 7, 30, and 60 days). SVM model 2 showed significantly higher satisfactory results than model 1, for both training and test phases and for all the considered soil levels. SVM models trends were implemented in a methodology of slope safety factor assessment. For a real event occurred in the tested slope, the triggering time was correctly predicted using data estimated by SVM model based on antecedent meteorological conditions. This confirms the necessity of including these predictors for building a SVM technique able to estimate correctly soil moisture dynamics in time. The results of this paper show a promising potential application of the SVM methodologies for modeling soil moisture required in slope stability analysis.     

An Alternative Process Model of Preferential Contaminant Travel Times in the Unsaturated Zone: Application to Rainier Mesa and Shoshone Mountain, Nevada

Simulating contaminant transport in unsaturated zones with sparse hydraulic property information is a difficult, yet common, problem. When contaminant transport may occur via preferential flow, simple modeling approaches can provide predictions of interest, such as the first arrival of contaminant, with minimal site characterization. The conceptual model for unsaturated zone flow at Rainier Mesa and Shoshone Mountain, Nevada National Security Site, establishes the possibility of preferential flow through lithologies between potential radionuclide sources and the saturated zone. After identifying preferential flow as a possible contaminant transport process, we apply a simple model to estimate first arrival times for conservatively transported radionuclides to reach the saturated zone. Simulated preferential flow travel times at Rainier Mesa are tens to hundreds of years for non-ponded water sources and 1 to 2 months for continuously ponded water sources; first arrival times are approximately twice as long at Shoshone Mountain. These first arrival time results should then be viewed as a worst-case scenario but not necessarily as a timescale for a groundwater-contamination hazard, because concentrations may be very low. The alternative approach demonstrated here for estimating travel times can be useful in situations where predictions are needed by managers for the fastest arrival of contaminants, yet budgetary or time constraints preclude more rigorous analysis, and when additional model estimates are needed for comparison (i.e., model abstraction).     

Classification of moisture and aeration regimes in sub-boreal forest soils

One hundred and two white spruce (Picea glauca (Moench) Voss) stands were studied in the Sub-boreal Spruce zone of British Columbia and were quantitatively classified into seven soil moisture regimes (moderately dry, slightly dry, fresh, moist, very moist, wet, and very wet) according to actual/potential evapotranspiration ratio, depth to gleyed layer or prominent mottling, and depth to groundwater table. The delineated soil moisture regimes demonstrated strong relationships with the composition of understory vegetation and white spruce foliar nutrients and site index. These relationships implied that the three differentiating characteristics used in the classification provided a good estimation of growing-season soil water supply. In addition to soil moisture regimes, three soil aeration regimes (adequate, restricted, and deficient) were delineated according to presence or absence of gleyed horizons and groundwater table, slope gradient, and soil texture. These soil aeration regimes helped in explaining the variation in white spruce site index, especially on water-surplus sites. Thus, an integrated classification of soil moisture-aeration regimes was proposed to explain the effect of soil moisture and aeration on white spruce productivity.The results of this study gave further evidence that soil moisture and aeration regimes, differentiated on the basis of climatic data and soil morphological properties, are useful measures of soil moisture and aeration conditions in sub-boreal forest soils.     

Optimization of fertirrigation efficiency in strawberry crops by application of fuzzy logic techniques

A high level of price support has favoured intensive agriculture and an increasing use of fertilisers and pesticides. This has resulted in the pollution of water and soils and damage to certain eco-systems. The target relationship that must be established between agriculture and environment can be called "sustainable agriculture". In this work we aim at relating strawberry total yield with nitrate concentration in water at different soil depths. To achieve this objective, we have used the Predictive Fuzzy Rules Generator (PreFuRGe) tool, based on fuzzy logic and data mining, by means of which the dose that allows a balance between yield and environmental damage minimization can be determined. This determination is quite simple and is done directly from the obtained charts. This technique can be used in other types of crops permitting one to determine in a precise way at which depth the appropriate dose of nitrate fertilizer must be correctly applied, on the one hand providing the maximum yield but, on the other hand, with the minimum loss of nitrates that leachate through the saturated zone polluting aquifers.     

Computation of groundwater resources and recharge in Chithar River Basin, South India

Groundwater recharge and available groundwater resources in Chithar River basin, Tamil Nadu, India spread over an area of 1,722 km² have been estimated by considering various hydrological, geological, and hydrogeological parameters, such as rainfall infiltration, drainage, geomorphic units, land use, rock types, depth of weathered and fractured zones, nature of soil, water level fluctuation, saturated thickness of aquifer, and groundwater abstraction. The digital ground elevation models indicate that the regional slope of the basin is towards east. The Proterozoic (Post-Archaean) basement of the study area consists of quartzite, calc-granulite, crystalline limestone, charnockite, and biotite gneiss with or without garnet. Three major soil types were identified namely, black cotton, deep red, and red sandy soils. The rainfall intensity gradually decreases from west to east. Groundwater occurs under water table conditions in the weathered zone and fluctuates between 0 and 25 m. The water table gains maximum during January after northeast monsoon and attains low during October. Groundwater abstraction for domestic/stock and irrigational needs in Chithar River basin has been estimated as 148.84 MCM (million m³). Groundwater recharge due to monsoon rainfall infiltration has been estimated as 170.05 MCM based on the water level rise during monsoon period. It is also estimated as 173.9 MCM using rainfall infiltration factor. An amount of 53.8 MCM of water is contributed to groundwater from surface water bodies. Recharge of groundwater due to return flow from irrigation has been computed as 147.6 MCM. The static groundwater reserve in Chithar River basin is estimated as 466.66 MCM and the dynamic reserve is about 187.7 MCM. In the present scenario, the aquifer is under safe condition for extraction of groundwater for domestic and irrigation purposes. If the existing water bodies are maintained properly, the extraction rate can be increased in future about 10 % to 15 %.     

Soil solution extraction techniques for microbial ecotoxicity testing: a comparative evaluation

The suitability of two different techniques (centrifugation and Rhizon sampler) for obtaining the interstitial pore water of soil (soil solution), integral to the ecotoxicity assessment of metal contaminated soil, were investigated by combining chemical analyses and a luminescence-based microbial biosensor. Two different techniques, centrifugation and Rhizon sampler, were used to extract the soil solution from Insch (a loamy sand) and Boyndie (a sandy loam) soils, which had been amended with different concentrations of Zn and Cd. The concentrations of dissolved organic carbon (DOC), major anions (F- , CI-, NO3, SO4(2-)) and major cations (K+, Mg2+, Ca2+) in the soil solutions varied depending on the extraction technique used. Overall, the concentrations of Zn and Cd were significantly higher in the soil solution extracted using the centrifugation technique compared with that extracted using the Rhizon sampler technique. Furthermore, the differences observed between the two extraction techniques depended on the type of soil from which the solution was being extracted. The luminescence-based biosensor Escherichia coli HB101 pUCD607 was shown to respond to the free metal concentrations in the soil solutions and showed that different toxicities were associated with each soil, depending on the technique used to extract the soil solution. This study highlights the need to characterise the type of extraction technique used to obtain the soil solution for ecotoxicity testing in order that a representative ecotoxicity assessment can be carried out.     

Comparison of soil solution sampling techniques to assess metal fluxes from contaminated soil to groundwater

The unsaturated zone plays a major role in elemental fluxes in terrestrial ecosystems. A representative chemical analysis of soil pore water is required for the interpretation of soil chemical phenomena and particularly to assess Trace Elements (TEs) mobility. This requires an optimal sampling system to avoid modification of the extracted soil water chemistry and allow for an accurate estimation of solute fluxes. In this paper, the chemical composition of soil solutions sampled by Rhizon® samplers connected to a standard syringe was compared to two other types of suction probes (Rhizon®?+?vacuum tube and Rhizon®?+?diverted flow system). We investigated the effects of different vacuum application procedures on concentrations of spiked elements (Cr, As, Zn) mixed as powder into the first 20 cm of 100-cm columns and non-spiked elements (Ca, Na, Mg) concentrations in two types of columns (SiO₂ sand and a mixture of kaolinite?+?SiO₂ sand substrates). Rhizon® was installed at different depths. The metals concentrations showed that (i) in sand, peak concentrations cannot be correctly sampled, thus the flux cannot be estimated, and the errors can easily reach a factor 2; (ii) in sand?+?clay columns, peak concentrations were larger, indicating that they could be sampled but, due to sorption on clay, it was not possible to compare fluxes at different depths. The different samplers tested were not able to reflect the elemental flux to groundwater and, although the Rhizon®?+?syringe device was more accurate, the best solution remains to be the use of a lysimeter, whose bottom is kept continuously at a suction close to the one existing in the soil.     

Saturated hydraulic conductivity of soils in the Horqin Sand Land of Inner Mongolia, northern China

Water is a limiting factor to plant growth in Horqin Sand Land of China. Knowledge of soil saturated hydraulic conductivity (K _sat) is of importance because K _sat influences soil evaporation and water cycling at various scales. In order to analyze the variation of K _sat along with sand types and soil depths, and its relationship with soil physiochemical properties, six typical lands were chosen, including mobile dune, fixed dune, pine woodland, poplar woodland, grassland, and cropland, and K _sat was measured in situ by Guelph Permeameter at each type of land. Soil bulk density, organic matter content, and soil particle size distribution were determined in parallel with K _sat measurement. The results showed that (1) The averaged K _sat was decreased in the order: mobile dune > fixed dune > pine woodland > poplar woodland > grassland > cropland; changes in K _sat varied considerably as soil depth increased, e.g., the changes of K _sat along with soil depth in fixed dune was fitted by exponential model, but it was fitted by parabola model in the pine woodland and grassland. (2) The K _sat values of fixed dune and mobile dune were varied considerably among three slope positions (dune top, windward slope, and leeward slope). (3) The relationships of K _sat and soil physiochemical property revealed that soil bulk density, organic matter content, and coarse sand fraction (2～0.1 mm) were the key factors affecting K _sat in Horqin Sand Land. Compared with clay and silt content proportion, sand fraction in this region showed a more significant positive correlation with K _sat.     

Accumulation of Tobacco mosaic virus (TMV) at different depths clay and loamy sand textural soils due to tobacco waste application

The effects of tobacco waste (TW) application to the soil surface on the accumulation of Tobacco mosaic virus (TMV) in clay and loamy sand textural soils at various depths were investigated in two different fields. The tobacco waste had been found to be infected with TMV. Eighteen months after TW application to the soil surface, soils were sampled at 20 cm intervals through to 80 cm depth. The DAS-ELISA method was performed to determine infection of soil with TMV. The viruses persisted in clay soil for a long period compared with loamy sand soil. There was no accumulation of TMV at any depth of loamy sand soil in Experimental Field 2. TMV adsorption to soil particles in 0-60 cm depth of clay soil was determined in all TW treatments in Experimental Field 1. The highest ELISA Absorbance (A405) values in all treatments were determined in the 20-40 cm soil depth that had the highest clay content. ELISA A405 values of TMV at different depths of clay soil gave significant correlations with clay content (r = 0.793**), EC values (r = 0.421**) and soil pH (r = -0.405**). Adsorption of TMV to net negatively charged clay particle surfaces increased with increasing EC values of soil solution. Decreasing soil pH and infiltration rate increased adsorption of TMV to clay particles. Higher infiltration rate and lower clay content in loamy sand soil caused leaching of TMV from the soil profile.