### 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 Source** – IRIS SEED 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.