Use SeisCloud with SeisClient

SeisClient Tutorial

Request Greens function and the synthetic waveform from point sources, multiple-point sources, and finite faults for various applications such as earthquake source characterization, induced seismicity monitoring, and Earthquake Early Warning (EEW).

SeisCloud provides the remote data server that hosts the 3D SGT database.

SeisClient can request precise 3D Greens’ functions and theoretical seismograms from SeisCloud.

1
2
3
4
# load packages
from seisclient.seisclient import SeisClient
from obspy.core.util.attribdict import AttribDict
import numpy as np

Check models and stations in SeisCloud

1
2
3
4
5
6
7
8
9
# request model parameters
req = SeisClient()

# specify the model name
model = 'SOCAL3D'

# request the model parameters
info_dict, _ = req.request_model_info(model)
print(info_dict)

Request synthetic waveform from point source

  • Set station, source, saving_dir
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
# station (eg: station SLA in Network CI)
station = AttribDict({
    'latitude': 35.890,
    'longitude': -117.283,
    'network': 'CI',
    'station': 'SLA',
    'location': '',
    'id': 'SLA'})

# source
origin = AttribDict({'latitude': 35.601,
                     'longitude': -117.597,
                     'depth_in_m': 2810.0,
                     'id': 'evt1'})

# the folder storing the request data. 
saving_dir = './data/'
  • Request synthetic waveform with CMTSOLUTION
1
2
3
4
5
6
7
8
# CMTSOLUTION
# set moment tensor eg: CMT in dyn.cm -> n.m
mt_cmt = np.array([-2.900000e+25, -4.120000e+26, 4.410000e+26, 8.580000e+25, 4.470000e+25, 1.010000e+26])/1E7

# request synthetic waveform 
syn_cmt = req.get_synthetic_mtsolution(model=model, station=station, origin=origin, save_dir=saving_dir, mt_RTP=mt_cmt, b_RTZ=True)
print(syn_cmt)
syn_cmt.plot(outfile="syn_cmt.png", show=False)
  • Request synthetic waveform with FORCESOLUTION
1
2
3
4
5
6
7
8
# FORCESOLUTION
# set force solution in ENZ.
force_enz = [1e15, 2e15, 3e15]

# request synthetic waveform
syn_f = req.get_synthetic_forcesolution(model=model, station=station, origin=origin, save_dir=saving_dir, force_enz=force_enz, b_RTZ=True)
print(syn_f)
syn_f.plot(outfile="syn_f.png", show=False)

Request 3D Greens function

  • 3D Greens function
1
2
3
# 3D Greens function
greens_3d = req.get_greens_3DMT(model=model, station=station, origin=origin, save_dir=saving_dir)
print(greens_3d)

Request synthetic waveform from multiple point sources and finite fault

Examples

0%