ALBERT

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

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

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

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

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

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

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

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

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