Bertel Nilsson

Spatially varying hydraulic and solute transport characteristics of a fractured till determined by field tracer tests, Funen, Denmark

A natural gradient tracer test that simulated one-dimensional flow was conducted in a large (4 m × 4.8 m surface area) isolated block of fractured till in Funen, Denmark, to determine both bulk and spatially distributed hydraulic properties. Hydraulic conductivity (based on free-flow tests) was higher in the upper 2.5 m (5.81 × 10−5 m s−1) than in the upper 4 m (8.00 × 10−6 m s−1) of till, indicating the more highly fractured and bioturbated nature of the upper part of the till. Chloride was introduced at the surface of the block as NaCl solution (490 mg L−1 Cl−), and breakthrough of Cl− was monitored in 13 horizontal screens at depths of 2.5 and 4 m. Arrival times of half of the input Cl− concentrations (C/c0=0.5) ranged about sevenfold among various screens at each sampling depth; 5–35 hours at 2.5 m and 11–68 hours at 4 m. Chloride breakthrough data for the two depths were fitted to both an equivalent porous media (EPM) model, based on the convective-dispersive equation, and a discrete fracture model. These results indicate differences in transport characteristics of the till with depth. At the 2.5 m depth the EPM model fit the initial part of the breakthrough curves quite well for most samplers but did not predict the long tails adequately. The converse was generally true for the discrete fracture model. Thus it appears that transport in the upper 2.5 m of till is controlled by a complex interconnected system of fractures and macropores. Field measurements of fractures and macropores at this site substantiate these results. At the 4 m depth the discrete fracture model described solute transport in all samplers better than the EPM model, indicating that vertical fracture flow dominates in the lower portion of the block. Fracture apertures and spacings were estimated based on the cubic law and fitting the breakthrough data to the discrete fracture model. Fitted fracture spacings at the 4 m depth (0.025–0.08 m) were slightly smaller than the minimum spacing (0.09 m) mapped in the field. Calculated apertures at the 4 m depth were large compared with other studies and ranged from 65 to 127 μm (average value of 94 μm).

Nitrate reduction in geologically heterogeneous catchments--a framework for assessing the scale of predictive capability of hydrological models.

In order to fulfil the requirements of the EU Water Framework Directive nitrate load from agricultural areas to surface water in Denmark needs to be reduced by about 40%. The regulations imposed until now have been uniform, i.e. the same restrictions for all areas independent of the subsurface conditions. Studies have shown that on a national basis about 2/3 of the nitrate leaching from the root zone is reduced naturally, through denitrification, in the subsurface before reaching the streams. Therefore, it is more cost-effective to identify robust areas, where nitrate leaching through the root zone is reduced in the saturated zone before reaching the streams, and vulnerable areas, where no subsurface reduction takes place, and then only impose regulations/restrictions on the vulnerable areas. Distributed hydrological models can make predictions at grid scale, i.e. at much smaller scale than the entire catchment. However, as distributed models often do not include local scale hydrogeological heterogeneities, they are typically not able to make accurate predictions at scales smaller than they are calibrated. We present a framework for assessing nitrate reduction in the subsurface and for assessing at which spatial scales modelling tools have predictive capabilities. A new instrument has been developed for airborne geophysical measurements, Mini-SkyTEM, dedicated to identifying geological structures and heterogeneities with horizontal and lateral resolutions of 30-50 m and 2m, respectively, in the upper 30 m. The geological heterogeneity and uncertainty are further analysed by use of the geostatistical software TProGS by generating stochastic geological realisations that are soft conditioned against the geophysical data. Finally, the flow paths within the catchment are simulated by use of the MIKE SHE hydrological modelling system for each of the geological models generated by TProGS and the prediction uncertainty is characterised by the variance between the predictions of the different models.

Transport and biodegradation of creosote compounds in clayey till, a field experiment

The transport and biodegradation of 12 organic compounds (toluene, phenol, o-cresol, 2,6-, 3,5-dimethylphenol, naphthalene, 1-methylnaphthalene, benzothiophene, dibenzofuran, indole, acridine, and quinoline) were studied at a field site located on the island of Funen, Denmark, where a clayey till 10–15 m deep overlies a sandy aquifer. The upper 4.8 m of till is highly fractured and the upper 2.5 m contains numerous root and worm holes. A 1.5–2 m thick sand lens is encountered within the till at a depth of 4.8 m. Sampling points were installed at depths of 2.5 m, 4 m, and in the sand lens (5.5 m) to monitor the downward migration of a chloride tracer and the organic compounds. Water containing organic compounds and chloride was infiltrated into a 4 m×4.8 m basin at a rate of 8.8 m3 day−1 for 7 days. The mass of naphthalene relative to chloride was 0.39–0.98 for the sampling points located at a depth of 2.5 m, 0.11–0.61 for the sampling points located at a depth of 4 m, and 0–0.02 for the sampling points located in the sand lens. A similar pattern was observed for eight organic compounds for which reliable results were obtained (toluene, phenol, o-cresol, 2,6-, 3,5-dimethylphenol, 1-methylnaphthalene, benzothiophene, and quinoline). This shows that the organic compounds were attenuated during the downward migration through the till despite the high infiltration rate. The attenuation process may be attributed to biodegradation.

Modeling fine-scale geological heterogeneity-examples of sand lenses in tills

Sand lenses at various spatial scales are recognized to add heterogeneity to glacial sediments. They have high hydraulic conductivities relative to the surrounding till matrix and may affect the advective transport of water and contaminants in clayey till settings. Sand lenses were investigated on till outcrops producing binary images of geological cross-sections capturing the size, shape and distribution of individual features. Sand lenses occur as elongated, anisotropic geobodies that vary in size and extent. Besides, sand lenses show strong non-stationary patterns on section images that hamper subsequent simulation. Transition probability (TP) and multiple-point statistics (MPS) were employed to simulate sand lens heterogeneity. We used one cross-section to parameterize the spatial correlation and a second, parallel section as a reference: it allowed testing the quality of the simulations as a function of the amount of conditioning data under realistic conditions. The performance of the simulations was evaluated on the faithful reproduction of the specific geological structure caused by sand lenses. Multiple-point statistics offer a better reproduction of sand lens geometry. However, two-dimensional training images acquired by outcrop mapping are of limited use to generate three-dimensional realizations with MPS. One can use a technique that consists in splitting the 3D domain into a set of slices in various directions that are sequentially simulated and reassembled into a 3D block. The identification of flow paths through a network of elongated sand lenses and the impact on the equivalent permeability in tills are essential to perform solute transport modeling in the low-permeability sediments.

Field trials of active and multi-port sock samplers in gravel-packed wells

Four novel sampling systems were tested for their ability to obtain level-representative groundwater samples from two gravel-packed boreholes in different hydrogeological settings at Skrydstrup (Denmark) and Torgau (Germany). Samples obtained with the novel devices were compared with those obtained from nearby permanent dedicated reference systems. Results obtained using the novel samplers vary from excellent agreement to no agreement. Numerical analysis was carried out in an attempt to understand the variation in the sampling results. Investigation showed that a single perfect sampler does not exist for these sites. Level-determined groundwater samples in gravel-packed wells can be obtained using active as well as multi-port sock samplers. The choice of the sampler depends mainly on the specific hydrogeological conditions, well completion, the type of pollution and the time available for sampling.

Development and testing of active groundwater samplers

Abstract Active groundwater sampling techniques are methods where the aquifer is flushed by pumping. The methods developed and tested represent non-dedicated methods for use in existing water wells. This paper describes two different sampling techniques: the Separation Pumping Technique (SP) and the Packer Baffle System (PBS). The methodology and design of the two systems is presented and the operational application is demonstrated by examples from full-scale field experiments. The methods are validated and their sensitivity to the well construction and the hydrogeological environment is assessed based on numerical modelling and controlled laboratory experiments. Active groundwater sampling techniques can be used for remedial pumping optimization and in obtaining hydraulic data and represent a fast operational and reliable sampling tool, also under heterogeneous and low permeability conditions.

Combining steam injection with hydraulic fracturing for the in situ remediation of the unsaturated zone of a fractured soil polluted by jet fuel.

A steam injection pilot-scale experiment was performed on the unsaturated zone of a strongly heterogeneous fractured soil contaminated by jet fuel. Before the treatment, the soil was stimulated by creating sub-horizontal sand-filled hydraulic fractures at three depths. The steam was injected through one hydraulic fracture and gas/water/non-aqueous phase liquid (NAPL) was extracted from the remaining fractures by applying a vacuum to extraction wells. The injection strategy was designed to maximize the heat delivery over the entire cell (10 m × 10 m × 5 m). The soil temperature profile, the recovered NAPL, the extracted water, and the concentrations of volatile organic compounds (VOCs) in the gas phase were monitored during the field test. GC-MS chemical analyses of pre- and post-treatment soil samples allowed for the quantitative assessment of the remediation efficiency. The growth of the heat front followed the configuration of hydraulic fractures. The average concentration of total hydrocarbons (g/kg of soil) was reduced by ∼ 43% in the upper target zone (depth = 1.5-3.9 m) and by ∼ 72% over the entire zone (depth = 1.5-5.5 m). The total NAPL mass removal based on gas and liquid stream measurements and the free-NAPL product were almost 30% and 2%, respectively, of those estimated from chemical analyses of pre- and post-treatment soil samples. The dominant mechanisms of soil remediation was the vaporization of jet fuel compounds at temperatures lower than their normal boiling points (steam distillation) enhanced by the ventilation of porous matrix due to the forced convective flow of air. In addition, the significant reduction of the NAPL mass in the less-heated deeper zone may be attributed to the counter-current imbibition of condensed water from natural fractures into the porous matrix and the gravity drainage associated with seasonal fluctuations of the water table.

Comparison and evaluation of model structures for the simulation of pollution fluxes in a tile-drained river basin.

The European Union Water Framework Directive requires an integrated pollution prevention plan at the river basin level. Hydrological river basin modeling tools are therefore promising tools to support the quantification of pollution originating from different sources. A limited number of studies have reported on the use of these models to predict pollution fluxes in tile-drained basins. This study focused on evaluating different modeling tools and modeling concepts to quantify the flow and nitrate fluxes in the Odense River basin using DAISY-MIKE SHE (DMS) and the Soil and Water Assessment Tool (SWAT). The results show that SWAT accurately predicted flow for daily and monthly time steps, whereas simulation of nitrate fluxes were more accurate at a monthly time step. In comparison to the DMS model, which takes into account the uncertainty of soil hydraulic and slurry parameters, SWAT results for flow and nitrate fit well within the range of DMS simulated values in high-flow periods but were slightly lower in low-flow periods. Despite the similarities of simulated flow and nitrate fluxes at the basin outlet, the two models predicted very different separations into flow components (overland flow, tile drainage, and groundwater flow) as well as nitrate fluxes from flow components. It was concluded that the assessment on which the model provides a better representation of the reality in terms of flow paths should not only be based on standard statistical metrics for the entire river basin but also needs to consider additional data, field experiments, and opinions of field experts.

The role of groundwater for lake-water quality and quantification of N seepage

The heterogeneous nature of both groundwater discharge to a lake (inflow) and nitrate concentrations in groundwater can lead to significant errors in calculations of nutrient loading. Therefore, an integrated approach, combining groundwater flow and transport modelling with observed nitrate and ammonium groundwater concentrations, was used to estimate nitrate loading from a catchment via groundwater to an oligotrophic flow-through lake (Lake Hampen, Denmark). The transport model was calibrated against three vertical nitrate profiles from multi-level wells and 17 shallow wells bordering a crop field near the lake. Nitrate concentrations in groundwater discharging to the lake from the crop field were on average 70 times higher than in groundwater from forested areas. The crop field was responsible for 96% of the total nitrate loading (16.2 t NO3 /year) to the lake even though the field only covered 4.5% of the catchment area. Consequently, a small change in land use in the catchment will have a large effect on the lake nutrient balance and possible lake restoration. The study is the first known attempt to estimate the decrease of nitrate loading via groundwater to a seepage lake when an identified catchment source (a crop field) is removed.

Subsurface nitrate reduction under wetlands takes place in narrow superficial zones.

ABSTRACT This study aims to investigate the depth distribution of the Nitrate Reduction Potential (NRP) on a natural and a re-established wetland. The obtained NRP provides a valuable data of the driving factors affecting denitrification, the Dissimilatory Nitrate Reduction to Ammonium (DNRA) process and the performance of a re-established wetland. Intact soil cores were collected and divided in slices for the determination of Organic Matter (OM) through Loss of Ignition (LOI) as well as Dissolved Organic Carbon (DOC) and NRP spiking nitrate in batch tests. The Nitrate Reduction (NR) was fitted as a pseudo-first order rate constant (k) from where NRPs were obtained. NR took place in a narrow superficial zone showing a dropping natural logarithmic trend along depth. The main driving factor of denitrification, besides depth, was OM. Although, DOC and LOI could not express by themselves and absolute correlation with NRP, high amounts of DOC ensured enough quantity and quality of labile OM for NR. Besides, high concentration of LOI but a scarce abundance of DOC failed to drive NR. DNRA was only important in superficial samples with high contents of OM. Lastly, the high NRP of the re-established wetland confirms that wetlands can be restored satisfactorily.