ALBERT

Description: We present an up-to-date high resolution picture of the ongoing crustal deformation field of Italy, based on an extensive combination of permanent and non-permanent GPS observations carried out since 1994. In addition, we present an updated map of contemporary S Hmax orientations computed by a multidisciplinary data set of well-constrained stress indicators, including both published results and novel analyses. The comparison of stress and geodetic strain-rates directions reveals that both patterns are near-parallel over a large part of the investigated area, highlighting that crustal stress and surface deformation are driven by the same mechanism. The comparison of the azimuthal patterns of surface strain and mantle deformation shows a modest correlation on the Alps and a low correlation along the Apennines chain and the Calabro-Peloritan Arc. Along the Apennines chain, this feature suggests the occurrence of significant strain partitioning and crust–mantle mechanical decoupling. Along the Calabro-Peloritan Arc, the apparent low correlation reflects a different mantle–crust mechanism of deformation to the ongoing subduction and rollback of the Ionian slab. In addition, the superposition of regional/local effects related to second-order sources (crustal lateral density changes, strength contrasts), which at regional/local scale modulate the crustal stress/strain-rate pattern, cannot be ruled out.

Description: A series of linear analysis was performed on the onset of thermal convection of highly compressible fluids, in order to deepen the fundamental insights into the mantle convection of massive super-Earths in the presence of strong adiabatic compression. We consider the temporal evolution (growth or decay) of an infinitesimal perturbation superimposed to a highly compressible fluid which is in a hydrostatic (motionless) and conductive state in a basally heated horizontal layer. As a model of pressure-dependence in material properties, we employed an exponential decrease in thermal expansivity α and exponential increase in (reference) density with depth. The linearized equations for conservation of mass, momentum and internal (thermal) energy are numerically solved for the critical Rayleigh number as well as the vertical profiles of eigenfunctions for infinitesimal perturbations. The above calculations are repeatedly carried out by systematically varying (i) the dissipation number (Di), (ii) the temperature at the top surface and (iii) the magnitude of pressure-dependence in α and . Our analysis demonstrated that the onset of thermal convection is strongly affected by the adiabatic compression, in response to the changes in the static stability of thermal stratification in the fluid layer. For sufficiently large Di where a thick sublayer of stable stratification develops in the layer, for example, the critical Rayleigh number explosively increases with Di, together with drastic decreases in the length scales of perturbations both in vertical and horizontal directions. In particular, for very large Di, a thick ‘stratosphere’ occurs in the fluid layer where the vertical motion is significantly suppressed, resulting in a shrink of the incipient convection in a thin sublayer of unstable thermal stratification. In addition, when Di exceeds a threshold value above which a thermal stratification becomes stable in the entire layer, no perturbation is allowed to grow with time regardless of the Rayleigh number and/or the horizontal wavelength. We also found that the effect of adiabatic compression becomes prominent for higher temperature at the top surface of the fluid layer. These findings may imply the crucial importance of adiabatic compression in understanding the dynamics and evolution of the mantles of massive super-Earths, particularly for those orbiting their parent stars very closely.

Description: Between 2011 and 2013, two dense transects were deployed across the central and western Pyrenees to get better constraints on the deep lithospheric architecture and discriminate the competing models of the structure and formation of the Pyrenees. Each transect recorded the regional and global seismicity during a period of approximately 1 yr. Here, we exploit the records of teleseismic compressional waves and of their conversions to shear waves on internal discontinuities in order to map lithospheric interfaces beneath the two transects. The migrated sections, obtained by performing common conversion point stacks, are in remarkable agreement with the results of the ECORS-Pyrenees and ECORS-Arzacq deep seismic surveys. However, the migrations of converted waves reveal new details of the deep lithospheric architecture that could not be seen with the active source experiments. The new images provide clear and definite evidence for the subduction of a thinned Iberian crust down to at least ~70 km depth, a result that has important implications for the formation of the Pyrenees. The subduction of the Iberian lithosphere leads to reconsider the amount of convergence between Iberia and Europe during the Cenozoic. A recent regional P -wave tomography, relying on the data of the PYROPE and IBERARRAY temporary experiments, revealed the segmentation of lithospheric structures by inherited Hercynian NE–SW transfer faults that were reactivated during the Albian rifting. Our migration images are consistent with this model, and give further support to the idea that the Pyrenees were produced by the tectonic inversion of a segmented hyperextended rift that was buried by subduction beneath the European Plate.

Description: In this paper, we propose an approach to compute the coseismic Earth's volume change based on a spherical-Earth elastic dislocation theory. We present a general expression of the Earth's volume change for three typical dislocations: the shear, tensile and explosion sources. We conduct a case study for the 2004 Sumatra earthquake ( M w 9.3), the 2010 Chile earthquake ( M w 8.8), the 2011 Tohoku-Oki earthquake ( M w 9.0) and the 2013 Okhotsk Sea earthquake ( M w 8.3). The results show that mega-thrust earthquakes make the Earth expand and earthquakes along a normal fault make the Earth contract. We compare the volume changes computed for finite fault models and a point source of the 2011 Tohoku-Oki earthquake ( M w 9.0). The big difference of the results indicates that the coseismic changes in the Earth's volume (or the mean radius) are strongly dependent on the earthquakes’ focal mechanism, especially the depth and the dip angle. Then we estimate the cumulative volume changes by historical earthquakes ( M w ≥ 7.0) since 1960, and obtain an Earth mean radius expanding rate about 0.011 mm yr –1 .

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: A new model of the deglaciation history of Antarctica over the past 25 kyr has been developed, which we refer to herein as ICE-6G_C (VM5a). This revision of its predecessor ICE-5G (VM2) has been constrained to fit all available geological and geodetic observations, consisting of: (1) the present day uplift rates at 42 sites estimated from GPS measurements, (2) ice thickness change at 62 locations estimated from exposure-age dating, (3) Holocene relative sea level histories from 12 locations estimated on the basis of radiocarbon dating and (4) age of the onset of marine sedimentation at nine locations along the Antarctic shelf also estimated on the basis of 14 C dating. Our new model fits the totality of these data well. An additional nine GPS-determined site velocities are also estimated for locations known to be influenced by modern ice loss from the Pine Island Bay and Northern Antarctic Peninsula regions. At the 42 locations not influenced by modern ice loss, the quality of the fit of postglacial rebound model ICE-6G_C (VM5A) is characterized by a weighted root mean square residual of 0.9 mm yr –1 . The Southern Antarctic Peninsula is inferred to be rising at 2 mm yr –1 , requiring there to be less Holocene ice loss there than in the prior model ICE-5G (VM2). The East Antarctica coast is rising at approximately 1 mm yr –1 , requiring ice loss from this region to have been small since Last Glacial Maximum. The Ellsworth Mountains, at the base of the Antarctic Peninsula, are inferred to be rising at 5–8 mm yr –1 , indicating large ice loss from this area during deglaciation that is poorly sampled by geological data. Horizontal deformation of the Antarctic Plate is minor with two exceptions. First, O'Higgins, at the tip of the Antarctic Peninsula, is moving southeast at a significant 2 mm yr –1 relative to the Antarctic Plate. Secondly, the margins of the Ronne and Ross Ice Shelves are moving horizontally away from the shelf centres at an approximate rate of 0.8 mm yr –1 , in viscous response to the early Holocene unloading of ice from the current locations of the ice shelf centers. ICE-6G_C (VM5A) fits the horizontal observations well (wrms residual speed of 0.7 mm yr –1 ), there being no need to invoke any influence of lateral variation in mantle viscosity. ICE-6G_C (VM5A) differs in several respects from the recently published W12A model of Whitehouse et al. First, the upper-mantle viscosity in VM5a is 5 10 20 Pa s, half that in W12A. The VM5a profile, which is identical to that inferred on the basis of the Fennoscandian relaxation spectrum, North American relative sea level histories and Earth rotation constraints, when coupled with the revised ICE-6G_C deglaciation history, fits all of the available constraints. Secondly, the net contribution of Antarctica ice loss to global sea level rise is 13.6 m, 2/3 greater than the 8 m in W12A. Thirdly, ice loss occurs quickly from 12 to 5 ka, and the contribution to global sea level rise during Meltwater Pulse 1B (11.5 ka) is large (5 m), consistent with sedimentation constraints from cores from the Antarctica ice shelf. Fourthly, in ICE-6G_C there is no ice gain in the East Antarctica interior, as there is in W12A. Finally, the new model of Antarctic deglaciation reconciles the global constraint upon the global mass loss during deglaciation provided by the Barbados record of relative sea level history when coupled with the Northern Hemisphere counterpart of this new model.

Description: The Himalaya is the result of the on-going convergence and collision of India and Asia. The internal configuration and processes that govern the rise of the Himalayan Mountains and Tibetan Plateau are crucial to understand continental collision zones. However, knowledge of the prior configuration of the colliding plates is equally important, since inherited (pre-orogenic/basement) structures can undeniably influence the development of the orogenic architecture throughout the orogen's cycle of collision and eventual collapse. Three northeast-trending palaeotopographic ridges of faulted Precambrian Indian basement underlie the Ganga basin south of the Himalaya. Our paper illustrates a crustal-scale fault origin for these ridges and succeeds in determining how far north beneath the Himalayan system they extend and how they ultimately govern the location of upper crustal faults in southern Tibet. Spectrally filtered EGM2008 Bouguer gravity data and edges in its horizontal gradient at different source depths (‘gravity worms’) over northern Peninsular India, the Himalaya and southern Tibet reveal several continuous Himalayan cross-strike discontinuities interpreted to represent crustal faults. Gravity lineaments in Peninsular India coincide with edges of the Precambrian basement ridges and megakinks up to 100 km wide develop in foreland cover sequences between the interpreted basement faults. The interpreted basement faults project northward beneath the Himalayan system and southern Tibet. Our results suggest that several active Himalayan cross-strike faults, such as the ones related to many graben in southern Tibet, are rooted in the underplated Indian lower crust or step en échelon along interpreted basement faults. Our interpretation thus suggests that south Tibet graben are spatially related to deep-seated crustal-scale faults rooted in the underplated Indian crust. These major discontinuities partition the Himalayan range into distinct zones, and could ultimately contribute to lateral variability in tectonic evolution along the orogen's strike.

Description: We explore the impact of deep ductile shear zones on post-seismic deformation following a finite length strike-slip earthquake. We show that the pattern of post-seismic vertical surface deformation surrounding the fault is a discriminant for the existence of high viscosities immediately below the seismogenic layer, regardless of whether the model contains purely distributed creep or also includes a component of localized creep at subseismogenic depths. Post-seismic deformation characterized by initially fast relaxation followed by a slower relaxation is predicted by models that include both localized creep in a subseismogenic shear zone and distributed creep in the surrounding region, even if they only contain steady Maxwell viscoelasticity. This post-seismic deformation is similar to that in models that approximate the ductile lithosphere and/or asthenosphere with Burgers viscoelasticity. We find that the post-seismic deformation following the 1997 M w 7.6 Manyi, China, earthquake, is consistent with a post-seismic model composed of a lower Maxwell viscoelastic region with viscosity 10 19 Pa s and a 5 km wide, Maxwell viscoelastic shear zone with viscosity 10 18 Pa s beneath the fault.

Description: The Amundsen Sea Embayment of West Antarctica represents a key component in the tectonic history of Antarctic–New Zealand continental breakup. The region played a major role in the plate-kinematic development of the southern Pacific from the inferred collision of the Hikurangi Plateau with the Gondwana subduction margin at approximately 110–100 Ma to the evolution of the West Antarctic Rift System. However, little is known about the crustal architecture and the tectonic processes creating the embayment. During two ‘RV Polarstern’ expeditions in 2006 and 2010 a large geophysical data set was collected consisting of seismic-refraction and reflection data, ship-borne gravity and helicopter-borne magnetic measurements. Two P -wave velocity–depth models based on forward traveltime modelling of nine ocean bottom hydrophone recordings provide an insight into the lithospheric structure beneath the Amundsen Sea Embayment. Seismic-reflection data image the sedimentary architecture and the top-of-basement. The seismic data provide constraints for 2-D gravity modelling, which supports and complements P -wave modelling. Our final model shows 10–14-km-thick stretched continental crust at the continental rise that thickens to as much as 28 km beneath the inner shelf. The homogenous crustal architecture of the continental rise, including horst and graben structures are interpreted as indicating that wide-mode rifting affected the entire region. We observe a high-velocity layer of variable thickness beneath the margin and related it, contrary to other ‘normal volcanic type margins’, to a proposed magma flow along the base of the crust from beneath eastern Marie Byrd Land—West Antarctica to the Marie Byrd Seamount province. Furthermore, we discuss the possibility of upper mantle serpentinization by seawater penetration at the Marie Byrd Seamount province. Hints of seaward-dipping reflectors indicate some degree of volcanism in the area after break-up. A set of gravity anomaly data indicate several phases of fully developed and failed rift systems, including a possible branch of the West Antarctic Rift System in the Amundsen Sea Embayment.

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 estimate Eurasia-North America Plate motion rotations at ~1-Myr intervals for the past 20 Myr from more than 11 000 crossings of 21 magnetic reversals from Chron 1n (0.78 Ma) to C6no (19.72 Ma) and flow lines digitized from the Charlie Gibbs, Bight and Molloy fracture zones and transform faults. Adjusted for outward displacement, the 21 best-fitting rotations determined from a simultaneous inversion of the numerous kinematic data reconstruct the reversal crossings with weighted root mean square misfits of only 1–2 km and 0.2–7 km for the transform fault and fracture zone crossings. The new rotations clearly define a ~1000 km southward shift of the rotation pole and 20 per cent slowdown in seafloor spreading rates between 7 and 6 Ma, preceded by apparently steady plate motion from 19.7 to ~7 Ma. Data for times since C3An.2 (6.7 Ma) are well fit by a stationary pole of rotation and constant rate of angular opening, consistent with steady motion since 6.7 Ma. The southward shift of the rotation pole at 7–6 Ma implies that Eurasia-North America motion in northeastern Asia changed from slowly convergent before 7 Ma to slowly divergent afterward. Crossings of magnetic reversals C1n through C3An.1 (6.0 Ma) are well fit everywhere in the Arctic basin and south to the Azores triple junction, indicating that the Eurasia and North America plates have not deformed along their mutual boundary since at least 6.0 Ma. However, the new rotations systematically overrotate magnetic lineations older than C3An.1 (6.0 Ma) within 200 km of the Azores triple junction and also overrotate lineations older than C5n along the Gakkel Ridge in the Arctic Basin. Barring misidentifications of the magnetic anomalies in those areas, the pattern and magnitude of the systematic misfits imply that slow (~1 mm yr –1 ) distributed or microplate deformation occurred in one or both regions.

Description: Here we inverted the GPS data to infer the coseismic slip of the Tohoku-Oki earthquake and the time-dependent afterslip distribution in the 4 months following the main shock. The Tohoku-Oki earthquake showed an unexpected magnitude and a characteristic depth-dependent differentiation of seismic energy radiation. In this context the estimation and comparison of the distribution of the fault portions that slip coseismically and post-seismically contribute to a better understanding of the variation of frictional characteristics of the plate interface. The inferred coseismic slip extends in a relatively compact region located updip from the hypocentre and reaches its highest value (about 60 m) near the trench. Afterslip occurs mostly outside the coseismic rupture and is distributed in two main modal centres. It reaches its largest values in an area located downdip of the coseismic slip and extends to a depth of 80 km. In the depth range between 30 and 50 km afterslip overlaps the portion of the fault that experienced historical moderate earthquakes, high-frequency seismic radiation and thrust-type aftershocks. The behaviour of this area can be explained by a rheologically heterogeneous region made of a ductile fault matrix interspersed with compact brittle asperities. On the contrary, the region beneath 50–60 km depth is probably characterized by a fully velocity strengthening behaviour. Southern afterslip, located off-Chiba Prefecture, is probably related to the M w 7.9 Ibaraki-Oki aftershock. The northward extension of the afterslip stops at a latitude of about 40°N, just south of the off-Aomori region. This may be related to three large events occurred in this area during the last century and the consequent strong coupling or complete depletion of the accumulated strain that characterize this region.

Description: We present a revised interpretation of magnetic anomalies and fracture zones on the Southwest Indian Ridge (SWIR; Africa–Antarctica) and the Southeast Indian Ridge (SEIR; Capricorn–Antarctica) and use them to calculate 2-plate finite rotations for anomalies 34 to 20 (84 to 43 Ma). Central Indian Ridge (CIR; Capricorn–Africa) rotations are calculated by summing the SWIR and SEIR rotations. These rotations provide a high-resolution record of changes in the motion of India and Africa at the time of the onset of the Reunion plume head. An analysis of the relative velocities of India, Africa and Antarctica leads to a refinement of previous observations that the speedup of India relative to the mantle was accompanied by a slowdown of Africa. The most rapid slowdown of Africa occurs around Chron 32Ay (71 Ma), the time when India's motion relative to Africa notably starts to accelerate. Using the most recent Geomagnetic Polarity Timescale (GTS12) we show that India's velocity relative to Africa was characterized by an acceleration from roughly 60 to 180 mm yr –1 between 71 and 66 Ma, a short pulse of superfast motion (~180 mm yr –1 ) between 66 and 63 Ma, an abrupt slowdown to 120 mm yr –1 between 63 and 62 Ma, and then a long period (63 to 47 Ma) of gradual slowing, but still fast motion (~100 mm yr –1 ), which ends with a rapid slowdown after Chron 21o (47 Ma). Changes in the velocities of Africa and India with respect to the mantle follow a similar pattern. The fastest motion of India relative to the mantle, ~220 mm yr –1 , occurs during Chron 29R. The SWIR rotations constrain three significant changes in the migration path of the Africa–Antarctic stage poles: following Chron 33y (73 Ma), following Chron 31y (68 Ma), and following Chron 24o (54 Ma). The change in the migration path of the SWIR stage poles following Chron 33y is coincident with the most rapid slowdown in Africa's motion. The change in the migration path after Chron 31y, although coincident with the most rapid acceleration of India's northward motion, may be related to changes in ridge push forces on the SWIR associated with the onset of extension along the Bain transform fault zone. The initial slowdown in India's motion relative to Africa between 63 and 62 Ma is more abrupt than predictions based on published plume head force models, suggesting it might have been caused by a change in plate boundary forces. The abrupt change in the migration path of the SWIR stage poles after Chron 24o is not associated with major changes in the velocities of either Africa or India and may reflect Atlantic basin plate motion changes associated with the arrival at the Earth's surface of the Iceland plume head. The abruptness of India's slowdown after Chron 21o is consistent with a collision event.

Description: Thin plate flexure theory provides an accurate model for the response of the lithosphere to vertical loads on horizontal length scales ranging from tens to hundreds of kilometres. Examples include flexure at seamounts, fracture zones, sedimentary basins and subduction zones. When applying this theory to real world situations, most studies assume a locally uniform plate thickness to enable simple Fourier transform solutions. However, in cases where the amplitude of the flexure is prominent, such as subduction zones, or there are rapid variations in seafloor age, such as fracture zones, these models are inadequate. Here we present a computationally efficient algorithm for solving the thin plate flexure equation for non-uniform plate thickness and arbitrary vertical load. The iterative scheme takes advantage of the 2-D fast Fourier transform to perform calculations in both the spatial and spectral domains, resulting in an accurate and computationally efficient solution. We illustrate the accuracy of the method through comparisons with known analytic solutions. Finally, we present results from three simple models demonstrating the differences in trench outer rise flexure when 2-D variations in plate rigidity and applied bending moment are taken into account. Although we focus our analysis on ocean trench flexure, the method is applicable to other 2-D flexure problems having spatial rigidity variations such as seamount loading of a thermally eroded lithosphere or flexure across the continental–oceanic crust boundary.

Description: Theory has been long established for computing the elastic response of a spherically symmetric terrestrial planetary body to both body tide and surface loading forces. However, for a planet with laterally heterogeneous mantle structure, the response is usually computed using a fully numerical approach. In this paper, we develop a semi-analytic method based on perturbation theory to solve for the elastic response of a planetary body with lateral heterogeneities in its mantle. We present a derivation of the governing equations for our second-order perturbation method and use them to study the high-order tidal effects caused by mode coupling between degree-2 body tide forcing and the laterally heterogeneous elastic structure of the mantle. We test our method by applying it to the Moon in which small long-wavelength lateral heterogeneities are assumed to exist in the elastic moduli of the lunar mantle. The tidal response of the Moon is determined mode by mode, for lateral heterogeneities with different depth ranges within the mantle and different horizontal scales. Our perturbation method solutions are compared with numerical results, showing remarkable agreement between the two methods. We conclude that our perturbation method provides accurate results and can be adapted to address a variety of forward and inverse response problems.

Description: Teleseismic data recorded at 13 broad-band seismological stations across northwest part of the Tethyan Himalaya and eastern Ladakh are analysed to determine the seismic characteristics of the crust and upper mantle beneath the northwest India–Asia collision zone. The receiver functions computed from teleseismic P- waveform for a wide range of backazimuth show strong azimuthal variation in the Indus suture zone (ISZ), the zone which marks the collision and subsequent subduction of both the Tethyan oceanic plate and Indian continental plate beneath Eurasia. The teleseismic waves piercing the ISZ do not show clear P -to- S ( Ps ) converted phase at the depth of Moho. In contrast, the waves piercing the Karakoram zone, Ladakh batholith and the Tethyan Himalayan region south of the ISZ clearly show the Moho converted Ps phase and corresponding inverted models reveal variation of crustal thickness from ~60 km beneath the Tethyan Himalaya to ~80 km beneath the Karakoram fault zone. A prominent intracrustal low velocity zone (IC-LVZ) is detected in the shear wave velocity models within the depth range ~15–40 km. The IC-LVZ identified at the stations both north and south of the ISZ can be interpreted as due to presence of fluid/partial melt. Our study provides compelling evidence that the mid-crustal low velocity zone does extend across the suture zone, in to the Tethyan Himalaya. The contact between this serpentinized ultramafic rocks and the eclogitized Indian continental crust in the suture zone is identified at ~47–50 km depth.

Description: The Tien Shan is the largest active intracontinental orogenic belt on Earth. To better understand the processes causing mountains to form at great distances from a plate boundary, we analyse passive source seismic data collected on 40 broad-band stations of the MANAS project (2005–2007) and 12 stations of the permanent KRNET seismic network to determine variations in crustal thickness and shear wave speed across the range. We jointly invert P - and S -wave receiver functions with surface wave observations from both earthquakes and ambient noise to reduce the ambiguity inherent in the images obtained from the techniques applied individually. Inclusion of ambient noise data improves constraints on the upper crust by allowing dispersion measurements to be made at shorter periods. Joint inversion can also reduce the ambiguity in interpretation by revealing the extent to which various features in the receiver functions are amplified or eliminated by interference from multiples. The resulting wave speed model shows a variation in crustal thickness across the range. We find that crustal velocities extend to ~75 km beneath the Kokshaal Range, which we attribute to underthrusting of the Tarim Basin beneath the southern Tien Shan. This result supports the plate model of intracontinental convergence. Crustal thickness elsewhere beneath the range is about 50 km, including beneath the Naryn Valley in the central Tien Shan where previous studies reported a shallow Moho. This difference apparently is the result of wave speed variations in the upper crust that were not previously taken into account. Finally, a high velocity lid appears in the upper mantle of the Central and Northern part of the Tien Shan, which we interpret as a remnant of material that may have delaminated elsewhere under the range.

Description: Errors in the satellite orbits are considered to be a limitation for Interferometric Synthetic Aperture Radar (InSAR) time-series techniques to accurately measure long-wavelength (〉50 km) ground displacements. Here we examine how orbital errors propagate into relative InSAR line-of-sight velocity fields and evaluate the contribution of orbital errors to the InSAR uncertainty. We express the InSAR uncertainty due to the orbital errors in terms of the standard deviations of the velocity gradients in range and azimuth directions (range and azimuth uncertainties). The range uncertainty depends on the magnitude of the orbital errors, the number and time span of acquisitions. Using reported orbital uncertainties we find range uncertainties of less than 1.5 mm yr –1  100 km –1 for ERS, less than 0.5 mm yr –1  100 km –1 for Envisat and ~0.2 mm yr –1  100 km –1 for TerraSAR-X and Sentinel-1. Under a conservative scenario, we find azimuth uncertainties of better than 1.5 mm yr –1  100 km –1 for older satellites (ERS and Envisat) and better than 0.5 mm yr –1  100 km –1 for modern satellites (TerraSAR-X and Sentinel-1). We validate the expected uncertainties using LOS velocity fields obtained from Envisat SAR imagery. We find residual gradients of 0.8 mm yr –1  100 km –1 or less in range and of 0.95 mm yr –1  100 km –1 or less in azimuth direction, which fall within the 1 to 2 uncertainties. The InSAR uncertainties due to the orbital errors are significantly smaller than generally expected. This shows the potential of InSAR systems to constrain long-wavelength geodynamic processes, such as continent-scale deformation across entire plate boundary zones.

Description: A method for subsurface recognition of blind geological bodies is presented using combined surface constraints and 3-D structural modelling that incorporates constraints from detailed mapping, and potential-field inversion modelling. This method is applied to the Mount Painter Province and demonstrates that addition of low density material is required to reconcile the gravity signature of the region. This method may be an effective way to construct 3-D models in regions of excellent structural control, and can be used to assess the validity of surface structures with 3-D architecture. Combined geological and potential-field constrained inversion modelling of the Mount Painter Province was conducted to assess the validity of the geological models of the region. Magnetic susceptibility constrained stochastic property inversions indicates that the northeast to southwest structural trend of the relatively magnetic meta-sedimentary rocks of the Radium Creek Group in the Mount Painter Inlier is reconcilable with the similar, northeast to southwest trending positive magnetic anomalies in the region. Radium Creek Group packages are the major contributor of the total magnetic response of the region. However field mapping and the results of initial density constrained stochastic property inversion modelling do not correlate with a large residual negative gravity anomaly central to the region. Further density constrained inversion modelling indicates that an additional large body of relatively low density material is needed within the model space to account for this negative density anomaly. Through sensitivity analysis of multiple geometrical and varied potential-field property inversions, the best-fitting model records a reduction in gravity rms misfit from 21.9 to 1.69 mGal, representing a reduction from 56 to 4.5 per cent in respect to the total dynamic range of 37.5 mGal of the residual anomaly. This best-fitting model incorporates a volumetrically significant source body of interpreted felsic, low density material (10 12 m 3 ) impinging on the central-west of the Mount Painter Inlier and overlying Neoproterozoic sequences, and the emplacement of more mafic affinities in the northeast and east. The spatial association and circular geometry of these granitoid bodies suggests an affinity with the Palaeozoic ~460–440 Ma British Empire Granite that outcrops in the Mount Painter Inlier. The intrusion of this additional material in the Palaeozoic could either be the product of; or contributed to, an increased local geotherm and heat flow in the region during the Palaeozoic.

Description: The Pacific and Australian plates in the South Island, New Zealand (NZ) converge at a rate of about 4 cm yr –1 . Accommodation of the continental part of this convergence in the lithospheric mantle is both poorly understood and currently controversial yet it is a problem of fundamental importance for understanding lithospheric thickening. End-member possibilities range from the classical model of asymmetric subduction to symmetric viscous thickening. Seismic tomography has the potential to image this process. However, tomographic images to date are poorly constrained due to the lack of appropriate earthquakes. Improved teleseismic tomography of the region has been achieved by increasing data coverage and applying a novel scheme of correcting for crustal structure by ray tracing through a newly created model of shallow shear wave velocity derived from the inversion of noise-based dispersion measurements. Our resulting models suggest the lithospheric mantle high velocities at the continental plate boundary extend no deeper than approximately 125 km, evidence against both previous models of viscous drip and typical asymmetric subduction zones. This high velocity core extends from north to south along the axis of South Island suggesting that mantle convergence is accommodated along the older, mid-Cenozoic, plate boundary. West of South Island, a high velocity west dipping zone may define the remnant Cretaceous subduction zone that has been distorted by Cenozoic transcurrent deformation. We present our new 3-D seismic velocity models together with a compatible tectonic model and discuss their implications for the nature of lithospheric evolution at this convergent boundary.

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.

Description: We have derived a shallow subsurface 2-D tomographic P -wave velocity image of the Deccan Volcanic Province (DVP) of India using first-arrival traveltime data along a 90-km-long N–S trending seismic profile in the Deccan Syneclise region. The tomographic image depicts smooth velocity variations of Quaternary and Tertiary (2.0–3.0 km s –1 ) sediments, basalts/traps (5.0–5.5 km s –1 ), sub-trappean Mesozoic sediments (4.3–4.5 km s –1 ) as well as the basement (5.9–6.1 km s –1 ) geometry down to a maximum depth of 5.0 km. Due to Late Cretaceous volcanism and outpouring of basaltic lava flows, this region is affected by numerous dyke intrusions and thick basaltic trap (2–3 km) exposed on the surface and surrounded by graben structures due to deep basinal faults forming a large igneous province. Although sub-basalt imaging is a major challenge for the oil industry, with the help of tomographic imaging technique of first-arrival seismic refraction data, we were able to image sub-trappean Mesozoic sediments (〈0.75 km) deposited below the two sequences of thick basaltic flows above the basement. The imaged Mesozoic sediments are expected to contain hydrocarbon because of their wide extension in this sedimentary basin with suitable trapping mechanism due to basalts. The robustness of the velocity image is assessed through numerous tests like velocity perturbations, 2 estimates, rms residuals of traveltime fit, uncertainty estimates through computation of ray-density or hits and series of checkerboard resolution tests with velocity anomalies having different cell size. The thickness of the basalt and the sub-trappean Mesozoic sediments along with the basement geometry obtained from tomography are constrained through ray-trace modelling and pre-stack depth migration (PSDM) of the wide-angle reflection phases for different shot gathers along the profile.

Description: Numerical experiments of convection with grain-damage are used to develop scaling laws for convective heat flow, mantle velocity and plate velocity across the stagnant lid and plate-tectonic regimes. Three main cases are presented in order of increasing complexity: a simple case wherein viscosity is only dependent on grain size, a case where viscosity depends on temperature and grain size, and finally a case where viscosity is temperature and grain size sensitive, and the grain-growth (or healing) is also temperature sensitive. In all cases, convection with grain-damage scales differently than Newtonian convection; whereas the Nusselt number (Nu), typically scales with the reference Rayleigh number, Ra 0 , to the 1/3 power, for grain-damage this exponent is larger because increasing Ra 0 also enhances damage. In addition, Nu, mantle velocity, and plate velocity are also functions of the damage to healing ratio, ( D / H ); increasing D / H increases Nu because more damage leads to more vigorous convection. For the fully realistic case, numerical results show stagnant lid convection, fully mobilized convection that resembles the temperature-independent viscosity case, and partially mobile or transitional convection, depending on D / H , Ra 0 , and the activation energies for viscosity and healing. Applying our scaling laws for the fully realistic case to Earth and Venus we demonstrate that increasing surface temperature dramatically decreases plate speed and heat flow, essentially shutting down plate tectonics, due to increased healing in lithospheric shear zones, as proposed previously. Contrary to many previous studies, the transitional regime between the stagnant lid and fully mobilized regimes is large, and the transition from stagnant lid to mobile convection is gradual and continuous. Thus planets could exhibit a full range of surface mobility, as opposed to the bimodal distribution of fully mobile lid planets and stagnant lid planets that is typically assumed.

Description: Three magnetotelluric (MT) profiles in northwestern Canada cross the central and western segments of Great Slave Lake shear zone (GSLsz), a continental scale strike-slip structure active during the Slave-Rae collision in the Proterozoic. Dimensionality analysis indicates that (i) the resistivity structure is approximately 2-D with a geoelectric strike direction close to the dominant geological strike of N45°E and that (ii) electrical anisotropy may be present in the crust beneath the two southernmost profiles. Isotropic and anisotropic 2-D inversion and isotropic 3-D inversions show different resistivity structures on different segments of the shear zone. The GSLsz is imaged as a high resistivity zone (〉5000  m) that is at least 20 km wide and extends to a depth of at least 50 km on the northern profile. On the southern two profiles, the resistive zone is confined to the upper crust and pierces an east-dipping crustal conductor. Inversions show that this dipping conductor may be anisotropic, likely caused by conductive materials filling a network of fractures with a preferred spatial orientation. These conductive regions would have been disrupted by strike-slip, ductile deformation on the GSLsz that formed granulite to greenschist facies mylonite belts. The pre-dominantly granulite facies mylonites are resistive and explain why the GSLsz appears as a resistive structure piercing the east-dipping anisotropic layer. The absence of a dipping anisotropic/conductive layer on the northern MT profile, located on the central segment of the GSLsz, is consistent with the lack of subduction at this location as predicted by geological and tectonic models.

Description: The thickness and equivalent global sea level contribution of an improved model of the central and northern Laurentide Ice Sheet is constrained by 24 relative sea level histories and 18 present-day GPS-measured vertical land motion rates. The final model, termed Laur16, is derived from the ICE-5G model by holding the timing history constant and iteratively adjusting the thickness history, in four regions of northern Canada. In the final model, the last glacial maximum (LGM) thickness of the Laurentide Ice Sheet west of Hudson Bay was ~3.4–3.6 km. Conversely, east of Hudson Bay, peak ice thicknesses reached ~4 km. The ice model thicknesses inferred for these two regions represent, respectively, a ~30 per cent decrease and an average ~20–25 per cent increase to the load thickness relative to the ICE-5G reconstruction, which is generally consistent with other recent studies that have focussed on Laurentide Ice Sheet history. The final model also features peak ice thicknesses of 1.2–1.3 km in the Baffin Island region, a modest reduction relative to ICE-5G and unchanged thicknesses for a region in the central Canadian Arctic Archipelago west of Baffin Island. Vertical land motion predictions of the final model fit observed crustal uplift rates well, after an adjustment is made for the elastic crustal response to present-day ice mass changes of regional ice cover. The new Laur16 model provides more than a factor of two improvement of the fit to the RSL data ( 2 measure of misfit) and a factor of nine improvement to the fit of the GPS data (mean squared error measure of fit), compared to the ICE-5G starting model. Laur16 also fits the regional RSL data better by a factor of two and gives a slightly better fit to GPS uplift rates than the recent ICE-6G model. The volume history of the Laur16 reconstruction corresponds to an up to 8 m reduction in global sea level equivalent compared to ICE-5G at LGM.

Description: The 2011 October 23 M W 7.1 Van earthquake in eastern Turkey caused ~600 deaths and caused widespread damage and economic loss. The seismogenic rupture was restricted to 10–25 km in depth, but aseismic surface creep, coincident with outcrop fault exposures, was observed in the hours to months after the earthquake. We combine observations from radar interferometry, seismology, geomorphology and Quaternary dating to investigate the geological slip rate and seismotectonic context of the Van earthquake, and assess the implications for continuing seismic hazard in the region. Transient post-seismic slip on the upper Van fault started immediately following the earthquake, and decayed over a period of weeks; it may not fully account for our long-term surface slip-rate estimate of ≥0.5 mm yr –1 . Post-seismic slip on the Bostanici splay fault initiated several days to weeks after the main shock, and we infer that it may have followed the M W  5.9 aftershock on the 9th November. The Van earthquake shows that updip segmentation can be important in arresting seismic ruptures on dip-slip faults. Two large, shallow aftershocks show that the upper 10 km of crust can sustain significant earthquakes, and significant slip is observed to have reached the surface in the late Quaternary, so there may be a continuing seismic hazard from the upper Van fault and the associated splay. The wavelength of folding in the hanging wall of the Van fault is dominated by the structure in the upper 10 km of the crust, masking the effect of deeper seismogenic structures. Thus, models of subsurface faulting based solely on surface folding and faulting in regions of reverse faulting may underestimate the full depth extent of seismogenic structures in the region. In measuring the cumulative post-seismic offsets to anthropogenic structures, we show that Structure-from-Motion can be rapidly deployed to create snapshots of post-seismic displacement. We also demonstrate the utility of declassified Corona mission imagery (1960s–1970s) for geomorphic mapping in areas where recent urbanization has concealed the geomorphic markers.

Description: Post-seismic deformation is commonly attributed to viscoelastic relaxation and/or afterslip, although discerning between the two driving mechanisms can be difficult. A major complication in modeling post-seismic deformation is that forward models can be computationally expensive, making it difficult to adequately search model space to find the optimal fault slip distribution and lithospheric viscosity structure that can explain observable post-seismic deformation. We propose an inverse method which uses coseismic and early post-seismic deformation to rapidly and simultaneously estimate a fault slip history and an arbitrarily discretized viscosity structure of the lithosphere. Our method is based on an approximation which is applicable to the early post-seismic period and expresses surface deformation resulting from viscoelastic relaxation as a linearized function with respect to lithospheric fluidity. We demonstrate this approximation using two-dimensional earthquake models. We validate the approximation and our inverse method using two three-dimensional synthetic tests. The success of our synthetic tests suggests that our method is capable of distinguishing the mechanisms driving early post-seismic deformation and recovering an effective viscosity structure of the lithosphere.

Description: Within a slab experiencing present-day lateral break-off, a particular type of earthquakes is expected to cluster at the detachment horizon tip: namely, events generated by reverse faulting, with the approximately horizontal compression involved acting along the strike of the slab. Such a cluster of moderate magnitude earthquakes (4.7 ≤  m b ≤ 5.0) was identified in this study at the 160–175 km depth range of the Vrancea seismogenic body, in the Southeast Carpathians mountains collision environment. The corresponding cluster epicentres were systematically positioned at the boundary between a region being subject ( cf. published GPS records), to present-day upward movements, and another one that underwent present-day subsidence. Such an overall setting seems to suggest that a lateral break-off is currently developing at the indicated depth within the Vrancea slab, leading to topographic uplift above the already detached slab section, and to enhanced subsidence above the section to which the gravitational slab pull was being transferred. In addition, by taking into account some systematic time correspondence which we documented between moderate magnitude events of the 160–175 km depth cluster and the subsequent strong Vrancea shocks ( M w ≥ 6.9), it appears that the latter, although occurring at much shallower depths (roughly, in the 80–140 km range), were also controlled by the break-off progress.

Description: Bajada del Diablo is located in the Northern Patagonian Massif, Chubut, Argentina. The study area includes several circular structures found in Miocene olivine basalts of the Quiñelaf Eruptive Complex and in the Late Pliocene/Early Pleistocene Pampa Sastre conglomerates. An impact origin has been proposed for these circular structures. With the aim of further investigate the proposed impact origin, topographic, gravimetric, magnetic, resistivity, palaeomagnetic and electromagnetic surveys in two circular structures (‘8’ and ‘G’) located in Pampa Sastre conglomerates and in basalts of the Quiñelaf Eruptive Complex were carried out. The new geophysical results support the hypothesis of an impact origin. However, the confirmation of such an origin through the findings of shock metamorphism evidences and/or the recovery of meteorites remains elusive.

Description: Estimating in situ stress based on hydraulic fracturing data typically depends on interpretation of the breakdown, secondary breakdown (‘reopening’) and shut-in pressure. While it has been recognized that the near-wellbore stress field should be taken into account and that the compressibility of the injection system and the viscous flow of the fluid can diminish the accuracy of stress estimates, these issues have not been well quantified. A coupled numerical model that includes the compressibility of the injection system and the flow of a viscous fluid in a plane-strain hydraulic fracture extending from a wellbore, in an impermeable rock, and in the presence of a non-isotropic in situ stress field provides a basic tool for estimating the order of the error associated with hydraulic fracturing stress measurements under non-ideal conditions. The main findings of this work are model-based guidelines on the values of relevant dimensionless parameter groups to ensure sufficient accuracy of stress estimates that use idealized models. When these guidelines cannot be met under field conditions, the model can be further applied to obtain first-order corrections that account for compressibility, viscosity and near-wellbore effects.

Description: We consider the problem of seismic velocity change estimation using ambient noise recordings. Motivated by Zhan et al. , we study how the velocity change estimation is affected by seasonal fluctuations in the noise sources. More precisely, we consider a numerical model and introduce spatio-temporal seasonal fluctuations in the noise sources. We show that indeed, as pointed out by Zhan et al. , the stretching method is affected by these fluctuations and produces misleading apparent velocity variations which reduce dramatically the signal to noise ratio of the method. We also show that these apparent velocity variations can be eliminated by an adequate normalization of the cross-correlation functions. Theoretically we expect our approach to work as long as the seasonal fluctuations in the noise sources are uniform, an assumption which holds for closely located seismic stations. We illustrate with numerical simulations in homogeneous and scattering media that the proposed normalization significantly improves the accuracy of the velocity change estimation. Similar behaviour is also observed with real data recorded in the Aegean volcanic arc. We study in particular the volcano of Santorini during the seismic unrest of 2011–2012 and observe a decrease in the velocity of seismic waves which is correlated with GPS measured elevation.

Description: Passive seismic experiment was carried out at the SW contact of the Dinarides and Pannonian basin to determine the crustal structure and velocity discontinuities. The aim of the experiment was to define the relationship between the Adriatic microplate and the Pannonian segment as a part of the European plate. Most of the temporary seismic stations were deployed in Croatia along the Alp07 profile—a part of the active-source ALP 2002 project. About 300-km-long profile stretches from Istra peninsula to the Drava river, in a WSW–ESE direction. Teleseismic events recorded on 13 temporary seismic stations along the profile were analysed by P -receiver function method. Two types of characteristic receiver functions (RF) have been identified, belonging to Dinaridic and Pannonian crusts as defined on the Alp07 profile, while in transitional zone there are both types. Three major crustal discontinuities can be identified for the Dinaridic type: sedimentary basement, intracrustal discontinuity and Mohorovičić discontinuity, whereas the Pannonian type revealed only two discontinuities. The intracrustal discontinuity was not observed in the Pannonian type, thus pointing to a single-layered crust in the Pannonian basin. Two interpretation methods were applied: forward modelling of the receiver functions and H – stacking method, and the results were compared with the active-source seismic data at deep refraction profile Alp07. The receiver function modelling has given reliable results of the Moho depths that are in accordance with the seismic refraction results at the end of the Alp07 profile, that is in the area of Pannonian crust characterized by simple crustal structure and low seismic velocities ( Vp between 5.9 and 6.2 km s –1 ). In the Dinarides and its peripheral parts, receiver function modelling regularly gives greater Moho depths, up to +15 per cent, due to more complex crustal structure. The depths of the Moho calculated by the H – stacking method vary within wide limits (±13 km), due to band limited data of short-period stations. The results at five stations have to be rejected because of huge deviations in comparison with all previous results, while at the other seven stations the Moho depths vary within ±15 per cent around the Moho discontinuity of the Alp07 profile.

Description: Megathrust earthquakes of magnitude close to 9 are followed by large-scale (thousands of km) and long-lasting (decades), significant crustal and mantle deformation. This deformation can be observed at the surface and quantified with GPS measurements. Here we report on deformation observed during the 5 yr time span after the 2010 M w 8.8 Maule Megathrust Earthquake (2010 February 27) over the whole South American continent. With the first 2 yr of those data, we use finite element modelling (FEM) to relate this deformation to slip on the plate interface and relaxation in the mantle, using a realistic layered Earth model and Burgers rheologies. Slip alone on the interface, even up to large depths, is unable to provide a satisfactory fit simultaneously to horizontal and vertical displacements. The horizontal deformation pattern requires relaxation both in the asthenosphere and in a low-viscosity channel along the deepest part of the plate interface and no additional low-viscosity wedge is required by the data. The vertical velocity pattern (intense and quick uplift over the Cordillera) is well fitted only when the channel extends deeper than 100 km. Additionally, viscoelastic relaxation alone cannot explain the characteristics and amplitude of displacements over the first 200 km from the trench and aseismic slip on the fault plane is needed. This aseismic slip on the interface generates stresses, which induce additional relaxation in the mantle. In the final model, all three components (relaxation due to the coseismic slip, aseismic slip on the fault plane and relaxation due to aseismic slip) are taken into account. Our best-fit model uses slip at shallow depths on the subduction interface decreasing as function of time and includes (i) an asthenosphere extending down to 200 km, with a steady-state Maxwell viscosity of 4.75 x 10 18 Pa s; and (ii) a low-viscosity channel along the plate interface extending from depths of 55–135 km with viscosities below 10 18 Pa s.

Description: We provide an updated present-day stress map for the Italian territory. Following the World Stress Map (WSM) Project guidelines, we list the different stress indicators, explaining the criteria used to select data. We discuss the data, which will also be included in the 2016 release of the WSM, highlighting the areas for which we have added stress information. Our map displays the minimum horizontal stress orientations inferred from crustal stress indicators down to 40 km depth using data of A–C quality, updated for earthquakes until December 2015. We have completely reviewed all data, and the data set now contains 855 entries, in contrast to the previous 715. The number of data with A–C quality of 630 corresponds to an increase of 26 per cent relative to the previous data set. In particular, the new data set contains the results of the analysis of borehole breakouts, critically reviewed data from earthquake focal mechanisms, data concerning active faults, formal inversions of focal mechanisms of seismic sequences or of restricted areas and one stress determination from overcoring. The new data set defines the stress field in areas not well covered by the previous data: the region north to the Po Plain and the central Adriatic sea, both characterized by a thrust- and strike-faulting regime, the northern Sicilian belt with a prevailing normal-faulting regime, and the Ionian sea with a strike-slip regime.

Description: Measurements of the velocity field associated with plumes rising through a viscous fluid are performed using stereoscopic Particle-Image Velocimetry in the Rayleigh number range 4.4 x 10 5 –6.4 x 10 5 . The experimental model is analogous to a mantle plume rising from the core–mantle boundary to the base of the lithosphere. The behaviour of the plume is studied throughout its life cycle, which is broken up into four stages; the Formation Stage, when the plume forms; the Rising Stage, when the plume rises through the fluid; the Spreading Stage, when the plume reaches the surface and spreads; and finally the Declining Stage, when the heat source has been removed and the plume weakens. The latter three stages are examined in terms of the Finite-Time Lyapunov Exponent fields and the advection of passive tracers throughout the flow. The temperature at the heater and near the fluid surface are measured using thermocouples to infer how the presence of a mantle plume would produce excess temperature near the lithosphere throughout the various stages of its life cycle. In all experiments, a time lag is observed between the removal of the heat source and the decline in the excess temperature near the surface, which is proportional to the rise time. A simple analytical model is presented, which suggests that under mantle conditions (i.e. negligible thermal diffusion), the relationship between the time lag and the rise time is robust and independent of the Rayleigh number; however, the constant of proportionality is closer to unity in the absence of diffusion. Once the heat source is removed, the excess temperature near the surface declines exponentially at a rate that is inversely proportional to the rise time. The implications of this result are discussed in terms of the decline in volcanism in the Louisville hotspot chain over the past 20 Ma. The rise velocity of material in the plume is examined; the rise velocity is found to vary significantly with the plume height in a manner that is inconsistent with many of the common semi-analytical models of thermal plumes in the literature. It is also argued that this height dependency will cause estimates of the rise velocity based on the decay series of uranium isotopes to significantly underestimate the true value.

Description: The Makran subduction zone is one of the last convergent margins to be investigated using space-based geodesy. While there is a lack of historical and modern instrumentation in the region, a sparse sampling of continuous and campaign measurements over the past decade has allowed us to make the first estimates of convergence rates. We combine GPS measurements from 20 stations located in Iran, Pakistan and Oman along with hypocentral locations from the International Seismological Centre to create a preliminary 3-D estimate of the geometry of the megathrust, along with a preliminary fault-coupling model for the Makran subduction zone. Using a convergence rate which is strongly constrained by measurements from the incoming Arabia plate along with the backslip method of Savage, we find the Makran subduction zone appears to be locked to a depth of at least 38 km and accumulating strain.We also find evidence for a segmentation of plate coupling, with a 300 km long section of reduced plate coupling. The range of acceptable locking depths from our modelling and the 900 km along-strike length for the megathrust, makes the Makran subduction zone capable of earthquakes up to M w  = 8.8. In addition, we find evidence for slow-slip-like transient deformation events on two GPS stations. These observations are suggestive of transient deformation events observed in Cascadia, Japan and elsewhere.

Description: In this study, we present a new synthesis of GPS velocities for tectonic deformation within the Tibetan Plateau and its surrounding areas, a combined data set of ~1854 GPS-derived horizontal velocity vectors. Assuming that crustal deformation is localized along major faults, a block modelling approach is employed to interpret the GPS velocity field. We construct a 30-element block model to describe present-day deformation in western China, with half of them located within the Tibetan Plateau, and the remainder located in its surrounding areas. We model the GPS velocities simultaneously for the effects of block rotations and elastic strain induced by the bounding faults. Our model yields a good fit to the GPS data with a mean residual of 1.08 mm a –1 compared to the mean uncertainty of 1.36 mm a –1 for each velocity component, indicating a good agreement between the predicted and observed velocities. The major strike-slip faults such as the Altyn Tagh, Xianshuihe, Kunlun and Haiyuan faults have relatively uniform slip rates in a range of 5–12 mm a –1 along most of their segments, and the estimated fault slip rates agree well with previous geologic and geodetic results. Blocks having significant residuals are located at the southern and southeastern Tibetan Plateau, suggesting complex tectonic settings and further refinement of accurate definition of block geometry in these regions.

Description: We present numerical models of mantle dynamics forced by plate velocities history in the last 450 Ma. The lower-mantle rheology and the thickness of a dense basal layer are systematically varied and several initial procedures are considered for each case. For some cases, the dependence on the mantle convection vigour is also examined. The resulting evolution of the CMB heat flux is analysed in terms of criteria to promote or inhibit reversals inferred from numerical dynamos. Most models present a rather dynamic lower mantle with the emergence of two thermochemical piles towards present-day. Only a small minority of models present two stationary piles over the last 450 Myr. At present-day, the composition field obtained in our models is found to correlate better with tomography than the temperature field. In addition, the temperature field immediately at the CMB (and thus the heat flux pattern) slightly differs from the average temperature field over the 100-km thick mantle layer above it. The evolution of the mean CMB heat flux or of the amplitude of heterogeneity seldom presents the expected correlation with the evolution of the palaeomagnetic reversal frequency suggesting these effects cannot explain the observations. In contrast, our analysis favours ‘inertial control’ on the geodynamo associated with polar cooling and in some cases break of Taylor columns in the outer core as sources of increased reversal frequency. Overall, the most likely candidates among our mantle dynamics models involve a viscosity increase in the mantle equal or smaller than 30: models with a discontinuous viscosity increase at the transition zone tend to agree better at present-day with observations of seismic tomography, but models with a gradual viscosity increase agree better with some of the criteria proposed to affect reversal frequency.

Description: Lateral viscosity variations (LVVs) in the mantle influence geodynamic processes and their surface expressions. With the observed long-wavelength geoid, free-air anomaly, gravity gradient in three directions and discrete, high-accuracy residual topography, we invert for depth- and temperature-dependent and tectonically regionalized mantle viscosity with a mantle flow model. The inversions suggest that long-wavelength gravitational and topographic signals are mainly controlled by the radial viscosity profile; the pre-Cambrian lithosphere viscosity is slightly (~ one order of magnitude) higher than that of oceanic and Phanerozoic lithosphere; plate margins are substantially weaker than plate interiors; and viscosity has only a weak apparent, dependence on temperature, suggesting either a balancing between factors or a smoothing of actual higher amplitude, but short wavelength, LVVs. The predicted large-scale lithospheric stress regime (compression or extension) is consistent with the world stress map (thrust or normal faulting). Both recent compiled high-accuracy residual topography and the predicted dynamic topography yield ~1 km amplitude long-wavelength dynamic topography, inconsistent with recent studies suggesting amplitudes of ~100 to ~500 m. Such studies use a constant, positive admittance (transfer function between topography and gravity), in contrast to the evidence which shows that the earth has a spatially and wavelength-dependent admittance, with large, negative admittances between ~4000 and ~10 4 km wavelengths.

Description: We present a new, computationally efficient numerical method to simulate global seismic wave propagation in realistic 3-D Earth models. We characterize the azimuthal dependence of 3-D wavefields in terms of Fourier series, such that the 3-D equations of motion reduce to an algebraic system of coupled 2-D meridian equations, which is then solved by a 2-D spectral element method (SEM). Computational efficiency of such a hybrid method stems from lateral smoothness of 3-D Earth models and axial singularity of seismic point sources, which jointly confine the Fourier modes of wavefields to a few lower orders. We show novel benchmarks for global wave solutions in 3-D structures between our method and an independent, fully discretized 3-D SEM with remarkable agreement. Performance comparisons are carried out on three state-of-the-art tomography models, with seismic period ranging from 34 s down to 11 s. It turns out that our method has run up to two orders of magnitude faster than the 3-D SEM, featured by a computational advantage expanding with seismic frequency.

Description: We report on a study to explore the deep electrical conductivity structure of the Dead Sea Basin (DSB) using magnetotelluric (MT) data collected along a transect across the DSB where the left lateral strike-slip Dead Sea transform (DST) fault splits into two fault strands forming one of the largest pull-apart basins of the world. A very pronounced feature of our 2-D inversion model is a deep, subvertical conductive zone beneath the DSB. The conductor extends through the entire crust and is sandwiched between highly resistive structures associated with Precambrian rocks of the basin flanks. The high electrical conductivity could be attributed to fluids released by dehydration of the uppermost mantle beneath the DSB, possibly in combination with fluids released by mid- to low-grade metamorphism in the lower crust and generation of hydrous minerals in the middle crust through retrograde metamorphism. Similar high conductivity zones associated with fluids have been reported from other large fault systems. The presence of fluids and hydrous minerals in the middle and lower crust could explain the required low friction coefficient of the DST along the eastern boundary of the DSB and the high subsidence rate of basin sediments. 3-D inversion models confirm the existence of a subvertical high conductivity structure underneath the DSB but its expression is far less pronounced. Instead, the 3-D inversion model suggests a deepening of the conductive DSB sediments off-profile towards the south, reaching a maximum depth of approximately 12 km, which is consistent with other geophysical observations. At shallower levels, the 3-D inversion model reveals salt diapirism as an upwelling of highly resistive structures, localized underneath the Al-Lisan Peninsula. The 3-D model furthermore contains an E–W elongated conductive structure to the northeast of the DSB. More MT data with better spatial coverage are required, however, to fully constrain the robustness of the above-mentioned off-profile features.

Description: Estimating how topography is maintained provides insights into the different factors responsible for surface deformations and their relative roles. Here, we develop a new and simple approach to assess the degree of isostatic compensation of continental topography at regional scale from GOCE gravity gradients. We calculate the ratio between the radial gradient observed by GOCE and that calculated from topography only. From analytical and statistical formulations, simple relationships between this ratio and the degree of compensation are obtained under the Airy–Heiskanen isostasy hypothesis. Then, a value of degree of compensation at each point of study area can be easily deduced. We apply our method to the Alaska-Canada Cordillera and validate our results by comparison with a standard isostatic gravity anomaly model and additional geophysical information for this area. Both our GOCE-based results and the isostatic anomaly show that Airy–Heiskanen isostasy prevails for the Yukon Plateau whereas additional mechanisms are required to support topography below the Northwest Territories Craton and the Yakutat collision zone.

Description: Bayesian sampling based inversions require many thousands or even millions of forward models, depending on how nonlinear or non-unique the inverse problem is, and how many unknowns are involved. The result of such a probabilistic inversion is not a single ‘best-fit’ model, but rather a probability distribution that is represented by the entire model ensemble. Often, a geophysical inverse problem is non-unique, and the corresponding posterior distribution is multimodal, meaning that the distribution consists of clusters with similar models that represent the observations equally well. In these cases, we would like to visualize the characteristic model properties within each of these clusters of models. However, even for a moderate number of inversion parameters, a manual appraisal for a large number of models is not feasible. This poses the question whether it is possible to extract end-member models that represent each of the best-fit regions including their uncertainties. Here, I show how a machine learning tool can be used to characterize end-member models, including their uncertainties, from a complete model ensemble that represents a posterior probability distribution. The model ensemble used here results from a nonlinear geodynamic inverse problem, where rheological properties of the lithosphere are constrained from multiple geophysical observations. It is demonstrated that by taking vertical cross-sections through the effective viscosity structure of each of the models, the entire model ensemble can be classified into four end-member model categories that have a similar effective viscosity structure. These classification results are helpful to explore the non-uniqueness of the inverse problem and can be used to compute representative data fits for each of the end-member models. Conversely, these insights also reveal how new observational constraints could reduce the non-uniqueness. The method is not limited to geodynamic applications and a generalized MATLAB code is provided to perform the appraisal analysis.

Description: Inference of the mantle viscosity from observations for glacial isostatic adjustment (GIA) process has usually been conducted through the analyses based on the simple three-layer viscosity model characterized by lithospheric thickness, upper- and lower-mantle viscosities. Here, we examine the viscosity structures for the simple three-layer viscosity model and also for the two-layer lower-mantle viscosity model defined by viscosities of 670, D (670- D km depth) and D ,2891 ( D -2891 km depth) with D -values of 1191, 1691 and 2191 km. The upper-mantle rheological parameters for the two-layer lower-mantle viscosity model are the same as those for the simple three-layer one. For the simple three-layer viscosity model, rate of change of degree-two zonal harmonics of geopotential due to GIA process (GIA-induced 2 ) of –(6.0–6.5)  x  10 –11 yr –1 provides two permissible viscosity solutions for the lower mantle, (7–20)  x  10 21 and (5–9)  x  10 22  Pa s, and the analyses with observational constraints of the 2 and Last Glacial Maximum (LGM) sea levels at Barbados and Bonaparte Gulf indicate (5–9)  x  10 22  Pa s for the lower mantle. However, the analyses for the 2 based on the two-layer lower-mantle viscosity model only require a viscosity layer higher than (5–10)  x  10 21  Pa s for a depth above the core–mantle boundary (CMB), in which the value of (5–10)  x  10 21  Pa s corresponds to the solution of (7–20)  x  10 21  Pa s for the simple three-layer one. Moreover, the analyses with the 2 and LGM sea level constraints for the two-layer lower-mantle viscosity model indicate two viscosity solutions: 670,1191  〉 3  x  10 21 and 1191,2891  ~ (5–10)  x  10 22  Pa s, and 670,1691  〉 10 22 and 1691,2891  ~ (5–10)  x  10 22  Pa s. The inferred upper-mantle viscosity for such solutions is (1–4)  x  10 20  Pa s similar to the estimate for the simple three-layer viscosity model. That is, these analyses require a high viscosity layer of (5–10)  x  10 22  Pa s at least in the deep mantle, and suggest that the GIA-based lower-mantle viscosity structure should be treated carefully in discussing the mantle dynamics related to the viscosity jump at ~670 km depth. We also preliminarily put additional constraints on these viscosity solutions by examining typical relative sea level (RSL) changes used to infer the lower-mantle viscosity. The viscosity solution inferred from the far-field RSL changes in the Australian region is consistent with those for the 2 and LGM sea levels, and the analyses for RSL changes at Southport and Bermuda in the intermediate region for the North American ice sheets suggest the solution of 670, D  〉 10 22 , D ,2891  ~ (5–10)  x  10 22  Pa s ( D  = 1191 or 1691 km) and upper-mantle viscosity higher than 6  x  10 20  Pa s.

Description: Mid-crustal conductors are a common phenomenon in magnetotelluric studies. In the Andean Cordillera of southern Chile, they appear to concentrate along major fault zones. A high-resolution, broad-band magnetotelluric survey including 31 stations has been carried out along two profiles perpendicular to (1) the Liquiñe-Ofqui Fault Systems (LOFS) and (2) the Villarrica-Quetrupillán-Lanín volcanic lineament running parallel to the Mocha-Villarrica Fault Zone (MVFZ). The survey aimed at tracing one of the known conductors from mid-crustal depth to near-surface along these faults. Directionality and dimensionality were analysed using tensor decomposition. Phase tensors and induction arrows reveal two major geoelectric strike directions following the strike of LOFS and MVFZ. 2-D inversion shows low resistivity zones along both fault systems down to a depth of 〉10 km, where the brittle-ductile transition is expected. Along the LOFS, the two anomalies are linked to (1) Lake Caburgua, where the LOFS broadens to about 2 km of lateral extension and seems to represent a pull-apart structure, and (2) the intersection with the Villarrica-Quetrupillán-Lanín volcanic lineament, where seismic activity was observed during the latest eruption in March 2015. A connection of the mid-crustal conductor to the ESE-WNW-striking fault zones is indicated from the presented data.

Description: We present new rotations that describe the relative positions and velocities of the Pacific and North America plates at 22 times during the past 19.7 Myr, offering 1-Myr temporal resolution for studies of the geotectonic evolution of western North America and other plate boundary locations. Derived from 18 000 magnetic reversal, fracture zone and transform fault identifications from the Pacific–Antarctic–Nubia–North America plate circuit and the velocities of 935 GPS sites on the Pacific and North America plates, the new rotations and GPS-derived angular velocity indicate that the rate of motion between the two plates increased by 70 per cent from 19.7 to 9±1 Ma, but changed by less than 2 per cent since 8 Ma and even less since 4.2 Ma. The rotations further suggest that the relative plate direction has rotated clockwise for most of the past 20 Myr, with a possible hiatus from 9 to 5 Ma. This conflicts with previously reported evidence for a significant clockwise change in the plate direction at 8–6 Ma. Our new rotations indicate that Pacific plate motion became obliquely convergent with respect to the San Andreas Fault of central California at 5.2–4.2 Ma, in agreement with geological evidence for a Pliocene onset of folding and faulting in central California. Our reconstruction of the northern Gulf of California at 6.3 Ma differs by only 15–30 km from structurally derived reconstructions after including 3–4 km Myr –1 of geodetically measured slip between the Baja California Peninsula and Pacific plate. This implies an approximate 15–30 km upper bound for plate non-rigidity integrated around the global circuit at 6.3 Ma. A much larger 200±54 km discrepancy between our reconstruction of the northern Gulf of California at 12 Ma and that estimated from structural and marine geophysical observations suggests that faults in northwestern Mexico or possibly west of the Baja California Peninsula accommodated large amounts of obliquely divergent dextral shear from 12–6.3 Ma. Pacific–North America plate motion since 16 Myr estimated with our new rotations agrees well with structurally summed deformation along two transects of western North America between the Colorado Plateau and western California, with a difference as small as 40 km out of 760 km of margin-parallel motion. A strong resemblance between a 20-Myr-to-present flow line reconstructed with our new rotations and the traces of the 700-km-long Queen Charlotte Fault and continental slope west of Canada suggests that the plate margin geometry was influenced by the passage of the Pacific plate and Yakutat block. The new rotations also suggest that (1) oblique convergence west of Canada initiated at 12–11 Ma, 5–8 Myr earlier than previously estimated, (2) no significant margin-normal shortening has occurred in areas of Canada located east of the Haida Gwaii archipelago since 20 Ma and (3) Pacific plate underthrusting of Haida Gwaii has accommodated the margin-normal component of plate motion since 12–11 Ma. Our rotations suggest an 70 per cent increase in the rate that the Pacific plate has been consumed by subduction beneath the Aleutian arc since 19.7 Ma, with still-unknown consequences for the rate of arc magmatism.

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: Among the outstanding tectonic questions regarding the convergence between the Tien Shan and Tarim basin in northwestern China are the manner in which deformation is accommodated within their lithospheres, and the extent that the Tarim lithosphere underthrusts the Tien Shan. In particular, the amount and type of deformation within the Tarim basin is poorly understood. It is also uncertain if the convergence between the Tarim and the Tien Shan takes place mainly along a discrete boundary, or if the Tarim lithosphere simply indents into the Kazach shield, forming the Tien Shan through crustal thickening accommodated by a distributed series of thrust faults. In this study we use hypocentres from published earthquake catalogues and waveforms recorded by regional seismic networks to determine earthquake source parameters through regional centroid moment tensor inversion. The entire dataset consists of 160 earthquakes that occurred between 1969 and 2009 and with moment magnitudes between 3.5 and 7 distributed throughout the central Tien Shan and northwestern Tarim Basin. The estimated focal depths of these earthquakes range from the near-surface to about 44 km. Focal mechanisms throughout much of the Tien Shan indicate active deformation accommodated by thrust faults from at least the upper crust to 30 km depth. South of the Tien Shan, the Jia-shi earthquake sequence within the Tarim basin suggests that both crustal shortening and localized flexure are part of a complicated process involving rotational convergence. Inside the Tarim basin, two earthquakes with thrust faulting mechanisms near the crust–mantle boundary beneath the Bachu uplift imply a brittle rheology of the lower crust. High-angle thrust events occur broadly across the Tien Shan, suggesting that the Tarim lithosphere as a whole is strong and indents into the Kazach shield to create the mountain range.

Description: The subduction of lithospheric plates is partitioned between subducting plate motion and lateral slab migration (i.e. trench retreat and advance). We use 3-D, dynamic models of subduction to address the role of a power-law mantle viscosity on subduction dynamics and, in particular, rates of trench retreat. For all numerical models tested, we find that a power-law rheology results in reduced rates of trench retreat, and elevated slab dip angles, relative to the equivalent isoviscous mantle model. We analyse the asthenospheric pressure distribution and the style of mantle flow, which exhibits only limited variability as a function of mantle rheology, in order to compute estimates of the mantle forces associated with subduction. The inclusion of a power-law rheology reduces the mantle shear force (which resists subducting plate motion) to a greater degree than it reduces the dynamic pressure gradient across the slab (which resists trench retreat). Therefore, the inclusion of a power-law mantle rheology favours a shift towards a subduction mode with a reduced trench retreat component, typically a relative reduction of order 25 per cent in our 3-D models. We suggest that this mechanism may be of importance for reducing the high trench retreat rates observed in many previous models to levels more in line with the average subduction partitioning observed on Earth at present (i.e. trench velocity ≤ plate velocity), for most absolute plate motion reference frames.

Description: Pulverized rocks (PR) are extremely incohesive and highly fractured rocks found within the damage zones of several large strike-slip faults around the world. They maintain their crystal structure, show little evidence of shearing or chemical alteration, and are believed to be produced by strong tensile forces. Several mechanisms for pulverization have been proposed based on simple qualitative analyses or laboratory experiments under simplified loading conditions. Numerical modelling, however, can offer new insights into what is needed to produce PR and likely conditions of formation. We perform dynamic rupture simulations of different earthquakes, varying the magnitude, the slip distribution, and the rupture speed (supershear and subshear), while measuring the stresses produced away from the fault. To contextualize our results, a basic threshold of 10 MPa is set as the tensile strength of the rock mass and recordings are made of where, when, and by how much this threshold is exceeded for each earthquake type. Guided by field observations, we discern that a large (〉 M w 7.1) subshear earthquake along a bimaterial fault produces a pulverized rock distribution most consistent with observations. The damage is asymmetric with the majority on the stiffer side of the fault extending out for several hundred metres. Within this zone there is a large and sudden volumetric expansion in all directions as the rupture passes. We propose that such an extreme tensile stress state, repeated for every earthquake, eventually produces the PR seen in the field.

Description: Kinematics of divergent boundaries and Rift-Rift-Rift junctions are classically studied using long-term geodetic observations. Since significant magma-related displacements are expected, short-term deformation provides important constraints on the crustal mechanisms involved both in active rifting and in transfer of extensional deformation between spreading axes. Using InSAR and GPS data, we analyse the surface deformation in the whole Central Afar region in detail, focusing on both the extensional deformation across the Quaternary magmato-tectonic rift segments, and on the zones of deformation transfer between active segments and spreading axes. The largest deformation occurs across the two recently activated Asal-Ghoubbet (AG) and Manda Hararo-Dabbahu (MH-D) magmato-tectonic segments with very high strain rates, whereas the other Quaternary active segments do not concentrate any large strain, suggesting that these rifts are either sealed during interdyking periods or not mature enough to remain a plate boundary. Outside of these segments, the GPS horizontal velocity field shows a regular gradient following a clockwise rotation of the displacements from the Southeast to the East of Afar, with respect to Nubia. Very few shallow creeping structures can be identified as well in the InSAR data. However, using these data together with the strain rate tensor and the rotations rates deduced from GPS baselines, the present-day strain field over Central Afar is consistent with the main tectonic structures, and therefore with the long-term deformation. We investigate the current kinematics of the triple junction included in our GPS data set by building simple block models. The deformation in Central Afar can be described by adding a central microblock evolving separately from the three surrounding plates. In this model, the northern block boundary corresponds to a deep EW-trending trans-tensional dislocation, locked from the surface to 10–13 km and joining at depth the active spreading axes of the Red Sea and the Aden Ridge, from AG to MH-D rift segments. Over the long-term, this plate configuration could explain the presence of the en-échelon magmatic basins and subrifts. However, the transient behaviour of the spreading axes implies that the deformation in Central Afar evolves depending on the availability of magma supply within the well-established segments.

Description: Volcanic crises are often preceded and accompanied by volcano deformation caused by magmatic and hydrothermal processes. Fast and efficient model identification and parameter estimation techniques for various sources of deformation are crucial for process understanding, volcano hazard assessment and early warning purposes. As a simple model that can be a basis for rapid inversion techniques, we present a compound dislocation model (CDM) that is composed of three mutually orthogonal rectangular dislocations (RDs). We present new RD solutions, which are free of artefact singularities and that also possess full rotational degrees of freedom. The CDM can represent both planar intrusions in the near field and volumetric sources of inflation and deflation in the far field. Therefore, this source model can be applied to shallow dikes and sills, as well as to deep planar and equidimensional sources of any geometry, including oblate, prolate and other triaxial ellipsoidal shapes. In either case the sources may possess any arbitrary orientation in space. After systematically evaluating the CDM, we apply it to the co-eruptive displacements of the 2015 Calbuco eruption observed by the Sentinel-1A satellite in both ascending and descending orbits. The results show that the deformation source is a deflating vertical lens-shaped source at an approximate depth of 8 km centred beneath Calbuco volcano. The parameters of the optimal source model clearly show that it is significantly different from an isotropic point source or a single dislocation model. The Calbuco example reflects the convenience of using the CDM for a rapid interpretation of deformation data.

Description: Tehran, the capital city of Iran with a population of over 12 million, is one of the largest urban centres within the seismically active Alpine–Himalayan orogenic belt. Although several historic earthquakes have affected Tehran, their relation to individual faults is ambiguous for most. This ambiguity is partly due to a lack of knowledge about the locations, geometries and seismic potential of structures that have been obscured by dramatic urban growth over the past three decades, and which have covered most of the young geomorphic markers and natural exposures. Here we use aerial photographs from 1956, combined with an ~1 m DEM derived from stereo Pleiades satellite imagery to investigate the geomorphology of a growing anticline above a thrust fault—the Pardisan thrust—within central Tehran. The topography across the ridge is consistent with a steep ramp extending from close to the surface to a depth of ~2 km, where it presumably connects with a shallow-dipping detachment. No primary fault is visible at the surface, and it is possible that the faulting dissipates in the near surface as distributed shearing. We use optically stimulated luminescence to date remnants of uplifted and warped alluvial deposits that are offset vertically across the Pardisan fault, providing minimum uplift and slip-rates of at least 1 mm yr –1 . Our study shows that the faults within the Tehran urban region have relatively rapid rates of slip, are important in the regional tectonics, and have a great impact on earthquake hazard assessment of the city and surrounding region.

Description: The 3-D shear velocity structure beneath South India's Dharwar Craton determined from fundamental mode Rayleigh waves phase velocities reveals the existence of anomalously high velocity materials in the depth range of 50–100 km. Tomographic analysis of seismograms recorded on a network of 35 broad-band seismographs shows the uppermost mantle shear wave speeds to be as high as 4.9 km s –1 in the northwestern Dharwar Craton, decreasing both towards the south and the east. Below ~100 km, the shear wave speed beneath the Dharwar Craton is close to the global average shear wave speed at these depths. Limitations of usable Rayleigh phase periods, however, have restricted the analysis to depths of 120 km, precluding the delineation of the lithosphere–asthenosphere boundary in this region. However, pressure–temperature analysis of xenoliths in the region suggests a lithospheric thickness of at least ~185 km during the mid-Proterozoic period. The investigations were motivated by a search for seismic indicators in the shallow mantle beneath the distinctly different parts of the Dharwar Craton otherwise distinguished by their lithologies, ages and crustal structure. Since the ages of cratonic crust and of the associated mantle lithosphere around the globe have been found to be broadly similar and their compositions bimodal in time, any distinguishing features of the various parts of the Dharwar shallow mantle could thus shed light on the craton formation process responsible for stabilizing the craton during the Meso- and Neo-Archean.

Description: Frontier hydrocarbon development projects in the deepwater slopes of the Gulf of Mexico Basin, Santos Basin and Lower Congo Basin all require wells to cross ductile layers of autochthonous or allochthonous salt moving at peak rates of 100 mm yr –1 . The Couette–Poiseuille number is introduced here to help pinpoint the depth of shear stress reversal in such salt layers. For any well-planned through salt, the probable range of creep forces of moving salt needs to be taken into account when designing safety margins and load-factor tolerance of the well casing. Drag forces increase with wellbore diameter, but more significantly with effective viscosity and speed of the creeping salt layer. The potential drag forces on cased wellbores in moving salt sheets are estimated analytically using a range of salt viscosities (10 15 –10 19 Pa s) and creep rates (0–10 mm yr –1 ). Drag on perfectly rigid casing of infinite strength may reach up to 13 Giga Newton per meter wellbore length in salt having a viscosity of 10 19 Pa s. Well designers may delay stress accumulations due to salt drag when flexible casing accommodates some of the early displacement and strain. However, all creeping salt could displace, fracture and disconnect well casing, eventually. The shear strength of typical heavy duty well casing (about 1000 MPa) can be reached due to drag by moving salt. Internal flow of salt will then fracture the casing near salt entry and exit points, but the structural damage is likely to remain unnoticed early in the well-life when the horizontal shift of the wellbore is still negligibly small (at less than 1 cm yr –1 ). Disruption of casing and production flow lines within the anticipated service lifetime of a well remains a significant risk factor within distinct zones of low-viscosity salt which may reach ultrafast creep rates of 100 mm yr –1 .

Description: The Antarctic Peninsula (AP) consists of a long lived and uniquely well preserved magmatic arc system. The broad tectonic structure of the AP arc is well understood. However, magmatic processes occurring along the arc are only constrained by regional geophysical and relatively sparse geological data. Key questions remain about the timing, volume, and structural controls on magma emplacement. We present new high resolution aeromagnetic data across Adelaide Island, on the western margin of the AP revealing the complex structure of the AP arc/forearc boundary. Using digital enhancement, 2-D modelling and 3-D inversion we constrain the form of the magnetic sources at the arc/forearc boundary. Our interpretation of these magnetic data, guided by geological evidence and new zircon U-Pb dating, suggests significant Palaeogene to Neogene magmatism formed ~25 per cent of the upper crust in this region (~7500 km 3 ). Significant structural control on Neogene magma emplacement along the arc/forearc boundary is also revealed. We hypothesize that this Neogene magmatism reflects mantle return flow through a slab window generated by Late Palaeogene cessation of subduction south of Adelaide Island. This mantle process may have affected the final stages of arc magmatism along the AP margin.

Description: The intrusion mechanism and internal structure of sills are still under debate. We present a detailed magnetic study, including anisotropy of magnetic susceptibility and rock magnetic analyses of a Cretaceous (94 Ma), 7-m-thick sill from the Lusitanian Basin in Portugal, the Foz da Fonte sill. The results, from both the top surface and a vertical profile, allow us to propose a model for the magmatic flow pattern and sense of flow. According to their location in the vertical profile, three magnetic fabric domains are identified: (1) at the borders, qualified as chilled margins (~0–50 cm), low anisotropies suggest that low velocity gradients and heterogeneous flow paths occurred during the initial emplacement stages; (2) in the centre of the sill, where low anisotropies are observed, low shear gradients and magma displacement close to pure translation is inferred and (3) in the intermediate zones, high anisotropy values are ascribed to zones having undergone high shear gradients. The mean magnetic lineations from the top surface and basal contact indicate an almost horizontal and NW–SE orientation (azimuth: 310°) which agrees with the preferred orientation of iron oxide grain clusters and with the elongation of vesicles considered as coaxial with the magma flow direction. Moreover, the magnetic foliation planes and the lineations show both a mirror imbrication relative to the average upper and lower border surfaces of the sill, pointing to a flow direction towards the SE. Based on these results and on the interpretation of two seismic reflection lines, we show that the Cabo Raso magnetic anomaly, located 25 km to NW of the FF-sill, is associated to Cretaceous magmatic intrusions from which the sill likely originated. This tectono-magmatic setting is discussed with respect to the West Iberia Late Cretaceous magmatism, integrating magnetic anomalies, isotope chronology and tectonics.

Description: Southern Mendoza and northern Neuquén Provinces, south of the Pampean Shallow Subduction region in western Argentina, are host to the 〈2 Myr Payunia Basaltic Province, which covers ~39 500 km 2 with primarily basaltic intraplate volcanism. This backarc igneous province can be explained by extension due to trench roll-back following steepening of a flat slab that existed in the middle to late Miocene. Magnetotelluric data collected at 37 sites from 67°W to 70°W and 35°S to 38°S are used to probe the source of the Payún Matrú basalts. These data, which require significantly 3-D structure, are inverted with a 3-D non-linear conjugate gradient algorithm that minimizes structure for a given data misfit. We identify two significant electrically conductive structures. One, called the SWAP (shallow western asthenospheric plume) approaches the surface beneath the Payún Matrú Caldera and the Trómen Volcano and dips westward towards the subducted Nazca slab. The second, called the DEEP (deep eastern plume) approaches the surface ~100 km to the southeast of Payún Matrú and dips steeply east to ~400 km depth while remaining above the subducted Nazca slab. We use a variety of model assessment techniques including forward modelling and constrained inversion to test the veracity of these features. We interpret the SWAP as the source of the 〈2 Myr intraplate volcanism. Our model assessment permits but does not require the SWAP to connect to the Nazca slab. The SWAP and DEEP are electrically connected only in the shallow crust, which is likely due to the Neuquén sedimentary basin and not a magmatic process. We propose that the SWAP and DEEP may have been more robustly connected in the past, but that the DEEP was decapitated to form the SWAP when shallow northwestward mantle flow resumed during steepening of the slab. The ~2 Myr basaltic volcanism is the result of this decapitated DEEP magma that had ponded below the crust until extension allowed eruption. The westward dip of the SWAP is interpreted to be the result of shear in the renewed mantle corner flow—this explains why the SWAP and Nazca slab can appear connected, yet there is no recent arc-signature magma in this region.

Description: This study presents the results of a deep seismic survey across the north Algerian margin, based on the combination of 2-D multichannel and wide-angle seismic data simultaneously recorded by 41 ocean bottom seismometers deployed along a north–south line extending 180 km off Jijel into the Algerian offshore basin, and 25 land stations deployed along a 100-km-long line, cutting through the Lesser Kabylia and the Tellian thrust-belt. The final model obtained using forward modelling of the wide-angle data and pre-stack depth migration of the seismic reflection data provides an unprecedented view of the sedimentary and crustal structure of the margin. The sedimentary layers in the Algerian basin are 3.75 km thick to the north and up to 4.5–5 km thick at the foot of the margin. They are characterized by seismic velocities from 1.9 to 3.8 km s –1 . Messinian salt formations are about 1 km thick in the study area, and are modelled and imaged using a velocity between 3.7 and 3.8 km s –1 . The crust in the deep sea basin is about 4.5 km thick and of oceanic origin, presenting two distinct layers with a high gradient upper crust (4.7–6.1 km s –1 ) and a low gradient lower crust (6.2–7.1 km s –1 ). The upper-mantle velocity is constrained to 7.9 km s –1 . The ocean–continent transition zone is very narrow between 15 and 20 km wide. The continental crust reaches 25 km thickness as imaged from the most landward station and thins to 5 km over a less than 70 km distance. The continental crust presents steep and asymmetric upper- and lower-crustal geometry, possibly due to either asymmetric rifting of the margin, an underplated body, or flow of lower crustal material towards the ocean basin. Present-time deformation, as imaged from three additional seismic profiles, is characterized by an interplay of gravity-driven mobile-salt creep and active thrusting at the foot of the tectonically inverted Algerian margin.

Description: Melt generation and migration are an important link between surface processes and the thermal and chemical evolution of the Earth's interior. However, their vastly different timescales make it difficult to study mantle convection and melt migration in a unified framework, especially for 3-D global models. And although experiments suggest an increase in melt volume of up to 20 per cent from the depth of melt generation to the surface, previous computations have neglected the individual compressibilities of the solid and the fluid phase. Here, we describe our extension of the finite element mantle convection code ASPECT that adds melt generation and migration. We use the original compressible formulation of the McKenzie equations, augmented by an equation for the conservation of energy. Applying adaptive mesh refinement to this type of problems is particularly advantageous, as the resolution can be increased in areas where melt is present and viscosity gradients are high, whereas a lower resolution is sufficient in regions without melt. Together with a high-performance, massively parallel implementation, this allows for high-resolution, 3-D, compressible, global mantle convection simulations coupled with melt migration. We evaluate the functionality and potential of this method using a series of benchmarks and model setups, compare results of the compressible and incompressible formulation, and show the effectiveness of adaptive mesh refinement when applied to melt migration. Our model of magma dynamics provides a framework for modelling processes on different scales and investigating links between processes occurring in the deep mantle and melt generation and migration. This approach could prove particularly useful applied to modelling the generation of komatiites or other melts originating in greater depths. The implementation is available in the Open Source ASPECT repository.

Description: We use recently published, high-resolution reconstructions of the Southwest Indian Ridge to test whether a previously described systematic difference between Global Positioning System (GPS) and 3.16-Myr-average estimates of seafloor spreading rates between Antarctica and Africa is evidence for a recent slowdown in Southwest Indian Ridge seafloor spreading rates. Along the Nubia-Antarctic segment of the ridge, seafloor opening rates that are estimated with the new, high-resolution reconstructions and corrected for outward displacement agree well with geodetic rate estimates and reduce previously reported, highly significant non-closure of the Nubia-Antarctic-Sur plate circuit. The observations are inconsistent with a slowdown in spreading rates and instead indicate that Nubia-Antarctic plate motion has been steady since at least 5.2 Ma. Lwandle-Antarctic seafloor spreading rates that are estimated from the new high-resolution reconstructions differ insignificantly from a GPS estimate, thereby implying steady Lwandle-Antarctic plate motion since 5.2 Ma. Between the Somalia and Antarctic plates, the new Southwest Indian Ridge reconstructions eliminate roughly half of the systematic difference between the GPS and MORVEL spreading rate estimates.We interpret the available observations as evidence that Somalia-Antarctic spreading rates have been steady since at least 5.2 Ma and postulate that the remaining difference is attributable to random and/or systematic errors in the plate kinematic estimates and the combined effects of insufficient geodetic sampling of undeforming areas of the Somalia plate, glacial isostatic adjustment in Antarctica and transient deformation triggered by the 1998 M w = 8.2 Antarctic earthquake, the 2004 M w = 9.3 Sumatra earthquake, or possibly other large historic earthquakes.

Description: Geophysical data from the MEDOC experiment across the Northern Tyrrhenian backarc basin has mapped a failed rift during backarc extension of cratonic Variscan lithosphere. In contrast, data across the Central Tyrrhenian have revealed the presence of magmatic accretion followed by mantle exhumation after continental breakup. Here we analyse the MEDOC transect E–F, which extends from Sardinia to the Campania margin at 40.5°N, to define the distribution of geological domains in the transition from the complex Central Tyrrhenian to the extended continental crust of the Northern Tyrrhenian. The crust and uppermost mantle structure along this ~400-km-long transect have been investigated based on wide-angle seismic data, gravity modelling and multichannel seismic reflection imaging. The P -wave tomographic model together with a P -wave-velocity-derived density model and the multichannel seismic images reveal seven different domains along this transect, in contrast to the simpler structure to the south and north. The stretched continental crust under Sardinia margin abuts the magmatic crust of Cornaglia Terrace, where accretion likely occurred during backarc extension. Eastwards, around Secchi seamount, a second segment of thinned continental crust (7–8 km) is observed. Two short segments of magmatically modified continental crust are separated by the ~5-km-wide segment of the Vavilov basin possibly made of exhumed mantle rocks. The eastern segment of the 40.5°N transect E–F is characterized by continental crust extending from mainland Italy towards the Campania margin. Ground truthing and prior geophysical information obtained north and south of transect E–F was integrated in this study to map the spatial distribution of basement domains in the Central Tyrrhenian basin. The northward transition of crustal domains depicts a complex 3-D structure represented by abrupt spatial changes of magmatic and non-magmatic crustal domains. These observations imply rapid variations of magmatic activity difficult to reconcile with current models of extension of continental lithosphere essentially 2-D over long distances.

Description: Northern Attica in Greece is characterized by a set of north dipping, subparallel normal faults. These faults were considered to have low tectonic activity, based on historical earthquake reports, instrumental seismicity and slip rate estimates. This study presents new data for one of these faults, the Milesi Fault. We run GIS based geomorphological analyses on fault offset distribution, field mapping of postglacial fault scarps and ground penetrating radar profiling to image hangingwall deformation. The first palaeoseismological trenching in this part of Greece allowed obtaining direct data on slip rates and palaeoearthquakes. The trenching revealed downthrown and buried palaeosols, which were dated by radiocarbon. The results of our investigations show that the slip rates are higher than previously thought and that at least four palaeoearthquakes with magnitudes of around M 6.2 occurred during the last 4000–6000 yr. We calculate an average recurrence interval of 1000–1500 yr and a maximum throw rate of ~0.4–0.45 mm a –1 . Based on the new geological earthquake data we developed a seismic hazard scenario, which also incorporates geological site effects. Intensities up to IX must be expected for Northern Attica and the southeastern part of Evia. Earthquake environmental effects like liquefaction and mass movements are also likely to occur. This scenario is in contrast to the official Greek seismic hazard zonation that is based on historical records and assigns different hazard zones for municipalities that will experience the same intensity by earthquakes on the Milesi Fault. We show that the seismic hazard is likely underestimated in our study area and emphasize the need to incorporate geological information in such assessments.

Description: Numerical models of mantle convection typically employ a temperature- or pressure-dependent viscous or viscoplastic rheology and a free slip upper boundary condition. The Earth, however, has a stress-free rather than a free slip surface condition. In addition, with decreasing temperature, the viscosity of rocks increases, which might induce a change from viscous to elastic behaviour (depending on the timescale of deformation). Here, we study the effects of both a Maxwell viscoelastic rheology and a free surface upper boundary condition on viscoelastic convection with a strongly temperature dependent rheology. We particularly focus on the effect of elasticity on the stress state of the lithosphere. Results show that convection vigor and heat transport are not significantly altered by the upper boundary condition or by elasticity. However, the stress state of the lithosphere is significantly affected by both factors. If elasticity is unimportant, a free surface upper boundary condition results in significantly elevated surface stresses (which are up to two magnitudes larger than in the free slip case). Elasticity counteracts this effect and significantly reduces the surface stresses, but distributes stresses over a thicker layer than in the case of a purely viscous rheology. At Earth-like conditions, this effect is significant. While it is warranted to use a free slip upper boundary condition and neglecting elasticity when studying mantle convection and its effect on the thermal state of the Earth, both factors are significant when one wants to predict the stress state of the lithosphere and related questions. Additional 2-D simulations of a plume impinging on a constant thickness and constant viscosity lithosphere show that reasonable parameters might induce lithospheric stress levels that are on the order of a GPa or larger, for viscous free surface models, and that these stresses are several orders of magnitude larger than stresses that occur for free slip models. This suggests that the effect of a free surface should not be ignored in models where the rheology is stress-dependent, such as in viscoplastic models of self-consistent plate tectonics.

Description: On 2012 August 11, a pair of large, damaging earthquakes struck the Varzaghan–Ahar region in northwest Iran, in a region where there was no major mapped fault or any well-documented historical seismicity. To investigate the active tectonics of the source region we applied a combination of seismological methods (local aftershock network, calibrated multiple event relocation and focal mechanism studies), field observations (structural geology and geomorphological) and inversions for the regional stress field. The epicentral region is north of the North Tabriz Fault. The first main shock is characterized by right-lateral strike-slip motion on an almost E–W fault plane of about 23 km length extending from the surface to a depth of about 14 km. The second main shock occurred on an ENE-striking fault that dips at 60–70° to the NW. Independent inversions of focal mechanisms and geologically determined fault kinematic data for the active stress state yield a transpressional tectonic regime with 1 oriented N132E. For the region northeast of the North Tabriz Fault, the presence of rigid lithosphere of the South Caspian Basin implies the kinematic adjustment by northward transferring of the contracted masses through both distributed deformation and structural deflections. Our results suggest that the kinematic adjustment inside a contracting wedge may occur along interacting crosswise or conjugate faults to accommodate low rates of internal deformation. At a global scale, our results indicate that despite the basic assumption of ‘rigid blocks’ in geodetic plate modelling, internal deformation of block-like regions could control the kinematics of deformation and the level of seismic hazard within and around such regions of low deformation rate.

Description: Many salt diapirs are thought to have formed as a result of down-building, which implies that the top of the diapir remained close to the surface during syn-halokinetic sediment deposition. Down-building is largely a 3-D process and in order to better understand what controls the patterns of the diapirs that form by this process, we here perform 3-D numerical models of down-built diapirs initiated by the gravity instability in linear viscous materials and compare the results with analytical models. We vary several parameters of the numerical models such as initial salt thickness, sedimentation rate, salt viscosity, salt-sediment viscosity ratio as well as the density of sediments. Down-building of 3-D diapirs only occurs for a certain range of parameters and is favoured by lower sediment/salt viscosity contrasts and sedimentation rates in agreement with analytical predictions and findings from previous 2-D models. However, the models show that the sedimentation rate has an additional effect on the formation and evolution of 3-D diapir patterns. At low sedimentation rates, salt ridges that form during early model stages remain preserved at later stages as well. For higher sedimentation rates, the initial salt ridges are covered up and finger-like diapirs form at their junctions, which results in different salt exposure patterns at the surface. Once the initial pattern of diapirs is formed, higher sedimentation rate can also result in covered diapirs if the diapir extrusion velocity is insufficiently large. We quantify the effect of sedimentation rate on the number of diapirs exposed at the surface as well as on their spacing and we explain the observations with analytical predictions using thick-plate analytical models. In some cases, this final pattern is distinctly different from the initial polygonal pattern.

Description: We present a generalized formalism for computing gravitationally self-consistent sea level changes driven by the combined effects of dynamic topography, geoid perturbations due to mantle convection, ice mass fluctuations and sediment redistribution on a deforming Earth. Our mathematical treatment conserves mass of the surface (ice plus ocean) load and the solid Earth. Moreover, it takes precise account of shoreline migration and the associated ocean loading. The new formalism avoids a variety of approximations adopted in previous models of sea level change driven by dynamic topography, including the assumption that a spatially fixed isostatic amplification of ‘air-loaded’ dynamic topography accurately accounts for ocean loading effects. While our approach is valid for Earth models of arbitrary complexity, we present numerical results for a set of simple cases in which a pattern of dynamic topography is imposed, the response to surface mass loading assumes that Earth structure varies only with depth and that isostatic equilibrium is maintained at all times. These calculations, involving fluid Love number theory, indicate that the largest errors in previous predictions of sea level change driven by dynamic topography occur in regions of shoreline migration, and thus in the vicinity of most geological markers of ancient sea level. We conclude that a gravitationally self-consistent treatment of long-term sea level change is necessary in any effort to use such geological markers to estimate ancient ice volumes.

Description: The Gibraltar arc and surrounding areas are a complex tectonic region and its tectonic evolution since Miocene is still under debate. Knowledge of its lithospheric structure will help to understand the mechanisms that produced extension and westward motion of the Alboran domain, simultaneously with NW–SE compression driven by Africa–Europe plates convergence. We perform a P -wave receiver function analysis in which we analyse new data recorded at 83 permanent and temporary seismic broad-band stations located in the South of the Iberian peninsula. These data are stacked and combined with data from a previous study in northern Morocco to build maps of thickness and average v P / v S ratio for the crust, and cross-sections to image the lithospheric discontinuities beneath the Gibraltar arc, the Betic and Rif Ranges and their Iberian and Moroccan forelands. Crustal thickness values show strong lateral variations in the southern Iberia peninsula, ranging from ~19 to ~46 km. The Variscan foreland is characterized by a relatively flat Moho at ~31 km depth, and an average v P / v S ratio of ~1.72, similar to other Variscan terranes, which may indicate that part of the lower crustal orogenic root was lost. The thickest crust is found at the contact between the Alboran domain and the External Zones of the Betic Range, while crustal thinning is observed southeastern Iberia (down to 19 km) and in the Guadalquivir basin where the thinning at the Iberian paleomargin could be still preserved. In the cross-sections, we see a strong change between the eastern Betics, where the Iberian crust underthrusts and couples to the Alboran crust, and the western Betics, where the underthrusting Iberian crust becomes partially delaminated and enters into the mantle. The structures largely mirror those on the Moroccan side where a similar detachment was observed in northern Morocco. We attribute a relatively shallow strong negative-polarity discontinuity to the lithosphere-asthenosphere boundary. This means relatively thin lithosphere ranging from ~50 km thickness in southeastern Iberia and northeastern Morocco to ~90–100 km beneath the western Betics and the Rif, with abrupt changes of ~30 km under the central Betics and northern Morocco. Our observations support a geodynamic scenario where in western Betics oceanic subduction has developed into ongoing continental subduction/delamination while in eastern Betics this process is inactive.

Description: A series of our linear analysis on the onset of thermal convection was applied to that of highly compressible fluids in a planar layer whose thermal conductivity and viscosity vary in space, in order to study the influences of spatial variations in physical properties expected in the mantles of massive terrestrial planets. The thermal conductivity and viscosity are assumed to exponentially depend on depth and temperature, respectively, while the variations in thermodynamic properties (thermal expansivity and reference density) with depth are taken to be relevant for the super-Earths with 10 times the Earth's. Our analysis demonstrated that the nature of incipient thermal convection is strongly affected by the interplay between the adiabatic compression and spatial variations in physical properties of fluids. Owing to the effects of adiabatic compression, a ‘stratosphere’ can occur in the deep mantles of super-Earths, where a vertical motion is insignificant. An emergence of ‘stratosphere’ is greatly enhanced by the increase in thermal conductivity with depth, while it is suppressed by the decrease in thermal expansivity with depth. In addition, by the interplay between the static stability and strong temperature dependence in viscosity, convection cells tend to be confined in narrow regions around the ‘tropopause’ at the interface between the ‘stratosphere’ of stable stratification and the ‘troposphere’ of unstable stratification. We also found that, depending on the variations in physical properties, two kinds of stagnant regions can separately develop in the fluid layer. One is well-known ‘stagnant-lids’ of cold and highly viscous fluids, and the other is ‘basal stagnant regions’ of hot and less viscous fluids. The occurrence of ‘basal stagnant regions’ may imply that convecting motions can be insignificant in the lowermost part of the mantles of massive super-Earths, even in the absence of strong increase in viscosity with pressure (or depth).

Description: We present a 2-D subsurface image of the Paganica Fault from a high-resolution refraction tomography and detailed geological investigation carried out across part of the northwestern segment of the 20-km-long Paganica–San Demetrio fault-system, and which was responsible of the 2009 April 6 M w 6.1 L'Aquila earthquake (central Italy). We acquired two seismic profiles crossing the Paganica basin with a dense-wide aperture configuration. More than 30 000 P wave first-arrival traveltimes were input to a non-linear tomographic inversion. The obtained 250–300 m deep 2-D Vp images illuminate the shallow portion of the Paganica Fault, and depict additional unreported splays defining a complex half-graben structure. We interpret local thickening of low- Vp (〈2400 m s –1 ) and intermediate- Vp (2600–3400 m s –1 ) regions as syn-tectonic clastic wedges above a high- Vp (3800–5000 m s –1 ) carbonate basement. These results are condensed in a 4.2-km-long section across the Paganica basin, clearly indicating that the Paganica Fault is a mature normal fault cutting the whole upper ~10 km of the crust. We evaluate a minimum cumulative net displacement of 650 ± 90 m and a total heave of 530 ± 65 m accomplished by the Paganica Fault, respectively. In the conservative hypothesis that the extension started during the Gelasian (1.80–2.59 Ma), we obtain a minimum long-term slip-rate of 0.30 ± 0.07 mm yr –1 and an extension-rate of 0.25 ± 0.06 mm yr –1 , respectively. Considering the regional averaged extensional field of ~1 mm yr –1 obtained from geodetic and geological analyses at 10 4  yr timescale, we infer that the Paganica Fault accounts for ~20 per cent of the NE-extension affecting this zone of the central Apennines axis due to the concurrent activity of other parallel normal fault-systems nearby (e.g. the Liri, Velino-Magnola, L'Aquila-Celano and Gran Sasso fault-systems).

Description: The Mozambique Ridge, a prominent basement high in the southwestern Indian Ocean, consists of four major geomorphological segments associated with numerous phases of volcanic activity in the Lower Cretaceous. The nature and origin of the Mozambique Ridge have been intensely debated with one hypothesis suggesting a Large Igneous Province origin. High-resolution seismic reflection data reveal a large number of extrusion centres with a random distribution throughout the southern Mozambique Ridge and the nearby Transkei Rise. Intrabasement reflections emerge from the extrusion centres and are interpreted to represent massive lava flow sequences. Such lava flow sequences are characteristic of eruptions leading to the formation of continental and oceanic flood basalt provinces, hence supporting a Large Igneous Province origin of the Mozambique Ridge. We observe evidence for widespread post-sedimentary magmatic activity that we correlate with a southward propagation of the East African Rift System. Based on our volumetric analysis of the southern Mozambique Ridge we infer a rapid sequential emplacement between ~131 and ~125 Ma, which is similar to the short formation periods of other Large Igneous Provinces like the Agulhas Plateau.

Description: Using an up-to-date global plate rotation model, applied to the endpoints of preserved major spreading ridge isochrons, we have calculated the explicitly reconstructable length-weighted mean global half-spreading rate (HSR), ridge length and area production as a function of time since the end of the Cretaceous Normal Superchron at 83.0 Ma. Our calculations integrate uncertainties in rotation parameters and chron boundary ages with the partial sampling uncertainties arising from progressive subduction of older oceanic lithosphere and its preserved spreading record. This record of directly reconstructable oceanic ridge production provides a well-constrained baseline that can be compared to reconstructions that include the largely unconstrained extrapolated histories of entirely subducted oceanic plates. The directly reconstructable global mean HSR has not varied by more than ±15 per cent about an average rate of 28.4 ± 4.6 mm a –1 since 83 Ma. No long-term secular trend is evident: a maximum global mean half-rate of 32 ± 6 mm a –1 occurred from 33.1 Ma to about 25.8 Ma, with minima of 26 ± 5 mm a –1 between about 56 and 40.2 Ma, and 24 ± 1 mm a –1 since 3.2 Ma. Only this most recent interval has a rate that differs significantly (at ±2) from the long-term mean. The global, reconstructable ridge length at 56 Ma decreases by less than 15 per cent relative to the modern ridge system; by 83 Ma it has decreased by 38 per cent. These relatively high preserved ridge fractions mean that the estimated uncertainty due to partial sampling stays roughly equivalent to the estimated rotation model uncertainties, allowing long-term spreading rate variations of 〉20 per cent since the Late Cretaceous to be ruled out. In contrast, prior to 83 Ma too little oceanic lithosphere is preserved to reliably reconstruct global spreading rates.

Description: Three 2-D Deep Electrical Resistivity Tomography (ERT) transects, up to 6.36 km long, were obtained across the Paganica-San Demetrio Basin, bounded by the 2009 L'Aquila M w 6.1 normal-faulting earthquake causative fault (central Italy). The investigations allowed defining for the first time the shallow subsurface basin structure. The resistivity images, and their geological interpretation, show a dissected Mesozoic-Tertiary substratum buried under continental infill of mainly Quaternary age due to the long-term activity of the Paganica-San Demetrio normal faults system (PSDFS), ruling the most recent deformational phase. Our results indicate that the basin bottom deepens up to 600 m moving to the south, with the continental infill largely exceeding the known thickness of the Quaternary sequence. The causes of this increasing thickness can be: (1) the onset of the continental deposition in the southern sector took place before the Quaternary, (2) there was an early stage of the basin development driven by different fault systems that produced a depocentre in the southern sector not related to the present-day basin shape, or (3) the fault system slip rate in the southern sector was faster than in the northern sector. We were able to gain sights into the long-term PSDFS behaviour and evolution, by comparing throw rates at different timescales and discriminating the splays that lead deformation. Some fault splays exhibit large cumulative throws (〉300 m) in coincidence with large displacement of the continental deposits sequence (〉100 m), thus testifying a general persistence in time of their activity as leading splays of the fault system. We evaluate the long-term (3–2.5 Myr) cumulative and Quaternary throw rates of most of the leading splays to be 0.08–0.17 mm yr –1 , indicating a substantial stability of the faults activity. Among them, an individual leading fault splay extends from Paganica to San Demetrio ne’ Vestini as a result of a post-Early Pleistocene linkage of two smaller splays. This 15 km long fault splay can explain the Holocene surface ruptures observed to be larger than those occurred during the 2009 L'Aquila earthquake, such as revealed by palaeoseismological investigations. Finally, the architecture of the basin at depth suggests that the PSDFS can also rupture a longer structure at the surface, allowing earthquakes larger than M 6.5, besides rupturing only small sections, as it occurred in 2009.

Description: Previous statistical studies showed that there was a correlation between the ultralow frequency (ULF) seismo-magnetic phenomena and local seismicity in the Kakioka region, Japan. In this study, utilizing Molchan's error diagram, we evaluate whether these phenomena contain precursory information and discuss how they can be used in short-term forecasting of sizable earthquakes. In practice, for given series of precursory signals and related earthquake events, each prediction strategy is characterized by the leading time of alarms ( ) and the length of alarm window ( L ). The leading time is the time length between a detected anomaly and its following alarm, and the alarm window is the duration that an alarm lasts. A modified area skill score measuring the area between actual prediction curve and random prediction line in Molchan's error diagram is used to assess the efficiency of different prediction strategies. The results indicate that predictions based on ULF magnetic data in Kakioka, Japan perform better than random prediction when is around 1 week and L is less than 4 d or is 13–14 d and L is less than 1 week. The optimal strategy of short-term forecasts has been established by setting at 8 d and L at 1 d. The methodology proposed in this study could also be useful in evaluating the prediction policy and optimizing other kinds of measurements for short-term earthquake forecasting.

Description: Inner core convection, and the corresponding variations in grain size and alignment, has been proposed to explain the complex seismic structure of the inner core, including its anisotropy, lateral variations and the F-layer at the base of the outer core. We develop a parametrized convection model to investigate the possibility of convection in the inner core, focusing on the dominance of the plume mode of convection versus the translation mode. We investigate thermal and compositional convection separately so as to study the end-members of the system. In the thermal case the dominant mode of convection is strongly dependent on the viscosity of the inner core, the magnitude of which is poorly constrained. Furthermore recent estimates of a large core thermal conductivity result in stable thermal stratification, hindering convection. However, an unstable density stratification may arise due to the pressure dependant partition coefficient of certain light elements. We show that this unstable stratification leads to compositionally driven convection, and that inner core translation is likely to be the dominant convective mode due to the low compositional diffusivity. The style of convection resulting from a combination of both thermal and compositional effects is not easy to understand. For reasonable parameter estimates, the stabilizing thermal buoyancy is greater than the destabilizing compositional buoyancy. However we anticipate complex double diffusive processes to occur given the very different thermal and compositional diffusivities.

Description: Relative sea level histories from previously glaciated areas have been used to study Earth rheology and ice sheet evolution during the last glacial cycle. The analysis of postglacial decay times has been used to place estimates on Earth viscosity structure that are relatively independent of uncertainty in the local ice history. Reconstructed sea levels from Ångermanland, Sweden have been commonly adopted for this purpose. We have assessed and compiled an updated relative sea level curve for this region, combining both varve-dated and radiocarbon dated index points. We fitted an exponential curve to the observations, taking into account estimates of eustatic sea level rise, elevation uncertainties as well as the geographical spread of the data sites to arrive at a decay time range (2 ) of 4.2–4.9 kyr for the whole record length (0–8 kyr) and 4.2–6.2 when 0–7 kyr fits are included. We computed model decay times using a large suite of over 900 ice and earth model combinations based on over 400 three-layer Earth viscosity models and more than 30 ice history reconstructions. Based on these extensive results, we confirm that decay time estimates are relatively independent of the regional ice model (at least within the range of ice chronology uncertainties) and so this data parametrization provides a relatively robust measure of Earth viscosity structure. We find that the observational constraints listed above are satisfied by 29 (8 kyr record) and 52 (7 and 8 kyr record) of the viscosity models considered. These subsets define uncertainty ranges in upper and lower mantle viscosity that are interdependent (Fig. 5). Consistent with previous analyses, we find that the observational decay time estimate does not provide useful constraints on model lithospheric thickness (within the range explored; 46–146 km).

Description: Seismological, geological and geodetic data have been integrated to characterize the seismogenic structure of the late 2013-early 2014 moderate energy (maximum local magnitude M Lmax = 4.9) seismic sequence that struck the interior of the Matese Massif, part of the Southern Apennines active extensional belt. The sequence, heralded by a M L = 2.7 foreshock, was characterized by two main shocks with M L = 4.9 and M L = 4.2, respectively, which occurred at a depth of ~17–18 km. The sequence was confined in the 10–20 km depth range, significantly deeper than the 1997–1998 sequence which occurred few km away on the northeastern side of the massif above ~15 km depth. The depth distribution of the 2013–14 sequence is almost continuous, albeit a deeper (16–19 km) and a shallower (11–15 km) group of events can be distinguished, the former including the main shocks and the foreshock. The epicentral distribution formed a ~10 km long NNW–SSE trending alignment, which almost parallels the surface trace of late Pliocene–Quaternary southwest-dipping normal faults with a poor evidence of current geological and geodetic deformation. We built an upper crustal model profile for the eastern Matese massif through integration of geological data, oil exploration well logs and seismic tomographic images. Projection of hypocentres on the profile suggests that the seismogenic volume falls mostly within the crystalline crust and subordinately within the Mesozoic sedimentary cover of Apulia, the underthrust foreland of the Southern Apennines fold and thrust belt. Geological data and the regional macroseismic field of the sequence suggest that the southwest-dipping nodal plane of the main shocks represents the rupture surface that we refer to here as the Matese fault. The major lithological discontinuity between crystalline and sedimentary rocks of Apulia likely confined upward the rupture extent of the Matese fault. Repeated coseismic failure represented by the deeper group of events in the sequence, activated in a passive fashion the overlying ~11–15 km deep section of the upper crustal normal faults. We consider the southwest-dipping Matese fault representative of a poorly known type of seismogenic structures in the Southern Apennines, where extensional seismogenesis and geodetic strain accumulation occur more frequently on NE-dipping, shallower-rooted faults. This is the case of the Boiano Basin fault located on the northern side of the massif, to which the 1997–1998 sequence is related. The close proximity of the two types of seismogenic faults at the Matese Massif is related to the complex crustal architecture generated by the Pliocene–early Pleistocene contractional and transpressional tectonics.

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: For the first time, a deep seismic data set acquired in the frame of the Algerian–French SPIRAL program provides new insights regarding the origin of the westernmost Algerian margin and basin. We performed a tomographic inversion of traveltimes along a 100-km-long wide-angle seismic profile shot over 40 ocean bottom seismometers offshore Mostaganem (Northwestern Algeria). The resulting velocity model and multichannel seismic reflection profiles show a thin (3–4 km thick) oceanic crust. The narrow ocean–continent transition (less than 10 km wide) is bounded by vertical faults and surmounted by a narrow almost continuous basin filled with Miocene to Quaternary sediments. This fault system, as well as the faults organized in a negative-flower structure on the continent side, marks a major strike-slip fault system. The extremely sharp variation of the Moho depth (up to 45 ± 3°) beneath the continental border underscores the absence of continental extension in this area. All these features support the hypothesis that this part of the margin from Oran to Tenes, trending N65–N70°E, is a fossil subduction-transform edge propagator fault, vestige of the propagation of the edge of the Gibraltar subduction zone during the westward migration of the Alborán domain.

Description: New Zealand straddles the boundary between the Australian and Pacific Plate. Cenozoic relative plate motion has resulted in a complex pattern of faulting and block rotation in a zone of continental lithosphere up to 250 km wide. I investigate the implications of the short-term kinematics for the strength of the deforming lithosphere. I use a compilation of seismic reflection/refraction studies and high quality receiver function analyses to determine both the regional structure of the crust, which ranges from 20 to 50 km thick, and fields of buoyancy stress (or GPE per unit volume). Deformation over thousands of years is quantified in terms of velocity and strain rate fields, based on an inversion of neotectonic fault slip and palaeomagnetic data, in the context of the short-term relative plate motions. Forces on the subduction megathrust, as well as deviatoric stresses in the behind subduction region, are calculated from simple 2-D force balances across the Hikurangi Margin, given negligible deviatoric stresses at the along-strike transition between backarc extension and compression. Average megathrust shear stresses are in the range 6–15 MPa, and average lithospheric stresses 〈20 MPa in the overriding plate. The regional lithospheric strength of the plate boundary zone, assuming a viscous rheology (Newtonian or power law), is determined from an inversion of the field of gradients of buoyancy stress (averaged over either the top 25 km of the crust, or 100-km-thick lithosphere) and strain rate, using the thin sheet stress balance equations, calibrated with the subduction force balance analysis. Effective viscosities for the deforming lithosphere and/or crust are in the range 0.1–5 x 10 21 Pa s, with marked weakening in zones of high strain rate, and an abrupt transition to viscosities 〉10 22 Pa s at the margins of the rigid plates. If lateral variations in effective viscosity are only due to non-Newtonian behaviour, these data indicate a bulk power law rheology, with exponent n in the range 2–6. Average lithospheric or crustal deviatoric stresses 〈30 MPa. Such low driving stresses for the deforming crust are likely to be the result of a combination of pore fluid pressures much greater than hydrostatic (〉〉40 per cent lithostatic) and low coefficients of friction (〈〈0.6) on crustal faults.

Description: On 1580 April 6 one of the most destructive earthquakes of northwestern Europe took place in the Dover Strait ( Pas de Calais ). The epicentre of this seismic event, the magnitude of which is estimated to have been about 6.0, has been located in the offshore continuation of the North Artois shear zone, a major Variscan tectonic structure that traverses the Dover Strait. The location of this and two other moderate magnitude historical earthquakes in the Dover Strait suggests that the North Artois shear zone or some of its fault segments may be presently active. In order to investigate the possible fault activity in the epicentral area of the AD 1580 earthquake, we have gathered a large set of bathymetric and seismic-reflection data covering the almost-entire width of the Dover Strait. These data have revealed a broad structural zone comprising several subparallel WNW–ESE trending faults and folds, some of them significantly offsetting the Cretaceous bedrock. The geophysical investigation has also shown some indication of possible Quaternary fault activity. However, this activity only appears to have affected the lowermost layers of the sediment infilling Middle Pleistocene palaeobasins. This indicates that, if these faults have been active since Middle Pleistocene, their slip rates must have been very low. Hence, the AD 1580 earthquake appears to be a very infrequent event in the Dover Strait, representing a good example of the moderate magnitude earthquakes that sometimes occur in plate interiors on faults with unknown historical seismicity.

Description: The normal fault-system responsible of the 2009 M w 6.1 L'Aquila earthquake (Paganica-San Demetrio fault-system) comprises several narrow, fault-parallel valleys of controversial origin. We investigated a key section of the southeastern portion of this fault network along the small Verupola Valley. In order to characterize its nature and possible tectonic activity, we applied multiple-geosciences techniques able to image at depth the structure associated to this peculiar landform. We integrated magnetometry, 2-D P wave and resistivity tomography, surface waves and seismic noise analysis coupled with field mapping, shallow boreholes and trenching. According to our results, the Verupola Valley is a ~30–40-m-deep graben controlled by a SW-dipping master fault and synthetic splays paired with an antithetic NE-dipping fault. The SW-dipping splays are active and cut very shallow (〈2 m deep) Late Pleistocene sediments. The small amount of cumulated vertical offset (~15 m) across the conjugated system may indicate a young fault inception or very low Quaternary slip-rates. Due to its structural continuity with the adjacent mapped strands of the Paganica–San Demetrio fault network, we relate the Verupola Valley to the recent activity of the southeastern segment of this fault system. We also suggest that other fault-parallel valleys can have the same tectonic origin and setting of the Verupola Valley. This latter represents a scale-independent analogue from metric scale (exposed in the palaeoseismological trenches) to the Middle Aterno Basin scale (seen from seismic profiles and fault mapping). Overall, the imaged structural style is coherent with the regional tectonic setting due to Quaternary crustal extension.