Velocity models

This is a list of the velocity models currently implemented in TieBeNN. As further local, dedicated velocity models are published/discovered, they can be integrated in TieBeNN as specified further below.

List of available velocity models

Number assigned

Description

0

IASP91 (Kennett et al. 1995)

1

AK135 (Montagner & Kennett 1996)

2

Layer-over-halfspace BGR velocity model (Schlittenhardt 1999)

3

WET, velocity model for Germany and adjacent regions (reference missing)

4

WB2012, Vogtland/West Bohemia velocity model (1D, P+S), averaged from the 3D models of Růžek & Horálek (2013)

5

DEU, 3-layer model for Germany with a Moho depth of 28.5 km. A vp/vs ratio of 1.68 was used to generate S-wave velocities (reference missing)

6

Crust1.0 model based on epicenter location. Layers from Leeman & Ammon

7

Crust1.0 + AK135: Crustal layers from Crust1.0, deeper structure from AK135

8

Insheim model (1D, P+S) from Küperkoch et al. (2018)

9

Regional model for Landau-Insheim (by J. Borns)

10

Landau industry model (1D, P+S) for northern Upper Rhine Graben

11

1D model for Central Alps (Diehl et al. 2021)

12

(Deactivated) 3D model for Central Alps (Diehl et al. 2021) – implementation under review

13

(Deactivated) 3D WEG model (shallow, high-res industry model around Rotenburg/Söhlingen/Soltau)

14

AlpsLocPS, 1D model for the Greater Alpine Region (Brazsus et al. 2024)

15

BENS, 1D model for Northern Rhine Area (Reamer & Hinzen 2004)

16

ASZmod1, 1D model for Albstadt Shear Zone (Mader et al. 2021)

17

(Deactivated) 3D P+S velocity model for the Upper Rhine Graben (Lengliné et al. 2023)

18

PO1 (1D P+S) for Novy-Kostel swarm region, constrained (Málek et al. 2023)

19

PO1 + WB2012 (from 11 km depth)

20

KIT6 (1D P+S) for the Eastern Eifel Volcanic Field (Ritter et al. 2024)

When using --vel_mode manual, the parameter --velmod must be set to one of the model numbers above.

When the Crust1.0 model is selected, TieBeNN automatically extracts a 1D velocity profile at the epicenter coordinates.

Structure of NonLinLoc 1D velocity model files

Each crustal velocity layer corresponds to one line in the text file using this format:

LAYER   DepthKm   VelP   GradVelP   VelS   GradVelS   Dens   GradDens

  • VelP, VelS in km/s

  • GradVelP, GradVelS in km/s/km (linear gradients)

  • Dens in kg/m³, GradDens in kg/m³/km

Note

The density is a carry-over field originally used in a waveform modeling code. It is not used in NonLinLoc programs, so any convenient value can be used (A. Lomax, pers. comm.).

How to add a new 1D velocity model

  1. Navigate to the folder: <repo>/tiebenn/data/velocity_models/ Create a new file named v_<number> (e.g. v21), with one LAYER line per layer.

  2. (Optional) List your model in the function docstring of velmods() in tools/velocity_models.py.

  3. Update/reinstall TieBeNN to integrate your new model (in the repository main directory):

    pip install .

  4. (Optional / recommended) Spread the love! Let the maintainer know via an issue or merge request so the model can be added to the official repo 🌍