Hydraulic conductivities of Boom Clay measured through various testing techniques in the laboratory exhibit similar K values in the order of 10 m/s. In situ measurements obtained from various scales at the HADES underground research facility yield K values that are very similar to values measured in the laboratory on samples of a few centimetres. Spatial analysis of K values across the Boom Clay at the Mol site reveals a typical profile with a very homogenous 61-m thick central part (the Putte and Terhagen Members), which is also the least permeable part of the Boom Clay. The geometric mean of the vertical (Kv) and horizontal (Kh) hydraulic conductivities for the Putte and Terhagen Members at the Mol site are 1.7×10 and 4.4×10 m/s, respectively, with a vertical anisotropy Kh/Kv of about 2.5. Higher K values, but still low (10 to 10 m/s), are observed in the more silty zones above and below. A regional analysis of vertical K variability of the Boom Clay in the northeast of Belgium confirms this typical K profile of the Boom Clay based on test results from other four regional boreholes. Furthermore, an increase in overall Boom Clay vertical hydraulic conductivity from the east towards the west in north-east Belgium was observed.
surements, in the framework of the environmental impact assessment of the near surface repository for low- and intermediate-level short-lived waste, realized by ONDRAF/NIRAS at Dessel, Belgium. Horizontal flow directions, horizontal and vertical gradients for the entire area of the groundwater model were estimated from measurements at shallow monitoring wells within the groundwater flow model domain, and compared to the flow directions and vertical gradients predicted by the model. For obtaining horizontal flow directions and gradients, triangulation of groundwater levels was performed for combinations of three neighboring hydraulic head observations in the same hydrogeological layer within the model. The simulated equivalents at the same monitoring wells were used to repeat the same methodology, and calculate flow direction components. This analysis was performed in SAGA GIS and was visualized through QGIS. Comparison of the flow directions and flow gradients obtained from measurements and simulations gives an indication on the model performance. The calculations were performed for three sandy hydrogeological units used in the model. A similar procedure was performed for the vertical hydraulic head gradients, where any combination of two hydraulic head observations at the same location but at different levels within the aquifer were used to validate the vertical gradients predicted by the model. Besides model validation on average hydraulic heads, the variability of flow direction and hydraulic gradients in time was checked, by using the actually measured monthly time series, to verify the applicability of the steady-state modelling approach. This basic assessment of flow directions and gradients using open source GIS can be used to identify potential areas of interest, were more detailed investigations would be recommended.
In the framework of the disposal of short-lived low- and intermediate-level radioactive waste in a near-surface disposal facility in Dessel (Belgium), extensive characterization of the hydraulic conductivity (K) in the shallow Neogene aquifer has been performed at a regional scale. In the last few years the small-scale heterogeneity has been additionally characterized by outcrop analogue, hydraulic direct push, and borehole core air permeameter studies. The gathered data now include a) more than 350 hydraulic conductivity measurements on samples from 8 cored boreholes, mostly reaching depths of 50 m and data at 2 m intervals, b) more than 5000 air permeability measurements on the same borehole cores, c) more than 250 cone penetration tests (CPTs) with depths down to 40 m and data at 2 cm intervals, d) over 100 dissipation tests performed during the CPT campaigns, e) 17 direct push injections loggings, 6 hydraulic profiling tool logs, and 6 direct push slug tests, f) several hundreds of air permeability measurements on outcrop analogues of the aquifer sediments, and g) numerous grain size analyses. The current study aims to quantify the heterogeneity of K from the centimetre- to the kilometre-scale and to check the compatibility of the spatial variability revealed by the different datasets. This is achieved through gathering all K values (either direct measurements, calibrated relative K values, or K estimates from secondary data), and the use of variography to quantify spatial variability in terms of two-points geostatistics. The results are discussed, and the main differences between the different data sources are explained. In a final step, different multi-scale variogram models are proposed for capturing the main characteristics of multi-scale variability within the shallow Neogene aquifer in Belgium.
Abstract. Lowland rivers and shallow aquifers are closely coupled and their interactions are crucial for maintaining healthy stream ecological functions. In order to explore river–aquifer interactions and lowland hydrological system in three Belgian catchments, we apply a combined approach of baseflow separation, impulse response modelling and time series analysis over a 30–year study period at catchment scale. Baseflow from hydrograph separation shows that the three catchments are groundwater-dominated. The recursive digital filter methods generate a smoother baseflow time series than the graphical methods, and yield more reliable results than the graphical ones. Impulse response modelling is applied with a two–step procedure. The first step where groundwater level response is modelled shows that groundwater level in shallow aquifers reacts fast to the system input, with most of the wells reaching their peak response during the first day. There is an overall trend of faster response time and higher response magnitude in the wet (October–March) than the dry (April–September) periods. The second step of baseflow response modelling shows that the system response is also fast and that simulated baseflow can capture some variations but not the peaks of the separated baseflow time series. The time series analysis indicates that components such as interflow and overland flow, contribute significantly to stream flow. They are somehow included as part of the separated baseflow, which is likely to be overestimated from hydrograph separation. The impulse response modelling approach from the groundwater flow perspective can be an optional method to estimate the baseflow, since it considers some level of the physical connection between river and aquifer in the subsurface. Further research is however recommended to improve the simulation, such as giving more weight to wells close to the river and adding more drainage dynamics to the model input.