Ambient Noise Probability Density Functions



Contents

Introduction

Characterizing Sources of Noise and Signal in the PDFs

Processing Methods

References

Acknowledgements



Introduction


A new system for analyzing data quality is now available to the seismology community allowing users to evaluate the long-term seismic noise levels for any broadband seismic data channel streaming into the buffer of uniform data (BUD) within the Incorporated Research Institutions for Seismology (IRIS) data management system (DMS). BUD is the IRIS DMS's acronym for the online data cache from which the DMC distributes near-real time miniSEED data holdings prior to formal archiving.

The new noise processing software uses a probability density function (PDF) to display the distribution of seismic power spectral density (PSD) and has been implemented
against the entire continuous data-stream available within the BUD utilizing the QUACK framework. QUACK is the software package at the IRIS DMC, responsible for managing the quality control (QC) of the real-time seismic data flowing into the BUD (see http://www.iris.washington.edu/servlet/quackquery/).

This noise processing system is unique in that there is no need to screen the data for  earthquakes, system glitches or general data artifacts, as is commonly done in seismic noise analysis. Instead with this new analysis, system transients map into a low-level background probability while ambient noise conditions reveal themselves as high probability occurrances.  In fact, examination of artifacts related to station operation and episodic cultural noise allows us to estimate both the overall station quality and a baseline level of earth noise at each site.

PDF noise plots are useful for characterizing the current and past performance of existing broadband sensors, for detecting operational problems within the recording system, and for evaluating the overall quality of data for a particular station. The advantages of this new approach include:
       
        1) provides an analytical view representing the true ambient noise levels rather than a simple absolute minimum;
       
        2) provides an assessment of the overall health of the instrument/station; and,
       
        3) provides an assessment of the health of recording and telemetry systems.


Below are PDF examples, with some artifacts and signals identified.

HLID - Typical ANSS station ~10km from  Hailey Idaho.
LTX - Lajitas TX. This station was instrumental in the orginal Peterson Low Noise Model, however due to increased cultural noise (0.1-1s, 1-10Hz) the highest probabilty power levels (i.e. mode, black line) are now significantly higher than the LNM. The minimum (red line) will approach the LNM <2% of the time.

Additional features and artifacts observed in the PDFs will be discussed below.  (For a more detailed discussion of the PDFs see McNamara and Buland (2004).)

Back to top


Characterizing Sources of Noise and Signal in the PDFs

Cultural noise

The most common source of seismic noise is from the actions of human beings at or near the surface of the Earth. This is often referred to as “cultural noise” and originates primarily from the coupling of traffic and machinery energy into the earth. Cultural noise propagates mainly as high-frequency surface waves (>1-10Hz,1-0.1s) that attenuate within several kilometers in distance and depth. For this reason cultural noise will generally be significantly reduced in boreholes, deep caves and tunnels. Cultural noise shows very strong diurnal variations and has characteristic frequencies depending on the source of the disturbance.

        Examples:
HLID - automobile traffic along a dirt road only 20 meters from station HLID creates a 20-30dB increase in power at about 0.1 sec period (10Hz). This type of cultural noise is observable in the PDFs as a region of low probability at high frequencies (1-10Hz, 0.1-1s).

ANMO - Borehole station in Albuquerque NM, with very low noise levels. For periods >1s the minimum closely tracks the Peterson LNM.

DWPF - Borehole station at Disney World Florida. Cultural noise levels are high due to human activity and microseism are high due to proximity to the Atlantic ocean.

ISCO - Abandoned mine tunnel station in Idaho Springs CO. Noise levels are low with minimal variations in cultural noise and microseisms.
Earthquakes

Our approach differs from many previous noise studies in that we make no attempt to screen the continuous waveforms to eliminate body and surface waves from naturally occurring earthquakes. Earthquake signals are included in our processing because they are generally low probability occurrences even at low power levels (small magnitude events). We are interested in the true noise that a given station will experience, thus we include all signals. For example, including events tells us something about the probability of teleseismic signals being obscured by small local events as well as various noise sources. Large teleseismic earthquakes can produce powers above ambient noise levels across the entire spectrum and are dominated by surface waves >10s, while small events dominate the short period, <1s. Earthquakes are observed in the PDFs as low probability smeared signal at short and long periods.

        Examples:
HLID - Body waves occur as low probabily signal in the 1sec range while surface waves are generally higher power at longer periods.

DWPF - Noise levels are high enough (0.1-1s) that body waves from smaller earthquakes are not readily detectable.

System artifacts

Since we make no attempt to screen waveforms for system transients such as data gaps and sensor glitches, the PDF plots contain numerous system generated artifacts that can be very useful for network quality control purposes. We have attempted to determine the source of several coherent, high power, low probability noise artifacts in the PDF plots. Several artifacts in the PDFs are easily explained and may be useful to the network operator. For example, data-gaps (due to telemetry drop-outs) and automatic mass re-centers (necessitated by “drift” in sensor mass position) are easily identifiable in the PDFs.  Should the probability of mass re-centering and/or telemetry drop-outs drastically increase, a remote network operator could readily diagnose the problem.

        Examples:
HLID  - Automatic mass re-centering and calibration pulses show up as low probability occurances in the PDF.

BINY - Binghamton NY, BHN component. Numerous mass recenters due to unstable sensor pad. We can determine that mass re-centers are automatically issued aproximately 1% of the time.

LTX - Lajitas TX.  Telemetry dropouts are somewhat common (~1-2%).

Back to top


Processing Methods

Power Spectral Density

Employing the algorithm used to develop the USGS Albuquerque Seismological Laboratory (ASL) low noise model (LNM; Peterson, 1993), we compute the power spectral density (PSD) for broadband stations streaming into the IRIS DMC BUD system.


In this algorithm, hour-long, continuous, and over-lapping (50%) time series segments are processed. (N.B. There is no removal of earthquakes, system transients and/or data glitches.) The instrument transfer function is removed from each segment, yielding ground acceleration (for easy comparison to the LNM).  Each hour-long time series is divided into 13 segments, each about 15 minutes long and overlapping by 75%, with each segment processed by:        
        1) removing the mean;
        2) removing the long period trend;
        3) tapering using a 10% sine function; and
        4) transforming using an FFT algorithm (Bendat and Piersol, 1971).

Segments are then averaged to provide a PSD for each one-hour time series segment. 


Probability Density Functions

For each channel, raw frequency distributions are constructed by gathering individual PSDs in the following manner:
       1) binning periods in 1/8 octave intervals; and
       2) binning power in 1 dB intervals.

Each raw frequency distribution bin is then normalized by the total number of PSDs to construct a Probability Density Function (PDF). The probability of occurrence of a given power at a particular period is plotted for direct comparison to the Peterson high and low noise models (HNM, LNM).

We also compute and plot the minimum, mode, and maximum powers for each period bin. A wealth of seismic noise information can be obtained from this statistical view of broadband seismic noise. For a more detailed discussion of the noise processing methods see McNamara and Buland (2004).

Click on the following links for sample PDFs with comments:

HLID  LTX  ANMO  BINY  DWPF  ISCO

or Return to IRIS DMC PDFs

Back to top


References

Bendat, J.S. and A.G. Piersol (1971). Random data: analysis and measurement procedures. John Wiley and Sons, New York, 407p.

McNamara, D.E. and R.P. Buland, Ambient Noise Levels in the Continental United States, Bull. Seism. Soc. Am., 94, 4, 1517-1527, 2004.

Peterson, J., Observation and modeling of seismic background noise, U.S. Geol. Surv. Tech. Rept., 93-322, 1-95, 1993.

Back to top



Acknowledgements


The algorithm and initial software were first developed by Dan McNamara and Ray Buland, at the United States Geological Survey (USGS) as a part of the data and network quality control (QC) system for the Advanced National Seismic System (ANSS).  Further development, supported by IRIS through funds it receives from the National Science Foundation (NSF) allowed for Richard Boaz to develop the system for implementation against the BUD dataset. Additonal installation and implementation support was provided by Bruce Weertman and Tim Ahern at the IRIS DMC.

Webpage authors: D. McNamara (USGS ANSS) and R. Boaz (www.boazconsultancy.com)

Back to top