Spatial measurements of stream baseflow,a relevant method for aquifer characterization and permeability evaluation. Application to a hard‐rock aquifer,the Oman ophiolite

The structure, functioning and hydrodynamic properties of aquifers can be determined from an analysis of the spatial variability of baseflow in the streams with which they are associated. Such analyses are based on simple low‐cost measurements. Through interpreting the hydrological profiles (Q = f(A)) it is possible to locate the aquifer(s) linked to the stream network and to determine the type of interrelated flow, i.e. whether the stream drains or feeds the aquifer. Using an analytical solution developed for situations with a positive linear relationship, i.e. where the baseflow increases linearly with increasing catchment size, it is also possible to estimate the permeability of the aquifer(s) concerned at catchment scale. Applied to the hard‐rock aquifers of the Oman ophiolite, this method shows that the ‘gabbro’ aquifer is more permeable than the ‘peridotite’ aquifer. As a consequence the streams drain the peridotites and ‘leak’ into the gabbro. The hydrological profiles within the peridotite are linear and positive, and indicate homogeneity in the hydrodynamic properties of these formations at the kilometre scale. The permeability of the peridotite is estimated at 5 · 10^?7 to 5 · 10^?8 m/s. Copyright © 2004 John Wiley & Sons, Ltd.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Late Quaternary submarine bedforms and ice‐sheet flow in Gerlache Strait and on the adjacent continental shelf,Antarctic Peninsula

Geophysical data from Gerlache Strait, Croker Passage, Bismarck Strait and the adjacent continental shelf reveal streamlined subglacial bedforms that were produced at the bed of the Antarctic Peninsula Ice Sheet (APIS) during the last glaciation. The spatial arrangement and orientation of these bedforms record the former drainage pattern and flow dynamics of an APIS outlet up‐flow, and feeding into, a palaeo‐ice stream in the Western Bransfield Basin. Evidence suggests that together, they represent a single ice‐flow system that drained the APIS during the last glaciation. The ice‐sheet outlet flowed north/northeastwards through Gerlache Strait and Croker Passage and converged with a second, more easterly ice‐flow tributary on the middle shelf to form the main palaeo‐ice stream. The dominance of drumlins with low elongation ratios suggests that ice‐sheet outlet draining through Gerlache Strait was comparatively slower than the main palaeo‐ice stream in the Western Bransfield Basin, although the low elongation ratios may also partly reflect the lack of sediment. Progressive elongation of drumlins further down‐flow indicates that the ice sheet accelerated through Croker Passage and the western tributary trough, and fed into the main zone of streaming flow in the Western Bransfield Basin. Topography would have exerted a strong control on the development of the palaeo‐ice stream system but subglacial geology may also have been significant given the transition from crystalline bedrock to sedimentary strata on the inner–mid‐shelf. In the broader context, the APIS was drained by a number of major fast‐flowing outlets through cross‐shelf troughs to the outer continental shelf during the last glaciation. Copyright © 2004 John Wiley & Sons, Ltd.     

Evolution of 2:1 layered silicates in low-grade metamorphosed Liassic shales of Central Switzerland

White mica from the Liassic black shales and slates in Central Switzerland was analysed by transmission electron microscopy (TEM) and electron microprobe to determine its textural and compositional evolution during very low-grade prograde metamorphism. Samples were studied from the diagenetic zone, anchizone and epizone (T ≈100°–450 °C). Phyllosilicate minerals analysed include illite/smectite (I/S), phengite, muscovite, brammallite, paragonite, margarite and glauconite. Textural evolution primarily is towards larger, more defect-free grains with compositions that approach those of their respective end-members. The smectite-to-illite transformation reduced the amounts of the exchange components SiK_?1Al_?1, MgSiAl_?2, and Fe³⁺Al_?1. These trends continue to a lesser degree in the anchizone and epizone. Correlations between the proportion of smectite in I/S and the composition of I/S indicate that smectite layers may contain a high layer charge. Illite in I/S bears a compositional resemblance to macrocrystalline phengite in some samples, but is different in others. Paragonite first appears in the upper diagenetic zone or lower anchizone as an interlayer-deficient brammallite, and it may be mixed with muscovite on the nanometre scale. Owing to the small calculated structure factor for paragonite-muscovite superstructures, conventional X-ray powder diffraction cannot distinguish between mixed-layer structures and a homogeneous compositionally intermediate solid solutions. However, indirect TEM evidence shows that irregularly shaped domains of Na- and K-rich mica exist below 10 nm. Subsequent coarsening of domains at higher grades produced discrete paragonite grains at the margins of muscovite crystals or in laths parallel to the basal plane of the host muscovite. Margarite appears in the epizone and follows a textural evolution similar to paragonite in that mixtures of margarite, paragonite, and muscovite may initially occur on the nanometre scale. However, no evidence of interlayer-poor margarite has been found.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

High-pressure eclogites in northern Jiangsu – southern Shandong province, eastern China

Abstract Eclogites are distributed for more than 500 km along a major tectonic boundary between the Sino-Korean and Yangtze cratons in central and eastern China. These eclogites usually have high-P assemblages including omphacite + kyanite and/or coesite (or its pseudomorph), and form a high-P eclogite terrane. They occur as isolated lenses or blocks 10 cm to 300 m long in gneisses (Type I), serpentinized garnet peridotites (Type II) and marbles (Type III). Type I eclogites were formed by prograde metamorphism, and their primary metamorphic mineral assemblage consists mainly of garnet [pyrope (Prp) = 15–40 mol%], omphacite [jadeite (Jd) = 34–64 mol%], pargasitic amphibole, kyanite, phengitic muscovite, zoisite, an SiO₂ phase, apatite, rutile and zircon. Type II eclogites characteristically contain no SiO₂ phase, and are divided into prograde eclogites and mantle-derived eclogites. The prograde eclogites of Type II are petrographically similar to Type I eclogites. The mantle-derived eclogites have high MgO/(FeO + Fe₂O₃) and Cr₂O₃ compositions in bulk rock and minerals, and consist mainly of pyrope-rich garnet (Prp = 48–60 mol%), sodic augite (Jd = 10–27 mol%) and rutile. Type III eclogites have an unusual mineral assemblage of grossular-rich (Grs = 57 mol%) garnet + omphacite (Jd = 30–34 mol%) + pargasite + rutile. Pargasitic and taramitic amphiboles, calcic plagioclase (An₆₈), epidote, zoisite, K-feldspar and paragonite occur as inclusions in garnet and omphacite in the prograde eclogites. This suggests that the prograde eclogites were formed by recrystallization of epidote amphibolite and/or amphibolite facies rocks with near-isothermal compression reflecting crustal thickening during continent–continent collision of late Proterozoic age. Equilibrium conditions of the prograde eclogites range from P > 26 kbar and T= 500–750°C in the western part to P > 28 kbar and T= 810–880°C in the eastern part of the high-P eclogite terrane. The prograde eclogites in the eastern part are considered to have been derived from a deeper position than those in the western part. Subsequent reactions, manifested by (1) narrow rims of sodic plagioclase or paragonite on kyanite and (2) symplectites between omphacite and quartz are interpreted as an effect of near-isothermal decompression during the retrograde stage. The conditions at which symplectites re-equilibrated tend to increase from west (P < 10 kbar and T < 580°C) to east (P > 9 kbar and T > 680°C). Equilibrium temperatures of Type II mantle-derived eclogites and Type III eclogite are 730–750°C and 680°C, respectively.     

Exposure dating and reinterpretation of coarse debris accumulations (‘rock glaciers’) in the Cairngorm Mountains,Scotland

Relict rock glaciers have considerable potential for contributing to palaeoclimatic reconstruction, but this potential is often undermined by lack of dating control and problems of interpretation. Here we reinvestigate and date four proposed ‘rock glaciers’ in the Cairngorm Mountains and show that the morphology of only one of these appears consistent with that of a true rock glacier produced by creep of underlying ice or ice‐rich sediment. All four features comprise rockslide or rock avalanche runout debris, and the possibility that all four represent unmodified runout accumulations cannot be discounted. Surface exposure dating of the four debris accumulations using cosmogenic ¹⁰Be produced uncertainty‐weighted mean ages of 15.4 ± 0.8 ka, 16.2 ± 1.0 ka, 12.1 ± 0.6 ka and 12.7 ± 0.8 ka. All four ages imply emplacement under cold stadial conditions, two prior to the Windermere Interstade of ca. 14.5–12.9 cal. ka BP and two during the Loch Lomond Stade of ca. 12.9–11.5 cal. ka BP. The above ages indicate that paraglacial rock‐slope failure on granite rockwalls occurred within a few millennia after deglaciation. The mean exposure ages obtained for runout debris at two sites – Strath Nethy (16.2 ± 1.0 ka) and Lairig Ghru (15.4 ± 0.8 ka) – are consistent with basal radiocarbon ages from Loch Etteridge, 22 km to the southwest (mean = 15.6 ± 0.3 cal. ka BP) and imply widespread deglaciation of the Cairngorms and adjacent valleys before 15 ka and possibly 16 ka. Copyright © 2008 John Wiley & Sons, Ltd.     

Late Quaternary landscapes in Central Australia: sedimentary history and palaeoecology of Puritjarra rock shelter

Puritjarra rock shelter provides a long record of late Quaternary vegetation in the Australian arid zone. Analysis of the sedimentary history of this rock shelter is combined with reanalysis of charcoal and phytolith records to provide a first‐order picture of changing landscapes in western Central Australia. These show a landscape responding to increasing aridity from 45 ka with deflation of clay‐rich red palaeosols (<45 ka) and sharp declines in grassland and other vegetation at 40–36 ka, and at the beginning of the Last Glacial Maximum (LGM) (24 ka). Vegetation in the catchment of the rock shelter recovered after 15 ka with expansion of both acacia woodland and spinifex grasslands, registering stronger summer rainfall in the interior of the continent. By 8.3 ka re‐vegetation of local palaeosols and dunes had choked off sediment supply to the rock shelter and the character of the sediments changed abruptly. Poaceae values peaked at 5.8 ka, suggesting the early–mid Holocene climatic optimum in Central Australia is bracketed between 8.3 and 5.8 ka. Local vegetation was disrupted in the late Holocene with a sharp decline in Poaceae at 3.8 ka, coinciding with an abrupt intensification of ENSO. Local grasslands recovered over the next two millennia and by 1.5 ka the modern vegetation appears to have become established. Copyright © 2009 John Wiley & Sons, Ltd.