How to consider strike, and rake to green functions

PyLith
1

Hello everyone,
I calculated the green functions for the northern Chilean subduction zone in order to obtain the surface displacements for any slip model and compare it with GNSS observations.
I got the following, where red arrows are my predictions and blue arrows are the observations (black are the residuals, the magnitude of the residuals arrow is wrong so ignore them):

horizontal_carrasco
horizontal_carrasco1920×981 304 KB

As you can see, the arrows of the residuals tend to have an orientation with respect to some angle, which suggests to me that strike or rake (or both) might be playing a role. To be honest I dont really know if there is a way to set the strike or rake in the step.cfg. So, is there a way to set the strike and rake at the moment of compute the green functions? or in order to consider this I have to rotate the GFs later?

Another question I have is how pylith delivers the green functions. If I take Okada’s formulation as a reference, I have to calculate the contributions of both strike-slip, dip-slip, and tensile, then the G matrix will have the form:

G = G_ss + G_ds + G_ts

With this in mind, it is easy to assign strike, rake and dip to any sub-fault, but as far as I understood, what pylith gives me in this case is G right?. Is there a way to get G_ss, G_ds???

Here I upload you the files that im using for the simulation, Im using pylith 2.2.2.

step07b.cfg (4.5 KB)
pylithapp.cfg (6.7 KB)
greensfns.cfg (2.2 KB)

Here I show you the mesh I am using. I tried to make it as realistic as possible, considering topography, bathymetry, slab 2.0, a realistic moho
geometry and a detailed tomography for the elastic properties of the crust.

mesh_sized
mesh_sized1160×731 206 KB

2

In step07b.cfg you set the impulse_dof to [1], which means PyLith is computing the Green’s functions only for the reverse component of slip.

[pylithapp.problem.interfaces.slab]
# Impulses for reverse slip.
# Note that it is possible to apply both left-lateral and updip slip
# (impulse_dof = [0,1]), but we separate the impulses into two problems.
impulse_dof = [1]

In the subduction-3d example, we use step07a.cfg to compute the Green’s functions for the left-lateral component. See the subduction-3d example for how to make use of both the left-lateral and reverse components in an inversion.