An analytical solution for predicting transient seepage into ditch drains from a ponded field

An analytical solution is provided for predicting time dependent seepage into an array of equally spaced parallel ditch drains in a homogeneous and anisotropic soil medium underlain by an impervious layer and receiving water from a ponded horizontal field of infinite extent. The solution can account for both unequal levels of water in the adjacent drains and variable depths of ponding at the soil surface. The validity of the developed model is tested by first reducing it to a steady state solution and then comparing predictions obtained from it for a few flow situations with corresponding predictions obtained from the analytical works of others. A numerical comparison of the developed model for a flow situation is also carried out using MODFLOW. The surface discharge distribution is found to show relatively greater uniformity at the early stages of simulation but with the progress of time, the extent of uniformity is found to reduce particularly for cases where the soil is subjected to a uniform depth of ponding. However, even when a soil surface is subjected to a constant depth of ponding, a high anisotropy ratio (ratio of horizontal to vertical hydraulic conductivity of soil) of the soil alone may lead to a marked improvement on the uniformity of the surface discharge distribution at all times in comparison to a soil having a low anisotropy ratio. A better uniformity of surface discharge may also be achieved by suitably adjusting the depths of ponding over the surface of the soil – regions close to the ditches be provided with zero or negligible depths of ponding and the ponding depths may be made to progressively increase with the increase in distance from the ditch faces. As the developed analytical model is of a general nature, it is hoped that the solution provided herein will lead to a better and realistic design of ditch drainage networks for controlling waterlogged areas and in reclaiming salt affected soils.     

STUDY OF INFILTRATION INTO A HETEROGENEOUS SOIL USING MAGNETIC RESONANCE IMAGING

In this paper, the feasibility of using magnetic resonance imaging (MRI) to study water infiltration into a heterogeneous soil is examined, together with its difficulties and limitations. MRI studies of ponded water infiltration into an undisturbed soil core show that the combination of one- and two-dimensional imaging techniques provides a visual and non-destructive means of monitoring the temporal changes of soil water content and the moisture profile, and the movement of the wetting front. Two-dimensional images show air entrapment in repetitive ponded infiltration experiments. During the early stages of infiltration, one-dimensional images of soil moisture profiles clearly indicate preferential flow phenomena. The observed advance of wetting fronts can be described by a linear relationship between the square root of infiltration time (√t) and the distance of the wetting front from the soil surface. Similarly, the cumulative infiltration is also directly proportional to √t. Furthermore, from the MRI infiltration moisture profiles, it is possible to estimate the parameters that feature in infiltration equations. © 1997 by John Wiley & Sons, Ltd.     

Occurrence of water ponding on soil surfaces depending on infiltration rates on Mongolian rangeland

The occurrence of water ponding on soil surfaces during and after heavy rainfall produces surface run‐off or surface water accumulation in low‐lying areas, which might reduce the water supply to soils and result in a reduction of the soil water that plants can use, especially in arid climates. On Mongolian rangeland, we observed ponded water on the surface of a specific soil condition subjected to a heavy rainfall of 30 mm/hr. By contrast, ponded water was not observed for the same type of soil where livestock grazing had been removed for 6–8 years via a fence or for nearby soil containing less clay. We measured the infiltration rate (the saturated hydraulic conductivity of the surface soil, K_s) of the three sites by applying ponded water on the soil surface (an intake rate test). The results showed that K_s in the rangeland was lower than the rainfall intensity in the site where water ponded on the soil surface; however, K_s of the soil inside of the fence has recovered to 3 times that of the soil outside of the fence to exceed the rainfall intensity. Heavy rainfall that exceeds the infiltration rate occurs several times a year at the livestock grazing site where we observed ponded water. Slight water repellency of the soil reduces rain infiltration to increase the possibility of surface ponding for the soil.     

DOC budgets of drained peat catchments: implications for DOC production in peat soils

The blocking of drainage ditches in peat has been proposed as a possible mitigation strategy for the widely observed increases in dissolved organic carbon (DOC) concentrations from northern peatlands. This study tested the hypothesis that drain‐blocking could lead to lower DOC concentrations by measuring the DOC export from a series of small peat‐covered catchments over a period of 2 years. Six catchments were chosen: two were pristine that had never been drained; three where drains had been blocked (one in 1995, and two in 2003); and a control peat drain catchment where the drain was left unblocked throughout the study. In the case where drains were blocked as part of thus study, the drains were observed for 2 months before blocking and 2 years after blocking. The results show that: (i) high concentrations of DOC can come from water ponded in the drain; (ii) the DOC export (flux of DOC per area of catchment) from the six study catchments shows a high degree of positive correlation with both catchment size and water yield; (iii) distinctly lower DOC export with water yield was observed for the catchments containing higher‐order channels (>27 500 m²) as opposed to single drain catchments (>7500 m²); (iv) drain‐blocking resulted in a statistically significant decrease in DOC export (average was 39%) but the effect upon DOC concentration explained only 1% of the variance in the data. The results suggest that drain blocking works by decreasing the flow from the drain, not by changing the production of DOC in the peat. The change in export with catchment size implies a considerable removal of DOC from large catchments. Copyright © 2009 John Wiley & Sons, Ltd.     

LE TRITIUM DANS LES PRÉCIPITATIONS DE L'ASIE DU SUD-EST ET LE PHÉNOMÉNE DE LA MOUSSON EN 1966 ET EN 1967

Abstract

The problem of non-steady flow of water in a soil-plant system can be described by adding a sink term to the continuity equation for soil water flow. In this paper the sink term is defined in two different ways. Firstly it is considered to be dependent on the hydraulic conductivity of the soil, on the difference in pressure head between the soil and the root-soil interface and some root effectiveness function. Secondly the sink is taken to be a prescribed function of the soil water content. The partial differential equation applying to the first problem is solved by both a finite difference (FD 1) and a finite element (FE 1) technique, that applying to the second problem by a finite difference approach (FD 2). The purpose of this paper is to verify the numerical models against field measurements, to compare the results obtained by the three numerical methods and to show how the finite element method can be applied to complex but realistic two-dimensional flow situations. Two examples are given. The first concerns one-dimensional flow and it compares numerical results with those obtained experimentally in the field from water balance studies on red cabbage (Brassica oleracea L. ‘Rode Herfst’) grown on a clay soil in the presence of a water table. The second example describes two-dimensional flow in a complex field situation in the Netherlands where flow takes place under cropped field conditions through five anisotropic layers. Water is supplied to the system by infiltration from two unlined ditches and is withdrawn from the system by evapotranspiration and by leakage to an underlying pumped aquifer.     

打桩荷载作用下自由场土体振动衰减规律研究

为研究打桩荷载作用下自由场土体振动衰减规律，建立了考虑桩-土相互作用的二维有限元数值模型，并通过Lamb问题解析解验证了数值模型的有效性。通过分析打桩深度、土体阻尼比、打桩荷载等级和土质条件等因素的影响，研究了土体表面振动特性及振动衰减规律。参数分析表明，打桩深度对微振动的影响较小，在距振源一定距离处的土体表面振动响应基本保持一致；土体阻尼比对土体表面振动的影响显著，阻尼比越小，土体表面振动响应越剧烈；不同场地软硬条件影响微振动的限制距离，在一定距离范围内，土质越软，土体表面振动响应越显著，防振距离越长。基于参数分析结果，对峰值速度衰减曲线进行拟合，拟合公式计算结果与模拟结果较吻合，可为振动敏感建筑场地的选择提供参考。     

Study on wall-type gravel drains as liquefaction countermeasure for underground structures

One of the most dramatic causes of damage to engineering structures during earthquakes has been the development of soil liquefaction beneath and around the structures. In order to dissipate the excess pore water pressures near structures, gravel drains are usually employed. In this study, the use of recycled concrete crushed stones as gravel drain materials is addressed. In order to investigate the performance of wall-type gravel drains, two series of shaking table tests were performed. The test results showed that gravel drains, when appropriate grain size distribution is considered, effectively dissipate the excess pore water pressure underneath the structure, and consequently reduce the magnitude of uplift. To supplement the laboratory tests, finite element analyses were also performed. For specified structure, ground and earthquake conditions, there is a critical width of gravel drain at which no uplift of structure will occur. The results of the model tests and the finite element analyses were then employed in developing design charts for determining the critical width of gravel drain to prevent buoyant rise of structure when the surrounding soil mass liquefies.     

城镇化下流域不透水面扩张对洪峰的影响——以南京秦淮河为例

长江下游秦淮河流域近年来由于城市化崛起导致不透水面迅速扩张,改变了流域水文过程,导致暴雨洪水灾害风险增大.本文以南京秦淮河流域为例,基于1988—2015年间下垫面和水文气象资料建立了流域水文模型,通过不透水面扩张情景分析,探讨了 1988—2015年间不透水面空间扩张及对流域洪水过程的影响.研究结果表明:(1)秦淮河全流域1988—2015年不透水率从3.92%增长到19.11%,且不同区域扩张速度有所差异;(2)2006—2015年不透水面情景下的洪峰流量平均涨幅大于城市化初期;受流域上下游位置和下垫面地形条件的影响,流域溧水河和句容河两河源处的不透水面变化对洪峰的影响较流域下游出口处更显著;(3)秦淮河流域及不同位置的不透水面扩张情景下,小洪水的洪峰响应均大于大洪水,且不透水面扩张发生在下游主干河流域时的大、小洪水洪峰涨幅差距略大于河源流域.     

The possible hydrologic effects of the proposed lignite open‐cast mining in Drama lignite field,Greece

The present study investigates the possible hydrologic effects of the proposed lignite open‐cast mining in Drama lignite field (north Greece). Recent years have seen a rapid increase in surface mining. This activity has generated a growing concern for the potential environmental impacts associated with large scale surface mining. In order to achieve a safe mine operation and allow extraction of lignite to considerable depths, extensive dewatering by pumping will be necessary, while at the same time it is desirable to avoid presence of overpumping conditions in the broader area. Based on stratigrafic, hydrologic and hydrogeologic data, a three‐dimensional finite difference model was developed in order to simulate the dewatering process of the western part of the lignite open‐cast mine in Drama and to predict both spatially and temporally the decline of ground water level down to the lignite surface. The dewatering of the part of the aquifer which underlies the mine area will influence the hydrological conditions of the broader region. The most important anticipated effects will be the abandonment of shallow wells as well as the decrease of ground water pumping rates of deep wells. Aquifer discharge towards the ditches of the study area will cease and there will be an inversion of ground water flow from the ditches towards the underlying aquifer. Dewatering activities will probably result in minor subsidence of the nearby peat deposits of Drama Philippi marshes. Moreover, sand pumping as well as the presence of gasses is likely to cause local subsidence phenomena, mainly in the pit slopes. Copyright © 2007 John Wiley & Sons, Ltd.     

Mitigation of railway induced ground vibration by heavy masses next to the track

The effectiveness of heavy masses next to the track as a measure for the reduction of railway induced ground vibration is investigated by means of numerical simulations. It is assumed that the heavy masses are placed in a continuous row along the track forming a wall. Such a continuous wall could be built as a gabion wall and also used as a noise barrier. Since the performance of mitigation measures on the transmission path strongly depends on local ground conditions, a parametric study is performed for a range of possible designs in a set of different ground types. A two-and-a-half dimensional coupled finite element–boundary element methodology is used, assuming that the geometry of the problem is uniform in the direction along the track. It is found that the heavy masses start to be effective above the mass–spring resonance frequency which is determined by the dynamic stiffness of the soil and the mass of the wall. At frequencies above this resonance frequency, masses at the soil׳s surface hinder the propagation of surface waves. It is therefore beneficial to make the footprint of the masses as large and stiff as possible. For homogeneous soil conditions, the effectiveness is nearly independent of the distance behind the wall. In the case of a layered soil with a soft top layer, the vibration reduction strongly decreases with increasing distance from the wall.     

The impact of ditch blocking on fluvial carbon export from a UK blanket bog

We investigated the effects of ditch blocking on fluvial carbon concentrations and fluxes at a 5‐year, replicated, control‐intervention field experiment on a blanket peatland in North Wales, UK. The site was hydrologically instrumented, and run‐off via open and blocked ditches was analysed for dissolved organic carbon (DOC), particulate organic carbon, dissolved carbon dioxide, and dissolved methane. DOC was also analysed in peat porewater and overland flow. The hillslope experiment was embedded within a paired control‐intervention catchment study, with 3 years of preblocking and 6 years of postblocking data. Results from the hillslope showed large reductions in discharge via blocked ditches, with water partly redirected into hillslope surface and subsurface flows, and partly into remaining open ditches. We observed no impacts of ditch blocking on DOC, particulate organic carbon, dissolved carbon dioxide or methane in ditch waters, DOC in porewaters or overland flow, or stream water DOC at the paired catchment scale. Similar DOC concentrations in ditch water, overland flow, and porewater suggest that diverting flow from the ditch network to surface or subsurface flow had a limited impact on concentrations or fluxes of DOC entering the stream network. The subdued response of fluvial carbon to ditch blocking in our study may be attributable to the relatively low susceptibility of blanket peatlands to drainage, or to physical alterations of the peat since drainage. We conclude that ditch blocking cannot be always be expected to deliver reductions in fluvial carbon loss, or improvements in the quality of drinking water supplies.     

Hydrological impact of roadside ditches in an agricultural watershed in Central New York: implications for non‐point source pollutant transport

Nonpoint source pollution and hydromodification are the leading causes of impairment to our nation's rivers and streams. Roadside ditch networks, ubiquitous in both rural and urban landscapes, intercept and shunt substantial quantities of overland runoff and shallow groundwater to stream systems. By altering natural flowpaths, road ditches contribute not only to hydromodification but also potentially to nonpoint‐source (NPS) pollution by acting as hydrological links between agricultural fields and natural streams. Unfortunately, the impacts of these alterations on watershed hydrology and water quality are not well understood. Through a series of field measurements, including field surveys and discharge monitoring, this study examined the effect of road ditch networks on basin morphometry, field‐ and watershed‐scale hydrology, and pollutant transport in a 38 km² agricultural watershed in south‐central NY. Salient findings include the following: (i) 94% of road ditches discharged to natural streams, effectively doubling the drainage density; (ii) on average, road ditches increased peak and total event flows in their receiving streams by 78% and 57%, respectively, but displayed significant variation across ditches; and (iii) ditches intercepted large quantities of surface and subsurface runoff from agricultural fields and therefore represent efficient conduits for the transport of agricultural NPS pollutants to sensitive receiving waterbodies. Our results provide useful information for hydrologists who wish to further understand how artificial drainage may be affecting watershed hydrology and for managers and engineers tasked with designing appropriate flood and NPS pollution control measures. Copyright © 2012 John Wiley & Sons, Ltd.     

Seasonal and event‐based drivers of runoff and phosphorus export through agricultural tile drains under sandy loam soil in a cool temperate region

Frequent algal blooms in surface water bodies caused by nutrient loading from agricultural lands are an ongoing problem in many regions globally. Tile drains beneath poorly and imperfectly drained agricultural soils have been identified as key pathways for phosphorus (P) transport. Two tile drains in an agricultural field with sandy loam soil in southern Ontario, Canada were monitored over a 28‐month period to quantify discharge and the concentrations and loads of dissolved reactive P (DRP) and total P (TP) in their effluent. This paper characterizes seasonal differences in runoff generation and P export in tile drain effluent and relates hydrologic and biogeochemical responses to precipitation inputs and antecedent soil moisture conditions. The generation of runoff in tile drains was only observed above a clear threshold soil moisture content (~0.49 m³·m^?3 in the top 10 cm of the soil; above field capacity and close to saturation), indicating that tile discharge responses to precipitation inputs were governed by the available soil‐water storage capacity of the soil. Soil moisture content approached this threshold throughout the non‐growing season (October – April), leading to runoff responses to most events. Concentrations of P in effluent were variable throughout the study but were not correlated with discharge (p > 0.05). However, there were significant relationships between discharge volume (mm) and DRP and TP loads (kg ha^?1) for events occurring over the study period (R² ≥ 0.49, p ≤ 0.001). This research has shown that the hydrologic and biogeochemical responses of tile drains in a sandy loam soil can be predicted to within an order of magnitude from simple hydrometric data such as precipitation and soil moisture once baseline conditions at a site have been determined. Copyright © 2016 John Wiley & Sons, Ltd.     

Sediment delivery from agricultural land to rivers via subsurface drainage

Diffuse sources of sediment and sediment‐associated nutrients are of increasing environmental concern because of their impacts on receiving water courses. The aim of the research reported here was to monitor the outflow from four field (land) drains at two farms in the English Midlands in order to estimate the quantity of sediment delivered to the local rivers and the most likely sources and processes involved. A multiparameter sediment unmixing model was employed, using environmental magnetic, geochemical and radionuclide tracers in order to determine the most likely origin of sediments transported through the drains. Results demonstrated that there was a generally linear relationship between drainflow sediment loss and drainflow volume and that the majority (>70%) of the sediment exported from the drains was derived from topsoil. Macropore flow through heavily cracked soils is supported by the data to be the most likely means of sediment delivery to the drains. In one catchment, drains contributed over 50% of the annual sediment budget. Spatial and temporal variations in the sources of sediment reaching one drain outlet were investigated in detail. A link between soil moisture deficit (SMD) and the frequency of high‐intensity rainfall events was used to explain the appearance and persistence of a new sediment source in this drain after October 1998. It is concluded that field drains have the potential to be significant conduits of sediment and agrochemicals in a wide variety of environments in the UK. It is also suggested that this potential may increase if projected climate change leads to more intense rainfall events and increases in SMD across a greater area of the UK. Copyright © 2005 John Wiley & Sons, Ltd.     

Peat slides in Teesdale and Weardale,northern pennines,july 1983: Description and failure mechanisms

Peat landslides that occurred in Upper Teesdale and Weardale in 1983 are described. The slides occurred after thunderstorm rainfall of a very rare intensity and return period. The slope geomorphological features and the vegetation associated with the slide areas are detailed. Baseline geotechnical data are used in a simple retrospective stability analysis, treating the peat slides as shallow translational slips. The analysis demonstrates the probable important role played by open moor drains, natural soil pipes, and prior cracking of the peat (owing to dessication or slow mass movement) in the rapid transfer of water to the failure zone within an impervious sensitive clay beneath the peat. The frequency of peat landslides and the visual persistence of sites is considered briefly in relation to substratum and peat stability within a context of climatic change. It is suggested that there may be a gradation of type from the classic Irish bog-burst, through failures of intermediate character, to peat slips involving rigid peat blocks. Antecedent conditions and the mechanism of failure need not be similar for mass movements at either end of this spectrum.     

FREEZE/THAW EFFECTS ON RILL AND GULLY EROSION IN THE NORTHWESTERN WHEAT AND RANGE REGION

Hydrology of the Northwestern Wheat and Range Region （NWRR） of the Pacific Northwest USA is dominated by winter events. Fifty-five to sixty-five percent of the precipitation occurs from November through March, and formation of impermeable frost, intensified by excessive tillage and tillage pans, has long been a major factor in rill and gully formation in the region. Saturated zones can form above tillage pans, freeze solid, and significantly reduce infiltration. Frost heaved surface soils thaw and weaken with warming temperatures or rain. Under these conditions, runoff from rain or snowmelt or a combination of the two is inevitable, concentrating in rills and channels, and carrying with it the loosened soil. Slopes are frequently quite steep and there may be little deposition above the toe slope. Classical over-fall head-cut gullies are uncommon. Concentrated flow channels form because of collection of water from impervious areas such as conventionally tilled fall seeded small grains and bare grass seed fields. In some soils, gullies are created by seepage from saturated layers above permanent restrictive layers in the soil. Gullies can also result from terrace failures due to rodent burrows and low compaction at the time of construction. In naturally unconsolidated soil, rodent activity can lead to piping failures that remove large quantities of subsurface soil and can become gullies. Rill measurements on conventionally tilled fields in Oregon, Washington, and Idaho were used to determine coefficients for relationships between slope length and steepness and size of rills. This paper describes these studies and the results as well as analysis of long-term erosion and weather records from southeastern Washington.     

A mathematical model of soil moisture spatial distribution on the hill slopes of the Loess Plateau

Based on important factors that affect soil moisture spatial distribution, such as the slope gradients, land use, vegetation cover, and surface water diffusion characteristics together with field measurements of soil moisture data obtained from the surface soil under different land use structures, a soil moisture spatial distribution model was established. The diffusion degree coefficient of surface water for different vegetations was estimated from soil moisture values obtained from field measurements. The model can be solved using the finite unit method. The soil moisture spatial distribution on the hill slopes in the Loess Plateau were simulated by the model. A comparison of the simulated values with measurement data shows that the model is a good fit.     

Surface flow boundary conditions in modeling land subsidence due to fluid withdrawal

Land subsidence due to subsurface fluid (water, gas, oil) withdrawal is often predicted by either finite element or finite difference numerical models based on coupled poroelastic theory, where the soil is represented as a semi-infinite medium bounded by the traction-free (ground) surface. One of the variables playing a most important role on the final outcome is the flow condition used on the traction-free boundary, which may be assumed as either permeable or impermeable. Although occasionally justified, the assumption of no-flow surface seems to be in general rather unrealistic. A permeable boundary where the fluid pressure is fixed to the external atmospheric pressure appears to be more appropriate. This paper addresses the response, in terms of land subsidence, obtained with a coupled poroelastic finite element model that simulates a distributed pumping from a horizontal aquifer confined between two relatively impervious layers, and takes either a permeable boundary surface, i.e., constant hydraulic potential, or an impermeable boundary, i.e., a zero Neumann flow condition. The analysis reveals that land subsidence is rather sensitive to the flow condition implemented on the traction-free boundary. In general, the no-flow condition leads to an overestimate of the predicted ground surface settlement, which could even be 1 order of magnitude larger than that obtained with the permeable boundary.     

A 2.5D coupled FE-BE model for the prediction of railway induced vibrations

Ground vibrations induced by railway traffic at grade and in tunnels are often studied by means of two-and-half dimensional (2.5D) models that are based on a Fourier transform of the coordinate in the longitudinal direction of the track. In this paper, the need for 2.5D coupled finite element-boundary element models is demonstrated in two cases where the prediction of railway induced vibrations is considered. A recently proposed novel 2.5D methodology is used where the finite element method is combined with a boundary element method, based on a regularized boundary integral equation. In the formulation of the boundary integral equation, Green's functions of a layered elastic halfspace are used, so that no discretization of the free surface or the layer interfaces is required. In the first case, two alternative models for a ballasted track on an embankment are compared. In the first model, the ballast and the embankment are modelled as a continuum using 2.5D solid elements, whereas a simplified beam representation is adopted in the second model. The free field vibrations predicted by both models are compared to those measured during a passage of the TGVA at a site in Reugny (France). A very large difference is found for the free field response of both models that is due to the fact that the deformation of the cross section of the embankment is disregarded in the simplified representation. In the second case, the track and free field response due to a harmonic load in a tunnel embedded in a layered halfspace are considered. A simplified methodology based on the use of the full space Green's function in the tunnel–soil interaction problem is investigated. It is shown that the rigorous finite element-boundary element method is required when the distance between the tunnel and the free surface and the layer interfaces of the halfspace is small compared to the wavelength in the soil.     

A model to couple overland flow and infiltration into macroporous vadose zone

Most vegetated land surfaces contain macropores that may have a significant effect on the rate of infiltration of water under ponded conditions on the ground surface. Owing to the small-scale variations of the land topography (microtopography), only portions of the land area may get ponded during the process of overland flow. As the macropores transmit water at much higher rates than the primary soil matrix, higher macropore activation in ponded areas produces larger effective infiltration rates into the soil. Therefore, overland flow and infiltration into the macroporous vadose zone are interrelated. Representing the microtopographic variation of the land surface by a simple sine wave function, a method was developed to relate the ponding area to the average ponding depth which was determined by overland flow. A numerical model coupling overland flow and infiltration into the macroporous vadose zone was developed. Overland flow was simulated using the St. Venant equations with the inertia terms neglected. A single macropore model was used to simulate the infiltration into the macroporous vadose zone. The interaction between overland flow and the infiltration into the macroporous vadose zone was analyzed for a hypothetical watershed. The sensitivity analysis revealed that the interaction of macropore flow and overland flow is significant. For the conditions tested, the macropore flow and the overland flow were found to be more sensitive to the macroporosity and less sensitive to the microtopographic surface variation.