Influence of the Neotethys rifting on the development of the Dampier Sub-basin (North West Shelf of Australia), highlighted by subsidence modelling

During the Late Palaeozoic and the Mesozoic, the development and evolution of the North West Shelf of Australia have been mostly driven by rifting phases associated with the break-up of Gondwana. These extensional episodes, which culminated in the opening of the Neotethys Ocean during the Permo-Carboniferous and a series of abyssal plains during the Jurassic-Cretaceous, are characterised by different stress regimes and modes of extension, and therefore had distinctive effects on the margin, and particularly on the Northern Carnarvon Basin.Interpretation of 3D and 2D seismic data enables a structural and stratigraphic analysis of the Late Palaeozoic sediments deposited in the proximal part of the Dampier Sub-basin (Mermaid Nose). Based on their seismic characters, stratigraphic relationship, internal patterns, lateral continuity, and architecture, these units are associated here with the Pennsylvanian?–Early Sakmarian glaciogenic Lyons Group and the Sakmarian–Artinskian Callytharra Formation. The former were deposited in a half-graben whose development is associated with the onset of the Neotethys rifting, and the latter is characterised by restricted deposition, inversion of prograding patterns, and uplift.The integration of seismo-stratigraphic characterisation of the Late Palaeozoic sequences and Mesozoic data from one exploration well (Roebuck-1) enables the construction of subsidence curves for the Mermaid Nose and the interpretation of its geohistory.The tectonic subsidence curves show a striking Permo-Carboniferous rifting phase related to the Neotethys rifting and a discrete Late Jurassic–Early Cretaceous event coeval with the opening and the spreading of the Argo Abyssal Plain.This result points out the predominance of the effects of the Permo-Carboniferous Neotethys episode, whereas the extension related to the Argo Abyssal Plain rifting that occurred later and closer to the studied area, had only limited effects on the subsidence of the proximal Dampier Sub-basin. Therefore, it supports a tectonic model with two distinct modes of extension for the Late Palaeozoic (widespread) and the Mesozoic (localised) rifting phases.     

Some considerations concerning seismic geomorphology and paleoseismology

Seismic geomorphology studies landforms which developed in connection with earthquakes. Among them, two different end members may be distinguished: 1) seismo-tectonic landforms, including surface faults and fractures, land uplift and subsidence at different scales, surface bulges, elongate ridges, and any other permanent ground deformations directly related to tectonic stress, and 2) seismo-gravitational landforms, such as landslides, deep-seated gravitational slope deformations, sinkholes, and fissures due to sediment compaction or liquefaction and sand blows, connected with both seismic shaking and gravitational stress.A clear-cut distinction between the two categories of landforms is not always easy to make (and in many instances not really useful), while there are, in many cases, ground effects that might be (and should be) considered as simultaneous combinations of seismo-tectonic and seismo-gravitational processes. This applies especially to surface fracturing and faulting which could be the combined result of tectonic stress, stress produced by seismic shaking, and gravitational stress.The objective of this paper is to review selected case histories mainly from Italy and the Mediterranean region, in order to show the importance of a comprehensive study of earthquake-generated landforms for understanding the seismicity level of the area under investigation. We argue that in earthquake prone areas, seismic landforms often constitute typical patterns (seismic landscapes) whose recognition, mapping and paleoseismological analysis may help in the evaluation of seismic hazards.     

Anomalous lithospheric structure of Northern Juan de Fuca plate — a consequence of oceanic rift propagation?

Anomalous crustal and upper mantle structure of northern Juan de Fuca plate is revealed from wide-angle seismic and gravity modelling. A 2-D velocity model is produced for refraction line II of the 1980 Vancouver Island Seismic Project (VISP80). The refraction data were recorded on three ocean bottom seismometers (OBSs) deployed at the ends and middle of a 110 km line oriented parallel to the North American continental margin. The velocity model is constructed via ray tracing and conforms to first-arrival amplitude observations and travel time picks of direct, converted and reflected phases. Between sub-sediment depths of 3 to 11 km, depths normally associated with the lower crust and upper oceanic mantle, the final model shows that compressional-wave velocities decrease significantly from southeast to northwest along the profile. At sub-sediment depths of 11 km at the northwestern end of the profile, P-wave velocities are as low as 7.2 km/s. A complementary 2-D gravity model using the geometry of the velocity model and velocity–density relationships characteristic of oceanic crust is produced. The high densities required to match the gravity field indicate the presence of peridotites containing 25–30% serpentine by volume, rather than excess gabbroic crust, within the deep low velocity zone. Anomalous travel time delays and unusual reflection characteristics observed from proximal seismic refraction and reflection experiments suggest a broader zone of partially serpentinized peridotites coincident with the trace of a pseudofault. We propose that partial serpentinization of the upper mantle is a consequence of slow spreading at the tip of a propagating rift.     

Directional dependence of P- and S-wave propagation and polarization in foliated rocks from the Kola superdeep well: Evidence from laboratory measurements and calculations based on TOF neutron diffraction

We have measured P- and S-wave velocities on two amphibolite and two gneiss samples from the Kola superdeep borehole as a function of pressure (up to 600 MPa) and temperature (up to 600 °C). The velocity measurements include compressional (Vp) and shear wave velocities (Vs1, Vs2) propagating in three orthogonal directions which were in general not parallel to inherent rock symmetry axes or planes. The measurements are accompanied by 3D-velocities calculations based on lattice preferred orientation (LPO) obtained by TOF (Time Of Flight) neutron diffraction analysis which allows the investigation of bulk volumes up to several cubic centimetres due to the high penetration depth of neutrons. The LPO-based numerical velocity calculations give important information on the different contribution of the various rock-forming minerals to bulk elastic anisotropy and on the relations of seismic anisotropy, shear wave splitting, and shear wave polarization to the structural reference frame (foliation and lineation). Comparison with measured velocities obtained for the three propagation directions that were not in accordance with the structural frame of the rocks (foliation and lineation) demonstrate that for shear waves propagating through anisotropic rocks the vibration directions are as important as the propagation directions. The study demonstrates that proper measurement of shear wave splitting by means of two orthogonal polarized sending and receiving shear wave transducers is only possible when their propagation and polarization directions are parallel and normal to foliation and lineation, respectively.     

Seismic velocities of granulite-facies xenoliths from Central Ireland: Implications for lower crustal composition and anisotropy

V_p and V_s values have been measured experimentally and calculated for granulite-facies lower crustal xenoliths from central Ireland close to the Caledonian Iapetus suture zone. The xenoliths are predominantly foliated and lineated metapelitic (garnet–sillimanite–K-feldspar) granulites. Their metapelitic composition is unusual compared with the mostly mafic composition of lower crustal xenoliths world-wide. Based on thermobarometry, the metapelitic xenoliths were entrained from depths of c. 20–25 ± 3.5 km and rare mafic granulites from depths of 31–33 ± 3.4 km. The xenoliths were emplaced during Lower Carboniferous volcanism and are considered to represent samples of the present day lower crust.V_p values for the metapelitic granulites range between 6.26 and 7.99 km s^{− 1} with a mean value of 7.09 ± 0.4 km s^{− 1}. Psammite and granitic orthogneiss samples have calculated V_p values of 6.51 and 6.23 km s^{− 1}, respectively. V_s values for the metapelites are between 3.86 and 4.34 km s^{− 1}, with a mean value of 4.1 ± 0.15 km s^{− 1}. The psammite and orthogneiss have calculated V_s values of 3.95 and 3.97 km s^{− 1}, respectively.The measured seismic velocities correlate with density and with modal mineralogy, especially the high content of sillimanite and garnet. V_p anisotropy is between 0.15% and 13.97%, and a clear compositional control is evident, mainly in relation to sillimanite abundance. Overall V_s anisotropy ranges from 1% to 11%. Poisson's ratio (σ) lies between 0.25 and 0.35 for the metapelitic granulites, mainly reflecting a high V_p value due to abundant sillimanite in the sample with the highest σ. Anisotropy is probably a function of deformation associated with the closure of the Iapetus ocean in the Silurian as well as later extension in the Devonian. The orientation of the bulk strain ellipsoid in the lower crust is difficult to constrain, but lineation is likely to be NE–SW, given the strike-slip nature of the late Caledonian and subsequent Acadian deformation.When corrected for present-day lower crustal temperature, the experimentally determined V_p values correspond well with velocities from the ICSSP, COOLE I and VARNET seismic refraction lines. Near the xenolith localities, the COOLE I line displays two lower crustal layers with in situ V_p values of 6.85–6.9 and 6.9–8.0 km s^{− 1}, respectively. The upper (lower velocity) layer corresponds well with the metapelitic granulite xenoliths while the lower (higher velocity) layer matches that of the basic granulite xenoliths, though their metamorphic pressures suggest derivation from depths corresponding to the present-day upper mantle.     

Clustering and visualization of earthquake data in a grid environment

We present a web client-server service WEB-IS, which we have developed for remote analysis and visualization of seismic data consisting of both small magnitude events and large earthquakes. We show that the problem-solving environment (PSE) intended for prediction of large magnitude earthquakes can be based on this WEB-IS idea. The clustering schemes, feature generation, feature extraction techniques and rendering algorithms form a computational framework of this environment. On the other hand, easy and fast access both to the seismic data distributed among distant computing resources and to computational and visualization resources can be realized in a GRID framework. We discuss the usefulness of NaradaBrokering (iNtegrated Asynchronous Real-time Adaptive Distributed Architecture) as a middleware, allowing for flexibility and high throughput for remote visualization of geophysical data. The WEB-IS functionality was tested both on synthetic and the actual earthquake catalogs. We consider the application of similar methodology for tsunami alerts.     

Liquefaction potential evaluations by energy-based method and stress-based method for various ground motions

An energy-based liquefaction potential evaluation method (EBM) previously developed was applied to a uniform sand model shaken by seismic motions recorded at different sites during different magnitude earthquakes. It was also applied to actual liquefaction case histories in Urayasu city during the 2011 M9.0 Tohoku earthquake and in Tanno-cho during the 2003 M8.0 Tokachi-oki earthquake. In all these evaluations, the results were compared with those by the currently used stress-based method (SBM) under exactly the same seismic and geotechnical conditions. It was found that EBM yields similar results with SBM for several ground motions of recent earthquakes but has easier applicability without considering associated parameters. In Urayasu city, the two methods yielded nearly consistent results by using an appropriate coefficient in SBM for the M9.0 earthquake, though both overestimated the actual liquefaction performance, probably because effects of plasticity and aging on in situ liquefaction strength were not taken into account. In Tanno-cho, EBM could evaluate actual liquefaction performance due to a small-acceleration motion during a far-field large magnitude earthquake while SBM could not.     

Simplified multi-block constitutive model predicting earthquake-induced landslide triggering and displacement along slip surfaces of saturated sand

Slopes consisting of saturated sand have recently moved down-slope tens or hundreds of meters under the action of earthquakes. This paper presents a simplified but accurate method predicting the triggering and displacement of such landslides. For this purpose, a simplified constitutive model simulating soil response of saturated sands along slip surfaces is proposed and validated. Then, this constitutive model is coupled with the multi-block sliding system model to predict the triggering and displacement of such slides. The multi-block model considers a general mass sliding on a trajectory which consists of n linear segments. The steps needed to apply this method are described in detail. The method was applied successfully to predict the triggering, the motion and the final configuration of the well-documented (a) Higashi Takezawa, (b) Donghekou and (c) Nikawa earthquake-induced slides.     

Geodynamic evolution of the Salinas de Añana diapir in the Basque-Cantabrian Basin,Western Pyrenees

The Salinas de Añana diapir is located in the Basque-Cantabrian basin part of the great evaporite basin, along with the Gulf of Mexico and the Central European basin, when the fragmentation of Pangea started. The evolution of these basins can only be achieved by understanding the control of salt in the sedimentary and tectonic evolution of these basins.Sedimentation began with clastic Buntsandstein sediments and minor Muschelkalk limestones. Subsequent Keuper evaporites are the bottom of sedimentary cover constituted by Jurassic limestones and marls, a clastic Lower Cretaceous and an alternant limestone and marl Upper Cretaceous, whose deposition has been conditioned by salt tectonics. The emplacement of salt extends from the Aptian until now, favored by the duplication of the salt thickness associated with the thrust of Sierra Cantabria, so it is an excellent example to study changes in the regime of intrusion along the time. The geodynamic evolution of the Salinas de Añana diapir was determined through the interpretation of 45 reprocessed seismic lines, along with information from three wells. Migration of the salt in this diapir, conditioned by N120E and N30E pre-Alpine basement lineations, was determined using time isopach maps of the various rock layers. Vertical evolution of the diapir was determined through the reconstruction of a north-south section at various geologic times by flattening the respective seismic horizons. A minimum of salt flow into the diapir coincides with a minimum rate of sedimentation during the Turonian. Similarly, maximum flows of salt into the diapir occurred during the Coniacian and Lower Santonian and again from the end of the Lower Miocene to the present, coinciding with maximum rates of sedimentation during these times. In the Tertiary, probably during the Oligocene, the diapir was displaced to the south by the Sierra Cantabria thrust, maintaining the contact between the evaporites of diapir and the same evaporites of the lower block. Since the Oligocene, the salts of the lower block migrated towards and into the diapir, deforming the trace of the overthrust.