Virtual Quake is a boundary element code designed to investigate long term fault system behavior and interactions between faults through stress transfer.
We currently have the ability to prescribe aftershocks for each simulated earthquake, by sampling the ETAS aftershock distribution functions.
Turning on this behavior is as simple as setting the simulation parameter sim.bass.max_generations = 1
HOWEVER: Currently there is a bug of some kind in the aftershock code. The simulation seems to have essentially an infinite number of earthquakes, with the simulation amassing hundreds of earthquakes in a tenth of a year interval.
This project first requires VALIDATION of the aftershock code, then a simple implementation of a simulation with and without aftershocks and quantifying the statistics.
Aftershock code is very short and simple, and is located in src/simulation/RunEvent.cpp in the processAftershock() function. There is also some code in src/simulation/BASSAftershocks.h and src/simulation/BASSAftershocks.cpp
We currently have the ability to prescribe aftershocks for each simulated earthquake, by sampling the ETAS aftershock distribution functions. Turning on this behavior is as simple as setting the simulation parameter
sim.bass.max_generations = 1HOWEVER: Currently there is a bug of some kind in the aftershock code. The simulation seems to have essentially an infinite number of earthquakes, with the simulation amassing hundreds of earthquakes in a tenth of a year interval.
This project first requires VALIDATION of the aftershock code, then a simple implementation of a simulation with and without aftershocks and quantifying the statistics.
Aftershock code is very short and simple, and is located in src/simulation/RunEvent.cpp in the
processAftershock()function. There is also some code in src/simulation/BASSAftershocks.h and src/simulation/BASSAftershocks.cpp