Modeling suggests that oblique extension facilitates rifting and continental break-up

[1] In many cases the initial stage of continental break-up was and is associated with oblique rifting. That includes break-up in the Southern and Equatorial Atlantic, separation from eastern and western Gondwana as well as many recent rift systems, like Gulf of California, Ethiopia Rift and Dead Sea fault. Using a simple analytic mechanical model and advanced numerical, thermomechanical modeling techniques we investigate the influence of oblique extension on the required tectonic force in a three-dimensional setting. While magmatic processes have been already suggested to affect rift evolution, we show that additional mechanisms emerge due to the three-dimensionality of an extensional syst…

A Framework to Assess the Information Dynamics of Source EEG Activity and Its Application to Epileptic Brain Networks

This study introduces a framework for the information-theoretic analysis of brain functional connectivity performed at the level of electroencephalogram (EEG) sources. The framework combines the use of common spatial patterns to select the EEG components which maximize the variance between two experimental conditions, simultaneous implementation of vector autoregressive modeling (VAR) with independent component analysis to describe the joint source dynamics and their projection to the scalp, and computation of information dynamics measures (information storage, information transfer, statistically significant network links) from the source VAR parameters. The proposed framework was tested on…

Constraining effective rheology through parallel joint geodynamic inversion

Abstract The dynamics of crust and lithosphere is to a large extent controlled by its effective viscosity. Unfortunately, extrapolation of laboratory experiments indicates that viscosity is likely to vary over many orders of magnitude. Additional methods are thus required to constrain the effective viscosity of the present-day lithosphere using more direct geophysical observations. Here we discuss a method, which couples 3D geodynamic models with observations (surface velocities and gravity anomalies) and with a Bayesian inversion scheme on massively parallel high performance computers. We illustrate that the basic principle of a joint geodynamic and gravity inversion works well with a simp…

Subduction initiation triggered the Caribbean large igneous province

AbstractSubduction provides the primary driving force for plate tectonics. However, the mechanisms leading to the formation of new subduction zones remain debated. An example is the Lesser Antilles Arc in the Atlantic. Previous initiation mechanisms have implied the transmission of subduction from the Pacific Ocean or the impact of a plume head. Here, we use geodynamic models to simulate the evolution of the Caribbean region during the Cretaceous, where the eastern Pacific subduction triggered the formation of a new subduction zone in the Atlantic. The simulations show how the collision of the old Caribbean plateau with the Central America margin lead to the formation of a new Atlantic subd…

Lower Crustal Rheology Controls the Development of Large Offset Strike‐Slip Faults During the Himalayan‐Tibetan Orogeny

International audience

Inferring rheology and geometry of subsurface structures by adjoint-based inversion of principal stress directions

SUMMARY Imaging subsurface structures, such as salt domes, magma reservoirs or subducting plates, is a major challenge in geophysics. Seismic imaging methods are, so far, the most precise methods to open a window into the Earth. However, the methods may not yield the exact depth or size of the imaged feature and may become distorted by phenomena such as seismic anisotropy, fluid flow, or compositional variations. A useful complementary method is therefore to simulate the mechanical behaviour of rocks on large timescales, and compare model predictions with observations. Recent studies have used the (non-linear) Stokes equations and geometries from seismic studies in combination with an adjoi…

Mutual Information Analysis of Brain-Heart Interactions in Epileptic Children

In this work we apply the network physiology paradigm to retrieve information from central and autonomic nervous systems before focal epileptic seizure, represented respectively by electroencephalogram (EEG) signals and R-R intervals (RRI), and investigate on the presence and strength of brain-heart interactions by computing mutual information (MI) measures. Statistical significance of MI values was tested through surrogate time series generated with the random shuffle approach. Our results suggest that the proposed method for aligning signals representing brain and heart activity measured with different sampling rates, is capable of revealing coupling between RRI representing heart system,…

Nonlinear brain-heart interactions in children with focal epilepsy assessed by mutual information of EEG and heart rate variability

Network physiology is a recent approach describing the human body as an integrated network composed of several organ systems which continuously interact to produce healthy and diseased states. In this work, we apply the network physiology paradigm to study dynamical interactions between EEG activity and heart rate variability in children suffering from focal epilepsy. We aim to study the characteristics of brainheart coupling between, before, and after seizures to better understand the physiological mechanisms underlying seizure onset in the pre-ictal phase and the recovery of normal autonomic function in the post-ictal phase. In perspective, linking the dynamic information of brain-heart c…

Deriving scaling laws in geodynamics using adjoint gradients

Abstract Whereas significant progress has been made in modelling of lithospheric and crustal scale processes in recent years, it often remains a challenge to understand which of the many model parameters is of key importance for a particular simulation. Determining this is usually done by manually changing the model input parameters and performing new simulations. For a few cases, such as for folding or Rayleigh-Taylor instabilities, one can use thick-plate stability analysis to derive scaling laws to obtain such insights. Yet, for more general cases, it is not straightforward to do this (apart from running many simulations). Here, we discuss a numerically cheaper approach to compute scalin…

Synergistic and Redundant Brain-Heart Information in Patients with Focal Epilepsy

In this work, partial information decomposition (PID) was applied to the time series of heart rate and EEG amplitude variability to investigate the dynamical interactions in brain-heart coupling before and after epileptic seizures. From ECG and EEG signals collected on 23 children suffering from focal epilepsy, the RR intervals and the EEG variance at ipsilateral and contralateral temporal electrodes were computed in four different time windows before and after the seizures. Static PID was used to obtain redundant, unique and synergistic components of the total information shared between the series of RR and EEG variance. Results highlight, in the progression from preictal to postictal stat…

Quantifying the thermo-mechanical impact of plume arrival on continental break-up

Abstract The arrival of a plume head at Earth's continental lithosphere is often considered to be an important factor for continental break-up. However, the impact of plume impingement on strength and duration of a rift remains unclear. In this study, we quantify the mechanical and thermal influence of a plume (i.e. lithosphere erosion) on continental break-up. To do that we apply the three-dimensional numerical code SLIM3D that features realistic elasto-visco-plastic rheology. We model the thermo-mechanical response of a segment of Earth's lithosphere that is affected both by extension as well as plume-related lithosphere erosion in order to evaluate the influence on the overall force budg…

Quantification and visualization of finite strain in 3D viscous numerical models of folding and overthrusting

Abstract Finite strain analysis and three-dimensional (3D) numerical modeling are important methods to understand the deformation history of rocks. Here, we analyze finite strain in 3D numerical simulations of power-law viscous folding and overthrusting. Simulations with different and laterally varying detachment geometries cause a lateral transition from folding (for thicker detachments) to overthrusting. We compute the 3D finite strain tensor, the principal strain values and their orientations. We compute the Nadai strain, e S , and the Lode’s ratio, ν , representing the strain symmetry (constriction or flattening). We apply Hsu diagrams to visualize strain distribution in e S - ν space, …

Equivalent continuum-based upscaling of flow in discrete fracture networks: The fracture-and-pipe model

Abstract. Predicting effective permeabilities of fractured rock masses is a key component of reservoir modelling. This is often realized with the discrete fracture network (DFN) method, where single-phase incompressible fluid flow is modelled in discrete representations of individual fractures in a network. Depending on the overall number of fractures, this can result in significant computational costs. Equivalent continuum models (ECM) provide an alternative approach by subdividing the fracture network into a grid of continuous medium cells, over which hydraulic properties are averaged for fluid flow simulations. While this has the advantage of lower computational costs and the possibility…

Comparison of continuous and discontinuous Galerkin approaches for variable-viscosity Stokes flow

We describe a Discontinuous Galerkin (DG) scheme for variable-viscosity Stokes flow which is a crucial aspect of many geophysical modelling applications and conduct numerical experiments with different elements comparing the DG approach to the standard Finite Element Method (FEM). We compare the divergence-conforming lowest-order Raviart-Thomas (RT0P0) and Brezzi-Douglas-Marini (BDM1P0) element in the DG scheme with the bilinear Q1P0 and biquadratic Q2P1 elements for velocity and their matching piecewise constant/linear elements for pressure in the standard continuous Galerkin (CG) scheme with respect to accuracy and memory usage in 2D benchmark setups. We find that for the chosen geodynami…

Information flow in EEG source networks in epileptic children with focal seizure activity

Scalp electroencephalographic (EEG) signals are influenced by several factors, including volume conduction and low spatial resolution, which can jeopardize the validity of brain connectivity analysis performed on the raw recordings. One possible solution is to identify, starting from scalp EEG signals, the underlying cortical source activations, and to apply connectivity metrics on the reconstructed source time series. In this work, the dynamics of information flow between cortical EEG signals obtained after source reconstruction were assessed in children suffering from focal epilepsy. In a group of 10 children with focal seizures, 5-second windows of the 19-channel EEG were obtained in the…

The effect of rheological approximations in 3-D numerical simulations of subduction and collision

Abstract Subduction and collision zones evolve differently from one another due to different rheological properties, different amounts of regional isostatic compensation, and the different mechanisms by which forces are applied to the convergent plates. The rheology of mantle and lithosphere is known to have the largest influence on the dynamics of subduction and continental collision. However, previous 3-D geodynamic models of subduction/collision processes have used various rheological approximations, making their results difficult to compare, since there is no clear understanding on the extent of these approximations on the dynamics. Here, we test the effect of rheological approximations…

Modeling evolution of the San Andreas Fault system in northern and central California

[1] We present a three-dimensional finite element thermomechanical model idealizing the complex deformation processes associated with evolution of the San Andreas Fault system (SAFS) in northern and central California over the past 20 Myr. More specifically, we investigate the mechanisms responsible for the eastward (landward) migrationof the San Andreas plate boundary over time, a process thathas largely determined the evolution and present structure of SAFS. Two possible mechanisms had been previously suggested. One mechanism suggests that the Pacific plate first cools and captures uprising mantle in the slab window, subsequently causing accretion of the continental crustal blocks. An alt…

Time, frequency and information domain analysis of short-term heart rate variability before and after focal and generalized seizures in epileptic children

OBJECTIVE In this work we explore the potential of combining standard time and frequency domain indexes with novel information measures, to characterize pre- and post-ictal heart rate variability (HRV) in epileptic children, with the aim of differentiating focal and generalized epilepsy regarding the autonomic control mechanisms. APPROACH We analyze short-term HRV in 37 children suffering from generalized or focal epilepsy, monitored 10 s, 300 s, 600 s and 1800 s both before and after seizure episodes. Nine indexes are computed in time (mean, standard deviation of normal-to-normal intervals, root mean square of the successive differences (RMSSD)), frequency (low-to-high frequency power rati…

Lithospheric stresses in Rayleigh–Bénard convection: effects of a free surface and a viscoelastic Maxwell rheology

Fe and Zn co-substituted beta-tricalcium phosphate (β-TCP): Synthesis, structural, magnetic, mechanical and biological properties

This work was supported by the European Social Fund under the No. 09.3.3- LMT-K-712 “Development of Competences of Scientists, other Researchers and Students through Practical Research Activities” measure. AK would like to express sincere gratitude for Fellowship administrated by The Japan Society for the Promotion of Science (JSPS). Fellow’s ID No.: L12546. Authors are grateful to R. Vargalis (Vilnius University) for taking SEM images. © 2020. This work is licensed under a CC BY-NC-ND license.

Entropy characteristics of heart rate wavelet multiscale components in epileptic children before and after seizures

In this work, we analyze the information content of the multiple time scale components of heart rate variability (HRV) in children with focal epilepsy. HRV components are extracted from 30 pediatric patients, monitored 10 min and 10 s before and after focal epileptic seizures, using wavelet multiscale decomposition (with 5, 15, 30, 60, 120, 180 s time scale), and then characterized computing Entropy (E), permutation entropy (PE), conditional entropy (CE) and information storage (IS). Moving from preictal to postictal windows, we find statistically significant differences in the CE and IS values of HRV components at short time scales, which reflect autonomic imbalance and appear as potential…

Pre- and post-ictal brain activity characterization using combined source decomposition and connectivity estimation in epileptic children

In this research, the study of functional connectivity between sources of electroencephalogram (EEG) activity assessed for different classes (well before seizure, preictal and post-ictal) was performed. EEG recordings were acquired from 12 subjects with focal epilepsy. Then, ten common spatial patterns (CSP) were obtained for EEG segments describing 95% of Riemannian distance between pairs of classes, followed by estimation of multivariate autoregressive (MVAR) models’ coefficients. The MVAR models were further used to extract coherence as a functional connectivity measures. Our results show that the coherence between CSP sources differs between baseline and pre-ictal segments: it has the l…