UNAVCO Logo
 
 
Projects Help with Projects Project Support & Services Borehole Services Siting Permitting Drilling BSM Data Processing Overview BSM Literature 2012 Strainmeter Workshop Report
Borehole Strainmeter (BSM) Data Processing Overview

How To Process Borehole Strainmeter Data

The basic steps involved in processing PBO borehole strainmeter data are outlined in the flowchart below. Full details are given in How PBO Strain Data Are Processed and summarized in Strainmeter Products Summary. The PBO Data Critical Design Review describes the processing techniques used for all PBO data sets.

For PBO strainmeters extensional strain is positive.

1. Download Raw Data

The raw data from each of the four gauges within a PBO strainmeter are available in 20 Hz, 1 Hz and 5-minute interval data rates. All BSM data are available in Bottle and SEED format (see PBO SEED channel codes) plus the low frequency data set is available in ASCII format.

2. Linearize Raw Data

The raw measurements from each gauge are outputs from capacitance bridges. To convert the gauge output, in counts, to strain along the axis of the gauge the following equation is used.

et=((Rt/1e+8)/(1-(Rt/1e+8))) - (Ro/1e+8)/(1-(Ro/1e+8)))*G/D

where et is the linearized strain at time t, Rt is the raw data point, Ro is a reference data point, G is the reference gap within the strainmeter (100 or 200 microns) and D is the instrument diameter (87 mm). All Level 2 time series are contain the linearized strain data.

3. Edit Data

Before tides and long term trends can be extracted from the gauge data non-tectonic and tectonic steps must be removed. Examples of such signals are coseismic offsets, instrument restarts or lightning strikes. UNAVCO's BSMAC isolates and finds such signals and flags them in the Level 2 Data Product.

4. Estimate the Earth and Ocean Load Tidal Signal

Once large steps have been removed from the time series the clean data can be input to a tidal analysis program such as BAYTAP-G and the amplitude and phase for the main tidal constituents determined. Each gauge is analyzed separately. When the tidal signal is being estimated for the first time it is often necessary to repeat steps three and four as small data steps will become apparent once the tides have been removed. These amplitudes and phases are calculated by UNAVCO's BSMAC and provided as a Level 2 Data Product both in an XML file and as a simplified ASCII file for each BSM.

5. Estimate The Barometric Pressure Response Coefficient

The strain changes are assumed to be linearly related to change in barometric pressure response. One linear response coefficient is estimated for each gauge . Typical responses range from a few to 10 nanostrain per millibar. These responses are also provided as a Level 2 Data Product.

6. Estimate the Long Term Trends

Estimating the long-term trends is one of the most difficult steps in processing BSM data as care must be taken to not remove the signal of interest. The expected long-term trend is one of compression as the borehole tends towards closure. For simple cases the trends can be forward modeled using an exponential plus linear trend. These forward modeled trends are provided as Level 2 Data Products. Alternative methods for detrending are high, low or band pass filtering to isolate the signal for interest.

7. Generate Areal and Shear Strains

The final step in BSM data processing is the calculation of areal and shear strains. In this step the gauge measurements are combined to form areal (Eee+Enn), differential extension (Eee-Enn) and tensor shear strain 2Ene . Currently the areal and shear strain measurements are generated using the manufacturer's calibrations and the tidal calibrations as published by Roeloffs (2010). These matrices are listed in the Level 2 XML files and a simplified ASCII version.

Requesting High Rate Processed Strain Data

High rate (1 sps) processed strain data is generated by UNAVCO after a M7 event any where in the world, any geophysical event of interest within the Plate Boundary Observatory and upon user request. The processed data set contains all the corrections that are normally seen in the 5 minute Level 2 data set : tidal, barometric pressure and borehole trend corrections plus offsets and data quality flags. In addition the files will contain tiltmeter, pore pressure and high rate barometric pressure data if they are collected at the site.

The data will be placed in an anonymous FTP directory for download and an email sent to the data requestor. Data sets generated following events may be updated with additional data in the hours following the event if required. The format is that of a space delimited ASCII file.

High Rate Data Format

For each event there is one tar file per strainmeter. Within each tar are five gzipped files, one for each gauge (CH0,CH1,CH2 and CH3) and one for tensor strain.

File Name Convention:

  • BNUM.CH0.event.txt.gz
  • BNUM.CH1.event.txt.gz
  • BNUM.CH2.event.txt.gz
  • BNUM.CH3.event.txt.gz
  • BNUM.tensor.LAB.event.txt.gz
  • BNUM.tensor.PBO_01.event.txt.gz

BNUM refers to the 4-character ID, e.g., B004 and "event" refers to the event name. E.g., B004.CH0.20110311Tohoku.txt.gz.

File Format

The format for the CH0, CH1, CH2 and CH3 files is:

Column

Description

Units

1

Date Time

UTC

2

Gauge measurement

Digital counts

3

Gauge strain

Microstrain

4

Interpolated missing values

Microstrain

5

Tide

Microstrain

6

Barometric correction

Microstrain

7

Barometric pressure

Millibars

8

Pore pressure

Millibars

9

Version

Generation date

The format for the tensor strain file is:

Column

Description

Units

1

DateTime

UTC

2

Eee+Enn Strain

Microstrain

3

Eee+Enn Tide

Microstrain

4

Eee+Enn Interpolated missing values

Microstrain

5

Eee+Enn barometric correction

Microstrain

6

Eee-Enn Strain

Microstrain

7

Eee-Enn Interpolated missing values

Microstrain

8

Eee-Enn Tide

Microstrain

9

Eee-Enn Barometric correction

Microstrain

10

2Ene Strain

Microstrain

11

2Ene Interpolated missing values

Microstrain

12

2Ene Tide

Microstrain

13

2Ene Barometric correction

Microstrain

15

Barometric Pressure

millibar

16

Pore Pressure

millibar

17

Version

Generation date

18

Tilt X

Microradian

19

Tilt Y

Microradian

The null value is 999999. If 1-sps barometric data are available it will be given otherwise there will only be 30 minute interval data. Pore Pressure and tiltmeter data are only available for a subset of sites. The "Interpolation" columns contained interpolated values to replace the null values. This column should be substracted from the linear strain to remove the null values. Barometric, pore and tiltmeter time-stamps have been rounded to the nearest second.

Last modified: Sunday, 11-May-2014 17:56:39 UTC

 

Sponsored by

National Science Foundation Logo National Aeronautics and Space Administration Logo