Earthquake Modeling using InSAR

Dear all,

I got the InSAR results in both ascending and descending orbits for my study area and an earthquake happened in that area. I would like to model it using Inversion and Okada. Could anybody please suggest me a good paper with a clear description of this model? I got some papers and read them but unfortunately, they are not very clear to implement and program by myself.

Best regards,

1 Like

May I ask what data you used to get these interferograms?
They look quite nice!

1 Like

I used S1-A images to get these interferograms in my study area using SNAP software


Sorry for not answering your question but I would like to comment that this is probably the best looking interferogram I have ever seen.


@qglaude Really?! :blush: thank you guys for your nice speaking :slightly_smiling_face:


Yes, very beautiful interferograms. I assume this is the M6 earthquake in Iran last year.

This GitHub site has code implementing the Okada dislocation models:


Note that the GitHub code of Ben Thompson has only the MATLAB and Python wrappers for the Okada programs. You have to also download the Fortran code from the link at the GitHub page. The wrappers make it easier to run the programs.

1 Like


Thank you so much for your great tips. I will try to understand and apply in my case study. Can I ask you to suggest some papers or references which help me to understand earthquake modeling especially using InSAR results? Because I would like to improve my knowledge in this field and also implement in my research.

Best regards,

This paper might be helpful:
Funning, G. J., B. Parsons, T. J. Wright, J. A. Jackson, and E. J. Fielding (2005), Surface displacements and source parameters of the 2003 Bam (Iran) earthquake from Envisat advanced synthetic aperture radar imagery, Journal of Geophysical Research B: Solid Earth , 110 (9), B09406, doi:10.1029/2004JB003338.

This paper is also a good place to start and includes a program implementation:
Feigl, K. L., and E. Dupré (1999), RNGCHN: a program to calculate displacement components from dislocations in an elastic half-space with applications for modeling geodetic measurements of crustal deformation, Computers and Geosciences , 25 (6), 695-704, doi:10.1016/S0098-3004(99)00003-5.

@EJFielding Thank you in advance for your suggestion. I will download them and read as soon as possible. I know that you are very busy :see_no_evil: but I really appreciate if you let me ask here my questions if I have some misunderstanding after reading them.

Sorry for my request

Dear all,

I calculated the amplitude image of two S1 images, which are used to generate that interferogram. I added them to related products as you can see in the following. Whenever I do the coregistration using S-1 TOPS Coregistration option, It ignore the amplitude as a output while I also need to coregister and match the amplitude of two images.

Could anybody please guide me about which option I can use to do it using SNAP?

I’m not sure why you need the amplitude but in my case, the coregistration produces the intensity of each image.


The intensity is nothing but the amplitude at the power of 2. Just take the square root of it to get the amplitude. If you plan to use it eslewhere of SNAP, don’t forget to convert the band (so that it is not a band created on-the-fly) and save the product.

@qglaude Thank you so much. I need to have an amplitude to do the following analysis based on the suggested references by Prof. Fielding.


Yes, I know. I generated the amplitude using Intensity and band math in SNAP, after applying TOPSAR-Split and Apply-Orbit-File for both two images. I need to coregister or match two amplitude products in any way to do that analysis.

isn’t the Offset-Tracking Tool exactly doing this? Offset tracking

@ABraun I don’t know. I am reading this reference now.

Michel, R., J.-P. Avouac, and J. Taboury (1999), Measuring ground displacements from SAR amplitude images: Application to the Landers earthquake, Geophys. Res. Lett., 26, 875 – 878.

Whenever I get the result and knowledge about this topic, I will share here with you. I found this paper and that technique interesting :slightly_smiling_face:

I think it is quite similar to the method of the Offset Trackign Tool. There is also a tutorial on it:
I would give it a try.

The tool takes GRD products.

@ABraun Yes, I agree with you based on this sentence from that paper:


Thank you so much for the tutorial. Does it only use the GRD products? Is it not possible with SLC images?

As far as I tested it the tool takes only GRD. But the coregistration is able to make good use of the resolution.

@ABraun Thank you so much. I will try all ways and share the results here. As the last question, this tool is only for S1 satellite or we can use for others like Envisat?