High-T, low-P metamorphism in the Palaeoproterozoic Halls Creek Orogen, northern Australia: the middle crustal response to a mantle-related transient thermal pulse

High‐T, low‐P metamorphic rocks of the Palaeoproterozoic central Halls Creek Orogen in northern Australia are characterised by low radiogenic heat production, high upper crustal thermal gradients (locally exceeding 40 °C km^?1) sustained for over 30 Myr, and a large number of layered mafic‐ultramafic intrusions with mantle‐related geochemical signatures. In order to account for this combination of geological and thermal characteristics, we model the middle crustal response to a transient mantle‐related heat pulse resulting from a temporary reduction in the thickness of the mantle lithosphere. This mechanism has the potential to raise mid‐crustal temperatures by 150–400 °C within 10–20 Myr following initiation of the mantle temperature anomaly, via conductive dissipation through the crust. The magnitude and timing of maximum temperatures attained depend strongly on the proximity, duration and lateral extent of the thermal anomaly in the mantle lithosphere, and decrease sharply in response to anomalies that are seated deeper than 50–60 km, maintained for <5 Myr in duration and/or have half‐widths <100 km. Maximum temperatures are also intimately linked to the thermal properties of the model crust, primarily due to their influence on the steady‐state (background) thermal gradient. The amplitudes of temperature increases in the crust are principally a function of depth, and are broadly independent of crustal thermal parameters. Mid‐crustal felsic and mafic plutonism is a predictable consequence of perturbed thermal regimes in the mantle and the lowermost crust, and the advection of voluminous magmas has the potential to raise temperatures in the middle crust very quickly. Although pluton‐related thermal signatures significantly dissipate within <10 Myr (even for very large, high‐temperature intrusive bodies), the interaction of pluton‐ and mantle‐related thermal effects has the potential to maintain host rock temperatures in excess of 400–450 °C for up to 30 Myr in some parts of the mid‐crust. The numerical models presented here support the notion that transient mantle‐related heat sources have the capacity to contribute significantly to the thermal budget of metamorphism in high‐T, low‐P metamorphic belts, especially in those characterised by low surface heat flow, very high peak metamorphic geothermal gradients and abundant mafic intrusions.     

The geological and geodynamic evolution of the eastern Black Sea basin: insights from 2-D and 3-D tectonic modelling

Subsidence mechanisms that may have controlled the evolution of the eastern Black Sea have been studied and simulated using a numerical model that integrates structural, thermal, isostatic and surface processes in both two- (2-D) and three-dimensions (3-D). The model enables the forward modelling of extensional basin evolution followed by deformation due to subsequent extensional and compressional events. Seismic data show that the eastern Black Sea has evolved via a sequence of interrelated tectonic events that began with early Tertiary rifting followed by several phases of compression, mainly confined to the edges of the basin. A large magnitude (approximately 12 km) of regional subsidence also occurred in the central basin throughout the Tertiary. Models that simulate the magnitude of observed fault controlled extension (β=1.13) do not reproduce the total depth of the basin. Similarly, the modelling of compressional deformation around the edges of the basin does little to enhance subsidence in the central basin. A modelling approach that quantifies lithosphere extension according to the amount of observed crustal thinning and thickening across the basin provides the closest match to overall subsidence. The modelling also shows that deep crustal and mantle–lithosphere processes can significantly influence the rate and magnitude of syn- to post-rift subsidence and shows that such mechanisms may have played an important role in forming the anomalously thin syn-rift and thick Miocene–Quaternary sequences observed in the basin. It is also suggested that extension of a 40–45 km thick pre-rift crust is required to generate the observed magnitude of total subsidence when considering a realistic bathymetry.     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

Cation distribution and structure modelling of spinel solid solutions

 This paper presents an improved generalisation of cation distribution determination based on an accurate fit of all crystal-chemical parameters. Cations are assigned to the tetrahedral and octahedral sites of the structure according to their scattering power and a set of bond distances optimised for spinel structure. A database of 295 spinels was prepared from the literature and unpublished data. Selected compositions include the following cations: Mg²⁺, Al³⁺, Si⁴⁺, Ti⁴⁺, V³⁺, Cr³⁺, Mn²⁺, Mn³⁺, Fe²⁺, Fe³⁺, Co²⁺, Ni²⁺, Zn²⁺ and vacancies. Bond distance optimisation reveals a definite lengthening in tetrahedral distance when large amounts of Fe³⁺ or Ni²⁺ are present in the octahedral site. This means that these cations modify the octahedral angle and hence the shared octahedral edge, causing an increase in the tetrahedral distance with respect to the size of the cations entering it. Some applications to published data are discussed, showing the capacity and limitations of the method for calculating cation distribution, and for identifying inconsistencies and inaccuracies in experimental data. Received: 19 February 2001 / Accepted: 1 June 2001     

Post-mid-Cretaceous eastern North Sea evolution inferred from 3D thermo-mechanical modelling

The Late Cretaceous–Cenozoic evolution of the eastern North Sea region is investigated by 3D thermo-mechanical modelling. The model quantifies the integrated effects on basin evolution of large-scale lithospheric processes, rheology, strength heterogeneities, tectonics, eustasy, sedimentation and erosion.

The evolution of the area is influenced by a number of factors: (1) thermal subsidence centred in the central North Sea providing accommodation space for thick sediment deposits; (2) 250-m eustatic fall from the Late Cretaceous to present, which causes exhumation of the North Sea Basin margins; (3) varying sediment supply; (4) isostatic adjustments following erosion and sedimentation; (5) Late Cretaceous–early Cenozoic Alpine compressional phases causing tectonic inversion of the Sorgenfrei–Tornquist Zone (STZ) and other weak zones.

The stress field and the lateral variations in lithospheric strength control lithospheric deformation under compression. The lithosphere is relatively weak in areas where Moho is deep and the upper mantle warm and weak. In these areas the lithosphere is thickened during compression producing surface uplift and erosion (e.g., at the Ringkøbing–Fyn High and in the southern part of Sweden). Observed late Cretaceous–early Cenozoic shallow water depths at the Ringkøbing–Fyn High as well as Cenozoic surface uplift in southern Sweden (the South Swedish Dome (SSD)) are explained by this mechanism.

The STZ is a prominent crustal structural weakness zone. Under compression, this zone is inverted and its surface uplifted and eroded. Contemporaneously, marginal depositional troughs develop. Post-compressional relaxation causes a regional uplift of this zone.

The model predicts sediment distributions and paleo-water depths in accordance with observations. Sediment truncation and exhumation at the North Sea Basin margins are explained by fall in global sea level, isostatic adjustments to exhumation, and uplift of the inverted STZ. This underlines the importance of the mechanisms dealt with in this paper for the evolution of intra-cratonic sedimentary basins.     

Lithosphere structure of European sedimentary basins from regional three-dimensional gravity modelling

A three-dimensional (3D) density model, approximated by two regional layers—the sedimentary cover and the crystalline crust (offshore, a sea-water layer was added), has been constructed in 1° averaging for the whole European continent. The crustal model is based on simplified velocity model represented by structure maps for main seismic horizons—the “seismic” basement and the Moho boundary. Laterally varying average density is assumed inside the model layers. Residual gravity anomalies, obtained by subtraction of the crustal gravity effect from the observed field, characterize the density heterogeneities in the upper mantle. Mantle anomalies are shown to correlate with the upper mantle velocity inhomogeneities revealed from seismic tomography data and geothermal data. Considering the type of mantle anomaly, specific features of the evolution and type of isostatic compensation, the sedimentary basins in Europe may be related into some groups: deep sedimentary basins located in the East European Platform and its northern and eastern margins (Peri-Caspian, Dnieper–Donets, Barents Sea Basins, Fore–Ural Trough) with no significant mantle anomalies; basins located on the activated thin crust of Variscan Western Europe and Mediterranean area with negative mantle anomalies of −150 to −200×10⁻⁵ ms⁻² amplitude and the basins associated with suture zones at the western and southern margins of the East European Platform (Polish Trough, South Caspian Basin) characterized by positive mantle anomalies of 50–150×10⁻⁵ ms⁻² magnitude. An analysis of the main features of the lithosphere structure of the basins in Europe and type of the compensation has been carried out.     

Propriétés fractales de la fragmentation et processus stochastiques de fracturation : approche géométrique 3D à l'aide du modèle ObsifracFractal properties of fragmentation and stochastic fracturing processes: geometrical 3D approach by OBSIFRAC model

A numerical rock fragmentation model was elaborated, producing a 3D puzzle of convex polyhedra, geometrically described in a database. In the first scenario, a constant proportion of blocks are fragmented at each step of the process and leads to fractal distribution. In the second scenario, division affects one random block at each stage of the process, and produces a Weibull volume distribution law. Imposing a minimal distance between the fractures, the third scenario reveals a power law. The inhibition of new fractures in the neighbourhood of existing discontinuities could be responsible for fractal properties in rock mass fragmentation. To cite this article: L. Empereur-Mot, T. Villemin, C. R. Geoscience 334 (2002) 127–133.     

Spatial variation of seismic motion induced by propagation of body waves

Spatial variation of earthquake ground motion is an important phenomenon that cannot be ignored in the design and safety of strategic structures. Several models have been developed to describe this variation using statistical, mathematical or physical approaches. The latter approach is not specific to an event. A recent contribution, which uses such an approach and called complete stochastic deamplification approach (CSDA), was developed [1]. The aim of this paper is to analyze the spatial variation of earthquake motion induced by the propagation of body waves using the CSDA. Coherency functions are evaluated for the cases of SH–SV–P waves propagating through stratified soil. Results obtained show that the variation of the coherency function is not the same for vertical and horizontal components and that the motion is more coherent at depth than at the free surface. In fact, we found that the rate of decrease with frequency and distance is not the same if P–SV waves propagate through stratified soil.     

Optimal placement of dampers and actuators based on stochastic approach

A general method is developed for optimal application of dampers and actuators by installing them at optimal location on seismic-resistant structures. The study includes development of a statistical criterion, formulation of a general optimization problem and establishment of a solution procedure. Numerical analysis of the seismic response in time-history of controlled structures is used to verify the proposed method for optimal device application and to demonstrate the effectiveness of seismic response control with optimal device location. This study shows that the proposed method for the optimal device application is simple and general, and that the optimally applied dampers and actuators are very efficient for seismic response reduction.     

河流水质风险评价的灰色-随机风险率方法

提出了量化影响河流水质的随机不确定性与灰色不确定性的水质超标灰色－随机风险率概念，建立了水质超标灰色－随机风险率评价模型。在水质单项参数评价模型中，将河流污染物浓度变量的分布处理成灰色概率分布，将污染物浓度超过水质类别标准值的风险率处理成灰色概率，即水质超标灰色－随机风险率。在水质综合评价模型中，河流水环境系统被考虑为担任某一使用可能的可靠性系统，而任意一种水质参数超标意味着河流水体使用功能不能得到应有的保证，也即表明水体综合评价超标，最后借鉴系统可靠性分析的理论和方法计算水质综合超标率。该方法应用于黄河花园口断面重金属污染风险评价。