GeoNet / help

An issues repo for technical help questions.
6 stars 3 forks source link

RIZ issue with broadband vs. accelerometer comparison #85

Closed elidana closed 2 years ago

elidana commented 2 years ago

Jesse Hutchinson (University of Canterbury) has identified an issue with RIZ seismic data when comparing data from the broadband and strong motion sensors installed at the station.

His procedure is explained in https://github.com/GeoNet/help/issues/84#issuecomment-931798699

processed Power Spectral Amplitude plots show a consistent difference between RINZ strong motion and broadband data for earthquakes with M>5 in the time period 2019-2020 pSA_diff_RIZ_datetime

It is likely that this difference is due to one of the sensors going flat during that time period. Further investigations to confirm the cause are needed.

/cc @SquirrelKnight, to help pinpoint time period and possible cause, can you please post here the list of events used to produce this figure?

SquirrelKnight commented 2 years ago

Here are the events used for the above measurements:

'2012p553582', '2012p595069', '2012p647755', '2012p681111', '2012p733490', '2012p766923', '2012p767122', '2020p263458', '2020p269386', '2020p332789', '2020p489837', '2020p498617', '2020p525251', '2020p548309', '2020p830671', '2020p926977', '2020p979274'

salichon commented 2 years ago

RIZ is Raoul island Remotely accessble stations

salichon commented 2 years ago
salichon commented 2 years ago

There is an - obvious- PSA/PGA difference in 2020 by RIZ Raoul island between HH/HN See belpw (no need to compare the PSAs Ratio is ~time 4

Event 2020p498617

<?xml version="1.0" encoding="UTF-8" standalone="yes"?>         <?xml version="1.0" encoding="UTF-8" standalone="yes"?>
<!DOCTYPE earthquake SYSTEM "stationlist.dtd">                  <!DOCTYPE earthquake SYSTEM "stationlist.dtd">
<stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml"   <stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml"
  <station code="RIZ" name="RIZ" netid="NZ" insttype="Trilliu |   <station code="RIZ" name="RIZ" netid="NZ" insttype="FBA-ES-
    <comp name="10.HHE">                                      |     <comp name="20.HNE">
      <acc value="3.2882998009" flag="0"/>                    |       <acc value="0.7794464103" flag="0"/>
      <vel value="1.2413733920" flag="0"/>                    |       <vel value="0.3010102158" flag="0"/>
      <psa03 value="5.2516239164" flag="0"/>                  |       <psa03 value="1.3026295004" flag="0"/>
      <psa10 value="0.5281061817" flag="0"/>                  |       <psa10 value="0.1306030042" flag="0"/>
      <psa30 value="0.0349679485" flag="0"/>                  |       <psa30 value="0.0086528750" flag="0"/>
    </comp>                                                         </comp>
    <comp name="10.HHN">                                      |     <comp name="20.HNN">
      <acc value="5.3235685639" flag="0"/>                    |       <acc value="1.4449198090" flag="0"/>
      <vel value="1.5349515077" flag="0"/>                    |       <vel value="0.3551027246" flag="0"/>
      <psa03 value="5.6611978434" flag="0"/>                  |       <psa03 value="1.4561612866" flag="0"/>
      <psa10 value="0.3342218803" flag="0"/>                  |       <psa10 value="0.0880081106" flag="0"/>
      <psa30 value="0.0292224447" flag="0"/>                  |       <psa30 value="0.0076309110" flag="0"/>
    </comp>                                                         </comp>
    <comp name="10.HHZ">                                      |     <comp name="20.HNZ">
      <acc value="3.5350992197" flag="0"/>                    |       <acc value="0.4533266631" flag="0"/>
      <vel value="0.6877658960" flag="0"/>                    |       <vel value="0.1281977146" flag="0"/>
      <psa03 value="1.2282583205" flag="0"/>                  |       <psa03 value="0.2997003001" flag="0"/>
      <psa10 value="0.1445904646" flag="0"/>                  |       <psa10 value="0.0345638400" flag="0"/>
      <psa30 value="0.0115920347" flag="0"/>                  |       <psa30 value="0.0028616972" flag="0"/>
    </comp>                                                         </comp>
  </station>                                                      </station>
</stationlist>                                                  </stationlist>

event 2020p979274

> <?xml version="1.0" encoding="UTF-8" standalone="yes"?>         <?xml version="1.0" encoding="UTF-8" standalone="yes"?>
> <!DOCTYPE earthquake SYSTEM "stationlist.dtd">                  <!DOCTYPE earthquake SYSTEM "stationlist.dtd">
> <stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml"   <stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml"
>   <station code="RIZ" name="RIZ" netid="NZ" insttype="Trilliu |   <station code="RIZ" name="RIZ" netid="NZ" insttype="FBA-ES-
>     <comp name="10.HHE">                                      |     <comp name="20.HNE">
>       <acc value="0.4941437803" flag="0"/>                    |       <acc value="0.1170454285" flag="0"/>
>       <vel value="0.1115322176" flag="0"/>                    |       <vel value="0.0282150623" flag="0"/>
>       <psa03 value="0.4723697473" flag="0"/>                  |       <psa03 value="0.1200957546" flag="0"/>
>       <psa10 value="0.0250678201" flag="0"/>                  |       <psa10 value="0.0062046699" flag="0"/>
>       <psa30 value="0.0022416171" flag="0"/>                  |       <psa30 value="0.0005787207" flag="0"/>
>     </comp>                                                         </comp>
>     <comp name="10.HHN">                                      |     <comp name="20.HNN">
>       <acc value="0.3657536328" flag="0"/>                    |       <acc value="0.0811721977" flag="0"/>
>       <vel value="0.0933475966" flag="0"/>                    |       <vel value="0.0241417859" flag="0"/>
>       <psa03 value="0.5231511377" flag="0"/>                  |       <psa03 value="0.1355195799" flag="0"/>
>       <psa10 value="0.0209135042" flag="0"/>                  |       <psa10 value="0.0053807146" flag="0"/>
>       <psa30 value="0.0019890443" flag="0"/>                  |       <psa30 value="0.0005128434" flag="0"/>
>     </comp>                                                         </comp>
>     <comp name="10.HHZ">                                      |     <comp name="20.HNZ">
>       <acc value="0.2302177579" flag="0"/>                    |       <acc value="0.0649731152" flag="0"/>
>       <vel value="0.0510248189" flag="0"/>                    |       <vel value="0.0134712057" flag="0"/>
>       <psa03 value="0.2288046422" flag="0"/>                  |       <psa03 value="0.0575988433" flag="0"/>
>       <psa10 value="0.0118311645" flag="0"/>                  |       <psa10 value="0.0029701454" flag="0"/>
>       <psa30 value="0.0008189927" flag="0"/>                  |       <psa30 value="0.0002130991" flag="0"/>
>     </comp>                                                         </comp>

---------------------------------------

GLKZ | RIZ HH psa computation shows consistent PSA values in scale Event 2020p498617

<!DOCTYPE earthquake SYSTEM "stationlist.dtd">
<stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml">
  <station code="GLKZ" name="GLKZ" netid="NZ" insttype="Trillium 120QA" lat="-29.262" lon="-177.92">
    <comp name="10.HHE">
      <acc value="2.1806454679" flag="0"/>
      <vel value="0.5180368734" flag="0"/>
      <psa03 value="3.8059693891" flag="0"/>
      <psa10 value="0.3447317382" flag="0"/>
      <psa30 value="0.0188931208" flag="0"/>
    </comp>
    <comp name="10.HHN">
      <acc value="1.9885158498" flag="0"/>
      <vel value="0.5018293198" flag="0"/>
      <psa03 value="3.1817822577" flag="0"/>
      <psa10 value="0.3104649784" flag="0"/>
      <psa30 value="0.0162603001" flag="0"/>
    </comp>
    <comp name="10.HHZ">
      <acc value="1.5119721137" flag="0"/>
      <vel value="0.2928724373" flag="0"/>
      <psa03 value="1.7520792673" flag="0"/>
      <psa10 value="0.1064284764" flag="0"/>
      <psa30 value="0.0088276304" flag="0"/>
    </comp>
  </station>
  <station code="RIZ" name="RIZ" netid="NZ" insttype="Trillium 120QA" lat="-29.2449" lon="-177.929">
    <comp name="10.HHE">
      <acc value="3.2882998009" flag="0"/>
      <vel value="1.2413733920" flag="0"/>
      <psa03 value="5.2516239164" flag="0"/>
      <psa10 value="0.5281061817" flag="0"/>
      <psa30 value="0.0349679485" flag="0"/>
    </comp>
    <comp name="10.HHN">
      <acc value="5.3235685639" flag="0"/>
      <vel value="1.5349515077" flag="0"/>
      <psa03 value="5.6611978434" flag="0"/>
      <psa10 value="0.3342218803" flag="0"/>
      <psa30 value="0.0292224447" flag="0"/>
    </comp>
    <comp name="10.HHZ">
      <acc value="3.5350992197" flag="0"/>
      <vel value="0.6877658960" flag="0"/>
      <psa03 value="1.2282583205" flag="0"/>
      <psa10 value="0.1445904646" flag="0"/>
      <psa30 value="0.0115920347" flag="0"/>
    </comp>
  </station>
</stationlist>
salichon commented 2 years ago

Obsidian datalogger related ?

FYI @SquirrelKnight @elidana WIP

salichon commented 2 years ago
salichon commented 2 years ago
salichon commented 2 years ago

and a consistent 1/4 ratio as identified by jesse between HH and BN on RIZ until a n strong motion issue along 2018 until a replacement in 2019 sept FYI @SquirrelKnight @elidana @ozym


Need to check if the HH Vs BN/HN ratio is still 4 in 2021 (after a Weak inst visit as 2020 Nov 10th)

event Sep 2021 2021p715493 M5.1 and August 2021 2021p638513 M 5.5

salichon commented 2 years ago

Results for 2021p715493 M5.5 GLKZ and RIZ Weak consistent output ~1-2%g Vs RIZ strong motion 0.3%g
Ratio RIZ HH Vs SM is 1/4 still

<stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml"   <stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml"
  <station code="GLKZ" name="GLKZ" netid="NZ" insttype="Trill     <station code="GLKZ" name="GLKZ" netid="NZ" insttype="Trill
    <comp name="10.HHE">                                            <comp name="10.HHE">
      <acc value="1.2434342470" flag="0"/>                            <acc value="1.2434342470" flag="0"/>
      <vel value="0.3226558428" flag="0"/>                            <vel value="0.3226558428" flag="0"/>
      <psa03 value="1.1799961255" flag="0"/>                          <psa03 value="1.1799961255" flag="0"/>
      <psa10 value="0.1873236385" flag="0"/>                          <psa10 value="0.1873236385" flag="0"/>
      <psa30 value="0.0122821212" flag="0"/>                          <psa30 value="0.0122821212" flag="0"/>
    </comp>                                                         </comp>
    <comp name="10.HHN">                                            <comp name="10.HHN">
      <acc value="1.0757444905" flag="0"/>                            <acc value="1.0757444905" flag="0"/>
      <vel value="0.2776559241" flag="0"/>                            <vel value="0.2776559241" flag="0"/>
      <psa03 value="1.1076307169" flag="0"/>                          <psa03 value="1.1076307169" flag="0"/>
      <psa10 value="0.1566439131" flag="0"/>                          <psa10 value="0.1566439131" flag="0"/>
      <psa30 value="0.0104541844" flag="0"/>                          <psa30 value="0.0104541844" flag="0"/>
    </comp>                                                         </comp>
    <comp name="10.HHZ">                                            <comp name="10.HHZ">
      <acc value="1.0311887968" flag="0"/>                            <acc value="1.0311887968" flag="0"/>
      <vel value="0.1782112942" flag="0"/>                            <vel value="0.1782112942" flag="0"/>
      <psa03 value="0.7133046028" flag="0"/>                          <psa03 value="0.7133046028" flag="0"/>
      <psa10 value="0.0506322085" flag="0"/>                          <psa10 value="0.0506322085" flag="0"/>
      <psa30 value="0.0068765641" flag="0"/>                          <psa30 value="0.0068765641" flag="0"/>
    </comp>                                                         </comp>
  </station>                                                      </station>
  <station code="RIZ" name="RIZ" netid="NZ" insttype="Trilliu |   <station code="RIZ" name="RIZ" netid="NZ" insttype="FBA-ES-
    <comp name="10.HHE">                                      |     <comp name="20.HNE">
      <acc value="1.5805146737" flag="0"/>                    |       <acc value="0.3822916356" flag="0"/>
      <vel value="0.6992911070" flag="0"/>                    |       <vel value="0.1721719665" flag="0"/>
      <psa03 value="3.9067508958" flag="0"/>                  |       <psa03 value="0.9890861015" flag="0"/>
      <psa10 value="0.3648450089" flag="0"/>                  |       <psa10 value="0.0904114417" flag="0"/>
      <psa30 value="0.0254664582" flag="0"/>                  |       <psa30 value="0.0064493938" flag="0"/>
    </comp>                                                         </comp>
    <comp name="10.HHN">                                      |     <comp name="20.HNN">
      <acc value="2.5137795811" flag="0"/>                    |       <acc value="0.5764419711" flag="0"/>
      <vel value="0.6970251287" flag="0"/>                    |       <vel value="0.1708662387" flag="0"/>
      <psa03 value="2.7988747423" flag="0"/>                  |       <psa03 value="0.7049039446" flag="0"/>
      <psa10 value="0.1880126050" flag="0"/>                  |       <psa10 value="0.0498220628" flag="0"/>
      <psa30 value="0.0116712685" flag="0"/>                  |       <psa30 value="0.0030854582" flag="0"/>
    </comp>                                                         </comp>
    <comp name="10.HHZ">                                      |     <comp name="20.HNZ">
      <acc value="1.1997534487" flag="0"/>                    |       <acc value="0.3079534936" flag="0"/>
      <vel value="0.2771769896" flag="0"/>                    |       <vel value="0.0701502819" flag="0"/>
      <psa03 value="1.4687793185" flag="0"/>                  |       <psa03 value="0.3726367741" flag="0"/>
      <psa10 value="0.0740232073" flag="0"/>                  |       <psa10 value="0.0187514684" flag="0"/>
      <psa30 value="0.0093899779" flag="0"/>                  |       <psa30 value="0.0023705178" flag="0"/>
    </comp>                                                         </comp>
  </station>                                                      </station>
</stationlist>                                                  </stationlist>

image


Check on event 2021p876265 <100Kmdistance M5.7 depth 30km OT 2021-11-21T05:12:29Z using MMKS and RIZ strong motion (measuer should consistent)

As shown Below There is a PGA ratio in the measure netween MMKS and RIZ strong motion

!DOCTYPE earthquake SYSTEM "stationlist.dtd">
<stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml">
  <station code="MMKS" name="MMKS" netid="NZ" insttype="Titan" lat="-29.2673" lon="-177.908">
    <comp name="20.HNE">
      <acc value="2.3375445330" flag="0"/>
      <vel value="1.5666094641" flag="0"/>
      <psa03 value="11.1049878581" flag="0"/>
      <psa10 value="1.4099762785" flag="0"/>
      <psa30 value="0.1446413040" flag="0"/>
    </comp>
    <comp name="20.HNN">
      <acc value="2.8647678642" flag="0"/>
      <vel value="2.6098697229" flag="0"/>
      <psa03 value="8.9490359984" flag="0"/>
      <psa10 value="3.5219782617" flag="0"/>
      <psa30 value="0.2265863963" flag="0"/>
    </comp>
    <comp name="20.HNZ">
      <acc value="0.8258273479" flag="0"/>
      <vel value="0.3882073192" flag="0"/>
      <psa03 value="2.9999018243" flag="0"/>
      <psa10 value="0.2759690486" flag="0"/>
      <psa30 value="0.0506806627" flag="0"/>
    </comp>
  </station>
  <station code="RIZ" name="RIZ" netid="NZ" insttype="FBA-ES-T-ISO" lat="-29.2449" lon="-177.929">
    <comp name="20.HNE">
      <acc value="0.5461324014" flag="0"/>
      <vel value="0.2547292779" flag="0"/>
      <psa03 value="1.6095682664" flag="0"/>
      <psa10 value="0.2454820776" flag="0"/>
      <psa30 value="0.0254289251" flag="0"/>
    </comp>
    <comp name="20.HNN">
      <acc value="0.6653214861" flag="0"/>
      <vel value="0.3193675578" flag="0"/>
      <psa03 value="2.2627113708" flag="0"/>
      <psa10 value="0.2506421204" flag="0"/>
      <psa30 value="0.0179910220" flag="0"/>
    </comp>
    <comp name="20.HNZ">
      <acc value="0.2648957490" flag="0"/>
      <vel value="0.0937906178" flag="0"/>
      <psa03 value="0.7140915188" flag="0"/>
      <psa10 value="0.0917864367" flag="0"/>
      <psa30 value="0.0093684920" flag="0"/>
    </comp>
  </station>
</stationlist>

---------- RIZ strong, Bias confirmed

salichon commented 2 years ago

FYI @elidana @mnaguit

Perhaps to be communicated into a delta meeting

salichon commented 2 years ago

Summary most likely scenarios (today)

FYI @SquirrelKnight @elidana

salichon commented 2 years ago

Output:

<?xml version="1.0" encoding="UTF-8" standalone="yes"?>
<!DOCTYPE earthquake SYSTEM "stationlist.dtd">
<stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml">
  <station code="RIZ" name="RIZ" netid="NZ" insttype="FBA-ES-T-ISO" lat="-29.2449" lon="-177.929">
    <comp name="20.HNE">
      <acc value="2.1806946883" flag="0"/>
      <vel value="1.0171284216" flag="0"/>
      <psa03 value="6.4269707508" flag="0"/>
      <psa10 value="0.9802045469" flag="0"/>
      <psa30 value="0.1015371409" flag="0"/>
    </comp>
    <comp name="20.HNN">
      <acc value="2.6566140876" flag="0"/>
      <vel value="1.2752276475" flag="0"/>
      <psa03 value="9.0349568291" flag="0"/>
      <psa10 value="1.0008084841" flag="0"/>
      <psa30 value="0.0718377561" flag="0"/>
    </comp>
    <comp name="20.HNZ">
      <acc value="1.0577229064" flag="0"/>
      <vel value="0.3745039464" flag="0"/>
      <psa03 value="2.8513517537" flag="0"/>
      <psa10 value="0.3665012273" flag="0"/>
      <psa30 value="0.0374081832" flag="0"/>
    </comp>
  </station>
</stationlist>

TO BE compared to MMKS strong for the same event about the same place for instance:cf https://github.com/GeoNet/help/issues/85#issuecomment-977450927 (this could be using also the colocated RIZ HH values too)

<stationlist created="" xmlns="ch.ethz.sed.shakemap.usgs.xml">
  <station code="MMKS" name="MMKS" netid="NZ" insttype="Titan" lat="-29.2673" lon="-177.908">
    <comp name="20.HNE">
      <acc value="2.3375445330" flag="0"/>
      <vel value="1.5666094641" flag="0"/>
      <psa03 value="11.1049878581" flag="0"/>
      <psa10 value="1.4099762785" flag="0"/>
      <psa30 value="0.1446413040" flag="0"/>
    </comp>
    <comp name="20.HNN">
      <acc value="2.8647678642" flag="0"/>
      <vel value="2.6098697229" flag="0"/>
      <psa03 value="8.9490359984" flag="0"/>
      <psa10 value="3.5219782617" flag="0"/>
      <psa30 value="0.2265863963" flag="0"/>
    </comp>
    <comp name="20.HNZ">
      <acc value="0.8258273479" flag="0"/>
      <vel value="0.3882073192" flag="0"/>
      <psa03 value="2.9999018243" flag="0"/>
      <psa10 value="0.2759690486" flag="0"/>
      <psa30 value="0.0506806627" flag="0"/>
    </comp>
  </station>

Results are consistent and makes sense - I don t think more checking is required from my side users and others coudl give it a go.

FYI @elidana @mnaguit @SquirrelKnight

...almost there ... :)

salichon commented 2 years ago

@SquirrelKnight This instrument related issue will be addressed soon while patching the metadata. it will encompass a series of sites as per https://github.com/GeoNet/tickets/issues/9943 over identified period of times. (related to Obsidian datalogger 4X response)

cheers FYI @elidana @mnaguit

SquirrelKnight commented 2 years ago

That is wonderful to hear. Thank you for all of the hard work!

On Wed, Feb 16, 2022 at 12:55 PM Jerome Salichon @.***> wrote:

@SquirrelKnight https://github.com/SquirrelKnight This instrument related issue will be addressed soon while patching the metadata. it will encompass a series of sites as per GeoNet/tickets#9943 https://github.com/GeoNet/tickets/issues/9943 over identified period of times. (related to Obsidian datalogger 4X response)

cheers FYI @elidana https://github.com/elidana @mnaguit https://github.com/mnaguit

— Reply to this email directly, view it on GitHub https://github.com/GeoNet/help/issues/85#issuecomment-1042302054, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJ5LPP67MF2Q4MZYIUQLUYDU3QFLLANCNFSM5FDXIQ4A . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

You are receiving this because you were mentioned.Message ID: @.***>

salichon commented 2 years ago

PR was merged correction applied today @SquirrelKnight documentaion to follow up with the sites affected

salichon commented 2 years ago

Re-open if additional issues