Skip to content

Added handling of missing SES in seismicity forecast

Cecilia Nievas requested to merge fix/skip_ses into main

The seismicity forecasts provided as input may not contain events for all realisations for which the forecast was run. This is the case when the software used to produce the seismicity forecast is set to output earthquakes above a certain threshold that results in some stochastic event sets being empty and, thus, non-existent in the forecast file. However, the user might want to take into account these "empty" or "missing" stochastic event sets when calculating the mean damage and losses associated with the whole forecast, as these stochastic event sets can be considered as not having produced any damage or losses.

The present MR addresses this issue by allowing the user to input two parameters in the config.yml file that control the IDs of the stochastic event sets to be taken into account:

  • continuous_ses_numbering: True/False (i.e. use IDs defined with consecutive numbering as specified in ses_range).
  • ses_range: ID_min, ID_max.

E.g., if continuous_ses_numbering=True and ses_range=1, 10, the software will consider 10 stochastic event sets with IDs [1, 2, 3, 4, 5, 6, 7, 8, 9, 10].

Closes https://git.gfz-potsdam.de/cnievas/real-time-loss-tools/-/issues/10.

Merge request reports