An overview of Task 6 of the Äspö Task Force: modelling groundwater and solute transport: improved understanding of radionuclide transport in fractured rock

An overview is presented of a 4-year study by the Äspö Task Force on Modelling of Groundwater Flow and Transport of Solutes, whose primary aim was to build a bridge between the approaches used for site characterisation (SC) and performance assessment (PA) associated with nuclear waste repositories. Eleven modelling teams representing six national radioactive waste organisations participated in eight modelling exercises whose objectives were: to assess simplifications used in PA models; to determine how, and to what extent, experimental tracer and flow experiments can constrain the range of parameters used in PA models; to support the design of SC programmes to assure that the results have optimal value for PA calculations; and to improve the understanding of site-specific flow and transport behaviour at different scales using SC models. The modelling tasks were concerned with flow and transport through single and multiple near-planar features on SC and PA timescales, including the diffusion of solutes into multiple immobile zones adjacent to fracture surfaces. In general, tracer tests provide only limited quantitative constraints on retention parameter values relevant to PA but nevertheless provide insight about the flow and transport processes, which is a key element of the bridge between SC and PA.     

Verification of net primary production estimation method in the Mongolian Plateau using landsat ETM+data

We plan to estimate global net primary production (NPP) of vegetation using the Advanced Earth Observing Satellite-II (ADEOS-II) Global Imager (GLI) multi-spectral data. We derive an NPP estimation algorithm from ground measurement data on temperate plants in Japan. By the algorithm, we estimate NPP using a vegetation index based on pattern decomposition (VIPD) for the Mongolian Plateau. The VIPD is derived from Landsat ETM+multi-spectral data, and the resulting NPP estimation is compared with ground data measured in a semi-arid area of Mongolia. The NPP estimation derived from satellite remote sensing data agrees with the ground measurement data within the error range of 15% when all above-ground vegetation NPP is calculated for different vegetation classifications.     

Anomalous change in atmospheric radon concentration sourced from broad crustal deformation: A case study of the 1995 Kobe earthquake

Before the Kobe earthquake, an anomalous increase in atmospheric Rn concentration was observed. By separating the measured concentration of atmospheric Rn into three components according to the distance from the monitoring station, the variation of Rn exhalation rate can be estimated for the respective area using the daily minimum and maximum concentrations. The mean rate of Rn exhalation gradually increased in an area of 20 km around the monitoring station, becoming five times higher than normal in the period between October 1994 and the date of the earthquake. This area had a large co-seismic displacement of up to 30 cm, which roughly corresponds to the crustal strain of 10⁻⁶-order, and it is considered the main source for the atmospheric Rn prior to the Kobe earthquake. Analyses revealed that the pre-seismic change in the atmospheric Rn concentration exhibited an anomalous pattern which would yield information on the spatial distribution of the mechanical response of the ground.     

Weathering mechanisms and their effects on landsliding in pelitic schist

Physical, chemical, and mineralogical analyses of undisturbed drill cores of pelitic schist from a landslide area in Japan clarified the mechanisms of chemical weathering of pelitic schist. Oxidizing surface water percolates downward and reaches an oxidation front, where chlorite is altered to Al‐vermiculite, graphite and pyrite are oxidized and depleted, and goethite precipitates. Oxidation of pyrite also occurs just below the oxidation front, probably by ferric iron. Pyrite oxidation yields sulphuric acid, which penetrates further downward, interacting with and weakening the rocks. In addition to this chemical weakening, stress release and shearing along schistosities form an incipient shear zone, which propagates to a sliding zone that forms the rupture surface of a landslide. Once a sliding zone has developed, it inhibits downward groundwater flow across it because of its low permeability, slowing the downward propagation of the weathering zone until this filtration barrier is broken by landslide movement. Copyright © 2010 John Wiley & Sons, Ltd.     

Investigation of the sliding behavior between steel and mortar for seismic applications in structures

The friction developed between a steel base plate and a mortar base contributes shear resistance to the building system during a seismic event. In order to investigate the possible sliding behavior between the base plate and the mortar, a shake table study is undertaken using a large rigid mass supported by steel contact elements which rest on mortar surfaces connected to the shake table. Horizontal input accelerations are considered at various magnitudes and frequencies. The results provide a constant friction coefficient during sliding with an average value of approximately 0.78. A theoretical formulation of the friction behavior is also undertaken. The theoretical equations show that the sliding behavior is dependent on the ratio of the friction force to the input force. The addition of vertical accelerations to the system further complicates the sliding behavior as a result of the varying normal force. This results in a variable friction resistance which is a function of the amplitude, phase, and frequency of the horizontal and vertical input motions. In general, this study showed a consistent and reliable sliding behavior between steel and mortar. Copyright © 2009 John Wiley & Sons, Ltd.     

SHRIMP zircon and EPMA monazite dating of granitic rocks from the Maizuru terrane, southwest Japan: Correlation with East Asian Paleozoic terranes and geological implications

Abstract The Maizuru terrane, distributed in the Inner Zone of southwest Japan, is divided into three subzones (Northern, Central and Southern), each with distinct lithological associations. In clear contrast with the Southern zone consisting of the Yakuno ophiolite, the Northern zone is subdivided into the western and eastern bodies by a high-angle fault, recognized mainly by the presence of deformed granitic rocks and pelitic gneiss. This association suggests an affinity with a mature continental block; this is supported by the mode of occurrence, and petrological and isotopic data. Newly obtained sensitive high mass-resolution ion microprobe (SHRIMP) zircon U–Pb ages reveal the intrusion ages of 424 ± 16 and 405 ± 18 Ma (Siluro–Devonian) for the granites from the western body, and 249 ± 10 and 243 ± 19 Ma (Permo–Triassic) for the granodiorites from the eastern body. The granites in the western body also show inherited zircon ages of around 580 and 765 Ma. In addition, electron probe microanalysis (EPMA) monazite U–Th–total Pb dating gives around 475–460 Ma. The age of intrusion, inherited ages, mode of occurrence, and geological setting of the Siluro–Devonian granites of the Northern zone all show similarities with those of the Khanka Massif, southern Primoye, Russia, and the Hikami granitic rocks of the South Kitakami terrane, Northeast Japan. We propose that both the Siluro–Devonian and Permo–Triassic granitic rocks of the Northern zone are likely to have been juxtaposed through the Triassic–Late Jurassic dextral strike-slip movement, and to have originated from the Khanka Massif and the Hida terrane, respectively. This study strongly supports the importance of the strike-slip movement as a mechanism causing the structural rearrangement of the Paleozoic–Mesozoic terranes in the Japanese Islands, as well as in East Asia.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

Control of Phosphate Concentration through Adsorption and Desorption Processes in Groundwater and Seawater Mixing at Sandy Beaches in Tokyo Bay, Japan

Field observation was conducted to monitor phosphate concentrations in groundwater and seawater mixing at two sandy beaches in Futtsu and Miura in Tokyo Bay, Japan. Dissolved phosphate concentrations were measured along transects from fresh groundwater aquifer to seawater adjacent the beaches. The concentrations were often high (up to 46 µM) in fresh groundwater samples (Cl^– < 0.2). Coastal seawater, on the other hand, exhibited low phosphate concentrations (1.5 µM or less). Along the transects, phosphate generally displayed non-conservative behavior during mixing of fresh and saline waters in the aquifer; concentrations as high as 100 µM were found around the upper limit of seawater intrusion (Cl^– = 2). Laboratory experiments were executed to identify the processes that control the phosphate behavior in the mixing processes. The results revealed that adsorption-desorption processes by the aquifer sand particles could significantly control the phosphate concentrations in the groundwater. Furthermore, the adsorption and/or desorption was found to be a function of salinity; the equilibrium concentration of dissolved phosphate in slurry of sand and water was the highest in freshwater and decreased considerably in saline water. The extreme concentration of phosphate may be caused by release from sand particles coinciding with the rapid change in salinity with tide.     

Formation of thermohaline front by cooling of the sea surface and inflow of the fresh water

A numerical study of characteristics of the front formed by cooling of the sea surface and inflow of the fresh water is made within a vertical two-dimensional plane without the rotation of the earth. The convective adjustment technique is employed to parameterize the small scale convective overturning process.A sharply edged shelf causing the marked difference of the heat capacity of the water columns between over the shelf (50 m depth) and in the deeper region (150 m depth), is responsible for the formation of a thermal front over the shelf edge and two cell circulations of the same sense with each other adjoining at the frontal region.Inflow of the fresh water at the upper region of the coast (which is equivalently replaced by the outflow of the salt from the same region) adds a circulation of the opposite sense to the above ones, which transports the water of low salinity and low temperature offshore and forms a sharp front both of temperature and salinity with the offshore water of high salinity and high temperature. This thermohaline front or Oceanic front has a remarkable character that the horizontal density gradient is minimal at the front due to the counteracting contributions of temperature and salinity to density. Response of this front to the sudden ceases of the surface cooling and the fresh water supply, is also studied in order to understand its transient behaviors.