swiss-seismological-service / scdetect

A computationally efficient earthquake detection module for SeisComP
https://scdetect.readthedocs.io
GNU Affero General Public License v3.0
15 stars 6 forks source link

Network Magnitudes #103

Closed mmesim closed 2 years ago

mmesim commented 2 years ago

It would be useful to get network magnitudes for each detection. Then during the post-processing with scolv we can screen magnitudes easily.

mmesim commented 2 years ago

I'm writing this here in order not to forget about it.

I couldn't see any Magnitudes on scolv event list, perhaps because I was looking at unassociated origins.

I did the following ->

1) Edit the scolv module, tab VisibleMagnitudes, and add MLhc, MLh, MRelative, etc..

2) Import in the db a scml file that contains Events

3) I got this

Screenshot from 2022-05-05 16-45-38

So my guess is that if we assign Network Magnitudes, then perhaps we will see those on scolv. It might need a tiny configuration on scevent to make sure it understands MRelative.

The Magnitudes tab on scolv, where you can see Station Magnitudes, might be a bit more complicated.

Thoughts @luca-s @FMassin ?

luca-s commented 2 years ago

If you are talking only about unassociated origins your understanding is correct, you won't be able to see the magnitude. That is because each origin might have multiple magnitudes and scevent task (one of many) is to select the preferred magnitude among all the origins and among all the magnitudes. You see that selected/preferred magnitude in the event list tab.However if an origin is not associated, it means scevent was not run or it was not able to create an event for that origin (because it doesn't satisfy the parameters configured for scevent). Without an event, there is no magnitude visualized. However you should be able to view all origin available magnitudes in the Magnitudes tab

damb commented 2 years ago

However you should be able to view all origin available magnitudes in the Magnitudes tab.

@mmesim, can you confirm that the StationMagnitudes are displayed after adding e.g. MRelative to the list of visibleMagnitudes?

mmesim commented 2 years ago

I cannot see StationMagnitudes.

mmesim commented 2 years ago

I had a short conversation with @luca-s yesterday and it is more clear to me what happened.

In brief, and @luca-s please correct me if I'm wrong,

To conclude, it is crucial to get network magnitudes.

luca-s commented 2 years ago

@damb basically in the scdetect-cc origin we found the stationMagnitude object but the magnitude one is missing. Please have a look how a magnitude looks like:

image

damb commented 2 years ago

@mmesim, which approach do you suggest to compute the NetworkMagnitude? Do you require it to be configurable?

@luca-s thanks.

mmesim commented 2 years ago

I think the median of the station magnitudes will do the trick.

No