Open calum-chamberlain opened 1 month ago
@calum-chamberlain Thanks very much for the query - helped us identify a bug in converting SCXML to QUAKEML in SeisComp. We will notify about this to Gempa.
In SCXML files, the uncertainty ellipsoid units are in meters, but when they are converted to QUAKEML using the stylesheet in SeisComp, they are multiplied by 1000, and hence the very large numbers you see.
We suspect this has been happening for events located with LOCSAT from ~2018-10-10 for using both NLL (https://service.geonet.org.nz/fdsnws/event/1/query?eventid=2018p760856) and LOCSAT (https://service.geonet.org.nz/fdsnws/event/1/query?eventid=2018p763141).
Whilst we figure this out with GEMPA, there are two options for you:
Hope this helps! Please let us know if you have any more information regarding this issue.
Cheers
Sent a ticket to Gempa to discuss and address this question
FYI @calum-chamberlain @pasansherath
a XSLT fix is proposed (conversionsheet) as in https://forum.seiscomp.de/uploads/short-url/vKuK8XEiraqGbQz2tkhznzmXv2C.tar
$SEISCOMP_ROOT/share/xml/
in seiscomp releases NB a change in the stylesheet ellipsoid units occurred in 08.2017 it looks like we revert it back. now
FYI @pasansherath
To be tested still 2024 Aug.
@pasansherath @salichon @calum-chamberlain
Further to this issue raised above, there is also another instance whereby the latitude and longitude uncertainties are inconsistent with the units given for the horizontal uncertainties.
I believe the origin uncertainties are too small (given in km) rather than in metres. This is different to the issue raised by Calum above where the uncertainties have been incorrectly multipled by 1000, rather some examples exist of them being divided by 1000. The issue seems to be widespread, though I haven't looked into the time period it affects.
An event example is: 2018p132507
from obspy.clients.fdsn import Client
client = Client("GEONET")
event = client.get_events(eventid="2018p132507")[0]
print(f"Max horizontal uncertainty: {event.preferred_origin().origin_uncertainty.max_horizontal_uncertainty}")
print(event.preferred_origin().origin_uncertainty.confidence_ellipsoid)
print(f"Latitude uncertainty: {event.origins[-1].latitude_errors.uncertainty}")
print(f"Longitude uncertainty: {event.origins[-1].longitude_errors.uncertainty}")
returns:
Max horizontal uncertainty: 7.434938502
ConfidenceEllipsoid
semi_major_axis_length: 11.48772057
semi_minor_axis_length: 4.294251743
semi_intermediate_axis_length: 7.526263248
major_axis_plunge: 142.2719735
major_axis_azimuth: -20.77884346
major_axis_rotation: 91.00429164
Latitude uncertainty: 2.724444892
Longitude uncertainty: 4.673446065
Here you can see that the horizontal uncertainties indicate (if given in the standard units of metres) an uncertainty of a few metres. This is suspiciously small, and is inconsistent with the latitude and longitude errors, which are listed as 2.7 and 4.7 km on the quakesearch webpage.
@emilyws1
Thanks very much for raising the issue. We have identified the issue here to be with the default style sheet not performing the expected unit conversions when going from seiscompxml to quakeml.
seiscompxml uses kilometer for origin depth uncertainty, origin uncertainty and confidence ellipsoid, quakeml uses meters for these.
We have raised this issue with Gempa and we will apply a fix once we receive a correct style sheet to use.
Please let us know if you have identified other issues in the quakemls.
Updated into Gempa community forum Fix to come
Kia ora team,
A student highlighted an inconsistency in the uncertainties provided for a couple of events:
Returns:
The units for these values should be SI (m), but the confidence ellipsoid seems inconsistent with the stated horizontal uncertainties. I wondered if something had gone wrong parsing from SC3ML to QuakeML or somewhere else in the chain?