Subduction dynamics as revealed by trench migration

International audienceNew estimates of trench migration rates allow us to address the dynamics of trench migration and back-arc strain. We show that trench migration is primarily controlled by the subducting plate velocity V-sub, which largely depends on its age at the trench. Using the hot and weak arc to back-arc region as a strain sensor, we define neutral arcs characterized by the absence of significant strain, meaning places where the forces (slab pull, bending, and anchoring) almost balance along the interface between the plates. We show that neutral subduction zones satisfy the kinematic relation between trench and subducting plate absolute motions: V-t = 0.5V(sub) - 2.3 (in cm a(-1)) in the HS3 reference frame. Deformation occurs when the velocity combination deviates from kinematic equilibrium. Balancing the torque components of the forces acting at the trench indicates that stiff (old) subducting plates facilitate trench advance by resisting bending

Some Results on Sprout

Abstract. Sprout is a lightweight stream cipher proposed by Armknecht and Mikhalev at FSE 2015. It has a Grain-like structure with two State Registers of size 40 bits each, which is exactly half the state size of Grain v1. In spite of this, the cipher does not appear to lose in security against generic Time-Memory-Data Tradeoff attacks due to the novelty of its design. In this paper, we first present improved results on Key Recovery with partial knowledge of the internal state. We show that if 50 of the 80 bits of the internal state are guessed then the remaining bits along with the Secret Key can be found in a reasonable time using a SAT solver. Thereafter we show that it is possible to perform a distinguishing attack on the full Sprout stream cipher in the multiple IV setting using around 240 randomly chosen IVs on an average. The attack requires around 248 bits of memory. Thereafter we will show that for every Secret Key, there exist around 230 IVs for which the LFSR used in Sprout enters the all zero state during the Keystream generating phase. Using this observation, we will first show that it is possible to enumerate Key-IV pairs that produce keystream bits with period as small as 80. We will then outline a simple Key recovery attack that takes time equivalent to 266.7 encryptions with negligible memory requirement. This although is not the best attack reported against this cipher in terms of the Time complexity, it is the best in terms of the memory required to perform the attack

Link-wise Artificial Compressibility Method

The Artificial Compressibility Method (ACM) for the incompressible Navier-Stokes equations is (link-wise) reformulated (referred to as LW-ACM) by a finite set of discrete directions (links) on a regular Cartesian mesh, in analogy with the Lattice Boltzmann Method (LBM). The main advantage is the possibility of exploiting well established technologies originally developed for LBM and classical computational fluid dynamics, with special emphasis on finite differences (at least in the present paper), at the cost of minor changes. For instance, wall boundaries not aligned with the background Cartesian mesh can be taken into account by tracing the intersections of each link with the wall (analogously to LBM technology). LW-ACM requires no high-order moments beyond hydrodynamics (often referred to as ghost moments) and no kinetic expansion. Like finite difference schemes, only standard Taylor expansion is needed for analyzing consistency. Preliminary efforts towards optimal implementations have shown that LW-ACM is capable of similar computational speed as optimized (BGK-) LBM. In addition, the memory demand is significantly smaller than (BGK-) LBM. Importantly, with an efficient implementation, this algorithm may be one of the few which is compute-bound and not memory-bound. Two- and three-dimensional benchmarks are investigated, and an extensive comparative study between the present approach and state of the art methods from the literature is carried out. Numerical evidences suggest that LW-ACM represents an excellent alternative in terms of simplicity, stability and accuracy.Comment: 62 pages, 20 figure

Machine Learning Can Predict the Timing and Size of Analog Earthquakes

Despite the growing spatiotemporal density of geophysical observations at subduction zones, predicting the timing and size of future earthquakes remains a challenge. Here we simulate multiple seismic cycles in a laboratory‐scale subduction zone. The model creates both partial and full margin ruptures, simulating magnitude M_w 6.2–8.3 earthquakes with a coefficient of variation in recurrence intervals of 0.5, similar to real subduction zones. We show that the common procedure of estimating the next earthquake size from slip‐deficit is unreliable. On the contrary, machine learning predicts well the timing and size of laboratory earthquakes by reconstructing and properly interpreting the spatiotemporally complex loading history of the system. These results promise substantial progress in real earthquake forecasting, as they suggest that the complex motion recorded by geodesists at subduction zones might be diagnostic of earthquake imminence

Towards higher order lattice Boltzmann schemes

In this contribution we extend the Taylor expansion method proposed previously by one of us and establish equivalent partial differential equations of DDH lattice Boltzmann scheme at an arbitrary order of accuracy. We derive formally the associated dynamical equations for classical thermal and linear fluid models in one to three space dimensions. We use this approach to adjust relaxation parameters in order to enforce fourth order accuracy for thermal model and diffusive relaxation modes of the Stokes problem. We apply the resulting scheme for numerical computation of associated eigenmodes and compare our results with analytical references

A lattice Boltzmann study of non-hydrodynamic effects in shell models of turbulence

A lattice Boltzmann scheme simulating the dynamics of shell models of turbulence is developed. The influence of high order kinetic modes (ghosts) on the dissipative properties of turbulence dynamics is studied. It is analytically found that when ghost fields relax on the same time scale as the hydrodynamic ones, their major effect is a net enhancement of the fluid viscosity. The bare fluid viscosity is recovered by letting ghost fields evolve on a much longer time scale. Analytical results are borne out by high-resolution numerical simulations. These simulations indicate that the hydrodynamic manifold is very robust towards large fluctuations of non hydrodynamic fields.Comment: 17 pages, 3 figures, submitted to Physica

Seamount Subduction and Megathrust Seismicity: The Interplay Between Geometry and Friction

Subducting seamounts are recognized as one of the key features influencing megathrust earthquakes. However, whether they trigger or arrest ruptures remains debated. Here, we use analog models to study the influence of a single seamount on megathrust earthquakes, separating the effect of topography from that of friction. Four different model configurations have been developed (i.e., flat interface, high and low friction seamount, low friction patch). In our models, the seamount reduces recurrence time, interseismic coupling, and fault strength, suggesting that it acts as a barrier: 80% of the ruptures concentrate in flat regions that surround the seamount and only smaller magnitude earthquakes nucleate above it. The low-friction zone, which mimics the fluid accumulation or the establishment of fracture systems in natural cases, seems to be the most efficient in arresting rupture propagation in our experimental setting

Efficient and accurate simulations of deformable particles immersed in a fluid using a combined immersed boundary lattice Boltzmann finite element method

The deformation of an initially spherical capsule, freely suspended in simple shear flow, can be computed analytically in the limit of small deformations [D. Barthes-Biesel, J. M. Rallison, The Time-Dependent Deformation of a Capsule Freely Suspended in a Linear Shear Flow, J. Fluid Mech. 113 (1981) 251-267]. Those analytic approximations are used to study the influence of the mesh tessellation method, the spatial resolution, and the discrete delta function of the immersed boundary method on the numerical results obtained by a coupled immersed boundary lattice Boltzmann finite element method. For the description of the capsule membrane, a finite element method and the Skalak constitutive model [R. Skalak et al., Strain Energy Function of Red Blood Cell Membranes, Biophys. J. 13 (1973) 245-264] have been employed. Our primary goal is the investigation of the presented model for small resolutions to provide a sound basis for efficient but accurate simulations of multiple deformable particles immersed in a fluid. We come to the conclusion that details of the membrane mesh, as tessellation method and resolution, play only a minor role. The hydrodynamic resolution, i.e., the width of the discrete delta function, can significantly influence the accuracy of the simulations. The discretization of the delta function introduces an artificial length scale, which effectively changes the radius and the deformability of the capsule. We discuss possibilities of reducing the computing time of simulations of deformable objects immersed in a fluid while maintaining high accuracy.Comment: 23 pages, 14 figures, 3 table