SAGE Facility MUSTANG metrics Web Service Documentation

pct_above_nhnm Percentage of PDF values above the New High Noise Model

Summary

Once hourly power spectral density (PSD) measurements are instrument-corrected, discrete period-power values for the day are binned into a probability density function (PDF) matrix. The number of PDF values that exceed the Peterson New High Noise Model (NHNM) (Peterson, 1993) for each discrete period are presented as a percentage of total PDF values.

Peterson, J, 1993, Observations and Modeling of Seismic Background Noise, U.S.G.S. OFR-93-322

Uses

An easy way to identify channels with noise problems is to search for pct_above_nhnm values larger than some threshold.

Data Analyzed

Traces – one N.S.L.C (Network.Station.Location.Channel) per measurement
Window – up to 24 hours depending on the completeness of data
Data SourceIRIS miniSEED archive or IRIS PH5 archive

SEED Channel Types – ?H?, ?L?, ?N?, ?G?, ?P? | High Gain, Low Gain, Accelerometer, Gravimeter, Geophone | excluding very long period channels

Algorithm

  • Request 24 hours of time series data for the current N.S.L.C.
  • Calculate PSDs (see noise-psd help)
  • Correct PSD curves for instrument response.
  • Bin PSD values into a PDF matrix using existing period values and 1 dB bins.
  • For each period in the matrix,
    • Make a running count of power values that exceed the NHNM,
    • Divide the final count by the total number of values in the matrix and multiply this ratio by 100 to get percent:
      pct_above_nhnm = 100 * ( PDF_matrix_values_above_NHNM / total_PDF_matrix_values )
      

Metric Values Returned

value – percentage of PDF values above the NHNM curve
target – the trace analyzed, labeled as N.S.L.C.Q (Network.Station.Location.Channel.Quality)
start – time of the first sample of time series data for that day in UTC
end – time of the last sample of time series data for that day in UTC
lddate – date/time the measurement was made and loaded into the MUSTANG database (UTC)

Notes

Instrument response has been removed.

Contact

See Also

Updated