E - A Realistic Time-independent Experiment

TL; DR

In a terminal, navigate to floatcsep/examples/case_e and type:

$ floatcsep run config.yml

After the calculation is complete, the results will be summarized in results/report.md.

Artifacts 

This example shows how to run a realistic testing experiment in Italy (based on https://doi.org/10.4401/ag-4844), summarizing the concepts of previous examples. The example has only a subset of the original models and evaluations. The input structure of the experiment is:

case_e
    └──  models
        ├── gulia-wiemer.ALM.italy.10yr.2010-01-01.xml
        ├── meletti.MPS04.italy.10yr.2010-01-01.xml
        └── zechar.TripleS-CPTI.italy.10yr.2010-01-01.xml
    ├── config.yml
    ├── models.yml
    └── tests.yml

Configuration 

Time 

The time configuration is manifested in the time-config inset.

time_config:
  start_date: 2010-1-1T00:00:00
  end_date: 2022-1-1T00:00:00
  intervals: 4
  growth: cumulative

Region 

The testing region is the official Italy CSEP Region obtained from csep.core.regions.italy_csep_region.
region_config:
  region: italy_csep_region
  mag_min: 5.0
  mag_max: 9.0
  mag_bin: 0.1
  depth_min: -2
  depth_max: 30

Catalog 

The catalog is obtained from an authoritative source, namely the Bollettino Sismico Italiano (http://terremoti.ingv.it/en/bsi ), using the function csep.query_bsi()
catalog: query_bsi

Models 

The models are set in models.yml. For simplicity, there are only three of the nineteen models originally submitted to the Italy Experiment.
- ALM:
    path: models/gulia-wiemer.ALM.italy.10yr.2010-01-01.xml
    forecast_unit: 10
    use_db: True
- MPS04:
    path: models/meletti.MPS04.italy.10yr.2010-01-01.xml
    forecast_unit: 10
    use_db: True
- TripleS-CPTI:
    path: models/zechar.TripleS-CPTI.italy.10yr.2010-01-01.xml
    forecast_unit: 10
    use_db: True
The .xml format is automatically detected and parsed by floatcsep readers.

Note

The forecasts are defined in [Earthquakes / 10-years], specified with the forecast_unit option.

Note

The use_db flag allows floatcsep to transform the forecasts into a database (HDF5), which speeds up the calculations.

Post-Process 

Additional options for post-processing can set using the postproc_config option.

postproc_config:
  plot_forecasts:
    region_border: True
    basemap: stock_img
    clabel_fontsize: 14
    cticks_fontsize: 12
    alpha_exp: 0.8
    catalog:
      legend: True
      legend_loc: 1
      legend_fontsize: 14
      legend_titlesize: 16
      markersize: 7
      markercolor: blue

See plot_spatial_dataset() for forecast plot options and plot_catalog() for the catalog placed on top.

Running the experiment 

The experiment can be run by simply navigating to the examples/case_a folder in the terminal and typing.
floatcsep run config.yml
This will automatically set all the calculation paths (testing catalogs, evaluation results, figures) and will create a summarized report in results/report.md.