The Bay Area Regional Deformation (BARD) network of continuously operating Global Positioning System (GPS) receivers monitors crustal deformation in the San Francisco Bay area (``Bay Area") and northern California (Murray et al., 1998). It is a cooperative effort of the BSL, the USGS, and several other academic, commercial, and governmental institutions. Started by the USGS in 1991 with 2 stations spanning the Hayward fault (King et al., 1995), BARD now includes 70 permanent stations (Figure 7.1). The principal goals of the BARD network are: 1) to determine the distribution of deformation in northern California across the wide Pacific-North America plate boundary from the Sierras to the Farallon Islands; 2) to estimate three-dimensional interseismic strain accumulation along the San Andreas fault (SAF) system in the Bay Area to assess seismic hazards; 3) to monitor hazardous faults and volcanoes for emergency response management; and 4) to provide infrastructure for geodetic data management and processing in northern California in support of related efforts within the surveying and other interested communities.
BARD currently includes 38 continuously operating stations in the Bay Area and northern California (Table 7.1), 14 near Parkfield along the central San Andreas fault, and 18 near the Long Valley caldera near Mammoth (Table 7.2). The BSL maintains 22 stations (including 2 with equipment provided by Lawrence Livermore National Laboratory (LLNL) and UC Santa Cruz). Other stations are maintained by the USGS (Menlo Park and Cascade Volcano Observatory), LLNL, Stanford University, UC Davis, UC Santa Cruz, and East Bay Municipal Utilities District, the City of Modesto, the National Geodetic Survey, Thales, Inc., and the Jet Propulsion Laboratory. Many of these stations are part of larger networks devoted to real-time navigation, orbit determination, and crustal deformation.
Between 1993 and 2001, the BSL acquired 29 Ashtech Z-12 and Micro-Z receivers from a variety of funding sources, including from federal (NSF and USGS), state (CLC), and private (EPRI) agencies. The network enhances continuous strain measurements in the Bay Area and includes several profiles between the Farallon Islands and the Sierra Nevada in order to better characterize the larger scale deformation field in northern California (Figure 7.1). Five more of the BSL receivers will be installed next year, 2 along the southern Hayward fault, and 3 as part of the NSF-funded Mini-PBO project establishing collocated GPS, and borehole strainmeter and seismometer observatories in the Bay Area (see Chapter 8).
The number of continuous GPS stations in northern California will dramatically increase over the next 5 years, with over 250 new site installations planned as part of the Plate Boundary Observatory (PBO) component of the NSF-funded Earthscope project. The BARD network will form the initial core of the northern California array, and the BSL recently received NSF funding to maintain 40 stations for an 18-month period. During this period, BARD and the other regional networks, such as SCIGN, BARGEN, and PANGA, will be developing plans to fold operation and maintenance of portions the existing networks into the PBO array at the end of 5 years. We are working closely with UNAVCO, Inc., who has primary responsibility for implementation of PBO, to facilitate this transition and are acting in an advisory role on siting issues for the new installations.
Today, raw and Rinex data files from the BSL stations and the other stations run by BARD collaborators are archived at the BSL/USGS Northern California Earthquake Data Center data archive maintained at the BSL (Romanowicz et al., 1994). The data are checked to verify their integrity, quality, completeness, and conformance to the RINEX standard, and are then made accessible, usually within 2 hours of collection, to all BARD participants and other members of the GPS community through Internet, both by anonymous FTP and by the World Wide Web (http://quake.geo.berkeley.edu/bard/).
Many of the BARD sites are classified as CORS stations by the NGS, which are used as reference stations by the surveying community. We coordinate efforts with surveying community at meetings of the Northern California GPS Users Group and the California Spatial Reference Center, and are currently developing plans to use the existing infrastructure at the NCEDC to provide a hub for a high-frequency real-time surveying network in the Bay Area. Data and ancillary information about BARD stations are also made compatible with standards set by the International GPS Service (IGS), which administers the global tracking network used to estimate precise orbits and has been instrumental in coordinating the efforts of other regional tracking networks. The NCEDC also retrieves data from other GPS archives, such as at SIO, JPL, and NGS, in order to provide a complete archive of all high-precision continuous GPS measurements collected in northern California.
The typical configuration of a BSL continuous GPS station installation has been described in detail in previous annual reports. We here provide a brief description and highlight some of the recent changes. During July 2002-June 2003, we performed maintenance on existing BARD stations, installed a new station, assisted collaborators with the installation of two new stations, and prepared for new stations near the Hayward fault, on the San Francisco peninsula, and north Bay area regions.
Each BSL BARD station uses a low-multipath choke-ring antenna, most of which are mounted to a reinforced concrete pillar approximately 0.5-1.0 meter above local ground level. The reinforcing steel bars of the pillar are drilled and cemented into rock outcrop to improve long-term monument stability. A low-loss antenna cable is used to minimize signal degradation on the longer cable setups that normally would require signal amplification. Low-voltage cutoff devices are installed to improve receiver performance following power outages. Most use Ashtech Z-12 receivers programmed to record data once every 30 seconds, observing up to 12 satellites simultaneously at elevations down to the horizon. The antennas are equipped with SCIGN antenna adapters and hemispherical domes, designed to provide security and protection from weather and other natural phenomenon, and to minimize differential radio propagation delays.
Data from most BSL-maintained stations are collected at 30-second intervals and transmitted continuously over serial connections (Table 7.1). Station TIBB uses a direct radio link to Berkeley, and MODB uses VSAT satellite telemetry. Nineteen stations use frame relay technology, either alone or in combination with radio telemetry. Thirteen GPS stations are collocated with broadband seismometers and Quanterra data collectors (Table 3.2). With the support of IRIS we developed software that converts continuous GPS data to MiniSEED opaque blockettes that are stored and retrieved from the Quanterra data loggers (Perin et al., 1998), providing more robust data recovery from onsite disks following telemetry outages.
Data from DIAB and MONB in the Bay Area, and 13 stations in the Parkfield regional (all but PKDB), are now being collected at 1 second intervals. Collecting at such high-frequency (for GPS) allows dynamic displacements due to large earthquakes to be better measured, such as was demonstrated by several studies following the 2002 Denali fault earthquake. However, this 30-fold increase in data can be limited by telemetry bandwidth issues. Data from the Parkfield stations are collected on an on-site computer, written to removable disk once per month, and sent to SOPAC for long-term archiving (decimated 30-sec data is acquired daily via the BSL frame relay circuit). In the Bay Area, we have converted the two stations that have sufficient bandwidth and are not collocated with seismic instrumentation. We are currently assessing bandwidth issues at other stations and are planning to convert to 1-second sampling where possible, such as the Mini-PBO stations, in the next year.
The BSL acquired 7 Ashtech MicroZ-CGRS (uZ) receivers with NSF funding for the Mini-PBO project. These receivers, designed for continuous station applications, use less power (5.6 W) than the Z-12 receivers due to the lack of an interactive screen, provide better remote receiver control, and can support serial telemetry in both native raw format and the receiver independent BINEX format. We installed a uZ at SBRN and replaced the Z-12 at OHLN with a uZ in May 2003 after the clock chip on Z-12 at the site began to malfunction. We are currently considering switching to the more compact BINEX format where possible, as this will reduce some of the bandwidth limitations and allow us to convert more stations to 1-second sampling.
The BSL also acquired several Wi-Lan VIP 110-24 VINES ethernet bridge radios. These 2.4 GHz spread spectrum radios use a tree structure to create a distributed ethernet backbone with speeds up to 11 Mbps. Each system uses a directional antenna to talk to its ``parent" in the tree, and an omni-directional antenna to talk to its children, if multiple, or a directional antenna if it has only 1 child. These radios offer several advantages over the Freewave radios used at other sites, including TCP/IP ethernet control, higher bandwidth, and greater flexibility for setting up networks. We installed a set of Wi-Lan radios at the SVIN Mini-PBO station to transmit data from the site to the frame relay circuit, and are assisting EBMUD in converting their continuous station to real-time telemetry using Wi-Lan radios.
In February 2003, telemetry flow of GPS data stopped at MUSB. Access to the site was initially limited by the winter snowpack, and then by the need to coordinate the visit with Southern California Edison engineers. During a visit to the site in August 2003 continuity tests revealed that the hardline antenna cable had apparently failed. This 70 m cable may have been damaged by repeated water freezing in the PVC conduit that houses it. We intend to replace the antenna cable and improve the drainage of the conduit in September 2003.
As part of the Plate Boundary Deformation project (Chapter 8), nine new continuous GPS sites were installed in the Parkfield area (see Table 7.2) in Summer 2001 by the USGS and SIO. These sites span about 25 km on either side of the San Andreas fault and are designed to link the BARD network in central and northern California to the SCIGN network in southern California. As part of this upgrade, the new station CARH was installed at Carr Hill near the original CARR station, which had been running since 1989. After allowing both stations to run side-by-side for nearly 2 years, CARR was turned off in April 2003. In February 2003, the NCEDC assumed responsibility from the USGS Pasadena for the telemetry download of these stations over their existing frame relay circuit at Parkfield. We installed an onsite LINUX computer that controlled the sequential download of data. In June 2003, the stations were upgraded to real-time streaming using Wi-Lan radios by SIO and the USGS.
Throughout the year, we have continued installations for the NSF-funded mini-PBO project establishing collocated GPS, and borehole strainmeter and seismometer observatories in the Bay Area. Completion of these sites have been hampered by problems with the original permitting at Ox Mt and Marin Headlands, and by the need to redesign the GPS antenna mount. The GPS system at the Mini-PBO site SBRN was installed in early March 2003. Two major changes were made since the installation of the first Mini-PBO site at OHLN. We increased the diameter of the upper part of the steel shroud, which protects the borehole casing, from 10 to 14 inches to ensure that the casing remains decoupled from the surface. We also used a new borehole adapter for the GPS mount that was machined from two stainless steel flanges. Should borehole access be needed, this adapter allows a very high level of horizontal accuracy when reinstalling the antenna. These GPS mounts will be installed at the remaining 3 sites in Fall 2003 after rainfall lessens the fire hazards posed by the required welding of the lower flange onto the casing. For more details about the Mini-PBO station installations, see Chapter 8).
We assisted Hat Creek Radio Observatory (HCRO), located in northeastern California near Mt Lassen (Figure 7.1), in designing and installing a continuous GPS station. The HCRO is installing the new Allen Telescope Array (ATA), which will consist of approximately 350 6.1-meter radio telescope dishes arrayed at the site, for both astrophysical and Search for Extraterrestrial Intelligence (SETI) studies. We previously assisted UC Berkeley astrophysicists in conducting an RTK survey of the HCRO site to determine the optimal locations for the 350 dishes using Trimble RTK equipment purchased for the project. After completion of the RTK survey, the base receiver was converted into the continuous station. The site is set amidst and underlain by extensive lava fields. After extensive reconnaissance of the site, we chose a monument location that is close to the main laboratory buildings, unlikely to be affected by future ATA dish placement, and on a reasonably stable lava flow. In June 2003, we assisted with the construction of a 12"-diameter concrete pier that is anchored to the lava flow outcrop. The Trimble Zephyr antenna was attached using a SCIGN adapter. We are currently establishing data acquisition procedures with the HCRO to archive the data at the NCEDC.
We also assisted Thales, Inc. (formerly Ashtech, Inc.) to establish a continuous GPS station on the roof of their Santa Clara office building. The chokering antenna is attached to a metal pin that was drilled and cemented into a corner of the roof's concrete parapet. Data is currently acquired daily by FTP from a server located at Thales, and we are investigating methods to acquire the data more rapidly using some of the TCP/IP capabilities of the recently developed Ashtech iCGRS receiver. Other agencies have also installed new continuous stations in the Bay Area, including an FAA site in Fremont (ZOA1) that will be used for Wide Area Augmentation System (WAAS) navigation control.
The BSL staff is evaluating the performance of the UNAVCO-designed L1 system in an urban setting. This single-frequency receiver is relatively inexpensive but is less accurate than dual-frequency receiver systems that can completely eliminate first-order ionospheric effects. Hence we expect the L1 system to be most useful for short baseline measurements where ionospheric effects tend to cancel due to similar propagation paths. The systems are self-contained, using solar power and integrated radio modems.
In April 2002, we installed 4 sites in a 10-km profile extending normal to the Hayward fault between the UC Berkeley campus and BARD station BRIB (Figure 7.2). Due to the topography of the East Bay hills, each site acts as a repeater for other sites. Data from WLDC passes through all the other stations, with its relay path being (in order) BDAM, VOLM, GRIZ, a repeater on the UC Berkeley Space Sciences Building, and then finally the master radio on the roof of McCone Hall where the BSL is located on campus. This profile, complemented by BRIB and EBMD to the west of the fault, will be most sensitive to variations in locking at 2-8 km depth. We expect that these systems will provide useful constraints on relative displacements near the Hayward fault in 3-5 years, and should help to resolve variations in creeping and locked portions of the fault (e.g., Bürgmann et al, 2000).
Between April 2002 and January 2003, the L1 system operated reasonably well, although problems with faulty batteries solar power regulators caused some loss of data. The Freewave radio at the repeater site SPSC was replaced with an Intuicom system. The original radio was sent in for routine maintenance and was found to have a frequency crystal that was beyond its normal operating range. In mid-January 2003, the solar panel at GRIZ was stolen, which resulted in damage to the cables located outside of the protective metal enclosure. The replacement solar panel was installed in a steel channel frame welded to the vertical steel post that forms the monument base. A 0.5"-thick Plexiglas layer was inserted to protect the surface of the solar panel. Acquisition of all data failed not long after this repair. Initial tests suggested a problem at the repeater site SPSC, but subsequent efforts failed to resolve the problem. In August 2003 we isolated the problem to bad cable connections at the GRIZ sites and re-established operations of the network.
We are processing the data using the GAMIT/GLOBK analysis package, which required modifications to handle L1-only observations. We corrected software provided by UNAVCO to synchronize the phase, pseudorange, and clock offset observables, which allows the data to be cleaned in an automatic fashion. Preliminary results suggest that repeatabilities of 1-2 mm in daily horizontal relative positions and 5 mm in the vertical on the shortest (several km) baselines can be achieved (Figure 7.3), but these degrade to 3-4 mm on the longer (10 km) baselines. We are investigating ways to simultaneously process the dual-frequency data from nearby BARD stations (e.g., BRIB, OHLN), with the single-frequency L1 data to improve these results. Currently data from second frequency on the BARD stations is not used, which degrades the definition of the local reference frame and repeatability of the baselines.
We use the GAMIT/GLOBK software developed at MIT and SIO to process data from the BARD and other nearby continuous GPS networks. We have recently modified our processing strategies to take better advantage of recent enhancements to the GAMIT software and automated scripts. These improvements include better accounting of ocean-tide effects, estimating gradients in atmospheric variations, and applying elevation-dependent weighting to the data observables. We process data from more than 70 stations within hours of the completion of the day using rapid or predicted orbits and are reprocessing older data from the present to 1991 using improved orbits, which we expect to be completed by Fall 2003. Data from 5 primary IGS fiducial sites located in North America and Hawaii are included in the solutions to help define a global reference frame. For long-term velocity estimates, we combine these solutions with global and regional solutions provided by SOPAC to better define a stable North America reference frame.
The estimated relative baseline determinations typically have 2-4 mm long-term scatter in the horizontal components (Figure 7.4) and the 10-20 mm scatter in the vertical. Average velocities for the longest running BARD stations during 1993-2003 are shown in Figure 7.5, with 95% confidence regions assuming only white noise. The velocities are relative to stable North America, as defined by the IGS and CORS fiducial stations. Together with students in the department who are now using the GAMIT software to process survey-mode data in the San Francisco Bay area, we are working to combine the survey-mode and continuous GPS solutions into a self-consistent velocity field for northern California.
Most of the Sierra Nevada sites (CMBB, QUIN, and ORVB), as well as SUTB in the Central Valley, show little relative motion, indicating that the northern Sierra Nevada-Central Valley is tectonically stable. The motion of these sites relative to North America differs from the inferred motion of the western Basin and Range Province, suggesting 3 mm/yr right-lateral shear across the Walker Lane-Mt. Shasta seismicity trend. Deformation along the coast in central California is dominated by the active SAF system, which accommodates about 35 mm/yr of right-lateral shear. The Farallon Island site (FARB) off the coast of San Francisco is moving at nearly the rate predicted by the NUVEL-1A Pacific-North America Euler pole. Two-dimensional modeling of the observed fault-parallel strain accumulation (Murray and Segall, 2001) predicts deep slip rates for the San Andreas, Hayward, and Calaveras/Concord faults are 19.31.8, 11.31.9, and 7.41.6 mm/yr, respectively, in good agreement with estimated geologic rates (174, 92, and 53 mm/yr, respectively). Most of the 46 mm/yr of relative motion is accommodated within a 100-wide zone centered on the SAF system and a broader zone in the Basin and Range Province in Nevada.
We are developing real-time analysis techniques that will enable rapid determinations (within minutes) of deformation following major earthquakes to complement seismological information. We use GAMIT/GLOBK processing techniques to estimate independent hourly solutions at the several cm-level horizontal precision and during the past year established an extension of the REDI system where estimates of postseismic positions are attempted when 10 minutes of data become available following an earthquake (Murray et al., 2002).
We currently process 1 hour data batches available within 20 minutes of measurement from more than 20 continuously telemetered BSL and other stations providing hourly data. The hourly solutions have higher scatter than the 24-hour solutions (3-10 mm in the horizontal and 10-30 mm in the vertical), but our simulations suggest that displacements 3-5 times these levels should be reliably detected, and that the current network should be able to resolve the finite dimensions and slip magnitude of a M7 earthquake on the Hayward fault. Due to the poor ability of GAMIT to resolve ambiguities from short data spans, estimates of coseismic displacements within minutes of an event have high (decimeter-level) uncertainty. We are testing a relatively new component of GAMIT that uses Kalman filtering techniques and improved ambiguity resolution methods to provide higher-precision kinematic positions. This method works well for networks with small interstation distances (e.g., near the 1999 Hector Mine earthquake), which aids ambiguity resolution, but has less success on more widely spaced networks, such as the continuous GPS stations in the vicinity of the 2002 Denali earthquake. The August 1998 M=5.1 San Juan Bautista earthquake (Uhrhammer et al., 1999) is the only event to have produced a detectable earthquake displacement signal (of 4 mm) at a BARD GPS receiver.
Mark Murray oversees the BARD program. André Basset, Bill Karavas, John Friday, Dave Rapkin, and Doug Neuhauser contribute to the operation of the BARD and L1 networks. Mark Murray and André Basset contributed to the preparation of this chapter.
Bürgmann, R., D. Schmidt, R.M. Nadeau, M. d'Alessio, E. Fielding, D. Manaker, T. V. McEvilly, and M. H. Murray, Earthquake potential along the northern Hayward fault, California, Science, 289, 1178-1182, 2000.
King, N. E., J. L. Svarc, E. B. Fogleman, W. K. Gross, K. W. Clark, G. D. Hamilton, C. H. Stiffler, and J. M. Sutton, Continuous GPS observations across the Hayward fault, California, 1991-1994, J. Geophys. Res., 100, 20,271-20,283, 1995.
Murray, M. H., and P. Segall, Continuous GPS measurement of Pacific-North America plate boundary deformation in northern California and Nevada, Geophys. Res. Lett., 28, 4315-4318, 2001.
Murray, M. H., R. Bürgmann, W. H. Prescott, B. Romanowicz, S. Schwartz, P. Segall, and E. Silver, The Bay Area Regional Deformation (BARD) permanent GPS network in northern California, EOS Trans. AGU, 79(45), Fall Meeting Suppl., F206, 1998.
Murray, M., Neuhauser D., Gee, L., Dreger, D., Basset, A., and Romanowicz, B., Combining real-time seismic and geodetic data to improve rapid earthquake information , EOS. Trans. AGU, 83(47), G52A-0957, 2002.
Perin, B. J., C. M. Meertens, D. S. Neuhauser, D. R. Baxter, M. H. Murray, and R. Butler, Institutional collaborations for joint seismic and GPS measurements, Seismol. Res. Lett., 69, 159, 1998.
Romanowicz, B., B. Bogaert, D. Neuhauser, and D. Oppenheimer, Accessing northern California earthquake data via Internet, EOS Trans. AGU, 75, 257-260, 1994.
Uhrhammer, R., L. S. Gee, M. Murray, D. Dreger, and B. Romanowicz, The 5.1 San Juan Bautista, California earthquake of 12 August 1998, Seismol. Res. Lett., 70, 10-18, 1999.
Berkeley Seismological Laboratory
215 McCone Hall, UC Berkeley, Berkeley, CA 94720-4760
Questions or comments? Send e-mail: email@example.com
© 2004, The Regents of the University of California