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.     

Hydrogeochemistry of the Simav geothermal field, western Anatolia, Turkey

Thermal waters hosted by Menderes metamorphic rocks emerge along fault lineaments in the Simav geothermal area. Thermal springs and drilled wells are located in the Eynal, Çitgöl and Na a locations, which are part of the Simav geothermal field. Studies were carried out to obtain the main chemical and physical characteristics of thermal waters. These waters are used for heating of residences and greenhouses and for balneological purposes. Bottom temperatures of the drilled wells reach 163°C with total dissolved solids around 2225 mg/kg. Surface temperatures of thermal springs vary between 51°C and 90°C. All the thermal waters belong to Na–HCO₃–SO₄ facies. The cold groundwaters are Ca–Mg–HCO₃ type. Dissolution of host rock and ion-exchange reactions in the reservoir of the geothermal system shift the Ca–Mg–HCO₃ type cold groundwaters to the Na–HCO₃–SO₄ type thermal waters. Thermal waters are oversaturated at discharge temperatures for aragonite, calcite, quartz, chalcedony, magnesite and dolomite minerals giving rise to a carbonate-rich scale. Gypsum and anhydrite minerals are undersaturated with all of the thermal waters. Boiling during ascent of the thermal fluids produces steam and liquid waters resulting in an increase of the concentrations of the constituents in discharge waters. Steam fraction, y, of the thermal waters of which temperatures are above 100°C is between 0.075 and 0.119. Reservoir pH is much lower than pH measured in the liquid phase separated at atmospheric conditions, since the latter experienced heavy loss of acid gases, mainly CO₂. Assessment of the various empirical chemical geothermometers and geochemical modelling suggest that reservoir temperatures vary between 175°C and 200°C.     

华北板内深部构造

华北板块的形成经历了早前寒武纪、燕山期及喜马拉雅期3个主要构造发展期，由于华北板块自身运动及所受应力场的作用，加之上地幔岩石圈的不均一性等因素，在中、新生代形成许多特殊的板内构造块。综合应用地质、地球物理和地球化学的成果，对华北板内深部结构进行了研究。从深部构造角度划分出6个金及金金属成矿带、4个金刚石成矿带，并对华北地区的地震及地热资源与新生事大陆裂谷的关系进行了探讨。     

模拟酸雨对福建四种山地土壤的淋溶与风化作用

在不同pH值 (2 .0、3.5、5 .0、6 .0及蒸馏水 )的模拟酸雨作用下 ,对福建四种类型山地土壤进行了为期 12 0d的淋溶试验。结果表明 ,盐基离子的淋溶量随模拟酸雨pH值的降低而增加 ,尤其当pH≤ 3.5时 ,增加最明显。 4种类型土壤中各离子淋失量大小均为 :Ca2 >Mg2 >K >Na ,这与原土壤中交换性盐基含量一致。在pH≥ 3.5的酸雨作用下 ,土壤中铝的淋失较少 ,而当酸雨pH =2 .0时 ,铝的淋失骤然升高。pH2 .0的模拟酸雨不仅增加了土壤铝离子的释放 ,而且加速了矿物的风化。在pH2 .0的模拟酸雨淋溶下 ,红壤、黄壤、山地草甸土和紫色土中由矿物风化输入的 5种阳离子总量分别为 2 .94、4 .33、4 .35和 6 .0 4 (cmol kg)土 ,显著大于其它酸雨处理。土壤对酸沉降的缓冲作用由阳离子交换、氢氧化铝水解以及原生矿物风化缓冲作用三部分组成 ,其相对重要性不仅决定于土壤酸度、盐基饱和度和矿物组成 ,而且与酸雨的pH值密切相关