A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Uniform and lateral preferential flows under flood irrigation at field scale

Flood irrigation is globally one of the most used irrigation methods. Typically, not all water that is applied during flood irrigation is consumed by plants or lost to evaporation. Return flow, the portion of applied water from flood irrigation that returns back to streams either via surface or subsurface flow, can constitute a large part of the water balance. Few studies have addressed the connection between vertical and lateral subsurface flows and its potential role in determining return flow pathways due to the difficulty in observing and quantifying these processes at plot or field scale. We employed a novel approach, combining induced polarization, time‐lapse electrical resistivity tomography, and time‐lapse borehole nuclear magnetic resonance, to identify flow paths and quantify changes in soil hydrological conditions under nonuniform application of flood irrigation water. We developed and tested a new method to track the wetting front in the subsurface using the full range of inverted resistivity values. Antecedent soil moisture conditions did not play an important role in preferential flow path activation. More importantly, boundaries between lithological zones in the soil profile were observed to control preferential flow pathways with subsurface run‐off occurring at these boundaries when saturation occurred. Using the new method to analyse time‐lapse resistivity measurements, we were able to track the wetting front and identify subsurface flow paths. Both uniform infiltration and preferential lateral flows were observed. Combining three geophysical methods, we documented the influence of lithology on subsurface flow processes. This study highlights the importance of characterizing the subsurface when the objective is to identify and quantify subsurface return flow pathways under flood irrigation.     

Soil sealing and soil water content under no‐tillage and conventional tillage in irrigated corn: Effects on grain yield

The main objective of this research was to analyse the effect of soil management on soil sealing and on soil water content under contrasting tillage practices and its influence on corn yield. The experimental research was carried out in a field cultivated with irrigated corn differentiated into three zones representing a gradient of soil texture (Z1, Z2, and Z3, i.e., increasingly coarser). Two plots under different soil management practices (conventional intensive tillage, CT, and no‐tillage, NT) were selected in each zone. The susceptibility to sealing of each soil and the steady infiltration rates were evaluated in the laboratory subjecting the soils to rainfall simulation applied at an intensity of 25 mm h^?1. In addition, soil porosity under each treatment was quantified. Soil water content (0–90 cm depth) was determined gravimetrically at the beginning and the end of the growing cycle and at the surface (0–5 cm) during three growing seasons and continuously at two depths (5–15 and 50–60 cm) during the last growing cycle. Soil water content was simulated using the SIMPEL model, which was calibrated for the experimental conditions. Corn yield and above‐ground biomass were also analysed. Significant differences in soil sealing among zones, with decreasing soil sealing for coarser textures, and treatments were observed with infiltration rates that were near twice in NT than in CT, being the effect of soil cover significant in the reduction of soil detachment and soil losses. NT showed higher soil water content than CT, especially in the surface layers. Above‐ground biomass production was smaller in CT than in NT, and in the areas with higher sealing susceptibility was 30% to 45% smaller than in other zones, reaching the smallest values in Z1. A similar reduction in corn yield was observed between treatments being smaller in CT than in NT. No‐tillage has been confirmed as an effective technique that benefits soil physical properties as well as crop yields in relation to CT, being its impact greater in soils susceptible to sealing.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Analytical solution to 1D coupled water infiltration and deformation in two‐layer unsaturated soils

An analytical solution to 1D coupled water infiltration and deformation in layered soils is derived using a Laplace transformation. Coupling between seepage and deformation, and initial conditions defined by arbitrary continuous pore‐water pressure distributions are considered. The analytical solutions describe the transient pore‐water pressure distributions during 1D, vertical infiltration toward the water table through two‐layer unsaturated soils. The nonlinear coupled formulations are first linearized and transformed into a form that is solvable using a Laplace transformation. The solutions provide a reliable means of comparing the accuracy of various numerical methods. Parameters considered in the coupled analysis include the saturated permeability (k_s), desaturation coefficient (α), and saturated volumetric water content (θ_s) of each soil layer, and antecedent and subsequent rainfall infiltration rates. The analytical solution demonstrates that the coupling of seepage and deformation plays an important role in water infiltration in layered unsaturated soils. A smaller value of α or a smaller absolute value of the elastic modulus of the soil with respect to a change in soil suction (H) for layered unsaturated soils means more marked coupling effect. A smaller absolute value of H of the upper layer soil also tends to cause more marked coupling effect. A large difference between the saturated coefficients of permeability for the top and bottom soil layers leads to reduced rainfall infiltration into the deep soil layer. The initial conditions also play a significant role in the pore‐water pressure redistribution and coupling effect. Copyright © 2011 John Wiley & Sons, Ltd.     

Melt infiltration into quartzite during partial melting in the Little Cottonwood Contact Aureole (UT,USA): implication for xenocryst formation

Melt infiltration into quartzite took place due to generation and migration of partial melts within the high‐grade metamorphic rocks of the Big Cottonwood (BC) formation in the Little Cottonwood contact aureole (UT, USA). Melt was produced by muscovite and biotite dehydration melting reactions in the BC formation, which contains pelite and quartzite interlayered on a centimetre to decimetre scale. In the migmatite zone, melt extraction from the pelites resulted in restitic schollen surrounded by K‐feldspar‐enriched quartzite. Melt accumulation occurred in extensional or transpressional domains such as boudin necks, veins and ductile shear zones, during intrusion‐related deformation in the contact aureole. The transition between the quartzofeldspathic segregations and quartzite shows a gradual change in texture. Here, thin K‐feldspar rims surround single, round quartz grains. The textures are interpreted as melt infiltration texture. Pervasive melt infiltration into the quartzite induced widening of the quartz–quartz grain boundaries, and led to progressive isolation of quartz grains. First as clusters of grains, and with increasing infiltration as single quartz grains in the K‐feldspar‐rich matrix of the melt segregation. A 3D–μCT reconstruction showed that melt formed an interconnected network in the quartzites. Despite abundant macroscopic evidence for deformation in the migmatite zone, individual quartz grains found in quartzofeldspathic segregations have a rounded crystal shape and lack quartz crystallographic orientation, as documented with electron backscatter diffraction (EBSD). Water‐rich melts, similar to pegmatitic melts documented in this field study, were able to infiltrate the quartz network and disaggregate grain coherency of the quartzites. The proposed mechanism can serve as a model to explain abundant xenocrysts found in magmatic systems.     

Infiltration,runoff and soil loss in Andisols affected by forest fire (Canary Islands,Spain)

Depending on the severity of the fire, forest fires may modify infiltration and soil erosion processes. Rainfall simulations were used to determine the hydrological effects of fire on Andisols in a pine forest burned by a wildfire in 2007. Six burned zones with different fire severities were compared with unburned zones. Infiltration, runoff and soil loss were analysed on slopes of 10% and 30%. Forest floor and soil properties were evaluated. Unburned zones exhibited relatively low infiltration (23 and 16 mm h^?1 on 10% and 30% slope angles, respectively) and high average runoff/rainfall ratios (43% and 50% on 10% and 30% slope angles, respectively), which were associated with the extreme water repellency of the forest floor. Nonetheless, this layer seems to provide protection against raindrop impact and soil losses were found to be low (8 and 16 g m^?2 h^?1 for 10% and 30% slope angles, respectively). Soil cover, soil structure and water repellency were the main properties affected by the fire. The fire reduced forest floor and soil repellency, allowing rapid infiltration. Moreover, a significant decrease was noted in soil aggregate stabilities in the burned zones, which limited the infiltration rates. Consequently, no significant differences in infiltration and runoff were found between the burned and the unburned zones. The decrease in post‐fire soil cover and soil stability resulted in order‐of‐magnitude increases in erosion. Sediment rates were 15 and 31 g m^?2 h^?1 on the 10% and 30% slope angles, respectively, in zones affected by light fire severity. In the moderate fire severity zones, these values reached 65 and 260 g m^?2 h^?1 for the 10% and 30% slope angles, respectively. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydro‐mechanical response of collapsible soils under different infiltration events

The behavior of a partially saturated soil during surface‐water infiltration is analyzed by means of an elasto‐plastic constitutive model formulated in terms of effective stress and extended to unsaturated conditions. The model is calibrated considering laboratory‐scale experimental results under suction‐controlled conditions. The wetting process in two collapsing soils, initially loaded at in situ stresses, is simulated by imposing two different boundary conditions: surface ponding and water flow. The stress paths resulting from the imbibition process are analyzed at different points inside the layer. Copyright © 2015 John Wiley & Sons, Ltd.     

河水入渗过程中近岸带地下水中砷的形成与影响因素研究

沈阳黄家水源地是我国北方地区典型的傍河地下水水源地,近岸带地下水中铁(Fe)、锰(Mn)、砷(As)含量严重超标。为查明地下水中As的来源与影响因素,对研究区河水、地下水以及土壤样品进行采集与测试,分析了水样常规指标与碳硫稳定同位素、土样中典型矿物、砷的含量及赋存形态。结果表明,研究区河水中As含量很低,而地下水中As含量普遍超标。河水入渗初期,氧化性河水使部分含As矿物发生氧化而释放As;随着河水入渗,地下水向还原环境转变,含As的Fe/Mn矿物发生还原性溶解,地下水中As含量逐渐升高。研究区典型矿物有黄铁矿、菱铁矿、软锰矿、赤铁矿、针铁矿、菱锰矿等,通过可交换态砷解吸、有机质结合态砷氧化、铁锰氧化物结合态砷还原性溶解等,介质中的As释放至地下水中。地下水中As含量与酸碱度(pH)、氧化还原电位(Eh)呈一定负相关,与溶解有机碳(DOC)、 {\mathrm{HCO}}_3^{-}、Fe、Mn含量呈正相关。     

物质组成特征对碎石土斜坡变形破坏影响研究

为研究物质组成空间分布特征对降雨下碎石土斜坡变形破坏的影响,以不同成因碎石土斜坡的特点为依照,进行了4种物质组成空间分布的碎石土斜坡模型及1种均质斜坡模型试验。结果表明,物质组成空间分布特征对降雨下碎石土斜坡变形破坏所发挥的作用与降雨入渗方向和应力调整方向有重要关系。当土料变化方向与入渗方向大角度相交且与应力调整方向相同时,斜坡渗透性为土料各自的渗透性,斜坡力学性质偏向于组成土料力学性质的最大值。当土料变化方向与入渗方向相同且与应力调整方向大角度相交时,斜坡渗透性接近组成土料渗透性的最小值,斜坡变形破坏受控于其中的软弱土料。降雨下斜坡各处均发生变形,且变形向周围调整并产生影响。破裂面更易发于软弱土料中。该研究有助于斜坡灾害勘察、预测评价及防治向精细化方向发展。