SCEC Workshop on Advancing Simulations of Sequences of Earthquakes and Aseismic Slip
Conveners: Brittany Erickson (U of Oregon), Junle Jiang (Oklahoma), Valere Lambert (UCSC)
Dates: November 2, 2021 (online)
SCEC Award: 21139
SUMMARY: Our group is working to advance computational methods for simulating Sequences of Earthquakes and Aseismic Slip (SEAS) by conducting a suite of code verification exercises. In SEAS models, the goal is to capture the interplay of aseismic fault slip—that ultimately lead to earthquake nucleation—and dynamic earthquake ruptures themselves, in an effort to understand which physical factors control the full range of observables such as aseismic deformation, earthquake initiation locations, ground shaking during dynamic rupture, recurrence times and magnitudes of earthquakes.
A total of 60 people participated in the online SCEC SEAS Workshop on November 2, 2021, including scientists from the U.S.A., Sweden, France, Switzerland, China, Germany, England, Japan, Canada, Brazil, and New Zealand. Almost half of our workshop participants were graduate students or postdocs.
The goal of this workshop was to share some recent progress in SEAS modeling that considers fluid effects and increased frictional and geometric complexity and to introduce our next benchmark problems BP6 and BP7, which focus on fluid and 3D effects, respectively. Results were summarized from the most recent benchmarks BP3-QD, BP1-FD, and BP5 which focused on dipping fault geometries, full dynamic and 3D effects, respectively. The group also discussed research targets and plans for SEAS in the near and longer term for the group, following scientific talks presented by SEAS modelers.
SESSION 1. Brittany Erickson kicked off the workshop with a brief history of SEAS activities since its initiation in 2018 and an overview of the workshop schedule and goals. All participants were then invited to introduce themselves and briefly describe their interest in SEAS and the workshop. Brittany then provided an update on the current status of model results for the two 2D benchmarks: BP1-FD (the first fully dynamic problem) and BP3-QD (the first plane strain problem with a dipping fault). Quantitative agreements have been made for BP1-FD (Figure 1) for several modeling groups, while other participating modeling groups continue to update/upload new results. For BP3-QD, the on-fault time series matched well across codes, and discrepancies in off-fault surface stations have been addressed by resolving an ambiguity in the benchmark description. A draft of the paper presenting these results will be circulated among co-authors and modelers by the end of 2021.
Next, Junle Jiang summarized highlights and lessons from the code comparison exercises for 3D SEAS benchmarks. A total of 10 modeling groups together explored how various computational and physical factors affect simulation outputs for all phases of multiple seismic cycles. The excellent agreement among best-resolved simulations lend confidence to participating numerical codes, offering best-practice examples to improve numerical simulations (Figure 2). These comparisons also reveal different sensitivities of model observables, which are important for integrating models with geophysical observations. The manuscript for 3D benchmarks is in review as of December 2021.
SESSION 2-3. Five science talks presented by Stacy Larochelle, Yuyun Yang, Sylvain Barbot, Pierre Romanet, and Duo Li (4 of which are early career scientists) focused on using SEAS modeling to understand the effects of fluids, friction, complex fault geometries and full dynamics on fault slip.
SESSION 4. During the afternoon session, Valère Lambert gave an overview of two new sets of benchmark problems targeted for 2022 (Figure 3). The first, BP6-QD-A/S/C, is a 2D quasi-dynamic problem of a 1-D fault with either constant or velocity-strengthening friction exhibiting an aseismic slip transient arising due to changes in effective normal stress from fluid injection and along-fault pore fluid diffusion. The second, BP7-QD/FD-A/S, is a 3D problem focusing on the interplay between aseismic slip and earthquake nucleation on a circular velocity-weakening asperity. BP6 and 7 have been modified based on feedback from group discussions, which includes using a smooth nucleation procedure for the initial rupture in BP7 and suggesting a finite domain size for volume-based methods.
SESSION 5. To conclude the meeting, the workshop conveners led the planning discussions for SEAS activities in 2022, including possible future benchmark problems, how to best maximize participation in SEAS benchmarks, and the importance of high-performance computing in the advancement of code. The group also discussed possible collaborations between SEAS and other groups (e.g. RSQSim) in future activities.
The final discussion highlighted three general problems to address for future benchmarks: 1) a 2D antiplane problem like BP1 but including bulk viscoelastic deformation, 2) a 2D dipping fault problem like BP3 with depth-variable normal stress, and 3) a 2D problem including a low-velocity damage zone. Looking forward, the group discussed options for revamp of the online platform to improve specific tools and functionality.
Presentation videos and presentations may be viewed by clicking the links below. PLEASE NOTE: Files are the author’s property. They may contain unpublished or preliminary information and should only be used while viewing the talk. Only the presentations for which SCEC has received permission to post publicly are included below.
TUESDAY, NOVEMBER 2, 2021
All times below are Pacific Daylight Time (PDT or UTC-7).
08:00 - 08:50 | Session 1: Workshop Goals, Introductions, Results to Date Moderator: Tran Huynh |
|
08:00 - 08:10 | Workshop Goals and Science Targets (PDF, 550KB) | Brittany Erickson |
08:10 - 08:30 | Group Introductions (PDF, 2MB) | All |
08:30 - 08:50 | Review of SEAS Group Benchmark Results (PDF, 2.8MB) | Brittany Erickson and Junle Jiang |
08:50 -09:00 | Break | |
09:00 - 10:00 | Session 2: Fluid Effects and Friction Moderator: Valère Lambert |
|
09:00 - 09:20 | Fluid-induced slip and earthquake nucleation on a rate-and-state fault (VIDEO) | Stacy Larochelle |
09:20 - 09:40 | Fault slip coupled to fluid migration with porosity and permeability evolution (VIDEO) | Yuyun Yang |
09:40 - 10:00 | Fault dynamics in non-isothermal conditions (VIDEO) | Sylvain Barbot |
10:00 - 10:10 | Break | |
10:10 - 11:00 | Session 3: Complex Geometries and Physics Moderator: Junle Jiang |
|
10:10 - 10:30 | Fully-dynamic earthquake cycle on a non planar fault using the spectral boundary element method (VIDEO) | Pierre Romanet |
10:30 - 10:50 | Spatiotemporal complexity of slow slip events on 3D non-planar subduction faults | Duo Li |
10:50 - 11:00 | Group Discussion | |
11:00 - 12:00 | Break | |
12:00 - 13:00 | Session 4: New Benchmarks Moderators: Brittany Erickson and Junle Jiang |
|
12:00 - 12:20 | New Benchmarks BP6 and 7 (PDF, 2.0MB) | Valère Lambert |
12:20 - 13:00 | Group Discussion | |
13:00 - 14:00 | Session 5: Future Directions for SEAS Moderators: Brittany Erickson and Valère Lambert |
|
13:00 - 13:30 | Group Planning for 2022 (PDF, 445KB) | |
13:30 - 14:00 | Group Discussion: Recommendations for Future Research Beyond 2022 | |
14:00 | Adjourn |
It is SCEC policy to foster harassment-free environments wherever our science is conducted. By accepting an invitation to participate in a SCEC-supported event, by email or online registration, participants agree to abide by the SCEC Activities Code of Conduct.
PARTICIPANTS
Kali Allison (UC Davis)
Martin Almquist (Uppsala Univ)
Pablo Ampuero (Uni Côte d'Azur)
Michael Barall (USGS)
Sylvain Barbot (USC)
Yehuda Ben-Zion (USC)
Segun Bodunde (U Oklahoma)
Alexandre Chen (U Oregon)
Luca Dal Zilio (ETHZ)
Luis Dalguer (3Q-Lab)
Shuai Deng (Sichuan Univ)
Ben Duan (TAMU)
Eric Dunham (Stanford)
Brittany Erickson (U Oregon)
Alice Gabriel (Munich)
Christine Goulet (USC)
Ruth Harris (USGS)
Tobias Harvey (U Oregon)
Jessica Hawthorne (Oxford)
Yihe Huang (Michigan)
Tran Huynh (USC)
Junle Jiang (U Oklahoma)
Yoshi Kaneko (GNS)
Taeho Kim (Caltech)
Jeremy Kozdon (NPS)
Christos Kyriakopoulos (U Memphis)
Valère Lambert (UCSC)
Stacy Larochelle (Caltech)
Duo Li (LMU Munich)
Dunyu Liu (UT Austin)
Yajing Liu (McGill)
Amy Lu (McGill)
Betsy Madden (Brasilia)
Phil Maechling (USC)
Rishav Mallick (EOS/Caltech)
Brendan Meade (Harvard)
Enrico Milanese (MIT)
Shiying Nie (USC)
So Ozawa (U Tokyo)
Edric Pauk (USC)
Andrea Perez (Victoria)
Fred Pollitz (USGS)
Casper Pranger (LMU Munich)
Sohom Ray (Dalhousie)
Pierre Romanet (NIED)
Nico Schliwa (LMU Munich)
Paul Segall (Stanford)
Md Shumon Mia (UIUC)
Yudong Sun (MIT)
Yuval Tal (Caltech)
Prithvi Thakur (Michigan)
Ben Thompson
Yuan Tian (Stanford)
Carsten Uphoff (LMU Munich)
Yongfei Wang (USC)
Meng Wei (URI)
Joseph Wick (UCSC)
Yuyun Yang (Stanford)
Wenqiang Zhang (McGill)