How HPC is Shaking Up Modeling of Mysterious Earthquakes

By Oliver Peckham

July 6, 2021

Earth modeling is extraordinarily complex – and when the earth quakes, it’s no different. The mechanisms of many earthquakes consistently flummox seismologists, making it difficult to understand why those earthquakes occurred, let alone perform the most important and elusive function in seismology: earthquake prediction. In a talk during ISC21, Alice-Agnes Gabriel, a professor of seismology at the Ludwig Maximilian University of Munich, explained how HPC-enabled earthquake codes have worked to close these gaps and take advantage of quickly accelerating supercomputing capabilities.

Dr. Alice-Agnes Gabriel

“Computational seismology has been a pioneering field for high-performance computing,” Gabriel said, “and [it] has itself been pioneered by HPC in order to image Earth’s interior to understand the dynamics of the metal, for example, and also to track down energy resources.” Computational seismology, she explained, just made sense for HPC: seismology is data-rich, and many of its dynamics could often be treated as linear systems.

On the other hand: “If you’re thinking about computational earthquake seismology, we’re facing a very different picture. Computational earthquake seismology is trying to understand earthquakes in a physics-based manner – we’re not imaging earthquakes, but we’re trying to model earthquakes based on physical first-order principles. We’re solving for spontaneous dynamic earthquake rupture as non-linear interaction of frictional failure across prescribed zones of weakness that are geological faults, that can take complicated shapes and interact with each other … and these frictional failure processes are linked to seismic wave propagation in a non-linear manner.”

So, in layperson’s terms: hunting for oil deposits, for instance, is easier than figuring out how many earthquakes occurred. And, indeed, there are many baffling earthquakes: earthquakes that do more damage than seismologists expect, earthquakes that occur along seemingly low-tension faults and even earthquakes that jump.

To model earthquakes, Gabriel explained, seismologists combine knowledge from a wide range of fields: seismology, geodesy, geology, tectonophysics, hydrology, numerical computing, data science, machine learning, applied mathematics, continuum mechanics, tribology, rock mechanics, materials science, and engineering. They amalgamate these processes in bootstrapped models, developing solvers that are able to read data from fault stresses and output rupture dynamics, ground deformation and synthetic observables (like simulated seismograms).

These models are robust, but three major challenges remain: first, earthquake source processes are very ill-constrained and highly non-linear; second, researchers are still working to identify which physical processes are dominant and relevant in real earthquakes and how to computationally budget for those processes’ inclusion in models; and third, teams struggle to assimilate all available knowledge in a manner both hardware- and software-suitable for modern HPC systems.

Progress on these challenges has, thankfully, been rapidly accelerating. Gabriel cited some of her team’s prior research from 2014 – itself a Gordon Bell Prize finalist – that modeled a large 1992 earthquake in California that confused researchers by seemingly jumping from fault to fault without a continuous connection. Ever since, she said, researchers had been modeling earthquakes at a petaflop scale, linking spontaneous dynamic rupture simulations with fault geometry and a range of other physical fields.

In the intervening seven years, of course, the scale is increasing: researchers are working to port those models to ever-larger spatiotemporal ranges in order to capture megathrust earthquakes (which represent all recorded earthquakes over magnitude 9.0) and tsunamis. “The extreme scale of megathrust earthquakes and tsunamis,” Gabriel said, “also requires extreme-scale runs.”

The investments required for such runs are rapidly decreasing. A “hero run” simulation of the 2004 Sumatra earthquake and tsunami, Gabriel said, took 14 hours on LRZ’s SuperMUC Phase 2 (2.8 Linpack petaflops) in 2017; now, typical runs take under six hours on just 16 nodes of LRZ’s SuperMUC-NG (19.5 Linpack petaflops).

Gabriel’s team uses a seismic solver called – appropriately – SeisSol. Using SeisSol, the researchers can solve seismic wave equations on unstructured tetrahedral meshes, combining complex geometries with heterogeneous inputs and high resolutions. Initially, SeisSol was simply a Fortran-based solver with MPI parallelization. Over time, it evolved: in 2014, for simulations of the 1992 California quake, the team implemented hybrid MPI/OpenMP parallelization, parallel I/O and a multiphysics offload function for many-core architectures, allowing the simulations to reach 96 billion degrees of freedom at 200,000 time steps.

For the Sumatra quake, SeisSol was further optimized through a bevy of features, including asynchronous I/O, overlapping computation and communication, and a feature called cluster-based local time-stepping that allowed each element to have its own time step. These optimizations allowed SeisSol to expand to 111 billion degrees of freedom at 3.3 million time steps, offering a 14-fold speedup and allowing the team to complete Sumatra simulations in hours rather than days. Essentially, Gabriel said, these speedups were the difference between a feasible simulation and an impossible simulation.

Image courtesy of Alice-Agnes Gabriel.

“We’ve been using this approach to shed light on puzzling earthquake observations,” Gabriel said. Lately, for instance, the researchers have been studying the 2016 Kaikōura earthquake, a magnitude 7.8 event that struck New Zealand to devastating – and befuddling – effect, perhaps most notably when the earthquake leapt from one fault to another fault 15km away, three times the previously thought maximum distance for a “jumping” earthquake. Other questions stumped researchers, too: who did the earthquake move so slowly? Why did the Hope Fault, which was thought to pose the greatest seismic hazard in the region, barely react?

A representation of the Kaikōura earthquake simulation. Image courtesy of Alice-Agnes Gabriel.

SeisSol, it turns out, was able to answer many of these questions that would have dead-ended researchers a decade ago. Using the solver, the team was able to show, for instance, that another fault – the Point Kean fault – acted as a crucial link between the two endpoints of the jump, and that unfavorable dynamic stresses on the Hope Fault led to its lack of activity.

The researchers are excited to develop even more functionality, exploring the crevices of earthquake modeling. Gabriel elaborated on how her team is modeling manmade earthquakes, such as the 2017 magnitude 5.5 earthquake in South Korea, and modeling the 2018 Sulawesi earthquake that devastated Indonesia – which was complicated by the unavailability of close seismometers recording the event, leading Gabriel and her colleagues to turn to sparse data from satellites and social media.

Improvements in the code, of course, are complemented by improvements in supercomputing. Gabriel’s team has made extensive use of top-tier supercomputers, including modern leaders like SuperMUC-NG, the 5.5-Linpack petaflops Shaheen II at KAUST, the 5.4-Linpack petaflops Mahti at CSC and the 23.5-Linpack petaflops Frontera at TACC.

All of this will be amplified by the forthcoming GPU port of SeisSol – and, of course, the imminent shift to the exascale era will open up a wide range of further possibilities. Gabriel mused on what they could learn from running the models a thousand times faster: higher resolutions, uncertainty quantification, rapid earthquake response simulations, more multi-physics components – a paradigm shift toward more holistic earthquake system modeling. “The interplay of advances in high-performance computing and dense observations,” Gabriel said, “will allow us to go beyond scenario-based analysis.”

Subscribe to HPCwire's Weekly Update!

Be the most informed person in the room! Stay ahead of the tech trends with industy updates delivered to you every week!

White House Scientific Integrity Report Addresses AI and ML Ethics

January 26, 2022

Earlier this month, the White House Office of Science and Technology Policy (OSTP) Scientific Integrity Task Force released a report titled “Protecting the Integrity of Government Science.” While broad-based and over Read more…

IBM Quantum Debuts Classical Entanglement Forging to Expand Simulation Capabilities

January 26, 2022

IBM last week reported a new technique – entanglement forging – that uses both quantum and classical computing resources to double the size of select simulation problems that can be solved on current quantum computer Read more…

Lenovo Launches Its TruScale HPC as a Service Offering

January 26, 2022

Lenovo today announced TruScale High Performance Computing as a Service (HPCaaS), which it says will offer a “cloud-like experience” to HPC organizations of all sizes. The new HPC-as-a-Service is part of the TruScale Read more…

Ceremorphic Touts Its HPC/AI Silicon Technology as It Exits Stealth

January 25, 2022

In a market still filling with fledging silicon chips, Ceremorphic, Inc. has exited stealth and is telling the world about what it calls its patented new ThreadArch multi-thread processor technology that is intended to h Read more…

Quantum Watch: Neutral Atoms Draw Growing Attention as Promising Qubit Technology

January 25, 2022

Currently, there are many qubit technologies vying for sway in quantum computing. So far, superconducting (IBM, Google) and trapped ion (IonQ, Quantinuum) have dominated the conversation. Microsoft’s proposed topologic Read more…

AWS Solution Channel

Register for the AWS “Speeds n’ Feeds” event on Feb. 9th

Since the debut of the first ‘Beowulf’ cluster in 1994, HPC has been a race between technologists squeezing as much performance as possible from hardware, and scale economics driving mass-production prices to affordable levels. Read more…

Meta’s Massive New AI Supercomputer Will Be ‘World’s Fastest’

January 24, 2022

Fresh off its rebrand last October, Meta (née Facebook) is putting muscle behind its vision of a metaversal future with a massive new AI supercomputer called the AI Research SuperCluster (RSC). Meta says that RSC will b Read more…

Lenovo Launches Its TruScale HPC as a Service Offering

January 26, 2022

Lenovo today announced TruScale High Performance Computing as a Service (HPCaaS), which it says will offer a “cloud-like experience” to HPC organizations of Read more…

Ceremorphic Touts Its HPC/AI Silicon Technology as It Exits Stealth

January 25, 2022

In a market still filling with fledging silicon chips, Ceremorphic, Inc. has exited stealth and is telling the world about what it calls its patented new Thread Read more…

Quantum Watch: Neutral Atoms Draw Growing Attention as Promising Qubit Technology

January 25, 2022

Currently, there are many qubit technologies vying for sway in quantum computing. So far, superconducting (IBM, Google) and trapped ion (IonQ, Quantinuum) have Read more…

Meta’s Massive New AI Supercomputer Will Be ‘World’s Fastest’

January 24, 2022

Fresh off its rebrand last October, Meta (née Facebook) is putting muscle behind its vision of a metaversal future with a massive new AI supercomputer called t Read more…

IBM Watson Health Finally Sold by IBM After 11 Months of Rumors

January 21, 2022

IBM has sold its underachieving IBM Watson Health unit for an undisclosed price tag to a global investment firm after almost a year’s worth of rumors that sai Read more…

Supercomputer Analysis Shows the Atmospheric Reach of the Tonga Eruption

January 21, 2022

On Saturday, an enormous eruption on the volcanic islands of Hunga Tonga and Hunga Haʻapai shook the Pacific Ocean. The explosion, which could be heard six tho Read more…

NSB Issues US State of Science and Engineering 2022 Report

January 20, 2022

This week the National Science Board released its biannual U.S. State of Science and Engineering 2022 report, as required by the NSF Act. Broadly, the report presents a near-term view of S&E based mostly on 2019 data. To a large extent, this year’s edition echoes trends from the last few reports. The U.S. is still a world leader in R&D spending and S&E education... Read more…

Multiverse Targets ‘Quantum Computing for the Masses’

January 19, 2022

The race to deliver quantum computing solutions that shield users from the underlying complexity of quantum computing is heating up quickly. One example is Multiverse Computing, a European company, which today launched the second financial services product in its Singularity product group. The new offering, Fair Price, “delivers a higher accuracy in fair price calculations for financial... Read more…

IonQ Is First Quantum Startup to Go Public; Will It be First to Deliver Profits?

November 3, 2021

On October 1 of this year, IonQ became the first pure-play quantum computing start-up to go public. At this writing, the stock (NYSE: IONQ) was around $15 and its market capitalization was roughly $2.89 billion. Co-founder and chief scientist Chris Monroe says it was fun to have a few of the company’s roughly 100 employees travel to New York to ring the opening bell of the New York Stock... Read more…

US Closes in on Exascale: Frontier Installation Is Underway

September 29, 2021

At the Advanced Scientific Computing Advisory Committee (ASCAC) meeting, held by Zoom this week (Sept. 29-30), it was revealed that the Frontier supercomputer is currently being installed at Oak Ridge National Laboratory in Oak Ridge, Tenn. The staff at the Oak Ridge Leadership... Read more…

AMD Launches Milan-X CPU with 3D V-Cache and Multichip Instinct MI200 GPU

November 8, 2021

At a virtual event this morning, AMD CEO Lisa Su unveiled the company’s latest and much-anticipated server products: the new Milan-X CPU, which leverages AMD’s new 3D V-Cache technology; and its new Instinct MI200 GPU, which provides up to 220 compute units across two Infinity Fabric-connected dies, delivering an astounding 47.9 peak double-precision teraflops. “We're in a high-performance computing megacycle, driven by the growing need to deploy additional compute performance... Read more…

Intel Reorgs HPC Group, Creates Two ‘Super Compute’ Groups

October 15, 2021

Following on changes made in June that moved Intel’s HPC unit out of the Data Platform Group and into the newly created Accelerated Computing Systems and Graphics (AXG) business unit, led by Raja Koduri, Intel is making further updates to the HPC group and announcing... Read more…

Nvidia Buys HPC Cluster Management Company Bright Computing

January 10, 2022

Graphics chip powerhouse Nvidia today announced that it has acquired HPC cluster management company Bright Computing for an undisclosed sum. Unlike Nvidia’s bid to purchase semiconductor IP company Arm, which has been stymied by regulatory challenges, the Bright deal is a straightforward acquisition that aims to expand... Read more…

D-Wave Embraces Gate-Based Quantum Computing; Charts Path Forward

October 21, 2021

Earlier this month D-Wave Systems, the quantum computing pioneer that has long championed quantum annealing-based quantum computing (and sometimes taken heat fo Read more…

Killer Instinct: AMD’s Multi-Chip MI200 GPU Readies for a Major Global Debut

October 21, 2021

AMD’s next-generation supercomputer GPU is on its way – and by all appearances, it’s about to make a name for itself. The AMD Radeon Instinct MI200 GPU (a successor to the MI100) will, over the next year, begin to power three massive systems on three continents: the United States’ exascale Frontier system; the European Union’s pre-exascale LUMI system; and Australia’s petascale Setonix system. Read more…

Three Chinese Exascale Systems Detailed at SC21: Two Operational and One Delayed

November 24, 2021

Details about two previously rumored Chinese exascale systems came to light during last week’s SC21 proceedings. Asked about these systems during the Top500 media briefing on Monday, Nov. 15, list author and co-founder Jack Dongarra indicated he was aware of some very impressive results, but withheld comment when asked directly if he had... Read more…

Leading Solution Providers

Contributors

Lessons from LLVM: An SC21 Fireside Chat with Chris Lattner

December 27, 2021

Today, the LLVM compiler infrastructure world is essentially inescapable in HPC. But back in the 2000 timeframe, LLVM (low level virtual machine) was just getting its start as a new way of thinking about how to overcome shortcomings in the Java Virtual Machine. At the time, Chris Lattner was a graduate student of... Read more…

2021 Gordon Bell Prize Goes to Exascale-Powered Quantum Supremacy Challenge

November 18, 2021

Today at the hybrid virtual/in-person SC21 conference, the organizers announced the winners of the 2021 ACM Gordon Bell Prize: a team of Chinese researchers leveraging the new exascale Sunway system to simulate quantum circuits. The Gordon Bell Prize, which comes with an award of $10,000 courtesy of HPC pioneer Gordon Bell, is awarded annually... Read more…

Meta’s Massive New AI Supercomputer Will Be ‘World’s Fastest’

January 24, 2022

Fresh off its rebrand last October, Meta (née Facebook) is putting muscle behind its vision of a metaversal future with a massive new AI supercomputer called t Read more…

Nvidia Defends Arm Acquisition Deal: a ‘Once-in-a-Generation Opportunity’

January 13, 2022

GPU-maker Nvidia is continuing to try to keep its proposed acquisition of British chip IP vendor Arm Ltd. alive, despite continuing concerns from several governments around the world. In its latest action, Nvidia filed a 29-page response to the U.K. government to point out a list of potential benefits of the proposed $40 billion deal. Read more…

Julia Update: Adoption Keeps Climbing; Is It a Python Challenger?

January 13, 2021

The rapid adoption of Julia, the open source, high level programing language with roots at MIT, shows no sign of slowing according to data from Julialang.org. I Read more…

Top500: No Exascale, Fugaku Still Reigns, Polaris Debuts at #12

November 15, 2021

No exascale for you* -- at least, not within the High-Performance Linpack (HPL) territory of the latest Top500 list, issued today from the 33rd annual Supercomputing Conference (SC21), held in-person in St. Louis, Mo., and virtually, from Nov. 14–19. "We were hoping to have the first exascale system on this list but that didn’t happen," said Top500 co-author... Read more…

TACC Unveils Lonestar6 Supercomputer

November 1, 2021

The Texas Advanced Computing Center (TACC) is unveiling its latest supercomputer: Lonestar6, a three peak petaflops Dell system aimed at supporting researchers Read more…

10nm, 7nm, 5nm…. Should the Chip Nanometer Metric Be Replaced?

June 1, 2020

The biggest cool factor in server chips is the nanometer. AMD beating Intel to a CPU built on a 7nm process node* – with 5nm and 3nm on the way – has been i Read more…

  • arrow
  • Click Here for More Headlines
  • arrow
HPCwire