6-C Wave parameter estimation: gulf of alaska earthquake

[sysdy]

Hello, why is my six-component seismic data parameter estimation inaccurate based on your tutorial "6-C Wave Parameter Estimation: Gulf

===============

This tutorial will teach you how to use TwistPy to process six-component seismic data recorded by collocated translational and rotational seismometers. The data used in this tutorial was recorded at the large ring laser observatory ROMY in Fürstenfeldbruck, Germany, after the January 2

In my parameter estimation, for the value of Vr, I set up the EstimatorConfiguration as follows:
est = EstimatorConfiguration(
wave_types=["R"],
method="DOT",
scaling_velocity=scaling_velocity,
use_ml_classification=False,
vr=[3000, 4000, 100],
xi=[-90, 90, 2],
phi=[150, 210, 1],
eigenvector=0,
)
No matter what range I set, it always displays the minimum value of this range, and all computed Vr values are the same. For example, vr=[3000, 4000, 100] ultimately results in a value of 3000. Do you know what might be causing this issue?

[solldavid]

Hi @sysdy

Thanks for reaching out! I am not sure I fully understand your question. I just ran the tutorial and it gives the exact same output as shown here.

Are you running this with your own data? If so, please note that setting phi=[150,210,1] limits the validity of the wave parameter estimation to Rayleigh waves traveling southward (with a propagation direction azimuth between 150 and 210 degrees). This is the predominant propagation direction expected for teleseismic waves traveling from Alaska to Germany, but it might not be a reasonable range for your own data. Please also make sure that your data complies with the coordinate system expected by TwistPy (North, East, Down).

[sysdy]

Thank you for your response. I am running the program using my own field-measured six-component data, but I am unable to estimate the parameters for the six components, especially the Rayleigh wave velocity Vr. No matter what range I set for Vr, the program always computes the lowest value within the range. For example, when I set vr=[3000, 4000, 100], the computed result is vr': array([[3000, 3000, 3000, 3000, 3000, 3000, 3000, 3000,..., 3000]]). Do you know how to solve this issue?