ALBERT

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

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

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

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

Description: The 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: 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: 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.