Skip to content

Addresses Large M Rupture Problem - Changes to Stochastic Rupture Generator

Graeme Weatherill requested to merge bug/large_m_fix into main

For some events in the OELFs, large magnitudes (M > 8.0) are occaisionally being produced. For these cases the stochastic rupture generator was generating surfaces that could not be validated by OpenQuake. This MR refactors the inputs and handling of the regional properties to constrain large ruptures including:

  1. Adds optional user control of the rupture regional properties (magnitude scaling relation, aspect ratio limits, largest area magnitude, upper seismogenic depth and lower seismogenic depth) to make them tectonic region dependent (as defined in the source model file). If not specified then this returns to a default set of properties, some of which will be over-ridden in the source model.
  2. Introduces the variable "largest area magnitude (area_mmax)", which sets a cap in the rupture generator such that the rupture area will not exceed that corresponding to the specified magnitude for any earthquake exceeding that magnitude
  3. Changes the handling of the random number generator such that the same random samples can be drawn according to a specific seed that can be set in the configuration file

Merge request reports