Questions Regarding Ground Motion Field Residuals in Event-Based Calculator

[Luis Guillermo Álvarez]

Hello OQ team,

I hope you're all doing well!

I've recently been working with the event-based calculator and I'm trying to extract specific information from the simulated ground motion fields. In particular, I'm interested in the inter-event and intra-event residuals for each site-earthquake pair.

Upon inspecting the HDF5 file, I noticed that the inter-event residuals are available. However, I couldn’t find the intra-event residuals. Could you please confirm whether these are stored somewhere in the file?

Additionally, I’d like to associate the inter-event residuals with the corresponding events or ruptures. I assume this information might be under ['gmf_data']['rup_info'], but I’m not sure how to interpret the data in that section. Is there any documentation available that explains the structure of the HDF5 file? Or could you kindly clarify what the contents of ['gmf_data']['rup_info'] represent?

I had posted a similar question on the forum some time ago but haven’t received a response yet. I’d really appreciate any guidance you can provide.

Thanks in advance for your help!

Best regards,
Luis

[Marco Pagani]

Hi Luis,

The event-based calculator does not store the within-event variability.

It would be too expensive in terms of storage to also save this information (since it’s per-site).

Best wishes,

Marco

--
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 visit https://groups.google.com/d/msgid/openquake-users/59be0e60-3d5b-4e0c-8442-1aef66e6ed6fn%40googlegroups.com.

[Luis Guillermo Álvarez]

Hello Marco

Thank you for your quick response. I understand that the amount of information can become overwhelming. It might be useful to give users the option to activate this feature depending on the specific application.

Regarding the second part of my question, I was wondering if there is any documentation on how the information is stored. Since the inter_residual data does not include an ID column linking it directly to the events, I assume that its rows correspond to those in the eid key. That is, the first element in inter_residual would represent the residual associated with the earthquake at index 0 of eid—which, coincidentally, matches both the key and the event ID.

Best regards,

Luis

[Marco Pagani]

Hi Luis,

Sorry for the late reply. We are currently holding the semi-annual Governing Board meeting.

Regarding the first part, if you need this information for a few ruptures, the most efficient approach is to do it using a Python script.

Regarding the second part, please have a look here https://docs.openquake.org/oq-engine/manual/latest/user-guide/advanced/working-with-the-dstore.html

Otherwise, reach out to me.

Marco

To view this discussion visit https://groups.google.com/d/msgid/openquake-users/9c682686-f977-4a5b-a2fb-4374adcd4547n%40googlegroups.com.