A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. 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     

A qualitative and quantitative evaluation of experimental model catchment evolution

Due to the geological time scales required for observation of catchment evolution, surrogates or analogues of field data are necessary to understand long‐term processes. To investigate long‐term catchment behaviour, two experimental model catchments that developed without rigid boundaries under controlled conditions are examined and a qualitative and quantitative analysis of their evolution is presented. Qualitatively, the experimental catchments have the visual appearance of field scale data. Observation demonstrates that changes in catchment shape and network form are conservative. Quantitative analysis suggests that the catchments reach an equilibrium form while a reduction in the channel network occurs. While the catchments are laboratory scale models, the results provide insights into field scale behaviour. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Multistep matched-field inversion for broad-band data from ASIAEX2001

This paper discusses the results of geoacoustic inversion carried out using explosive charge data from the Asian Seas International Acoustic Experiment (ASIAEX) East China Sea (ECS) Experiment. A multifrequency incoherent matched-field inversion processor and a genetic algorithm (GA) are used for the inversion. A multistep matched field inversion approach is presented, which makes use of the varying sensitivities of wave fields at various frequencies to reduce the inversion problem into a sequence of smaller inversions with fewer unknowns to estimate at each stage. Different parameters are estimated using data at different frequencies according to their sensitivities. Inversion results for different areas in the ECS region are summarized and compared with core data.     

Structural development along the Billefjorden Fault Zone in the area between Kjellströmdalen and Adventdalen/Sassendalen, central Spitsbergen

The Billefjorden Fault Zone represents a major lineament on Spitsbergen with a history of tectonic activity going back into the Devonian and possibly earlier. Recent structural, sedimcntological and stratigraphical investigations indicate that most of the stratigraphic thickness variations within the Mesozoic strata along the Billefjorden Fault Zone south of Isfjordcn are due to Tertiary compressional tectonics related to the transpressive Eocene West-Spitsbergen Orogeny. No convincing evidence of distinct Mesozoic extensional events, as suggested by previous workers, has been recognized. Tertiary compressional tectonics are characterized by a combined thin-skinned/thick-skinned structural style. Decollement zones arc recognized in the Triassic Sassendalen Group (tower Décollement Zone) and in the Jurassic/Cretaceous Janusfjellet Subgroup (Upper Décollement Zone). East-vergent folding and reverse faulting associated with these decollement' zones have resulted in the development of compressional structures, of which the major arc the Skolten and Tronfjellct Anticlines and the Advcntelva Duplex. Movements on one or more high angle east-dipping reverse faults in the pre-Mesozoic basement have resulted in the development of the Juvdalskampcn Monocline, and are responsible for out-of-sequence thrusting and thinning of the Mesozoic sequence across the Billefjorden Fault Zone. Preliminary shortening calculations indicate an eastward displacement of minimum 3-4 km, possibly as much as 10 km for the Lower Cretaceous and younger rocks across the Billefjorden Fault Zone.     

Deep-ocean piezometer probe technology for geotechnical investigations

Two single-sensor piezometer probes, 8 mm in diameter, were developed for deep-ocean geotechnical investigations. These probes were tested in a hyperbaric chamber pressurized to 55 MPa (8000 psi). Testing was performed for a period of five weeks under high hydrostatic pressure with the probes inserted in reconstituted illitic marine sediment. Small differential pore-water pressures were generated in response to both mechanically and thermally generated forcing functions. During deep-ocean simulated pressure tests, the sensors exhibited excellent sensitivity and stability. These developments in piezometer-probe technology provide a quantitative means of assessing important geotechnical parameters of fine-grained seabed deposits.     

A one-dimensional model for the interaction between continental-scale ice sheets and atmospheric stationary waves

 The great continental ice sheets of the Pleistocene represented a significant topographic obstacle to the westerly winds in northern midlatitudes. This work explores how consequent changes in the atmospheric stationary wave pattern might have affected the shape and growth of the ice sheets themselves. A one dimensional (1-D) model is developed which permits an examination of the types and magnitudes of the feedbacks that might be expected. When plausible temperature perturbations are introduced at the ice-sheet margin which are proportional to the stationary wave amplitude, the equilibrium shape of the ice sheet is significantly altered, and depends on the sign of the perturbation. The proposed feedback also affects the response of the ice sheet to time-varying climate forcing. The results suggest that the evolution of a continental-scale ice sheet with a land-based margin may be significantly determined by the changes it induces in the atmospheric circulation. Received: 1 October 1999 / Accepted: 17 July 2000     

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.