The Bay Area Regional Deformation (BARD) network of permanent, continuously operating Global Positioning System (GPS) receivers monitors crustal deformation in the San Francisco Bay area ("Bay Area") and northern California (Murray et al., 1998a). It is a cooperative effort of the Berkeley Seismological Laboratory at UC Berkeley (BSL), the U.S. Geological Survey (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 45 permanent stations and will expand to about 55 stations by July 2001 (Figure 6.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 BARD Consortium and with surveying, meteorological, and other interested communities.
BARD presently includes 45 permanent, continuously operating stations, 8 of which monitor the Long Valley caldera near Mammoth. The remaining 37 stations (Table 6.1) include 19 maintained by the BSL (including 3 with equipment provided by Lawrence Livermore National Laboratory (LLNL), the Omnistar Corporation, and UC Santa Cruz), 6 by the USGS, 2 by Trimble Navigation, and one each by LLNL, Stanford University, UC Davis, and East Bay Municipal Utilities District. Other stations are maintained by institutions outside of northern California, such as the National Geodetic Survey, the Jet Propulsion Laboratory, and the Scripps Institution of Oceanography, as part of larger networks devoted to real-time navigation, orbit determination, and crustal deformation.
Between 1993 and 1996, the BSL acquired 5 Ashtech Z-12 receivers from UC Berkeley and private (EPRI) funding, which together with 2 USGS receivers, formed the nucleus of the initial BARD network. Since 1996, the BSL has acquired additional Ashtech Z-12 receivers with Dorne-Margolin design choke ring antennas: 13 in 1996 from a combination of federal (NSF), state (CLC), and private (EPRI) funding, and 4 in the last year from federal (USGS) funding. All but 5 of these receivers have been permanently installed to enhance continuous strain measurements in the Bay Area and to consolidate the regional geodetic network. The network 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 6.1). Three more of the BSL receivers will be installed next year, and 7 additional receivers will installed during 2000-2002 as part of the NSF-funded mini-PBO project, which will establish collocated GPS/seismometer/borehole strainmeter observatories in the Bay Area (see Chapter 7).
In 1996, researchers from the BSL, the USGS, Stanford University, LLNL, UC Davis, and UC Santa Cruz formed a consortium of institutions involved in studies of the tectonic deformation in the San Francisco Bay area and northern California. Members of the BARD consortium agreed to pool existing resources and coordinate development of new ones in order to advance an integrated strategy for improving the temporal and spatial resolution of the strain field. They agreed in principle to the continued development of the network of permanently deployed GPS receivers, to the development and maintenance of a pool of GPS receivers for survey-mode operations that may be deployed in semi-permanent mode in the Bay Area when not otherwise in use, to archiving of all data at the NCEDC, and to the development of a coordinated data analysis facility that will process permanent, semi-permanent, and survey data.
In the remainder of this section, we describe BARD-related activities the BSL has performed over the last year, including the installation of a new station, remonumentation and other upgrades to existing stations, continued assessment of an experimental single-frequency receiver, improvement in data archive and processing methods, and analysis of the data to estimate deformation signals monitored by the network.
During July 1999-June 2000, one new BARD station (MODB) was installed near Alturas in Modoc county in the northeast corner of California. It is collocated with the new broadband seismic station MOD. This station should provide important kinematic constraints on the poorly understood tectonics of the region, which is roughly where the Basin and Range extensional province in Nevada transitions into the backarc of the Cascadia subduction zone in southern Oregon. This region may also accommodate deformation associated with the eastern California shear zone as it extends north from Owens Valley.
The design of MODB GPS station is similar to other BSL permanent stations. It uses a low-multipath choke-ring antenna mounted to a reinforced concrete pillar approximately 0.5 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. It uses a low-loss antenna cable 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. The Ashtech Z-12 receiver is programmed to record data once every 30 seconds, observing up to 12 satellites simultaneously at elevations down to the horizon.
Unlike most other stations, which use a combination of radio modem or frame-relay technologies to provide continuous telemetry, MODB uses the collocated Quanterra seismic datalogger to package the GPS data into blockettes, which are then sent using VSAT satellite telemetry to Berkeley via Golden, Colorado. Initial problems with this telemetry method, primarily caused by long forwarding times of the data packets, have recently been resolved and we are preparing to make the data publicly available.
New monuments were installed at the existing stations BRIB and CMBB during early 2000. These sites were the first GPS installations performed by the BSL staff, and the antennas were mounted directly to bedrock near the ground surface, which was considered an optimal configuration at that time (late 1993). However, subsequent tests performed by UNAVCO on low antenna mounts revealed that estimates of tropospheric water vapor from the GPS data are strongly correlated with signal multipath errors, which can degrade the precision of the vertical position estimates. The new concrete pier monuments at BRIB and CMBB, similar to most of the other BSL GPS stations, elevate the antennas 0.5-1.0 m above the ground surface, which helps to minimize the correlations between multipath and tropospheric parameters. At the same time, a low-multipath choke-ring antenna was installed at BRIB, and SCIGN-designed adapters and protective domes were installed at both stations to improve their security and make them more consistent with the other BSL GPS stations that were upgraded during 1999.
The BSL staff helped to improve data retrieval from the UCD1 antenna located on the roof of the UC Davis Geology Building. Data had previously been downloaded manually at weekly intervals. We helped to configure a new telemetry path using a serial connection to a computer where data is continuously made available via ftp. We also helped UC Santa Cruz install the same system for a new BARD station (UCSC) that was established on the roof of the Earth Sciences Building on the UC Santa Cruz campus. This station uses a permanently installed Trimble antenna, but will operate in semi-permanent mode depending on receiver availability, primarily when they are not in use on survey-mode projects.
Significant disruptions were experienced at several stations. The NUNE station ceased operation in December 1999 due to water damage sustained during a winter storm. This station has repeatedly experienced problems due to landsliding of the near surface during severe storms that displace the receiver and antenna enclosure relative to the more deeply anchored antenna. We have decided to abandon this site in favor of a new site on nearby Pt. Wilson to be installed as part of the mini-PBO project.
The MOLA station, located at a former U.S. Navy fuel supply depot, has experienced significant power and telemetry disruption during the last 2 years as the site undergoes hazard waste cleanup. We established radio telemetry to the station via Yerba Buena Island to replace the phone line that had been became permanently inoperable in 1998, and installed a temporary power line to replace the previous line that been cut in mid-2000. We are also considering moving this station if a suitable site on nearby Pt. San Pablo can be established as part of the mini-PBO project.
The BSL staff completed the antenna dome upgrade project begun last year. We purchased SCIGN-designed hemispherical domes using federal (USGS) funding. Domes cover the antennas to provide security and protection from the weather and other natural phenomenon. Previously the BSL stations had a mixture of dome types or none at all, adding a potential non-uniformity to signal delays and antenna phase patterns. The new SCIGN dome is designed for the Dorne-Margolin antennas and minimizes differential radio propagation delays by being hemispherical about the phase center and uniform in thickness at the 0.1 mm level. It is also very resistant to damage and, in its tall form in combination with the SCIGN-designed antenna adapter, can completely cover the dome and cable connections for added protection. Tall domes and adapters were installed at the new and upgraded stations, as well as 6 stations that require the most security due to public accessibility: BRIB, CMBB, DIAB, MODB, MUSB, PTRB, SAOB, TIBB, and YBHB. Short domes were installed at 8 less accessible stations: HOPB, LUTZ, MHCB, MONB, ORVB, PKDB, POTB, and SODB.
Data from all BSL-maintained stations are collected at 30-second intervals and transmitted continuously over serial connections. Two sites (MOLA and TIBB) have direct radio links to Berkeley, and MODB uses VSAT satellite telemetry. The remaining 17 sites use frame relay technology, either alone or in combination with radio telemetry. Twelve GPS stations are collocated with broadband seismometers and Quanterra data collectors. With the support of IRIS we have developed software that will allow continuous GPS data to be converted to MiniSEED opaque blockettes, which can be stored and retrieved from the Quanterra dataloggers (Perin et al., 1998). This approach preserves GPS data during telemetry outages. During the past year, we installed a new version of the Quanterra system software that recognizes the new MiniSEED blockettes, established generally reliable connections at 10 of the 12 collocated sites, and worked to resolve problems at the remaining 2 stations (BRIB and POTB). We are currently developing procedures to preferentially use the MiniSEED data where available, based on our experiences with data from MODB, which is only available using that method.
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 (NCEDC) 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/).
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.
Many of the BARD sites are classified as CORS stations by the NGS, which are used as reference stations by the surveying community. All continuous stations operating in July 1998 were included in a statewide adjustment of WGS84 coordinates for this purpose; a more recent adjustment is currently underway. Members of the BARD project regularly discuss these and other common issues with the surveying community at meetings of the Northern California GPS Users Group.
In the past 2 years the BARD Project and the NCEDC have collaborated with UNAVCO and other members of the GPS community to define database schema and file formats for the GPS Seamless Archive Centers (GSAC) project. When completed this project will allow a user to access the most current version of GPS data and metadata from distributed GSAC locations. The NCEDC will participate at several levels in the GSAC project: as a primary provider of data collected from BSL-maintained stations, as a wholesale collection point for other data collected in northern California, and as a retail provider for the global distribution of all data archived within the GSAC system. We have produced monumentation files describing the data sets that are produced by the BARD project or archived at the NCEDC, and are during the last year began creating incremental files describing changes to the holdings of the NCEDC so that other members of the GSAC community can provide up-to-date information about our holdings.
The data from the BARD sites generally are of high quality and measure relative horizontal positions at the 2-4 mm level. The 24-hour RINEX data files are processed daily with an automated system using high-precision IGS orbits. Final IGS orbits, available within 7-10 days of the end of a GPS week, are used for final solutions. Preliminary solutions for network integrity checks and rapid fault monitoring are also estimated from Predicted IGS orbits (available on the same day) and from Rapid IGS orbits (available within 1 day). 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. Average station coordinates are estimated from 24 hours of observations using the GAMIT software developed at MIT and SIO, and the solutions are output with weakly constrained station coordinates and satellite state vectors.
Processing of data from the BARD and other nearby networks is split into 7 geographical subregions: the Bay Area, northern California, Long Valley caldera, Parkfield, southern and northern Pacific Northwest, and the Basin and Range Province. Each subnet includes the 5 IGS stations and 3 stations in common with another subnet to help tie the subnets together. The weakly constrained solutions are combined using the GLOBK software developed at MIT, which uses Kalman filter techniques and allows tight constraints to be imposed a posteriori. This helps to ensure a self-consistent reference frame for the final combined solution. The subnet solutions for each day are combined assuming a common orbit to estimate weakly constrained coordinate-only solutions. These daily coordinate-only solutions are then combined with tight coordinate constraints to estimate day-to-day coordinate repeatabilities, temporal variations, and site velocities.
The estimated relative baseline determinations typically have 2-4 mm WRMS scatter about a linear fit to changes in north and east components and the 10-20 mm WRMS scatter in the vertical component. Average velocities for the longest running BARD stations during 1993-2000 are shown in Figure 6.2, with 95% confidence regions. We have allowed randow-walk variations in the site positions in order to more accurate characterization of the long-term stability of the site monuments and day-to-day correlations in position. The velocities are relative to stable North America, as defined by the IGS fiducial stations, which we assume have relative motions given by Kogan et al., (2000).
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 in the Pacific Northwest is generally consistent with interseismic strain accumulation along the Cascadia megathrust, the interface between the Juan de Fuca and North America plates, particularly in Washington where the velocity vectors are nearly parallel to the oblique convergence direction. Greater arc-parallel motion in Oregon and northern California may be due to the influence of the SAF system to the south and clockwise rotation of the southern Oregon forearc ( Savage et al., 2000).
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 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. Additional details about the deformation results and modeling are include in a later chapter, entitled "Continuous GPS Measurement of Pacific-North America Plate Boundary Deformation in Northern California".
These BARD results are being combined with older VLBI and spatially dense Geodolite EDM and survey GPS measurements collected by the USGS, Stanford University, and UC Berkeley (Murray et al., 1998b). The combined velocity map will provide significantly improved constraints on three-dimensional locking depth and deep-slip models of strain accumulation, which will be used for seismic hazard assessment along the SAF system.
We are also developing real-time analysis techniques that will enable rapid determinations (minutes) of deformation following major earthquakes to complement seismological information and aid determinations of earthquake location, magnitude, geometry, and strong motion (Murray et al., 1998c). We currently process data available within 1 hour of measurement from the 18 continuous telemetry BSL stations, and several other stations that make their data available on an hourly basis. The data are binned into 1 hour files and processed simultaneously. The scatter of these hourly solutions is much higher than the 24-hour solutions: 10 mm in the horizontal and 30-50 mm in the vertical. 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 M=7 earthquake on the Hayward fault. We are currently investigating other analysis techniques that should improve upon these results, such as using a Kalman filter that can combine the most recent data with previous data in near real-time. 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 at a BARD GPS receiver.
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 polars and an integrated radio modems. During 1999, the BSL borrowed 2 receivers and a master radio from UNAVCO to perform the evaluation, but persistent hardware and software problems limited progress on this project. UNAVCO subsequently resolved many of the problems and in summer 2000, we received new, improved equipment and software for 4 systems and a master radio. We are currently working closely with UNAVCO engineers to finishing testing and securing land permits in order to deploy the systems on a 10-km profile extending normal to the Hayward fault between the UC Berkeley campus and the permanent BRIB site.
We have also resolved processing problems reported in the previous Annual Report. These problems were caused by inappropriate input parameters to the automatic editing routines, which significantly degraded ambiguity resolution. We currently achieve 2-4 mm repeatability in horizontal components on 1.5-km long baselines, and 3-5 mm repeatability on 10-km long baselines. These results are similar to those obtained using Bernese processing by UNAVCO, and may be improved by incorporating ionospheric models using data from the nearby dual-frequency permanent receivers. 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).
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.
Kogan, M. G., G. M. Steblov, R. W. King, T. A. Herring, D. I. Frolov, S. G. Egorov, V. Y. Levin, A. Lerner-Lam, A. Jones, Geodetic constraints on the rigidity and relative motion of Eurasia and North America, Geophys. Res. Lett., 27, 2041-2044, 2000.
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, 1998a.
Murray, M. H., W. H. Prescott, R. Bürgmann, J. T. Freymueller, P. Segall, J. Svarc, S. D. J. Williams, M. Lisowski, and B. Romanowicz, The deformation field of the Pacific-North America plate boundary zone in northern California from geodetic data, 1973-1989, EOS Trans. AGU, 79(45), Fall Meeting Suppl., F192, 1998b.
Murray, M. H., D. S. Dreger, D. S. Neuhauser, D. R. Baxter, L. S. Gee, and B. Romanowicz, Real-time earthquake geodesy, Seismol. Res. Lett., 69, 145, 1998c.
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.
Savage, J. C., J. L. Svarc, W. H. Prescott, and M. H. Murray, Deformation across the forearc of the Cascadia subduction zone at Cape Blanco, Oregon, J. Geophys. Res., 105, 3095-3102, 2000.
Uhrhammer, R., L. S. Gee, M. Murray, D. Dreger, and B. Romanowicz, The Mw 5.1 San Juan Bautista, California earthquake of 12 August 1998, Seismol. Res. Lett., 70, 10-18, 1999.