Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. Copyright © 2011 John Wiley & Sons, Ltd.     

Using a H∞ filter assimilation procedure to estimate root zone soil water content

Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.     

The response of CO2 flux to rain pulses at a saline desert

As the substantial component of the ecosystem respiration, soil CO₂ flux is strongly influenced by infrequent and unpredictable precipitation in arid region. In the current study, we investigated the response of soil CO₂ flux to rain pulses at a saline desert in western China. Soil CO₂ flux was measured continuously during the whole growing season of 2009 at six sites. We found that there were remarkable changes in amplitude or diurnal patterns of soil CO₂ flux induced by rainfall events: from bimodal before rain to a single peak after that. Further analysis indicated that there is a significant linear relationship (P < 0.001) between soil CO₂ flux and soil temperature (T_soil). However, a hysteresis between the waveform of diurnal course of CO₂ flux and T_soil was observed: with soil CO₂ flux always peaked earlier than T_soil. Furthermore, a double exponential decay function was fitted to the soil CO₂ flux after rainfall, and total carbon (C) releases were estimated by numerical integration for rainfall events. The relative enhancement and total C release, in association with the rain pulses, was linearly related to the amount of precipitation. According to the size and frequency of rainfall events, the total amount of C release induced by rain pulses was computed as much as 7.88 g C·m^–2 in 2009, equivalent to 10.25% of gross primary production. These results indicated that rain pulses played a significant role in the carbon budget of this saline desert ecosystem, and the size of them was a good indicator of rain‐induced flux enhancement. Copyright © 2012 John Wiley & Sons, Ltd.     

Subsurface lateral flow generation in aspen and conifer‐dominated hillslopes of a first order catchment in northern Utah

Mountain headwater catchments in the semi‐arid Intermountain West are important sources of surface water because these high elevations receive more precipitation than neighboring lowlands. This study examined subsurface runoff in two hillslopes, one aspen dominated, the other conifer dominated, adjacent to a first order stream in snow‐driven northern Utah. Snow accumulation, soil moisture, trenchflow and streamflow were examined in hillslopes and their adjacent stream. Snow water equivalents (SWEs) were greater under aspen stands compared to conifer, the difference increasing with higher annual precipitation. Semi‐variograms of shallow spatial soil moisture patterns and transects of continuous soil moisture showed no increase in soil moisture downslope, suggesting the absence of subsurface flow in shallow (～12 cm) soil layers of either vegetation type. However, a clear threshold relationship between soil moisture and streamflow indicated hillslope–stream connectivity, deeper within the soil profile. Subsurface flow was detected at ～50 cm depth, which was sustained for longer in the conifer hillslope. Soil profiles under the two vegetation types varied, with deep aspen soils having greater water storage capacity than shallow rocky conifer soils. Though SWEs were less under the conifers, the soil profile had less water storage capacity and produced more subsurface lateral flow during the spring snowmelt. Copyright © 2010 John Wiley & Sons, Ltd.     

Impact of grid resolution on the integrated and distributed response of a coupled surface–subsurface hydrological model for the des Anglais catchment,Quebec

Digital elevation models (DEMs) at different resolutions (180, 360, and 720 m) are used to examine the impact of different levels of landscape representation on the hydrological response of a 690‐km² catchment in southern Quebec. Frequency distributions of local slope, plan curvature, and drainage area are calculated for each grid size resolution. This landscape analysis reveals that DEM grid size significantly affects computed topographic attributes, which in turn explains some of the differences in the hydrological simulations. The simulations that are then carried out, using a coupled, process‐based model of surface and subsurface flow, examine the effects of grid size on both the integrated response of the catchment (discharge at the main outlet and at two internal points) and the distributed response (water table depth, surface saturation, and soil water storage). The results indicate that discharge volumes increase as the DEM is coarsened, and that coarser DEMs are also wetter overall in terms of water table depth and soil water storage. The reasons for these trends include an increase in the total drainage area of the catchment for larger DEM cell sizes, due to aggregation effects at the boundary cells of the catchment, and to a decrease in local slope and plan curvature variations, which in turn limits the capacity of the watershed to transmit water downslope and laterally. The results obtained also show that grid resolution effects are less pronounced during dry periods when soil moisture dynamics are mostly controlled by vertical fluxes of evaporation and percolation. Copyright © 2010 John Wiley & Sons, Ltd.     

Seasonal rainfall–runoff relationships in a lowland forested watershed in the southeastern USA

Hydrological processes of lowland watersheds of the southern USA are not well understood compared to a hilly landscape due to their unique topography, soil compositions, and climate. This study describes the seasonal relationships between rainfall patterns and runoff (sum of storm flow and base flow) using 13 years (1964–1976) of rainfall and stream flow data for a low‐gradient, third‐order forested watershed. It was hypothesized that runoff–rainfall ratios (R/P) are smaller during the dry periods (summer and fall) and greater during the wet periods (winter and spring). We found a large seasonal variability in event R/P potentially due to differences in forest evapotranspiration that affected seasonal soil moisture conditions. Linear regression analysis results revealed a significant relationship between rainfall and runoff for wet (r² = 0·68; p < 0·01) and dry (r² = 0·19; p = 0·02) periods. Rainfall‐runoff relationships based on a 5‐day antecedent precipitation index (API) showed significant (r² = 0·39; p < 0·01) correspondence for wet but not (r² = 0·02; p = 0·56) for dry conditions. The same was true for rainfall‐runoff relationships based on 30‐day API (r² = 0·39; p < 0·01 for wet and r² = 0·00; p = 0·79 for dry). Stepwise regression analyses suggested that runoff was controlled mainly by rainfall amount and initial soil moisture conditions as represented by the initial flow rate of a storm event. Mean event R/P were higher for the wet period (R/P = 0·33), and the wet antecedent soil moisture condition based on 5‐day (R/P = 0·25) and 30‐day (R/P = 0·26) prior API than those for the dry period conditions. This study suggests that soil water status, i.e. antecedent soil moisture and groundwater table level, is important besides the rainfall to seasonal runoff generation in the coastal plain region with shallow soil argillic horizons. Copyright © 2011 John Wiley & Sons, Ltd.     

Soil saturation index of salt marshes subjected to spring‐neap tides: a new variable for describing marsh soil aeration condition

Organized spatial distribution of plants (plant zonation) in salt marshes has been linked to the soil aeration condition in the rhizosphere through simplistic tidal inundation parameters. Here, a soil saturation index (ratio of saturation period to tidal period at a soil depth) is introduced to describe the soil aeration condition. This new index captures the effects of not only the tidal inundation period and frequency but also the flow dynamics of groundwater in the marsh soil. One‐dimensional numerical models based on saturated flow with the Boussinesq approximations and a two‐dimensional variably saturated flow model were developed to explore the behaviour of this new soil aeration variable under the influence of spring‐neap tides. Simulations revealed two characteristic zones of soil aeration across the salt marsh: a relatively well aerated near‐creek zone and a poorly aerated interior zone. In the near‐creek zone, soils undergo periodic wetting and drying as the groundwater table fluctuates throughout the spring‐neap cycle. In the interior, the soil remains largely water saturated except for neap tide periods when limited drainage occurs. Although such a change of soil aeration condition has been observed in previous numerical simulations, the soil saturation index provides a clear delineation of the zones that are separated by an ‘inflexion point’ on the averaged index curve. The results show how the saturation index represents the effects of soil properties, tidal parameters and marsh platform elevation on marsh soil aeration. Simulations of these combined effects have not been possible with traditional tidal inundation parameters. The saturation index can be easily derived using relatively simple models based on five non‐dimensional variables. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling the impacts of subsurface drainage on surface runoff and sediment yield in the Le Sueur Watershed,Minnesota, USA

Abstract

Tile drainage influences infiltration and surface runoff and is thus an important factor in the erosion process. Tile drainage reduces surface runoff, but questions abound on its influence on sediment transport through its dense network and into the stream network. The impact of subsurface tiling on upland erosion rates in the Le Sueur River watershed, USA, was assessed using the Water Erosion Prediction Project (WEPP) model. Six different scenarios of tile drainage with varying drainage coefficient and management type (no till and autumn mulch-till) were evaluated. The mean annual surface runoff depth, soil loss rate and sediment delivery ratio (SDR) for croplands, based on a 30-year simulation for the watershed with untiled autumn mulch-till (Scenario 1), were estimated to be 83.5 mm, 0.27 kg/m² and 86.7%, respectively; on no-till management systems (Scenario 4), the respective results were 72.3 mm, 0.06 kg/m² and 88.2%. Tile drains reduced surface runoff, soil loss and SDR estimates for Scenario 1 by, on average, 14.5, 8.1 and 7.9%, respectively; and for Scenario 4 by an estimated 31.5, 22.1 and 20.2%, respectively. The impact of tile drains on surface runoff, soil loss and SDR was greater under the no-till management system than under the autumn mulch-till management system. Comparison of WEPP outputs with those of the Soil Water Assessment Tool (SWAT) showed differences between the two methods.

Editor Z.W. Kundzewicz

Citation Maalim, F.K. and Melesse. A.M., 2013. Modelling impacts of subsurface drainage on surface runoff and sediment yield in the Le Sueur Watershed in Minnesota, USA. Hydrological Sciences Journal, 58 (3), 570–586.     

Microtopography,water storage and flow patterns in a fine‐textured soil under agricultural use

Agricultural use of soils implies tillage and often compaction and therefore influences processes on soil surface and affects infiltration of water into the subsoil. Although many studies on soil surface processes or flow patterns in soils exist, works relating both are rare in literature. We did two tracer experiments with Brilliant Blue FCF on a tilled and compacted plot and a non‐tilled one to investigate water storage on the soil surface during simulated rainfall and changes of soil microtopography, to analyse the associated flow patterns in the soil and to relate both to tillage and compaction. Our results show that storage was larger on the tilled and compacted plot than on the non‐tilled one. After tillage, transport processes above the plough pan were partly disconnected from those underneath because macropores were disrupted and buried by the tillage operation. However, preferential flow along cracks occurred on both plots and the macropores buried below the tillage pan still functioned as preferential flow paths. Therefore, we conclude that the studied soil is susceptible to deep vertical solute propagation at dry conditions when cracks are open, irrespective of tillage and compaction. Copyright © 2012 John Wiley & Sons, Ltd.     

STUDIES ON ENVIRONMENTAL BACKGROUND VALUES IN CHINA

An investigation on environmental background values was made in an area of about 1,140,000 km2, which included temperate and subtropical zones of China. The environmental background values of 142 soil environment units, 18 main soil types, 87 surface water environment units, 8 aquatic organism environment units and 20 underground water environment units were obtained. The rules, causes and effecting factors of regional differentiation of the environmental background values were deduced from over 200,000 various data.