Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Mezen Bay—a key area for understanding Weichselian glaciations in northern Russia

Sediment successions in coastal cliffs around Mezen Bay, southeastern White Sea, record an unusually detailed history of former glaciations, interstadial marine and fluvial events from the Weichselian. A regional glaciation model for the Weichselian is based on new data from the Mezen Bay area and previously published data from adjacent areas. Following the Mikulinian (Eemian) interglacial a shelf‐centred glaciation in the Kara Sea is reflected in proglacial conditions at 100–90 ka. A local ice‐cap over the Timan ridge existed between 75 and 65 ka. Renewed glaciation in the Kara Sea spread southwestwards around 60 ka only, interrupted by a marine inundation, before it advanced to its maximum position at about 55–50 ka. After a prolonged ice‐free period, the Scandinavian ice‐sheet invaded the area from the west and terminated east of Mezen Bay about 17 ka. The previously published evidence of a large ice‐dammed lake in the central Arkhangelsk region, Lake Komi, finds no support in this study. Copyright © 2003 John Wiley & Sons, Ltd.     

Glacial rebound and sea-level change in the British Isles

Observations of sea levels around the coastline of the British Isles for the past 10,000–15,000 years exhibit a major regional variation and provide an important data base for testing models of glacial rebound as well as models of the Late Devensian ice sheet. A high-resolution rebound model has been developed which is consistent with both the spatial and temporal patterns of sea-level change and which demonstrates that the observations are the result of (i) the glacio-isostatic crustal rebound in response to the unloading of the ice sheet over Britain and, to a lesser degree, of the ice sheet over Fennoscandia, and (ii) the rise in sea-level from the melting Late Pleistocene ice sheets, including the response of the crust to the water loading (the hydro-isostatic effect). The agreement between model and observations is such that there is no need to invoke vertical crustal movements for Great Britain and Ireland of other than glacio-hydro-isostatic origin. The rebound contributions are important throughout the region and nowhere is it sufficiently small for the sea-level change to approximate the eustatic sea-level rise. The observational data distribution around the periphery as well as from sites near the centre of the former ice sheet is sufficient to permit constraints to be established on both earth model parameters specifying the mantle viscosity and lithospheric thickness and the extent and volume of the ice sheet at the time of the last glaciation. Preliminary solutions are presented which indicate an upper mantle viscosity of (3–5)10²⁰ Pas, a lithospheric thickness of about 100 km or less, and an ice model that was not confluent with the Scandinavian ice sheet during the last glaciation and whose maximum thickness over Scotland is unlikely to have exceeded about 1500 m.     

Ultra-Deep Optical Spectroscopy with PMAS. Using the Nod-and-Shuffle Technique

PMAS, the Potsdam Multi-Aperture Spectrophotometer, is a new integral field spectrograph in the optical, which is optimized for good transmission and high image quality from 350 nm to 1 μm. We present our plan to implement a CCD charge-shuffle mode to allow for beam switching with a very high degree of sky subtraction accuracy for faint object 3-D spectroscopy. This revised version was published online in July 2006 with corrections to the Cover Date.     

Electric Discharges in the Protoplanetary Nebula as a Source of Impulse Magnetic Fields to Promote Dust Aggregation

Lightning discharge generated in the protoplanetary nebula is viewed as a temporally isolated surge in the flow of electrically charged particles, similar to that of terrestrial lightning. If the current is intense enough, a powerful circular impulse magnetic field is generated around the instantaneous virtual electric conductor. Such magnetic field is capable of magnetizing dust grains containing ferromagnetic components present in its vicinity to their saturation levels. As a result, dust grains attract one another, forming the aggregates. This magnetically driven attraction suggests an important process possibly operational at an early stage of the planetary accretion. Based on both a classical model for electric conductor, and the theory of Lienard–Wiechert electromagnetic potentials, our calculations show that the magnetic impulse due to a discharge channel of a few cm in diameter transferring a charge of about 10⁴ electrons reaches as high as 10 T. At these magnetic fields, the ferromagnetic dust grains, and possibly the already-formed larger aggregates as well, are easily magnetized to the saturation levels, producing compact clusters exhibiting permanent magnetic moments.     

Invariance and scaling properties in the distributions of contributing area and energy in drainage basins

The cumulative probability distributions for stream order, stream length, contributing area, and energy dissipation per unit length of channel are derived, for an ordered drainage system, from Horton's laws of network composition. It is shown how these distributions can be related to the fractal nature of single rivers and river networks. Finally, it is shown that the structure proposed here for these probability distributions is able to fit the observed frequency distributions, and their deviations from straight lines in a log-log plot.     

Fluid entrapment and reequilibration during subduction and exhumation: A case study from the high-grade Nestos shear zone, Central Rhodope, Greece

A fluid inclusion study on metamorphic minerals of successive growth stages was performed on highly deformed paragneisses from the Nestos Shear Zone at Xanthi (Central Rhodope), in which microdiamonds provide unequivocal evidence for ultrahigh-pressure (UHP) metamorphism. The correlation of fluid inclusion density isochores and fluid inclusion reequilibration textures with geothermobarometric data and the relative chronology of micro- and macro-scale deformation stages allow a better understanding of both the fluid and metamorphic evolution along the P–T–d path. Textural evidence for subduction towards the NE is recorded by the orientation of intragranular NE-oriented fluid inclusion planes and the presence of single, annular fluid inclusion decrepitation textures. These textures occur within quartz “foam” structures enclosed in an earlier generation of garnets with prolate geometries and rarely within recrystallized matrix quartz, and reequilibrated both in composition and density during later stages of exhumation. No fluid inclusions pertaining to the postulated ultrahigh-pressure stage for microdiamond-bearing garnet–kyanite–gneisses have yet been found. The prolate shape of garnets developed during the earliest stages of exhumation that is recorded structurally by (L  S) tectonites, which subsequently accommodated progressive ductile SW shearing and folding up to shallow crustal levels. The majority of matrix kyanite and a later generation of garnet were formed during SW-directed shear under plane-strain conditions. Fluid inclusions entrapped in quartz during this stage of deformation underwent density loss and transformed to almost pure CO₂ inclusions by preferential loss of H₂O. Those inclusions armoured within garnet retained their primary 3-phase H₂O–CO₂ compositions. Reequilibration of fluid inclusions in quartz aggregates is most likely the result of recrystallization along with stress-induced, preferential H₂O leakage along dislocations and planar lattice defects which results in the predominance of CO₂ inclusions with supercritical densities. Carbonic fluid inclusions from adjacent kyanite–corundum-bearing pegmatoids and, the presence of shear-plane-parallel fluid inclusion planes within late quartz boudin structures consisting of pure CO₂-fluid inclusions with negative crystal shapes, bear witness of the latest stage of deformation by NE-directed extensional shear.This study shows that the textures of early fluid inclusions that formed already during the prograde metamorphic path can be preserved and used to derive information about the kinematics of subduction that is difficult to obtain from other sources. The textures of early inclusions, together with later generations of unaltered primary and secondary inclusions in metamorphic index minerals that can be linked to specific deformation stages and even P–T conditions, are a welcome supplement for the reconstruction of a rather detailed P–T–d path.     

Morphological analysis of the drainage system in the Eastern Alps

We study the morphology of the major rivers draining the Eastern Alps to test whether the active tectonics of this part of the orogen is reflected in the shape of channel profiles of the river network. In our approach we compare channel profiles measured from digital elevation models with numerically modelled channel profiles using a stream power approach. It is shown that regions of high stream power coincide largely with regions of highest topography and largest uplift rates, while the forelands and the Pannonian Basin are characterised by a significantly lower stream power. From stream power modelling we conclude that there is young uplift at the very east of the Eastern Alps, in the Bohemian Massif and in the Pohorje Range. The impact of the Pleistocene glaciations is explored by comparing properties of rivers that drain in proximal and distal positions relative to the ice sheet during the last glacial maximum. Our analysis shows that most knick points, wind gaps and other non-equilibrium features of catchments covered by ice during the last glaciations (Salzach, Enns) can be correlated with glacial processes. In contrast the ice free catchments of the Mur and Drava are characterized by channels in morphological equilibrium at the first approximation and are showing only weak evidence of the strong tectonic activity within these catchments. Finally, the channel profiles of the Adige and the divide between the upper Rhine and Danube catchments differ significantly from the other catchments. We relate this to the fact that the Adige and the Rhine respond to different base levels from the remainder of the Eastern Alps: The Adige may preserve a record from the Messininan base level change and the Rhine is subject to the base level lowering in the Rhine Graben.     

Use of sinkhole and specific capacity distributions to assess vertical gradients in a karst aquifer

The carbonate-rock aquifer in the Great Valley, West Virginia, USA, was evaluated using a database of 687 sinkholes and 350 specific capacity tests to assess structural, lithologic, and topographic influences on the groundwater flow system. The enhanced permeability of the aquifer is characterized in part by the many sinkholes, springs, and solutionally enlarged fractures throughout the valley. Yet, vertical components of subsurface flow in this highly heterogeneous aquifer are currently not well understood. To address this problem, this study examines the apparent relation between geologic features of the aquifer and two spatial indices of enhanced permeability attributed to aquifer karstification: (1) the distribution of sinkholes and (2) the occurrence of wells with relatively high specific capacity. Statistical results indicate that sinkholes (funnel and collapse) occur primarily along cleavage and bedding planes parallel to subparallel to strike where lateral or downward vertical gradients are highest. Conversely, high specific capacity values are common along prominent joints perpendicular or oblique to strike. The similarity of the latter distribution to that of springs suggests these fractures are areas of upward-convergent flow. These differences between sinkhole and high specific capacity distributions suggest vertical flow components are primarily controlled by the orientation of geologic structure and associated subsurface fracturing.