ALBERT

Description: While the kinematics of Anatolia plate and the North Anatolian Fault System (NAFS) has been studied extensively, the slip rate and locking depth along the NAFS are usually assumed constant in the analyses due to the lack of sufficient data. This is also partly due to the reasonably good fit of Euler small circle and partly due to the lack of spatial resolution of observations to determine slip rates independently from locking depths. On the other hand, recent geodetic studies show a contrast for locking depth between Marmara and other parts of the NAFS, implying a non-uniform locking depth across the NAFS. In this study, we analyse new GPS data and homogenously combine available data sets covering the eastern part of the NAFS to form the most complete data set. In particular, we incorporate the first results of Turkish Real-Time Kinematic GPS Network (CORS-TR) into our data set. A detailed analysis of three profiles within the NAFS reveals an increase of locking depth in the middle profile to 19.1 ± 3.4 km from 11.9 ± 3.5 km in the easternmost profile while the slip rate is nearly constant (20–22 mm yr –1 ), which implies a variation of strain rate of ~100 nanostrain yr –1 . Assuming a constant locking depth throughout whole NAFS gives an average locking depth of 14.3 ± 1.7 km. Our best estimates of slip rates in block modelling which takes the variation of locking depths into account are in the range between 22.5 and 22.8 mm yr –1 over eastern part of the NAFS.

Description: The Taupo Volcanic Zone (TVZ) is one of the world's most productive regions of rhyolitic volcanism and contains the highly active Okataina Volcanic Centre (OVC). Within the TVZ, intra-arc extension is expressed as normal faulting within a zone known as the Taupo Rift. The OVC is located within a complex part of the rift, where volcanism and deformation is considered influenced by rift structure and kinematics. There has been significant research on the structural, volcanic and geophysical properties of the rift and OVC, but less focus on deformation using geodetic data. The limited studies that have utilized geodetic data do not clearly resolve the distribution of deformation and strain rates within the rift and OVC. This is essential to ensure that deformation signals from volcanic processes at the OVC are correctly identified and distinguished from those related to regional tectonic or local hydrothermal processes within the rift. In this paper, we present a picture of contemporary deformation at the OVC and within the surrounding rift in detail, using existing and new GPS campaign and continuous GPS (cGPS) data collected between 1998 and 2011. The results show a highly heterogeneous deformation and strain rate field (both extension and shortening) through the study area, partitioned into different parts of the rift. Our results agree well with earlier geodetic studies, as well as identify new features, but some deformation patterns conflict with long-term geological observations. In the OVC, we observe a locally rotated horizontal velocity field, significant vertical deformation and variable strain rates across the caldera. In the Tarawera Rift, we identify elevated extension and shear rates, which may have significant implications for volcanism there. A shortening pattern is identified through the central rift, which is unexpected in an intra-arc rifting environment. We attempt to explain the source/s of shortening and extension and discuss their implications for geodetic monitoring efforts in the OVC.

Description: The pre- and current collision of the Juan Fernández Ridge with the central Chilean margin at 31°–33°S is characterized by large-scale crustal thinning and long-term subsidence of the submarine forearc caused by subduction erosion processes. Here, we study the structure of the central Chilean margin in the ridge–trench collision zone by using wide-angle and multichannel seismic data. The transition from the upper to middle continental slope is defined by a trenchward dipping normal scarp with variable offsets of 500–2000 m height. Beneath the scarp, the 2-D velocity–depth models show a prominent lateral velocity contrast of 〉1 s –1 that propagates deep into the continental crust defining a major lateral seismic discontinuity. The discontinuity is interpreted as the lithological contact between the subsided/collapsed outermost forearc (composed of eroded and highly fractured volcanic rocks) and the seaward part of the uplifted Coastal Cordillera (made of less fractured metamorphic/igneous rocks). Extensional faults are abundant in the collapsed outermost forearc, however, landward of the continental slope scarp, both extensional and compressional structures are observed along the uplifted continental shelf that forms part of the Coastal Cordillera. Particularly, at the landward flank of the Valparaíso Forearc Basin (32°–33.5°S), shallow crustal seismicity has been recorded in 2008–2009 forming a dense cluster of thrust events of M w 4–5. The estimated hypocentres spatially correlate with the location of the fault scarp, and they highlight the upper part of the seismic crustal discontinuity.

Description: Geodetic observations of interseismic deformation in the Western United States provide constraints on microplate rotations, earthquake cycle processes, and slip partitioning across the Pacific–North America Plate boundary. These measurements may be interpreted using block models, in which the upper crust is divided into microplates bounded by faults that accumulate strain in a first-order approximation of earthquake cycle processes. The number and geometry of microplates are typically defined with boundaries representing a limited subset of the large number of potentially seismogenic faults. An alternative approach is to include a large number of potentially active faults bounding a dense array of microplates, and then algorithmically estimate the boundaries at which strain is localized. This approach is possible through the application of a total variation regularization (TVR) optimization algorithm, which simultaneously minimizes the L 2 norm of data residuals and the L 1 norm of the variation in the differential block motions. Applied to 3-D spherical block models, the TVR algorithm can be used to reduce the total variation between estimated rotation vectors, effectively grouping microplates that rotate together as larger blocks, and localizing fault slip on the boundaries of these larger block clusters. Here we develop a block model comprised of 137 microplates derived from published fault maps, and apply the TVR algorithm to identify the kinematically most important faults in the western United States. This approach reveals that of the 137 microplates considered, only 30 unique blocks are required to approximate deformation in the western United States at a residual level of 〈2 mm yr –1 .

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

Description: We explore thermal convection of a fluid with a temperature-dependent viscosity in a basally heated 3-D spherical shell using linear stability analyses and numerical experiments, while considering the application of our results to terrestrial planets. The inner to outer radius ratio of the shell f assumed in the linear stability analyses is in the range of 0.11–0.88. The critical Rayleigh number R c for the onset of thermal convection decreases by two orders of magnitude as f increases from 0.11 to 0.88, when the viscosity depends sensitively on the temperature, as is the case for real mantle materials. Numerical simulations carried out in the range of f  = 0.11–0.55 show that a thermal boundary layer (TBL) develops both along the surface and bottom boundaries to induce cold and hot plumes, respectively, when f is 0.33 or larger. However, for smaller f values, a TBL develops only on the bottom boundary. Convection occurs in the stagnant-lid regime where the root mean square velocity on the surface boundary is less than 1 per cent of its maximum at depth, when the ratio of the viscosity at the surface boundary to that at the bottom boundary exceeds a threshold that depends on f . The threshold decreases from 10 6.5 at f  = 0.11 to 10 4 at f  = 0.55. If the viscosity at the base of the convecting mantle is 10 20 –10 21  Pa s, the Rayleigh number exceeds R c for Mars, Venus and the Earth, but does not for the Moon and Mercury; convection is unlikely to occur in the latter planets unless the mantle viscosity is much lower than 10 20  Pa s and/or the mantle contains a strong internal heat source.

Description: Subducting oceanic lithosphere is an example of a thin sheet-like object whose characteristic lateral dimension greatly exceeds its thickness. Here we exploit this property to derive a new hybrid boundary-integral/thin sheet (BITS) representation of subduction that combines in a single equation all the forces acting on the sheet: gravity, internal resistance to bending and stretching, and the tractions exerted by the ambient mantle. For simplicity, we limit ourselves to 2-D. We solve the BITS equations using a discrete Lagrangian approach in which the sheet is represented by a set of vertices connected by edges. Instantaneous solutions for the sinking speed of a slab attached to a trailing flat sheet obey a scaling law of the form V / V Stokes  = fct(St), where V Stokes is a characteristic Stokes sinking speed and St is the sheet's flexural stiffness. Time-dependent solutions for the evolution of the sheet's shape and thickness show that these are controlled by the viscosity ratio between the sheet and its surroundings. An important advantage of the BITS approach is the possibility of generalizing the sheet's rheology, either to a viscosity that varies along the sheet or to a non-Newtonian shear-thinning rheology.

Description: The aftershock distribution of the 1933 Sanriku-oki outer trench earthquake is estimated by using modern relocation methods and a newly developed velocity structure to examine the spatial extent of the source-fault and the possibility of a triggered interplate seismicity. In this study, we first examined the regional data quality of the 1933 earthquake based on smoked-paper records and then relocated the earthquakes by using the 3-D velocity structure and double-difference method. The improvements of hypocentre locations using these methods were confirmed by the examination of recent earthquakes that are accurately located based on ocean bottom seismometer data. The results show that the 1933 aftershocks occurred under both the outer- and inner-trench-slope regions. In the outer-trench-slope region, aftershocks are distributed in a ~280-km-long area and their depths are shallower than 50 km. Although we could not constrain the fault geometry from the hypocentre distribution, the depth distribution suggests the whole lithosphere is probably not under deviatoric tension at the time of the 1933 earthquake. The occurrence of aftershocks under the inner trench slope was also confirmed by an investigation of waveform frequency difference between outer and inner trench earthquakes as recorded at Mizusawa. The earthquakes under the inner trench slope were shallow (depth less double equals30 km) and the waveforms show a low-frequency character similar to the waveforms of recent, precisely located earthquakes in the same area. They are also located where recent activity of interplate thrust earthquakes is high. These suggest that the 1933 outer-trench-slope main shock triggered interplate earthquakes, which is an unusual case in the order of occurrence in contrast with the more common pairing of a large initial interplate shock with subsequent outer-slope earthquakes. The off-trench earthquakes are distributed about 80 km width in the trench perpendicular direction. This wide width cannot be explained from a single high-angle fault confined at a shallow depth (depth less double equals50 km). The upward motion of the 1933 tsunami waveform records observed at Sanriku coast also cannot be explained from a single high-angle west-dipping normal fault. If we consider additional fault, involvement of high-angle, east-dipping normal faults can better explain the tsunami first motion and triggering of the aftershock in a wide area under the outer trench slope. Therefore multiple off-trench normal faults may have activated during the 1933 earthquake. We also relocated recent (2001–2012) seismicity by the same method. The results show that the present seismicity in the outer-trench-slope region can be divided into several groups along the trench. Comparison of the 1933 rupture dimensions based on our aftershock relocations with the morphologies of fault scarps in the outer trench slope suggest that the rupture was limited to the region where fault scarps are largely trench parallel and cross cut the seafloor spreading fabric. These findings imply that bending geometry and structural segmentation of the incoming plate largely controls the spatial extent of the 1933 seismogenic faulting. In this shallow rupture model for this largest outer trench earthquake, triggered seismicity in the forearc and structural control of faulting represent an important deformation styles for off-trench and shallow megathrust zones.

Description: Geodetic, geologic and palaeomagnetic data reveal that Oregon (western USA) rotates clockwise at 0.3 to 1.0° Ma –1 (relative to North America) about an axis near the Idaho–Oregon–Washington border, while northeast Washington is relatively fixed. This rotation has been going on for at least 15 Ma. The Yakima fold and thrust belt (YFTB) forms the boundary between northern Oregon and central Washington where convergence of the clockwise-rotating Oregon block is apparently accommodated. North–south shortening across the YFTB has been thought to occur in a fan-like manner, increasing in rate to the west. We obtained high-accuracy, high-density geodetic GPS measurements in 2012–2014 that are used with earlier GPS measurements from the 1990s to characterize YFTB kinematics. The new results show that the deformation associated with the YFTB starts at the Blue Mountains Anticline in northern Oregon and extends north beyond the Frenchman Hills in Washington, past the epicentre of the 1872 M w 7.0 Entiat earthquake to 49°N. The north–south strain rate across the region is 2 to 3 x 10 –9 yr –1 between the volcanic arc and the eastern edge of the YFTB (241.0°E); east of there it drops to about 10 –9 yr –1 . At the eastern boundary of the YFTB, faults and earthquake activity are truncated by a north-trending, narrow zone of deformation that runs along the Pasco Basin and Moses Lake regions near 240.9°E. This zone, abutting the Department of Energy Hanford Nuclear Reservation, accommodates about 0.5 mm yr –1 of east to northeast shortening. A similar zone of N-trending transpression is seen along 239.9°E where there is a change in the strike of the Yakima folds. The modern deformation of the YFTB is about 600 km wide from south to north and internally may be controlled by pre-existing crustal structure.

Description: The Central Anatolian orogenic plateau is represented by young volcanism, rapid plateau uplift and distinctive (past and active) tectonic deformation. In this study, we consider observational data in terms of regional present-day geodynamics in the region. The residual topography of Central Anatolia was derived to define the regional isostatic conditions according to Airy isostasy and infer the potential role of ‘dynamic topography’. 2-D thermomechanical forward models for coupled mantle-lithosphere flow/deformation were conducted along an N–S directional profile through the region (e.g. northern/Pontides, interior and southern/Taurides). These models were based on seismic tomography data that provide estimates about the present-day mantle thermal structure beneath the Anatolian plate. We compare the modelling results with calculated residual topography and independent data sets of geological deformation, gravity and high surface heat flow/widespread geothermal activity. Model results suggest that there is ~1 km of mantle flow induced dynamic topography associated with the sublithospheric flow driven by the seismically inferred mantle structure. The uprising mantle may have also driven the asthenospheric source of volcanism in the north (e.g. Galatia volcanic province) and the Cappadocia volcanic province in the south while elevating the surface in the last 10 Myr. Our dynamic topography calculations emphasize the role of vertical forcing under other orogenic plateaux underlain by relatively thin crust and low-density asthenospheric mantle.

Description: In many geological systems, inversion of density stratification sets in Rayleigh–Taylor (RT) instabilities, leading to an ascent of relatively low-density materials through the high-density overburden in the form of diapirs. These diapirs often originate from dipping low-density layers. This study aims to show how the initial tilt of such source layers can control the ascent behaviour of diapirs initiated by RT instabilities. Using two-layer viscous models we produced RT instabilities in physical experiments, and investigated the effects of source-layer tilts ( β ). Our experiments suggest that these diapirs ascend with contrasting lateral spreading rates in the up and down slope directions, resulting in their axi-asymmetric geometry. However, their heads retain a circular outline on the horizontal top surface, where the upwelling axis is located away from their geometric centre in the upslope direction. In this paper, we present a series of experimental models to demonstrate the spectrum of axi-symmetric to -asymmetric geometrical transitions with increasing β . Our experiments also reveal that when β is large (〉4°) the diapirs become unstable, resulting in a continuous migration of their stems in the upslope direction. Using the volume of fluid method we ran computational fluid dynamic (CFD) simulations to study the underlying hydrodynamics of axi-asymmetric diapiric growth. The CFD simulations show that β 〉 0° conditions develop stronger flow vortices on the downslope side of an ascending diapir, leading to a pressure difference between the up- and downslope flanks. Such a differential pressure causes the diapir head to spread at a faster rate in the tilt direction. An estimate of the asymmetric spreading rates is given as a function of β . Our present study provides a fundamental understanding of the hydrodynamic flow structure responsible for the asymmetric growth of RT instabilities on tilted source layers, as applicable to a wide range of large-scale geological settings, such as sedimentary basins and subduction zones.

Description: The geomechanical analysis of a highly compartmentalized reservoir is performed to simulate the seafloor subsidence due to gas production. The available observations over the hydrocarbon reservoir consist of bathymetric surveys carried out before and at the end of a 10-yr production life. The main goal is the calibration of the reservoir compressibility c M , that is, the main geomechanical parameter controlling the surface response. Two conceptual models are considered: in one (i) c M varies only with the depth and the vertical effective stress (heterogeneity due to lithostratigraphic variability); in another (ii) c M varies also in the horizontal plane, that is, it is spatially distributed within the reservoir stratigraphic units. The latter hypothesis accounts for a possible partitioning of the reservoir due to the presence of sealing faults and thrusts that suggests the idea of a block heterogeneous system with the number of reservoir blocks equal to the number of uncertain parameters. The method applied here relies on an ensemble-based data assimilation (DA) algorithm (i.e. the ensemble smoother, ES), which incorporates the information from the bathymetric measurements into the geomechanical model response to infer and reduce the uncertainty of the parameter c M . The outcome from conceptual model (i) indicates that DA is effective in reducing the c M uncertainty. However, the maximum settlement still remains underestimated, while the areal extent of the subsidence bowl is overestimated. We demonstrate that the selection of the heterogeneous conceptual model (ii) allows to reproduce much better the observations thus removing a clear bias of the model structure. DA allows significantly reducing the c M uncertainty in the five blocks (out of the seven) characterized by large volume and large pressure decline. Conversely, the assimilation of land displacements only partially constrains the prior c M uncertainty in the reservoir blocks marginally contributing to the cumulative seafloor subsidence, that is, blocks with low pressure.

Description: This work aims to explore the ongoing tectonic activity of structures in the outermost sector of the Northern Apennines, which represents the active leading edge of the thrust belt and is dominated by compressive deformation. We have applied the Persistent Scatterer Interferometry (PSI) technique to obtain new insights into the present-day deformation pattern of the frontal area of the Northern Apennine. PSI has proved to be effective in detecting surface deformation of wide regions involved in low tectonic movements. We used 34 Envisat images in descending geometry over the period of time between 2004 and 2010, performing about 300 interferometric pairs. The analysis of the velocity maps and of the PSI time-series has allowed to observe ground deformation over the sector of the Po Plain between Piacenza and Reggio Emilia. The time-series of permanent GPS stations located in the study area, validated the results of the PSI technique, showing a good correlation with the PS time-series. The PS analysis reveals the occurrence of a well-known subsidence area on the rear of the Ferrara arc, mostly connected to the exploitation of water resources. In some instances, the PS velocity pattern reveals ground uplift (with mean velocities ranging from 1 to 2.8 mm yr –1 ) above active thrust-related anticlines of the Emilia and Ferrara folds, and part of the Pede-Apennine margin. We hypothesize a correlation between the observed uplift deformation pattern and the growth of the thrust-related anticlines. As the uplift pattern corresponds to known geological features, it can be used to constrain the seismo-tectonic setting, and a working hypothesis may involve that the active Emilia and Ferrara thrust folds would be characterized by interseismic periods possibly dominated by aseismic creep.

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

Description: An important real world application of doublet flow occurs in well design of both geothermal and hydrocarbon reservoirs. A guiding principle for fluid management of injection and extraction wells is that mass balance is commonly assumed between the injected and produced fluid. Because the doublets are considered closed loops, the injection fluid is assumed to eventually reach the producer well and all the produced fluid ideally comes from stream tubes connected to the injector of the well pair making up the doublet. We show that when an aquifer background flow occurs, doublets will rarely retain closed loops of fluid recirculation. When the far-field flow rate increases relative to the doublet's strength, the area occupied by the doublet will diminish and eventually vanishes. Alternatively, rather than using a single injector (source) and single producer (sink), a linear array of multiple injectors separated by some distance from a parallel array of producers can be used in geothermal energy projects as well as in waterflooding of hydrocarbon reservoirs. Fluid flow in such an arrangement of parallel source-sink arrays is shown to be macroscopically equivalent to that of a line doublet. Again, any far-field flow that is strong enough will breach through the line doublet, which then splits into two vortices. Apart from fundamental insight into elementary flow dynamics, our new results provide practical clues that may contribute to improve the planning and design of doublets and direct line drives commonly used for flow management of groundwater, geothermal and hydrocarbon reservoirs.

Description: We estimate fluid sources around a subducted seamount along the northern Hikurangi subduction margin of New Zealand, using thermomechanical numerical modelling informed by wedge structure and porosities from multichannel seismic data. Calculated fluid sources are input into an independent fluid-flow model to explore the key controls on overpressure generation to depths of 12 km. In the thermomechanical models, sediment transport through and beneath the wedge is calculated assuming a pressure-sensitive frictional rheology. The change in porosity, pressure and temperature with calculated rock advection is used to compute fluid release from compaction and dehydration. Our calculations yield more precise information about source locations in time and space than previous averaged estimates for the Hikurangi margin. The volume of fluid release in the wedge is smaller than previously estimated from margin-averaged calculations (~14 m 3  yr –1  m –1 ), and is exceeded by fluid release from underlying (subducting) sediment (~16 m 3  yr –1  m –1 ). Clay dehydration contributes only a small quantity of fluid by volume (~2 m 3  yr –1  m –1 from subducted sediment), but the integrated effect is still significant landward of the seamount. Fluid source terms are used to estimate fluid pressures around a subducting seamount in the fluid-flow models, using subducted sediment permeability derived from porosity, and testing two end-members for décollement permeability. Models in which the décollement acts as a fluid conduit predict only moderate fluid overpressure in the wedge and subducting sediment. However, if the subduction interface becomes impermeable with depth, significant fluid overpressure develops in subducting sediment landward of the seamount. The location of predicted fluid overpressure and associated dehydration reactions is consistent with the idea that short duration, shallow, slow slip events (SSEs) landward of the seamount are caused by anomalous fluid pressures; alternatively, it may result from frictional effects of changing clay content along the subduction interface.

Description: A method is outlined by means of which it is possible to estimate high-resolution vertical displacements due to an earthquake even in the case where high-resolution topography is lacking before the earthquake. This result can be achieved by combining a highly accurate, post-event digital elevation model (DEM), for example lidar, with archived satellite imagery. The method is illustrated by calculating vertical displacements for the 2010 El Mayor-Cucapah earthquake. For this earthquake, there are both pre- and post-event lidar DEMs from which vertical displacements may also be estimated after correcting for the lateral advection of topography due to horizontal displacements. A comparison between the two means of deriving vertical displacements shows generally good agreement, with the displacements obtained using satellite imagery performing better in high relief areas. As a result of this property, we are able to trace the vertical offsets due to the El Mayor-Cucapah earthquake as the rupture jumped from the Pescadores fault to the Borrego fault in propagating through the high relief of the Sierra Cucapah.

Description: We present high-resolution tomographic images in source areas of 26 large crustal earthquakes ( M 6.0–7.2) which occurred in Northeast Japan (Tohoku) during the past 120 yr from 1894 to 2014. Prominent low-velocity (low- V ) and high Poisson's ratio (high- ) anomalies are revealed in the crust and mantle wedge under the source areas. Beneath the volcanic front and backarc areas, the low- V and high- zones reflect arc-magma related high-temperature anomalies which are produced by joint effects of corner flow in the mantle wedge and fluids from dehydration of the subducting Pacific slab. The hot anomalies cause locally thinning and weakening of the brittle seismogenic layer above them. Low-frequency micro-earthquakes are observed in the lower crust and uppermost mantle in or around the low- V zones, which reflect ascending of arc magma and fluids from the mantle wedge to the crust. No volcano and magma exist in the forearc area due to low temperature there, hence the low- V zones in the forearc reflect fluids from the slab dehydration. The ascending fluids may have produced a ‘water wall’ in the mantle wedge and crust beneath the forearc area. When the water enters active faults in the crust, the fault-zone friction is reduced and so large earthquakes can be induced. These results indicate that the nucleation of a large earthquake is not entirely a mechanical process, but is closely associated with subduction dynamics and physical and chemical properties of rocks in the crust and upper mantle. In particular, arc magma and fluids play an important role in the seismogenesis.

Description: We present observations and models of the Sulaiman Range of western Pakistan that shed new light on the evolution and deformation of fold-thrust belts. Earthquake source inversions show that the seismic deformation in the range is concentrated in the thick pile of sediments overlying the underthrusting lithosphere of the Indian subcontinent. The slip vectors of the earthquakes vary in strike around the margin of the range, in tandem with the shape of the topography, suggesting that gravitational driving forces arising from the topography play an important role in governing the deformation of the region. Numerical models suggest that the active deformation, and the extreme plan-view curvature of the range, are governed by the presence of weak sediments in a pre-existing basin on the underthrusting Indian Plate. These sediments affect the stress-state in the over-riding mountain range and allow for the rapid propagation of the nose of the range and the development of extreme curvature and laterally varying surface gradients.

Description: The mechanical damage characteristics of sandstone subjected to cyclic loading is very significant to evaluate the stability and safety of deep excavation damage zones. However to date, there are very few triaxial experimental studies of sandstone under cyclic loading. Moreover, few X-ray micro-computed tomography (micro-CT) observations have been adopted to reveal the damage mechanism of sandstone under triaxial cyclic loading. Therefore, in this research, a series of triaxial cyclic loading tests and X-ray micro-CT observations were conducted to analyse the mechanical damage characteristics of sandstone with respect to different confining pressures. The results indicated that at lower confining pressures, the triaxial strength of sandstone specimens under cyclic loading is higher than that under monotonic loading; whereas at confining pressures above 20 MPa, the triaxial strength of sandstone under cyclic loading is approximately equal to that under monotonic loading. With the increase of cycle number, the crack damage threshold of sandstone first increases, and then significantly decreases and finally remains constant. Based on the damage evolution of irreversible deformation, it appears that the axial damage value of sandstone is all higher than the radial damage value before the peak strength; whereas the radial damage value is higher than the axial damage value after the peak strength. The evolution of Young's modulus and Poisson's ratio of sandstone can be characterized as having four stages: (i) Stage I: material strengthening; (ii) Stage II: material degradation; (iii) Stage III: material failure and (iv) Stage IV: structure slippage. X-ray micro-CT observations demonstrated that the CT scanning surface images of sandstone specimens are consistent with actual surface crack photographs. The analysis of the cross-sections of sandstone supports that the system of crack planes under triaxial cyclic loading is much more complicated than that under triaxial monotonic loading. More axial and lateral tensile cracks were observed in the specimens under cyclic loading than under monotonic loading.

Description: We determine the deep structure of the eastern Algerian basin and its southern margin in the Annaba region (easternmost Algeria), to better constrain the plate kinematic reconstruction in this region. This study is based on new geophysical data collected during the SPIRAL cruise in 2009, which included a wide-angle, 240-km-long, onshore–offshore seismic profile, multichannel seismic reflection lines and gravity and magnetic data, complemented by the available geophysical data for the study area. The analysis and modelling of the wide-angle seismic data including refracted and reflected arrival travel times, and integrated with the multichannel seismic reflection lines, reveal the detailed structure of an ocean-to-continent transition. In the deep basin, there is an ~5.5-km-thick oceanic crust that is composed of two layers. The upper layer of the crust is defined by a high velocity gradient and P -wave velocities between 4.8 and 6.0 km s –1 , from the top to the bottom. The lower crust is defined by a lower velocity gradient and P -wave velocity between 6.0 and 7.1 km s –1 . The Poisson ratio in the lower crust deduced from S -wave modelling is 0.28, which indicates that the lower crust is composed mainly of gabbros. Below the continental edge, a typical continental crust with P -wave velocities between 5.2 and 7.0 km s –1 , from the top to the bottom, shows a gradual seaward thinning of ~15 km over an ~35-km distance. This thinning is regularly distributed between the upper and lower crusts, and it characterizes a rifted margin, which has resulted from backarc extension at the rear of the Kabylian block, here represented by the Edough Massif at the shoreline. Above the continental basement, an ~2-km-thick, pre-Messinian sediment layer with a complex internal structure is interpreted as allochthonous nappes of flysch backthrusted on the margin during the collision of Kabylia with the African margin. The crustal structure, moreover, provides evidence for Miocene emplacement of magmatic intrusions in both the deep basin and the continental margin. Based on the crustal structure, we propose that the eastern Algerian basin opened during the southeastward migration of the European forearc before the collision, along a NW–SE elongated spreading centre that ran perpendicular to the subduction trend. Such an atypical geometry is explained by the diverging directions of the subduction rollback during the backarc opening: eastward for the Corsica–Sardinia block, and southward for the Kabylian blocks. This geometry of the forearc can be interpreted as the surface expression of a slab tear at depth, which is responsible for atypical magmatism in the overlying backarc oceanic basin.

Description: In November 2010, intense seismic activity including 29 events with a magnitude above 5.0, started in the western part of the Gulf of Aden, where the structure of the oceanic spreading ridge is characterized by a series of N115°-trending slow-spreading segments set within an EW-trending rift. Using signals recorded by permanent and temporary networks in Djibouti and Yemen, we located 1122 earthquakes, with a magnitude ranging from 2.1 to 5.6 from 2010 November 1 to 2011 March 31. By looking in detail at the space–time distribution of the overall seismicity, and both the frequency and the moment tensor of large earthquakes, we re-examine the chronology of this episode. In addition, we also interpret the origin of the activity using high-resolution bathymetric data, as well as from observations of seafloor cable damage caused by high temperatures and lava flows. The analysis allows us to identify distinct active areas. First, we interpret that this episode is mainly related to a diking event along a specific ridge segment, located at E044°. In light of previous diking episodes in nearby subaerial rift segments, for which field constraints and both seismic and geodetic data exist, we interpret the space–time evolution of the seismicity of the first few days. Migration of earthquakes suggests initial magma ascent below the segment centre. This is followed by a southeastward dike propagation below the rift immediately followed by a northwestward dike propagation below the rift ending below the northern ridge wall. The cumulative seismic moment associated with this sequence reaches 9.1 x 10 17 Nm, and taking into account a very low seismic versus geodetic moment, we estimate a horizontal opening of ~0.58–2.9 m. The seismic activity that followed occurred through several bursts of earthquakes aligned along the segment axis, which are interpreted as short dike intrusions implying fast replenishment of the crustal magma reservoir feeding the dikes. Over the whole period, the opening is estimated to be ~1.76–8.8 m across the segment. A striking feature of this episode is that the seismicity remained confined within one individual segment, whereas the adjacent en-echelon segments were totally quiescent, suggesting that the magma supply system of one segment is disconnected from those of the neighbouring segments. Second, we identify activity induced by the first intrusion with epicentres aligned along an N035°E-trending, ~30 km long at the northwestern end of the active opening segment. This group encompasses more than seven earthquakes with magnitude larger than 5.0, and with strike-slip focal mechanisms consistent with the faults identified in the bathymetry and the structural pattern of the area. We propose that a transform fault is currently in formation which indicates an early stage of the ridge segmentation, at the locus of the trend change of the spreading ridge, which also corresponds to the boundary between a clear oceanic lithosphere and the zone of transform between continental and oceanic crust.

Description: We challenge the perspective that seismicity could contribute to polar motion by arguing quantitatively that, in first approximation and on the average, interseismic deformations can compensate for it. This point is important because what we must simulate and observe in Earth Orientation Parameter time-series over intermediate timescales of decades or centuries is the residual polar motion resulting from the two opposing processes of coseismic and interseismic deformations. In this framework, we first simulate the polar motion caused by only coseismic deformations during the longest period available of instrumental seismicity, from 1900 to present, using both the CMT and ISC-GEM catalogues. The instrumental seismicity covering a little longer than one century does not represent yet the average seismicity that we should expect on the long term. Indeed, although the simulation shows a tendency to move the Earth rotation pole towards 133°E at the average rate of 16.5 mm yr –1 , this trend is still sensitive to individual megathrust earthquakes, particularly to the 1960 Chile and 1964 Alaska earthquakes. In order to further investigate this issue, we develop a global seismicity model (GSM) that is independent from any earthquake catalogue and that describes the average seismicity along plate boundaries on the long term by combining information about present-day plate kinematics with the Anderson theory of faulting, the seismic moment conservation principle and a few other assumptions. Within this framework, we obtain a secular polar motion of 8 mm yr –1 towards 112.5°E that is comparable with that estimated from 1900 to present using the earthquake catalogues, although smaller by a factor of 2 in amplitude and different by 20° in direction. Afterwards, in order to reconcile the idea of a secular polar motion caused by earthquakes with our simplest understanding of the seismic cycle, we adapt the GSM in order to account for interseismic deformations and we use it to quantify, for the first time ever, their contribution to polar motion. Taken together, coseismic and interseismic deformations make the rotation pole wander around the north pole with maximum polar excursions of about 1 m. In particular, the rotation pole moves towards about Newfoundland when the interseismic contribution dominates over the coseismic ones (i.e. during phases of low seismicity or, equivalently, when most of the fault system associated with plate boundaries is locked). When megathrust earthquakes occur, instead, the rotation pole is suddenly shifted in an almost opposite direction, towards about 133°E.

Description: While elasticity is a defining characteristic of the Earth's lithosphere, it is often ignored in numerical models of long-term tectonic processes in favour of a simpler viscoplastic description. Here we assess the consequences of this assumption on a well-studied geodynamic problem: the growth of normal faults at an extensional plate boundary. We conduct 2-D numerical simulations of extension in elastoplastic and viscoplastic layers using a finite difference, particle-in-cell numerical approach. Our models simulate a range of faulted layer thicknesses and extension rates, allowing us to quantify the role of elasticity on three key observables: fault-induced topography, fault rotation, and fault life span. In agreement with earlier studies, simulations carried out in elastoplastic layers produce rate-independent lithospheric flexure accompanied by rapid fault rotation and an inverse relationship between fault life span and faulted layer thickness. By contrast, models carried out with a viscoplastic lithosphere produce results that may qualitatively resemble the elastoplastic case, but depend strongly on the product of extension rate and layer viscosity U x L . When this product is high, fault growth initially generates little deformation of the footwall and hanging wall blocks, resulting in unrealistic, rigid block-offset in topography across the fault. This configuration progressively transitions into a regime where topographic decay associated with flexure is fully accommodated within the numerical domain. In addition, high U x L favours the sequential growth of multiple short-offset faults as opposed to a large-offset detachment. We interpret these results by comparing them to an analytical model for the fault-induced flexure of a thin viscous plate. The key to understanding the viscoplastic model results lies in the rate-dependence of the flexural wavelength of a viscous plate, and the strain rate dependence of the force increase associated with footwall and hanging wall bending. This behaviour produces unrealistic deformation patterns that can hinder the geological relevance of long-term rifting models that assume a viscoplastic rheology.

Description: Steady plate subduction elastically brings about permanent lithospheric deformation in island arcs, though this effect has been neglected in most studies based on elastic dislocation theory. We investigate the characteristics of the permanent lithospheric deformation using a kinematic model, in which steady slip motion is given along a plate interface in the elastic lithosphere overlying the viscoelastic asthenosphere under gravity. As a rule of thumb, long-term lithospheric deformation can be understood as a bending of an elastic plate floating on non-viscous fluid, because the asthenosphere behaves like water on the long term. The steady slip below the lithosphere–asthenosphere boundary does not contribute to long-term lithospheric deformation. Hence, the key parameters that control the lithospheric deformation are only the thickness of the lithosphere and the geometry of the plate interface. Slip on a plate interface generally causes substantial vertical displacement, and gravity always tries to retrieve the original gravitational equilibrium. For a curved plate interface gravity causes convex upward bending of the island-arc lithosphere, while for a planar plate interface gravity causes convex downward bending. Larger curvature and thicker lithosphere generally results in larger deformation. When the curvature changes along the plate interface, internal deformation is also involved intrinsically, which modifies the deformation field due to gravity. Because the plate interface generally has some curvature, at least near the trench, convex upward bending of the island-arc lithosphere, which involves uplift of island-arc and subsidence around the trench, is always realized. On the other hand, the deformation field of the island-arc lithosphere sensitively depends on lithospheric thickness and plate interface geometry. These characteristics obtained by the numerical simulation are consistent with observed topography and free-air gravity anomalies in subduction zones: a pair of topography and gravity anomalies, high in the arc and low around the trench, is observed without exceptions all over the world, while there are large variety in the amplitude and horizontal scale of the topography and gravity anomalies.

Description: We investigate the influence of the glacial isostatic adjustment (GIA) on the deformation at the surface and at seismogenic depths in Fennoscandia. The surface strain rate field, derived from geodetic data, is controlled by GIA which causes NW–SE extension of up to 4  x  10 –9  yr –1 in most of mainland Fennoscandia, surrounded by regions of radial shortening towards the centre of uplift. The seismic deformation field, derived from a new compilation of focal mechanisms, shows consistent NW–SE compression on the Norwegian continental margin and a tendency towards tension in mainland Fennoscandia. The seismic moment rate is at least two orders of magnitude smaller than the geodetic moment rate. We propose that the low level of seismicity and the tendency towards tensional focal mechanisms in mainland Fennoscandia may be explained by the destructive interference of the regional stress from ridge push with the flexural stress due to GIA.

Description: New Guinea is a region characterized by rapid oblique convergence between the Pacific and Australian tectonic plates. The detailed tectonics of the region, including the partitioning of relative block motions and fault slip rates within this complex boundary plate boundary zone are still not well understood. In this study, we quantify the distribution of the deformation throughout the central and western parts of Papua New Guinea (PNG) using 20 yr of GPS data (1993–2014). We use an elastic block model to invert the regional GPS velocities as well as earthquake slip vectors for the location and rotation rates of microplate Euler poles as well as fault slip parameters in the region. Convergence between the Pacific and the Australian plates is accommodated in northwestern PNG largely by the New Guinea Trench with rates exceeding 90 mm yr –1 , indicating that this is the major active interplate boundary. However, some convergent deformation is partitioned into a shear component with ~12 per cent accommodated by the Bewani-Torricelli fault zone and the southern Highlands Fold-and-Thrust Belt. New GPS velocities in the eastern New Guinea Highlands region have led to the identification of the New Guinea Highlands and the Papuan Peninsula being distinctly different blocks, separated by a boundary through the Aure Fold-and-Thrust Belt complex which accommodates an estimated 4–5 mm yr –1 of left-lateral and 2–3 mm yr –1 of convergent motion. This implies that the Highlands Block is rotating in a clockwise direction relative to the rigid Australian Plate, consistent with the observed transition to left-lateral strike-slip regime observed in western New Guinea Highlands. We find a better fit of our block model to the observed velocities when assigning the current active boundary between the Papuan Peninsula and the South Bismark Block to be to the north of the city of Lae on the Gain Thrust, rather than on the more southerly Ramu-Markham fault as previously thought. This may indicate a temporary shift of activity onto out of sequence thrusts like the Gain Thrust as opposed to the main frontal thrust (the Ramu-Markham fault). In addition, we show that the southern Highlands Fold-and-Thrust Belt is the major boundary between the rigid Australian Plate and the New Guinea Highlands Block, with convergence occurring at rates between ~6 and 13 mm yr –1 .

Description: We examine the geodetically derived rotational variations for the rate of change of degree-two harmonics of Earth's geopotential, $\skew5\dot J_2$ , and true polar wander, combining a recent melting model of glaciers and the Greenland and Antarctic ice sheets taken from the IPCC 2013 Report (AR5) with two representative GIA ice models describing the last deglaciation, ICE5G and the ANU model developed at the Australian National University. Geodetically derived observations of $\skew4\dot J_2$ are characterized by temporal changes of –(3.7 ± 0.1)  x 10 –11 yr –1 for the period 1976–1990 and –(0.3 ± 0.1)  x 10 –11 yr –1 after ~2000. The AR5 results make it possible to evaluate the recent melting of the major ice sheets and glaciers for three periods, 1900–1990, 1991–2001 and after 2002. The observed $\skew4\dot J_2$ and the component of $\skew4\dot J_2$ due to recent melting for different periods indicate a long-term change in $\skew4\dot J_2$ —attributed to the Earth's response to the last glacial cycle—of –(6.0–6.5)  x 10 –11 yr –1 , significantly different from the values adopted to infer the viscosity structure of the mantle in most previous studies. This is a main conclusion of this study. We next compare this estimate with the values of $\skew4\dot J_2$ predicted by GIA ice models to infer the viscosity structure of the mantle, and consequently obtain two permissible solutions for the lower mantle viscosity ( lm ), ~10 22 and (5–10)  x 10 22 Pa s, for both adopted ice models. These two solutions are largely insensitive to the lithospheric thickness and upper mantle viscosity as indicated by previous studies and relatively insensitive to the viscosity structure of the D '' layer. The ESL contributions from the Antarctic ice sheet since the last glacial maximum (LGM) for ICE5G and ANU are about 20 and 30 m, respectively, but glaciological reconstructions of the Antarctic LGM ice sheet have suggested that its ESL contribution may have been less than ~10 m. The GIA-induced $\skew4\dot J_2$ for GIA ice models with an Antarctic ESL component of ~10 m suggests two permissible lower mantle viscosity solutions of lm ~ 2  x 10 22 and ~5  x 10 22 Pa s or one solution with (2–5)  x 10 22 Pa s. These results suggest that the effective lower mantle viscosity is larger than ~10 22 Pa s regardless of the uncertainties for an Antarctic ESL component. We also examine the polar wander due to recent melting and GIA processes, suggesting that the observed polar wander may be significantly attributed to convection motions in the mantle and/or another cause, particularly for permissible lower mantle viscosity solution of (5–10)  x 10 22 Pa s.

Description: In this work, we test the analysis of anisotropy of magnetic susceptibility (AMS) for the study of diapirs in the Naval salt wall that is located in the western end of the South Pyrenean Central Unit. Two hundred eighty-six specimens from 30 sites were sampled in the mobilized Middle-Upper Triassic evaporites and shales that constitute the core of this diapir. Resulting magnetic ellipsoids are mostly oblate and scalar parameters are similar to those found in weakly deformed, sedimentary rocks. From the observation of samples and thin sections, magnetic fabrics can be related to petrofabric: magnetic foliation corresponds to flattening or shear planes and magnetic lineation to flow/shear direction, elongation of grains or tight microfolds. The magnetic foliation shows a roughly circular pattern with shallow dips and some deflections due to Pyrenean structures. Magnetic lineations show a slightly dominant NE–SW trend and shallow plunges. Magnetic lineations and foliations are consistent with the theoretical pattern expected at the roof of a vertically grown salt wall, thus supporting that magnetic fabrics registered evaporitic flow at late stages of diapiric growth. We demonstrate that AMS is an efficient and promising tool to define diapiric flow where tectonic fabric cannot be directly measured. Knowledge of the geometry and internal deformation of the diapiric body helps to minimize ambiguities in the interpretation of magnetic fabrics in terms of internal deformation.

Description: Several large earthquakes in the Hellenic subduction zone have been documented in historical records from around the eastern Mediterranean, but the relative seismic quiescence of the region over the period of instrumental observation means that the exact locations of these earthquakes and their tectonic significance are not known. We present AMS radiocarbon dates from uplifted late Holocene palæoshorelines from the island of Rhodes, showing that uplift is most consistent with a single large ( M W  ≥ 7.7) reverse-faulting earthquake between about 2000 BC and 200 BC. Analysis of the uplift treating the earthquake as a dislocation in an elastic half-space shows a predominantly reverse-faulting event with a slip vector oblique to the direction of convergence between Rhodes and Nubia. We suggest that the fault responsible for the uplift dips at an angle of 30–60° above the more gently dipping oblique subduction interface. The highly oblique convergence across the eastern Hellenic plate boundary zone appears to be partitioned into reverse slip on faults that strike parallel to the boundary and strike-parallel or oblique slip on the subduction interface. Hydrodynamical simulation of tsunami propagation from a range of tectonically plausible sources suggests that earthquakes on the fault uplifting Rhodes represent a significant tsunami hazard for Rhodes and SW Turkey, and also possibly for Cyprus and the Nile Delta.

Description: Reconstructions of the history of convergence between the Nubia and Eurasia plates constitute an important part of a broader framework for understanding deformation in the Mediterranean region and the closing of the Mediterranean Basin. Herein, we combine high-resolution reconstructions of Eurasia-North America and Nubia-North America Plate motions to determine rotations that describe Nubia-Eurasia Plate motion at ~1 Myr intervals for the past 20 Myr. We apply trans-dimensional hierarchical Bayesian inference to the Eurasia-North America and Nubia-North America rotation sequences in order to reduce noise in the newly estimated Nubia-Eurasia rotations. The noise-reduced rotation sequences for the Eurasia-North America and Nubia-North America Plate pairs describe remarkably similar kinematic histories since 20 Ma, consisting of relatively steady seafloor spreading from 20 to 8 Ma, ~20 per cent opening-rate slowdowns at 8–6.5 Ma, and steady plate motion from ~7 Ma to the present. Our newly estimated Nubia-Eurasia rotations predict that convergence across the central Mediterranean Sea slowed by ~50 per cent and rotated anticlockwise after ~25 Ma until 13 Ma. Motion since 13 Ma has remained relatively steady. An absence of evidence for a significant change in motion immediately before or during the Messinian Salinity Crisis at 6.3–5.6 Ma argues against a change in plate motion as its causative factor. The detachment of the Arabian Peninsula from Africa at 30–24 Ma may have triggered the convergence rate slowdown before 13 Ma; however, published reconstructions of Nubia-Eurasia motion for times before 20 Ma are too widely spaced to determine with confidence whether the two are correlated. A significant discrepancy between our new estimates of Nubia-Eurasia motion during the past few Myr and geodetic estimates calls for further investigation.

Description: The Yellowstone–East Snake River Plain hotspot track has been intensely studied since several decades and is widely considered to result from the interaction of a mantle plume with the North American plate. An integrated conclusive geodynamic interpretation of this extensive data set is however presently still lacking, and our knowledge of the dynamical processes beneath Yellowstone is patchy. It has been argued that the Yellowstone plume has delaminated the lower part of the thick Wyoming cratonic lithosphere. We derive an original dynamic model to quantify delamination processes related to mantle plume–lithosphere interactions. We show that fast (~300 ka) lithospheric delamination is consistent with the observed timing of formation of successive volcanic centres along the Yellowstone hotspot track and requires (i) a tensile stress regime within the whole lithosphere exceeding its failure threshold, (ii) a purely plastic rheology in the lithosphere when stresses reach this yield limit, (iii) a dense lower part of the 200 km thick Wyoming lithosphere and (iv) a decoupling melt horizon inside the median part of the lithosphere. We demonstrate that all these conditions are verified and that ~150 km large and ~100 km thick lithospheric blocks delaminate within 300 ka when the Yellowstone plume ponded below the 200 km thick Wyoming cratonic lithosphere. Furthermore, we take advantage of the available extensive regional geophysical and geological observation data sets to design a numerical 3-D upper-mantle convective model. We propose a map of the ascending convective sheets contouring the Yellowstone plume. The model further evidences the development of a counter-flow within the lower part of the lithosphere centred just above the Yellowstone mantle plume axis. This counter-flow controls the local lithospheric stress field, and as a result the trajectories of feeder dykes linking the partial melting source within the core of the mantle plume with the crust by crosscutting the lithospheric mantle. This counter-flow further explains the 50 km NE shift observed between the mantle plume axis and the present-day Yellowstone Caldera as well as the peculiar shaped crustal magma chambers.

Description: We present source models for the August 2014 Murmuri (Dehloran) earthquake sequence in the Zagros Mountains of Iran. An M w 6.2 mainshock was followed by an aftershock sequence containing five events of M w ≥ 5.4. Models of P and SH waveforms show that all events had dominantly thrust-faulting mechanisms, and had centroid depths that place them within the thick sedimentary sequence, above the crystalline basement. The combination of our estimated focal mechanisms, relative relocations of the event hypocentres and the surface displacement patterns observed using InSAR imply that the mainshock and largest aftershock ruptured different fault planes and both contributed to the surface deformation. The fault planes both slipped in horizontally elongated patches, possibly due to rheological layering limiting the updip and downdip extent of rupture. The slip vector of the Murmuri mainshock implies that the decollement beneath the Lorestan Arc is weaker than any such feature beneath the Dezful Embayment, providing an explanation for the plan-view sinuosity of the range-front of the Zagros Mountains.

Description: Displacements and stress-field changes associated with earthquakes, volcanoes, landslides and human activity are often simulated using numerical models in an attempt to understand the underlying processes and their governing physics. The application of elastic dislocation theory to these problems, however, may be biased because of numerical instabilities in the calculations. Here, we present a new method that is free of artefact singularities and numerical instabilities in analytical solutions for triangular dislocations (TDs) in both full-space and half-space. We apply the method to both the displacement and the stress fields. The entire 3-D Euclidean space $\mathbb {R}^{3}$ is divided into two complementary subspaces, in the sense that in each one, a particular analytical formulation fulfils the requirements for the ideal, artefact-free solution for a TD. The primary advantage of the presented method is that the development of our solutions involves neither numerical approximations nor series expansion methods. As a result, the final outputs are independent of the scale of the input parameters, including the size and position of the dislocation as well as its corresponding slip vector components. Our solutions are therefore well suited for application at various scales in geoscience, physics and engineering. We validate the solutions through comparison to other well-known analytical methods and provide the MATLAB codes.

Description: Using dynamic rupture models of a right-lateral fault embedded in an elastic or elastoplastic 3-D medium, we investigate elastic and inelastic responses of compliant fault zones to nearby earthquake ruptures. We particularly examine effects of fault zone depth, width, shape and rigidity reduction on the surface displacement field. Our results from elastic models show that deeper and wider fault zones generally result in larger residual displacements. However, for shallow fault zones, the vertical residual displacement is insensitive to or even decreases with fault zone width. The width of horizontal displacement anomalies across a fault zone is only indicative of the fault zone width near the Earth's surface. There are trade-off effects among fault zone depth, width, shape and rigidity reduction on the amplitude of surface residual displacements. Our elastoplastic models show that plastic strain can occur along the entire fault zone near the Earth's surface and in the extensional quadrant at depth, if fault zone rocks are close to failure before a nearby earthquake happens. Compared with results from elastic models, plastic strain near the Earth's surface generally enhances surface displacements of a fault zone and does not change the trend of effects of fault zone depth and width, while plastic strain at depth can result in reduced retrograde motion or sympathetic motion across the fault zone, and introduce complexities in effects of fault zone depth and width. Sympathetic horizontal motion more likely occurs across a narrow fault zone with inelastic response at depth. Vertical motion in the extensional quadrant may actually decrease with fault zone width in elastoplastic models. Sympathetic horizontal motion, or small retrograde horizontal motion in conjunction with large vertical motion above a fault zone is indicative of inelastic response of a fault zone at depth.

Description: Models of the glacial isostatic adjustment process, which is dominated by the influence of the Late Pleistocene cycle of glaciation and deglaciation, depend on two fundamental inputs: a history of ice-sheet loading and a model of the radial variation of mantle viscosity. These models may be tested and refined by comparing their local predictions of relative sea level history to geological inferences based upon appropriate sea level indicators. The U.S. Atlantic coast is a region of particular interest in this regard, due to the fact that data from the length of this coast provides a transect of the forebulge associated with the former Laurentide ice sheet. High-quality relative sea level histories from this region are employed herein to explore the ability of current models of mantle viscosity to explain the inferred evolution of relative sea level that have accompanied forebulge collapse following deglaciation. Existing misfits are characterized, and alternatives are explored for their reconciliation. It is demonstrated that a new model of mantle viscosity, referred to herein as VM6, when coupled with the latest model of deglaciation history ICE-6G_C, is able to eliminate the majority of these misfits, while continuing to reconcile a wide range of other important geophysical observables, as well as additional relative sea level data from the North American. West coast which also record the collapse of the forebulge but which have not been employed in tuning the viscosity profile to enable ICE-6G_C (VM6) to fit the East coast data set.

Description: In the context of the 2014 realization of the International Terrestrial Reference Frame, the International DORIS (Doppler Orbitography Radiopositioning Integrated by Satellite) Service (IDS) has delivered to the IERS a set of 1140 weekly SINEX files including station coordinates and Earth orientation parameters, covering the time period from 1993.0 to 2015.0. From this set of weekly SINEX files, the IDS combination centre estimated a cumulative DORIS position and velocity solution to obtain mean horizontal and vertical motion of 160 stations at 71 DORIS sites. The main objective of this study is to validate the velocities of the DORIS sites by comparison with external models or time-series. Horizontal velocities are compared with two recent global plate models (GEODVEL 2010 and NNR-MORVEL56). Prior to the comparisons, DORIS horizontal velocities were corrected for Global Isostatic Adjustment from the ICE-6G (VM5a) model. For more than half of the sites, the DORIS horizontal velocities differ from the global plate models by less than 2–3 mm yr –1 . For five of the sites (Arequipa, Dionysos/Gavdos, Manila and Santiago) with horizontal velocity differences with respect to these models larger than 10 mm yr –1 , comparisons with GNSS estimates show the veracity of the DORIS motions. Vertical motions from the DORIS cumulative solution are compared with the vertical velocities derived from the latest GPS cumulative solution over the time span 1995.0–2014.0 from the University of La Rochelle solution at 31 co-located DORIS-GPS sites. These two sets of vertical velocities show a correlation coefficient of 0.83. Vertical differences are larger than 2 mm yr –1 at 23 percent of the sites. At Thule, the disagreement is explained by fine-tuned DORIS discontinuities in line with the mass variations of outlet glaciers. Furthermore, the time evolution of the vertical time-series from the DORIS station in Thule show similar trends to the GRACE equivalent water height.

Description: We have used ensemble averages of satellite-derived free-air gravity anomaly data, together with inverse modelling techniques, to determine the effective elastic thickness, T e , of circum-Pacific subducting oceanic lithosphere and its relationship to plate age. Synthetic modelling tests show that T e can be recovered best using gravity anomaly, rather than bathymetry, data and profiles that are at least 750 km long. Inverse modelling based on a uniform T e elastic plate suggests that T e increases with age of the subducting oceanic lithosphere and is given approximately by the depth to the 390 ± 10 °C oceanic isotherm based on a cooling plate model. Misfits between the observed and calculated gravity anomalies are significantly improved if a mechanically weak zone is included between the trench axis and the outer rise. This weak zone is coincident with observations of bend-faulting and seismicity. Inverse modelling shows that T e landward of the outer rise is generally 40–65 per cent less than the T e seaward of the outer rise. Both landward and seaward T e increases with age of the lithosphere and are given by the depth to the 342–349 °C and 671–714 °C oceanic isotherm, respectively. A dependence of T e on age is consistent with models for the cooling of oceanic lithosphere as it moves away from a mid-ocean ridge and the temperature-dependent ductile creep of oceanic lithospheric minerals such as olivine. By comparing the observed T e to the predicted T e based on laboratory-derived yield strength envelopes and an assumption of elastic-perfectly plastic deformation, we have attempted to constrain the rheology of oceanic lithosphere. Regardless of the assumed friction coefficient, the dry-olivine low-temperature plasticity flow laws of Goetze, Evans & Goetze, Raterron et al . and Mei et al . all provide quite a good fit to the observed T e at circum-Pacific subduction zones. This result contrasts with the Hawaiian Islands, where these flow laws are generally too strong to fit the observations. The discrepancy in rheology within Pacific plate may be caused by differences in the timescale of loading and therefore the amount of viscoelastic stress relaxation that has occurred. Other possibilities include thermal rejuvenation and magma-assisted flexure at the Hawaiian Islands.

Description: Large gaps and inconsistencies remain in published estimates of Nubia–Somalia plate motion based on reconstructions of seafloor spreading data around Africa. Herein, we use newly available reconstructions of the Southwest Indian Ridge at ~1-Myr intervals since 20 Ma to estimate Nubia–Somalia plate motion farther back in time than previously achieved and with an unprecedented degree of temporal resolution. At the northern end of the East African rift, our new estimates of Nubia–Somalia motion for six times from 0.78 Ma to 5.2 Ma differ by only 2 per cent from the rift-normal component of motion that is extrapolated from a recently estimated GPS angular velocity. The rate of rift-normal extension thus appears to have remained steady since at least 5.2 Ma. Our new rotations indicate that the two plates have moved relative to each other since at least 16 Ma and possibly longer. Motion has either been steady since at least 16 Ma or accelerated modestly between 6 and 5.2 Ma. Our Nubia–Somalia rotations predict 42.5 ± 3.8 km of rift-normal extension since 10.6 Ma across the well-studied, northern segment of the Main Ethiopian Rift, consistent with 40–50 km estimates for extension since 10.6 Myr based on seismological surveys of this narrow part of the plate boundary. Nubia–Somalia rotations are also derived by combining newly estimated Somalia–Arabia rotations that reconstruct the post-20-Ma opening of the Gulf of Aden with Nubia–Arabia rotations estimated via a probabilistic analysis of plausible opening scenarios for the Red Sea. These rotations predict Nubia–Somalia motion since 5.2 Myr that is consistent with that determined from Southwest Indian Ridge data and also predict 40 ± 3 km of rift-normal extension since 10.6 Ma across the Main Ethiopian Rift, consistent with our 42.5 ± 3.8 km Southwest Indian Ridge estimate. Our new rotations exclude at high confidence level previous estimates of 12 ± 13 and 123 ± 14 km for rift-normal extensions across the Main Ethiopian Rift since 10.6 Ma based on reconstructions of Chron 5n.2 along the Southwest Indian Ridge. Sparse coverage of magnetic reversals older than 16 Ma along the western third of the Southwest Indian Ridge precludes reliable determinations of Nubia–Somalia plate motion before 16 Ma, leaving unanswered the key question of when the motion between the two plates began.

Description: The TAMMAR segment of the Mid-Atlantic Ridge forms a classic propagating system centred about two degrees south of the Kane Fracture Zone. The segment is propagating to the south at a rate of 14 mm yr –1 , 15 per cent faster than the half-spreading rate. Here, we use seismic refraction data across the propagating rift, sheared zone and failed rift to investigate the crustal structure of the system. Inversion of the seismic data agrees remarkably well with crustal thicknesses determined from gravity modelling. We show that the crust is thickened beneath the highly magmatic propagating rift, reaching a maximum thickness of almost 8 km along the seismic line and an inferred (from gravity) thickness of about 9 km at its centre. In contrast, the crust in the sheared zone is mostly 4.5–6.5 km thick, averaging over 1 km thinner than normal oceanic crust, and reaching a minimum thickness of only 3.5 km in its NW corner. Along the seismic line, it reaches a minimum thickness of under 5 km. The PmP reflection beneath the sheared zone and failed rift is very weak or absent, suggesting serpentinisation beneath the Moho, and thus effective transport of water through the sheared zone crust. We ascribe this increased porosity in the sheared zone to extensive fracturing and faulting during deformation. We show that a bookshelf-faulting kinematic model predicts significantly more crustal thinning than is observed, suggesting that an additional mechanism of deformation is required. We therefore propose that deformation is partitioned between bookshelf faulting and simple shear, with no more than 60 per cent taken up by bookshelf faulting.

Description: Knowledge and understanding of the ocean–continent transition (OCT) structure, continent–ocean boundary (COB) location and crustal type are of critical importance in evaluating rifted continental margin formation and evolution. OCT structure, COB location and magmatic type also have important implications for the understanding of the geodynamics of continental breakup and in the evaluation of petroleum systems in deep-water frontier oil and gas exploration at rifted continental margins. Mapping the distribution of thinned continental crust and lithosphere, its distal extent and the start of unequivocal oceanic crust and hence determining the OCT structure and COB location at rifted continental margins is therefore a generic global problem. In order to assist in the determination of the OCT structure and COB location, we present methodologies using gravity anomaly inversion, residual depth anomaly (RDA) analysis and subsidence analysis, which we apply to the west Iberian rifted continental margin. The west Iberian margin has one of the most complete data sets available for deep magma-poor rifted margins, so there is abundant data to which the results can be calibrated. Gravity anomaly inversion has been used to determine Moho depth, crustal basement thickness and continental lithosphere thinning; subsidence analysis has been used to determine the distribution of continental lithosphere thinning; and RDAs have been used to investigate the OCT bathymetric anomalies with respect to expected oceanic bathymetries at rifted continental margins. These quantitative analytical techniques have been applied to the west Iberian rifted continental margin along profiles IAM9, Lusigal 12 (with the TGS-extension) and ISE-01. Our predictions of OCT structure, COB location and magmatic type (i.e. the volume of magmatic addition, whether the margin is ‘normal’ magmatic, magma-starved or magma-rich) have been tested and validated using ODP wells (Legs 103, 149 and 173), which provide observational constraints on the west Iberian margin.

Description: The large-scale geological evolution of the North Atlantic Realm during the past 450 Myr is largely understood, but crucial elements remain uncertain. These involve the Caledonian orogeny, the formation of the North Atlantic and accompanying igneous activity, and the present-day high topography surrounding the North Atlantic. Teleseismic receiver function interpretation in the Central Fjord Region of East Greenland recently suggested the presence of a fossil Caledonian subduction complex, including a slab of eclogitised mafic crust and an overlying wedge of serpentinised mantle peridotite. Here we further investigate this topic using inverse receiver functions modelling. The obtained velocity models are tested with regard to their consistency with the regional gravity field and topography. We find that the obtained receiver function model is generally consistent with gravity and isostasy. The western part of the section, with topography of 〉1000 m, is clearly supported by the 40-km-thick crust. The eastern part requires additional buoyancy as provided by the hydrated mantle wedge. The geometry, velocities and densities are consistent with interpretation of the lithospheric structure as a fossil subduction zone complex. The spatial relations with Caledonian structures suggest a Caledonian origin. The results indicate that topography is isostatically compensated by density variations within the lithosphere, and that significant dynamic topography is not required at the present-day.

Description: We present a high-resolution shear wave tomography of the Japan subduction zone down to a depth of 700 km, which is determined by inverting a large number of high-quality S -wave arrival-time data from local earthquakes and teleseismic events. The subducting Pacific and Philippine Sea (PHS) slabs are revealed clearly as high-velocity (high- V ) zones, whereas low-velocity (low- V ) anomalies are revealed in the mantle wedge above the two slabs. The PHS slab has subducted aseismically down to a depth of 480 km under the Japan Sea and to a depth of 540 km under the Tsushima Strait. A window is revealed within the aseismic PHS slab, being consistent with P -wave tomography. Prominent low- V and high-Poisson's ratio ( ) anomalies exist below the PHS slab and above the Pacific slab, which reflect hot and wet mantle upwelling caused by the joint effect of deep dehydration of the Pacific slab and convective circulation process in the mantle wedge above the Pacific slab. The hot and wet mantle upwelling has caused the complex geometry and structure of the PHS slab in SW Japan, and contributed to the Quaternary volcanism along the Japan Sea coast. In eastern Japan, low- V zones are revealed at depths of 200–700 km below the Pacific slab, which may reflect hot upwelling from the lower mantle or even the core–mantle boundary.

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

Description: To investigate the coupling relationship between great earthquake and ionosphere, the GPS-derived total electron contents (TECs) by the Center for Orbit Determination in Europe and the foF2 data from the Space Weather Prediction Center were used to analyse the impending ionospheric anomalies before the Iquique Mw8.2 earthquake in Chile on 2014 April 1. Eliminating effects of the solar and geomagnetic activities on ionosphere by the sliding interquartile range with the 27-day window, the TEC analysis results represent that there were negative anomalies occurred on 15th day prior to the earthquake, and positive anomalies appeared in 5th day before the earthquake. The foF2 analysis results of ionosonde stations Jicamarca, Concepcion and Ramey show that the foF2 increased by 40, 50 and 45 per cent, respectively, on 5th day before the earthquake. The TEC anomalous distribution indicates that there was a widely TEC decrement over the epicentre with the duration of 6 hr on 15th day before the earthquake. On 5th day before the earthquake, the TEC over the epicentre increased with the amplitude of 15 TECu, and the duration exceeded 6 hr. The anomalies occurred on the side away from the equator. All TEC anomalies in these days were within the bounds of equatorial anomaly zone where should be the focal area to monitor ionospheric anomaly before strong earthquakes. The relationship between ionospheric anomalies and geomagnetic activity was detected by the cross wavelet analysis, which implied that the foF2 was not affected by the magnetic activities on 15th day and 5th day prior to the earthquake, but the TECs were partially affected by anomalous magnetic activity during some periods of 5th day prior to the earthquake.

Description: We present the first estimates of Southwest Indian Ridge (SWIR) plate motions at high temporal resolution during the Quaternary and Neogene based on nearly 5000 crossings of 21 magnetic reversals out to C6no (19.72 Ma) and the digitized traces of 17 fracture zones and transform faults. Our reconstructions of this slow-spreading mid-ocean ridge reveal several unexpected results with notable implications for regional and global plate reconstructions since 20 Ma. Extrapolations of seafloor opening distances to zero-age seafloor based on reconstructions of reversals C1n (0.78 Ma) through C3n.4 (5.2 Ma) reveal evidence for surprisingly large outward displacement of 5 ± 1 km west of 32°E, where motion between the Nubia and Antarctic plates occurs, but 2 ± 1 km east of 32°E, more typical of most mid-ocean ridges. Newly estimated SWIR seafloor spreading rates are up to 15 per cent slower everywhere along the ridge than previous estimates. Reconstructions of the numerous observations for times back to 11 Ma confirm the existence of the hypothesized Lwandle plate at high confidence level and indicate that the Lwandle plate's western and eastern boundaries respectively intersect the ridge near the Andrew Bain transform fault complex at 32°E and between ~45°E and 52°E, in accord with previous results. The Nubia–Antarctic, Lwandle–Antarctic and Somalia–Antarctic rotation sequences that best fit many magnetic reversal, fracture zone and transform fault crossings define previously unknown changes in the Neogene motions of all three plate pairs, consisting of ~20 per cent slowdowns in their spreading rates at 7.2 $^{+0.9 }_{ -1.4}$ Ma if we enforce a simultaneous change in motion everywhere along the SWIR and gradual 3°–7° anticlockwise rotations of the relative slip directions. We apply trans-dimensional Bayesian analysis to our noisy, best-fitting rotation sequences in order to estimate less-noisy rotation sequences suitable for use in future global plate reconstructions and geodynamic studies. Notably, our new Nubia–Antarctic reconstruction of C5n.2 (11.0 Ma) predicts 20 per cent less opening than do two previous estimates, with important implications for motion that is estimated between the Nubia and Somalia plates. A Nubia–Somalia rotation determined from our Nubia–Antarctic and Somalia–Antarctic plate rotations for C5n.2 (11.0 Ma) predicts cumulative opening of 45 ± 4 km (95 per cent uncertainty) across the northernmost East Africa rift since 11.0 Ma, 70 per cent less than a recent 129 ± 62 km opening estimate based on a now-superseded interpretation of Anomaly 5 along the western portion of the SWIR.

Description: Geodynamic models predict that rifting of thick, ancient continental lithosphere should not occur unless it is weakened by heating and magmatic intrusion. Therefore, the processes occurring along sections of the western branch of the East African Rift, where ~150 km thick, Palaeoproterozoic lithosphere is rifting with no surface expression of magmatism, are a significant challenge to understand. In an attempt to understand the apparently amagmatic extension we probed the regional uppermost mantle for signatures of thermal alteration using compressional ( Vp ) and shear ( Vs ) wave speeds derived from Pn and Sn tomography. Pervasive thermal alteration of the uppermost mantle and possibly the presence of melt can be inferred beneath the Rungwe volcanic centre, but no signatures on a similar scale were discerned beneath amagmatic portions of the western rift branch encompassing the southern half of the Lake Tanganyika rift and much of the Rukwa rift. In this region, Vp and Vs wave speeds indicate little, if any, heating of the uppermost mantle and no studies have reported dyking. Vp / Vs ratios are consistent with typical, melt-free, olivine-dominated upper mantle. Although our resolution limit precludes us from imaging potential localised magmatic intrusions with dimensions of tens of kilometres, the absence of surface volcanism, the amagmatic upper crustal rupture known to have occurred at disparate locations on the western branch, the presence of lower crustal seismicity and the low temperatures implied by the fast seismic wave speeds in the lower crust and uppermost mantle in this region suggests possible amagmatic extension. Most dynamic models predict that this should not happen. Indeed even with magmatic intrusion, rifting of continental lithosphere 〉100 km thick is considered improbable under conditions found on Earth. Yield strength envelopes confirm that currently modelled stresses are insufficient to produce the observed deformation along these portions of the rift system. Stresses arising from the gravitational force related to the uplift of the East African Plateau provide only one-eighth of the minimum stress necessary to produce observed lower crustal earthquakes in the western branch. We expect that some of this disparity may be accounted for by considering smaller scale bending stresses and dynamic feedbacks between brittle and elastic deformation and between faulting, topography and weathering that are not currently included in models of the East African Rift.

Description: A Neogene magmatic reactivation of the Manihiki Plateau, a large igneous province (LIP) in the central Pacific, is studied using seismic reflection data. Igneous diapirs have been identified exclusively within a narrow WNW–ESE striking corridor in the southern High Plateau (HP), which is parallel to the Neogene Pacific Plate motion and overlaps with an extrapolation of the Society Islands Hotspot (SIH) path. The igneous diapirs are characterized by a narrow width (〉5 km), penetration of the Neogene sediments, and they become progressively younger towards the East (23–10 Ma). The magmatic source appears to be of small lateral extent, which leads to the conclusion that the diapirs represent Neogene hotspot volcanism within a LIP, and thus may be an older, previously unknown extension of the SIH track (〉4.5 Ma). Comparing hotspot volcanism within oceanic and continental lithosphere, we further conclude that hotspot volcanism within LIP crust has similarities to tectonically faulted continental crust.

Description: Many mechanisms control the state of stress within Earth plates. First-order well-known mechanisms include stresses induced by lateral variations of lithospheric density structure, sublithospheric tractions, ridge push and subduction pull. In this study, we attempt to quantify the influence of these mechanisms to understand the origin of stresses in the lithosphere, choosing the African plate (TAP) as an example. A finite-element based suite, Proshell, was developed to combine several data sets, to estimate the gravitational potential energy (GPE) of the lithosphere and to calculate stresses acting on the real (non-planar) geometry of TAP. We introduce several quantitative parameters to measure the degree of fit between the model and observations. Our modelling strategy involves nine series of numerical experiments. We start with the simplest possible model and then, step by step, build it up to be a more physically realistic model, all the while discussing the influence of each additional component. The starting (oversimplified) model series (1) is based on the CRUST2 data set for the crust and a half-space-cooling approximation of the lithospheric mantle. We then describe models (series 2–5) that account for lithospheric mantle density heterogeneities to build a more reliable GPE model. The consecutive series involve basal traction from the convective mantle (series A, C) and the rheological heterogeneity of the TAP via variations in its effective elastic thickness (series B, C). The model quality reflects the increase in complexity between series with an improving match to observed stress regimes and directions. The most complex model (series D) also accounts for the bending stresses in the elastic lithosphere and achieves a remarkably good fit to observations. All of our experiments were based on the iteration of controlling parameters in order to achieve the best fit between modelled and observed stresses, always considering physically feasible values. This gives us confidence that our methodology appropriately models the stress pattern of Africa, and that it may be further applied to other plates on Earth. Our modelling approach allows us to quantify several important features controlling the lithospheric stress pattern. Even though the initial oversimplified model does not fit the observations satisfactorily, it shows how ridge push may create significant compressive stresses in the lithosphere. More complex models show the importance of the density structure of the lithosphere, specifically in the subcrustal lithosphere. The stress regime within the TAP mainly results from a global balance of masses and mass moments between continental and oceanic parts of the plate. The orientation of stresses, in turn, is influenced more by local features expressed by topographic and crustal density variations, whereas existent subcrustal density variations appear to be smoothed by the crust above. The models show that accounting separately for either basal tractions or rheological heterogeneities brings moderate improvement, but the combination of these two mechanisms results in a substantially better match between model and observations. The bending stresses caused by isostatical re-adjustment improve the model match, but they have to be analysed with caution because of their depth-dependent nature.

Description: The Red Sea and Gulf of Aden represent two young basins that formed between Africa and Arabia since the early Oligocene, floored by oceanic crust or by transitional and thinned continental crust. While in the easternmost Gulf of Aden, the rift–drift transition can be dated chron C6 (~20.1 Ma), here we show that in the Red Sea the first pulse of seafloor spreading occurred during chron C3n.2n (~4.6 Ma) around ~17.1°N (present-day coordinates) and propagated southwards from this location, separating the Danakil microplate from Arabia. It is also shown that seafloor spreading between Arabia and Nubia started later, around chron 2A (~2.58 Ma), and propagated northwards. At present, there is no magnetic evidence for the existence of a linear spreading centre in the northern Red Sea at latitudes higher than ~24°N and in the southern Red Sea below ~14.8°N. The present-day plate kinematics of this region can be described with high accuracy by a network of five interacting plates (Nubia, Arabia, Somalia, Sinai and Danakil) and six triple junctions. For times older than anomaly 2A (~2.58 Ma) and up to anomaly 3, the absence of marine magnetic anomalies between Arabia and Nubia prevents a rigorous kinematic description of the five-plates system. However, there is strong evidence that the unique changes in plate motions during the last 5 Myr were a dramatic slowdown at chron C2 (~1.77 Ma) in the spreading or extension rates along the ridge and rift axes, thereby a good representation of the real plate motions can be obtained anyway by backward extension of the oldest Arabia–Nubia and Arabia–Danakil stage rotations determined on the basis of marine magnetic anomalies, respectively, C2–C2A and C2A–C3. The proposed kinematic reconstructions are accompanied by a geodynamic explanation for the genesis of large continent–continent fracture zones at the rift–drift transition and by an analysis of the strain associated with plate motions in Afar, northeastern Egypt and Sinai.

Description: On 2016 February 6 the South Taiwan earthquake ( M w 6.4) occurred in the Meinong District of Kaohsiung, southern Taiwan, at a depth of 17 km. It caused 117 fatalities and widespread damage to infrastructures, especially in the Tainan city. To clarify the generating mechanism of this damaging earthquake, we determined high-resolution 3-D images of P - and S -wave velocity ( V P , V S ) and Poisson's ratio ( ) in the epicentral area. We used 105 712 P - and 61 250 S -wave arrival times of 8279 local earthquakes (1.5 ≤ M ≤ 6.4) recorded at 41 seismic stations in South Taiwan during 2000–2011. In the upper crust (depth ≤ 10 km), the most remarkable feature is low- V P , low- V S and high- anomalies in areas with known active faults in the southwestern and easternmost parts of Taiwan. In contrast, high- V P , high- V S and low- anomalies become dominant in the lower crust. The hypocentre of the 2016 South Taiwan earthquake is located in a boundary zone where seismic velocity and Poisson's ratio change drastically in both the horizontal and vertical directions. Furthermore, the hypocentre is underlain by a vertically elongated high- anomaly at depths of 23–40 km, which may reflect ascending fluids from the upper (or uppermost) mantle. The low- V and high- anomalies in the upper crust coincide with areas of low heat flow, negative Bouguer gravity anomaly, and low magnetotelluric resistivity, which may reflect crustal fluids contained in the young fold-and-thrust belt. These results suggest that the 2016 South Taiwan earthquake was triggered by ascending fluids from dehydration of the subducting Eurasian slab, invading into active faults with a high loading rate.

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

Description: In this work, we analyse in detail how a large earthquake could cause stress changes on volcano plumbing systems and produce possible positive feedbacks in promoting new eruptions. We develop a sensitivity analysis that considers several possible parameters, providing also new constraints on the methodological approach. The work is focus on the M w 8.8 2010 earthquake that occurred along the Chile subduction zone near 24 historic/Holocene volcanoes, located in the Southern Volcanic Zone. We use six different finite fault-slip models to calculate the static stress change, induced by the coseismic slip, in a direction normal to several theoretical feeder dykes with various orientations. Results indicate different magnitudes of stress change due to the heterogeneity of magma pathway geometry and orientation. In particular, the N–S and NE–SW-striking magma pathways suffer a decrease in stress normal to the feeder dyke (unclamping, up to 0.85 MPa) in comparison to those striking NW–SE and E–W, and in some cases there is even a clamping effect depending on the magma path strike. The diverse fault-slip models have also an effect (up to 0.4 MPa) on the results. As a consequence, we reconstruct the geometry and orientation of the most reliable magma pathways below the 24 volcanoes by studying structural and morphometric data, and we resolve the stress changes on each of them. Results indicate that: (i) volcanoes where post-earthquake eruptions took place experienced earthquake-induced unclamping or very small clamping effects, (ii) several volcanoes that did not erupt yet are more prone to experience future unrest, from the point of view of the host rock stress state, because of earthquake-induced unclamping. Our findings also suggest that pathway orientation plays a more relevant role in inducing stress changes, whereas the depth of calculation (e.g. 2, 5 or 10 km) used in the analysis, is not key a parameter. Earthquake-induced magma-pathway unclamping might contribute to promote new eruptions at volcanoes as far as 450 km from the epicentre.

Description: Seismic studies indicate that the Earth's inner core has a complex structure and exhibits a strong elastic anisotropy with a cylindrical symmetry. Among the various models which have been proposed to explain this anisotropy, one class of models considers the effect of the Lorentz force associated with the magnetic field diffused within the inner core. In this paper, we extend previous studies and use analytical calculations and numerical simulations to predict the geometry and strength of the flow induced by the poloidal component of the Lorentz force in a neutrally or stably stratified growing inner core, exploring also the effect of different types of boundary conditions at the inner core boundary (ICB). Unlike previous studies, we show that the boundary condition that is most likely to produce a significant deformation and seismic anisotropy is impermeable, with negligible radial flow through the boundary. Exact analytical solutions are found in the case of a negligible effect of buoyancy forces in the inner core (neutral stratification), while numerical simulations are used to investigate the case of stable stratification. In this situation, the flow induced by the Lorentz force is found to be localized in a shear layer below the ICB, whose thickness depends on the strength of the stratification, but not on the magnetic field strength. We obtain scaling laws for the thickness of this layer, as well as for the flow velocity and strain rate in this shear layer as a function of the control parameters, which include the magnitude of the magnetic field, the strength of the density stratification, the viscosity of the inner core and the growth rate of the inner core. We find that the resulting strain rate is probably too small to produce significant texturing unless the inner core viscosity is smaller than about 10 12 Pa s.

Description: We analyse Bouguer anomaly data and previously published Moho depths estimated from receiver functions in order to determine the amount of isostatic compensation or uncompensation of the Rif topography in northern Morocco. We use Moho depth variations extracted from receiver function analyses to predict synthetic Bouguer anomalies that are then compared to observed Bouguer anomaly. We find that Moho depth variations due to isostatic compensation of topographic and/or intracrustal loads do not match Moho depth estimates obtained from receiver function analyses. The isostatic misfit map evidences excess crustal root as large as 10 km in the western part of the study area, whereas a ‘missing’ crustal root of ~5 km appears east of 4.3°E. This excess root/missing topography correlates with the presence of a dense mantle lid, the noticeable southwestward drift of the Western Rif area, and with a current surface uplift. We propose that a delaminated mantle lid progressively detaching westward or southwestward from the overlying crust is responsible for viscous flow of the ductile lower crust beneath the Rif area. This gives rise to isostatic uplift and westward drift due to viscous coupling at the upper/lower crust boundary. At the same time, the presence of this dense sinking mantle lid causes a negative dynamic topography, which explains why the observed topography is too low compared to the crustal thickness.

Description: Flow-like landslides are rapidly moving fluid–solid mixtures that can cause significant destruction along paths that run far from their original sources. Existing models for run out prediction and motion simulation of flow-like landslides have many limitations. In this paper, we develop a new method named ‘Tsunami Squares’ to simulate the generation, propagation and stoppage of flow-like landslides based on conservation of volume and momentum. Landslide materials in the new method form divisible squares that are displaced, then further fractured. The squares move under the influence of gravity-driven acceleration and suffer decelerations due to basal and dynamic frictions. Distinctively, this method takes into account solid and fluid mechanics, particle interactions and flow regime transitions. We apply this approach to simulate the 1982 El Picacho landslide in San Salvador, capital city of El Salvador. Landslide products from Tsunami Squares such as run out distance, velocities, erosion and deposition depths and impacted area agree well with field investigated and eyewitness data.

Description: Numerous regional plate reorganizations and the coeval ages of the Hawaiian Emperor bend (HEB) and Louisville bend of 50–47 Ma have been interpreted as a possible global tectonic plate reorganization at ~chron 21 (47.9 Ma). Yet for a truly global event we would expect a contemporaneous change in Africa absolute plate motion (APM) reflected by physical evidence distributed on the Africa Plate. This evidence has been postulated to take the form of the Réunion-Mascarene bend which exhibits many HEB-like features, such as a large angular change close to ~chron 21. However, the Réunion hotspot trail has recently been interpreted as a sequence of continental fragments with incidental hotspot volcanism. Here we show that the alternative Réunion-Mascarene Plateau trail can also satisfy the age progressions and geometry of other hotspot trails on the Africa Plate. The implied motion, suggesting a pivoting of Africa from 67 to 50 Ma, could explain the apparent bifurcation of the Tristan hotspot chain, the age reversals seen along the Walvis Ridge, the sharp curve of the Canary trail, and the diffuse nature of the St. Helena chain. To test this hypothesis further we made a new Africa APM model that extends back to ~80 Ma using a modified version of the Hybrid Polygonal Finite Rotation Method. This method uses seamount chains and their associated hotspots as geometric constraints for the model, and seamount age dates to determine APM through time. While this model successfully explains many of the volcanic features, it implies an unrealistically fast global lithospheric net rotation, as well as improbable APM trajectories for many other plates, including the Americas, Eurasia and Australia. We contrast this speculative model with a more conventional model in which the Mascarene Plateau is excluded in favour of the Chagos-Laccadive Ridge rotated into the Africa reference frame. This second model implies more realistic net lithospheric rotation and far-field APMs, but fails to explain key details of the Atlantic Ocean volcanic chains. Both models predict a Canary plume influence beneath the Madeiras. Neither model, when projected via the global plate circuit into the Pacific, predicts any significant change in plate motion around chron 21. Consequently, Africa APM models do not appear to provide independent support for a chron 21 global reorganization.

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

Description: The source parameters and slip distribution of the 2013 May 11 M w 6.1 Minab earthquake are studied using seismology, geodesy and field observations. We observe left-lateral strike-slip motion on a fault striking ENE–WSW; approximately perpendicular to previously studied faults in the Minab–Zendan–Palami fault zone. The fault that ruptured in 2013 is one of a series of ~E–W striking left-lateral faults visible in the geology and geomorphology. These accommodate a velocity field equivalent to right-lateral shear on ~N–S striking planes by clockwise rotations about vertical axes. The presence of these faults can reconcile differences in estimates of fault slip rates in the western Makran from GPS and Quaternary dating. The longitudinal range of shear in the western Makran is likely to be controlled by the distance over which the underthrusting Arabian lithosphere deepens in the transition from continent–continent collision in the Zagros to oceanic subduction in the Makran.

Description: The Earth's outer core is a rotating ellipsoidal shell of compressible, stratified and self-gravitating fluid. As such, in the treatment of geophysical problems a realistic model of this body needs to be considered. In this work, we consider compressible and stratified fluid core models with different stratification parameters, related to the local Brunt-Väisälä frequency, in order to study the effects of the core's density stratification on the frequencies of some of the inertial-gravity modes of this body. The inertial-gravity modes of the core are free oscillations with periods longer than 12 hr. Historically, an incompressible and homogeneous fluid is considered to study these modes and analytical solutions are known for the frequencies and the displacement eigenfunctions of a spherical model. We show that for a compressible and stratified spherical core model the effects of non-neutral density stratification may be significant, and the frequencies of these modes may change from model to model. For example, for a spherical core model the frequency of the spin-over mode, the (2, 1, 1) mode, is unaffected while that of the (4, 1, 1) mode is changed from –0.410 for the Poincaré core model to –0.434, –0.447 and –0.483 for core models with the stability parameter β = –0.001, –0.002 and –0.005, respectively, a maximum change of about 18 per cent when β = –0.005. Our results also show that for small stratification parameter, |β| ≤ 0.005, the frequency of an inertial-gravity mode is a nearly linear function of β but the slope of the line is different for different modes, and that the effects of density stratification on the frequency of a mode is likely related to its spatial structure, which remains the same in different Earth models. We also compute the frequencies of some of the modes of the ‘PREM’ (spherical shell) core model and show that the frequencies of these modes may also be significantly affected by non-zero β.

Description: Reaching sub-salt hydrocarbon targets in the deeper part of the Gulf of Mexico requires drilling through a salt canopy. The suture zones in the salt canopy are potential drilling hazards due to anomalous pressure behaviour of entrapped sediments. The Pólya vector field of coalescing salt sheets inside the canopy is used to explain suture formation and distinguish between upright and inclined suture contacts. Our analytical models, based on complex potentials, provide exact solutions for multiple source flows as they compete for space when spreading into the viscous continuum of the salt canopy. The velocity gradient tensor yields the strain rate tensor, which is used to map the principal strain rate magnitude inside the canopy. Quantification of one of the principal strain rates is sufficient because the plane deformation assumption ensures the two principal strain rates are equal in magnitude (but of opposite sign); the third principal dimension can have neither strain nor deviatoric stress. Visualization of the locations where the principal stress vanishes or peaks (with highs and lows) is useful for pre-drilling plans because such peaks must be avoided and the stress-free locations provide the safer drilling sites. A case study—of the Walker Ridge region—demonstrates the practical application of our new method.

Description: The Indochina block is important to our understanding of the extrusion model as a consequence of the Indo-Eurasia collision. The lithospheric structure of this block, however, remains obscured due to a lack of sufficient instrumentation for high resolution seismic imaging. We present a shear velocity model derived from Rayleigh wave phase velocity tomography using data from recently deployed seismic networks in this region. Our inversion results for lithospheric structure show strong correlations with tectonic history in this block. A prominent slow-velocity anomaly (5 per cent) is observed in northern Indochina along the Ailao Shan-Red River (ASRR) shear zone including Chuxiong basin, Lanping-Simao fold belt and Thailand rift basin, which has seen extensive deformation events since Eocene. The Khorat Plateau basin is characterized by thick continental keel type lithosphere, consistent with palaeomagnetic and geological observations indicating this basin has experienced much less deformation than the surrounding regions. Additionally, our inversion imaged a sharp, lithospheric-scale velocity contrast across the southeastern segment of ASRR, indicative of a thin and thus relatively weak lithosphere southwest of Red River Fault. The thin lithosphere, low asthenospheric seismic velocities we observe and the average crustal thicknesses in the region suggest that the topography high is dynamically supported by upwelling asthenosphere rather than thickening of the crust/lithosphere. Based on the occurrence of Palaeogene volcanism and its timing, we prefer an explanation of thinning of the lithosphere and allowing a throughgoing fault rather than emplacement of a thin terrane to explain the thin lithosphere. Therefore, the anomalously thin lithosphere between Khorat Plateau and the ASRR in conjunction with other geological observations is generally consistent with the extrusion model for Indochina, which requires localization of lithospheric deformation around tectonic blocks.

Description: On Earth, oceanic plates subduct beneath a variety of overriding plate (OP) styles, from relatively thin and negatively buoyant oceanic OPs to thick and neutrally/positively buoyant continental OPs. The inclusion of an OP in numerical models of self-consistent subduction has been shown to reduce the rate that subducting slabs roll back relative to the equivalent single plate models. We use dynamic, 2-D subduction models to investigate how the mechanical properties, namely viscosity, thickness, and density, of the OP modify the slab rollback rate and the state of stress of the OP. In addition, we examine the role of the subducting plate (SP) viscosity. Because OP deformation accommodates the difference between the slab rollback rate and the far-field OP velocity, we find that the temporal variations in the rollback rate results in temporal variations in OP stress. The slabs in our models roll back rapidly until they reach the lower mantle viscosity increase, at which point the rollback velocity decreases. Concurrent with this reduction in rollback rate is a switch from an OP dominated by extensional stresses to a compressional OP. As in single plate models, the viscosity of the SP exerts a strong control on subducting slab kinematics; weaker slabs exhibit elevated sinking velocities and rollback rates. The SP viscosity also exerts a strong control on the OP stress regime. Weak slabs, either due to reduced bulk viscosity or stress-dependent weakening, have compressional OPs, while strong slabs have dominantly extensional OPs. While varying the viscosity of the OP alone does not substantially affect the OP stress state, we find that the OP thickness and buoyancy plays a substantial role in dictating the rate of slab rollback and OP stress state. Models with thick and/or negatively buoyant OPs have reduced rollback rates, and increased slab dip angles, relative to slabs with thin and/or positively buoyant OPs. Such elevated trench rollback for models with positively buoyant OPs induces extensional stresses in the OP, while OPs that are strongly negatively buoyant are under compression. While rollback is driven by the negative buoyancy of the subducting slab in such models of free subduction, we conclude that the physical properties of the OP potentially play a significant role in modulating both rollback rates and OP deformation style on Earth.

Description: In the central Andes, the Nazca plate displays large along strike variations in dip with a near horizontal subduction angle between 28 and 32°S referred to the Pampean flat slab segment. The upper plate above the Pampean flat slab has high rates of crustal seismicity and active basement cored uplifts. The SIEMBRA experiment, a 43-broad-band-seismic-station array was deployed to better characterize the Pampean flat slab region around 31°S. In this study, we explore the lithospheric structure above the flat slab as a whole and its relation to seismicity. We use the SIEMBRA data to perform a joint inversion of teleseismic receiver functions and Rayleigh wave phase velocity dispersion to constrain the shear wave velocity variations in the lithosphere. Our joint inversion results show: (1) the presence of several upper-plate mid-crustal discontinuities and their lateral extent that are probably related to the terrane accretion history; (2) zones of high shear wave velocity in the upper-plate lower crust associated with a weak Moho signal consistent with the hypothesis of partial eclogitization in the lower crust; (3) the presence of low shear-wave velocities at ~100 km depth interpreted as the subducting oceanic crust. Finally, in order to investigate the relation of the lithospheric structure to seismicity, we determine an optimal velocity–depth model based on the joint inversion results and use it to perform regional moment tensor inversions (SMTI) of crustal and slab earthquakes. The SMTI for 18 earthquakes that occurred between 2007 and 2009 in the flat slab region below Argentina, indicates systematically shallower focal depths for slab earthquakes (compared with inversions using previous velocity models). This suggests that the slab seismicity is concentrated mostly between 90 and 110 km depths within the subducting Nazca plate's oceanic crust and likely related to dehydration. In addition, the slab earthquakes exhibit extensional focal mechanisms suggesting new faulting at the edges of the flat portion of the slab. SMTI solutions for upper-plate crustal earthquakes match well the geological observations of reactivated structures and agree with crustal shortening. Our new constraints on flat slab structure can aid earthquake characterization for regional seismic hazard assessment and efforts to help understand the mechanisms for slab flattening in the central Andes.

Description: The elastic displacement and stress fields due to rectangular faults and opening-mode fractures within an anisotropic homogeneous half-space are derived in this paper. The solution is expressed in terms of the mathematically elegant and computationally powerful Stroh formalism and can be applied to the generally anisotropic half-space or a transversely isotropic half-space with any oriented isotropic plane. For any flat fault or opening-mode fracture of polygonal shape, one needs only to carry out a simple line integral from 0 to in order to express the fault-induced response. Numerical examples are presented to demonstrate the effect of the anisotropy and fault orientation on the internal and surface responses of the half-space. Our results prove that both rock anisotropy and fault orientation could dramatically change the fields in the domain and one needs to consider these properties as accurately as possible to be able to predict the response in the domain precisely. Anisotropy of the rock mass may alter the dominant displacement and stress components at observation points in the model domain as compared to the isotropic case.

Description: Shear deformation of partially molten rock in laboratory experiments causes the emergence of melt-enriched sheets (bands in cross-section) that are aligned at about 15°–20° to the shear plane. Deformation and deviatoric stress also cause the coherent alignment of pores at the grain scale. This leads to a melt-preferred orientation that may, in turn, give rise to an anisotropic permeability. Here we develop a simple, general model of anisotropic permeability in partially molten rocks. We use linearized analysis and nonlinear numerical solutions to investigate its behaviour under simple-shear deformation. In particular, we consider implications of the model for the emergence and angle of melt-rich bands. Anisotropic permeability affects the angle of bands and, in a certain parameter regime, it can give rise to low angles consistent with experiments. However, the conditions required for this regime have a narrow range and seem unlikely to be entirely met by experiments. Anisotropic permeability may nonetheless affect melt transport and the behaviour of partially molten rocks in Earth's mantle.

Description: Eastern Anatolia region between north–south colliding Arabian and Eurasian plates has no significant crustal root and shallow (upper) mantle flow beneath seems to be vertically supporting its high topography. It has a high surface heat flow and the underlying mantle is characterized by low seismic velocity zones. Using a mantle density/temperature variation field derived from P -wave seismic velocity, current shallow mantle flow and resultant dynamic topography of Eastern Anatolia and adjacent Arabian foreland and Caucasus areas were calculated along a vertical section. The section crosses the tectonic boundaries interrelated with slab bodies (high seismic velocity/cold regions) and the low velocity zones above the slabs. According to the modelling experiments, the surface topography of Eastern Anatolia seems to be supported by shallow mantle flow dynamics. On the other hand, residual topography for the region was calculated using high resolution crustal thickness data. Positive residual topography that suggests an undercompensated state of Eastern Anatolia is in concordance with the dynamic topography anomaly. The modelled local shallow mantle flow support due to the density contrast between hot (low velocity) zones and underlying cold slab bodies beneath the area may be the present-day snapshot of the mantle flow uplift in Eastern Anatolia presence of which was previously suggested.

Description: Crystallization experiments in the dendritic regime have been carried out in hypergravity conditions (from 1 to 1300 g) from an ammonium chloride solution (NH 4 Cl and H 2 O). A commercial centrifuge was equipped with a slip ring so that electric power (needed for a Peltier device and a heating element), temperature and ultrasonic signals could be transmitted between the experimental setup and the laboratory. Ultrasound measurements (2–6 MHz) were used to detect the position of the front of the mushy zone and to determine attenuation in the mush. Temperature measurements were used to control a Peltier element extracting heat from the bottom of the setup and to monitor the evolution of crystallization in the mush and in the liquid. A significant increase of solid fraction and attenuation in the mush is observed as gravity is increased. Kinetic undercooling is significant in our experiments and has been included in a macroscopic mush model. The other ingredients of the model are conservation of energy and chemical species, along with heat/species transfer between the mush and the liquid phase: boundary-layer exchanges at the top of the mush and bulk convection within the mush (formation of chimneys). The outputs of the model compare well with our experiments. We have then run the model in a range of parameters suitable for the Earth's inner core. This has shown the role of bulk mush convection for the inner core and the reason why a solid fraction very close to unity should be expected. We have also run melting experiments: after crystallization of a mush, the liquid has been heated from above until the mush started to melt, while the bottom cold temperature was maintained. These melting experiments were motivated by the possible local melting at the inner core boundary that has been invoked to explain the formation of the anomalously slow F-layer at the bottom of the outer core or inner core hemispherical asymmetry. Oddly, the consequences of melting are an increase in solid fraction and a decrease in attenuation. It is hence possible that surface seismic velocity and attenuation of the inner core are strongly affected by melting.

Description: The coupling that exists between surface processes and deformation within both the shallow crust and the deeper mantle-lithosphere has stimulated the development of computational geodynamic models that incorporate a free surface boundary condition. We introduce a treatment of this boundary condition that is suitable for staggered grid, finite difference schemes employing a structured Eulerian mesh. Our interface capturing treatment discretizes the free surface boundary condition via an interface that conforms with the edges of control volumes (e.g. a ‘staircase’ representation) and requires only local stencil modifications to be performed. Comparisons with analytic solutions verify that the method is first-order accurate. Additional intermodel comparisons are performed between known reference models to further validate our free surface approximation. Lastly, we demonstrate the applicability of a multigrid solver to our free surface methodology and demonstrate that the local stencil modifications do not strongly influence the convergence of the iterative solver.

Description: Observations of the temporal variations in the volume flux of a plume can provide useful constraints on geodynamic models of plumes and plume-plate interactions. Furthermore, they can be compared with observations at other plumes and may be analysed further to understand the nature and cause of the variations. The volume plume flux is typically derived from a sum of edifice and compensation root volumes. The former can be obtained via the application of regional–residual separation procedures that split the observed relief into regional (swell) and residual (edifice) components, while the latter is generally inferred from the former using local (Airy) or regional (flexural) compensation models. Most regional–residual techniques used in past studies give non-unique results and provide no estimates of the uncertainty in the separation, which impacts the significance of the results. Here, the optimal robust separator (ORS) method achieves a unique separation for the swell and edifice components of the Hawaiian Ridge and furthermore obtain confidence bounds on the total volume flux. A fast spectral method for plate flexure with different edifice and infill densities is used to determine compensation volumes. Although my flux estimates have assigned confidence bounds, these are much smaller than the flux estimates themselves. A comparison of my new results to published volume flux curves shows that my revised flux estimates are lower by a factor of 2–3. Reproducing the prior higher results demonstrates that these discrepancies appear to be related to shortcomings in the implementation of the methodology used in the separation. The variability in the Hawaiian plume flux occurs at two different time scales: A short (1–2 Myr) periodicity related to the spacing of islands and seamounts, which ultimately is related to plume-plate flexural interactions, and a much longer (10–15 Myr) periodicity that may be related to plate kinematic changes. Superimposed on these trends may be an exponential increase towards more recent times, but this trend may also be explained by a higher flux level during the period when the plume was positioned beneath the relatively younger lithospheric segment bracketed by the Murray and Molokai fracture zones. Landslides and erosion of the edifice may imply an underestimation of total volumes by 5–10 per cent. The main uncertainty facing studies of plume flux is related to the unknown quantity of magmatic underplating.

Description: Ellesmere Island in Arctic Canada displays a complex geological evolution. The region was affected by two distinct orogenies, the Palaeozoic Ellesmerian orogeny (the Caledonian equivalent in Arctic Canada and Northern Greenland) and the Palaeogene Eurekan orogeny, related to the opening of Baffin Bay and the consequent convergence of the Greenland plate. The details of this complex evolution and the present-day deep structure are poorly constrained in this remote area and deep geophysical data are sparse. Receiver function analysis of seven temporary broad-band seismometers of the Ellesmere Island Lithosphere Experiment complemented by two permanent stations provides important data on the crustal velocity structure of Ellesmere Island. The crustal expression of the northernmost tectonic block of Ellesmere Island (~82°–83°N), Pearya, which was accreted during the Ellesmerian orogeny, is similar to that at the southernmost part, which is part of the Precambrian Laurentian (North America-Greenland) craton. Both segments have thick crystalline crust (~35–36 km) and comparable velocity–depth profiles. In contrast, crustal thickness in central Ellesmere Island decreases from ~24–30 km in the Eurekan fold and thrust belt (~79.7°–80.6°N) to ~16–20 km in the Hazen Stable Block (HSB; ~80.6°–81.4°N) and is covered by a thick succession of metasediments. A deep crustal root (~48 km) at ~79.6°N is interpreted as cratonic crust flexed beneath the Eurekan fold and thrust belt. The Carboniferous to Palaeogene sedimentary succession of the Sverdrup Basin is inferred to be up to 1–4 km thick, comparable to geologically-based estimates, near the western margin of the HSB.

Description: Joint analysis of the seismic velocities and geoid, gravity and gravity gradients are used to constrain the viscosity profile within the mantle as well as the lateral density variations. Recent ESA's Gravity field and steady-state Ocean Circulation Explorer measurements of the second-order derivatives of the Earth's gravity potential give new possibilities to determine these mantle properties. Using a simple mantle model and seismic tomography results, we investigate how the gravitational potential, the three components of the gravity vector and the gravity gradients can bring information on the radial viscosity profile and on the mantle mass anomalies. We start with lateral density variations in the Earth's mantle based either on slab history or deduced from seismic tomography. The main uncertainties are: for the latter case, the relationship between seismic velocity and density—the so-called density/velocity scaling factor—and for the former case, the variation with depth of the density contrast between the cold slabs and the surrounding mantle. We perform a Monte Carlo search for the viscosity and the density/velocity scaling factor profiles within the mantle, which allows to fit the observed geoid, gravity and gradients of gravity. We compute the posterior probability distribution of the unknown parameters, and find that the gravity gradients improve the estimate of the scaling factor within the upper mantle, because of their sensitivity to the masses within the upper mantle, whereas the geoid and the gravity better constrain the scaling factor in the lower mantle. In the upper mantle, it is less than 0.02 in the upper part and about 0.08–0.14 in the lower part, and it is significantly larger for depths greater than 1200 km (about 0.32–0.34). In any case, the density/velocity scaling factor between 670 and 1150 km depth is not well constrained. We show that the viscosity of the upper part of the mantle is strongly correlated with the viscosity of the lower part of the mantle and that the viscosity profile is characterized by a decrease in the lower part of the upper mantle (about 10 20 –2  x  10 20 Pa s) and by an increase (about 10 23 –2  x  10 23 Pa s) at the top of the lower mantle (between 670 and 1150 km). The viscosity of the mantle below 1150 km depth is well estimated in our Monte Carlo search and is about 10 22 –4  x  10 22 Pa s.

Description: Large-scale topography may be due to several causes, including (1) variations in crustal thickness and density structure, (2) oceanic lithosphere age differences, (3) subcrustal density variations in the continental lithosphere and (4) convective flow in the mantle beneath the lithosphere. The last contribution in particular may change with time and be responsible for continental inundations; distinguishing between these contributions is therefore important for linking Earth's history to its observed geological record. As a step towards this goal, this paper aims at such distinction for the present-day topography: the approach taken is deriving a ‘model’ topography due to contributions (3) and (4), along with a model geoid, using a geodynamic mantle flow model. Both lithosphere thickness and density anomalies beneath the lithosphere are inferred from seismic tomography. Density anomalies within the continental lithosphere are uncertain, because they are probably due to variations in composition and temperature, making a simple scaling from seismic to density anomalies inappropriate. Therefore, we test a number of different assumptions regarding these. As a reality check, model topography is compared, in terms of both correlation and amplitude ratio, to ‘residual’ topography, which follows from observed topography after subtracting contributions (1) and (2). The model geoid is compared to observations as well. Comparatively good agreement is found if there is either an excess density of 0.2 per cent in the lithosphere above 150 km depth, with anomalies below as inferred from tomography, or if the excess density is 0.4 per cent in the entire lithosphere. Further, a good fit is found for viscosity 10 20 Pa s in the asthenosphere, increasing to 10 23 Pa s in the lower mantle above D'. Results are quite dependent on which tomography models they are based on; for some recent ones, topography correlation is 0.6, many smaller scale features are matched, topography amplitude is less than 30 per cent too large, while geoid variance reduction exceeds 70 per cent—overall a considerable improvement compared to previous models. Correlation becomes less if smaller scale features (corresponding to spherical harmonic degrees 15 and higher), which are probably largely due to anomalies in the lithosphere, are also considered. Comparison of results with different viscosity structures, and a regional comparison of amplitude ratios, indicates that lateral viscosity variations can be quite strong, but only leading to moderate variations in model topography of a factor probably less than two.

Description: Oblique convergence across Tibet leads to slip partitioning with the coexistence of strike-slip, normal and thrust motion on major fault systems. A key point is to understand and model how faults interact and accumulate strain at depth. Here, we extract ground deformation across the Haiyuan Fault restraining bend, at the northeastern boundary of the Tibetan plateau, from Envisat radar data spanning the 2001–2011 period. We show that the complexity of the surface displacement field can be explained by the partitioning of a uniform deep-seated convergence. Mountains and sand dunes in the study area make the radar data processing challenging and require the latest developments in processing procedures for Synthetic Aperture Radar interferometry. The processing strategy is based on a small baseline approach. Before unwrapping, we correct for atmospheric phase delays from global atmospheric models and digital elevation model errors. A series of filtering steps is applied to improve the signal-to-noise ratio across high ranges of the Tibetan plateau and the phase unwrapping capability across the fault, required for reliable estimate of fault movement. We then jointly invert our InSAR time-series together with published GPS displacements to test a proposed long-term slip-partitioning model between the Haiyuan and Gulang left-lateral Faults and the Qilian Shan thrusts. We explore the geometry of the fault system at depth and associated slip rates using a Bayesian approach and test the consistency of present-day geodetic surface displacements with a long-term tectonic model. We determine a uniform convergence rate of 10 [8.6–11.5] mm yr –1 with an N89 [81–97]°E across the whole fault system, with a variable partitioning west and east of a major extensional fault-jog (the Tianzhu pull-apart basin). Our 2-D model of two profiles perpendicular to the fault system gives a quantitative understanding of how crustal deformation is accommodated by the various branches of this thrust/strike-slip fault system and demonstrates how the geometry of the Haiyuan fault system controls the partitioning of the deep secular motion.

Description: In this paper, we study stability of the rate, state and temperature friction (RSTF) model. The Segall and Rice approach is used to model heat transfer at the sliding interface with its surroundings. The effect of pore pressure is not considered in the model to avoid the complex expression for critical stiffness. Linear stability analysis of the spring-mass sliding system is carried out with the ageing law under the quasistatic conditions in order to determine the critical stiffness above which sliding behaviour changes from unstable to stable or vice versa. Our numerical simulations establish that critical stiffness of the heated surface may increase or decrease from corresponding to the critical stiffness of the unheated surface depending on the relative values of two contradictory parameters related with velocity effect and temperature effect. Parametric studies are also carried out to understand shear velocity and temperature of the sliding surface dependence of steady friction. The RSTF model is also used to study the gravity induced failure of a creeping rock slope and the results are justified.

Description: The Siple Coast region of Antarctica contains a number of fast-flowing ice streams, which control the dynamics and mass balance of the region. These ice streams are known to undergo stagnation and reactivation cycles, which lead to ice thickness changes that may be sufficient to excite a viscous solid Earth response (glacial isostatic adjustment; GIA). This study aims to quantify Siple Coast ice thickness changes during the last 2000 yr in order to determine the degree to which they might contribute to GIA and associated present-day bedrock uplift rates. This is important because accurate modelling of GIA is necessary to determine the rate of present-day ice-mass change from satellite gravimetry. Recently-published reconstructions of ice-stream variability were used to create a suite of kinematic models for the stagnation-related thickening of Kamb Ice Stream since ~1850 AD, and a GIA model was used to predict present-day deformation rates in response to this thickening. A number of longer-term loading scenarios, which include the stagnation and reactivation of ice streams across the Siple Coast over the past 2000 yr, were also constructed, and used to investigate the longer term GIA signal in the region. Uplift rates for each of the ice loading histories, based on a range of earth models, were compared with regional GPS-observed uplift rates and an empirical GIA estimate. We estimate Kamb Ice Stream to have thickened by 70–130 m since stagnation ~165 years ago. Modelled present-day vertical motion in response to this load increase peaks at –17 mm yr –1 (i.e. 17 mm yr –1 subsidence) for the weakest earth models tested here. Comparison of the solid Earth response to ice load changes throughout the last glacial cycle, including ice stream stagnation and reactivation across the Siple Coast during the last 2000 yr, with an empirical GIA estimate suggests that the upper mantle viscosity of the region is greater than 1 x 10 20 Pa s. When upper mantle viscosity values of 1 x 10 20  Pa s or smaller are combined with our suite of ice-load scenarios we predict uplift rates across Siple Coast that are at least 4 mm yr –1 smaller than those predicted by the empirical GIA estimate. GPS data are unable to further constrain model parameters due to the distance of the GPS sites from the study area. Our results demonstrate that Late Holocene ice load changes related to the stagnation and reactivation of ice streams on the Siple Coast may play a dominant role in defining the present-day uplift signal. However, both the detailed Earth structure and deglacial history of the region need to be better constrained in order to reduce uncertainties associated with the GIA signal of this region.

Description: The Castrovillari scarps (Cfs) are located in northern Calabria (Italy) and consist of three main WSW-dipping fault scarps resulting from multiple rupture events. At the surface, these scarps are defined by multiple breaks in slope. Despite its near-surface complexity, the faults likely merge to form a single normal fault at about 200 m depth, which we refer to as the Castrovillari fault. We present the results of a multidisciplinary and multiscale study at a selected site of the Cfs with the aim to (i) characterize the geometry at the surface and at depth and (ii) obtain constraints on the fault slip history. We investigate the site by merging data from quantitative geomorphological analyses, electrical resistivity and ground penetrating radar surveys, and palaeoseismological trenching along a ~40 m high scarp. The closely spaced investigations allow us to reconstruct the shallow stratigraphy, define the fault locations, and measure the faulted stratigraphic offsets down to 20 m depth. Despite the varying resolutions, each of the adopted approaches suggests the presence of sub-parallel fault planes below the scarps at approximately the same location. The merged datasets permit the evaluation of the fault array (along strike for 220 m within a 370-m-wide zone). The main fault zone consists of two closely spaced NW–SE striking fault planes in the upper portion of the scarp slope and another fault at the scarp foot. The 3-D image of the fault surfaces shows west to southwest dipping planes with values between 70° and 80°; the two closely spaced planes join at about 200 m below the surface. The 8-to-12-m-high upper fault, which shows the higher vertical displacements, accommodated most of the deformation during the Holocene. Results from the trenching analysis indicate a minimum slip per event of 0.6 m and a maximum short-term slip rate of 0.6 mm yr –1 for the Cf. The shallow subsurface imaging techniques are particularly helpful in evaluating the possible field uncertainties related to postfaulting modification by erosional/depositional/human processes, such as within stream valleys and urbanized zones.

Description: With the progress of mantle convection modelling over the last decade, it now becomes possible to solve for the dynamics of the interior flow and the surface tectonics to first order. We show here that tectonic data (like surface kinematics and seafloor age distribution) and mantle convection models with plate-like behaviour can in principle be combined to reconstruct mantle convection. We present a sequential data assimilation method, based on suboptimal schemes derived from the Kalman filter, where surface velocities and seafloor age maps are not used as boundary conditions for the flow, but as data to assimilate. Two stages (a forecast followed by an analysis) are repeated sequentially to take into account data observed at different times. Whenever observations are available, an analysis infers the most probable state of the mantle at this time, considering a prior guess (supplied by the forecast) and the new observations at hand, using the classical best linear unbiased estimate. Between two observation times, the evolution of the mantle is governed by the forward model of mantle convection. This method is applied to synthetic 2-D spherical annulus mantle cases to evaluate its efficiency. We compare the reference evolutions to the estimations obtained by data assimilation. Two parameters control the behaviour of the scheme: the time between two analyses, and the amplitude of noise in the synthetic observations. Our technique proves to be efficient in retrieving temperature field evolutions provided the time between two analyses is 10 Myr. If the amplitude of the a priori error on the observations is large (30 per cent), our method provides a better estimate of surface tectonics than the observations, taking advantage of the information within the physics of convection.

Description: In this study, the Moho depth is estimated based on the refined spherical Bouguer gravity disturbance and DTM2006 topographic data using the Vening Meinesz-Moritz gravimetric-isostatic hypothesis. In this context, we compute the refined spherical Bouguer gravity disturbances in a set of 1° x 1° blocks. The spherical terrain correction, a residual correction to each Bouguer shell, is computed using rock heights and ice sheet thicknesses from the DTM2006 and Earth2014 models. The study illustrates that the defined simple Bouguer gravity disturbance corrected for the density variations of the oceans, ice sheets and sediment basins and also the non-isostatic effects needs a significant terrain correction to become the refined Bouguer gravity disturbance, and that the isostatic gravity disturbance is significantly better defined by the latter disturbance plus a compensation attraction. Our study shows that despite the fact that the lateral variation of the crustal depth is rather smooth, the terrain affects the result most significantly in many areas. The global numerical results show that the estimated Moho depths by the simple and refined spherical Bouguer gravity disturbances and the seismic CRUST1.0 model agree to 5.6 and 2.7 km in RMS, respectively. Also, the mean value differences are 1.7 and 0.2 km, respectively. Two regional numerical studies show that the RMS differences between the Moho depths estimated based on the simple and refined spherical Bouguer gravity disturbance and that using CRUST1.0 model yield fits of 4.9 and 3.2 km in South America and yield 3.2 and 3.4 km in Fennoscandia, respectively.

Description: The rifting episode that occurred in Dabbahu–Manda–Hararo (Ethiopia) between 2005 and 2010 during which 14 dyke intrusions were emitted, was a unique opportunity to study interactions between tectonic deformation and magmatic processes. While magmatism has been shown to control primarily the spatial and temporal distribution of dyke intrusions during this accretion sequence, the role of faults in accommodating plate spreading in rift segments is poorly understood. During interdyking periods, transient ground deformation due to magma movement is generally observed. Investigating such a small-scale deformation and in particular the movement along faults during these periods will help understanding the factors that trigger fault movement in magmatic rifts. We analyse fault activity during three interdyking periods: 2006 December–June (d0–d1), 2007 January–July (d5–d6) and 2009 November–January (d10–d11). The time–space evolution of surface displacements along ~700 faults is derived from pairs of ascending and descending SAR interferograms. Surface slip distributions are then compared with codyking ground deformation fields. The results show that faults are mainly activated above the areas affected by magma emplacement during interdyking periods. A detailed analysis of brittle deformation during the six months following the 2005 September intrusion shows asymmetric deformation on the rift shoulders, with significant opening on faults located to the west of the dyke. We explain this feature by the activation of westward dipping pre-existing faults, with block rotations in between. In addition, we observe that the strip encompassing the activated faults narrows by 30 per cent from co- to interdyking period. This suggests that magma keeps migrating to shallower depths after the dyke intrusion. During a rifting episode, activation of faults in a pre-existing fracture network therefore seems to be mainly controlled by deep magma processes.

Description: We investigate how uncertainties in seismic and density structure of the crust propagate to uncertainties in mantle density structure. The analysis is based on interpretation of residual upper-mantle gravity anomalies which are calculated by subtracting (stripping) the gravitational effect of the crust from the observed satellite gravity field data (GOCE Direct release 3). Thus calculated residual mantle gravity anomalies are caused mainly by a heterogeneous density distribution in the upper mantle. Given a relatively small range of expected compositional density variations in the lithospheric mantle, knowledge on uncertainties associated with incomplete information on crustal structure is of utmost importance for progress in gravity modelling. Uncertainties in the residual upper-mantle gravity anomalies result chiefly from uncertainties in (i) seismic V P velocity–density conversion for the crust and (ii) uncertainties in the seismic crustal structure (thickness and average V P velocities of individual crustal layers, including the sedimentary cover). We examine the propagation of these uncertainties into determinations of lithospheric mantle density and analyse both sources of possible uncertainties by applying different velocity-to-density conversions and by introducing variations into the crustal structure which correspond to typical resolution of high-quality and low-quality seismic models. We apply our analysis to Siberia (the West Siberian Basin and the Siberian Craton) for which a new regional seismic crustal model, SibCrust, has recently become available. For the same region, we also compute upper-mantle gravity and density anomalies based on three global crustal models (CRUST 5.1, CRUST 2.0 and CRUST 1.0) and compare the results based on four different crustal models. A large uncertainty in the V P -to-density conversion may result in the uncertainty in lithospheric mantle density anomalies of ca. 0.02–0.03 g cm –3 (i.e. 0.5–1 per cent, which is comparable to compositional density anomalies expected for continental lithosphere mantle). Similar values of uncertainties may be caused by a 0.2 km s –1 error in average crustal V P velocities or by a 2 km uncertainty in the Moho depth. One of the largest uncertainties is caused by errors in thickness of the sedimentary layer, and a 2 km error leads to ca. 0.03 g cm –3 error in lithospheric mantle densities. Large deviations (locally ±10 km) of the Moho depth in global crustal models (CRUST 5.1, CRUST2.0 and CRUST1.0) from the high-resolution regional seismic model of the crust, SibCrust, may produce artefact residual mantle gravity anomalies of up to ±150 mGal locally, caused by large errors in crustal gravity corrections. These errors in gravity anomalies produce up to ca. 0.04 g cm –3 ( ca. 1.2 per cent) errors in density of the lithospheric mantle, which may well correspond to the amplitude of real density anomalies in the mantle. Our results demonstrate that gravity modelling alone cannot reliably constrain the crustal structure, including the Moho depth and thickness of sediments.

Description: For a period of about 1 yr between the summers of 2010 and 2011, 25 broad-band seismographs were deployed in a roughly linear array across the eastern end of the Qaidam basin and the Qilian Shan in the northeastern Tibetan plateau. This region is probably the most suitable place to study the ongoing convergence interaction between the high Tibetan plateau and the main Asian continental plate. Low-frequency P receiver function analysis of the data provides an image of the crust and mantle down to 700 km depth. In addition to the Moho at 45–65 km depth beneath the profile, the 410 and 660 km discontinuities bounding the mantle transition zone can be identified at 400–410 and 650–660 km depths, respectively. A possible increase in temperature in the upper mantle thought to exist beneath the northern part of the high Tibetan plateau is thus confined to this part of the plateau and lower upper-mantle temperatures similar to those beneath southern Tibet occur beneath the Qaidam basin and Qilian Shan. When higher frequencies are included in the P receiver function analysis, a positive Ps converter dipping down to the south from 70–75 km depth at 37.9°N to about 110 km depth at 36°N is imaged. As this feature is only seen in high-frequency images and not in the low-frequency image, it is modelled as the positive Ps conversion from the base of an approximately 5-km-thick anisotropic layer at the top of the Asian mantle lithosphere which is currently subducting. This south-dipping converter continues to the south on the INDEPTH IV profile. S receiver function analysis completes the image of the structure below the Qilian Shan profile with the identification of the lithosphere–asthenosphere boundary (LAB). The LAB of the Asian Plate is identified for a reference slowness of 6.4 s deg –1 at 12–14 s (105–125 km depth) between 38 and 41°N below the northern part of the S receiver function profile. To the south it increases in depth such that it is at about 19 s (170 km depth) between 34 and 35°N at the southern end of the profile. The LAB of the Asian Plate occurs at similar depths on the INDEPTH IV profile at the latitudes where the INDEPTH IV and Qilian Shan profiles overlap. As on the INDEPTH IV profile to the south, between 34 and 35°N at the southern end of the Qilian Shan profile there is evidence from the S receiver functions for the LAB of a separate Tibetan Plate.

Description: Lower and upper bounds for present deformation rates across faults in central California between the San Andreas Fault and Pacific coast are estimated from a new Global Positioning System (GPS) velocity field for central, western California in light of geodetic evidence presented in a companion paper for slow, but significant deformation within the Pacific Plate between young seafloor in the eastern Pacific and older seafloor elsewhere on the plate. Transects of the GPS velocity field across the San Andreas Fault between Parkfield and San Juan Buatista, where fault slip is dominated by creep and the velocity field thus reveals the off-fault deformation, show that GPS sites in westernmost California move approximately parallel to the fault at an average rate of 3.4 ± 0.4 mm yr –1 relative to the older interior of the Pacific Plate, but only 1.8 ± 0.6 mm yr –1 if the Pacific Plate frame of reference is corrected for deformation within the plate. Modelled interseismic elastic deformation from the weakly coupled creeping segment of the San Andreas Fault is an order-of-magnitude too small to explain the southeastward motions of coastal sites in western California. Similarly, models that maximize residual viscoelastic deformation from the 1857 Fort Tejon and 1906 San Francisco earthquakes mismatch both the rates and directions of GPS site motions in central California relative to the Pacific Plate. Neither thus explains the site motions southwest of the San Andreas fault, indicating that the site motions measure deformation across faults and folds outboard of the San Andreas Fault. The non-zero site velocities thus constitute strong evidence for active folding and faulting outboard from the creeping segment of the San Andreas Fault and suggest limits of 0–2 mm yr –1 for the Rinconada Fault slip rate and 1.8 ± 0.6 to 3.4 ± 0.4 mm yr –1 for the slip rates integrated across near-coastal faults such as the Hosgri, San Gregorio and San Simeon faults.

Description: We combine new, well-determined GPS velocities from Clarion, Guadalupe and Socorro islands on young seafloor in the eastern Pacific basin with newly estimated velocities for 26 GPS sites from older seafloor in the central, western and southern parts of the Pacific Plate to test for deformation within the interior of the Pacific Plate and estimate the viscosity of the asthenosphere below the plate. Relative to a Pacific Plate reference frame defined from the velocities of the 26 GPS sites in other areas of the Pacific Plate, GPS sites on Clarion and Guadalupe islands in the eastern Pacific move 1.2 ± 0.6 mm yr –1 (1) towards S09°W ± 38° and 1.9 ± 0.3 mm yr –1 towards S19°E ± 10°, respectively. The two velocities, which are consistent within their 95 per cent uncertainties, both differ significantly from Pacific Plate motion. Transient volcanic deformation related to a 1993–1996 eruption of the Socorro Island shield volcano renders our GPS velocity from that island unreliable for the tectonic analysis although its motion is also southward like those of Clarion and Guadalupe islands. We test but reject the possibilities that drift of Earth's origin in ITRF2008 or unmodelled elastic offsets due to large-magnitude earthquakes around the Pacific rim since 1993 can be invoked to explain the apparent slow southward motions of Clarion and Guadalupe islands. Similarly, corrections to the Pacific Plate GPS velocity field for possible viscoelastic deformation triggered by large-magnitude earthquakes since 1950 also fail to explain the southward motions of the two islands. Viscoelastic models with prescribed asthenospheric viscosities lower than 1  x 10 19 Pa s instead introduce statistically significant inconsistencies into the Pacific Plate velocity field, suggesting that the viscosity of the asthenosphere below the plate is higher than 1  x 10 19 Pa s. Elastic deformation from locked Pacific–North America Plate boundary faults is also too small to explain the southward motions of the two islands. Horizontal thermal contraction of the plate interior may explain the motion observed at Clarion and Guadalupe islands, as might long-term tectonic deformation of the plate interior.

Description: In this study, a new method for computing the sensitivity of the glacial isostatic adjustment (GIA) forward solution with respect to the Earth's mantle viscosity, the so-called the forward sensitivity method (FSM), and a method for computing the gradient of data misfit with respect to viscosity parameters, the so-called adjoint-state method (ASM), are presented. These advanced formal methods complement each other in the inverse modelling of GIA-related observations. When solving this inverse problem, the first step is to calculate the forward sensitivities by the FSM and use them to fix the model parameters that do not affect the forward model solution, as well as identifying and removing redundant parts of the inferred viscosity structure. Once the viscosity model is optimized in view of the forward sensitivities, the minimization of the data misfit with respect to the viscosity parameters can be carried out by a gradient technique which makes use of the ASM. The aim is this paper is to derive the FSM and ASM in the forms that are closely associated with the forward solver of GIA developed by Martinec. Since this method is based on a continuous form of the forward model equations, which are then discretized by spectral and finite elements, we first derive the continuous forms of the FSM and ASM and then discretize them by the spectral and finite elements used in the discretization of the forward model equations. The advantage of this approach is that all three methods (forward, FSM and ASM) have the same matrix of equations and use the same methodology for the implementation of the time evolution of stresses. The only difference between the forward method and the FSM and ASM is that the different numerical differencing schemes for the time evolution of the Maxwell and generalized Maxwell viscous stresses are applied in the respective methods. However, it requires only a little extra computational time for carrying out the FSM and ASM numerically. An straightforward approach to compute the gradient of the data misfit is the brute-force method, whereby the partial derivatives of the misfit with respect to model parameters are approximated by the centred difference of two forward model runs. Although the brute-force method is useful for computing the gradient of the data misfit with respect to a small number of model parameters, it becomes expensive for a viscosity model with a large number of parameters. The ASM offers an efficient alternative for computing the gradient of the misfit since the computational time of the ASM is independent of the number of viscosity parameters. The ASM is thus highly efficient for calculating the gradient of the misfit for models with large numbers of parameters. However, the forward-model solution for each time step must be stored, hence the memory demands scale linearly with the number of time steps. This is the main drawback of the ASM.

Description: Relative to the gravitational potential energy of the Earth's monopole, the multipole energy has received far less attention. In this paper, we recapitulate the basic physics from first principles and derive the formulas for multipole energies in analogy to classical electrostatic theory. We focus on the zonal quadrupole energy associated with the Earth's oblateness, the dominant term in Earth's gravity field apart from the monopole. We find the gravitational energy E oblateness 10 –6 | E monopole | = +2.5 x 10 26 J. We examine the implications of E oblateness and its changes associated with long-term ‘secular’ decreases in the oblateness parameter J 2 . We find the rate of loss of E oblateness due to the Earth rounding induced by the present-day GIA is about –200 GW, an amount quite significant in the kinetic energy budget of the mantle heat engine that drives the plate tectonics that has been estimated to be ~1 TW. We also assert that the tidal braking and the global earthquake dislocations, both resulting in Earth rounding on long-term geological timescales, are accompanied with a secular decrease of E oblateness at nearly the same rate of several GW.

Description: Large-scale chemical lateral heterogeneities are inferred in the Earth's lowermost mantle by seismological studies. We explore the model space of thermochemical convection that can maintain reservoirs of dense material for a long period of time, by using similar analysis in 3-D spherical geometry. In this study, we focus on the parameters thought to be important in controlling the stability and structure of primordial dense reservoirs in the lower mantle, including the chemical density contrast between the primordial dense material and the regular mantle material (buoyancy ratio), thermal and chemical viscosity contrasts, volume fraction of primordial dense material and the Clapeyron slope of the phase transition at 660 km depth. We find that most of the findings from the 3-D Cartesian study still apply to 3-D spherical cases after slight modifications. Varying buoyancy ratio leads to different flow patterns, from rapid upwelling to stable layering; and large thermal viscosity contrasts are required to generate long wavelength chemical structures in the lower mantle. Chemical viscosity contrasts in a reasonable range have a second-order role in modifying the stability of the dense anomalies. The volume fraction of the initial primordial dense material does not effect the results with large thermal viscosity contrasts, but has significant effects on calculations with intermediate and small thermal viscosity contrasts. The volume fraction of dense material at which the flow pattern changes from unstable to stable depends on buoyancy ratio and thermal viscosity contrast. An endothermic phase transition at 660 km depth acts as a ‘filter’ allowing cold slabs to penetrate while blocking most of the dense material from penetrating to the upper mantle.

Description: Relative sea level curves contain coupled information about absolute sea level change and vertical lithospheric movement. Such curves may be constructed based on, for example tide gauge data for the most recent times and different types of geological data for ancient times. Correct account for vertical lithospheric movement is essential for estimation of reliable values of absolute sea level change from relative sea level data and vise versa. For modern times, estimates of vertical lithospheric movement may be constrained by data (e.g. GPS-based measurements), which are independent from the relative sea level data. Similar independent data do not exist for ancient times. The purpose of this study is to test two simple inversion approaches for simultaneous estimation of lithospheric uplift rates and absolute sea level change rates for ancient times in areas where a dense coverage of relative sea level data exists and well-constrained average lithospheric movement values are known from, for example glacial isostatic adjustment (GIA) models. The inversion approaches are tested and used for simultaneous estimation of lithospheric uplift rates and absolute sea level change rates in southwest Scandinavia from modern relative sea level data series that cover the period from 1900 to 2000. In both approaches, a priori information is required to solve the inverse problem. A priori information about the average vertical lithospheric movement in the area of interest is critical for the quality of the obtained results. The two tested inversion schemes result in estimated absolute sea level rise of ~1.2/1.3 mm yr –1 and vertical uplift rates ranging from approximately –1.4/–1.2 mm yr –1 (subsidence) to about 5.0/5.2 mm yr –1 if an a priori value of 1 mm yr –1 is used for the vertical lithospheric movement throughout the study area. In case the studied time interval is broken into two time intervals (before and after 1970), absolute sea level rise values of ~0.8/1.2 mm yr –1 (before 1970) and ~2.0 mm yr –1 (after 1970) are found. The uplift patterns resulting from the different inversions suggest that the lithospheric post-GIA response changes near the border between the Danish Basin and the Fennoscandian Shield. The obtained patterns of vertical lithospheric movement rates are comparable to results from other studies based on different and similar data types. Main differences between the inversion results and the results from other studies are caused by factors such as the simplifications included in the inversion approach, such as neglecting local sea level variation caused by the dominant wind patterns, and the a priori values chosen for the vertical uplift rates. The tests of the inversion schemes reveal that realistic values of absolute sea level rise and lithospheric uplift may be simultaneously estimated provided that reliable prior knowledge regarding the overall lithospheric uplift in the study area is available beforehand. In the presented parametrizations, only one absolute sea level change rate value is estimated for each studied time interval while several vertical movement rates are found, and the inverse estimate of absolute sea level change rate is practically insensitive with respect to the choice of a priori value of absolute sea level change, as long as the uncertainty assigned to this a priori value is kept sufficiently high.

Description: In 1356, a magnitude 6–7 earthquake occurred near Basel, in Switzerland. But recent compilations of GPS measurements reveal that measured horizontal deformation rates in northwestern continental Europe are smaller than error bars on the measurements, proving present tectonic activity, if any, is very small in this area. We propose to reconcile these apparently antinomic observations with a mechanical model of the lithosphere that takes into account the geometry of the lithosphere–asthenosphere boundary, assuming that the only loading mechanism is gravity. The lithosphere is considered to be an elastoplastic material satisfying a Von Mises plasticity criterion. The model, which is 400 km long, 360 km wide and 230 km thick, is centred near Belfort in eastern France, with its width oriented parallel to the N145°E direction. It also takes into account the real topography of both the ground surface and that of the Moho discontinuity. Not only does the model reproduce observed principal stress directions orientations, it also identifies a plastic zone that fits roughly the most seismically active domain of the region. Interestingly, a somewhat similar stress map may be produced by considering an elastic lithosphere and an ad-hoc horizontal ‘tectonic’ stress field. However, for the latter model, examination of the plasticity criterion suggests that plastic deformation should have taken place. It is concluded that the present-day stress field in this region is likely controlled by gravity and rheology, rather than by active Alpine tectonics.

Description: We have investigated variations in transition zone thickness under the Borborema Province of NE Brazil by migrating and stacking teleseismic P -wave receiver functions at 32 seismic stations in the region. The Borborema Province represents the western portion of a larger Neoproterozoic mobile belt that occupied much of northern Gondwana, where extensional processes in the Mesozoic lead to the formation of a number of intracontinental basins and ultimately continental breakup. Episodes of intraplate volcanism and uplift marked the evolution of the Province during the Cenozoic, but it is unclear whether those episodes originated from shallow or deep-seated magmatic sources. On one hand, the elliptical shape of the uplifted area, the stress pattern of the Cenozoic deformation and the time overlap between uplift and volcanism suggest doming from thermal activation due to a deep-seated mantle plume. On the other hand, geochronological dates of volcanic bodies in the Province are better understood if resulting from lithospheric erosion by a shallow, small-scale convection cell. Large temperature anomalies are expected to be associated with mantle upwellings, and constraints on the depth extent of the upwellings can be obtained from transition zone thickness. Thinning of the transition zone with respect to its nominal 250 km value is considered diagnostic for positive temperature anomalies, while thickening is considered diagnostic for negative anomalies. Our results show that transition zone thickness is normal, around 250 km, throughout the Province and suggest that thermal perturbations—if present—are confined to the upper mantle. We argue that our results are consistent with a local, shallow magmatic source for the Cenozoic intraplate volcanism of the Borborema Province, although other proposed scenarios—such as channeling of upwelling plume material along lithospheric thin spots—cannot be ruled out with our analysis.

Description: Standard techniques for computed tomography imaging are not directly applicable to a carbonate rock because of the geometric complexity of its pore space. In this study, we first characterized the pore structure in Majella limestone with 30 per cent porosity. Microtomography data acquired on this rock was partitioned into three distinct domains: macropores, solid grains, and an intermediate domain made up of voxels of solid embedded with micropores below the resolution. A morphological analysis of the microtomography images shows that in Majella limestone both the solid and intermediate domains are interconnected in a manner similar to that reported previously in a less porous limestone. We however show that the macroporosity in Majella limestone is fundamentally different, in that it has a percolative backbone which may contribute significantly to its permeability. We then applied for the first time 3-D-volumetric digital image correlation (DIC) to characterize the mode of mechanical failure in this limestone. Samples were triaxially deformed over a wide range of confining pressures. Tomography imaging was performed on these samples before and after deformation. Inelastic compaction was observed at all tested pressures associated with both brittle and ductile behaviors. Our DIC analysis reveals the structure of compacting shear bands in Majella limestone deformed in the transitional regime. It also indicates an increase of geometric complexity with increasing confinement—from a planar shear band, to a curvilinear band, and ultimately to a diffuse multiplicity of bands, before shear localization is inhibited as the failure mode completes the transition to delocalized cataclastic flow.

Description: The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995–1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr –1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N–S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the ‘r’ factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2 uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.

Description: The resistivity structure of the lithospheric mantle beneath the Proterozoic Grenville Province in southern Ontario, Canada is investigated using 84 magnetotelluric (MT) sites divided into four profiles. Depth-based regional geoelectric dimensionality analyses of the MT responses indicate that the mantle lithosphere north of Lake Ontario can be subdivided into upper (45–150 km) and deeper (〉200 km) lithospheric mantle layers with regional strike azimuths of N85°E (±5°) and N65°E (±5°), respectively. MT responses from the Grenville Front and the northwest part of the Central Gneiss Belt are compatible with the presence of 2-D resistivity structures but farther to the southeast, in the southeast part of the Central Gneiss Belt and Central Metasedimentary Belt, they suggest the presence of localized 3-D structures. 2-D inversion of distortion-free MT responses images a large scale very resistive (〉20 000  m) region that extends 300 km southeast of the Grenville Front and for at least 800 km along-strike in the lithospheric mantle beneath the Grenville Province. This feature is interpreted to be Superior Province lithosphere and the corresponding N85°E geoelectric strike to be associated with the fabric of the Superior Province. The base of the resistor reaches depths of 280 km on two of the three MT profiles north of Lake Ontario and this depth is interpreted to be the base of the lithosphere. A large region of enhanced conductivity in the lower lithosphere, spatially correlated with decreased seismic velocity, is bounded to the northwest by a subvertical resistivity anomaly located near the Kirkland Lake and Cobalt kimberlite fields. The enhanced conductivity in the lower lithosphere is attributed to refertilization by fluids associated with Cretaceous kimberlite magmatism and can be explained by water content in olivine of 50 wt ppm in background areas with higher values in a localized anomaly beneath the kimberlite fields. Farther to the southeast the resistivity models include a lithospheric conductor between 100 and 150 km depth beneath the Central Metasedimentary Belt. The enhanced conductivity is attributed to grain boundary graphite films, associated with the Cretaceous kimberlitic magmatic process, or to water and carbon, introduced into the mantle during the pre-Grenvillian tectonism.

Description: Displacements and stress-field changes associated with earthquakes, volcanoes, landslides and human activity are often simulated using numerical models in an attempt to understand the underlying processes and their governing physics. The application of elastic dislocation theory to these problems, however, may be biased because of numerical instabilities in the calculations. Here, we present a new method that is free of artefact singularities and numerical instabilities in analytical solutions for triangular dislocations (TDs) in both full-space and half-space. We apply the method to both the displacement and the stress fields. The entire 3-D Euclidean space $\mathbb {R}^{3}$ is divided into two complementary subspaces, in the sense that in each one, a particular analytical formulation fulfils the requirements for the ideal, artefact-free solution for a TD. The primary advantage of the presented method is that the development of our solutions involves neither numerical approximations nor series expansion methods. As a result, the final outputs are independent of the scale of the input parameters, including the size and position of the dislocation as well as its corresponding slip vector components. Our solutions are therefore well suited for application at various scales in geoscience, physics and engineering. We validate the solutions through comparison to other well-known analytical methods and provide the MATLAB codes.

Description: Models of the glacial isostatic adjustment process, which is dominated by the influence of the Late Pleistocene cycle of glaciation and deglaciation, depend on two fundamental inputs: a history of ice-sheet loading and a model of the radial variation of mantle viscosity. These models may be tested and refined by comparing their local predictions of relative sea level history to geological inferences based upon appropriate sea level indicators. The U.S. Atlantic coast is a region of particular interest in this regard, due to the fact that data from the length of this coast provides a transect of the forebulge associated with the former Laurentide ice sheet. High-quality relative sea level histories from this region are employed herein to explore the ability of current models of mantle viscosity to explain the inferred evolution of relative sea level that have accompanied forebulge collapse following deglaciation. Existing misfits are characterized, and alternatives are explored for their reconciliation. It is demonstrated that a new model of mantle viscosity, referred to herein as VM6, when coupled with the latest model of deglaciation history ICE-6G_C, is able to eliminate the majority of these misfits, while continuing to reconcile a wide range of other important geophysical observables, as well as additional relative sea level data from the North American. West coast which also record the collapse of the forebulge but which have not been employed in tuning the viscosity profile to enable ICE-6G_C (VM6) to fit the East coast data set.

Description: The origin of shear wave splitting (SWS) in the mantle beneath Iceland is examined using numerical models that simulate 3-D mantle flow and the development of seismic anisotropy due to lattice-preferred orientation (LPO). Using the simulated anisotropy structure, we compute synthetic SKS waveforms, invert them for fast polarization directions and split times, and then compare the predictions with the results from three observational studies of Iceland. Models that simulate a mantle plume interacting with the Mid-Atlantic Ridge in which the shallow-most mantle has a high viscosity due to the extraction of water with partial melting, or in which C-type olivine LPO fabric is present due to high water content in the plume, produce the largest chi-squared misfits to the SWS observations and are thus rejected. Models of a low-viscosity mantle plume with A-type olivine fabric everywhere, or with the added effects of E-type fabric in the plume below the solidus produce lower misfits. The lowest misfits are produced by models that include a rapid (~50 km Myr –1 ) northward regional flow (NRF) in the mid-upper mantle, either with or without a plume. NRF was previously indicated by a receiver function study and a regional tomography study, and is shown here to be a major cause of the azimuthal anisotropy beneath Iceland. The smallest misfits for the models with both a plume and NRF are produced when LPO forms above depths of 300–400 km, which, by implication, also mark the depths above which dislocation creep dominates over diffusion creep. This depth of transition between dislocation and diffusion creep is greater than expected beneath normal oceanic seafloor, and is attributed to the unusually rapid strain rates associated with an Iceland plume and the NRF.

Description: The resistivity structure of the lithospheric mantle beneath the Proterozoic Grenville Province in southern Ontario, Canada is investigated using 84 magnetotelluric (MT) sites divided into four profiles. Depth-based regional geoelectric dimensionality analyses of the MT responses indicate that the mantle lithosphere north of Lake Ontario can be subdivided into upper (45–150 km) and deeper (〉200 km) lithospheric mantle layers with regional strike azimuths of N85°E (±5°) and N65°E (±5°), respectively. MT responses from the Grenville Front and the northwest part of the Central Gneiss Belt are compatible with the presence of 2-D resistivity structures but farther to the southeast, in the southeast part of the Central Gneiss Belt and Central Metasedimentary Belt, they suggest the presence of localized 3-D structures. 2-D inversion of distortion-free MT responses images a large scale very resistive (〉20 000  m) region that extends 300 km southeast of the Grenville Front and for at least 800 km along-strike in the lithospheric mantle beneath the Grenville Province. This feature is interpreted to be Superior Province lithosphere and the corresponding N85°E geoelectric strike to be associated with the fabric of the Superior Province. The base of the resistor reaches depths of 280 km on two of the three MT profiles north of Lake Ontario and this depth is interpreted to be the base of the lithosphere. A large region of enhanced conductivity in the lower lithosphere, spatially correlated with decreased seismic velocity, is bounded to the northwest by a subvertical resistivity anomaly located near the Kirkland Lake and Cobalt kimberlite fields. The enhanced conductivity in the lower lithosphere is attributed to refertilization by fluids associated with Cretaceous kimberlite magmatism and can be explained by water content in olivine of 50 wt ppm in background areas with higher values in a localized anomaly beneath the kimberlite fields. Farther to the southeast the resistivity models include a lithospheric conductor between 100 and 150 km depth beneath the Central Metasedimentary Belt. The enhanced conductivity is attributed to grain boundary graphite films, associated with the Cretaceous kimberlitic magmatic process, or to water and carbon, introduced into the mantle during the pre-Grenvillian tectonism.

Description: Using dynamic rupture models of a right-lateral fault embedded in an elastic or elastoplastic 3-D medium, we investigate elastic and inelastic responses of compliant fault zones to nearby earthquake ruptures. We particularly examine effects of fault zone depth, width, shape and rigidity reduction on the surface displacement field. Our results from elastic models show that deeper and wider fault zones generally result in larger residual displacements. However, for shallow fault zones, the vertical residual displacement is insensitive to or even decreases with fault zone width. The width of horizontal displacement anomalies across a fault zone is only indicative of the fault zone width near the Earth's surface. There are trade-off effects among fault zone depth, width, shape and rigidity reduction on the amplitude of surface residual displacements. Our elastoplastic models show that plastic strain can occur along the entire fault zone near the Earth's surface and in the extensional quadrant at depth, if fault zone rocks are close to failure before a nearby earthquake happens. Compared with results from elastic models, plastic strain near the Earth's surface generally enhances surface displacements of a fault zone and does not change the trend of effects of fault zone depth and width, while plastic strain at depth can result in reduced retrograde motion or sympathetic motion across the fault zone, and introduce complexities in effects of fault zone depth and width. Sympathetic horizontal motion more likely occurs across a narrow fault zone with inelastic response at depth. Vertical motion in the extensional quadrant may actually decrease with fault zone width in elastoplastic models. Sympathetic horizontal motion, or small retrograde horizontal motion in conjunction with large vertical motion above a fault zone is indicative of inelastic response of a fault zone at depth.

Description: The Algerian margin formed through back-arc opening of the Algerian basin (Mediterranean Sea) resulting from the roll-back of the Tethyan slab. Recent geophysical data acquired along the Algerian margin showed evidence of active or recent compressive deformation in the basin due to the ongoing Africa–Eurasia convergence. Published data from four wide-angle seismic profiles have allowed imaging the deep structure of the Algerian margin and its adjacent basins. In this study, we converted these velocity models into density models, then into isostatic anomalies. This allowed us to image an isostatic disequilibrium (relative to a local isostasy model) reaching a maximum amplitude at the margin toe. Converting isostatic anomalies into Moho depth variations shows that the Moho extracted from wide-angle seismic data is deeper than the one predicted by a local isostasy model in the oceanic domain, and shallower than it in the continental domain. These anomalies can be interpreted by opposite flexures of two plates separated by a plate boundary located close to the margin toe. We use a finite element model to simulate the lithospheric flexure. The amplitude of the equivalent vertical Moho deflection is larger in the central part of the study area (6–7 km) than on the easternmost and westernmost profiles (3 km). The effective elastic thickness used to best match the computed deflection is always extremely low (always less than 10 km) and probably reflects the relatively low strength of the lithosphere close to the plate boundary. Comparison with other wide-angle seismic profiles across an active and a passive margin show that the North Algerian margin displays isostatic anomalies close to that of an active margin. Finally, plate flexure is highest at the southern tip of the ocean-continent transition, possibly indicating that a former passive margin detachment is reactivated as a crustal scale reverse fault pre-dating a future subduction.

Description: Wide-angle reflection/refraction seismic profiles were recorded across the Cyprus Arc, the plate boundary between the African Plate and the Aegean–Anatolian microplate, from the Eratosthenes Seamount to the Hecataeus Rise immediately south of Cyprus. The resultant models were able to resolve detail of significant lateral velocity variations, though the deepest crust and Moho are not well resolved from the seismic data alone. Conclusions from the modelling suggest that (i) Eratosthenes Seamount consists of continental crust but exhibits a laterally variable velocity structure with a thicker middle crust and thinner lower crust to the northeast; (ii) the Hecataeus Rise has a thick sedimentary rock cover on an indeterminate crust (likely continental) and the crust is significantly thinner than Eratosthenes Seamount based on gravity modelling; (iii) high velocity basement blocks, coincident with highs in the magnetic field, occur in the deep water between Eratosthenes and Hecataeus, and are separated and bounded by deep low-velocity troughs and (iv) one of the high velocity blocks runs parallel to the Cyprus Arc, while the other two appear linked based on the magnetic data and run NW–SE, parallel to the margin of the Hecataeus Rise. The high velocity block beneath the edge of Eratosthenes Seamount is interpreted as an older magmatic intrusion while the linked high velocity blocks along Hecataeus Rise are interpreted as deformed remnant Tethyan oceanic crust or mafic intrusives from the NNW–SSE oriented transform margin marking the northern boundary of Eratosthenes Seamount. Eratosthenes Seamount, the northwestern limit of rifted continental crust from the Levant Margin, is part of a jagged rifted margin transected by transform faults on the northern edge of the lower African Plate that is being obliquely subducted under the Aegean–Anatolian upper plate. The thicker crust of Eratosthenes Seamount may be acting as an asperity on the subducting slab, locally locking up subduction of the Cyprus Arc on its northern margin, while deformed Tethyan oceanic crust remains trapped between its northeastern margin and the Hecataeus Rise.