Point Source
SISMOWINE WP2_LOH1
SISMOWINE is intended as a long-term interactive web interface for verifying numerical modeling methods in seismology. Numerical-method developers and numerical modelers may compare their solutions with other solutions. SISMOWINE is a continuation of the original SPICE Code Validation interface established within the 6th Framework Programme project.
LOH1 is used as an example here to illustrate the implementation of source point for earthquake nucleation in SeisSol. The details of the LOH1 model can also be found at the point source example.
The model uses Right-handed Cartesian, x positive North, y positive East, z positive downward, all coordinates in meters. The source is buried at 2000 m in a half-space Earth (Figure [fig:loh1]. The top layer is 1000 m thick and the bottom layer is 33000 m. The material parameters are listed in Table [table:loh1].
Vp (m/s) |
Vs(m/s) |
density |
Qp |
Qs |
|
|---|---|---|---|---|---|
layer |
4000 |
2000 |
2600 |
Inf |
Inf |
half-space |
6000 |
3464 |
2700 |
Inf |
Inf |
Table: Material properties in LOH1.
Geometry of LOH1.
Geometry
The mesh is generated using Gmsh.
Geometry of LOH1 model (Gmsh). A 1 km layer of low velocity (Vp=4000 m/s, vs=2000 m/s) is at the top of high velocity (vp=6000 m/s, vs=3464 m/s).
Point source input
The point source needs to be turned on in parameter.par file.
&SourceType
Type = 50
FileName=’LOH1_source.dat’
/
The source input file can be generated using the script. Duration of the source is 4 seconds.
Results
Benchmark of x-component particle velocity at receiver point 1 (0.0, 693.0,0.1). Bule is 4-order SeisSol and orange is SISMOWINE result.