Hydrological control of water quality - Modelling base cation weathering and dynamics across heterogeneous boreal catchments.

Linking biogeochemical processes to water flow paths and solute travel times is important for understanding internal catchment functioning and control of water quality. Base cation weathering is a process closely linked to key factors affecting catchment functioning, including water pathways, soil contact time, and catchment characteristics, particularly in silicate-dominated areas. However, common process-based weathering models are often calibrated and applied for individual soil profiles, which can cause problems when trying to extrapolate results to catchment scale and assess consequences for stream water and groundwater quality. Therefore, in this work, base cation export was instead modelled using a fully calibrated 3D hydrological model (Mike SHE) of a boreal catchment, which was expanded by adding a relatively simple but still reasonably flexible and versatile weathering module including the base cations Na, K, Mg, and Ca. The results were evaluated using a comprehensive dataset of water chemistry from groundwater and stream water in 14 nested sub-catchments, representing different catchment sizes and catchment characteristics. The strongest correlations with annual and seasonal observations were found for Ca (r = 0.89-0.93, p < 0.05), Mg (r = 0.90-0.95, p < 0.05), and Na (r = 0.80-0.89, p < 0.05). These strong correlations suggest that catchment hydrology and landscape properties primarily control weathering rates and stream dynamics of these solutes. Furthermore, catchment export of Mg, Ca, and K was strongly connected to travel times of discharging stream water (r = 0.78-0.83). Conversely, increasing Na export was linked to a reduced areal proportion of mires (r = -0.79). The results suggest that a significant part (~45%) of the catchment stream export came from deep-soil weathering sources (>2.5 m). These results have implications for terrestrial and aquatic water quality assessments. If deep soils are present, focusing mainly on the shallow soil could lead to misrepresentation of base cation availability and the acidification sensitivity of groundwater and water recipients such as streams and lakes.

Microtomography-based Inter-Granular Network for the simulation of radionuclide diffusion and sorption in a granitic rock.

Field investigation studies, conducted in the context of safety analyses of deep geological repositories for nuclear waste, have pointed out that in fractured crystalline rocks sorbing radionuclides can diffuse surprisingly long distances deep into the intact rock matrix; i.e. much longer distances than those predicted by reactive transport models based on a homogeneous description of the properties of the rock matrix. Here, we focus on cesium diffusion and use detailed micro characterisation data, based on micro computed tomography, along with a grain-scale Inter-Granular Network model, to offer a plausible explanation for the anomalously long cesium penetration profiles observed in these in-situ experiments. The sparse distribution of chemically reactive grains (i.e. grains belonging to sorbing mineral phases) is shown to have a strong control on the diffusive patterns of sorbing radionuclides. The computed penetration profiles of cesium agree well with an analytical model based on two parallel diffusive pathways. This agreement, along with visual inspection of the spatial distribution of cesium concentration, indicates that for sorbing radionuclides the medium indeed behaves as a composite system, with most of the mass being retained close to the injection boundary and a non-negligible part diffusing faster along preferential diffusive pathways.

Identification and characterization of potential discharge areas for radionuclide transport by groundwater from a nuclear waste repository in Sweden.

This paper describes solute transport modeling carried out as a part of an assessment of the long-term radiological safety of a planned deep rock repository for spent nuclear fuel in Forsmark, Sweden. Specifically, it presents transport modeling performed to locate and describe discharge areas for groundwater potentially carrying radionuclides from the repository to the surface where man and the environment could be affected by the contamination. The modeling results show that topography to large extent determines the discharge locations. Present and future lake and wetland objects are central for the radionuclide transport and dose calculations in the safety assessment. Results of detailed transport modeling focusing on the regolith and the upper part of the rock indicate that the identification of discharge areas and objects considered in the safety assessment is robust in the sense that it does not change when a more detailed model representation is used.

Upscaling of hydraulic conductivity and telescopic mesh refinement.

Performance assessments of repositories for the underground disposal of nuclear fuel and waste include models of ground water flow and transport in the host rocks. Estimates of hydraulic conductivity, K, based on field measurements may require adjustment (upscaling) for use in numerical models, but the choice of upscaling approach can be complicated by the use of nested modeling, large-scale fracture zones, and a high degree of heterogeneity. Four approaches to upscaling K are examined using a reference case based on exhaustive site data and an application of nested modeling to evaluate performance assessment of a waste repository. The upscaling approaches are evaluated for their effects on the flow balance between nested modeling domains and on simple measures of repository performance. Of the upscaling approaches examined in this study, the greatest consistency of boundary flows was achieved using the observed scale dependence for the rock domains, measured values from the large-scale interference test for the conductor domain, and a semivariogram regularization based on the Moye model for packer test interpretation. Making the assumption that large fracture zones are two-dimensional media results in the greatest changes to the median of travel time and improves the flow balance between the nested models. The uncertainty of upscaling methods apparently has a small impact on median performance measures, but a significant impact on the variances and earliest arrival times.