pointsource.rst

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

The comparison with solution is shown in Figure [fig:compare_loh1].

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.