Slip rate during the Landers Earthquake from Hernandez et al.1999
Initiation of rupture from Campillo and Ionescu 1997
My research focuses on seismic waves, their sources, and how they propagate through complex media. Understanding these phenomena allows to create models of the Earth and its dynamics that provide new insights into various aspects of geophysics and improve our understanding of seismic hazards. My work began with the problem of modelling of the seismic radiation from dynamic cracks. These developments were applied to the imaging of the rupture process of large earthquakes. A new inversion approach in the frequency domain was proposed. It was successfully applied to describe the growth and propagation of rupture at depth from local seismic and geodetic data (GPS, InSAR).
An earthquake corresponding to a drastic change of magnitude of deformation velocity, the question of the nature of the instability is a key problem that was addressed through a novel mathematical approach based on a spectral method. It allows to describe precisely the phase of initiation and to characterize the transition stable/unstable in complex geometries. We introduced a concept of effective friction of a heterogeneous fault based on spectral arguments. My work on seismic source is still in progress. I worked on the inversion of slip history during silent earthquakes from permanent geodetic data and we analyzed the scaling properties of slip during earthquakes.
Scattering by a topography from Sanchez-Sesma and Campillo 1991
Understanding wave propagation is the central issue for the interpretation of seismic records. I studied extensively the regional propagation of short period waves using numerical simulation and data analysis. I showed the frequency dependence of the classical quality factor of rocks and the importance of scattering. I proposed maps of quality factors and discovered zones of extinction of guided waves across the mountain ranges of Alps and Pyrenees. A method of simulation of wave propagation based on boundary integral equations was developed to investigate such complex geological structures. This method was also applied to problems of seismic hazard assessment (effect of basin or topography) and of discrimination of earthquakes and explosions.
Modified from Margerin et al. 2000
The quantitative analysis of seismograms convinced me to turn my interests to the problem of scattering of seismic waves in the Earth. This represents a significant change of perspective in my career since I had to investigate some problems that are traditionally more relevant to physics than to seismology. We develop a model of coda decay based on the radiative transfer theory in which the traditional seismological parameter Qcoda is interpreted as a time of residence of diffuse waves in the Earth crust. In modern physics, the concept of ‘mesoscopic’ regime opened a new field of investigation of a world of phenomena between microscopic and macroscopic behaviours. We found that multiple scattering of seismic waves is relevant of this approach. We demonstrated the existence of the phenomena of equipartion of energies in the seismic coda, and we showed that energy ratios are in agreement with theory. We showed also the first demonstration of weak localisation of seismic waves, a typical mesoscopic effect in which diffuse average energy is exhibiting a behaviour associated with reciprocal properties of the wave field and interferences.
Comparison between the nine cross correlations and the nine components of the theoretical Green tensor from Campillo and Paul 2003
An interesting consequence of this preservation of phase effect in supposedly random fields is the possibility of passive reconstruction of the response of the Earth between two distant points from the correlation of the field recorded at these points. We demonstrated the pertinence of this concept with actual coda waves. Then we applied the same approach to records obtained in absence of earthquake, that is with ambient noise records. In this case, and in the period band we considered, the so-called seismic noise consists of the permanent agitation produced by the interaction between ocean waves and the solid Earth.
Using long continuous noise records, we have reconstructed the Earth response between stations and use these virtual seismograms to extract surface waves and to produce 3D tomographic images of the Earth interior with an unprecedented resolution.
The development of this technique was a breakthrough, which rapidly led to numerous applications ranging from near-surface exploration to deep Earth structures such as the mantle-core boundary.
Change seismic velocity associated with earthquake activity, from Brenguier et al. 2008
We studied carefully the precision and stability of this novel type of data. We found that the measurements can be remarkably stable. We therefore proposed to use it for monitoring temporal changes at depth with repeated noise-based measurements. We showed a new type of signal preceding the eruptions of a volcano. The long-term work on random fields and mesoscopy led to new approaches to investigate the changing Earth. My present priority is to develop new seismological tools to monitor slight changes in fault zones that could be associated with the evolution of the tectonic stress. This includes new imaging techniques based on the radiative transfer theory with, for the first time, the coupling of body and surface waves.
Space-time monitoring of groundwater fluctuations with passive seismic interferometry, from Mao et al. 2022
We also apply matricial methods developed in medical imaging and optics to build high resolution images of the scattering strength to complement the traditional seismic representations for volcanoe and fault zones (left: magmatic structure beneath a volcano from Blondel et al. 2018
Indication that apparently continuous slow slip consists of a cluster of slow transient, from Frank et al. 2018
In the last years, we studied the large slow slip events recently discovered in Mexico by using jointly GPS and seismic data. We developed techniques that allow for imaging the slip evolution during the events and for detecting low magnitude events. We have seen the triggering of slow slip by a distant mega-earthquake and the triggering of a large local event by a slow slip event.
One of our findings is the intermittent nature of the large slow slip events that can be decomposed in a cascade of short transient slips by using simultaneously geodetic and seismological data. This discovery challenges the present-day widely accepted models of friction.
Response of the seismic signal to surface temperature from Steinmann et al. 2022
The signature of the continuous dynamics of a volcano revealed in the 3C record of a single station, from Steinmann et al., 2024
Recently we started a project for applying deep learning to data analysis. A chair in the newly installed Interdisciplinary Institute for Artificial Intelligence of UGA (MIAI) supports this program. The goal of this project, based on multiscale representations of the signals, is to integrate geodetic and seismological data in the same classification. Scatnet is an analysis tool developed as part of the ERC F-IMAGE and used for a number of applications.
We recently obtained positive results: predicting GPS time series from seismic data. We are now reaching another of our goals: extract from raw seismograms a feature space in which it is possible to blindly separate the different sources of waves, but also to disentangle source and medium change effects.
This opens up the possibility of extracting data directly from temporal variations in rock properties (due to meteorological forcing or deep deformation) and from changes in weak signal sources at depth.