Next-day gridded forecasts

Grid-based forecast is one of the two forecasts types used for CSEP experiments and supported by PyCSEP. For a specific forecasting period, a grid-based forecast is specified by providing the expected number of events (or rate) for each square bin in a pre-defined space-magnitude grid. In this way, grid-based forecasts can be formatted as tables with each row representing a space-magnitude bin and each column reporting informations about the bin location (e.g. space coordinates, magnitude, depth). Grid-forecasts for multiple forecasting periods are then formatted as multiple tables, one for each period.

Here, we provide access to 27 models that, within CSEP, produced daily grid-based forecasts (between 00:00:01 and 23:59:59) of events with magnitude above 3.95 (or 4.95 depending on the model) and depth below 30 km, covering the period from August 1, 2007 to August 30, 2018 for a total of over 50,000 forecasts. These seismicity models includes several ensemble models, non-parametric models, and different flavors of the well-established Epidemic-Type Aftershock Sequence (ETAS), and Short-Term Earthquake Probability (STEP) models.

All the grid-based forecasts provided here are formatted accordingly to the PyCSEP default format . They are a collection of .dat files (one for each day covered by the model) where row represents space-magnitude bins and column contains all the information needed to identify the bin and the rate. More specifically the columns (in order) are lower longitude extreme, upper longitude extreme, lower latitude extreme, upper latitude extreme, lower depth extreme, upper depth extreme, lower magnitude extreme, upper magnitude extreme, rate, and a flag. A more detailed explanation can be found on the PyCSEP website. All the forecasts are based on the same space-magnitude grid which consists of 0.1° x 0.1° (lat x lon), 0.1 Mw units, and one bin for the depth going from 0 to 30km.

The forecasts are provided as a .hdf5 file storing a system of nested folders. The system of nested folders is structered in the following way: one folder per year, inside each year folder we have one folder per month, and inside each month folder we have the .dat file representing the forecasts for the days of the respective month and year. We provided a tutorial (link) on how to use the .hdf5 file in combination with PyCSEP to evaluate the forecasts against observed data. This comprehends testing the consistency of the forecasts with the observations and comparing forecasts from different models.

List of Models

Model nameStarting dateEnd dateMissing daysMin magnitudeLinks
ETAS01 August 200730 August 201803.951. Description
2. Download


pyCSEP v0.6.3 is now available on PyPI and conda-forge. Visit the GitHub page for more information.

Recent BlogS