ALBERT

Description: The Piano di Pezza fault is the central section of the 35 km long L'Aquila-Celano active normal fault-system in the central Apennines of Italy. Although palaeoseismic data document high Holocene vertical slip rates (~1 mm yr –1 ) and a remarkable seismogenic potential of this fault, its subsurface setting and Pleistocene cumulative displacement are still poorly known. We investigated for the first time the shallow subsurface of a key section of the main Piano di Pezza fault splay by means of high-resolution seismic and electrical resistivity tomography coupled with time-domain electromagnetic soundings (TDEM). Our surveys cross a ~5-m-high fault scarp that was generated by repeated surface-rupturing earthquakes displacing Holocene alluvial fans. We provide 2-D Vp and resistivity images, which show significant details of the fault structure and the geometry of the shallow basin infill material down to 50 m depth. Our data indicate that the upper fault termination has a sub-vertical attitude, in agreement with palaeoseismological trench evidence, whereas it dips ~50° to the southwest in the deeper part. We recognize some low-velocity/low-resistivity regions in the fault hangingwall that we relate to packages of colluvial wedges derived from scarp degradation, which may represent the record of some Holocene palaeo-earthquakes. We estimate a ~13–15 m throw of this fault splay since the end of the Last Glacial Maximum (~18 ka), leading to a 0.7–0.8 mm yr –1 throw rate that is quite in accordance with previous palaeoseismic estimation of Holocene vertical slip rates. The 1-D resistivity models from TDEM soundings collected along the trace of the electrical profile significantly match with 2-D resistivity images. Moreover, they indicate that in the fault hangingwall, ~200 m away from the surface fault trace, the pre-Quaternary carbonate basement is at ~90–100 m depth. We therefore provide a minimal ~150–160 m estimate of the cumulative throw of the Piano di Pezza fault system in the investigated section. We further hypothesize that the onset of the Piano di Pezza fault activity may date back to the Middle Pleistocene (~0.5 Ma), so this is a quite young active normal fault if compared to other mature normal fault systems active since 2–3 Ma in this portion of the central Apennines.

Description: We present the crustal resistivity structure of the Pamir and Southern Tian Shan orogenic belts at the northwestern promontory of the India–Asia collision zone. The magnetotelluric (MT) data were recorded along a roughly north–south trending, 350 km long corridor from the Pamir Plateau in southern Tajikistan across the Pamir frontal ranges, the Alai Valley and the southwestern Tian Shan to Osh in the Kyrgyz part of the Fergana Basin. In total, we measured at 178 sites, whereof 26 combine broad band and long period recordings. One of the most intriguing features of the 2-D and 3-D inversion results is a laterally extended zone of high electrical conductivity below the Pamir Plateau, with resistivities below 1 m, starting at a depth of ~10–15 km. The high conductivity can be explained with the presence of partially molten rocks at middle to lower crustal levels, possibly related to ongoing migmatization and/or middle/lower crustal flow underneath the Southern Pamir. This interpretation is consistent with a low velocity zone found from local earthquake tomography, relatively high v p / v s ratios, elevated surface heat flow, and thermomechanical modelling suggesting that melting temperatures are reached in the felsic middle crust. In the upper crust of the Pamir and Tian Shan, the Palaeozoic–Mesozoic suture zones appear as electrically conductive, whereas the compact metamorphic rocks of the Muskol-Shatput Dome of the Central Pamir are highly resistive. The intra-montane basin of the Alai Valley—sandwiched between the Pamir and Tian Shan—exhibits a generally conductive upper crust that bifurcates into two conductors at depth. One of them connects to the active Main Pamir Thrust, which is absorbing most of today's convergence between the Pamir and the Tian Shan. Several deeper zones of high conductivity in the middle and lower crust of Central and Northern Pamir likely record fluid release due to metamorphism associated with active continental subduction/delamination.

Description: Described are results of laboratory experiments which revealed regularities of gradual transition from stick-slip mode to aseismic creep. The behaviour of model gouge-filled fault was investigated with experimental setup of the spring-bock model. It was experimentally proven that small variations of a percentage of materials with velocity strengthening and velocity weakening properties in the fault principal slip zone may result in significant variation of the portion of seismic energy radiated during a fault slip event. The tests simulated different modes of interblock sliding whose characteristic values of scaled kinetic energy varied by several orders of magnitude, while differences in contact strength and shear stress drop remained relatively small. The obtained results led to the conclusion that the earthquake radiation efficiency and the fault slip mode are governed by the ratio of two parameters—maximum fault slip-weakening rate and shear stiffness of the enclosing massif. The ratio can be essentially changed by small variations of the material composition of the fault principal slip zone.

Description: The Scythian Platform (ScP) with a heterogeneous basement of Baikalian–Variscan–Cimmerian age is located between the East European Craton (EEC) on the north and the Crimean–Caucasus orogenic belt and the Black Sea (BS) Basin on the south. In order to get new constrains on the basin architecture and crustal structure of the ScP and a better understanding of the tectonic processes and evolution of the southern margin of the EEC during Mesozoic and Cenozoic time, a 630-km-long seismic wide-angle refraction and reflection (WARR) profile DOBRE-5 was acquired in 2011 October. It crosses in a W–E direction the Fore-Dobrudja Trough, the Odessa Shelf of the BS and the Crimean Plain. The field acquisition included eight chemical shot points located every 50 km and recorded by 215 stations placed every ~2.0 km on the land. In addition, the offshore data from existing profile 26, placed in the Odessa Shelf, were used. The obtained seismic model shows clear lateral segmentation of the crust within the study region on four domains: the Fore-Dobrudja Domain (km 20–160), an offshore domain of the Karkinit Trough at the Odessa Shelf of the BS (km 160–360), an onshore domain of the Central Crimean Uplift (Crimean Plain, km 360–520) and the Indolo-Kuban Trough at the Kerch Peninsula (km 520–620) that is the easternmost part of the Crimea. Two contrasting domains of the ScP within the central part of the DOBRE-5 profile, the Karkinit Trough and the Central Crimean Uplift, may represent different stages of the ScP formation. A deep Karkinit Trough with an underlying high-velocity (〉7.16 km s –1 ) lower crust body suggests its rifting-related origin during Early Cretaceous time. The Central Crimean Uplift represents a thick (up to 47 km) crustal domain consisting of three layers with velocities 5.8–6.4, 6.5–6.6 and 6.7–7.0 km s –1 , which could be evidence of this part of the ScP originating on the crust of Precambrian craton (EEC). The thick heterogeneous basement of the Central Crimean Uplift shows inclusions of granitic bodies associated with magmatic activity related with Variscan orogeny within the ScP. General bending and crustal scale buckling of the Central Crimean Uplift with a wavelength of 230 km could be an effect of the Alpine compressional tectonics in the adjacent Crimean Mountains. The extended/rifted continental margin of the ScP (EEC) at the Odessa Shelf and buckling/uplifted domain of the Central Crimean Uplift affected by compressional tectonics, are separated by the N–S oriented Western Crimean Fault. The crust of the southern margin of the EEC is separated from the ScP, which originated on the EEC crust tectonised and reworked during the Palaeozoic–Mesozoic, by the crustal fault of ~W–E orientation, which corresponds with the Golitsyn Fault observed at the surface between the EEC and the ScP. The Fore-Dobrudja Domain with a thick (〉10 km) heterogeneous basement and two subhorizontal layers in the crystalline crust (with velocities 6.2–6.3 and 6.4–6.65 km s –1 ) differs from the ScP crust and its origin could be very similar to that of the Trans-European Suture Zone and Palaeozoic West European Platform.

Description: We develop and validate a systematic approach to infer plate boundary strength and rheological parameters in models of mantle flow from surface velocity observations. Based on a realistic rheological model that includes yielding and strain rate weakening from dislocation creep, we formulate the inverse problem in a Bayesian inference framework. To study the distribution of parameters that are consistent with the observations, we compute the maximum a posteriori (MAP) point, Gaussian approximations of the parameter distribution around that MAP point, and employ Markov Chain Monte Carlo (MCMC) sampling methods. The computation of the MAP point and the Gaussian approximation require first and second derivatives of an objective function subject to non-linear Stokes equations; these derivatives are computed efficiently using adjoint Stokes equations. We set up 2-D numerical experiments with many of the elements expected in a global geophysical inversion. This setup incorporates three subduction zones with slab and weak zone (interplate fault) geometry consistent with average seismic characteristics. With these experiments, we demonstrate that when the temperature field is known, we can recover the strength of plate boundaries, the yield stress and strain rate exponent in the upper mantle. When the number of uncertain parameters increases, there are trade-offs between the inferred parameters. These trade-offs depend on how well the observational data represents the surface velocities, and on the weakness of plate boundaries. As the plate boundary coupling drops below a threshold, the uncertainty of the inferred parameters increases due to insensitivity of plate motion to plate coupling. Comparing the trade-offs between inferred rheological parameters found from the Gaussian approximation of the parameter distribution and from MCMC sampling, we conclude that the Gaussian approximation—which is significantly cheaper to compute—is often a good approximation, in particular locally around the MAP point. Thus, the method can be applied to the global problem of inferring non-linear constitutive parameters and plate coupling factors for each subduction zone in a global geophysical inversion with known slab structure.

Description: In Europe, common input data types for seismic hazard evaluation include earthquake catalogues, seismic zonation models and ground motion models, all with well-constrained epistemic uncertainties. In contrast, neotectonic deformation models and their related uncertainties are rarely considered in earthquake forecasting and seismic hazard studies. In this study, for the first time in Europe, we developed a seismic hazard model based exclusively on active fault and geodynamic deformation models. We applied it to the External Dinarides, a slow-deforming fold-and-thrust belt in the Central Mediterranean. The two deformation models furnish consistent long-term earthquake rates above the M w 4.7 threshold on a latitude/longitude grid with 0.2° spacing. Results suggest that the use of deformation models is a valid alternative to empirical-statistical approaches in earthquake forecasting in slow-deforming regions of Europe. Furthermore, we show that the variability of different deformation models has a comparable effect on the peak ground motion acceleration uncertainty as do the ground motion prediction equations.

Description: The Zambales Ophiolite Complex (ZOC) in the western portion of Luzon Island, Philippines represents a typical exposure of an emplaced crust-upper mantle section of an ancient lithosphere. The ZOC is divided into the Acoje and Coto Blocks based on petrological, geochemical and age disparities, thus implying diverse and complex origins for the archipelago's lithospheric sources. We used gravity and magnetic data to reveal differences in the subsurface characteristics of the two ophiolite blocks for the first time. Low Bouguer gravity (〈135 mGal) and magnetic (〈69 nT) anomalies characterize the ophiolitic units in the Acoje Block whereas high Bouguer gravity (〉150 mGal) and magnetic (〉110 nT) anomalies typify the Coto Block. Such contrasting signatures further extend at greater depths which reflect the disparities in the crustal density, the basement structure, and the depth to Moho. Petrophysical characteristics such as density, magnetic susceptibility and natural remanent magnetization (NRM) intensities correlate well with the interpreted Bouguer gravity and magnetic anomalies. Densities of gabbros and peridotites from the Acoje and Coto Blocks reveal mean values ranging from 2640 to 2810 kg m –3 and 2570 to 2690 kg m –3 , respectively. Magnetic susceptibility (〉4.43 x 10 –3 SI) and NRM (〉0.69 A m –1 ) data are also generally higher over the Coto Block. Both Bouguer gravity and magnetic anomaly maps reveal a prominent steep gradient that potentially marks the structural contact between Acoje and Coto Blocks. We further infer that the steep anomaly gradient validates the presence of the Lawis Fault Zone which separates the two ophiolite blocks in the Masinloc Massif. Recent field evidence from rock exposures in Coto reveals both right lateral and vertical displacements along the fault zone.

Description: The North Anatolian Fault Zone (NAFZ) starts to branch off in the western Bolu plain. The branches of the NAFZ in this location create the Almacık block which is surrounded by the latest surface ruptures of significant earthquakes that occurred between 1944 and 1999, but its northeastern part remains unruptured. The most recently formed rupture, that was a result of the 1999 November 12 Düzce earthquake, ended to the northwest of the Bakacak Fault. The connection between the Bakacak Fault and the main branch of the NAFZ via the Bolu plain has until now remained unknown. This paper establishes that the route of the missing link runs through the Dağkent, Kasaplar and Bürnük faults, a finding achieved with the help of seismic reflection studies. The paper also argues that the cross cutting nature of these newly determined faults and a stress analysis based on focal mechanism solutions of recent earthquakes demonstrate the termination of the suggested pull-apart nature of the Bolu plain.

Description: Geophysical data are the main source of information about the subsurface. Geophysical techniques are, however, highly non-unique in determining specific physical parameters and boundaries of subsurface objects. To obtain actual physical information, an inversion process is often applied, in which measurements at or above the Earth surface are inverted into a 2- or 3-D subsurface spatial distribution of the physical property. Interpreting these models into structural objects, related to physical processes, requires a priori knowledge and expert analysis which is susceptible to subjective choices and is therefore often non-repeatable. In this research, we implemented a recently introduced object-based approach to interpret the 3-D inversion results of a single geophysical technique using the available a priori information and the physical and geometrical characteristics of the interpreted objects. The introduced methodology is semi-automatic and repeatable, and allows the extraction of subsurface structures using 3-D object-oriented image analysis (3-D OOA) in an objective knowledge–based classification scheme. The approach allows for a semi-objective setting of thresholds that can be tested and, if necessary, changed in a very fast and efficient way. These changes require only changing the thresholds used in a so-called ruleset, which is composed of algorithms that extract objects from a 3-D data cube. The approach is tested on a synthetic model, which is based on a priori knowledge on objects present in the study area (Tanzania). Object characteristics and thresholds were well defined in a 3-D histogram of velocity versus depth, and objects were fully retrieved. The real model results showed how 3-D OOA can deal with realistic 3-D subsurface conditions in which the boundaries become fuzzy, the object extensions become unclear and the model characteristics vary with depth due to the different physical conditions. As expected, the 3-D histogram of the real data was substantially more complex. Still, the 3-D OOA-derived objects were extracted based on their velocity and their depth location. Spatially defined boundaries, based on physical variations, can improve the modelling with spatially dependent parameter information. With 3-D OOA, the non-uniqueness on the location of objects and their physical properties can be potentially significantly reduced.

Description: We use continuous GPS measurements from 31 stations in southern Mexico to model coseismic slip and post-seismic deformation from the 2012 March 20 M w  = 7.5 Ometepec earthquake, the first large thrust earthquake to occur below central Mexico during the modern GPS era. Coseismic offsets ranging from ~280 mm near the epicentre to 5 mm or less at sites far from the epicentre are fit best by a rupture focused between ~15 and 35 km depth, consistent with an independent seismological estimate. The corresponding geodetic moment of 1.4 10 20 N·m is within 10 per cent of two independent seismic estimates. Transient post-seismic motion recorded by GPS sites as far as 300 km from the rupture has a different horizontal deformation gradient and opposite sense of vertical motion than do the coseismic offsets. A forward model of viscoelastic relaxation as a result of our new coseismic slip solution incorrectly predicts uplift in areas where post-seismic subsidence was recorded and indicates that viscoelastic deformation was no more than a few per cent of the measured post-seismic deformation. The deformation within 6 months of the earthquake was thus strongly dominated by fault afterslip. The post-seismic GPS time-series are well fit as logarithmically decaying fault afterslip on an area of the subduction interface up to 10 times larger than the earthquake rupture zone, extending as far as 220 km inland. Afterslip had a cumulative geodetic moment of 2.0 10 20 N·m, ~40 per cent larger than the Ometepec earthquake. Tests for the shallow and deep limits for the afterslip require that it included much of the earthquake rupture zone as well as regions of the subduction interface where slow slip events and non-volcanic tremor have been recorded and areas even farther downdip on the flat interface. Widespread afterslip below much of central Mexico suggests that most of the nearly flat subduction interface in this region is conditionally stable and thus contributes measurable transient deformation to large areas of Mexico south of and in the volcanic belt.