No GMFs: GMFs below minimum_intensity threshold

132 views
Skip to first unread message

Jocelin AU TSANG

unread,
Aug 20, 2024, 3:03:13 AM8/20/24
to OpenQuake Users
Hi everyone,

I am new to OpenQuake as an honours student. I now encounter a runtime error when I tried to run scenario damage with condition under different GMPEs.
Runtime Error: No GMFs were generated, perhaps they were all below the minimum_intensity threshold.

I guess a possible reason is due to 'nan' nominal bias mean for GMPEs Bradley2013Volc, ChiouYounggs2014 and ChiouYoungsf2014Japan. When I ran scenario damage with ZhaoEtAl2016Asc, it ran through with a nominal bias -1.565.
[2024-08-13 16:04:56 #88 INFO] GSIM: [Bradley2013Volc], IMT: SA(0.6), Nominal bias mean: nan, Nominal bias stddev: 0.104
[2024-08-13 15:40:07 #83 INFO] GSIM: [ZhaoEtAl2016Asc], IMT: SA(0.6), Nominal bias mean: -1.565, Nominal bias stddev: 0.110

I am using OpenQuake version 3.20.1 and Python 3.11.7. I have attached my job_scenariodamge.ini as a photo.Screenshot 2024-08-15 at 23.46.18.png

Thanks!


Jocelin
August 2024

Anirudh Rao

unread,
Aug 27, 2024, 4:36:56 AM8/27/24
to OpenQuake Users
Hi Jocelin,

Probably the error message is misleading, but my guess is that there might be missing site parameters in your site model file. ZhaoEtAl2016Asc requires only vs30, but Bradley2013Volc, ChiouYoungs2014, and ChiouYoungs2014Japan require also vs30measured and z1pt0 in addition to vs30. The calculation might be failing silently due to the missing site paramaters.

Best regards,
Anirudh

Michele Simionato

unread,
Sep 1, 2024, 11:07:44 PM9/1/24
to OpenQuake Users
You need to use the command `oq zip job_scenariodamage.ini` and send to me the generated job.zip, otherwise it is impossible to debug the problem.

        Michele 
 
On Tuesday, August 20, 2024 at 9:03:13 AM UTC+2 jocelin....@gmail.com wrote:

Jocelin AU TSANG

unread,
Sep 12, 2024, 8:32:03 AM9/12/24
to OpenQuake Users
Hi Anirudh and Michele,

Thanks for your replies!

Anirudh, there are already measured, z1pt0 and z2pt5 in the input. Do you have any other suggestions on why there might be a problem?
Here is the first 15 lines of the [site para] file: 
<site lon="  106.564000" lat="   -5.441670" vs30="  356.31" vs30Type="measured" z1pt0="  150.66" z2pt5="    0.26"/>
    <site lon="  106.574000" lat="   -5.441670" vs30="  355.84" vs30Type="measured" z1pt0="  151.11" z2pt5="    0.26"/>
    <site lon="  106.534000" lat="   -5.451360" vs30="  355.10" vs30Type="measured" z1pt0="  151.82" z2pt5="    0.27"/>
    <site lon="  106.544000" lat="   -5.451360" vs30="  355.07" vs30Type="measured" z1pt0="  151.84" z2pt5="    0.27"/>
    <site lon="  106.554000" lat="   -5.451360" vs30="  355.38" vs30Type="measured" z1pt0="  151.55" z2pt5="    0.26"/>
    <site lon="  106.564000" lat="   -5.451360" vs30="  355.36" vs30Type="measured" z1pt0="  151.57" z2pt5="    0.26"/>
    <site lon="  106.574000" lat="   -5.451360" vs30="  355.35" vs30Type="measured" z1pt0="  151.58" z2pt5="    0.26"/>
    <site lon="  106.584000" lat="   -5.451360" vs30="  354.29" vs30Type="measured" z1pt0="  152.58" z2pt5="    0.27"/>
    <site lon="  106.524000" lat="   -5.461040" vs30="  353.72" vs30Type="measured" z1pt0="  153.13" z2pt5="    0.27"/>
    <site lon="  106.534000" lat="   -5.461040" vs30="  354.32" vs30Type="measured" z1pt0="  152.56" z2pt5="    0.27"/>
    <site lon="  106.544000" lat="   -5.461040" vs30="  354.26" vs30Type="measured" z1pt0="  152.62" z2pt5="    0.27"/>
    <site lon="  106.554000" lat="   -5.461040" vs30="  354.42" vs30Type="measured" z1pt0="  152.47" z2pt5="    0.27"/>
    <site lon="  106.564000" lat="   -5.461040" vs30="  354.40" vs30Type="measured" z1pt0="  152.48" z2pt5="    0.27"/>
    <site lon="  106.574000" lat="   -5.461040" vs30="  354.43" vs30Type="measured" z1pt0="  152.46" z2pt5="    0.27"/>
    <site lon="  106.584000" lat="   -5.461040" vs30="  353.60" vs30Type="measured" z1pt0="  153.25" z2pt5="    0.27"/>
...

Michele, please find the attached zip file. 

Thanks again!


Best,
Jocelin
Jocelin_job.zip

Anirudh Rao

unread,
Sep 12, 2024, 9:12:34 AM9/12/24
to openqua...@googlegroups.com
Hi Jocelin,

Converting your site model into a csv file, the attached job is working for me.

Anirudh

--
You received this message because you are subscribed to the Google Groups "OpenQuake Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to openquake-use...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/openquake-users/2a624cc6-f294-43aa-b9e7-c8aeec9b9284n%40googlegroups.com.
job_scenariodamage.ini
cipta_Vs30.csv

Jocelin AU TSANG

unread,
Sep 12, 2024, 10:13:44 AM9/12/24
to OpenQuake Users
Hi Anirudh,

Thanks for your prompt reply.

I have tried running oq with the files you provided me. It calculated nominal bias mean equals nan, same as my previous run. It then stopped at this line, and my computer crashed and needed to be restarted.

[2024-09-12 23:44:36 #140 INFO] GSIM: [Bradley2013Volc], IMT: SA(0.3), Nominal bias mean: nan, Nominal bias stddev: 0.129



Best,
Jocelin

Siti Raudha Kaisuku

unread,
Sep 22, 2024, 11:07:03 AM9/22/24
to openqua...@googlegroups.com

Hello everyone,

I have never used OpenQuake before for seismic hazard analysis using the probabilistic method. I would like to ask if anyone has a script for processing background sources, faults, and subduction zones so that I can process the data in OpenQuake.

Thank you all in advance!

Reply all
Reply to author
Forward
0 new messages