Virtues of modeling many faults: New method illuminates shape of Alaskan quake

In a new study published in Scientific Reports, a research team led by the University of Tsukuba developed a new method to model the details of complex earthquake rupture processes affecting systems of multiple faults. They then applied this method to the magnitude 7.9 earthquake that struck the Gulf of Alaska near Kodiak Island on January 23, 2018.

As study co-author Professor Yuji Yagi explains, "Our method uses a flexible finite-fault inversion framework with improved smoothness constraints. This approach allows us to analyze seismic P waves and estimate the focal mechanisms and rupture evolution of geometrically complex earthquakes involving rupture of multiple fault segments."

Based on the distribution of aftershocks within one week of the main shock of the Gulf of Alaska earthquake, this method was applied to represent slip along a horizontal plane at a depth of 33.6 km.

The main rupture stage of the earthquake, which lasted for 27 seconds, affected fault segments oriented both north-south and east-west.

"Our results confirm previous reports that this earthquake ruptured a conjugate fault system in a multi-shock sequence," says study first author Shinji Yamashita. "Our model further suggests that this rupture tended to occur along weak zones in the sea floor: fracture zones that extend east-west, as well as plate-bending faults that run parallel to north-south-oriented magnetic lineaments."

These features caused discontinuities in the fault geometry that led to irregular rupture behavior. "Our findings show that irregular rupture stagnation 20 kilometers north of the earthquake's epicenter may have been promoted by a fault step across the seafloor fracture zone," explains co-author Assistant Professor Ryo Okuwaki, "They also indicate a causal link between rupture evolution and pre-existing bathymetric features in the Gulf of Alaska."

This method represents a promising step forward in modeling earthquake rupture processes in complex fault systems based only on seismic body waves, which may improve modeling of seismic wave propagation and mapping of complex fault networks in tectonically active areas.