Unsaturated Zone Transport Research Articles

Identifying the nitrate transport and transformations under apple orchards in the loess depositusing stable isotopes of water and nitrate

Understanding the behavior of nitrate in the unsaturated zone is critical for predicting groundwater nitrate pollution. Although this issue has been intensively investigated for shallow soil profiles, it remains poorly understood in deep soils. We collected soil samples up to 20 m depth under farmland and apple orchards with trees that were 16 and 22 years old (A16 and A22). We then determined soil-water content, stable isotopes (δ2H and δ18O) and nitrate content under two land use types. The stable isotopes of nitrate under two apple orchards were measured to a depth of 16 m. The mean nitrate contents in the shallow soils (0–6 m) under apple orchards (78.8 ± 87.8 mg kg−1 for A16, 180 ± 155 mg kg−1 for A22) were significantly higher than that of farmland (2.40 ± 1.56 mg kg−1) (p < 0.01), whereas the nitrate contents in the deep soils (6–20 m) were almost similar between apple orchards and farmland. Ammonium or urea-based nitrogen (N) fertilizers (30–61%) were the predominant contributors to nitrate in the shallow soil layers, followed by soil organic N and manure & sewage, whereas the soil organic N (mean 34%) was the primary source of nitrate in deep soils. Isotopic data revealed that the N transformation was controlled by nitrification, and denitrification was not apparent throughout the soil profiles. The downward transport of nitrate lagged behind soil-water. Nitrate was leached by piston flow in the vadose zone, and the nitrate from N application since the 1980 s has only infiltrated to 6 m depth. The difficulty in recharging deep soils mitigates the transport of nitrate to lower depths. The combination of water-nitrate stable isotopes provides new insights for identifying nitrate sources, transformation, and transport in the thick unsaturated zone.

Transport of perfluorooctanoic acid in unsaturated porous media mediated by SDBS

Per– and polyfluoroalkyl substances (PFAS) contamination is ubiquitous and likely co–exist with hydrocarbon surfactants in the subsurface environment. The transport of PFAS in unsaturated porous media (quartz sand–packed columns) with an anionic hydrocarbon surfactant was investigated using perfluorooctanoic acid (PFOA) and sodium dodecyl benzene sulfonate (SDBS) as representative PFAS and hydrocarbon surfactant at 1.5 and 30 mM NaCl, respectively. Our results show that SDBS inhibited the retardation of PFOA during transport. The degree of the inhibited retardation was more significant at higher concentration of SDBS (CSDBS) for a given concentration of PFOA (CPFOA) and higher ionic strength. In contrast, PFOA only marginally affected the transport of SDBS under all tested conditions. This is mainly due to the greater surface activity of SDBS with respect to PFOA, favoring the competitive preferential adsorption of SDBS at the air–water interface and then the transport of PFOA in the unsaturated sand. This assumption is further supported by the greater air–water adsorption coefficient (Kai) and surface excess (Γ) for SDBS than those for PFOA, as well as the surface tension isotherms that fairly weighted toward to SDBS in the binary PFOA–SDBS system. Our findings indicate the critical role of anionic hydrocarbon surfactant in retarding PFOA transport in unsaturated zone. This should be considered when assessing the transport potential and environmental risk, and developing remedial strategies of PFAS in the subsurface environment.

Simulating contaminant transport in unsaturated and saturated groundwater zones

The demand for clean and adequate water is rising rapidly with increasing population. This growing demand for water necessitates the measurement of the quantity and quality of water. Simulation modeling has become increasingly popular in the last two decades largely because of their predictive ability. This paper reviews the approaches for simulation modeling in groundwater resources management, focusing on models that have been used to simulate contaminant transport through the aquifer system. Recent research papers that have integrated the models for unsaturated and saturated zones have also been studied and described. Integrated models require assessment of the complex interactions between the groundwater zones and the movement of water and solute through them. Due to this, integrated models provide a more accurate modeling approach than models that have been independently developed for saturated and unsaturated zones. Application of such models is encouraged at the regional level to arrive at the best groundwater management decisions. PRACTITIONER POINTS: In the past few decades, modeling of contaminant transport in groundwater systems has seen tremendous applications. A number of models exist that independently simulate flow and solute transport in unsaturated and saturated zones. Recently, focus has been given on developing advanced coupled modeling approaches that require less inputs and run times.

An analytical model of vapor intrusion with fluctuated water table

Intrusion of subsurface contaminants in vapor phase from groundwater into overlying buildings through cracks or other pathways poses potential risk to human health. Previous studies indicated that water table fluctuations may have significant impacts on vapor intrusion, resulting in fluctuated concentrations of indoor volatile organic compounds (VOCs). Analytical solutions for vapor intrusion processes are flexible tools in predicting VOC concentrations and assessing health risks of VOC-contaminated sites. However, an analytical solution considering the impacts of fluctuated water table is not available. In this study, an analytical model of the vapor intrusion with fluctuated water table was developed. The VOC transport in unsaturated zones was described by diffusion equation and the varying water table was treated as a moving boundary. The residual VOCs in the unsaturated zone induced by water table fluctuations were also considered and were described by a source term of the vapor diffusion. The analytical solution was tested with its numerical solution. The results indicate that the fluctuated water table affects the dynamics of VOCs mainly in two ways, changing the length of vapor migration path and delivering the residual VOCs into the unsaturated zone. The VOC concentration at the outlet varies inversely with the water table fluctuation owing to the release processes of the residual VOCs dominate. The analytical solution matches well to the laboratory experiment data in an existing study and could be a simple tool for the risk assessment of vapor intrusion when considering the impacts of water table fluctuations.

Suitability of common models to estimate hydrology and diffuse water pollution in North-eastern German lowland catchments with intensive agricultural land use

Various process-based models are extensively being used to analyze and forecast catchment hydrology and water quality. However, it is always important to select the appropriate hydrological and water quality modeling tools to predict and analyze the watershed and also consider their strengths and weaknesses. Different factors such as data availability, hydrological, hydraulic, and water quality processes and their desired level of complexity are crucial for selecting a plausible modeling tool. This review is focused on suitable model selection with a focus on desired hydrological, hydraulic and water quality processes (nitrogen fate and transport in surface, subsurface and groundwater bodies) by keeping in view the typical lowland catchments with intensive agricultural land use, higher groundwater tables, and decreased retention times due to the provision of artificial drainage. In this study, four different physically based, partially and fully distributed integrated water modeling tools, SWAT (soil and water assessment tool), SWIM (soil and water integrated model), HSPF (hydrological simulation program– FORTRAN) and a combination of tools from DHI (MIKE SHE coupled with MIKE 11 and ECO Lab), have been reviewed particularly for the Tollense River catchment located in North-eastern Germany. DHI combined tools and SWAT were more suitable for simulating the desired hydrological processes, but in the case of river hydraulics and water quality, the DHI family of tools has an edge due to their integrated coupling between MIKE SHE, MIKE 11 and ECO Lab. In case of SWAT, it needs to be coupled with another tool to model the hydraulics in the Tollense River as SWAT does not include backwater effects and provision of control structures. However, both SWAT and DHI tools are more data demanding in comparison to SWIM and HSPF. For studying nitrogen fate and transport in unsaturated, saturated, and river zone, HSPF was a better model to simulate the desired nitrogen transformation and transport processes. However, for nitrogen dynamics and transformations in shallow streams, ECO Lab had an edge due its flexibility for inclusion of user-desired water quality parameters and processes. In the case of SWIM, most of the input data and governing equations are similar to SWAT but it does not include water bodies (ponds and lakes), wetlands and drainage systems. In this review, only the processes that were needed to simulate the Tollense River catchment were considered, however the resulted model selection criteria can be generalized to other lowland catchments in Australia, North-western Europe and North America with similar complexity.

Numerical study on fingered flow and solute transport in unsaturated zone

Numerical study on fingered flow and solute transport in unsaturated zone

Chromium(VI) transport and fate in unsaturated zone and aquifer: 3D Sandbox results

The simulation of Cr(VI) behavior in an unsaturated zone and aquifer, using a 3D experimental set-up were performed to illustrate the distribution, transport and transformation of Cr(VI), and further to reveal the potential harm of Cr(VI) after entering the groundwater. The result indicated that chromium(VI) was transported in the vertical direction, meanwhile, was transported in the horizontal direction under the influence of groundwater flow. The direction and distance away from the pollution source zone had great effect on the chromium(VI) concentration. At the sampling sites near the pollution source zone, there was a sudden increase of chromium(VI) concentration. The concentration of chromium(III) concentration in some random effluent samples was not detected. Chromium had not only transported but also had fraction and specie transformation in the unsaturated zone and aquifer. The relative concentration of residue fraction chromium was decreased with time. The content of Fe–Mn oxide fraction chromium was increased with time. The relative content of exchangeable and carbonate-bound fraction chromium was lower and the content variations were not obvious. Chromium(VI) (91–98%) was first reduced to chromium(III) rapidly. The oxidation reaction occurred later and the relative content of chromium(VI) was increased again. The presence of manganese oxides under favorable soil conditions can promote the reoxidation of Cr(III) to Cr(VI).

Multimodel analysis of anisotropic diffusive tracer‐gas transport in a deep arid unsaturated zone

AbstractGas transport in the unsaturated zone affects contaminant flux and remediation, interpretation of groundwater travel times from atmospheric tracers, and mass budgets of environmentally important gases. Although unsaturated zone transport of gases is commonly treated as dominated by diffusion, the characteristics of transport in deep layered sediments remain uncertain. In this study, we use a multimodel approach to analyze results of a gas‐tracer (SF6) test to clarify characteristics of gas transport in deep unsaturated alluvium. Thirty‐five separate models with distinct diffusivity structures were calibrated to the tracer‐test data and were compared on the basis of Akaike Information Criteria estimates of posterior model probability. Models included analytical and numerical solutions. Analytical models provided estimates of bulk‐scale apparent diffusivities at the scale of tens of meters. Numerical models provided information on local‐scale diffusivities and feasible lithological features producing the observed tracer breakthrough curves. The combined approaches indicate significant anisotropy of bulk‐scale diffusivity, likely associated with high‐diffusivity layers. Both approaches indicated that diffusivities in some intervals were greater than expected from standard models relating porosity to diffusivity. High apparent diffusivities and anisotropic diffusivity structures were consistent with previous observations at the study site of rapid lateral transport and limited vertical spreading of gas‐phase contaminants. Additional processes such as advective oscillations may be involved. These results indicate that gases in deep, layered unsaturated zone sediments can spread laterally more quickly, and produce higher peak concentrations, than predicted by homogeneous, isotropic diffusion models.

Field‐Scale Sulfur Hexafluoride Tracer Experiment to Understand Long Distance Gas Transport in the Deep Unsaturated Zone

A natural gradient SF6 tracer experiment provided an unprecedented evaluation of long distance gas transport in the deep unsaturated zone (UZ) under controlled (known) conditions. The field‐scale gas tracer test in the 110‐m‐thick UZ was conducted at the U.S. Geological Survey's Amargosa Desert Research Site (ADRS) in southwestern Nevada. A history of anomalous (theoretically unexpected) contaminant gas transport observed at the ADRS, next to the first commercial low‐level radioactive waste disposal facility in the United States, provided motivation for the SF6 tracer study. Tracer was injected into a deep UZ borehole at depths of 15 and 48 m, and plume migration was observed in a monitoring borehole 9 m away at various depths (0.5–109 m) over the course of 1 yr. Tracer results yielded useful information about gas transport as applicable to the spatial scales of interest for off‐site contaminant transport in arid unsaturated zones. Modeling gas diffusion with standard empirical expressions reasonably explained SF6 plume migration, but tended to underpredict peak concentrations for the field‐scale experiment given previously determined porosity information. Despite some discrepancies between observations and model results, rapid SF6 gas transport commensurate with previous contaminant migration was not observed. The results provide ancillary support for the concept that apparent anomalies in historic transport behavior at the ADRS are the result of factors other than nonreactive gas transport properties or processes currently in effect in the undisturbed UZ.

Insights gained from geochemical studies in the Waterloo Moraine: Indications and implications for anthropogenic loading

The Waterloo Moraine is a historic and continuing groundwater resource for the growing Region of Waterloo in southern Ontario, Canada, providing hundreds of thousands of litres of fresh water daily to the regional water supply. The complex stratigraphy of the Moraine is characterized by a system of multiple aquifers separated by regionally discontinuous aquitards, complicating resource assessments. The Region’s extensive monitoring and production well network has been sampled repeatedly over the last four decades, providing insights into how long-term pumping and urban development have affected aquifer geochemistry. Various environmental tracers were used to identify groundwater residence time and mixing including 3H, 3H/3He, CFCs, SF6, 14C, and 4He. Studies on the Waterloo Moraine provided important contributions to the application of these environmental tracers by validating unsaturated zone transport theory and documenting the effects of localized atmospheric CFC contamination. These studies have also provided important information about groundwater flow and mixing in the Moraine. Both historical and recent hydrogeological studies identified urban applications of road deicers (salt) and rural agricultural practices (nitrate loadings) as potential problems for future use of groundwater as a domestic water supply for the region, with Cl and NO3-N concentrations in some areas approaching 900 and 17 mg/L, respectively. The effects of these urban and rural contaminant issues are explored, and efforts to mitigate and reduce the impact, including the effects of best (beneficial) management practices, are discussed. The potential impacts of these non-point source contaminants on the public water supply should be considered when planning new residential developments.

Electrical resistivity tomography as monitoring tool for unsaturated zone transport: an example of preferential transport of deicing chemicals

Non-invasive spatially resolved monitoring techniques may hold the key to observe heterogeneous flow and transport behavior of contaminants in soils. In this study, time-lapse electrical resistivity tomography (ERT) was employed during an infiltration experiment with deicing chemical in a small field lysimeter. Deicing chemicals like potassium formate, which frequently impact soils on airport sites, were infiltrated during snow melt. Chemical composition of seepage water and the electrical response was recorded over the spring period 2010. Time-lapse electrical resistivity tomographs are able to show the infiltration of the melt water loaded with ionic constituents of deicing chemicals and their degradation product hydrogen carbonate. The tomographs indicate early breakthrough behavior in parts of the profile. Groundtruthing with pore fluid conductivity and water content variations shows disagreement between expected and observed bulk conductivity. This was attributed to the different sampling volume of traditional methods and ERT due to a considerable fraction of immobile water in the soil. The results show that ERT can be used as a soil monitoring tool on airport sites if assisted by common soil monitoring techniques.

Links between climate change, water-table depth, and water chemistry in a mineralized mountain watershed

Recent studies suggest that climate change is causing rising solute concentrations in mountain lakes and streams. These changes may be more pronounced in mineralized watersheds due to the sensitivity of sulfide weathering to changes in subsurface oxygen transport. Specific causal mechanisms linking climate change and accelerated weathering rates have been proposed, but in general remain entirely hypothetical. For mineralized watersheds, a favored hypothesis is that falling water tables caused by declining recharge rates allow an increasing volume of sulfide-bearing rock to become exposed to air, thus oxygen. Here, we test the hypothesis that falling water tables are the primary cause of an increase in metals and SO4 (100–400%) observed since 1980 in the Upper Snake River (USR), Colorado. The USR drains an alpine watershed geologically and climatologically representative of many others in mineralized areas of the western U.S. Hydrologic and chemical data collected from 2005 to 2011 in a deep monitoring well (WP1) at the top of the USR watershed are utilized. During this period, both water table depths and groundwater SO4 concentrations have generally increased in the well. A numerical model was constructed using TOUGHREACT that simulates pyrite oxidation near WP1, including groundwater flow and oxygen transport in both saturated and unsaturated zones. The modeling suggests that a falling water table could produce an increase in metals and SO4 of a magnitude similar to that observed in the USR (up to 300%). Future water table declines may produce limited increases in sulfide weathering high in the watershed because of the water table dropping below the depth of oxygen penetration, but may continue to enhance sulfide weathering lower in the watershed where water tables are shallower. Advective air (oxygen) transport in the unsaturated zone caused by seasonally variable recharge and associated water table fluctuations was found to have little influence on pyrite oxidation rates near WP1. However, this mechanism could be important in the case of a shallow dynamic water table and more abundant/reactive sulfides in the shallow subsurface. Data from WP1 and numerical modeling results are thus consistent with the falling water table hypothesis, and illustrate fundamental processes linking climate and sulfide weathering in mineralized watersheds.

Comparisons of diffusive and advective fluxes of gas phase volatile organic compounds (VOCs) in unsaturated zones under natural conditions

Diffusive flux is traditionally treated as the dominant mechanism of gas transport in unsaturated zones under natural conditions, and advective flux is usually neglected. However, some researchers have found that pressure-driven and density-driven advective flux may also be significant under certain conditions. This article compares the diffusive, pressure-driven and density-driven advective fluxes of gas phase volatile organic compound (VOCs) in unsaturated zones under various natural conditions. The presence of a less or more permeable layer at ground surface in a layered unsaturated zone is investigated for its impact on the net contribution of advective and diffusive fluxes. Results show although the transient advective flux can be greater than the diffusive flux, under most of the field conditions the net contribution of the advective flux is one to three orders of magnitude less than the diffusive flux, and the influence of the density-driven flux is undetectable. The advective flux contributes comparably with the diffusive flux only when the gas-filled porosity is less than 0.05. The presence of a less permeable layer at ground surface slightly increases the total flux in the underlying layer, while the presence of a more permeable layer at ground surface significantly increases the total flux in it. When the magnitude of water table fluctuation is less than 1cm, and the period is greater than 0.5day, the fluctuation of the water table can be simulated by fixing the water table position and setting a fluctuating moving velocity at the water table.

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).

Radionuclide Transport in the Unsaturated Zone at Yucca Mountain, Nevada

This paper presents a comprehensive treatment of large‐scale unsaturated zone (UZ) transport below the repository at the Yucca Mountain, Nevada, reporting on studies investigating the UZ as a component of the Lower Natural Barrier to radionuclide transport from the site. Regardless of the disposition of Yucca Mountain as the site for a nuclear waste repository, these studies and models of radionuclide transport through the UZ have relevance to the study of contaminant transport through deep vadose zones in general. For many radionuclides, including most strongly sorbing species, the Yucca Mountain UZ was found to prevent or substantially reduce the rate of movement of radionuclides to the accessible environment. The UZ radionuclide transport models developed for Yucca Mountain simulate the features and processes that influence the capability of the UZ below the repository to reduce the movement of radionuclides. Transport of radionuclides away from the repository depends on the rate and pathways of flow within the UZ, the retardation of radionuclides, and the rate of transport by colloids. Several basic processes affect radionuclide transport in the UZ, including the quantity of recharge and deep percolation of water, flow partitioning between fractures and rock matrix, diffusion of radionuclides from water flowing in fractures into the pores of the rock matrix, sorption of radionuclides onto rock or mineral surfaces, and colloid filtration. This paper summarizes the treatment of these and other processes, describing the experimental and observational basis for the models, conceptualizations, and predictions, and then presents some representative results.

Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship.

Quasi 3D modeling of water flow in vadose zone and groundwater

Summary The complexity of subsurface flow systems calls for a variety of concepts leading to the multiplicity of simplified flow models. One habitual simplification is based on the assumption that lateral flow and transport in unsaturated zone are not significant unless the capillary fringe is involved. In such cases the flow and transport in the unsaturated zone above groundwater level can be simulated as a 1D phenomenon, whereas the flow and transport through groundwater are viewed as 2D or 3D phenomena. A new approach for a numerical scheme for 3D variably saturated flow using quasi 3D Richards’ equation and finite difference scheme is presented. The corresponding numerical algorithm and the QUASI-3D computer code were developed. Results of the groundwater level simulations were compared with transient laboratory experimental data for 2D data constant–flux infiltration, quasi-3D HYDRUS-MODFLOW numerical model and a FULL-3D numerical model using Richards’ equation. Hypothetical 3D examples of infiltration, pumping and groundwater mound dissipation for different spatial–time scales are presented. Water flow simulation for the Alto Piura aquifer (Peru) demonstrates the QUASI-3D model application at the regional scale. Computationally the QUASI-3D code was found to be more efficient by an order of 10–300%, while being accurate with respect to the benchmark fully 3D variable saturation code, when the capillary fringe was considered.

Evaluation of Land Surface Scheme SABAE-HW in Simulating Snow Depth, Soil Temperature and Soil Moisture within the BOREAS Site, Saskatchewan

The Soil Atmosphere Boundary, Accurate Evaluation of Heat and Water (SABAE-HW) model is a multilayered, one-dimensional, physically based version of the Canadian Land Surface Scheme (CLASS) and uses the same methodologies as CLASS, version 2.6. SABAE provides an improved interface for groundwater modelling to simulate soil moisture, soil temperature, energy fluxes and snow depth for a wide range of soil and vegetation. This paper reports the results of the first field comparison of SABAE-HW using an extensive ten-year dataset from the Boreal Ecosystem Atmosphere Study (BOREAS) and the Boreal Ecosystem Research and Monitoring Sites (BERMS) project, an area in central Saskatchewan, Canada, rich in terms of hydrological and meteorological data. The model is also independently tested and verified with the Simultaneous Heat and Water (SHAW) model, which is an unsaturated-zone transport model. Two boundary conditions are considered at the bottom of the soil profile: a water table boundary condition and a unit g...

AN ANALYSIS OF THE EFFECT OF HYDRAULIC PARAMETERS ON RADIONUCLIDE MIGRATION IN AN UNSATURATED ZONE

A One-Dimensional Water Flow and Contaminant Transport in Unsaturated Zone (FTUNS) code has been developed in order to interpret radionuclide migration in an unsaturated zone. The pore-size distribution index (n) and the inverse of the air-entry value (α) for an unsaturated zone were measured by KS M ISO 11275 method. The hydraulic parameters of the unsaturated soil are investigated by using soil from around a nuclear facility in Korea. The effect of hydraulic parameters on radionuclide migration in an unsaturated zone has been analyzed. The higher the value of the n-factor, the more the cobalt concentration was condensed. The larger the value of α-factor, the faster the migration of cobalt was and the more aggregative the cobalt concentration was. Also, it was found that an effect on contaminant migration due to the pore-size distribution index (n) and the inverse of the air-entry value (α) was minute. Meanwhile, migrations of cobalt and cesium are in inverse proportion to the Freundich isotherm coefficient. That is to say, the migration velocity of cobalt was about 8.35 times that of cesium. It was conclusively demonstrated that the Freundich isotherm coefficient was the most important factor for contaminant migration.

Regional transport modelling for nitrate trend assessment and forecasting in a chalk aquifer

Regional degradation of groundwater resources by nitrate has become one of the main challenges for water managers worldwide. Regulations have been defined to reverse observed nitrate trends in groundwater bodies, such as the Water Framework Directive and the Groundwater Daughter Directive in the European Union. In such a context, one of the main challenges remains to develop efficient approaches for groundwater quality assessment at regional scale, including quantitative numerical modelling, as a decision support for groundwater management. A new approach combining the use of environmental tracers and the innovative ‘Hybrid Finite Element Mixing Cell’ (HFEMC) modelling technique is developed to study and forecast the groundwater quality at the regional scale, with an application to a regional chalk aquifer in the Geer basin in Belgium. Tritium data and nitrate time series are used to produce a conceptual model for regional groundwater flow and contaminant transport in the combined unsaturated and saturated zones of the chalk aquifer. This shows that the spatial distribution of the contamination in the Geer basin is essentially linked to the hydrodynamic conditions prevailing in the basin, more precisely to groundwater age and mixing and not to the spatial patterns of land use or local hydrodispersive processes. A three-dimensional regional scale groundwater flow and solute transport model is developed. It is able to reproduce the spatial patterns of tritium and nitrate and the observed nitrate trends in the chalk aquifer and it is used to predict the evolution of nitrate concentrations in the basin. The modelling application shows that the global inertia of groundwater quality is strong in the basin and trend reversal is not expected to occur before the 2015 deadline fixed by the European Water Framework Directive. The expected time required for trend reversal ranges between 5 and more than 50 years, depending on the location in the basin and the expected reduction in nitrate application. To reach a good chemical status, nitrate concentrations in the infiltrating water should be reduced as soon as possible below 50 mg/l; however, even in that case, more than 50 years is needed to fully reverse upward trends.