Influence of temporal fluctuations and spatial heterogeneity on pollution transport in porous media

The combined influence of temporal fluctuations and spatial heterogeneity on non-reactive solute transport mechanisms in porous media can be understood by performing simulations of steady and unsteady flow and transport in heterogeneous media. The study focuses on issues such as the degree of heterogeneity, correlation length, separation of the combined effects of temporal and spatial variations, and ergodicity conditions under unsteady flow conditions. It is shown that the effect of temporal variations on solute transport is masked by the strong effect of spatial heterogeneity. There is no obvious difference in plume shape between steady and unsteady flow conditions; the first and the second spatial moments of the plume of the unsteady-state flow condition fluctuate around the steady-state flow condition with the same period of oscillations as the input signal at small storage coefficient (S????0.001). At a relatively high standard deviation in hydraulic conductivity and a small storage coefficient, the unsteady flow condition sharpens the temporal variations in macrodispersion coefficients. The magnitude of the longitudinal macrodispersion coefficient under unsteady flow condition is almost doubled at the maximum values. However, the transverse macrodispersion coefficient fluctuates around zero. The Kubo number and Peclet number ranges are 1.2?C64 and 10?C250, respectively.     

盆地地下水流系统形成与影响因素分析

目前区域（盆地）地下水流系统模拟研究中,常用的定水头与通量两种上边界条件刻画方法与实际条件存在差距。通过对比两种方法的差异和各自适用条件,采用解析法讨论地下水位的形成控制机制,提出了改进后的变通量上边界数值模型,并以鄂尔多斯盆地北部白垩系地下水流系统为例分析了盆地地下水流系统的形成与影响因素。研究表明,鄂尔多斯盆地北部白垩系水流系统地下水位受地形、补给条件和渗透系数三者共同控制,同时特有的气候、地形和岩性组合通过控制地下水位影响地下水流系统的发育演化。采用变通量上边界法探讨上边界条件改变对盆地水流系统的影响,对深刻认识区域地下水流系统形成演化机制,揭示地下水系统与上边界气候变化、植被生态变化之间的相互作用关系具有一定优势。     

氢氧同位素在地下水流系统的重分布：从高程效应到深度效应

大气降水的氢氧同位素含量具有高程效应,降水入渗后参与地下水循环,其高程效应如何受地下水流系统的影响转化为地下水氢氧同位素的深度效应？现有研究对于这个问题缺少定量认识。文章构建单向倾斜盆地和双峰波状盆地的稳态地下水循环理论模型,采用MODFLOW模拟剖面二维地下水流场、采用MT3DMS模拟重同位素分子的对流-弥散过程,得到地下水D和¹⁸O含量的空间分布,探讨了氢氧同位素高程效应在地下水流系统转化为深度效应的机理。结果表明：在单斜盆地,补给区大气降水D和¹⁸O含量的高程效应转化为排泄区地下水δD和δ¹⁸O值随埋深增大而指数型衰减的深度效应;在双峰波状盆地,当含水层渗透性相对入渗强度较大时(K₀/w=1 000),仅发育一个区域地下水流系统,在区域地下水的排泄区δD和δ¹⁸O随埋深增大呈现S形曲线分布;当含水层渗透性相对入渗强度较小时(K0/w=250),双峰波状盆地发育多个局部地下水流系统,区域地下水的排泄区δD和δ¹⁸O随埋深增大呈现S形曲线,而局部地下水排...     

The use of spectral analysis-based exact solutions to characterize topography-controlled groundwater flow

Spectral analysis enhances the ability to analyze groundwater flow at a steady state by separating the top boundary condition into its periodic forms. Specifically, spectral analysis enables comparisons of the impact of individual spatial scales on the total flow field. New exact spectral solutions are presented for analyzing 3D groundwater flow with an arbitrarily shaped top boundary. These solutions account for depth-decaying, anisotropic and layered permeability while utilizing groundwater flux or the phreatic surface as a top boundary condition. Under certain conditions, groundwater flow is controlled by topography. In areas where the groundwater flow is controlled by the topography, the unknown water table is often approximated by the topography. This approximation induces a systematic error. Here, the optimal resolution of digital elevation models (DEMs) is assessed for use as a top boundary in groundwater flow models. According to the analysis, the water-table undulation is smoother than the topography; therefore, there is an upper limit to the resolution of DEMs that should be used to represent the groundwater surface. The ability to represent DEMs of various spectral solutions was compared and the results indicate that the fit is strongly dependent on the number of harmonics in the spectral solution.     

Evidence for deep groundwater flow and convective heat transport in mountainous terrain, Delta County, Colorado, USA

The Tongue Creek watershed lies on the south flank of Grand Mesa in western Colorado, USA and is a site with 1.5 km of topographic relief, heat flow of 100 mW/m², thermal conductivity of 3.3 W m^–1 °C^–1, hydraulic conductivity of 10^-8 m/s, a water table that closely follows surface topography, and groundwater temperatures 3–15°C above mean surface temperatures. These data suggest that convective heat transport by groundwater flow has modified the thermal regime of the site. Steady state three-dimensional numerical simulations of heat flow, groundwater flow, and convective transport were used to model these thermal and hydrological data. The simulations provided estimates for the scale of hydraulic conductivity and bedrock base flow discharge within the watershed. The numerical models show that (1) complex three-dimensional flow systems develop with a range of scales from tens of meters to tens of kilometers; (2) mapped springs are frequently found at locations where contours of hydraulic head indicate strong vertical flow at the water table, and; (3) the distribution of groundwater temperatures in water wells as a function of surface elevation is predicted by the model.     

Natural spatial and temporal variations in groundwater chemistry in fractured,sedimentary rocks: scale and implications for solute transport

Natural tracers (major ions, δ¹⁸O, and O₂) were monitored to evaluate groundwater flow and transport to a depth of 20 m below the surface in fractured sedimentary (primarily shale and limestone) rocks. Large temporal variations in these tracers were noted in the soil zone and the saprolite, and are driven primarily by individual storm events. During nonstorm periods, an upward flow brings water with high TDS, constant δ¹⁸O, and low dissolved O₂ to the water table. During storm events, low TDS, variable δ¹⁸O, and high dissolved O₂ water recharges through the unsaturated zone. These oscillating signals are rapidly transmitted along fracture pathways in the saprolite, with changes occurring on spatial scales of several meters and on a time scale of hours. The variations decreased markedly below the boundary between the saprolite and less weathered bedrock. Variations in the bedrock units occurred on time scales of days and spatial scales of at least 20 m. The oscillations of chemical conditions in the shallow groundwater are hypothesized to have significant implications for solute transport. Solutes and colloids that adsorb onto aquifer solids can be released into solution by decreases in ionic strength and pH. The decreases in ionic strength also cause thermodynamic undersaturation of the groundwater with respect to some mineral species and may result in mineral dissolution. Redox conditions are also changing and may result in mineral dissolution/precipitation. The net result of these chemical variations is episodic transport of a wide range of dissolved solutes or suspended particles, a phenomenon rarely considered in contaminant transport studies.     

指示娘子关泉群水动力环境的水化学—同位素信息分析

本文依据宏量组分、微量元素和氢氧同位氧组成监测资料，分析了娘子关泉群水动力环境及其水文地质演化趋势。除国家和城西泉主要排泄局部流动系统的地下水外，其它各泉均为不同空间尺度流动系统地下水混合、排泄之产物。在人类活动影响下，娘子关泉群水质恶化，流量衰减，水帘洞、程家泉断流，城西、坡底泉也面临断流之危险。目前抽取的水帘洞泉的水资源量组成已与五龙泉和苇泽关泉相似，以区域和中间流动系统地下水为主。     

Concepts for the sustainable management of multi-scale flow systems: the groundwater system within the Laufen Basin, Switzerland

Many groundwater systems consist of multi-scale aquifer units. The exchange processes and rates between these aquifer units are complex. In order to manage such complex systems, a subdivision into different catchments, sub-catchments or groundwater bodies as manageable units is required. The sustainable management of water resources requires a comprehensive view of water-quality and water-quantity aspects not only for water supply issues, but generally also for flood protection and riverine ecosystem functions. Such transformations require an improved understanding of recharge and exchange processes between different aquifer units as well as aquifer-surface water interaction-processes at different spatiotemporal scales. The main objective of this study is to illustrate concepts by defining the geometry and scales of different aquifer units within a sedimentary basin. The Laufen Basin in the Jura Mountains represents a sub-catchment of the River Birs (Switzerland). Its structure is characterized by a pronounced local relief and a series of aquifer units which are typical for many complex groundwater systems in front of mountain chains such as the alpine foreland and the Jura Mountains of Central Europe. A combination of different concepts is required to understand multi-scale flow systems and to describe the various hydrogeological processes. Three concepts are proposed for the Laufen Basin, including: (1) a regional flow-system analysis, based on the concept of hierarchical groundwater flow systems; (2) the river-corridor concept for understanding aquifer-surface water interaction processes; and (3) the calculation of the dynamic vulnerability index and the aquifer base gradient approach for karst flow and fractured flow systems.     

Simulation of spatial and temporal trends in nitrate concentrations at the regional scale in the Upper Dyle basin,Belgium

Models are the only tools capable of predicting the evolution of groundwater systems at a regional scale, by taking into account a large amount of information. This study presents the association of a water balance model (WetSpass) with a groundwater flow and solute transport model (SUFT3D, saturated and unsaturated flow and transport in 3D) in order to simulate the present and future groundwater quality in terms of nitrate in the Upper Dyle basin (439 km²) Belgium. The HFEMC (hybrid finite element mixing cell) method implemented in the SUFT3D code is used to model groundwater flow and nitrate transport. Spatially distributed recharge, modelled with WetSpass, is considered for prescribing the recharge to the groundwater flow model. The feasibility of linking the WetSpass model with the finite-elements SUFT3D code is demonstrated. Time evolution and distribution of nitrate concentration are then simulated using the calibrated model. Nitrate inputs are spatially distributed according to land use. The spatial simulations and temporal trends are compared with previously published data on this aquifer and show good results.     

A new numerical method of considering local longitudinal dispersion in single fractures

The solutions of advection–dispersion equation in single fractures were carefully reviewed, and their relationships were addressed. The classic solution, which represents the resident or flux concentration within the semi‐infinite fractures under constant concentration or flux boundary conditions, respectively, describes the effluent concentration for a finite fracture. In addition, it also predicts the cumulative distribution of solute particle residence time passing through a single fracture under pulse injection condition, based on which a particle tracking approach was developed to simulate the local advection–dispersion in single fractures. We applied the proposed method to investigate the influence of local dispersion in single fractures on the macrodispersion in different fracture systems with relatively high fracture density. The results show that the effects of local dispersion on macrodispersion are dependent on the heterogeneity of fracture system, but generally the local dispersion plays limited roles on marodispersion at least in dense fracture network. This trend was in agreement with the macrodispersion in heterogeneous porous media. Copyright © 2013 John Wiley & Sons, Ltd.     

应用泉钙华环境记录和地下水流动系统探讨娘子关泉群演变历史

从研究娘子关泉群出露区堆积的不同时期泉钙华入手，对泉群自中更新世以来的发育阶段进行了划分，在此基础上，运用地下水流动系统理论分析了泉群的时空演化过程．这一过程包括４个阶段：Ｑ２时期的泉群发育雏形期，Ｑ３时期的泉群发育全盛期，Ｑ４时期的泉群发育相对稳定期和近代泉群流量衰减期．随着时间的推移，泉点出露位置有规律地发生横向迁移和垂向下移．近年来由于人类活动引起区域性地下水位降低，从根本上改变了岩溶水流动系统的动态平衡，导致有的泉已干涸，有的将要干涸，泉群流量总体上呈现持续衰减趋势．     

Groundwater as a geologic agent: An overview of the causes, processes, and manifestations

 The objective of the present paper is to show that groundwater is a general geologic agent. This perception could not, and did not, evolve until the system nature of basinal groundwater flow and its properties, geometries, and controlling factors became recognized and understood through the 1960s and 1970s. The two fundamental causes for groundwater's active role in nature are its ability to interact with the ambient environment and the systematized spatial distribution of its flow. Interaction and flow occur simultaneously at all scales of space and time, although at correspondingly varying rates and intensities. Thus, effects of groundwater flow are created from the land surface to the greatest depths of the porous parts of the Earth's crust, and from a day's length through geologic times. Three main types of interaction between groundwater and environment are identified in this paper, with several special processes for each one, namely: (1) Chemical interaction, with processes of dissolution, hydration, hydrolysis, oxidation-reduction, attack by acids, chemical precipitation, base exchange, sulfate reduction, concentration, and ultrafiltration or osmosis; (2) Physical interaction, with processes of lubrication and pore-pressure modification; and (3) Kinetic interaction, with the transport processes of water, aqueous and nonaqueous matter, and heat. Owing to the transporting ability and spatial patterns of basinal flow, the effects of interaction are cumulative and distributed according to the geometries of the flow systems. The number and diversity of natural phenomena that are generated by groundwater flow are almost unlimited, due to the fact that the relatively few basic types are modified by some or all of the three components of the hydrogeologic environment: topography, geology, and climate. The six basic groups into which manifestations of groundwater flow have been divided are: (1) Hydrology and hydraulics; (2) Chemistry and mineralogy; (3) Vegetation; (4) Soil and rock mechanics; (5) Geomorphology; and (6) Transport and accumulation. Based on such a diversity of effects and manifestations, it is concluded that groundwater is a general geologic agent. Received, December 1998 · Revised, January 1999 · Accepted, January 1999     

Modelling water-table depth in a primary aquifer to identify potential wetland hydrogeomorphic settings on the northern Maputaland Coastal Plain,KwaZulu-Natal,South Africa

The primary aquifer on the Maputaland Coastal Plain in northern KwaZulu-Natal, South Africa, is the principal source of water for rivers, lakes and most of the wetlands in dry periods, and is recharged by these systems in wet periods. Modelling hydrologic conditions that control regional water-table depth can provide insight into the spatial patterns of wetland occurrence and of the persistence of wet conditions that control their character. This project used a groundwater model (MODFLOW) to simulate 10-year water-table fluctuations on the Maputaland Coastal Plain from January 2000 to December 2010, to contrast the conditions between wet and dry years. Remote sensing imagery was used to map “permanent” and “temporary” wetlands in dry and wet years to evaluate the effectiveness of identifying the suitable conditions for their formation using numerical modelling techniques. The results confirm that topography plays an important role on a sub-regional and local level to support wetland formation. The wetlands’ extent and distribution are directly associated with the spatial and temporal variations of the water table in relation to the topographical profile. Groundwater discharge zones in the lowland (1–50 masl) areas support more permanent wetlands with dominantly peat or high organic soil substrates, including swamp forest and most of the permanent open water areas. Most temporary wetlands associated with low-percentage clay occurrence are through-flow low-lying interdune systems characterised by regional fluctuation of the water table, while other temporary wetlands are perched or partially perched. The latter requires a more sophisticated saturated-unsaturated modelling approach.     

Estimation of groundwater recharge using a GIS-based distributed water balance model in Dire Dawa, Ethiopia

Sustainable groundwater management requires knowledge of recharge. Recharge is also an important parameter in groundwater flow and transport models. Spatial variation in recharge due to distributed land-us.e, soil texture, topography, groundwater level, and hydrometeorological conditions should be accounted for in recharge estimation. However, conventional point-estimates of recharge are not easily extrapolated or regionalized. In this study, a spatially distributed water balance model WetSpass was used to simulate long-term average recharge using land-use, soil texture, topography, and hydrometeorological parameters in Dire Dawa, a semiarid region of Ethiopia. WetSpass is a physically based methodology for estimation of the long-term average spatial distribution of surface runoff, actual evapotranspiration, and groundwater recharge. The long-term temporal and spatial average annual rainfall of 626 mm was distributed as: surface runoff of 126 mm (20%), evapotranspiration of 468 mm (75%), and recharge of 28 mm (5%). This recharge corresponds to 817 l/s for the 920.12 km² study area, which is less than the often-assumed 1,000 l/s recharge for the Dire Dawa groundwater catchment.     

GIS-based spatial and temporal prediction system development for regional land subsidence hazard mitigation

An integrated GIS-based approach for establishing a spatial and temporal prediction system for groundwater flow and land subsidence is proposed and applied to a subsidence-progressed Japanese coastal plain. Various kinds of fundamental data relating to groundwater flow and land subsidence are digitized and entered into a GIS database. A surface water hydrological cycle simulation is performed using a GIS spatial data operation for the entire plain, and the spatial and temporal groundwater infiltration quantity is hereby obtained. Through the data transformation from the GIS database to a groundwater flow code (MODFLOW), a 3D groundwater flow model is established and unsteady groundwater flow simulation for the past 21 years is conducted with results which compare satisfactorily with observed results. Finally, a Visual Basic code is developed for land subsidence calculations considering aquifer and aquitard deformation. Future land subsidence in the plain is predicted assuming different water pumping scenarios, and the results provide important information for land subsidence mitigation decision-making.     

An analytical study on groundwater flow in drainage basins with horizontal wells

Analytical studies on release/capture zones are often limited to a uniform background groundwater flow. In fact, for basin-scale problems, the undulating water table would lead to the development of hierarchically nested flow systems, which are more complex than a uniform flow. Under the premise that the water table is a replica of undulating topography and hardly influenced by wells, an analytical solution of hydraulic head is derived for a two-dimensional cross section of a drainage basin with horizontal injection/pumping wells. Based on the analytical solution, distributions of hydraulic head, stagnation points and flow systems (including release/capture zones) are explored. The superposition of injection/pumping wells onto the background flow field leads to the development of new internal stagnation points and new flow systems (including release/capture zones). Generally speaking, the existence of n injection/pumping wells would result in up to n new internal stagnation points and up to 2n new flow systems (including release/capture zones). The analytical study presented, which integrates traditional well hydraulics with the theory of regional groundwater flow, is useful in understanding basin-scale groundwater flow influenced by human activities.     

Spatial and temporal constraints on regional-scale groundwater flow in the Pampa del Tamarugal Basin,Atacama Desert,Chile

Aquifers within the Pampa del Tamarugal Basin (Atacama Desert, northern Chile) are the sole source of water for the coastal city of Iquique and the economically important mining industry. Despite this, the regional groundwater system remains poorly understood. Although it is widely accepted that aquifer recharge originates as precipitation in the Altiplano and Andean Cordillera to the east, there remains debate on whether recharge is driven primarily by near-surface groundwater flow in response to periodic flood events or by basal groundwater flux through deep-seated basin fractures. In addressing this debate, the present study quantifies spatial and temporal variability in regional-scale groundwater flow paths at 20.5°S latitude by combining a two-dimensional model of groundwater and heat flow with field observations and δ¹⁸O isotope values in surface water and groundwater. Results suggest that both previously proposed aquifer recharge mechanisms are likely influencing aquifers within the Pampa del Tamarugal Basin; however, each mechanism is operating on different spatial and temporal scales. Storm-driven flood events in the Altiplano readily transmit groundwater to the eastern Pampa del Tamarugal Basin through near-surface groundwater flow on short time scales, e.g., 10⁰–10¹ years, but these effects are likely isolated to aquifers in the eastern third of the basin. In addition, this study illustrates a physical mechanism for groundwater originating in the eastern highlands to recharge aquifers and salars in the western Pampa del Tamarugal Basin over timescales of 10⁴–10⁵ years.     

Prediction in hydrogeology: two case histories

An inherent aspect of hydrogeology is the dynamic character of groundwater flow within the more passive lithospheric medium. Prediction therefore requires insight into the spatial pattern of aquifer properties as well as the dynamic character of groundwater flow systems. Both aspects together determine hydrogeological system behaviour in space and time. The spatial characteristics of an aquifer are mainly based on geological and geomorphological structure, whereas the dynamics of the system result predominantly from climate and topography; knowledge of paleoconditions are a prerequisite for sound understanding of present systems. The integration of exploration and forecast is illustrated by two case studies. The first study focuses on an assessment of groundwater resource and their sustainable development in semi-arid Botswana. The second study concerns the impact of land and water management in the Netherlands coastal area on the groundwater regime and its future development, including the redistribution of salt and fresh groundwater on a long time scale.     

A groundwater-basin approach to conceptualize and simulate post-Pleistocene subsurface flow in a semi-arid region, southeastern New Mexico and western Texas, USA

Numerical simulation was used to enhance conceptual understanding of the post-Pleistocene hydrogeology of a layered sequence of clastic and evaporite sediments. This work is part of an effort to evaluate the suitability of the Waste Isolation Pilot Plant (WIPP), New Mexico, USA, as a repository for transuranic waste. The numerical model is three-dimensional, extends laterally to topographic features that form the actual boundaries of a regional groundwater system, and uses a free surface with seepage face as an upper boundary condition to simulate the effect of change in recharge rate on the position of the water table. Simulation results suggest that the modern-day flow field is still adjusting to the drying of the climate that has occurred since the end of the Pleistocene Epoch. A wetter climate at the end of the Pleistocene resulted in a shallow water table, and patterns of groundwater flow were controlled by the intermediate features of the land-surface topography. As the climate became drier and the water table declined, groundwater flow began to increasingly reflect the land-surface topography at the scale of the entire groundwater basin. The modern-day flow pattern has not equilibrated with either the present recharge rate or the position of the water table. Received, November 1998/Revised, December 1999/Accepted, January 2000     

Groundwater pollution due to sugar-mill effluent, at Sonai, Maharashtra, India

 Analyses of 126 samples collected from 18 dug wells in the shallow basaltic aquifer over a period of 7 months have revealed spatial as well as temporal changes in the chemical properties of groundwater. While the temporal changes have been attributed to dilution and concentration phenomena governed by climatic factors, the spatial variations in the geochemical characteristics of groundwater appeared to be related to pollution due to effluents from the Mula Sugar Factory. The cause of groundwater pollution is the effluent carried by a stream flowing through the area. Fluctuations in the groundwater table, influent water quality character of the stream, less capacity to accommodate large volume of effluent and occurrence of zero base flow (under natural conditions) in the stream are the factors favoring infiltration of constituents of waste water into the underlying weathered basaltic aquifer. Pollutants have entered into the shallow aquifer by downward percolation through the zone of aeration to form a recharge mound at the water table and, further, lateral movement below the water table. The plume of polluted groundwater has a lateral extent of a few meters in the upstream area and more than 400 m on either side of the stream in the downstream part. The zone of polluted groundwater has an areal extent of more than 3.5 km². Groundwater is the only source available for drinking and agricultural purposes. It is recommended that the base of the lagoons and the stream used for release of plant effluent should be waterproofed for the protection of groundwater in the Sonai area. Received: 30 April 1997 · Accepted: 23 September 1997