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.     

Microtextures and crystal chemistry in P21/c pigeonites

Summary ?A single-crystal X-ray investigation was performed on crystals of P2₁/c natural pigeonite with varying Ca and Fe^* ( = Fe²⁺ + Mn²⁺) contents, in order to verify the effect of microtextural disorder on structure refinements and to constrain the crystal chemistry of pigeonite. Antiphase domains and exsolution lamellae affect differently the refinement results. In a crystal free of exsolution the structure obtained after refinement with all reflections is an average of that of the antiphase domains and of their boundaries, whereas in an exsolved crystal it represents only the structure of the prevailing pigeonite lamellae. The refinement using only h + k odd reflections seems to give the structure of the Ca-free pigeonite characteristic of the antiphase domains rather than that of Ca-rich domain walls. The ratio of the scale factors in refinements with all reflections and with only h + k odd reflections allows the ratios of the exsolved augite and pigeonite phases to be estimated. The crystal chemistry of the investigated samples follows the trends outlined by data on Ca-free and Fe-free synthetic samples. In particular, it is shown that Ca and Fe^* substitution for Mg induce similar changes in the average structure, i.e. both induce an expansion in the M1 polyhedron and decrease the difference between the M2–O3 distances. Received October 18, 2001; revised version accepted February 15, 2002     

Ionprobe U-Pb zircon ages from the high-pressure/low-temperature mélange of Syros, Greece: age diversity and the importance of pre-Eocene subduction

The Cycladic blueschist belt in the central Aegean Sea has experienced high‐pressure (HP) metamorphism during collisional processes between the Apulian microplate and Eurasia. The general geological and tectonometamorphic framework is well documented, but one aspect which is yet not sufficiently explored is the importance of HP mélanges which occur within volcano‐sedimentary successions. Unresolved issues concern the range in magmatic and metamorphic ages recorded by mélange blocks and the significance of eventual pre‐Eocene HP metamorphism. These aspects are here addressed in a U‐Pb zircon study focusing on the block–matrix association exposed on the island of Syros. Two gneisses from a tectonic slab of this mélange, consisting of an interlayered felsic gneiss‐glaucophanite sequence, yielded zircon ²⁰⁶Pb/²³⁸U ages of 240.1 ± 4.1 and 245.3 ± 4.9 Ma, respectively, similar to Triassic ages determined on zircon in meta‐volcanic rocks from structurally coherent sequences elsewhere in the Cyclades. This strongly suggests that parts of these successions have been incorporated in the mélanges and provides the first geochronological evidence that the provenance of mélange blocks/slabs is neither restricted to a single source nor confined to fragments of oceanic lithosphere. Zircon from a jadeitite and associated alteration zones (omphacitite, glaucophanite and chlorite‐actinolite rock) all yielded identical ²⁰⁶Pb/²³⁸U ages of c. 80 Ma. Similar Cretaceous U‐Pb zircon ages previously reported for mélange blocks have been interpreted by different authors to reflect magmatic or metamorphic ages. The present study adds a further argument in favour of the view that zircon formed newly in some rock types at c. 80 Ma, due to hydrothermal or metasomatic processes in a subduction zone environment, and supports the interpretation that the Cycladic blueschist belt records both Cretaceous and Eocene HP episodes and not only a single Tertiary HP event.     

Size and shape evolution of pores in a viscoplastic matrix under compression

The frequent use of soils and earth materials for hydraulic capping and for geo‐environmental waste containment motivated our interest in detailed modelling of changes in size and shape of macro‐pores to establish links between soil mechanical behaviour and concurrent changes in hydraulic and transport properties. The objective of this study was to use finite element analysis (FEA) to test and extend previous analytical solutions proposed by the authors describing deformation of a single macro‐pore embedded in linear viscoplastic soil material subjected to anisotropic remote stress. The FEA enables to consider more complex pore geometries and provides a detailed picture of matrix yield behaviour to explain shortcomings of approximate analytical solutions. Finite element and analytical calculations agreed very well for linear viscous as well as for viscoplastic materials, only limited for the case of isotropic remote stress due to the simplifications of the analytical model related to patterns and onset of matrix‐yielding behaviour. FEA calculations were compared with experimental data obtained from a compaction experiment in which pore deformation within a uniform modelling clay sample was monitored using CAT scanning. FEA predictions based on independently measured material properties and initial pore geometry provided an excellent match with experimentally determined evolution of pore size and shape hence lending credence to the potential use of FEA for more complex pore geometries and eventually connect macro‐pore deformation with hydraulic properties. Copyright © 2006 John Wiley & Sons, Ltd.     

Ocean-floor hydrothermal metamorphism in the Limousin ophiolites (western French Massif Central): evidence of a rare preserved Variscan oceanic marker

The Limousin ophiolite is located at the suture zone between two major thrust sheets in the western French Massif Central. This ophiolitic section comprises mantle‐harzburgite, mantle‐dunite, wehrlites, troctolites and layered gabbros. It has recorded a static metamorphic event transforming the gabbros into undeformed amphibolites and the magmatic ultramafites into serpentinites and/or pargasite‐bearing chloritites. With various thermobarometric methods, it is possible to show that the different varieties of amphibole have registered low‐P (c. 0.2 GPa) conditions with temperature ranging from high‐T, late‐magmatic conditions to greenschist–zeolite metamorphic facies. The abundance of undeformed metamorphic rocks (which is typical of the lower oceanic crust), the occurrence of Ca–Al (–Mg) metasomatism illustrated by the growth of Ca–Al silicates in veins or replacing the primary magmatic minerals, the P–T conditions of the metamorphism and the numerous similarities with oceanic crustal rocks from Ocean Drilling Program and worldwide ophiolites are the main arguments for an ocean‐floor hydrothermal metamorphism in the vicinity of a palaeo‐ridge. Among the West‐European Variscan ophiolites, the Limousin ophiolites constitute an extremely rare occurrence that has not been involved in any HP (subduction‐related) or MP (orogenic) metamorphism as observed in other ophiolite occurrences (i.e. France, Spain and Germany).     

Estimation of dispersion in estuaries by a self-consistent method

Longitudinal and transverse dispersion in estuaries may be determined in principle by measuring the change in time of the concentration field of a water soluble dye. In practice, measurements may take the form of several sampling runs crossing the dye patch. The coefficients of dispersion are usually estimated from the second moments of the concentration field; the use of second moments, in particular, introduces errors because measurements at long times or distance are unduly emphasized.A self-consistent method, free of the disadvantages inherent in the use of higher moments, is described in which the distribution of concentration in space is related to the zeroth spatial moment of concentration. In this formulation, based upon sampling runs across the dye patch, the concentration field in the estuary may be directly and simply related to the dispersion coefficients in both the intercept and gradient of a linear graph.The methods are illustrated by the characterization of a series of experiments in the Inner Estuary of Milford Haven.     

Identification of physical processes inherent in artificial neural network rainfall runoff models

The emergence of artificial neural network (ANN) technology has provided many promising results in the field of hydrology and water resources simulation. However, one of the major criticisms of ANN hydrologic models is that they do not consider/explain the underlying physical processes in a watershed, resulting in them being labelled as black‐box models. This paper discusses a research study conducted in order to examine whether or not the physical processes in a watershed are inherent in a trained ANN rainfall‐runoff model. The investigation is based on analysing definite statistical measures of strength of relationship between the disintegrated hidden neuron responses of an ANN model and its input variables, as well as various deterministic components of a conceptual rainfall‐runoff model. The approach is illustrated by presenting a case study for the Kentucky River watershed. The results suggest that the distributed structure of the ANN is able to capture certain physical behaviour of the rainfall‐runoff process. The results demonstrate that the hidden neurons in the ANN rainfall‐runoff model approximate various components of the hydrologic system, such as infiltration, base flow, and delayed and quick surface flow, etc., and represent the rising limb and different portions of the falling limb of a flow hydrograph. Copyright © 2004 John Wiley & Sons, Ltd.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

Theoretical accuracy of synthetic aperture sonar micronavigation using a displaced phase-center antenna

The Cramer-Rao lower bounds on the cross-track translation and rotation of a displaced phase-center antenna (DPCA) in the slant range plane between two successive pings (known as DPCA sway and yaw in what follows) are computed, assuming statistically homogeneous backscatter. These bounds are validated using experimental data from a 118-182-kHz sonar, showing an accuracy of the order of 20 microns on the ping-to-ping cross-track displacements. Next, the accuracy required on the DPCA sway and yaw in order to achieve a given synthetic aperture sonar (SAS) beampattern specification, specified by the expected SAS array gain, is computed as a function of the number P of pings in the SAS. Higher accuracy is required when P increases to counter the accumulation of errors during the integration of the elementary ping-to-ping estimates: the standard deviation must decrease as P/sup -1/2/ for the DPCA sway and P/sup -3/2/ for the yaw. Finally, by combining the above results, the lower bounds on DPCA micronavigation accuracy are established. These bounds set an upper limit to the SAS length achievable in practice. The maximum gain Q in cross-range resolution achievable by a DPCA micronavigated SAS is computed as a function of the key SAS parameters. These theoretical predictions are compared with simulations and experimental results.     

Metamorphism in Archaean greenstone belts: calculated fluid compositions and implications for gold mineralization

An inescapable consequence of the metamorphism of greenstone belt sequences is the release of a large volume of metamorphic fluid of low salinity with chemical characteristics controlled by the mineral assemblages involved in the devolatilization reactions. For mafic and ultramafic sequences, the composition of fluids released at upper greenschist to lower amphibolite facies conditions for the necessary relatively hot geotherm corresponds to those inferred for greenstone gold deposits (X_CO2= 0.2–0.3). This result follows from the calculation of mineral equilibria in the model system CaO–MgO–FeO–Al₂O₃–SiO₂–H₂O–CO₂, using a new, expanded, internally consistent dataset. Greenstone metamorphism cannot have involved much crustal over-thickening, because very shallow levels of greenstone belts are preserved. Such orogeny can be accounted for if compressive deformation of the crust is accompanied by thinning of the mantle lithosphere. In this case, the observed metamorphism, which was contemporaneous with deformation, is of the low-P high-T type. For this type of metamorphism, the metamorphic peak should have occurred earlier at deeper levels in the crust; i.e. the piezothermal array should be of the ‘deeper-earlier’type. However, at shallow crustal levels, the piezothermal array is likely to have been of ‘deeper-later’type, as a consequence of erosion. Thus, while the lower crust reached maximum temperatures, and partially melted to produce the observed granites, mid-crustal levels were releasing fluids prograde into shallow crustal levels that were already retrograde. We propose that these fluids are responsible for the gold mineralization. Thus, the contemporaneity of igneous activity and gold mineralization is a natural consequence of the thermal evolution, and does not mean that the mineralization has to be a consequence of igneous processes. Upward migration of metamorphic fluid, via appropriate structurally controlled pathways, will bring the fluid into contact with mineral assemblages that have equilibrated with a fluid with significantly lower X_CO2. These assemblages are therefore grossly out of equilibrium with the fluid. In the case of infiltrated metabasic rocks, intense carbonation and sulphidation is predicted. If, as seems reasonable, gold is mobilized by the fluid generated by devolatilization, then the combination of processes proposed, most of which are an inevitable consequence of the metamorphism, leads to the formation of greenstone gold deposits predominantly from metamorphic fluids.