New insights into the genesis of volcanic-hosted massive sulfide deposits on the seafloor from numerical modeling studies

Numerical computer simulations have been used to gain insight into the evolution of marine hydrothermal systems and the formation conditions of massive sulfide deposits in ancient and modern submarine volcanic terrains. Simulation results have been used to gain a better understanding of the formation of massive sulfide ore deposits, their location, zonation, size, and occurrence in various geotectonic settings.Most hydrothermal fluid discharging at the seafloor exhibits temperatures ranging from 200 °C to about 410 °C and average fluid discharge velocities of 1 to 2 m/s in agreement with seafloor observations. Mass calculations imply that average massive sulfide deposits may form in ~ 5000 years while giant deposits take longer than 5000 years to accumulate; supergiant deposits either need much longer time to form (> 35,000 years) or at least 100 ppm of metal in solution. Results indicate that supergiant deposits may only form in certain geotectonic environments where longevity and preservation potential of the hydrothermal system are high. An additional process (mineral precipitate cap) is proposed here to explain the zinc content of massive sulfide deposits. This cap would prevent the widespread dissolution of anhydrite and the ‘wash-out’ of zinc by subsequent hydrothermal fluid discharge.     

Advances in our understanding of the gem corundum deposits of the West Pacific continental margins intraplate basaltic fields

Recent discoveries over the last decade of new gemfields, exploitation of new and existing deposits, and application of relatively new techniques have greatly increased our knowledge of the basalt-derived gem sapphire–ruby–zircon deposits. In this paper we focus on the Late Mesozoic to Cenozoic intraplate basaltic fields of the West Pacific continental margins. We review advances made in understanding the genesis of these deposits, based on the application of newer techniques. We also critically review existing data on the gem corundum deposits, in order to further refine a model for their origin.In some of the intraplate basaltic fields, corundum-bearing xenoliths have been found showing a range of PT formation conditions from 790 °C at 0.85 GPa to as much as 1100 to 1200 °C at 1.0 to 2.5 GPa. Although most magmatic sapphires contain syngenetic inclusions of columbite-group phases, zircon, spinel and rutile, some contain additional nepheline and K-feldspar, suggesting crystallization from more undersaturated alkaline magma while the Weldborough field of NE Tasmania also contains molybdenite and beryl, suggesting at least some interaction with more fractionated ‘granitic-type’ magmas. There is a large range in PT conditions calculated for the metamorphic rubies (from 780 to 940 °C, through 800 to 1150 °C up to 1000 to 1300 °C). Fluid/melt inclusion studies on magmatic corundums generally suggest that they formed in a CO₂-rich environment from a ‘syenitic’ melt under a range of T conditions from 720 to 880 °C up to 1000 to 1200 °C. Oxygen isotope studies reveal that typical magmatic corundums have values of + 4.4 to 6.9‰, whereas metamorphic corundums from the same basaltic host have lower values of + 1.3 to 4.2‰.Geochronological studies have shown that some fields produced a simple eruptive and zircon/corundum crystallization event while others had multiple eruptive events but only one or two zircon crystallization events. For a few fields, some corundums/zircons crystallized in storage regions and then remained relatively inert for periods of 200 to 400 Ma without significant change before transport to the surface in the Cenozoic. Tectonic studies of the Australian region suggest that many of the corundums crystallized from magmas related to episodic basaltic volcanism in a tectonic regime of extension, associated with the opening of the Tasman and Coral Seas. For the Asian region, the magmatic–polygenetic corundums within the basaltic fields largely crystallized in a tectonic regime of distributed E–W extension, whereas the metamorphic-metasomatic corundums crystallised in a transpressional regime associated with the collision of the Indian Plate with the Eurasian Plate (e.g., [Garnier, V., Giuliani, G., Maluski, H., Ohnenstetter, D., Deloule, E., 2003. Ar–Ar and U–Pb ages of marble-hosted ruby deposits from Central and South-east Asia. Geophysical Research Abstracts 5, 03751; Garnier, V., Giuliani, G., Ohnenstetter, D., and Schwarz, D., 2004. Les gisements de corindon: classification et genese. Les placers a corindon gemme. Le Regne Mineral 55, 7-47; Garnier, V., Ohnenstetter, D., Giuliani, G., Maluski, H., Deloule, E., Phan Trong, T., Pham Van, L., Hoang Quang, V., 2005a. Age and significance of ruby-bearing marble from the Red River Shear Zone, Northern Vietnam. Canadian Mineralogist 43, 1315–1329]).     

Shear wave velocity and attenuation structure for the shallow crust of the southern Korean peninsula from short period Rayleigh waves

We analyzed the short period Rayleigh waves from the first crustal-scale seismic refraction experiment in the Korean peninsula, KCRUST2002, to determine the shear wave velocity and attenuation structure of the uppermost 1 km of the crust in different tectonic zones of the Korean peninsula and to examine if this can be related to the surface geology of the study area. The experiment was conducted with two large explosive sources along a 300-km long profile in 2002. The seismic traces, recorded on 170 vertical-component, 2-Hz portable seismometers, show distinct Rayleigh waves in the period range between 0.2 s and 1.2 s, which are easily recognizable up to 30–60 km from the sources. The seismic profiles, which traverse three tectonic regions (Gyeonggi massif, Okcheon fold belt and Yeongnam massif), were divided into five subsections based on tectonic boundaries as well as lithology. Group and phase velocities for the five subsections obtained by a continuous wavelet transform method and a slant stack method, respectively, were inverted for the shear wave models. We obtained shear wave velocity models up to a depth of 1.0 km. Overall, the shear wave velocity of the Okcheon fold belt is lower than that of the Gyeonggi and Yeongnam massifs by  0.4 km/s in the shallowmost 0.2 km and by 0.2 km/s at depths below 0.2 km. Attenuation coefficients, determined from the decay of the fundamental mode Rayleigh waves, were used to obtain the shear wave attenuation structures for three subsections (one for each of the three different tectonic regions). We obtained an average value of Q_β^− 1 in the upper 0.5 km for each subsection. Q_β^− 1 for the Okcheon fold belt ( 0.026) is approximately three times larger than Q_β^− 1 for the massif areas ( 0.008). The low shear wave velocity in the Okcheon fold belt is consistent with the high attenuation in this region.     

Palaeomagnetism of the Loch Doon Granite Complex, Southern Uplands of Scotland: The Late Caledonian palaeomagnetic record and an Early Devonian episode of True Polar Wander

The Southern Uplands terrane is an Ordovician–Silurian back-arc/foreland basin emplaced at the northern margin of the Iapetus Ocean and intruded by granite complexes including Loch Doon (408.3 ± 1.5 Ma) during Early Devonian times. Protracted cooling of this 130 km³ intrusion recorded magnetic remanence comprising a predominant (‘A’) magnetisation linked to initial cooling with dual polarity and mean direction D / I = 237 / 64° (α₉₅ = 4°, palaeopole at 316°E, 21°N). Subsidiary magnetisations include Mesozoic remanence correlating with extensional tectonism in the adjoining Irish Sea Basin (‘B’, D / I = 234/− 59°) and minority populations (‘C’, D / I = 106/− 2° and ‘D’, D / I = 199/1°) recording emplacement of younger ( 395 Ma) granites in adjoining terranes and the Variscan orogenic event. The ‘A’ directions have an arcuate distribution identifying anticlockwise rotation during cooling. A comparable rotation is identified in the Orthotectonic Caledonides to the north and the Paratectonic Caledonides to the south following closure of Iapetus. Continental motion from midsoutherly latitudes ( 40°S) at 408 Ma to equatorial palaeolatitudes by  395 Ma is identified and implies minimum rates of continental movement between 430 and 390 Ma of 30–70 cm/year, more than double maximum rates induced by plate forces and interpreted as a signature of true polar wander. Silurian–Devonian palaeomagnetic data from the British–Scandinavian Caledonides define a 430–385 Ma closed loop comparable to the distributed contemporaneous palaeomagnetic poles from Gondwana. They reconcile pre-430 Ma and post-380 Ma APW from this supercontinent and show that Laurentia–Baltica–Avalonia lay to the west of South America with a relict Rheic Ocean opening to the north which closed to produce Variscan orogeny by a combination of pivotal closure and right lateral transpression.     

Syntectonic crustal melting and high-grade metamorphism in a transpressional regime, Variscan Massif Central, France

Hot collisional orogens are characterized by abundant syn-kinematic granitic magmatism that profoundly affects their tectono-thermal evolutions. Voluminous granitic magmas, emplaced between 360 and 270 Ma, played a visibly important role in the evolution of the Variscan Orogen. In the Limousin region (western Massif Central, France), syntectonic granite plutons are spatially associated with major strike–slip shear zones that merge to the northwest with the South Armorican Shear Zone. This region allowed us to assess the role of magmatism in a hot transpressional orogen. Microstructural data and U/Pb zircon and monazite ages from a mylonitic leucogranite indicate synkinematic emplacement in a dextral transpressional shear zone at 313 ± 4 Ma. Leucogranites are coeval with cordierite-bearing migmatitic gneisses and vertical lenses of leucosome in strike–slip shear zones. We interpret U/Pb monazite ages of 315 ± 4 Ma for the gneisses and 316 ± 2 Ma for the leucosomes as the minimum age of high-grade metamorphism and migmatization respectively. These data suggest a spatial and temporal relationship between transpression, crustal melting, rapid exhumation and magma ascent, and cooling of high-grade metamorphic rocks.Some granites emplaced in the strike–slip shear zone are bounded at their roof by low dip normal faults that strike N–S, perpendicular to the E–W trend of the belt. The abundant crustal magmatism provided a low-viscosity zone that enhanced Variscan orogenic collapse during continued transpression, inducing the development of normal faults in the transpression zone and thrust faults at the front of the collapsed orogen.     

Deep-rooted piercement structures in deep sedimentary basins — Manifestations of supercritical water generation at depth?

Deep-rooted enigmatic piercement structures in sedimentary basins, including ‘mud volcanoes’, ‘shale diapirs’, ‘salt diapirs’, and ‘asphalt volcanoes’, range in size from less than 1 km², surface area, up to 64 km², and have often an unknown depth of penetration due to incomplete imaging. We propose that they form a family associated with fluid flow. Our argument is based partly on their inferred location (above deep faults) and on the chemical analysis of emitted products, which includes liquid clays, brines and other substances from salt diapirs, and asphalt and light oils from the asphalt volcanoes. We explain these compositions by chemical alteration caused partly by supercritical water, a phase of water existent at high pressure and temperature, locally and temporarily achieved at depths generally beyond 10 km below surface, i.e., at the sediment–crust boundary. Our hypothesis overcomes some of the problems with interpreting fluid flow products, which are otherwise very difficult to explain. In case this hypothesis can be further verified, the family could perhaps be called ‘hydrothermally associated piercement structures’.     

Le passage Ordovicien–Silurien et la partie inférieure du Silurien (Sud-Est du Massif armoricain, France)The Ordovician–Silurian transition and the lower part of the Silurian (southeastern Armorican Massif, France)

A new section in the Silurian graptolitic ‘phtanites’ (black cherts) of Les Fresnaies at Chalonnes-sur-Loire (SE Armorican Massif) shows for the first time that these rocks (1) succeed conformably to Uppermost Ordovician (Hirnantian) glaciomarine deposits and (2) contain successive graptolite assemblages that characterise the base of the Silurian, the whole Rhuddanian and Aeronian stages and the lower part of the Telychian. To cite this article: J.M. Piçarra et al., C. R. Geoscience 334 (2002) 1177–1183.     

Sustained turbidity currents and their interaction with debrite-related topography; Labuan Island, offshore NW Borneo, Malaysia

The Temburong Fm (Early Miocene), Labuan Island, offshore NW Borneo, was deposited in a lower-slope to proximal basin-floor setting, and provides an opportunity to study the deposits of sustained turbidity currents and their interaction with debrite-related topography. Two main gravity-flow facies are identified; (i) slump-derived debris-flow deposits (debrites) — characterised by ungraded silty mudstones in beds 1.5 to > 60 m thick which are rich in large (> 5 m) lithic clasts; and (ii) turbidity current deposits (turbidites) — characterised by medium-grained sandstone in beds up to 2 m thick, which contain structureless (T_a) intervals alternating with planar-parallel (T_b) and current-ripple (T_c) laminated intervals. Laterally discontinuous, cobble-mantled scours are also locally developed within turbidite beds. Based on these characteristics, these sandstones are interpreted to have been deposited by sustained turbidity currents. The cobble-mantled scours indicate either periods of intense turbidity current waxing or individual flow events. The sustained turbidity currents are interpreted to have been derived from retrogressive collapse of sand-rich mouth bars (breaching) or directly from river effluent (hyperpycnal flow). Analysis of the stratal architecture of the two facies indicates that routing of the turbidity currents was influenced by topographic relief developed at the top of the underlying debrite. In addition, turbidite beds are locally eroded at the base of an overlying debrite, possibly due to clast-related substrate ‘ploughing’ during the latter flow event. This study highlights the difficulty in constraining the origin of sustained turbidity currents in ancient sedimentary sequences. In addition, this study documents the importance large debrites may have in generating topography on submarine slopes and influencing routing of subsequent turbidity currents and the geometry of their associated deposits.     

Potential source rocks, organic geochemistry and thermal maturation in the southern depocenter (Kipourio–Grevena) of the Mesohellenic Basin, central Greece

The study area is the southern depocenter (depth > 4200 m) of the Mesohellenic Basin which extends between Kipourio and Grevena, central Greece. The Mesohellenic Basin is a Middle-Tertiary intramontane basin developed within the Hellenide orogen. Previous studies have focused on the depositional environments, configuration and hydrocarbon potential of the basin. In this paper we present additional geochemical and petrographic data from outcrop samples of the basin's southern depocenter, which is considered the most promising area, in terms of hydrocarbon prospectivity. A total number of thirty six samples were analysed: Rock-Eval pyrolysis, maceral analysis, vitrinite reflectance and thermal alteration index, bitumens extraction, liquid chromatography, and GC-MS. The samples were collected from deltaic deposits and submarine fan sediments of Late Eocene to Late Oligocene age. The TOC values of the analysed samples range between rich and very rich and the organic matter consists mainly of type III kerogen and the organic matter consider to be predominately gas prone. The thermal maturity assessed from T_max and vitrinite reflectance shows an immature stage of the organic matter along with the presence of layers having reached the very early mature stage. Vitrinite reflectance measurements and maturity calculations (applying the Lopatin modeling), reveal that the lower part of the depocenter sediments falls within the ‘oil window’. The extractable organic matter (EOM) (mg bitumens/g TOC) indicate the existence of samples (from deltaic deposits) with high ratio of transformation (EOM) (> 100 mg bitumen/g TOC). The GC and GC-MS analyses of the biomarkers indicate mainly the occurrence of terrestrial organic matter reflecting oxidizing conditions and both immature and very early mature stages. The results of the Rock-Eval pyrolysis and the distribution of the isoprenoids support the assumption of the input of an organic matter mixture.     

Late Quaternary activity and dextral strike-slip movement on the Karakax Fault Zone, northwest Tibet

Field observations and interpretations of satellite images reveal that the westernmost segment of the Altyn Tagh Fault (called Karakax Fault Zone) striking WNW located in the northwestern margin of the Tibetan Plateau has distinctive geomorphic and tectonic features indicative of right-lateral strike-slip fault in the Late Quaternary. South-flowing gullies and N–S-trending ridges are systematically deflected and offset by up to ~ 1250 m, and Late Pleistocene–Holocene alluvial fans and small gullies that incise south-sloping fans record dextral offset up to ~ 150 m along the fault zone. Fault scarps developed on alluvial fans vary in height from 1 to 24 m. Riedel composite fabrics of foliated cataclastic rocks including cataclasite and fault gouge developed in the shear zone indicate a principal right-lateral shear sense with a thrust component. Based on offset Late Quaternary alluvial fans, ¹⁴C ages and composite fabrics of cataclastic fault rocks, it is inferred that the average right-lateral strike-slip rate along the Karakax Fault Zone is ~ 9 mm/a in the Late Quaternary, with a vertical component of ~ 2 mm/a, and that a M 7.5 morphogenic earthquake occurred along this fault in 1902. We suggest that right-lateral slip in the Late Quaternary along the WNW-trending Karakax Fault Zone is caused by escape tectonics that accommodate north–south shortening of the western Tibetan Plateau due to ongoing northward penetration of the Indian plate into the Eurasian plate.     

Tephras in lacustrine sediments of the Sarliève marsh (French Massif Central): age and preservation

The Sarliève marsh sediments (Massif Central, France) contain two tephras. The first tephra [, ca. 12 000 BP], regionally well known, enables to date the beginning of lacustrine infill to the Lateglacial. The second tephra, the ‘tephra de Sarliève’, the emitting volcano of which is unknown, would be dated to around the Early Subboreal from pollen data. This occurrence, after the discovery of the ‘tephra de Beaunit’, emphasizes that volcanic eruption(s) occurred in the ‘Chaîne des Puys’ or in the volcanic Cézallier more than 1000 years after the last known eruption (Pavin) in the ‘Chaîne des Puys’ at around 6.6/6.7 ka (5800/5900 BP). In the Sarliève piles, these tephras, well preserved in thick and more silicated deposits of deltas, were not observed in carbonated basin sediments where they were altered. The abundance of authigenic zeolites formed during the Lateglacial in restricted depocentre lacustrine waters allows us to detect initial CF1 tephra occurrence. To cite this article: A. Fourmont et al., C. R. Geoscience 338 (2006).     

A cryo-SEM study of aggregate and floc structure changes during clay settling and raking processes

The destabilization of kaolinite suspension by anionic flocculant addition occurs in three zones; free settling, hindered settling and compression which usually includes a final bed raking process in mineral processing practice. This paper reports changes in the kaolinite aggregate and floc structures in the different settling and raking zones by cryo-vitrification/cryo-SEM techniques with image analysis combining micro- and macro-flocs. Cryo-SEM images indicate that, even during free settling, fine clay particles are bridged predominantly in edge–edge (E–E) with some face–face (F–F) configurations forming single, small flocs and some chain structures. When these small flocs and chains settle into the hindered settling zone, the collision between flocs and chains results in “honeycomb” network structures formed with lateral chain-like extension. The settled bed consists of these honeycomb structures with both inter-aggregate and intra-aggregate trapped water and has relatively low bed density (e.g. < 12 wt.% for a 2 wt.% slurry). The effect of the raking process in dramatically improving thickener underflow solids has been extensively studied but the structural changes in flocs and aggregates in this process are less well defined. Raking the compression zone for 1 h at 3 rpm can release some of the trapped water in the “honeycomb” structure and the bed density for 2 wt.% slurry improves dramatically to more than 36 wt.%. Cryo-SEM illustrates the extensive restructuring of flocs from predominantly E–E to predominantly F–F in many areas. The STructural IMage ANalysis (STIMAN) software is used to combine a series of images at magnifications from 1000× to 8000×, including both macro- and micro-flocs. This structural analysis comparing the un-raked and raked bed samples gives increases in total particle area of 30% and in relative particle area of 6%. The relatively low energy rake action of the shear stress results in the disruption of the E–E chains and the honeycomb structure, partly releasing the trapped water and inducing some E–E to F–F aggregate restructuring are clearly illustrated in these results.     

Fossil fuel energy resources of Ethiopia: Coal deposits

The gravity of Ethiopian energy problem has initiated studies to explore various energy resources in Ethiopia, one among this is the exploration for coal resources. Studies confirmed the presence of coal deposits in the country. The coal-bearing sediments are distributed in the Inter-Trappean and Pre-Trap volcanic geological settings, and deposited in fluvio-lacustrine and paludal environments in grabens and half-grabens formed by a NNE–SSW and NNW–SSE fault systems. Most significant coal deposits are found in the Inter-Trappean geological setting. The coal and coal-bearing sediments reach a maximum thickness of 4 m and 300 m, respectively. The best coal deposits were hosted in sandstone–coal–shale and mudstone–coal–shale facies. The coal formations of Ethiopia are quite unique in that they are neither comparable to the coal measures of the Permo-Carboniferous Karroo Formation nor to the Late Devonian–Carboniferous of North America or Northwestern Europe. Proximate analysis and calorific value data indicated that the Ethiopian coals fall under lignite to high volatile bituminous coal, and genetically are classified under humic, sapropelic and mixed coal. Vitrinite reflectance studies confirmed 0.3–0.64% Ro values for the studied coals. Palynology studies confirmed that the Ethiopian coal-bearing sediments range in age from Eocene to Miocene. A total of about 297 Mt of coal reserve registered in the country. The coal reserve of the country can be considered as an important alternative source of energy.     

Fashion and phases of late Pleistocene aggradation and incision in the Alaknanda River Valley, western Himalaya, India

We study the aggradation and incision of the Alaknanda River Valley during the late Pleistocene and Holocene. The morphostratigraphy in the river valley at Deoprayag shows the active riverbed, a cut terrace, and a fill terrace. The sedimentary fabric of the fill terrace comprises four lithofacies representing 1) riverbed accretion, 2) locally derived debris fan, 3) the deposits of waning floods and 4) palaeoflood records. The sedimentation style, coupled with geochemical analysis and Optically Stimulated Luminescence (OSL) dating, indicate that this terrace formed in a drier climate and the river valley aggraded in two phases during 21–18 ka and 13–9 ka. During these periods, sediment supply was relatively higher. Incision began after 10 ka in response to a strengthened monsoon and aided by increase of the tectonic gradient. The cut terrace formed at ~ 5 ka during a phase of stable climate and tectonic quiescence. The palaeoflood records suggest wetter climate 200–300 yr ago when the floods originated in the upper catchment of the Higher Himalaya and in the relatively drier climate ~ 1.2 ka when locally derived sediments from the Lesser Himalaya dominated flood deposits. Maximum and minimum limits of bedrock incision rate at Deoprayag are 2.3 mm/a and 1.4 mm/a.     

Antimony quartz and antimony–gold quartz veins from northern Portugal

Antimony- and Pb–Sb-quartz veins from the Bragança district, Portugal, are mainly hosted by Silurian phyllites. Antimony–Au-quartz veins from the Dúrico–Beirã region are mainly hosted by a Cambrian schist–metagraywacke complex, as well as Ordovician phyllites and quartzites. The deposits were mostly exploited in the late 19th Century. Mineralogical characteristics and chemical compositions of individual ore minerals are similar in the two areas. First and second generations of arsenopyrite precipitated at 390 and 300 °C, respectively. Berthierite and stibnite are the most abundant Sb-bearing minerals and precipitated between 225 and 128 °C, native antimony at < 200 °C. Drastic fluid cooling is the main cause of mineral precipitation. The Pb isotope compositions of stibnite suggest a homogeneous crustal source of lead, from the metasedimentary sequences, for Sb, Pb–Sb and Sb–Au deposits in both areas, which is consistent with the findings for comparable mineralizations elsewhere in Europe. Remobilization of Pb is related to Variscan metamorphism and deformation.     

Seismic microzonation for Penang using geospatial contour mapping

Shear wave velocity (V _s) and the fundamental site period of the subsurface condition are the primary parameters that affect seismic soil amplification in particular sites. Within the topmost layer of the soil, which measures 30 m, the average shear wave velocity V _s30 is commonly used to build codes for site classification for the design of earthquake-resistant structures and to conduct microzonation studies. In this study, the development of a microzonation map for V _s30 distribution, National Earthquake Hazard Reduction Program V _s30 site classification, and a fundamental site period for Penang are presented. The multichannel analysis of surface wave (MASW) test was conducted for more than 50 sites with available borehole data to develop the microzonation maps. The ten selected V _s profiles measured by MASW show a good correlation with the data obtained using empirical correlations in a previous study. The highest V _s values were identified at the northeastern and southeastern parts of Penang Island, corresponding to the shallow bedrock and the outcrop zone. Conversely, the lowest V _s values were found in the northwestern and southwestern parts of the Penang mainland owing to the thick layer of soft clay and silt deposits. The site period map shows the variation in site periods, with the highest value of 1.03 s at the western part of the Penang mainland and the lowest value of 0.02 s at the eastern part of the Penang Island. The microzonation maps developed in this study are vital to studies on seismic hazard and earthquake mitigation programs in Malaysia.     

Lithospheric flexure, uplift and expected horizontal strain rate in the Pannonian Carpathian region

The Apuseni Mountains are located between the Pannonian Basin and the Transylvanian Basin along a direction of SE convergence with the Carpathian belt. A flexural model based on the cylindrical bending of a semi-infinite, isostatically supported, thin elastic plate is here examined with the Apuseni playing the role of flexural bulge, and under the assumption that the plate is deforming under the action of a vertical shear force and a bending moment applied at the end of the plate, beneath the Carpathians. The model yields estimates of the plate thickness ranging between 13 and 14.5 km, depending on the assumed density contrast between crust/sediments and mantle providing buoyancy. The vertical shear force which is necessary to bend the plate is in the range between 60 and 300 × 10¹¹ N m^− 1, depending on the assumed density contrast. This force is shown to be modelled by a gravitational ‘slab pull’ force, using model parameters derived from seismic tomography. If the height of the flexural bulge, after correction for erosion, is allowed to increase, the model yields an estimate of the horizontal strain rate at the top of the bulge. For example, 5 mm/yr vertical change of the flexural bulge of a 14 km thick plate results in a horizontal deformation rate of approximately 7 nanostrain/yr at the top of the bulge, a value which is at the threshold of sensitivity of continuous GPS measurements. Different vertical rates will change the horizontal strain rate almost proportionally.     

Relationships between hydrodynamic parameters and grain size in two contrasting transitional environments: The Lagoons of Venice and Cabras, Italy

A comparison was made of shallow water sediments from the Lagoon of Venice (LV) and the Lagoon of Cabras (LC), comparing depositional environments and exploring the relationships between hydrodynamics and sedimentological parameters. The two water bodies are very different in size (LV: 360 km²; LC: 22 km²), and the sediments predominantly consist of silty-clay (LV: Mz ≈ 26 μm; LC: Mz ≈ 6 μm). However, there are large differences between the two lagoons with respect to sand (LV: mean 19%; LC: mean ~ 3%) and clay (LV: mean 20%; LC: mean 45%) contents. The Lagoon of Venice (mean depth ~ 1 m) can be considered a tidal basin, whereas the Lagoon of Cabras (mean depth ~ 2 m) has the character of a coastal lake in which wind is the main hydrodynamic forcing factor. A comparison of sediment grain-size distributions with water circulation patterns in different parts of the lagoons highlighted some interesting differences. Grain-size analyses of samples reveal a deficiency of particles around 8 μm in the LC, which is interpreted as reflecting the transition between cohesive flocs/aggregates and non-cohesive coarser silt particles, while the transition limit in the LV is ~ 20 μm. Thus, particles are cohesive below 8 μm in the LC and below ~ 20 μm in the LV. This is probably because of the differences in the clay/silt ratio, which is much lower in the LV (~ 0.3) than in LC (~ 1), conferring a “silt-dominated network structure” on most of the LV sediments.The hydrographical data used were root mean square velocity (RMSV) and water residence time (WRT), computed under the main forcing conditions. The results show a general correlation between RMSV and sortable silt in the LC, and between RMSV and coarser sediments (63–105 μm) in the LV. Some significant differences between the lagoons were detected in the degree of correlation between WRT and grain size. Root mean square velocity (~ 7 cm s^− 1 in the LV and ~ 3 cm s^− 1 in the LC) was a greater forcing factor in the LC than in the LV. Conversely, WRT, which is on average ~ 16 days in the LV and ~ 19 days in the LC, has more influence in the LV. This study highlights the usefulness of comparing environments with different hydrodynamic energies, e.g., tidal and/or wind-driven currents, to elucidate and thereby improve our understanding of the processes governing the spatial distribution of sedimentological features, the transport mechanisms of sediments, and the relationship between them. The results demonstrate that the approach outlined in this study has the potential to provide a universal hydro-sedimentological classification scheme.     

The Gibraltar Arc seismogenic zone (part 2): Constraints on a shallow east dipping fault plane source for the 1755 Lisbon earthquake provided by tsunami modeling and seismic intensity

The Great Lisbon earthquake has the largest documented felt area of any shallow earthquake and an estimated magnitude of 8.5–9.0. The associated tsunami ravaged the coast of SW Portugal and the Gulf of Cadiz, with run-up heights reported to have reached 5–15 m. While several source regions offshore SW Portugal have been proposed (e.g.— Gorringe Bank, Marquis de Pombal fault), no single source appears to be able to account for the great seismic moment as well as all the historical tsunami amplitude and travel time observations. A shallow east dipping fault plane beneath the Gulf of Cadiz associated with active subduction beneath Gibraltar, represents a candidate source for the Lisbon earthquake of 1755.Here we consider the fault parameters implied by this hypothesis, with respect to total slip, seismic moment, and recurrence interval to test the viability of this source. The geometry of the seismogenic zone is obtained from deep crustal studies and can be represented by an east dipping fault plane with mean dimensions of 180 km (N–S) × 210 km (E–W). For 10 m of co-seismic slip an Mw 8.64 event results and for 20 m of slip an Mw 8.8 earthquake is generated. Thus, for convergence rates of about 1 cm/yr, an event of this magnitude could occur every 1000–2000 years. Available kinematic and sedimentological data are in general agreement with such a recurrence interval. Tsunami wave form modeling indicates a subduction source in the Gulf of Cadiz can partly satisfy the historical observations with respect to wave amplitudes and arrival times, though discrepancies remain for some stations. A macroseismic analysis is performed using site effect functions calculated from isoseismals observed during instrumentally recorded strong earthquakes in the region (M7.9 1969 and M6.8 1964). The resulting synthetic isoseismals for the 1755 event suggest a subduction source, possibly in combination with an additional source at the NW corner of the Gulf of Cadiz can satisfactorily explain the historically observed seismic intensities. Further studies are needed to sample the turbidites in the adjacent abyssal plains to better document the source region and more precisely calibrate the chronology of great earthquakes in this region.     

Thermal maturity of a fold–thrust belt based on vitrinite reflectance analysis in the Western Foothills complex, western Taiwan

Burial depth, cumulative displacement, and peak temperature of frictional heat of a fault system are estimated by thermal analysis in the fold–thrust belt of the Western Foothills complex, western Taiwan based on the vitrinite reflectance technique. The regional thermal structure across the complex reveals that the rocks were exposed to maximum temperatures ranging from 100 °C to 180 °C, which corresponds to a burial depth of 3.7–6.7 km. A large thermal difference of 90 °C were observed at the Shuilikeng fault which make the eastern boundary of the fold–thrust belt where it is in contact with metamorphic rock of Hsuehshan Range. The large thermal difference corresponds to cumulative displacements on the Shuilikeng fault estimated to be in the range of 5.2–6.9 km. However, thermal differences in across the Shuangtung and Chelungpu faults cannot be determined apparently due to small vertical offsets. The large displacement observed across the Shuilikeng fault is absent at the other faults which are interpreted to be younger faults within the piggyback thrust system. Localized high temperatures adjacent to fault zones were observed in core samples penetrating the Chelungpu fault. Three major fracture zones were observed at core lengths of 225 m, 330 m, and 405 m and the two lower zones which comprise dark gray narrow shear zones. A value of vitrinite reflectance of 1.8%, higher than the background value of 0.8%, is limited at a narrow shear zone of 1 cm thickness at the fracture zone at 330 m. The estimated peak temperature in the range of 550–680 °C in the shear zone is far higher than the background temperature of 130 °C, and it is interpreted as due to frictional heating during seismic faulting.