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: We describe a multihomogeneity theory for source-parameter estimation of potential fields. Similar to what happens for random source models, where the monofractal scaling-law has been generalized into a multifractal law, we propose to generalize the homogeneity law into a multihomogeneity law. This allows a theoretically correct approach to study real-world potential fields, which are inhomogeneous and so do not show scale invariance, except in the asymptotic regions (very near to or very far from their sources). Since the scaling properties of inhomogeneous fields change with the scale of observation, we show that they may be better studied at a set of scales than at a single scale and that a multihomogeneous model is needed to explain its complex scaling behaviour. In order to perform this task, we first introduce fractional-degree homogeneous fields, to show that: (i) homogeneous potential fields may have fractional or integer degree; (ii) the source-distributions for a fractional-degree are not confined in a bounded region, similarly to some integer-degree models, such as the infinite line mass and (iii) differently from the integer-degree case, the fractional-degree source distributions are no longer uniform density functions. Using this enlarged set of homogeneous fields, real-world anomaly fields are studied at different scales, by a simple search, at any local window W , for the best homogeneous field of either integer or fractional-degree, this yielding a multiscale set of local homogeneity-degrees and depth estimations which we call multihomogeneous model. It is so defined a new technique of source parameter estimation (Multi-HOmogeneity Depth Estimation, MHODE), permitting retrieval of the source parameters of complex sources. We test the method with inhomogeneous fields of finite sources, such as faults or cylinders, and show its effectiveness also in a real-case example. These applications show the usefulness of the new concepts, multihomogeneity and fractional homogeneity-degree, to obtain valid estimates of the source parameters in a consistent theoretical framework, so overcoming the limitations imposed by global-homogeneity to widespread methods, such as Euler deconvolution.

Description: In this paper, two separate but related goals are tackled. The first one is to demonstrate that in some saturated rock textures the non-linear behaviour of induced polarization (IP) and the violation of Ohm's law not only are real phenomena, but they can also be satisfactorily predicted by a suitable physical-mathematical model, which is our second goal. This model is based on Fick's second law. As the model links the specific dependence of resistivity and chargeability of a laboratory sample to the injected current and this in turn to its pore size distribution, it is able to predict pore size distribution from laboratory measurements, in good agreement with mercury injection capillary pressure test results. This fact opens up the possibility for hydrogeophysical applications on a macro scale. Mathematical modelling shows that the chargeability acquired in the field under normal conditions, that is at low current, will always be very small and approximately proportional to the applied current. A suitable field test site for demonstrating the possible reliance of both resistivity and chargeability on current was selected and a specific measuring strategy was established. Two data sets were acquired using different injected current strengths, while keeping the charging time constant. Observed variations of resistivity and chargeability are in agreement with those predicted by the mathematical model. These field test data should however be considered preliminary. If confirmed by further evidence, these facts may lead to changing the procedure of acquiring field measurements in future, and perhaps may encourage the design and building of a new specific geo-resistivity meter. This paper also shows that the well-known Marshall and Madden's equations based on Fick's law cannot be solved without specific boundary conditions.

Description: We present a new 3-D traveltime tomography code (TOMO3D) for the modelling of active-source seismic data that uses the arrival times of both refracted and reflected seismic phases to derive the velocity distribution and the geometry of reflecting boundaries in the subsurface. This code is based on its popular 2-D version TOMO2D from which it inherited the methods to solve the forward and inverse problems. The traveltime calculations are done using a hybrid ray-tracing technique combining the graph and bending methods. The LSQR algorithm is used to perform the iterative regularized inversion to improve the initial velocity and depth models. In order to cope with an increased computational demand due to the incorporation of the third dimension, the forward problem solver, which takes most of the run time (~90 per cent in the test presented here), has been parallelized with a combination of multi-processing and message passing interface standards. This parallelization distributes the ray-tracing and traveltime calculations among available computational resources. The code's performance is illustrated with a realistic synthetic example, including a checkerboard anomaly and two reflectors, which simulates the geometry of a subduction zone. The code is designed to invert for a single reflector at a time. A data-driven layer-stripping strategy is proposed for cases involving multiple reflectors, and it is tested for the successive inversion of the two reflectors. Layers are bound by consecutive reflectors, and an initial velocity model for each inversion step incorporates the results from previous steps. This strategy poses simpler inversion problems at each step, allowing the recovery of strong velocity discontinuities that would otherwise be smoothened.

Description: Iterative substitution of the coupled Marchenko equations is a novel methodology to retrieve the Green's functions from a source or receiver array at an acquisition surface to an arbitrary location in an acoustic medium. The methodology requires as input the single-sided reflection response at the acquisition surface and an initial focusing function, being the time-reversed direct wavefield from the acquisition surface to a specified location in the subsurface. We express the iterative scheme that is applied by this methodology explicitly as the successive actions of various linear operators, acting on an initial focusing function. These operators involve multidimensional crosscorrelations with the reflection data and truncations in time. We offer physical interpretations of the multidimensional crosscorrelations by subtracting traveltimes along common ray paths at the stationary points of the underlying integrals. This provides a clear understanding of how individual events are retrieved by the scheme. Our interpretation also exposes some of the scheme's limitations in terms of what can be retrieved in case of a finite recording aperture. Green's function retrieval is only successful if the relevant stationary points are sampled. As a consequence, internal multiples can only be retrieved at a subsurface location with a particular ray parameter if this location is illuminated by the direct wavefield with this specific ray parameter. Several assumptions are required to solve the Marchenko equations. We show that these assumptions are not always satisfied in arbitrary heterogeneous media, which can result in incomplete Green's function retrieval and the emergence of artefacts. Despite these limitations, accurate Green's functions can often be retrieved by the iterative scheme, which is highly relevant for seismic imaging and inversion of internal multiple reflections.

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: Employing dynamic reciprocity can be an effective tool to simplify the calculation of elastic wavefields for borehole problems and to check the results. We analytically obtain the reciprocity relations for the elastodynamic fields generated by multipole sources in a fluid–solid configuration: if the multipole sources are located in the fluid, the particle displacement due to a dipole source is reciprocal to the particle acceleration due to a single force; the fluid pressure due to a dipole source is reciprocal to the particle acceleration due to a monopole source; the particle displacement due to a quadrupole source is reciprocal to the spatial partial derivative of the particle acceleration due to a single force; the fluid pressure due to a quadrupole source is reciprocal to the spatial partial derivative of the particle acceleration due to a monopole source. These relations are tested by numerical experiments for different borehole problems, including acoustic logging, single-well imaging and vertical seismic profiling. A reciprocity test can be used as a quick check of a finite-difference algorithm and the implementation of the sources, although it cannot detect errors due to improper discretization of the interfaces.

Description: The area of the 9.1-km-deep Continental Deep Drillhole (KTB) in Germany is used as a case study for a geothermal reservoir situated in folded and faulted metamorphic crystalline crust. The presented approach is based on the analysis of 3-D seismic reflection data combined with borehole data and hydrothermal numerical modelling. The KTB location exemplarily contains all elements that make seismic prospecting in crystalline environment often more difficult than in sedimentary units, basically complicated tectonics and fracturing and low-coherent strata. In a first step major rock units including two known nearly parallel fault zones are identified down to a depth of 12 km. These units form the basis of a gridded 3-D numerical model for investigating temperature and fluid flow. Conductive and advective heat transport takes place mainly in a metamorphic block composed of gneisses and metabasites that show considerable differences in thermal conductivity and heat production. Therefore, in a second step, the structure of this unit is investigated by seismic waveform modelling. The third step of interpretation consists of applying wavenumber filtering and log-Gabor-filtering for locating fractures. Since fracture networks are the major fluid pathways in the crystalline, we associate the fracture density distribution with distributions of relative porosity and permeability that can be calibrated by logging data and forward modelling of the temperature field. The resulting permeability distribution shows values between 10 –16 and 10 –19 m 2 and does not correlate with particular rock units. Once thermohydraulic rock properties are attributed to the numerical model, the differential equations for heat and fluid transport in porous media are solved numerically based on a finite difference approach. The hydraulic potential caused by topography and a heat flux of 54 mW m –2 were applied as boundary conditions at the top and bottom of the model. Fluid flow is generally slow and mainly occurring within the two fault zones. Thus, our model confirms the previous finding that diffusive heat transport is the dominant process at the KTB site. Fitting the observed temperature–depth profile requires a correction for palaeoclimate of about 4 K at 1 km depth. Modelled and observed temperature data fit well within 0.2 °C bounds. Whereas thermal conditions are suitable for geothermal energy production, hydraulic conditions are unfavourable without engineered stimulation.

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: This paper presents the results from the Deflo-hydroacoustic experiment in the Southern Indian Ocean using three autonomous underwater hydrophones, complemented by two permanent hydroacoustic stations. The array monitored for 14 months, from November 2006 to December 2007, a 3000 x 3000 km wide area, encompassing large segments of the three Indian spreading ridges that meet at the Indian Triple Junction. A catalogue of 11 105 acoustic events is derived from the recorded data, of which 55 per cent are located from three hydrophones, 38 per cent from 4, 6 per cent from five and less than 1 per cent by six hydrophones. From a comparison with land-based seismic catalogues, the smallest detected earthquakes are m b 2.6 in size, the range of recorded magnitudes is about twice that of land-based networks and the number of detected events is 5–16 times larger. Seismicity patterns vary between the three spreading ridges, with activity mainly focused on transform faults along the fast spreading Southeast Indian Ridge and more evenly distributed along spreading segments and transforms on the slow spreading Central and ultra-slow spreading Southwest Indian ridges; the Central Indian Ridge is the most active of the three with an average of 1.9 events/100 km/month. Along the Sunda Trench, acoustic events mostly radiate from the inner wall of the trench and show a 200-km-long seismic gap between 2 °S and the Equator. The array also detected more than 3600 cryogenic events, with different seasonal trends observed for events from the Antarctic margin, compared to those from drifting icebergs at lower (up to 50°S) latitudes. Vocalizations of five species and subspecies of large baleen whales were also observed and exhibit clear seasonal variability. On the three autonomous hydrophones, whale vocalizations dominate sound levels in the 20–30 and 100 Hz frequency bands, whereas earthquakes and ice tremor are a dominant source of ambient sound at frequencies 〈20 Hz.

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: Seismic reflections from the oceanic water column contain information about ocean temperature and salinity. Even though seismic waveform inversion is effective for studying oceanic structure, its application is limited in the absence of sufficient direct temperature/velocity measurements. Here, two methods are developed to invert pre-stack seismic waveform data for temperature and salinity when independent temperature/velocity data are sparse or unavailable, allowing estimation of water-column temperature/salinity from any marine seismic reflection data set. The first method combines a genetic algorithm (GA) with non-linear least squares inversion, and the second method is a parallel implementation of a GA. Both methods produce results to an accuracy between 0 and 0.1 °C in estimating temperature when applied to a field data set from the South China Sea. Although the second approach is superior, it is computationally demanding and requires large parallel computers. The first approach runs extremely fast on parallel computers and can even be run on much smaller machines to provide results in a reasonable runtime. While both methods are viable choices for estimating temperature and salinity, the choice of one over the other will largely depend upon the available computational resources and the time frame within which the inversion needs to be completed.

Description: A new approach of seismoelectric imaging has been recently proposed to detect saturation fronts in which seismic waves are focused in the subsurface to scan its heterogeneous nature and determine saturation fronts. Such type of imaging requires however a complete modelling of the seismoelectric properties of porous media saturated by two immiscible fluid phases, one being usually electrically insulating (for instance water and oil). We combine an extension of Biot dynamic theory, valid for porous media containing two immiscible Newtonian fluids, with an extension of the electrokinetic theory based on the notion of effective volumetric charge densities dragged by the flow of each fluid phase. These effective charge densities can be related directly to the permeability and saturation of each fluid phase. The coupled partial differential equations are solved with the finite element method. We also derive analytically the transfer function connecting the macroscopic electrical field to the acceleration of the fast P wave (coseismic electrical field) and we study the influence of the water content on this coupling. We observe that the amplitude of the co-seismic electrical disturbance is very sensitive to the water content with an increase in amplitude with water saturation. We also investigate the seismoelectric conversions (interface effect) occurring at the water table. We show that the conversion response at the water table can be identifiable only when the saturation contrasts between the vadose and saturated zones are sharp enough. A relatively dry vadose zone represents the best condition to identify the water table through seismoelectric measurements. Indeed, in this case, the coseismic electrical disturbances are vanishingly small compared to the seismoelectric interface response.

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: Rock permeability is an important parameter for the formation evaluation. In this paper, a new method with streaming current is proposed to determine the sample permeability based on the electrokinetic effects, and is proved by the experimental measurements. Corresponding to this method, we have designed an experimental setup and a test system, then performed the streaming current (potential) and electro-osmosis pressure experiments with 23 sandstone samples at 0.05 mol l –1 NaCl solution. The streaming current (potential) coefficient and electro-osmosis pressure coefficient are obtained, respectively, with the experimental data at low frequencies with AC lock-in technique. The electrokinetic permeabilities are further calculated with these coefficients. The results are consistent well with the gas permeability measured with Darcy's law, which verifies the current method for estimating rock permeability. Our measurements are also analysed and compared with previous measurements. The results indicate that our method can reflect the essence of electrokinetic effects better and simplify the electrokinetic measurements as well. In addition, we discuss the influences of experimental artefacts (core holder and confining pressure installation) on the electrokinetic data. The results show that the trough phenomenon, appeared in frequency curves of streaming current (potential) coefficients, is induced by the resonance of the core-holder/vibrator system. This is important for the design of electrokinetic setup and the analysis of low-frequency response of the electrokinetic coupling coefficients.

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: Multiparameter full waveform inversion (FWI) applied to an elastic orthorhombic model description of the subsurface requires in theory a nine-parameter representation of each pixel of the model. Even with optimal acquisition on the Earth surface that includes large offsets, full azimuth, and multicomponent sensors, the potential for trade-off between the elastic orthorhombic parameters are large. The first step to understanding such trade-off is analysing the scattering potential of each parameter, and specifically, its scattering radiation patterns. We investigate such radiation patterns for diffraction and for scattering from a horizontal reflector considering a background isotropic model. The radiation patterns show considerable potential for trade-off between the parameters and the potentially limited resolution in their recovery. The radiation patterns of C 11 , C 22 , and C 33 are well separated so that we expect to recover these parameters with limited trade-offs. However, the resolution of their recovery represented by recovered range of model wavenumbers varies between these parameters. We can only invert for the short wavelength components (reflection) of C 33 while we can mainly invert for the long wavelength components (transmission) of the elastic coefficients C 11 and C 22 if we have large enough offsets. The elastic coefficients C 13 , C 23 , and C 12 suffer from strong trade-offs with C 55 , C 44 , and C 66 , respectively. The trade-offs between C 13 and C 55 , as well as C 23 and C 44 , can be partially mitigated if we acquire P – SV and SV – SV waves. However, to reduce the trade-offs between C 12 and C 66 , we require credible SH – SH waves. The analytical radiation patterns of the elastic constants are supported by numerical gradients of these parameters.

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: Typically, seismic data are sparsely and irregularly sampled due to limitations in the survey environment and these cause problems for key seismic processing steps such as surface-related multiple elimination or wave-equation-based migration. Various interpolation techniques have been developed to alleviate the problems caused by sparse and irregular sampling. Among many interpolation techniques, matching pursuit interpolation is a robust tool to interpolate the regularly sampled data with large receiver separation such as crossline data in marine seismic acquisition when both pressure and particle velocity data are used. Multicomponent matching pursuit methods generally used the sinusoidal basis function, which have shown to be effective for interpolating multicomponent marine seismic data in the crossline direction. In this paper, we report the use of wavelet basis functions which further enhances the performance of matching pursuit methods for de-aliasing than sinusoidal basis functions. We also found that the range of the peak wavenumber of the wavelet is critical to the stability of the interpolation results and the de-aliasing performance and that the range should be determined based on Nyquist criteria. In addition, we reduced the computational cost by adopting the inner product of the wavelet and the input data to find the parameters of the wavelet basis function instead of using L-2 norm minimization. Using synthetic data, we illustrate that for aliased data, wavelet-based matching pursuit interpolation yields more stable results than sinusoidal function-based one when we use not only pressure data only but also both pressure and particle velocity together.

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: Deep sea pockmarks underlain by chimney-like or pipe structures that contain methane hydrate are abundant along the Norwegian continental margin. In such hydrate provinces the interaction between hydrate formation and fluid flow has significance for benthic ecosystems and possibly climate change. The Nyegga region, situated on the western Norwegian continental slope, is characterized by an extensive pockmark field known to accommodate substantial methane gas hydrate deposits. The aim of this study is to detect and delineate both the gas hydrate and free gas reservoirs at one of Nyegga's pockmarks. In 2012, a marine controlled-source electromagnetic (CSEM) survey was performed at a pockmark in this region, where high-resolution 3-D seismic data were previously collected in 2006. 2-D CSEM inversions were computed using the data acquired by ocean bottom electrical field receivers. Our results, derived from unconstrained and seismically constrained CSEM inversions, suggest the presence of two distinctive resistivity anomalies beneath the pockmark: a shallow vertical anomaly at the underlying pipe structure, likely due to gas hydrate accumulation, and a laterally extensive anomaly attributed to a free gas zone below the base of the gas hydrate stability zone. This work contributes to a robust characterization of gas hydrate deposits within subseafloor fluid flow pipe structures.

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: One of the main applications of magnetic field measurements in boreholes is the detection of unexploded ordnance or buried utility structures like pipes or tiebacks. Even though the advantage of fully oriented magnetic vector measurements have long been recognized and could significantly reduce costs and risks, the tools used for those purposes typically measure only the total magnetic field, the vertical and horizontal components or gradients thereof. The Göttingen Bohrloch Magnetometer uses three fibre optic gyros to record its orientation and thus enables us to compute high-quality three-component magnetic vector data regardless of borehole orientation. The measurements described in this paper were run in the scientific borehole Cuxhaven Lüdingworth 1/1A, which was drilled as a part of the ‘Coastal Aquifer Test field’ project to study the dynamics of the saltwater/freshwater interface. As the drill string got stuck during drilling of the first borehole, a second hole was drilled in the immediate vicinity. The drill string lies at a depth between 80 and 114 m at a distance of only 2.5 m southeast of the borehole used for the measurements, making it an ideal target to demonstrate the benefits of vector magnetic surveys. Although the theories to calculate magnetic fields of objects with different shapes is well established and do not need to be tested, they almost exclusively include approximations of the geometry. It is not obvious whether these approximations are suited to describe real data, or whether additional effects or refinements have to be considered. We use both a simplified monopole model and a cylinder model to fit the data and are able to determine the position of the drill string within a statistical error of approximately 10 cm. Additionally, we show that the location of the drill string could not have been determined by measurements of the total field or its horizontal and vertical component alone and that those methods would require the drilling of additional boreholes to obtain an unambiguous result.

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 analysed four newly retrieved tide gauge records of the 1998 July 17 Papua New Guinea (PNG) tsunami to study statistical and spectral properties of this tsunami. The four tide gauge records were from Lombrum (PNG), Rabaul (PNG), Malakal Island (Palau) and Yap Island (State of Yap) stations located 600–1450 km from the source. The tsunami registered a maximum trough-to-crest wave height of 3–9 cm at these gauges. Spectral analysis showed two dominant peaks at period bands of 2–4 and 6–20 min with a clear separation at the period of ~5 min. We interpreted these peak periods as belonging to the landslide and earthquake sources of the PNG tsunami, respectively. Analysis of the tsunami waveforms revealed 12–17 min delay in landslide generation compared to the origin time of the main shock. Numerical simulations including this delay fairly reproduced the observed tide gauge records. This is the first direct evidence of the delayed landslide source of the 1998 PNG tsunami which was previously indirectly estimated from acoustic T-phase records.

Description: For extensional wave propagation along a cylindrical pipe, there exists a natural stopband in the frequency range between the first and second modes. This study explores the feasibility and practicality of building a drill collar acoustic extensional-wave isolator by combining the stopbands of pipes of different thicknesses. Numerical modelling shows that this is indeed possible and a stopband of designated width can be obtained using an optimization procedure. Laboratory measurement on an optimized design further verified this concept. The result provides a viable approach for the acoustic isolation design of a logging while drilling acoustic tool.

Description: The Murray Ridge/Dalrymple Trough system forms the boundary between the Indian and Arabian plates in the northern Arabian Sea. Geodetic constraints from the surrounding continents suggest that this plate boundary is undergoing oblique extension at a rate of a few millimetres per year. We present wide-angle seismic data that constrains the composition of the Ridge and of adjacent lithosphere beneath the Indus Fan. We infer that Murray Ridge, like the adjacent Dalrymple Trough, is underlain by continental crust, while a thin crustal section beneath the Indus Fan represents thinned continental crust or exhumed serpentinized mantle that forms part of a magma-poor rifted margin. Changes in crustal structure across the Murray Ridge and Dalrymple Trough can explain short-wavelength gravity anomalies, but a long-wavelength anomaly must be attributed to deeper density contrasts that may result from a large age contrast across the plate boundary. The origin of this fragment of continental crust remains enigmatic, but the presence of basement fabrics to the south that are roughly parallel to Murray Ridge suggests that it separated from the India/Seychelles/Madagascar block by extension during early breakup of Gondwana.

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: In conventional marine seismic surveys, due to the time-delayed reflections from the sea surface on both source and receiver sides, ghosts present in recorded seismograms and lead to both phase and spectrum distortions (especially near certain frequency notches). To achieve a high-quality broad-band image/velocity model with conventional reverse time migration (RTM)/full waveform inversion (FWI) that adopts a synthetic zero-phase source wavelet and absorbing surface condition during wavefield modelling, marine seismic data have to be pre-processed to remove ghost effects. However, seismic deghosting is not a trivial task. Instead of employing an external deghosting process, we propose a strategy to compensate for ghost effects during FWI and RTM, which consists of two parts: first, to address phase distortions due to ghost effects by means of obtaining an accurate source wavelet estimation and adopting an appropriate surface boundary condition in both forward and backward wave propagation to appropriately generate ghosts; secondly, to build a compensation operator in the adjoint state computation to mitigate spectrum distortions caused by dominant ghost effects. To demonstrate the success and robustness of the proposed strategy, we present both synthetic experiments and field examples, which suggest that this strategy can lead to successful applications of FWI/RTM directly on marine seismic data without an extra deghosting process.

Description: Seismic full waveform inversion (FWI) has been applied to simple elastic problems with certain symmetries, such as isotropic, transverse isotropic or vertical transversely isotropic media. In this study, the FWI concept is extended to the most general anisotropic case with 21 independent elastic material parameters and no symmetry plane (triclinic). Beside a short description of the 3-D finite-difference scheme to solve the forward problem and the FWI optimization algorithm, we present a sensitivity study for a simple anisotropic medium. This test problem consists of a homogenous triclinic anisotropic full space, which contains 21 spatially separated spheres. In each sphere one component of the elastic tensor deviates by 5 per cent from the background medium. The resolution of the different spheres, ambiguities between the different elastic parameters, as well as the effect of the acquisition geometry can be systematically investigated. Due to the high computational costs of the triclinic forward problem a few compromises have to be made regarding the acquisition geometries. Point sources are replaced by plane wave sources which lead to a limitation of incidence angles and therefore a strong decrease in resolution of the nondiagonal elastic tensor components. It is shown that, despite these limitations, a tomographic acquisition geometry would be able to resolve to some extent a monoclinic symmetry via FWI. Restricting the acquisition geometries (e.g. VSP combined with reflection seismic or reflection seismic only) significantly reduces the number of resolvable tensor elements in strict dependence of the covered incidence angles.

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: Non-linear resonant coupling of edge waves can occur with tsunamis generated by large-magnitude subduction zone earthquakes. Earthquake rupture zones that straddle beneath the coastline of continental margins are particularly efficient at generating tsunami edge waves. Using a stochastic model for earthquake slip, it is shown that a wide range of edge-wave modes and wavenumbers can be excited, depending on the variability of slip. If two modes are present that satisfy resonance conditions, then a third mode can gradually increase in amplitude over time, even if the earthquake did not originally excite that edge-wave mode. These three edge waves form a resonant triad that can cause unexpected variations in tsunami amplitude long after the first arrival. An M ~ 9, 1100 km-long continental subduction zone earthquake is considered as a test case. For the least-variable slip examined involving a Gaussian random variable, the dominant resonant triad includes a high-amplitude fundamental mode wave with wavenumber associated with the along-strike dimension of rupture. The two other waves that make up this triad include subharmonic waves, one of fundamental mode and the other of mode 2 or 3. For the most variable slip examined involving a Cauchy-distributed random variable, the dominant triads involve higher wavenumbers and modes because subevents, rather than the overall rupture dimension, control the excitation of edge waves. Calculation of the resonant period for energy transfer determines which cases resonant coupling may be instrumentally observed. For low-mode triads, the maximum transfer of energy occurs approximately 20–30 wave periods after the first arrival and thus may be observed prior to the tsunami coda being completely attenuated. Therefore, under certain circumstances the necessary ingredients for resonant coupling of tsunami edge waves exist, indicating that resonant triads may be observable and implicated in late, large-amplitude tsunami arrivals.

Description: Although seismic attenuation measurements have great potential to enhance our knowledge of physical conditions and rock properties, their application is limited because robust methods for improving both the resolution and accuracy of attenuation estimates have not yet been established. We propose attenuation estimation methods for zero-offset vertical seismic profile (VSP) data by combining seismic interferometry (SI) and the modified median frequency shift (MMFS) method developed for attenuation estimation using sonic waveform data. The configuration of zero-offset VSP data is redatumed to that of the sonic logging measurement by adopting two types of SI: deconvolution interferometry and crosscorrelation interferometry (CCI). Then, we can apply the MMFS method to the redatumed VSP data. Although the amplitude information estimated from CCI is biased, we propose a correction method for this bias to correctly estimate attenuation. First, to investigate the performance both in resolution and accuracy, we apply different trace separations to synthetic data with random noise at different signal-to-noise ratio levels. Second, we estimate the influence of residual reflection events after wavefield separation on attenuation estimation. The proposed methods provide more stable attenuation estimates in comparison with the spectral ratio method because the mean-median procedure suppresses random events and characteristic features caused by residual reflection events in spectral domain. Our numerical experiments also imply that the proposed methods can estimate scattering attenuation values even if frequency components are not changed between the two receiver depths. Finally, by preliminarily applying the proposed methods to field VSP data, we demonstrate the advantages of the proposed method in the resolution and stability of attenuation estimates and these observations correlate with those of numerical tests.

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: One interpretation of a seismic image is the instantaneous scattered wave response of a colocated pseudo-source and pseudo-receiver at each point in the subsurface model. If there is no model perturbation at a point then there will be no instantaneous scattered wave so nothing will be imaged; if something is imaged then there must be a perturbation at that location. By extension, so-called extended images (EIs) represent the full spatio-temporal response between offset subsurface pseudo-sources and pseudo-receivers which can be used to constrain elastic properties around each image point. However, one-sided illumination of the subsurface (from the Earth's surface), errors in the initial velocity model estimate, and the use of a linearized, single-scattering assumption (as is usual in seismic imaging) cause errors in EI gathers such as missing events, incorrect amplitudes, and spurious energy. By creating elastic (P-to-P and P-to-S) EIs in a synthetic example of subsalt imaging, we demonstrate the advantages of incorporating multiply scattered waves correctly by non-linear imaging, and of including transmitted waves by using two-sided receiver arrays, and discuss how the recently developed autofocussing methods could provide us with the various required subsurface wavefields. Pre- and post-imaging f–k filtering procedures are introduced to further improve the quality of the EIs by (explicitly or implicitly) limiting the directions of waves arriving at the subsurface pseudo-source and receiver survey line. These filters suppress strong linear events that arise from the erroneous interaction of near-horizontally propagating waves which are not naturally accounted for due to the lack of sources and receivers on either side of the imaging target. Finally, we analyse the sensitivity of elastic P-to-P EIs to errors in the migration velocity models and show that events in the EI are shifted in opposite directions when constructed using reflection or transmission data. In other words, velocity errors are mapped into the EIs differently in the case of one-sided from two-sided illumination. This leads to the potential for new methods of migration velocity analysis when surface and borehole seismic data are jointly acquired.

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: Uncertainty estimation is a vital part of geophysical numerical modelling. There exist a variety of methods aimed at uncertainty estimation, which are often complicated and difficult to implement. We present an inversion technique that produces multiple solutions, based on bootstrap resampling, to create a qualitative uncertainty measure for 2-D magnetotelluric inversion models. The approach is easy to implement, can be used with almost any inversion code, and does not require access to the inversion software's source code. It is capable of detecting the effect of data uncertainties on the model result rather than just analysing the effect of model variations on the model response. To obtain uncertainty estimates for an inversion model, the original data set is resampled repeatedly and alternate data set realizations are created and inverted. This ensemble of solutions is then statistically analysed to determine the variability between the different solutions. The process yields interpretable uncertainty maps for the inversion model and we demonstrate its effectiveness to qualitatively quantify uncertainty in synthetic model tests and a case study.

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: Electromagnetic (EM) imaging methods are useful tools for monitoring subsurface changes in pore-fluid content and the associated changes in electrical permittivity and conductivity. The most common method for georadar tomography uses a high frequency ray-theoretic approximation that is valid when material variations are sufficiently small relative to the wavelength of the propagating wave. Georadar methods, however, often utilize EM waves that propagate within heterogeneous media at frequencies where ray theory may not be applicable. In this paper we describe EM wave propagation 3-D Fréchet sensitivity kernels that capture the data sensitivity to material perturbations for a given source–receiver combination. Various data functional types are formulated that consider all three components of the electric wavefield and incorporate near-, intermediate- and far-field contributions. We show that EM waves exhibit substantial variations for different relative source–receiver component orientations. The 3-D sensitivities also illustrate out of plane effects that are not captured in 2-D sensitivity kernels and can influence results obtained using 2-D inversion methods to image structures that are in reality 3-D.

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: This paper focuses on the modelling of fluid-filled poroelastic double porosity media under quasi-static and dynamic regimes. The double porosity model is derived from a two-scale homogenization procedure, by considering a medium locally characterized by blocks of poroelastic Biot microporous matrix and a surrounding system of fluid-filled macropores or fractures. The derived double porosity description is a two-pressure field poroelastic model with memory and viscoelastic effects. These effects result from the ‘time-dependent’ interaction between the pressure fields in the two pore networks. It is shown that this homogenized double porosity behaviour arises when the characteristic time of consolidation in the microporous domain is of the same order of magnitude as the macroscopic characteristic time of transient regime. Conversely, single porosity behaviours occur when both timescales are clearly distinct. Moreover, it is established that the phenomenological approaches that postulate the coexistence of two pressure fields in ‘instantaneous’ interaction only describe media with two pore networks separated by an interface flow barrier. Hence, they fail at predicting and reproducing the behaviour of usual double porosity media. Finally, the results are illustrated for the case of stratified media.

Description: The Okavango Delta of northern Botswana is one of the world's largest inland deltas or megafans. To obtain information on the character of sediments and basement depths, audiomagnetotelluric (AMT), controlled-source audiomagnetotelluric (CSAMT) and central-loop transient electromagnetic (TEM) data were collected on the largest island within the delta. The data were inverted individually and jointly for 1-D models of electric resistivity. Distortion effects in the AMT and CSAMT data were accounted for by including galvanic distortion tensors as free parameters in the inversions. By employing Marquardt–Levenberg inversion, we found that a 3-layer model comprising a resistive layer overlying sequentially a conductive layer and a deeper resistive layer was sufficient to explain all of the electromagnetic data. However, the top of the basal resistive layer from electromagnetic-only inversions was much shallower than the well-determined basement depth observed in high-quality seismic reflection images and seismic refraction velocity tomograms. To resolve this discrepancy, we jointly inverted the electromagnetic data for 4-layer models by including seismic depths to an interface between sedimentary units and to basement as explicit a priori constraints. We have also estimated the interconnected porosities, clay contents and pore-fluid resistivities of the sedimentary units from their electrical resistivities and seismic P -wave velocities using appropriate petrophysical models. In the interpretation of our preferred model, a shallow ~40 m thick freshwater sandy aquifer with 85–100 m resistivity, 10–32 per cent interconnected porosity and 〈13 per cent clay content overlies a 105–115 m thick conductive sequence of clay and intercalated salt-water-saturated sands with 15–20 m total resistivity, 1–27 per cent interconnected porosity and 15–60 per cent clay content. A third ~60 m thick sandy layer with 40–50 m resistivity, 10–33 per cent interconnected porosity and 〈15 per cent clay content is underlain by the basement with 3200–4000 m total resistivity. According to an interpretation of helicopter TEM data that cover the entire Okavango Delta and borehole logs, the second and third layers may represent lacustrine sediments from Paleo Lake Makgadikgadi and a moderately resistive freshwater aquifer comprising sediments of the recently proposed Paleo Okavango Megafan, respectively.

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: Interpolation and random noise removal is a pre-requisite for multichannel techniques because the irregularity and random noise in observed data can affect their performances. Projection Onto Convex Sets (POCS) method can better handle seismic data interpolation if the data's signal-to-noise ratio (SNR) is high, while it has difficulty in noisy situations because it inserts the noisy observed seismic data in each iteration. Weighted POCS method can weaken the noise effects, while the performance is affected by the choice of weight factors and is still unsatisfactory. Thus, a new weighted POCS method is derived through the Iterative Hard Threshold (IHT) view, and in order to eliminate random noise, a new adaptive method is proposed to achieve simultaneous seismic data interpolation and denoising based on dreamlet transform. Performances of the POCS method, the weighted POCS method and the proposed method are compared in simultaneous seismic data interpolation and denoising which demonstrate the validity of the proposed method. The recovered SNRs confirm that the proposed adaptive method is the most effective among the three methods. Numerical examples on synthetic and real data demonstrate the validity of the proposed adaptive method.

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: Capacitively Coupled Resistivity (CCR) is conventionally used to emulate DC resistivity measurements and may provide important information about the ice content of material in periglacial areas. The application of CCR theoretically enables the determination of both electrical parameters, that is, the resistivity and the electrical permittivity, by analysing magnitude and phase shift spectra. The electrical permittivity may dominate the impedance, especially in periglacial areas or regions of hydrogeological interest. However, previous theoretical work suggested that the phase shift may strongly depend on electrode height above ground, implying that electrode height must be known with great accuracy to determine electrical permittivity. Here, we demonstrate with laboratory test measurements, theoretical modelling and by analysing the Jacobian matrix of the inversion that the sensitivity towards electrode height is drastically reduced if the electrical permittivity is frequency dependent in a way that is typical for ice. For the first time, we used a novel broad-band CCR device ‘Chameleon’ for a field test located in one of the ridge galleries beneath the crest of Mount Zugspitze. A permanently ice covered bottom of a tunnel was examined. For the inversion of the measured spectra, the frequency dependence of the electrical parameters was parametrized in three different ways: A Debye Model for pure ices, a Cole–Cole Model for pure ices and a dual Cole–Cole Model including interfacial water additionally. The frequency-dependent resistivity and permittivity spectra obtained from the inversion, including low- and high-frequency limits, agree reasonably well with laboratory and field measurements reported in the literature.

Description: We evaluate the applicability and the effectiveness of GPR attribute analysis for high-resolution glacier imaging and characterization, testing this approach on 4-D GPR multifrequency data collected in a small glacier in the Eastern Alps, by repeating the acquisition along the same profiles in four different periods of the year 2013. The main objectives are to image and characterize the glacier's inner structure and to quantitatively monitor the seasonal thawing of near-surface frozen materials (snow/firn). A multiattribute approach is used to characterize the subsurface through different attribute categories, including instantaneous and textural attributes considering not only amplitude-, phase- and frequency-related attributes, but also other more complex and integrated parameters. We combine information from more than one attribute into a single image with composite displays, using overlays or mixed displays. The results demonstrate that the developed GPR attribute analysis can provide significant improvements in the discrimination of GPR signals, and obtain enhanced and more constrained data interpretations.

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: Migration velocity analysis aims at automatically updating the large-scale components of the velocity model, called macromodel. Extended Common Image Gathers are panels used to evaluate focusing after imaging and are constructed as a function of a spatial shift introduced in the imaging condition. We investigate how transmitted waves can also be used in migration velocity analysis: instead of back-propagating the residuals associated with reflected waves, we propose to back-propagate the full wavefield. The image function, equivalent to the migrated section for reflected data, does not exhibit localized events in space along horizons but is still sensitive to the choice of the background velocity model and can thus be coupled to the same objective function defined in the image domain. In order to enhance the benefits of direct waves, we consider a cross-well configuration. Direct waves provide a large illumination between two vertical wells. Associated Common Image Gathers present different characteristics than the ones associated with reflected waves in surface acquisition. In particular, energy is spread over up to the maximum penetration depth. We invert cross-well seismic data along two lines. In the first case, the input data contain the full wavefield dominated by transmitted waves. It demonstrates the possibility to handle transmitted waves to determine the velocity model. It appears that the misfit in the data domain is largely reduced after inversion. In the second case, we use the same algorithm, but with reflected observed data only, as in a classical approach. Most of velocity updates are localized around the reflectivity, leading to an incorrect final model. This demonstrates the benefit of transmitted waves for migration velocity analysis in a cross-well configuration.

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: An integrated approach to interpret Self-Potential (SP) anomalies based on spectral analysis and tomographic methods is presented. The Maximum Entropy Method (MEM) is used for providing accurate estimates of the depth of the anomaly source. The 2-D tomographic inversion technique, based on the underground charge occurrence probability (COP) function, is, then, used to fully characterize the anomalous body, as the MEM is not helpful in delineating the shape of the anomaly source. The proposed integrated approach is applied for the inversion of synthetic SP data generated by geometrically simple anomalous bodies, such as cylinders and inclined sheets. This numerical study has allowed the determination of mathematical relationships between zero lines of the COP distributions, the polarization angles and the positions along the profile of the causative sources, which have been of great help for interpreting the related SP anomalies. Finally, the analysis of field examples shows the high potential applicability of the proposed integrated approach for SP data inversion.

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: Volcanic and tectonic events are the main processes involved in the generation of the oceanic crust and responsible for the seismicity associated with seafloor spreading. To monitor this activity, usually not or poorly detected by land-based seismological stations, we deployed from February 2012 to February 2013 a network of autonomous hydrophones to compare the behaviour of the ultraslow-spreading Southwest Indian ridge (SWIR) with that of the intermediate-spreading Southeast Indian ridge (SEIR). The rate of seismicity is similar for both ridges, suggesting that there is no systematic relationship between seismicity and spreading rates. The along-axis distribution of the seismic events, however, does differ, reflecting the rate dependence of accretion modes. Earthquakes are sparse and regularly spaced and scattered along the SWIR, reflecting prevailing tectonic processes. By contrast, along the SEIR, events are irregularly distributed and focus at ridge-segment ends and transforms faults, reflecting the ridge segmentation; only two swarms occurred at a segment centre and are probably caused by a magmatic event. This seismicity distribution thus looks controlled by segment-scale crustal heterogeneities along the SEIR and by regional-scale contrasting accretion processes along the SWIR, probably driven by different lithospheric and asthenospheric dynamics on either side of the Melville fracture zone. The comparison of hydroacoustic and teleseismic catalogues shows that, along these spreading ridges, the background seismicity observed in 1 yr by a hydroacoustic network is representative of the seismicity observed over two decades by land-based networks.

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: Frequency-domain loop–loop electromagnetic (EM) methods are sensitive to the magnetic susceptibility of the Earth as well as its resistivity. Thus, inversion techniques have been used to simultaneously reconstruct both resistivity and susceptibility models from EM data. However, to take full advantage of inversion methods, calibration errors must be assessed and removed because ignoring them can result in misleading models. We present a multidimensional inversion method that jointly inverts EM and direct current (DC) resistivity data to derive offset errors as well as resistivity and susceptibility models, assuming that calibration errors can be represented by in-phase and quadrature offsets at each frequency. Addition of independent data such as DC data is effective for more accurately estimating the offsets, resulting in more reliable subsurface models. Synthetic examples involving small-loop EM data show that simultaneous inversion for resistivity and susceptibility is not stable, because of strong correlations between in-phase offset parameters and background susceptibility, but that the offsets are well determined when the data misfit is reduced rapidly in the early iteration step. Improvements achieved by joint inversion are mainly on the resistivity model. For airborne electromagnetic (AEM) data, the inversion process is stable, because AEM data are acquired using more loop–loop geometries and a wider range of frequencies. As a result, both the resistivity and susceptibility models are significantly improved by joint inversion.

Description: We have developed an efficient numerical scheme for 3-D electromagnetic (EM) simulations using an exponential finite-difference (FD) method with non-uniform grids. The method uses the set of exponential basis functions {1, exp [ ± ( x x + y y + z z )]}, where the exponents x , y and z must be chosen carefully depending on the simulation frequency and local node conductivity. The method achieves an approximation of the oscillatory and exponentially decaying EM fields that is better than that obtained via the low-degree polynomial fitting from standard FDs—and hence also leads to more accurate results. An important property of the exponential FD method is that it tends to the standard FD method when the exponents x , y and z tend to zero. We applied the standard and exponential FD methods to three marine controlled-source EM modelling scenarios: deep-water, shallow-water and intermediate water depth. For the deep-water scenario, we found that the proposed exponential FD method gave two to three times more accurate results as compared to the standard FD method on the same grid. For the shallow-water and intermediate water depth scenarios, the exponential FD method improved the accuracy of the upgoing fields; it gave 2–2.5 times more accurate results for the upgoing fields than the standard FD method on the same grid. Consequently, the method can achieve the same accuracy with a coarser grid and hence is faster than the standard FD method, as demonstrated using a frequency-domain iterative solver.

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: Rocks in nature are multiphase composites of minerals and fluids with varying elastic properties. Predicting the effects that geochemical and geological processes, which lead to removal, alteration or substitution of one or more phases, have on seismic velocities is of considerable practical interest in geophysics. We present a new theory that quantifies the effect of replacing ‘one' or ‘more' phases (fluid or solid) on the elastic properties of multimineralic isotropic rocks. These new exact solutions relax the assumption of unchanging rock microstructure upon substitution of the pore-filling material. By developing equations that do not depend on this assumption, which is fundamental to conventional substitution theories, for example, Biot–Gassmann theory, we expand their range of applicability. Since the new exact solutions require usually unknown parameters, we also develop approximations that depend only on initial effective stiffness, elastic properties and volume fractions of constituents.

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: Geothermal well GRT-1 (Rittershoffen, Alsace) was drilled in 2012. Its open-hole section (extending down to a depth of 2.6 km) penetrated fractured sandstones and granite. In 2013, the well was subjected to Thermal, Chemical and Hydraulic (TCH) stimulation, which improved the injectivity index fivefold. The goal of the study was to assess the impact of the stimulation by comparing pre- and post-stimulation well-logging (acoustic and temperature [ T ] logs) and mud-logging data. This comparison revealed modifications of almost all the natural fractures. However, not all of these fractures are associated with permeability enhancement, and the post-stimulation T logs are important for characterizing this enhancement. Chemical alteration due to mechanical erosion at the tops and bottoms of the fractures was observed in the sandstones. These zones display indications of very small new permeability after the TCH stimulation. Because a major fault zone caved extensively where it crosses the borehole, it was not imaged in the acoustic logs. However, this originally permeable zone was enhanced as demonstrated by the T logs. Based on the natural injectivity of this fault zone, hydraulic erosion and thermal microcracking of its internal quartz veins are associated with this permeability enhancement. Although local changes in the borehole wall observed in the acoustic images cannot be directly linked to the improved injectivity index, the comparison of the acoustic image logs allows for identification of fracture zones impacted by the TCH stimulation.

Description: We investigate the spatial coherence of underwater ambient noise using a yearlong time-series measured off Ascension Island. Qualitative agreement with observed cross-correlations is achieved using a simple range-dependent model, constrained by earlier, active tomographic studies in the area. In particular, the model correctly predicts the existence of two weakly dispersive normal modes in the microseism frequency range, with the group speed of one of the normal modes being smaller than the sound speed in water. The agreement justifies our interpretation of the peaks of the measured cross-correlation function of ambient noise as modal arrivals, with dispersion that is sensitive to crustal velocity structure. Our observations are consistent with Scholte to Moho head wave coupled propagation, with double mode conversion occurring due to the bathymetric variations between receivers. We thus demonstrate the feasibility of interrogating crustal properties using noise interferometry of moored hydrophone data at ranges in excess of 120 km.

Description: Recent tectonic reconstructions of the South Atlantic have partitioned the ocean basin into several segments based upon one or more proposed intraplate South American deformation zones. In several of these reconstructions, opening of the southern segment(s) by seafloor spreading prior to Aptian-Albian time is accompanied by contemporaneous strike-slip motion along an intraplate boundary extending southeastward from the Andean Cochabamba—Santa Cruz bend to the Rio Grande Fracture Zone (RGFZ). We have examined new magnetic data over the Pelotas, Santos and Campos Basins, offshore Argentina and Brazil, acquired by ION-GXT in tandem with long-offset, long record seismic reflection data, and identified seafloor spreading anomalies M4, M3, M2 and M0 (~131, ~129, ~128 and ~125 Ma). Integrating these results with our earlier work, we have been able to correlate magnetochrons M4, M3, M2 and M0 north and south of the RGFZ on the South American margin, and north and south of the Walvis Ridge on the African side. Our results are therefore inconsistent with diachronous opening models that involve substantial continental strike-slip motion north of RGFZ during M4 to M0 time. Although the ocean basin may have opened from south to north, our results indicate that seafloor spreading began north of the RGFZ earlier than previously proposed.

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: Determining thin layer thickness is very important for reservoir characterization and CO 2 quantification. Given its high time–frequency resolution and robustness, the complex spectral decomposition method was applied on time-lapse 3-D seismic data from the Ketzin pilot site for CO 2 storage to evaluate the frequency-dependent characteristics of thin layers at the injection level. Higher temporal resolution and more stratigraphic details are seen in the all-frequency and monochromatic reflectivity amplitude sections obtained by complex spectral decomposition compared to the stacked sections. The mapped geologic discontinuities within the reservoir are consistent with the preferred orientation of CO 2 propagation. Tuning frequency mapping shows the thicknesses of the reservoir sandstone and gaseous CO 2 is consistent with the measured thickness of the sandstone unit from well logging. An attempt to discriminate between pressure effects and CO 2 saturation using the extracted tuning frequency indicates that CO 2 saturation is the main contributor to the amplitude anomaly at the Ketzin site. On the basis of determined thickness of gaseous CO 2 in the reservoir, quantitative analysis of the amount of CO 2 was performed and shows a discrepancy between the injected and calculated CO 2 mass. This may be explained by several uncertainties, like structural reservoir heterogeneity, a limited understanding of the complex subsurface conditions, error of determined tuning frequency, the presence of ambient noise and ongoing CO 2 dissolution.

Description: Comprehensive analytical solutions to 3-D axisymmetric problems for static response of multilayered thermoelastic media subjected to surface loads and containing sources are presented in this study. The solution strategy employs Laplace and Hankel transforms to the field variables. The problem is formulated in cylindrical coordinate system and in this coordinate system vector surface harmonics and generalized propagator matrix are introduced to find the solution for the problem for the behaviour of thermoelastic multilayered media subject to surface loads and containing heat sources. A high-order adaptive Gaussian quadrature method with continued fraction expansions is employed to approximate the integral solutions expressed in terms of semi-infinite Hankel-type integrals. It is the first time to apply the proposed solution method to investigate the behaviour of repository for heat-emitting high-level nuclear waste (HLW) in a geological formation where the HLW can be regarded as a decaying with time point heat source.

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: We present 1D anisotropic inversion of magnetotelluric (MT) data as a potential tool for mapping structural permeability in sedimentary basins. Using 1D inversions of a 171 site, broadband MT data set from the Koroit region of the Otway Basin, Victoria, Australia, we have delineated an electrically anisotropic layer at approximately 2.5 to 3.5 km depth. The anisotropy strike is consistent between stations at approximately 160° east of north. The depth of anisotropy corresponds to the top depth of the Lower Cretaceous Crayfish Group, and the anisotropy factor increases from west to east. We interpret the anisotropy as resulting from north–northwest oriented, fluid-filled fractures resulting in enhanced electrical and hydraulic conductivity. This interpretation is consistent with permeability data from well formation tests. It is also consistent with the orientation of mapped faults in the area, which are optimally oriented for reactivation in the current stress field.

Description: A recent paper by Tarasov and Titov suggested that the Cole–Cole conductivity models should be preferred to the Pelton resistivity model in fitting induced polarization responses. Each model has four parameters: resistivity (or its inverse conductivity ), chargeability m , time parameter and frequency dependence c . Tarasov and Titov showed that in fitting experimental data, 3 of the parameters directly correspond across the two formulations, but that there is a difference between the Pelton and Cole–Cole model in that parameter is only the same at low frequency, but is an entangled function with m and c at high frequencies in the Pelton formulation. This claim of inconsistent was based on using the Pelton complex resistivity formulation and the Cole–Cole complex conductivity equation to analyse complex conductivity data as a function of frequency. However, if Pelton and Cole–Cole models are used to fit complex resistivity, rather than fitting conductivity, then it is the Cole–Cole model that has an entangled parameter at high chargeability m values. Simple testing shows that a Pelton model used to fit resistivity has all four model parameters directly corresponding with conductivity data fitted with a Cole–Cole model. The conclusion is that Cole–Cole and Pelton models in fact require identical independent parameters that can fit experimental or synthetic data, provided the Pelton resistivity formulation is restricted to fitting to resistivity data, and the Cole–Cole conductivity formulation is restricted to fitting conductivity results.

Description: Stratum quality factors ( P -wave and S -wave quality factors, Q p and Q s ) have gradually been utilized in the study of physical state of crust and uppermost mantle, tectonic evolution, hydrogeololgy, gas hydrates, petroleum exploration, etc. Different opinions of the seismic attenuation mechanism result in various approaches to estimate the P -wave and S -wave quality factors. Considering the viscoelasticity of the underground medium, the constitutive matrix of the Earth medium is written as the superposition of homogeneous background medium, elastic perturbation medium and viscoelastic perturbation medium. Under the hypothesis of Born integral and stationary phase approximation, the seismic reflectivity is initially raised in terms of P -wave and S -wave moduli, density, P -wave and S -wave quality factors. Furthermore, incorporating the complex seismic traces with the seismic wavelets at different offsets, a two-step inversion approach is proposed to estimate the P -wave and S -wave quality factors. The AVO/AVA Bayesian inversion approach is suggested to estimate the P -wave modulus and S -wave modulus with the real component of the pre-stack seismic data initially. Taking the estimated P -wave and S -wave moduli as prior information, the P -wave and S -wave quality factors are further estimated with the imaginary component of the complex pre-stack seismic data, which is the quadrature of the original data. Finally, synthetic examples demonstrate that the proposed approach is able to estimate P -wave and S -wave quality factors stably and properly, and two field data examples demonstrate that the proposed approach may work as an efficient approach to fluid identification.

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: In this paper, we modelled the electrical transport behaviour of bimodal carbonate rocks from a reservoir in China using dual-pore networks. One basic assumption, generally supported by experimental data and microstructure observations in the reservoir samples, was that the low porosity, monomodal rocks had the same properties and structure as the microporous matrix of the high porosity, bimodal samples. We assumed that the matrix was homogeneous and always interconnected but that the connectivity and the pore size distribution of macropore system was randomly variable. Both pore systems were supposed to act locally as ‘in parallel’ electrical conductors, an approach previously used by Bauer et al. Hence, the effect of matrix properties, macropore size distribution and connectivity on electrical properties of bimodal rocks could be modelled and investigated. We simulated electrical current through 3-D, simple cubic and body-centred cubic networks with different coordination numbers, different pipe radius distributions of macropore system and different matrix properties. The main result was that the formation factor of dual-pore network obeyed a ‘universal’ scaling relationship (i.e. independent of lattice type). Based on this result, we extended the power-law model derived by Bernabé et al. for monomodal porous media. We developed methods for evaluating the scale-invariant pore structure parameters in the model using conventional core analysis and satisfactorily tested the proposed model against experimental data from the Chinese reservoir as well as some other previously published data sets.