EGU22-10523, updated on 28 Mar 2022
EGU General Assembly 2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.

A Discontinuous-Galerkin approach to model non-classical nonlinearity observed from lab to global scales

Zihua Niu1, Alice-Agnes Gabriel1, Dave May1,2, Christoph Sens-Schönfelder3, and Heiner Igel1
Zihua Niu et al.
  • 1University of Munich, Institute of Geophysics, Department of Earth and Environmental Sciences, Munich, Germany (
  • 2Scripps Institution of Oceanography, UC San Diego, La Jolla, USA
  • 3Helmholtz Centre Potsdam, German Research Centre for Geosciences GFZ, Telegrafenberg, Section 2.4, D-14473 Potsdam, Germany

Under dynamic perturbations, it has been observed that materials like sedimentary rocks show complex mechanical behaviors. They include the simultaneous dependence of the elastic moduli and attenuation on strain at the same time scale of the perturbations, as well as the conditioning and recovery of the elastic moduli that may happen at time scales that are much larger. The latter cases were recently referred to as non-classical nonlinearity. Aside from laboratory experiments, comparable observations of the non-classical nonlinearity have been made in the field over the past two decades with the development of long-term continuous monitoring of the velocity field inside the Earth using methods such as ambient noise interferometry.


A variety of mathematical models that can potentially quantify the non-classical nonlinearity have already been proposed, e.g., the Damage–Breakage Rheology Model, the Internal Variable Model and the Godunov–Peshkov–Romenski model. However, implementing them in numerical schemes suitable to reproduce nonlinear effects in wave propagation on the local, regional, or global scale is challenging. This can be of interest for constraining a more realistic dynamic rheology for the Earth with the field observations.


In this work, wave propagation in different non-classical nonlinear models is implemented in FEniCS using the discontinuous Galerkin (DG) method in 1D. Behaviors of the different models are systematically studied and quantitatively compared against measurements. This work lays the foundation for an extension to the simulation of 2D/3D wave propagation in the Earth on the large-scale DG simulation frameworks, e.g., SeisSol and ExaHyPE.

How to cite: Niu, Z., Gabriel, A.-A., May, D., Sens-Schönfelder, C., and Igel, H.: A Discontinuous-Galerkin approach to model non-classical nonlinearity observed from lab to global scales, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-10523,, 2022.

Comments on the display material

to access the discussion