New imaging technique uses Earth's warped surface to reveal rocky interior

The new technique, described by researchers as "deformation imaging," provides results comparable to seismic imaging but offers direct information about the rigidity of the planet's crust and mantle. This property is essential for understanding how earthquakes and other large-scale geological processes work, said Simone Puel, who developed the method for a research project at the University of Texas Institute for Geophysics while in graduate school at the UT Jackson School of Geosciences.

"Material properties like rigidity are critical to understand the different processes that occur in a subduction zone or in earthquake science in general," Puel said.

"When combined with other techniques like seismic, electromagnetic or gravity, it should be possible to actually produce a much more comprehensive mechanical model of an earthquake in a way that has never been done before."

Puel, who is now a postdoctoral scholar at the California Institute of Technology, published the theory behind his method in Geophysical Journal International earlier this year. A recent study published in June in Science Advances shows it in action. It used GPS data recorded during Japan's 2011 Tohoku earthquake to image the subsurface down to about 100 kilometers underground.

The image revealed the tectonic plates and volcanic system beneath the Japanese portion of the Pacific Ring of Fire, including an area of low rigidity that's thought to be a deep magma reservoir feeding the system—the first time such a reservoir has been detected using only surface information.

The method relies on the fact that Earth's crust is a hodge-podge of rocky material with differing elastic properties. Some parts are more pliant, and other parts are more rigid. This causes the crust to contract and expand unevenly. During an earthquake, for example, the Earth vibrates in a way that reflects what it's made of, leaving the surface deformed in telltale ways.

To turn this uneven deformation into an image of the subsurface, the researchers constructed a computer model that treats the Earth as if it is a simplified elastic material, while allowing its elastic strength to vary in three dimensions.

The model then computed the subsurface rigidity based on how much the GPS sensors had moved in relation to one another during the earthquake. The result is a 3D picture of the Earth's interior based on changes on the surface.

An advantage of the new method is that it can use measurements made by satellites. These include NASA's upcoming NISAR spacecraft, a joint mission with the Indian Space Research Organization that will map the entire globe in very high resolution every 12 days.

Using the new technique, NISAR could offer important insights into some of the world's most geologically hazardous regions, said study co-author Thorsten Becker, a professor at the Jackson School. By continuously mapping the Earth's surface, the satellite will allow scientists to track structural changes in earthquake faults as they progress through their earthquake cycle.

Co-author Omar Ghattas, a professor at the UT Walker Department of Mechanical Engineering and UT Oden Institute for Computational Engineering and Sciences, said that the new method could be an important step to building digital twins of the Earth. These complex computer models perpetually improve themselves by identifying where to make new observations, then assimilating the new data.

More information: S Puel et al, An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data, Geophysical Journal International (2023). DOI: 10.1093/gji/ggad442

Simone Puel et al, Volcanic arc rigidity variations illuminated by coseismic deformation of the 2011 Tohoku-oki M9, Science Advances (2024). DOI: 10.1126/sciadv.adl4264

Journal information: Geophysical Journal International, Science Advances

Source: University of Texas at Austin