Groundwater flow modeling of periods with periglacial and glacial climate conditions for the safety assessment of the proposed high-level nuclear waste repository site at Forsmark,Sweden

The impact of periglacial and glacial climate conditions on groundwater flow in fractured crystalline rock is studied by means of groundwater flow modeling of the Forsmark site, which was recently proposed as a repository site for the disposal of spent high-level nuclear fuel in Sweden. The employed model uses a thermal-hydraulically coupled approach for permafrost modeling and discusses changes in groundwater flow implied by the climate conditions found over northern Europe at different times during the last glacial cycle (Weichselian glaciation). It is concluded that discharge of particles released at repository depth occurs very close to the ice-sheet margin in the absence of permafrost. If permafrost is included, the greater part discharges into taliks in the periglacial area. During a glacial cycle, hydraulic gradients at repository depth reach their maximum values when the ice-sheet margin passes over the site; at this time, also, the interface between fresh and saline waters is distorted the most. The combined effect of advances and retreats during several glaciations has not been studied in the present work; however, the results indicate that hydrochemical conditions at depth in the groundwater flow model are almost restored after a single event of ice-sheet advance and retreat.     

Geoscientific R&D for high level radioactive waste disposal in Sweden — current status and future plans

SKB (Svensk Kärnbränslehantering AB) is responsible for all handling, transport and storage of the nuclear wastes outside the Swedish nuclear power stations. According to Swedish law, SKB is responsible for an R&D-programme needed to take care of the radwastes. The programme comprises, among others, a general supportive geo-scientific R&D and the Äspö Hard Rock Laboratory (HRL) for more in-situ specific tasks.

Sweden is geologically located in the Fennoscandian shield which is dominated by gneisses and granitoids of Precambrian age. The Swedish reference repository concept thus considers an excavated vault at ca. 500 m depth in crystalline rocks. In this concept (KBS-3), copper canisters with high level waste will be emplaced in deposition holes from a system of tunnels. Blocks of highly compacted swelling bentonite clay are placed in the holes leaving ample space for the canisters. At the final closure of the repository, the galleries are backfilled with a mixture of sand and bentonite. This repository design aims to make the disposal system as redundant as possible. Although the KBS-3 concept is the reference concept, alternative concepts and/or repository lay-outs are also studied. The main alternative, currently under development at SKB, is disposal in boreholes with depths of 4–5 km. The geoscientific research will to a great extent be guided by the demands posed by the performance and safety assessments, as well as the constuctability issues. Some main functions of the geological barrier are fundamental for the long-term safety of a repository. These are: bedrock mechanical stability, a chemically stable environment as well as a slow and stable groundwater flux. The main time-table for the final disposal of long-lived radioactive waste in Sweden foresees the final selection of the disposal system and site during the beginning of next decade.     

Effect of transport-pathway simplifications on projected releases of radionuclides from a nuclear waste repository (Sweden)

The Swedish Nuclear Fuel and Waste Management Company has recently submitted an application for a license to construct a final repository for spent nuclear fuel, at approximately 500?m depth in crystalline bedrock. Migration pathways through the geosphere barrier are geometrically complex, with segments in fractured rock, deformation zones, backfilled tunnels, and near-surface soils. Several simplifications of these complex migration pathways were used in the assessments of repository performance that supported the license application. Specifically, in the geosphere transport calculations, radionuclide transport in soils and tunnels was neglected, and deformation zones were assumed to have transport characteristics of fractured rock. The effects of these simplifications on the projected performance of the geosphere barrier system are addressed. Geosphere performance is shown to be sensitive to how transport characteristics of deformation zones are conceptualized and incorporated into the model. Incorporation of advective groundwater travel time within backfilled tunnels reduces radiological dose from non-sorbing radionuclides such as I-129, while sorption in near-surface soils reduces radiological doses from sorbing radionuclides such as Ra-226. These results help quantify the degree to which geosphere performance was pessimistically assessed, and provide some guidance on how future studies to reduce uncertainty in geosphere performance may be focused.     

Comparison of vertical and horizontal deposition hole concept for disposal of radioactive waste based on rock mechanical in situ stress-strength analyses

In Finland, a spent fuel repository is being planned in hard, crystalline rock to the depth of 300–800 m. The repository will consist of parallel deposition tunnels which are interconnected by a central tunnel. The spent fuel canisters will be placed in holes drilled according to the tunnel floor—vertical deposition hole concept or in the tunnel wall-horizontal deposition hole concept. Based on the numerical stress-strength analyses, both concepts are found to be feasible. However, when the maximum horizontal stress is perpendicular to the axis of the deposition hole, the rock in the vicinity of the horizontal deposition hole seems to have more stable behavior than in the case of the vertical deposition hole.     

An empirical probabilistic approach for constraining the uncertainty of long-term solute transport predictions in fractured rock using in situ tracer experiments

Although radionuclide tracer tests have been carried out for over 30  years, the role of tracer tests in radioactive waste repository performance assessment (PA) has been questioned due to the differences between the time scales for tracer tests and PA. The possibility of using in situ tracer tests to constrain PA time scale (over 10,000  years) solute transport has been demonstrated using a systematic “microstructural model” approach to define advective and retentive materials. A series of simulations were conducted of the “TRUE-1” sorbing solute transport experiments in a well-characterized block of fractured granite. These experiments were then used to constrain the uncertainty of long-term transport on the same pathways, using a gradient several orders of magnitude smaller. The comparison of uncertainty of this long-term transport, with and without this conditioning step, provided a measure of the ability of tracer tests to reduce PA time-scale uncertainty. Although this approach for quantifying uncertainty reduction is somewhat empirical, it does indicate the potential usefulness of tracer experiments in reducing the uncertainty of the key PA time-scale transport parameters such as the flow wetted surface, provided that immobile zone properties such as sorption (Kd), porosity and diffusivity are available.     

Reactive transport modeling of the interaction between a high-pH plume and a fractured marl: the case of Wellenberg

In the context of the proposed low- and intermediate-level radioactive waste repository at Wellenberg (Switzerland), calculations simulating the interaction between hyperalkaline solutions and a fractured marl, at 25 °C, have been performed. The aim of these calculations is to evaluate the possible effects of mineral dissolution and precipitation on porosity and permeability changes in such a fractured marl, and their impact on repository performance. Solute transport and chemical reaction are considered in both a high-permeability zone (fracture), where advection is important, and the wall rock, where diffusion is the dominant transport mechanism. The mineral reactions are promoted by the interaction between hyperalkaline solutions derived from the degradation of cement (a major component of the engineered barrier system in the repository) and the host rock. Both diffusive/dispersive and advective solute transport are taken into account in the calculations. Mineral reactions are described by kinetic rate laws. The fluid flow system under consideration is a two-dimensional porous medium (marl, 1% porosity), with a high-permeability zone simulating a fracture (10% porosity) crossing the domain. The dimensions of the domain are 6 m per 1 m, and the fracture width is 10 cm. The fluid flow field is updated during the course of the simulations. Permeabilities are updated according to Kozeny's equation. The composition of the solutions entering the domain is derived from modeling studies of the degradation of cement under the conditions at the proposed underground repository at Wellenberg. Two different cases have been considered in the calculations. These 2 cases are representative of 2 different stages in the process of degradation of cement (pH 13.5 and pH 12.5). In both cases, the flow velocity in the fracture diminishes with time, due to a decrease in porosity. This decrease in porosity is caused by the precipitation of calcite (replacement of dolomite by calcite) and other secondary minerals (brucite, sepiolite, analcime, natrolite, tobermorite). However, the decrease in porosity and flow velocity is much more pronounced in the lower pH case. The extent of the zone of mineral alteration along the fracture is also much more limited in the lower pH case. The reduction of porosity in the fractures would be highly beneficial for repository performance, since it would mean that the solutions coming from the repository and potentially carrying radionuclides in solution would have to flow through low-conductive rock before they would be able to get to higher-conductive features. The biggest uncertainty in the reaction rates used in the calculations arises from the surface areas of the primary minerals. Additional calculations making use of smaller surface areas have also been performed. The results show that the smaller surface areas (and therefore smaller reaction rates) result in a smaller reactivity of the system and smaller porosity changes.     

Thermal constraints on disposal of heat-emitting waste in argillaceous rocks

Although the post-closure thermal phase of a repository containing heat-emitting radioactive wastes is likely to be of comparatively short duration when judged against performance assessment time-scales, the effects of temperature increase and associated thermal gradients within an argillaceous host-rock are likely to be wide-ranging and include: (a) geomechanical effects on the rock-mass and engineered structures (linings, emplacement boreholes, etc.), (b) perturbation of groundwater pressures and flow by ‘aquathermal mechanisms’, (c) possible near-field effects on mass transport properties and mechanisms (e.g. coupled flow phenomena), and (d) effects on the chemical evolution of the near-field environment. Maximum allowable temperature in a repository is likely to be constrained by host-rock and/or clay buffer/backfill thermal behaviour. Given the diverse range of the thermal responses, the identification of a primary constraint on near-field temperature and the specification of a maximum allowable temperature remains problematic.     

Mechanism and countermeasures of preceding tunnel distortion induced by succeeding EPBS tunnelling in close proximity

Twin tunnels are frequently used to address the increasing transportation demands in large cities. To ensure the safety of twin tunnels in close proximity, it is often necessary to take protective measures that have not been well studied. Field monitoring was conducted for a project of twin earth pressure balance shield (EPBS) tunnels in typical soft ground. The preceding tunnel was reinforced by various measures, including trailer bracing, compensation grouting, artificial freezing and scaffold bracing. The entire deformation of the reinforced tunnel was recorded during the succeeding tunnelling process. A three dimensional finite-element method (FEM) model was established to simulate the entire process of twin EPBS tunnelling, particularly the reinforcement measures. The computed deformations of the reinforced tunnel were consistent with the measured data. Furthermore, the stress history and pore pressure of the surrounding soil were analysed to investigate the deformation mechanism of the tunnel. Both the measured and computed results indicate that although the face pressure of the succeeding tunnel was smaller than the earth pressure at rest, the preceding tunnel could still experience an inward horizontal convergence and a deflection away from the succeeding tunnel. These distortion modes were caused by the squeezing effect of the horizontal soil arch in front of the succeeding tunnel face. Finally, convergence and deflection indices were proposed to quantify and assess the effectiveness of the reinforcement measures. The trailer bracing, as an “in-tunnel” reinforcement technique, was found to be the most effective method for controlling tunnel convergence. However, artificial freezing as an “out-tunnel” reinforcement technique led to the largest reductions in tunnel deflection. A combination of both “in-tunnel” and “out-tunnel” reinforcements was recommended.     

Applying numerical modeling for designing strategies of effective groundwater remediation

Selection of effective groundwater remediation scenarios is a complex issue that requires understanding of contaminants’ transport processes. The effectiveness of cleanup measures may be verified by fate and transport numerical modeling. The goal of this work was to present the usefulness of fate and transport modeling for planning, verification and fulfillment of effective groundwater remediation methods. Selection methodology was developed, which is based on results of numerical flow and transport modeling. A field site located in south-east Poland was selected as a case study, in which groundwater contamination of trichloroethene and tetrachloroethene was detected. The results indicated that “pump and treat” was the most effective among the studied remediation methods, followed by permeable reactive barrier and in situ chemical oxidation. Natural attenuation-based remediation was demonstrated to be the least suitable, as it requires the longest time to reach predefined remediation goals, principally due to low sorption capacity and unfavorable hydrogeochemical conditions for biodegradation. Fate and transport numerical modeling allowed simulating different remediation strategies, and thus the decision-making process was facilitated.     

Catchment-scale soil erosion and sediment yield simulation using a spatially distributed erosion model

Increasing rainfall intensity and frequency due to extreme climate change and haphazard land development are aggravating soil erosion problems in Korea. A quantitative estimate of the amount of sediment from the catchment is essential for soil and water conservation planning and management. Essential to catchment-scale soil erosion modeling is the ability to represent the fluvial transport system associated with the processes of detachment, transport, and deposition of soil particles due to rainfall and surface flow. This study applied a spatially distributed hydrologic model of rainfall–runoff–sediment yield simulation for flood events due to typhoons and then assessed the impact of topographic and climatic factors on erosion and deposition at a catchment scale. Measured versus predicted values of runoff and sediment discharge were acceptable in terms of applied model performance measures despite underestimation of simulated sediment loads near peak concentrations. Erosion occurred widely throughout the catchment, whereas deposition appeared near the channel network grid cells with a short hillslope flow path distance and gentle slope; the critical values of both topographic factors, providing only deposition, were observed at 3.5 (km) (hillslope flow path distance) and 0.2 (m/m) (local slope), respectively. In addition, spatially heterogeneous rainfall intensity, dependent on Thiessen polygons, led to spatially distinct net-erosion patterns; erosion increased gradually as rainfall amount increased, whereas deposition responded irregularly to variations in rainfall.     

动水条件下递增开放式涌水通道高效治理技术

陕西榆神矿区西湾露天煤矿东端帮一条爆破扩张裂隙沟通矿坑外侧烧变岩含水体,形成长度大于65 m的递增开放式涌水通道,水量稳定在400 m3/h左右,流速高达0.36 m/s,导致该矿面临生态、环保、经济、安全等一系列问题。为尽快封堵该涌水通道,化解以上问题,综合考量环保、工期、成本及施工条件等,针对裂隙主通道位置的确定提出逐步加密布置钻孔的快速探查方法,针对动水条件下递增开放式涌水通道常规灌注方案骨料留存率极低的问题,提出新型骨料充填技术。采用上述方法,成功探查出2个与主涌水通道连通的钻孔,明确通道位置;治理时,先在通道外端口布设“骨料拦阻网”,再采用镁渣核和钙质结核搭配,提高骨料留存率,实现通道的快速充填,经双液浆加固,水量剩余6.6 m3/h,仅用6 d即成功封堵该涌水通道。结果表明,本次提出的爆破裂隙探查方法和递增开放式涌水通道封堵方案科学、高效,该成功案例可为类似条件下涌水通道的探查与治理提供参考。      

Die Unterschiede in den Sedimentationsvorgängen zwischen dem Autunien und Saxonien von Mittel- und Westeuropa

After a comment on the stratigraphical classification and the subdivision of Permian strata, a review on their occurrences in Central and Western Europe including the Alps is given. Within this area, the strata can be summarized to a northern, a central, and a “Verrucano” province. The Autunian (Lower Rotliegend) has been typically developed in the central province where it is best preserved. It is characterized by many different sedimentation processes the most important of which are described. They lead to the conclusion that the sediment deposition mainly took place in fresh water lakes and through fluviatile transport. On the other hand, the deposition of the Saxonian (Upper Rotliegend) took place in a much drier climate. Finally, the differences between the two types of climate are described.     

Rough evaluation of the water-inflow discharge in abandoned mining tunnels using a simplified water balance model: the case of the Cogne iron mine (Aosta Valley, NW Italy)

One of the most complex hydrogeological problems in the design and maintenance of drainage systems in abandoned mining sites is quantifying the maximum water infiltration and, therefore, the amount that is potentially drainable by the tunnels. This problem is compounded when water-inflow data are limited or lacking. The aim of the study was to present a single but reliable model for making this evaluation; this model was applied to the case history of the abandoned Cogne iron ore mining complex (Western Alps, Aosta Valley Region, NW Italy). The study focused on quantifying the amount of water infiltrating into the mine drifts, using a water balance model in a Geographic Information System (GIS) environment. In the model, five different infiltration scenarios were calculated, including a detailed analysis of rainfall data, snow density and thickness (Snow Water Equivalent calculation), and melting periods. The maximum water discharge that could affect the mine tunnels was, therefore, determined under several scenarios of normal precipitation conditions and during heavy rainfall, including the case of the Cogne valley flood in October 2000, used as a reference for the limit conditions. Taking into account the various approximations considered, the results can be considered a good indication of the magnitude of the total amount of water that should be drained out through abandoned mine drifts and in the drainage network during implementation of final closure of the mine.     

Decadal-Timescale Estuarine Geomorphic Change Under Future Scenarios of Climate and Sediment Supply

Future estuarine geomorphic change, in response to climate change, sea-level rise, and watershed sediment supply, may govern ecological function, navigation, and water quality. We estimated geomorphic changes in Suisun Bay, CA, under four scenarios using a tidal-timescale hydrodynamic/sediment transport model. Computational expense and data needs were reduced using the morphological hydrograph concept and the morphological acceleration factor. The four scenarios included (1) present-day conditions; (2) sea-level rise and freshwater flow changes of 2030; (3) sea-level rise and decreased watershed sediment supply of 2030; and (4) sea-level rise, freshwater flow changes, and decreased watershed sediment supply of 2030. Sea-level rise increased water levels thereby reducing wave-induced bottom shear stress and sediment redistribution during the wind-wave season. Decreased watershed sediment supply reduced net deposition within the estuary, while minor changes in freshwater flow timing and magnitude induced the smallest overall effect. In all future scenarios, net deposition in the entire estuary and in the shallowest areas did not keep pace with sea-level rise, suggesting that intertidal and wetland areas may struggle to maintain elevation. Tidal-timescale simulations using future conditions were also used to infer changes in optical depth: though sea-level rise acts to decrease mean light irradiance, decreased suspended-sediment concentrations increase irradiance, yielding small changes in optical depth. The modeling results also assisted with the development of a dimensionless estuarine geomorphic number representing the ratio of potential sediment import forces to sediment export forces; we found the number to be linearly related to relative geomorphic change in Suisun Bay. The methods implemented here are widely applicable to evaluating future scenarios of estuarine change over decadal timescales.     

Controls on the alluviation of oxbow lakes by bed‐material load along the Sacramento River,California

Differences in the nature and quantity of sediment filling oxbow lakes have significant implications for the evolution of meandering rivers and the development of floodplains, influencing rates of meander migration and the valley width over which migration takes place. In an effort to identify the controls on the alluviation of oxbow lakes by coarse bed material, this study examined the sedimentary records stored within oxbow lakes of the Sacramento River of California, USA, and found that the volume of gravel in storage correlated negatively with the diversion angle separating flow between the river channel and the entrance into each lake. A method was devised for estimating the original channel bathymetry of the studied lakes and for modelling the hydraulic and sediment‐transport effects of the diversion angle within channels recently abandoned by meander cut‐off. The diversion angle determines the width of a flow separation within the abandoned‐channel entrance, reducing the discharge diverted from the river channel and thus limiting the ability of the abandoned channel to transport bed material. Aggradation rates are faster within entrances to abandoned channels with high diversion angles, resulting in the rapid isolation of lakes that store only a small volume of coarse‐grained sediment. Aggradation rates are slower within channel entrances where diversion angles are low, resulting in the slow transitioning of such channels into oxbow lakes with a larger and more extensive accumulation of coarse‐grained sediment. These findings compare well with observations in other natural settings and the mechanism which is described for the control of the diversion may explain why some oxbow lakes remain as open‐water environments for centuries, whereas others are filled completely within decades of cut‐off.     

Slope stabilization in difficult conditions: the case study of a debris slide in NW Italian Alps

A case study of a debris slide (estimated volume of about 35,000 m³) is described in this paper. This slide occurred in April 2009 in the North Western Italian Alps (Aosta valley) and damaged the SR25 road along the Valgrisenche valley. Ground investigations started with severe safety and logistic issues being posed. Given the need to open as soon as possible the road, the design of the landslide stabilization works was carried out using a “design as you go” approach. The stabilization measures were conceived to be flexible in order to allow for changes and integration during construction, in line with the progressive refinement of the geological–geotechnical slope model being developed. Back analysis by means of the limit equilibrium method (LEM) and the finite element method (FEM) was used. Groundwater level rise following heavy rainfall and spring snow melting was found to be the main cause of the debris slide. The stabilization works were designed by using both the LEM and FEM methods. The stability conditions of the engineered slope were assessed based on the available performance monitoring data.     

Numerical study of the performance of tunnel plugs

High axial flow rates within and along tunnels, excavated for deposition of high level nuclear waste, may increase the transport capacity of potential escape routes for radionuclides that have been released to the buffer surrounding the waste canisters, or to the tunnel backfill materials, from leaking fuel canisters. High flow rates may be found if the tunnel interior has a high permeability, either initially due to the composition of the backfill material or at later times due to degradation of low permeability backfill material components. If a disturbed rock zone, a DRZ, of increased permeability was created as a result of direct damage done to the wall region during tunnel excavation, this zone may contain additional possible escape routes. Low-permeability tunnel plugs, keyed into slots cut in the rock walls, have been suggested as means of reducing axial flow rates. This paper deals with different aspects of such tunnel-plug systems. A preliminary estimate of the potential for tunnel plugs to reduce axial flow rates is made using analytical expressions. A number of numerical techniques are employed to investigate the hydraulic, mechanical and hydromechanical performance. These include μFLOW (FEM flow calculations), FLAC (finite difference mechanical calculations) and 3DEC (distinct element mechanical calculations). The mechanical calculations concern the mechanical stability in the rock surrounding the slot and permeability changes in that rock caused by stress redistribution. The results show that the effects of tunnel plugs are very significant in most cases. A discussion on the validity of the results and the applicability of the used methods is provided.     

Modelling the evolution of hydrochemical conditions in the Fennoscandian Shield during Holocene time using multidisciplinary information

Three-dimensional, large-scale models for groundwater flow and solute transport are used for the low-temperature, fractured crystalline rock sites in Sweden that are being considered for the geological disposal of spent nuclear fuel. It has been suggested that comparisons between measured and simulated present-day hydrochemical data provide a means to constrain the complex influences of past climatic events and to improve the ability to understand the palaeohydrogeological evolution of the physical system studied. Here the authors demonstrate how the integration of multidisciplinary data and models from one of the sites in Sweden (Forsmark) can aid the appraisal of the hydrochemical conditions at 8000 BC, which is the selected starting point for the palaeohydrogeological modelling of the hydrochemical conditions in the Fennoscandian Shield during the Holocene (last 10 ka). Since a firm understanding of the evolution of the hydrochemical conditions is important for the long-term safety assessment, recognition of the initial hydrochemical conditions is essential for the overall build-up of confidence in the modelling process.