Institute for Crustal Studies
2000-2001 Annual Report

Next Page

Previous Page

Table of Contents

Awards Administered

Crustal Materials


Frank Spera

12/15/97 - 12/14/01

Department of Energy, DE-FG03-91ER14211

7/15/99 — 6/30/02

National Science Foundation, EAR - 9972922

Magma Rheology, Mixing of Rheological Fluids, Molecular Dynamics Simulation
and Lithospheric Dynamics

Experimental Rheology of Magmatic Emulsions and Foams: Effects of
Composition, Temperature, Shear Rate and Bubble Content

Construction of a high temperature rheometer based on a concentric cylinder design has been completed. This device can investigate the rheology of magma, a complex mixture of solids, melt and vapor bubbles, at temperatures from 700 °C to 1375 °C for shear rates in the range 10-5 to 1 s-1. The instrument has been calibrated and certified using NBS standard glasses and is accurate to within 0.02 log10 units in viscosity throughout the temperature and shear rate range. This 2-sigma uncertainty is better than the NBS quoted value of 0.029 log10 units. A paper was published in Review of Scientific Instruments where details may be found. Experiments are presently underway exploring the relative viscosity, defined as the ratio of the viscosity of the mixture to the viscosity of single phase melt for a rhyolite from the DOE sponsored drilling site at Obsidian Dome, near Long Valley California as well as melt-vapor emulsions of molten albite and molten orthoclase compositions. Results suggest that for shear rates exceeding a critical value, the relative viscosity decreases with increasing bubble content whereas for shear rates below the critical threshold, the relative viscosity increases as the bubble content increases in accord with the classical result derived by G.I. Taylor 60 years ago. The goal of this work is to develop comprehensive non-linear constitutive relations for melt-vapor emulsions in the range of shear rate and temperature relevant to crustal magmas. A short paper has been submitted to Earth and Planetary Science Letters detailing some of these results and presently another paper is being prepared for submission to Journal of Volcanological and Geothermal Research. Other work completed this year includes an extensive set of MD simulations applicable to molten and glassy CaAl2Si2O8 at temperatures and pressures characteristic of continental lithosphere, study of the glass transition in glassy anorthite and MD study of melts in the system NaAlO2-SiO2. This work was published in American Mineralogist. Other MD simulations for compositions in the system Na2O-SiO2-H2O are also underway. A final sub-project concerns porous media thermohaline convection. We have confirmed that in low porosity geologic media, a doubly-advective instability sets in that leads to chaotic behavior for Darcy Rayleigh numbers and buoyancy ratios in the range found in nature. These simulations take account of the temperature and composition dependence of fluid properties (H2O +NaCl) and allow for the importance of dispersion (rather than simple molecular diffusion) of solute. This doubly-advective instability has application to diagenesis, the formation of crustal ore deposits and the evolution of metamorphic terrains. This work is in press at Earth and Planetary Science Letters.

Frank Spera

2/1/01 — 1/31/03

National Science Foundation, Award EAR-0073932

Engergy-Constrained Geochemical Models for Open System Magma

Existing tools that model the geochemical evolution of open-system magma bodies suffer from several limitations, the most critical of which are lack of accounting of energy-conservation and lack of inclusion of the consequences of wallrock partial melting. A new tool has been developed that models the geochemical evolution of a magma body-wallrock system in which energy, mass and species are conserved as wallrock and magma approach thermal equilibrium. We are attempting to build upon this EC-AFC framework by accomplishing five new tasks. (1) Implementation of a self-consistent treatment of magma Replenishment to the existing model to develop a user-friendly EC-RAFC model. (2) More tightly couple phase equilibria and thermochemical constraints to EC-RAFC. This is accomplished using information from the MELTS code (e.g., Ghiorso, 1997) to implement a procedure whereby melt productivity functions fa(Ta) and fm(Tm), which are essential in EC-RAFC, are determined for user-defined compositions of initial magma, wallrock and replenishment magma. In addition, phase equilibria data will enable better approximations of bulk partition coefficients for magma and wallrock as a system undergoes RAFC. (3) Allow for the development and maintenance of compositionally distinct melt reservoirs within an EC-RAFC magma body. This will enable us to apply EC-RAFC to zoned ignimbrites. The prevalence of such deposits in the geological record (e.g., in continental hydrothermal-magmatic regions) and their associations with zoned magma bodies make this an essential step for developing an EC-RAFC model that is broadly applicable. (4) Allow for a melt extraction efficiency less than unity (i.e., where 100% of wallrock anatectic melt is not added to the standing magma body); this modification accommodates a scenario where mass and energy addition can be decoupled and addresses the physics of melt extraction in a quantitative manner. (5) Apply the enhanced version of the EC-RAFC simulator to six to a dozen carefully documented magmatic systems spanning a range of geologic environments in order to explore potential systematics in the relative roles of fractional crystallization, recharge and assimilation in magma bodies associated with particular tectono-magmatic settings (e.g., MOR, continental arcs). This research has the potential to impact the work of the large number of geoscientists who generate and interpret geochemical, geochronological, petrographic and petrologic data from extant and ancient magmatic systems and will parallel the dramatic improvement in the amount and quality of geochemical data describing magmatic systems.

Crustal Structure and Tectonics

James R. Boles

Jordan F. Clark,

Ira Leifer


University of California Energy Institute, UCSB 20000774

Temporal Variation of a Major Natural gas Seep,
Santa Barbara Channel, California

A large natural methane seep, occurring a kilometer offshore in 64m of water, has been monitored on an hourly basis for 10 months. About 16.4 cubic meters (580 MCF) of methane/day is captured from the seep by two large steel tents (each 30m by 30m) open at the bottom to the ocean. The gas is piped to shore where it is metered and processed. We have analyzed the variation in seep rate as a function of anthropogenic and natural factors. The affect of hydrocarbon production from the nearby South Ellwood oil field on this seep is not clear. There is no clear relation between seep rate and change in formation pressure in key wells over a 10 year period, and from a month long platform shutdown. Results, however, show a clear tidal influence on the gas flow rate; these results are the first documentation of the effect of tides on natural gas seeps.

High tide results in a reduced flow, whereas low tide correlates with an increase in flow rate. Time series analyses clearly shows all of the four principal diurnal tidal components at 12.00 hr, 12.42 hr, 23.93 hr, and 25.82 hr. The correlation indicates that each 30 cm (foot) increase of sea height results in a decrease of 3.0 to 3.8 cubic meter/hour or 0.5% of the hourly flow rate. A simple physical model indicates that the observed changes can best be accounted for by a pore activation model, whereby gas is released from small pores at low pressures but release is inhibited at higher pressure. Our study implies that sea level drops, even those as small as tidal variations, can result in significant increase in degassing from sediment. Major sea level changes might result in substantial degassing of the sediment.


Douglas W. Burbank

4/1/01 — 1/31/04

National Science Foundation, Award EAR-0117242

Scaling and Displacement Relationships for Thrust Faults


Through field studies in New Zealand and Kyrgyzstan, we will quantify the style of deformation generated by the initiation, propagation, and linkage of active thrust faults. A primary goal of the project is to document variations in slip along thrust faults with strike-lengths ranging from a few meters to >10 km. Although similar studies on normal faults have prompted a significant step forward in the perceived understanding of that type of fault, no such population studies exist for thrust faults. We propose to analyze the spectacularly well preserved deformation associated with the Ostler fault on the South Island of New Zealand and with a suite of thrusts that slice through bedrock and are well exposed in the central Tien Shan of Kyrgyzstan. Excellent structural markers (fluvial terraces and an exhumed, low-relief erosion surface) provide the critical reference datum for quantification of displacement-length scaling on thrusts, population statistics of thrust faults, linkage among thrusts, and transfer of displacement between them. In addition we will be able to document rates of thrust propagation at time scales of >103 yr, the timing of linkage among them, and variations in the fabric and character of the fault zone as slip increases or dies along a fault trace. An important focus will be on the magnitude and nature of "off-fault" deformation in order to determine how permanent strain is accumulated in the vicinity of active thrust faults. We will utilize differential GPS surveying, ground-penetrating radar, trenching, and micro- to macroscopic analysis of fault-zone rocks to document patterns of strain at several spatial and temporal scales. It is our expectation that research in these unique field settings will provide a robust statistical foundation for description of the geometric properties of thrust faults. Although generally unavailable until now, such research should underpin a significantly improved understanding of thrust faults.

Douglas W. Burbank


National Science Foundation, Award EAR-0196414

Strain rates, patterns, and partitioning during continental transtension

In recent years, considerable attention has been focused on the modes of accommodation of the stresses that arise due to oblique motion between crustal blocks. Numerous studies of forearcs illustrate strain partitioning in both transpressional and transtensional settings. It is less well known, however, whether similar models relating obliquity to partitioning can be applied to regimes of oblique continental extension outside of forearc settings. Nor has it been well documented how structures evolve when the stress regime changes from one of orthogonal extension to one of highly oblique extension or even pure strike-slip faulting.

We are currently testing several hypotheses regarding the nature of continental transtension and a transition from extension to transtension in the northern Owens Valley, White-Inyo Mountains and Fish Lake Valley (a part of the Eastern California Shear Zone):

  • Hypothesis 1: Strike-slip faults that are active during transtension represent reactivated normal faults formed under a previous extensional regime.
  • Hypothesis 2: When sub-parallel strike-slip faults display compensatory changes in mean slip rates (i.e., when the slip rate on one strike-slip fault decreases, the other increases), motion is effectively transferred between the faults by the development of pull-apart basins and linking faults.
  • Hypothesis 3: Under regimes of transtension, two sets of faults develop where one set shows predominantly strike-slip motion, and the other, predominantly dip slip.)
  • Hypothesis 4: The angle between faults and the regional slip vector controls whether strain is accommodated by strain partitioning on to separate faults or whether oblique-slip faults accommodate most strain.

Field work is being conducted in the White-Inyo Mountains of east-central California and is focussed around the following topics:

  • Geometry, timing, rates, and style of the Pliocene/Miocene structural development with special attention to variations along the range.
  • Rates and style of Plio-Pleistocene deformation along range-bounding faults.
  • Late Quaternary deformation in the study area, including both the sense and the rate of slip along faults and defining other strain accommodation, such as via folding.

NSF post-doc Eric Kirby and graduate student Jessica Jager are collecting data on kinematic slip indicators of fault motions, tilt and displacement of strain markers, provenance indicators and facies, Ar-Ar dating and geochemical correlation of tephra, surveyed and mapped offsets of geomorphic features (terraces, pediments, channels, etc.), and dating of these features (radiocarbon, cosmogenic nuclide, and scarp diffusion modeling) to determine the rate and style of deformation since Miocene times.

This study promises to bring a detailed chronologic perspective to bear on the nature of continental transtension. It will define the style and magnitude of strain partitioning and its relationship to regional strain fields, it will determine some mechanisms whereby slip is transferred between strike-slip faults, it will document the presence or absence of significant temporal variations in the along-strike growth of a mountain range, and it will yield a perspective on rates of extensional and transtensional deformation ranging from Miocene to late Quaternary times.

Douglas W. Burbank


National Aeronautics & Space Administration, Award Nag5-10520

Tectonic-climate interactions in active orogenic belts: Quantification

of dynamic topography with SRTM data

In a recently completed NASA mission, most of the terrestrial earth between 60°N and 60°S was imaged with a tandem radar system aboard the shuttle in order to create a high-resolution digital topographic data base. We propose to utilize the newly acquired topographic data in conjunction with data on the magnitude of erosion, recent patterns of deformation, the extent of modern and past glaciers, and variations in modern climate to address several key questions about the evolution of these ranges. At the present time, little of the shuttle-based data has been released to investigators. As a result, we are basing most of our analysis on older 90-m resolution data and on smaller structures (rather than whole ranges) for which we do have high-resolution data. Based on our previous studies in the NW Himalaya, we concluded that the mean topography was controlled by climatic gradients in glacial snowlines, rather than by the strong variations in rock-uplift rates that we have documented. Other recent results suggest that several rapidly deforming ranges (Southern Alps, Himalaya, Taiwan) have attained a long-term topographic steady state, whereby material influx and exflux are in balance within the range. Other equally active, but perhaps "younger" ranges are still in a pre-steady-state condition. These observations spur the following questions:

  • What progressive topographic changes occur within a rapidly deforming landscape before it attains a steady-state form?
  • To what extent can topographic characteristics be used to differentiate between pre-steady-state and steady-state orogenic belts?
  • Are there fundamental differences in the topography of active mountain belts between glaciated and non-glaciated terrains? If so, how are these differences manifested in pre-steady-state and steady-state mountains?
  • To what extent are fundamental controls on topographic characteristics of mountain belts attributable to climate and glacial erosion, or to spatial and temporal variations in rates of rock uplift rates?

In order to address these questions, we are studying an array of collisional ranges that span from clearly pre-steady state to probable steady state, that encompass regions of heavy glaciation to ones with no glaciation, and that represent contrasting climatic zones. These ranges include the Himalaya, the Tien Shan, the Transverse Ranges of California, the Central Ranges of Taiwan, the Olympic Mountains, the Chugach-St. Elias Range of SE Alaska, and the Southern Alps. In each of these ranges, we are merging data on cooling and erosion histories, present and past distribution of glaciers, histories of deformation, sediment fluxes out of ranges, geodetic data, and both existing and new topographic data. We are also finding it useful to study modest scale (5-10 km) active folds that are propagating laterally and being attacked by erosion as the emerge as positive topographic entities. The along-fold changes in topography provide useful insights on the way in which larger ranges could evolve toward a topographic steady state. Syntheses and comparisons of data from these larger ranges over the next year will permit the development of new models for the temporal growth and decay of deforming mountain belts and will provide insights on the climatic, tectonic, and erosional controls on the topographic evolution of ranges as they interact with past and present surface processes.

Cathy Busby

1/1/98 - 8/31/00

American Chemical Society - Petroleum Research Fund, 32624-AC8

High Resolution Sequence Stratigraphic Framework of a
Tectonically Active Forearc Basin

Mesozoic rocks of the Baja California Peninsula form one of the best-exposed and most areally extensive, well preserved convergent margin complexes in the world. Busby proposes that this convergent margin shows an evolutionary trend that is typical of arc systems facing large ocean basins: a progression from highly extensional through mildly extensional to compressional strain regimes. This proposal focuses on the southernmost end of the Peninsular Ranges forearc basin complex of southern California and northern Baja California, Mexico, the Rosario embayment. Busby suggests that the Peninsular Ranges forearc basin complex formed as forearc strike slip basins under a compressional convergent margin strain regime, where strong coupling promoted the development of strike slip faults during oblique subduction. Strike slip basins are known to be the most tectonically active and complex basin type of all, characteristically showing very rapidly alternating subsidence, uplift and tilting events of great magnitude. This makes them ideal for the kind of detailed "tectono-stratigraphic analysis" that is proposed here. Busby and a graduate student will spend 3 to 4 months in the field mapping the entire basin at a scale of 1:50,000 or better using aerial photographs, satellite imagery and topographic bases. Packages of strata will then be dated using foraminifera, nannofossils, macrofossils and magnetostratigraphic data.

Bradley Hacker


National Science Foundation, EAR-9809840

The Thermal, Petrological and Seismological Structure of Subducting
Oceanic Lithosphere

We compiled and assessed a database of the relevant thermodynamic properties of phases expected in mafic and ultramafic rocks at subsolidus pressures up to 8 GPa. These data are used to calculate mineral physical properties (e.g., density, Vp, Vs) at elevated P and T. Rock physical properties are derived using a Voigt-Reuss-Hill average of mineral physical properties weighted by mineral proportions.

A phase diagram for 5 ultramafic rock compositions was constructed using Holland and Powell's Thermocalc program and natural mineral compositions, and modified in light of recent experiments on hydrous phases at high P. A phase diagram for mafic rocks was constructed from Thermocalc and from petrological studies that reported mineral modes, mineral compositions, and MORB bulk-rock compositions. These phase diagrams permit direct assessment of phase proportions, reaction stoichiometries, H2O content, and, in combination with the rock physical property calculations described above, density, Vp, and Vs of mafic and ultramafic rocks in subduction zones.

Application of these principles to a thermal model results in a cross section detailing the mineralogy, H2O content, density, Vp, and Vs, of a specific subduction zone. These predictions are then compared to models of slab devolatilization related to arc magma generation, models of devolatilization-caused earthquakes, calculations of slab buoyancy, and observations of slab wave speed.

For example, the Philippine Sea plate beneath SW Japan shows seismic velocities of Vp=6.8 and Vs=4.0 km/s in the upper 60 km (Ohkura, 2000). These observations are well fit by calculated velocities for epidote-blueschist facies mafic rocks, the expected equilibrium mineral assemblage, implying that phase transformations are keeping pace with increasing P and T. In contrast, in at least 5 other colder Pacific subduction zones, high-frequency energy that travels up a slab waveguide from intermediate-depth earthquakes is delayed 5-9% relative to long-wavelength energy propagated through the surrounding mantle. This may reflect lawsonite + amphibole stability or partial transformation in the eclogite stability field.

Bradley Hacker

James Mattinson


National Science Foundation, EAR-9814889

Exhumation of Ultrahigh-Pressure Rocks in the
Scandinavian Caledonides

Unraveling the high-temperature exhumation history of Norwegian UHP eclogites is severely hampered by the paucity of Caledonian igneous rocks, which could otherwise provide opportunities for measuring P, T, deformation, and the ages of high-temperature minerals. Each of the few known plutons is therefore viewed with considerable interest. The Sognekollen pluton is particularly noteworthy because it is the only pluton involved with the Nordfjord-Sognfjord Detachment Fault--the premier, large-magnitude exhumation structure in the Western Gneiss Region of Norway. The crustal section in the Sognekollen area is a W-plunging synform that exposes eclogite-bearing orthogneiss overlain by metasedimentary rock, metavolcanic rock, and capping Devonian sedimentary rocks. We calculate equilibration conditions of c. 2.0 GPa and 500°C for eclogites reported by Chauvet et al. (1992). Kyanite-staurolite-chloritoid-garnet-(coesite pseudomorphs?) rocks within the metasedimentary section yield univariant equilibria at comparable pressures and temperatures, implying at least 70 km of exhumation. The Sognekollen pluton is a leucocratic garnet-epidote bearing body with associated felsic dikes that intruded along the metasedimentary/metavolcanic rock contact. Pooled U/Pb isotopic ratios from two of the dikes yield a 1.3 Ga upper intercept and a 437+/-19 Ma lower intercept, which we interpret as the crystallization age of the dikes. Zircons from the pluton itself are U rich and show significant Pb loss, but similarities in rock compositions imply that the pluton and dikes are consanguineous. The metasedimentary host rocks to the pluton underwent regional amphibolite-facies metamorphism at 600-1000 MPa. They document strong down-to-the west extension with a range of outcrop- to crystal-scale structures. Discrete mylonite zones, cataclastic faults, and outcrop-scale brittle faulting indicate that the same deformation regime continued during cooling through greenschist-facies temperatures until c. 395 Ma (hornblende and muscovite 40Ar/39Ar ages from Chauvet and Dallmeyer, 1992).


Bradley Hacker

James Mattinson


Woodhole Oceanographic Institute, Award A100169

Constraints on the Genesis of Continental Crust via Arc Magmatism: Geology, Geochemistry, Structure and Physical Properties of the Talkeetna Arc Crustal Section, South Central Alaska

This research is an interdisciplinary study of the exposed crustal section of the Talkeetna island arc to address the rate and mechanism of continental growth at convergent margins. The Talkeetna arc section, particularly in the Nelchina-Tonsina region, is a crustal section through an accreted, Jurassic subduction-related magmatic arc, from volcanic, volcaniclastic and sedimentary rocks at the top to residual mantle peridotites at the base. Detailed studies of exposed, island arc sections can provide a crucial link between geophysical observations and geochemical studies of volcanic rocks in active arcs. In the same way that ophiolite studies have provided an ideal counterpoint to both marine geophysics and analysis of mid-ocean ridge basalts in developing a complete picture of crustal accretion at oceanic spreading ridges, studies of arc sections will be essential to progress in understanding arc magmatism and crustal genesis over the coming decade.

Our study addresses a key question regarding arc processes: What is the rate and mechanism of continental growth at convergent margins? Arc crust is basaltic. Continental crust is formed mainly in arcs. Continental crust is andesitic. How can these three ideas be reconciled? Clearly, although each of these points is open to question, research funded by this grant will mainly address the first one. However, in our proposed study we will also determine compositional layering of the arc crust, and constrain proposed scenarios for lower crustal delamination during and after arc magmatism. Most importantly, work on the Talkeetna section can provide crucial constraints on the composition and timing of accretion of the middle and lower arc crust, which are generally not exposed in active arcs. Evidence from rare plutonic exposures in the Aleutian and Izu Bonin arcs indicates that different processes or different primary magmas may produce intrusive and extrusive rocks, so that inversion of volcanic rock compositions is not sufficient to determine the bulk crustal composition, nor the nature of primary melts passing from the mantle into the lower crust.

In our proposed study, we will address various aspects of crustal genesis in a practical way, applying a variety of analytical techniques. Our research goals are to (1) make an improved geological map and use a fully two -dimensional approach to better constrain the relative proportions of the different rock types in the Talkeetna arc crustal section; (2) determine the deformation history of lower crustal rocks, identify and date faults along which section may be missing or repeated, , make thermobarometric estimates from mineral compositions, and thereby better constrain how much of the section is missing and how to interpolate the data for the rocks that are exposed; (3) make extensive measurements of physical properties of rock samples for comparison with the growing database on seismic refraction and reflection in arcs; (4) conduct detailed investigations of residual mantle harzburgites and igneous ultramafic rocks just below the base of the gabbroic crust in order to determine the mode of melt transport from the mantle into the crust, and determine the relative importance of ultramafic "cumulates"; (5) conduct major trace and isotopic analysis on a comprehensive suite of samples from all the igneous rocks in the section, together with extensive new geochronology studies, to delineate how many liquid lines of descent are represented in the arc crustal section, and which are most volumetrically important; (6) use petrological and trace element modeling techniques to quantify the possible proportions of different rock types, for comparison with the proportions of rock compositions determined by geologic mapping, providing a further constraint on the bulk composition of the arc crust; and (7) to look at the P-T-t history of the lower crust, in order to constrain Moho temperatures in active arcs, and to provide constraints on the density and viscosity of arc lower crust during and after magmatism.

Jeff Lee

7/1/99 — 5/31/01

National Science Foundation, EAR-9902968

Collaborative Research: Dynamics of Intraplate Fault Systems in the
Northern Basin and Province

This collaborative project is investigating the dynamics of a system of intraplate faults in the northern Basin and Range province using paleoseismic techniques on range front faults in central and eastern Nevada. In light of ongoing research on contemporary strain accumulation across this system of faults using GPS, it is timely to obtain a precise chronology of late Quaternary strain release across the province. The two data sets will provide an unprecedented basis for modeling the dynamics of Basin and Range faulting.

To date our research efforts have on range front faults of the Cortez Mountains and Schell Creek Range. For the Cortez Mountains we have (1) completed detailed surficial mapping and fault scarp profiling of the 60-km long range front fault zone, (2) permitted, opened, and initiated logging of three trenches across an alluvial scarp at the mouth of Fourmile Canyon, (3) sampled various units in the alluvial gravels offset by and overlapping the scarp for 14C and 26Cl geochronology, and (4) sampled a 6-m high bedrock scarp for surface exposure age dating. For the Schell Creek Range we have (1) completed reconnaissance surficial mapping of a 70-km long range bounding fault segment and (2) examined a fluvial scarp that transects the main fault zone appropriate for application of alluvial scarp dating.

Bruce Luyendyk


University of California Energy Institute, UCSB 20000798

A New Device to Sample Marine Oil Slicks (ND-SOS)

Oil and tarballs continually wash up on the shores of Santa Barbara Channel (SBC) including both popular tourist beaches in Santa Barbara and those on the islands in the protected Channel Islands National Marine Sanctuary. The oil arises from both natural and anthropogenic sources and is transported by currents and winds to the beaches. The anthropogenic sources are related to oil production, boat transportation, and run-off of car motor oil. The environmental toll taken by an oil spill results from both the quantity of oil released and the changes it undergoes while being carried to shore. Although it is popularly assumed that oil on SBC beaches results from oil production, research suggests otherwise. However, convincing data linking marine slicks and oil on beaches to natural oil and gas seeps are lacking. Oil slicks evolve in time due to numerous processes. The knowledge of this process is imperfect because the source of the oil is often unknown and sampling has been spotty. Our goal is to provide detailed estimates of the quantity and composition of an oil slick as it varies in time and position. We will adapt a technique developed to sample the sea surface using a rotating glass drum suspended on a small catamaran. By towing the catamaran sampler through a slick a continuous sample record can be obtained. This record can then be later analyzed for oil composition and quantity. The sampler will be tested on a natural marine oil slick within the Coal Oil Point hydrocarbon seep field, with the goal of attempting to quantify the amount of oil released from a single oil seep. The proximity of the seeps to UCSB allows easy access to the sites. This UCEI-funded research will be used as matching funds in a companion project funded by the U.S. Minerals Management Service.


Bruce Luyendyk

8/15/97 - 7/31/02

National Science Foundation, OPP-9615281

Collaborative Research: Air-Ground Study of Tectonics at the Boundary
Between the Eastern Ross Embayment and Western Marie Byrd Land, Antarctica: Basement Geology and Structure, and Influences on West Antarctic Glaciation

During the 1998/99 season SOAR conducted an aerogeophysical survey in western Marie Byrd Land in the region of Edward VII Peninsula and the Ford Ranges. A total of 212 lines were flown amounting to about 36,000 km of gravity, magnetic, radar penetration and laser altimetry data over an area of 130,000 km2. All navigation was by precision differential GPS. During 1999/00 we received most of the data products from Lamont Doherty Earth Observatory (LDEO) and the University of Texas at Austin Institute of Geophysics (UTIG).

Efforts to date have focused on quality control and editing of the data. Processed and edited magnetics and gravity data were received from LDEO in February 2000. Because the technician doing the processing had not anticipated that we would want to verify his work, there was a delay until March and April, 2000, until complete raw data and corrections were obtained. Data quality and processing is on the whole very good to excellent. Our changes to the original processing will probably be limited to adding back data in areas of limited crossover coverage where conservative acceptance criteria employed by LDEO have produced significant data gaps.

Surface-elevation and ice-thickness data were received from UTIG in March and April, 2000, although the surface elevation data set is still not quite complete. Processing of these data involves fewer subjective decisions, so we do not anticipate significant changes to the editing, but gridding is more difficult because the ice and bedrock surfaces lack the inherent smoothness of potential fields that were measured at significant elevation. Current efforts are focused on choosing an optimal gridding strategy for the bedrock surface, which is needed for quantitative interpretation of the magnetic and gravity data.

Paleomagnetic samples taken during the complementary ground program have been processed at UCSB, UC Santa Cruz, and Caltech. In addition, AMS on all samples is complete. To aid in the interpretation of the aerogeophysical data we are compiling magnetic properties and densities of the rock samples. We are also in contact with Christine Smith Siddoway who is mapping brittle structure in the Ford Ranges. These features will guide our interpretation of potential fields’ data.

Very preliminary interpretation of the relation between the bedrock topography and the potential fields has begun. As expected, the gravity field shows high coherence with the bedrock surface over most of the area. The Fosdick Mountains area is one of few places where the magnetic fieldshows similar coherence, suggesting that the range has moderate strength of magnetization and uniform polarity and intensity. A few high-amplitude positive magnetic anomalies occur on the geographic southeast part of the survey area, suggesting young volcanic or shallow intrusive centers under the ice.

Paleomagnetic Characteristic directions from Cretaceous and Devonian intrusive rocks in the southern Ford Ranges are all normally magnetized and in a direction statistically identical to the synthetic mid-Cretaceous paleomagnetic pole for East Antarctica (EANT) of DiVenere, et al. (1994). The southern Ford Ranges pole is distinct from the northern Ford Ranges pole of Luyendyk et al. (1996). There are two possible interpretations: 1) There is an EANT/MBL (Marei Byrd Land) crustal boundary of Late Mesozoic or Cenozoic age located within the Ford Ranges beneath the Crevasse Valley glacier that separates the northern and southern Ford Ranges. Or 2), The relative motion of MBL with respect to EANT since the Cretaceous, as recorded in the northern and southern Ford Ranges, is purely the result of vertical-axis block rotations there and the southern Ford Ranges have not been rotated. The second interpretation is reinforced by the fact that none of the southern Ford Ranges VGPs have been tilted about a horizontal axis whereas tilting is seen at several locations in the northern Ford Ranges (Luyendyk et al., 196). Anisotropy of magnetic susceptibility (AMS) measurements have been made at UCSB and UC Santa Cruz on our samples. The measurements reveal a regional pattern of foliation and lineation that depends on the age of the intrusive rocks, as published by Adams (1987),suggesting the development of fabric during or shortly after initial intrusion.

Bruce Luyendyk


National Science Foundation, OPP-9725876

GPS Measurement of Isostatic Rebound and Tectonics
in Marie Byrd Land, West Antarctica

In the 1998/99 season three autonomous Global Positioning Satellite (GPS) receiver systems were installed in western Marie Byrd Land. The objective is to continually record positions at three bedrock points to monitor crustal deformation and postglacial rebound. Stations were installed in the Rockefeller Mountains on Edward VII Peninsula (MBL1; Mt. Patterson), the Phillips Mountains of the northern Ford Ranges (MBL2; Mt. Carbone), and the central Ford Ranges (MBL3; Clark Mountains). During this past season (99/00) these stations were visited twice. On the first visit in November 1999 it was found that all stations suffered wind damage to the wind generator systems. In addition, the disk drives in the data logging system had failed. Before revisiting the stations in January 2000 we designed a smaller wind generator modeled after a system used by Dr. R. Dibble on Mt. Erebus. We decided to switch from logging data on disc to solid state memory. During the January visit MBL2 and MBL3 were upgraded with flash Databridges and Multi-fan Generators. MBL2 was brought into operation, but we were unable to confirm it was logging due to weather and time constraints. MBL3 was operating fine after upgrading. MBL1 was removed except for the tower and antenna but not upgraded due to weather and time constraints. It is now set up and operating at McMurdo Station. We recovered data for sites MBL1, MBL2, and MBL3 during the 1998-1999 season and some data from the site MBL1 in November 1999, providing a time baseline of just under one year. While the velocity results are inconclusive they do suggest that spreading is occurring across the Ross Embayment. To strengthen our solution we used data available from the GPS stations Mount Coats and Mount Cocks from February and November 1999 (from Dr. C. A. Raymond). Our results show that MCM4 (McMurdo) is moving south 7+- 13 mm/yr, west 15+-24 mm/yr, and up 3+-44 mm/yr relative to MBL1, located in the Rockefeller mountains. The errors are scaled one-sigma errors or scaled by a factor of 1.8 such that the chi-squared per degree-of-freedom is one in the velocity solution. Examined in a different way, the baseline between MCM4 and MBL1 is lengthening at a rate of 17+- 22 mm/yr and moving transverse to each other at a negligible rate of 2+-16 mm/yr. The low vertical rate between the two stations lends confidence that the solution is fairly robust. Additional years of data will improve the velocity estimate while reducing the errors.

Christopher Sorlien


University of Southern California

Southern California Earthquake Center, USC 572726

Rapid Subsidence and South Propagation of the Active Santa Monica
Mountains-Channel Islands Thrust

Funding for this project arrived recently and work on this project is just getting underway. Preliminary work shows that the offshore Santa Monica fault, known as the Dume fault, has a 700 m-high sea floor escarpment west of its intersection with the Palos Verdes fault, and about 2 km of vertical separation of the youngest pre-folding horizon (not yet identified). North-dipping reflections beneath and south of this escarpment are interpreted to be fault-plane reflections. The Malibu Coast fault is identified in the hanging-wall north of the Dume fault, and mapped to be continuous with the Santa Cruz Island fault. The hanging-wall block of the Dume fault shows little net vertical motion, as there is little erosion of the mainland continental shelf southwest of Point Dume. Lack of uplift with respect to sea level is nearly irrelevant to activity on the Dume fault and underlying blind Santa Monica Mountains thrust; it is the structural relief that is important. Vertical motion is the sum of tectonic uplift, sediment compaction (can be ignored on shelf), and isostatic subsidence.



Ralph Archuleta

4/1/91 - 1/31/02

University of Southern California

Southern California Earthquake Center, USC 572726

See individual research summaries below for projects included in this grant.

Ralph Archuleta


University of Southern California

Southern California Earthquake Center, USC 572726

The Portable Broadband Instrument Center (PBIC)

The Portable Broadband Instrument Center (PBIC) provides seismic instrumentation to SCEC investigators for seismic research in southern California. Having control of instruments allows for rapid redeployment of the equipment in the event of a significant southern California earthquake. Past aftershock deployments have used PBIC equipment to supplement Southern California Seismic Network (SCSN) coverage and to obtain digital records at existing strong ground motion sites. PBIC instrumentation is compatible with Incorporated Research Institutes for Seismology (IRIS) PASSCAL equipment and has been used in several cooperative projects. In addition, the PBIC develops calibration and other quality control methods for use with the recording equipment and performs routine maintenance and repairs on seismic instrumentation for other SCEC institutions.

The PBIC provided equipment to nine different research experiments over this past year. Substantial field assistance was provided on several of the larger projects including the active phase of LARSE II and the Hector Mine deployment. The PBIC has completed the basic timing corrections and preliminary event associations of the Hector Mine data. Further processing steps are already underway and new data sets will be available in the next few months. This preliminary data set is available via the web at The Santa Barbara Array, delayed by the Hector Mine deployment and subsequent experiments, now has nine stations operational. The equipment was used heavily over the past year with the exception of a several month period following the completion of the Hector Mine deployment. This time was used to perform further maintenance on DASs and to upgrade all the GPS units to correct a bug introduced by the WNRO and Y2K upgrades performed in mid 1999. ARGOS satellite telemetry systems from PBIC and PASSCAL were used to monitor equipment status of the Hector Mine deployment and allowed researchers to better gauge battery replacement intervals and the rate of data collection.

The PBIC web page,, registered over 21000 hits from over 2000 different users and is being redesigned to provide easier and expanded access to user documentation and legacy information. Summaries and repair histories for all major pieces of PBIC equipment are now available online. Response information as well as the equipment information is being migrated to a web accessible database form.

Outreach programs continue to play an important role in the PBIC. ICS researchers and the PBIC participated in seismological demonstrations at local schools and hosted a booth at the Earth Day 2000 celebration in Santa Barbara at Santa Barbara Community College.

Ralph Archuleta

2/1/99 - 1/31/02

University of Southern California

Southern California Earthquake Center, USC 572726

Strong Motion Database (SMDB)

SMDB continues to see significant use from both the seismological and earthquake engineering communities. In 1999, more than 139,000 Web pages were accessed from the site; 3.7 GB of strong motion data were downloaded from SMDB; and 9,760 files were downloaded from off-site FTP and Web sites maintained by other agencies.

These numbers translate to an average daily access of 380 Web pages, download of 10 MB of strong-motion data from SMDB, and download of 27 files from external sites, representing increases of 83%, 23%, and 63% percent, respectively, over 1998 numbers (Figure 2).

Improvements to SMDB in 1999 include:

Data from thirty-four new earthquakes have been added with 440 station records. This is an increase of 32% for the number of station records over the beginning of the year. The new data includes records from the 10/16/99 Hector Mine earthquake (Figure 3), as well as from earthquakes in Alaska, Chile, Mexico, and Japan. Links to external FTP sites have been added from the University of Nevada Reno (for Guerrero Gap records) and from K-NET (for Japanese records). The database now contains 143 earthquakes and 1,796 site records.

Quality assessment has been conducted to improve the accuracy of important parameters in the database. All database values for epicentral distance, hypocentral distance, peak ground acceleration, response spectral amplitudes, and filter frequencies have been compared to values determined from the data and updated where discrepancies were found.

Query times have been reduced by up to 50% by optimizing code. For example, a query returning all data from the 1994 Northridge CA earthquake has been reduced from thirteen to seven seconds.

A new map access method has been added to the database ( This access method produce a map based on a user query showing all earthquakes and stations that resulted from that query. Clicking on an earthquake or station on the map produces a page showing a summary of data for that earthquake or station. The user can also zoom in or out on the map. An example is shown in Figure 3.

Navigation of the Web site has been improved by adding the major Web site links to the top of each page produced by the site.

Following requests from users, data summary pages now show sensor locations and component angles for each trace. These parameters are shown along with the station name, station owner, peak ground acceleration, and epicentral distance.

A new page has been added to the site listing all agencies that have contributed data to the database: ( This new page makes it easier for the user to identify the agency from which data has been collected and gives additional credit to the agencies involved.

A program to convert strong-motion data to the Seismic Analysis Code (SAC) format has been written. This program reads fourteen different ASCII strong-motion formats and converts them to one single format. The program is available from

To increase outreach to the engineering community, the database was presented at the EERI annual meeting in San Diego on February 5-7, 1999.

Daniel Lavallée


University of Southern California

Southern California Earthquake Center, USC 572726

Modeling of Nonlinear Strong Ground Motion During the 1994 Northridge
Earthquake at the Van Norman Dam Complex

SCEC sponsored studies and strong motion data recorded at the Van Norman Dam Complex (VNDC) during the 17 January 1994 Northridge earthquake (Mw 6.7) have been used to determine the presence of nonlinear effects in strong ground motion. The data available included borehole velocity profiles, weak to strong motion records, and dynamic soil laboratory tests. Measurements of ground motion observed at the Jensen Generator Building (JGB) site located on bedrock, have been coupled with the nonlinear model to generate scenarios of ground shaking at the Jensen Main Building (JMB) site. The nonlinear model formulation includes effects such as anelasticity, hysteretic behavior –also known as the memory effect– and cyclic degradation due to pore water pressure. In this study, two situations have been investigated; one that includes the effect of pore pressure and the second without it.

The results of the numerical experiments conducted with this site support the assumption that pore pressure cyclic mobility did contribute significantly to the ground shaking observed at the surface. A characteristic features of simulations including pore pressure effects for this site is the depletion of the high frequency in the signal in good agreement with the recorded acceleration. None of the scenarios generated without pore pressure were able to reproduce this feature. Also spiky waveforms in the tail can be detected in both synthetic and observed accelerograms. Similar conclusions applied to both horizontal components of the JMB site. We have to stress that in this study, conclusive finding of nonlinearity is based on direct simulations of nonlinear soil dynamics and the quantification of the nonlinearity in terms of the model parameters.

The model formulation includes nonlinear effects such as anelasticity, hysteretic behavior –also known as the memory effect– and cyclic degradation due to pore water pressure. The effect of pore pressure has been incorporated using a model developed by Towhata and Ishihara, 1985; and Iai, 1990a, b. The constitutive equation corresponds to a plane strain multiple mechanism model —or multi inelastic spring model— used to simulate cyclic mobility of sands under undrained conditions. In this model, the pore pressure development is correlated to the shear work.

Several hysteresis models have been developed and discussed in the literature (e.g., Pyke, 1979, Li and Liao, 1993; McCall, 1994; Muravskii and Frydman, 1998; Yoshida et al., 1998, Xu et al., 1998). Masing’s original formulation is inadequate to describe the hysteretic behavior of material under noncyclic loadings (see Pyke, 1979 and references therein; also Li and Liao, 1993). For instance, numerical simulations with noncyclic signals suggest that application of the Masing rules leads to an unphysical situation–such as the computed stress exceeding the strength of the material (see Figure 1).

To understand the behavior of the soil during strong shaking we have developed a general formulation of hysteresis based on the Masing rules (see Kramer, 1996). The generalized Masing rules provide a framework for understanding the nonuniform dilation and translation of stress-strain loops for a material subject to non-periodic stresses (or strains). This new hysteresis formulation has several interesting features. It has a functional representation and it includes the Cundall-Pyke hypothesis (Pyke, 1979) and Masing original formulation as special cases. In its most elementary implementation, the generalized Masing rule is even simpler than the Masing and extended Masing rules (Kramer, 1996). The model depends only on one free parameter named the fiducial point. This parameter controls the size of the loop in the stress-strain space and therefore can be related to the amount of energy dissipated through the nonlinear property of the material. We have derived a relationship between the anelastic damping of a stress-strain loop and the fiducial point for cyclic loadings (Archuleta et al., 1999). In other words, the generalized Masing rules provide a mean to introduce the effect of the damping ratio into nonlinear modeling independently of the other soil parameters. The theoretical formulation of the generalized Masing rules is presented in Archuleta et al., 1999; 2000. Incorporation of the pore pressure model with the Generalized Masing rule is discussed in Archuleta et al, 1999.

Peng-Cheng Liu


University of Southern California

Southern California Earthquake Center, USC 572726

Source Parameter Inversion for 1999 M 7.6 Chi-Chi Taiwan Earthquake
Using 3D Green's Functions

We will invert the near-field ground motion from the magnitude 7.6 Chi-Chi, Taiwan, earthquake on September 21, 1999 to determine the source parameters of finite fault by using 3D Green's functions. Because the 3D Green's functions play key role in this kind inversion, we should have the 3D velocity structure of the region around the 1999 Chi-Chi earthquake, as well as we should have technique to simulate wave propagation in 3D structure. Our study will base on the known 3D velocity structure for this region, instead of construct a new one, although, that is certain to be improved in the course of this study. We developed an efficient 3D viscoelastic finite difference (FD) algorithm in which we incorporated the improved modeling of absorbing boundaries and constant Q. In additional the 4th order, staggered-grid, 3D-FD algorithm is designed with discontinuous grid spacing so that computational and memory requirements can be significantly reduced. Having the 3D structure and numerical technique in hand, we can compute 3D synthetic Green's functions. A global hybrid search algorithm will then be employed to solve for the two components of fault slip, rise time, and rupture time for each subfault. For each time history we limit the time window to the direct P and S waves to mitigate the influence of the uncertainties in the 3D-velocity structure on later arrivals.


Kim Olsen

2/1/97 - 1/31/02

University of Southern California

Southern California Earthquake Center, USC 572726

Three-Dimensional Elastic Finite-Difference
Simulation of a Dynamic Rupture

Professor Raul Madariaga, Ecole Normale Superieure, and Dr. Eiichi Fukuyama, National Research Institute for Earth Science and Disaster Prevention, Japan, visited ICS during the summer of 2000, funded through OKP03.

In collaboration with Prof. Madariaga I have studied the conditions for rupture initiation and propagation along a fault surface subject to a heterogeneous stress distribution and/or inhomogeneous frictional parameters using a fourth-order finite-difference program. We have identified a non-dimensional parameter that controls rupture. This parameter measures the ratio of local available strain energy to fracture energy derived from the friction law. A bifurcation occurs when this parameter has values greater than a certain critical value that depends mildly on the geometry of the stress distribution on the fault. We have shown that the non-dimensional number can be used to describe dynamic rupture in the presence of a highly heterogeneous stress field for the 1992 Landers earthquake.

Together with Dr. Fukuyama and Prof. Madariaga I initiated the development of a hybrid method for flexible and efficient modeling of dynamic rupture propagation on curved fault geometry and its radiation in a heterogeneous three-dimensional medium. The dynamic rupture propagation is computed using the boundary integral equation (BIE) method. The computation of radiated waves away from the fault is carried out by an efficient fourth-order staggered-grid finite-difference (FD) method. The hybrid method enables dynamic modeling of rupture propagation on curved or multi-segmented faults in laterally and vertically heterogeneous earth models with an accurate free surface. In addition to the model of dynamic rupture on a single fault, the efficiency of the hybrid method can be used to model the statistics of recurrent ruptures on multiple, curved faults. We demonstrated the promise of the method by simulating spontaneous rupture with radiation on a curved fault.

Sophie Peyrat, a Ph.D. candidate at Ecole Normale Superieure in Paris visited ICS for a month during November and December 1999, funded through this project. During this time Sophie and I continued ongoing work together with Prof. Madariaga on simulated dynamic rupture for the 28 June 1992, {\bf M} 7.3 Landers earthquake. Our model used a single, vertical fault plane, combining the three segments of the Landers earthquake, a simple slip-weakening friction law, and a heterogeneous initial stress deduced from the slip inversion results of Wald and Heaton. While Olsen et al. (published in Science) obtained a fairly good fit between the dynamic radiation and strong motion records at four sites at various azimuths to the fault, we carried the work a step further by modifying the initial stress field derived from Wald and Heaton in a trial-and-error procedure, thereby obtaining a much improved fit between the dynamic radiation and strong motion records. In fact, we found that our dynamic radiation provided a much-improved fit to the strong motion data compared to that obtained by the Wald and Heaton kinematic model with a single fault plane involved.

During a visit at the National Research Institute for Earth Science and Disaster Prevention, Japan, during October-November 1999 I initiated a study on the variation of dynamic rupture velocity in homogeneous and heterogeneous media with Dr. Eiichi Fukuyama. This collaboration was continued during Dr. Fukuyama's visit at ICS July-August 2000. We use a long rectangular fault with primarily in-plane dynamics and the boundary integral and finite-difference methods for our study. On the homogeneous fault three regimes of rupture propagation appear. In regime I (large rupture resistance) the velocity increases slowly toward the Rayleigh wave speed. In regime II (intermediate rupture resistance) the velocity repeatedly jumps to super-shear values for a short distance but drops back to the Rayleigh wave speed. In regime III (low rupture resistance) velocity jumps to super-shear values and obtains values between the S- and P-velocity. Very narrow (critical) ranges of rupture resistance characterize the transitions between the three regimes. Such critical boundary is also observed in the presence of an asperity at the middle of the fault where the initial stress is raised by a certain amount. Here, the boundary includes two corners between which the critical boundary is located where the product of the width and height of the asperity becomes constant.

Kim Olsen


University of Southern California (PG&E), 050186

3D Modeling of the Santa Clara Basin

I have participated in the ground motion simulation code verification study in southern California. Four difference groups have compared synthetics for simulations in uniform halfspace and layered models with point and extended sources. The comparisons show excellent agreements between the synthetics from the different groups. Comparisons to be made include synthetics computed in the SCEC 3D-velocity model version 2.

Kim Olsen

2/1/99 - 1/31/02

University of Southern California

Southern California Earthquake Center, USC 572726

Ground Motion Modeling in Los Angeles

Ellen Gottschaemmer, a Ph.D. candidate at Karlsruhe University in Germany, visited ICS for a month during November and December 1999, funded through this project. During this time Ellen and I compared the accuracy of two implementations of the explicit planar free-surface boundary condition for 3-D fourth-order velocity-stress staggered-grid finite-differences, 1/2 grid apart vertically, in a uniform halfspace. Due to the staggered grid, the closest distance between the free surface and some wavefield components for both implementations is 1/2 grid spacing. When compared to a reflectivity solution computed at the staggered positions closest to the surface, the total misfit for all three components of the wavefield is generally found to be larger for the free surface co-located with the normal stresses, compared to that for the free-surface co-located with the xz and yz stresses. However, this trend is reversed when compared to the reflectivity solution exactly at the free surface. When the wavefield is averaged across the free surface, thereby centering the staggered wavefield exactly on the free surface, the free surface condition at the xz and yz stresses generates the smallest total misfit for increasing epicentral distance.

Sophie Peyrat, a Ph.D. candidate at Ecole Normale Superieure in Paris visited ICS for a month during November and December 1999. During this time Sophie and I simulated dynamic rupture for the 28 June 1992, M 7.3 Landers earthquake. Our model used a single, vertical fault plane, combining the three segments of the Landers earthquake, a simple slip-weakening friction law, and a heterogeneous initial stress deduced from the slip inversion results of Wald and Heaton. While Olsen et al. (published in Science) obtained a fairly good fit between the dynamic radiation and strong motion records at four sites at various azimuths to the fault, we carried the work a step further by obtaining a much improved fit between the dynamic radiation and strong motion records by modifying the initial stress field derived from Wald and Heaton in a trial-and-error procedure. In fact, Peyrat et al. found that their dynamic radiation provided a much-improved fit to the strong motion data compared to that obtained by the Wald and Heaton kinematic model with a single fault plane involved.


Jamison Steidl

2/1/98 - 1/31/02

University of Southern California

Southern California Earthquake Center, USC 572726

SCEC Borehole Instrumentation Initiative


Application of Site Characterization Studies to Ground Motion Prediction: Uphole/Downhole Modeling at the Van Norman Dam Complex

The 1997 Uniform Building Code (UBC) used in the design of structures by the engineering community places a great deal of emphasis on the average shear wave velocity in the upper 30 meters to classify sites and to assign site response correction factors. The emphasis on characterization of the near-surface properties in California, especially at sites with strong motion instrumentation, provides a wealth of new information for site response studies. Borehole geotechnical data coupled with downhole instrumentation allow for the direct estimation of the effects of surface geology on seismic ground motions and the ability to calibrate and improve our physical models of soil response for different levels of ground motion.

In March of 1997, a workshop was held to discuss the initiation of a borehole instrumentation program within SCEC to be coordinated with other ongoing drilling programs in Southern California. Shortly after the workshop the first year of the program was approved with three borehole sites planned for the first year, and three per year proposed for the following three years. The long term scientific objectives of this program are: to estimate the degree of nonlinearity for strong ground shaking on typical Southern California soils; to improve our capabilities in predicting the effects of near-surface soil conditions on ground motions; and to examine the details of the earthquake source process.

This marks the third year of the borehole instrumentation program. To date, seven boreholes have been drilled, logged and sampled for near-surface soil properties, and cased for deployment of the permanent borehole instrumentation. Four of the seven boreholes have been instrumented and are providing data real-time to the Caltech/USGS Southern California Seismic Network (SCSN). This data is available via the SCEC data center archive. These four boreholes are located on the exterior of the Los Angeles basin, three to the north within the Santa Monica Mtns. in a linear array with about 7 km spacing, and the fourth to the south-west on the Palos Verdes Peninsula. A fifth borehole located at the CDMG California strong motion instrumentation program (CSMIP) Obregon Park site in the LA Basin should be instrumented this fall and data will be streamed real-time to both the SCSN and CDMG in Sacramento. A 350 meter deep borehole has been drilled and cased in the LA Basin at the Long Beach water reclamation plant, and along with a shallow 30 meter borehole at the same location, will be instrumented this fall/winter. This vertical array will provide critical data on the soil behavior of a "typical" LA Basin site in both the small and large strain regimes as the waves propagate up through the soil column.

In addition to the Borehole Instrumentation Initiative, a second SCEC funded project uses portable borehole instrumentation at the Van Norman Dam site in the San Fernando Valley to help calibrate the current numerical wave-propagation methods for predicting soil response. This data can also be used for development of new nonlinear modeling techniques for the behavior of soils under large strains, as the instrumentation allows for ground motion up to 0.5-g downhole.

Using the borehole data from both SCEC projects as the input to our linear models, and the geotechnical site characterization data provided under collaboration with the ROSRINE project, we are able to reproduce the surface observations in the time, frequency, and response spectral domains for frequencies up to 10 Hz. The degree of non-linear behavior in California soils at large strain levels is a critical issue for determining the maximum plausible ground motions from large earthquakes. However, the lack of observed large-strain data from these new borehole sites in Southern California forces us to calibrate our nonlinear models with data from other regions. We use the Port Island vertical array observations from the Hyogo-ken Nambu (Kobe) earthquake and vertical array data from the Kushiro-Oki earthquake to test newly developed nonlinear soil behavior models.



Ralph Archuleta

11/26/95 - 11/25/01

Institute for Protection and Nuclear Safety, 4060-00001217

Institute for Protection and Nuclear Safety, 4060-9B016120/SH

Institute for Protection and Nuclear Safety, 4060-0A523570/SH

04/01/00 — 03/31/03

Nuclear Regulatory Commission, 04-00-038

Garner Valley Downhole Seismographic Array (GVDSA) Project

Project Background:

The Garner Valley downhole seismographic array (GVDSA) project, installed under the U.S. Nuclear Regulatory Commission contract NRC-04-87-108 in cooperation with the French Commissariat a l'Energie Atomique (CEA), has two main scientific objectives. The first is to understand the effects of the near-surface soil conditions on seismic ground motion to improve the ground motion prediction capabilities for design, seismic hazard assessment, and hazard mitigation. The second is to understand the effect of earthquake ground motion on the hydraulic conductivity of ground water systems for the deep storage of nuclear waste.

The near-surface geological site conditions have been shown to be the dominant factor in controlling the amplitude and variation of strong ground motion, and the damage patterns that result from large earthquakes. A unique set of data collected from the Garner Valley project makes it possible to advance two major areas of engineering seismology. The first problem is how weak motion scales to strong motion. The second one is how the recordings at different soil types scale to each other, especially with respect to a competent rock ("reference") site. The understanding of competing effects of amplification and attenuation (including non-linearity) is of a vital importance for seismic design studies. The site is located near the Anza segment of the seismically active San Jacinto fault in Southern California, which has a relatively high probability for a large earthquake of magnitude 6.5 or greater in the near future.

At Garner Valley we measure the ground motion during earthquakes in the bedrock 500 and 220 meters below the surface, at 50 meters below the surface in a zone of weathered granite, and at 22, 15, and 6 meters below the surface in a layer of soft alluvium. The ground motion is also measured at the surface above these borehole instruments by 5 stations in a linear array, one of which is directly above the borehole instruments.

In the deepest borehole (500 meter), downhole pressure transducers are located within sealed off fracture zones. The effects of earthquake ground motion and the rock-mass hydraulic response to ground motion are important factors in the short- and long-term performance of a high-level nuclear waste repository. This part of the project was designed and undertaken to provide fundamental data regarding the influence of earthquake ground motion on dynamic and static changes to the pore pressure in the rock mass. In addition to the downhole pressure transducers, static and dynamic changes in the pore pressure at different levels within the bedrock borehole are also measured via tubes (sampling lines) which are connected to pressure transducers at the surface and extend into the borehole to various depths.

Located 3 km from the GVDA main station is a remote rock station, which consists of surface and 30 meter borehole accelerometers connected to a Kinemetrics K2 data recorder. The data are stored on the K2 and transmitted via radio modem to the Garner Valley main station. The 30 meter borehole instrument and surface instrument directly above provide critical data on the effects of the weathered rock layer on the earthquake ground motion and help to define the "reference" or input motion for use in site-specific hazard analysis.

For more information on the Garner Valley Downhole Array project please see our website at

Current Research:

One of the highlights of the past year has been another significant donation to ICS. The operation and maintenance of a third vertical array observatory has been added to the engineering seismology test site program at ICS. The Borrego Valley Downhole Array (BVDA), donated to ICS in Fall of1999, is another state-of-the-art engineering seismology test site installed by Agbabian Associates for Kajima Corporation in 1993. This donation, while not officially valued yet by the Development Office at UCSB, is an even larger donation than the Hollister Earthquake Observatory (HEO), which was donated to ICS in 1998 and valued at 1.25 million dollars. It is the data from vertical arrays like GVDA, BVDA, and HEO that provide the calibration data for the extensive site characterization efforts of our geotechnical engineering colleagues, and continue to demonstrate to the earthquake engineering community the first class nature of engineering seismology research conducted at ICS. BVDA and HEO, along with GVDA, are all maintained under the NRC and IPSN funding agreements.

Another highlight of the last year’s research has been the data recorded from the M7.1 Hector Mine Earthquake. The vertical array data from GVDA and BVDA are sure to be the subject of future publications on engineering seismology research by ICS personnel. While the earthquake was located just over 100 km from the GVDA site, it was the first event which demonstrated clear evidence of pore pressure build up in the near-surface soil layers, even at relatively low strain levels. This kind of data is critical to our understanding of nonlinear soil behavior and the processes, which lead to soil liquefaction. In addition to the pore pressure build up in the near-surface soils, changes in the deep pore pressure were observed within the 500-meter borehole, which penetrates 400 meters into the country rock. These pore pressure changes are only just being analyzed and will also no doubt be the subject of future publications.

Extensive work on the development of new methods for incorporating nonlinear behavior into the dynamics of wave propagation through near-surface soils has taken place in the past years and have culminated in the successful defense of Fabian Bonilla’s Ph.D. Thesis. While simple linear models of soil response are well suited for modeling the low-strain observations we have seen to date at GVDA. HEO and BVDA, the degree of non-linear behavior in California soils at large strain levels is still a critical unresolved issue for determining the maximum plausible ground motions from large earthquakes. The data from the Hector Mine earthquake at the GVDA site which include pore pressure observations will be critical in validation and calibration of the nonlinear effective stress modeling of the GVDA site at the onset of the nonlinear soil behavior. The lack of observed large-strain data from borehole sites in California forces us to continue to calibrate our nonlinear models at large strain levels with data from other regions, such as Japan..


Ralph Archuleta


National Science Foundation, EAR 97-25709

Dynamic Earthquake Rupture Simulation on Dipping Faults

Our main goal in this project was to investigate the dynamics of dip-slip faulting, and specifically to characterize the effect of asymmetric geometry on the rupture propagation, fault slip, and resultant ground motion. For compressive tectonic regimes such as the Los Angeles area, Japan, and Central and South America, and in extensional regimes such as the Mediterranean and the Great Basin of Nevada, Utah, and Idaho, most seismic hazard lies in such non-vertical (dipping) faults. A general goal was to show that an important difference between a vertical and a non-vertical fault is the breakdown of symmetry with respect to the free surface. Because of this geometrical asymmetry in the non-vertical case, seismic waves radiated by an earthquake rupture can bounce off the free surface and hit the fault again and thus modify the stress field on the fault. This interaction causes variations in the normal stress on the fault, and these variations can affect the friction and hence the rupture and slip dynamics of the earthquake. The net result is that the time dependent normal stress produces asymmetric ground motion in the proximity of the fault. Using a numerical simulation technique, we sought to document this effect and explain its physical origin. We performed a number of dynamic simulations of earthquakes on non-vertical faults, using both two and three-dimensional finite element techniques. In the process, we developed a fault boundary condition for implementation in the finite element codes. These research activities were the bulk of David Oglesby’s dissertation research in the Department of Geological Sciences at UC Santa Barbara. The results have been presented in a variety of professional and education/outreach presentations, as well as in an article in Science magazine. Additional scientific journal publications will appear soon.

The asymmetry of the fault with respect to the free surface manifests itself chiefly by causing the normal stress on the fault to change with time. Ahead of the crack tip, due to the shear stress increase, the normal stress change is tensional for a normal fault and compressional for a thrust fault. This effect makes it more difficult for a thrust fault to break through to the surface, but this effect also causes a much higher energy release if the fault does break through. Behind the crack tip, in the slipping region of the fault, the stress changes are of opposite sign due to the drop from static to sliding friction on the fault. Therefore, the effect of the free surface on normal stress also changes sign: In the slipping region near the free surface, the normal stress on a normal fault is increased, while it is decreased for a thrust fault. After starting to slip, a normal fault will have a stronger frictional force holding it back, and have decreased particle motion. Conversely, a thrust fault will have lower friction, a greater stress drop, and increased particle motion. Additionally, for both thrust and normal faults, the hanging wall has higher ground motion than the footwall. This effect is due to the fact that the hanging wall has less volume and thus less mass than the footwall near the free surface.

The results of our simulation show that 1) in all cases the thrust fault produces higher ground motion than the normal fault on the free surface above the fault and 2) there is a large discontinuity in particle velocity and displacement as one crosses from the footwall to the hanging wall. The results of our simulations may explain some observations in the vicinity of non-vertical dip-slip faults, such as increased ground motion in the hanging wall compared to the footwall and larger ground motion for thrust faults than normal faults given the same site and source geometry. An important caveat is that the effects shown above exist only for faults that intersect the free surface of the Earth, not for buried faults.

We have also simulated the dynamics of dip-slip faults that have an abrupt change in dip with depth (a "bend" in the fault). We found that the pure dynamic effects of this geometry are small compared to the effects of the free surface mentioned before. The implications are that the deep structure of faults may not have a great effect on the near-fault strong ground motion.

Ralph Archuleta

Kim Olsen


US Geological Survey, 99HQGR0077

Earthquake Simulations and Ground Motion Prediction
for the San Francisco Bay Region

Dr. Ralph Archuleta, Dr. Peng-Cheng Liu and Dr. Kim Olsen have initiated a study to examine the effects of 3D crustal structure on kinematic slip inversion. They computed the ground motion using a finite-difference method and the reciprocity principle in 1D and 3D earth models of the 1979 Coyote Lake earthquake and inverted the tractions on the fault plane. Good comparisons of finite-difference and reflectivity synthetics validated the use of the finite-difference method. The inversion using the 3D Greens functions showed a better fit to data than obtained from the 1D Greens functions, suggesting that the forward and inverse modeling methodologies were appropriate, and that the accuracy of kinematic slip inversion may be improved by including the 3D earth structure.

Edward Keller

James Kennett


National Science Foundation, EAR-9803115

Investigation of a Very Rapid Tectonic Process: Direction and Rates of Lateral Propagation of Reverse Faulting and Folding

Major objectives during the past year included: 1) development of geomorphic criteria to recognize the process of lateral propagation of active folds in the Santa Barbara Fold Belt; and 2) initiation of work on the mid-channel trend associated with the Oak Ridge and mid-channel faults. Preliminary evidence suggests that the fold ,which has the largest topographic expression of any anticlinal fold in the Santa Barbara Channel is propagating laterally to the west..

Daniel Lavallée


National Science Foundation, 08990536

Modeling of Hysteretc and Anelastic Soils: From Laboratory Data of
Earthquake Strong Ground Motion

The main objective of this research project is to study nonlinear effects in earthquake strong ground motion. The premise is that nonlinear soil dynamics are essential to a full understanding of earthquake shaking. The model includes nonlinear effects such as anelasticity, hysteretic behavior (also known as the memory effect), and lost of stiffness due to pore water pressure.

To understand and reproduce the behavior of the soil during strong shaking we have developed a new formulation of hysteresis based on the Masing rules. The generalized Masing rules provide a framework for understanding the nonuniform dilation and translation of stress-strain loops for a material subject to non-periodic stresses (or strains). This new hysteresis formulation has several interesting features. It has a functional representation and in limiting cases it reduces to standard models of hysteresis. In its most elementary implementation, the generalized Masing rule is even simpler than the Masing and extended Masing rules. The model depends only on one free parameter named the fiducial point. This parameter controls the size of the loop in the stress-strain space and therefore can be related to the amount of energy dissipated through the nonlinear property of the material. We have derived a relationship between the anelastic damping of a stress-strain loop and the fiducial point for cyclic loadings. In other words, the generalized Masing rules provide a mean to introduce the effect of the damping ratio into nonlinear modeling independently of the other soil parameters.

Using the in situ observations from the Garner Valley downhole seismographic array, we have modeled scenarios of ground motions at the surface for this site. The simulations show amplitude reduction as well as the shift of the fundamental frequency to lower frequencies as observed on vertical arrays. The synthetic accelerograms show the development of intermittent behavior–high frequency peaks riding on low frequency carrier–as observed in the acceleration records mentioned above. Inclusion of pore pressure effects produces large strains in the soil with large amplification in the low-frequency band of the ground motion. Manifestation of spiky waveforms in the synthetic accelerogram is enhanced by the degradation of the soil parameters that is synchronous with a build up of pore pressure. Comparisons between the nonlinear model predictions with those computed with the equivalent linear model for the identical physical situation demonstrate that the latter model does not capture some essential elements of nonlinear soil response.

Daniel Lavallee


U.S. Geological Survey Award 01HQGR0049

Modeling and Estimates of Nonlinear Effects in Strong Earthquake Motion for the Los Angeles Area

SCEC sponsored studies and strong motion data recorded during the 17 January 1994 Northridge earthquake (Mw 6.7) have been used to determine the presence of nonlinear effects in strong ground motion. The data available included borehole velocity profiles, weak to strong motion records, and dynamic soil laboratory tests. Several sedimentary stations have responded with sharp resonance peaks. Signals recorded at these sites are good candidates to investigate nonlinear effects. Measurements of ground motion observed at the bedrock site Jensen Generator Building (JGB) have been coupled with the nonlinear model to generate scenarios of ground shaking at the Jensen Main Building (JMB) site. While measurements of ground motion observed at the bedrock site Pacoima Dam Downstream (PCD or PAC), have been coupled with the nonlinear model to generate scenarios of ground shaking at the Newhall Fire Station (NWH) site. The same procedure was used to generate scenarios of nonlinear motion at the Rinaldi Receiving Station (RIN) site using the ground motion recorded at the Los Angeles Dam (LAD) site. The results of this study confirm the presence of nonlinear effects in the mainshock observed at the JMB, NWH and RIN sites.

Kim Olsen


Los Alamos National Laboratory, F42200017-3Z

General Support on Computations on Near-Surface Wave Amplification in
the LA Basin

During 1997-1999, Dr. Eric Jones, Los Alamos Natl. Laboratories, and I have developed a hybrid method to model non-linear soil amplification using finite-fault ground motion. We compute low-frequency source-time functions at a datum plane below the area of interest using a 3-D elastic finite-difference method, and high-frequency source-time functions by a stochastic approach developed by Alexei Tumarkin. The deterministic and stochastic source-time functions are then merged and propagated up to the surface through a 1D-soil column using a full non-linear method. We have tested the method on the site NWH (Newhall Firestation) in the San Fernando Valley, where borehole velocity logs and laboratory strength data are available. We find linear/nonlinear ratios similar to those obtained from ratios of weak and strong motion recordings from the Northridge earthquake, supporting the proposed pervasive nonlinearity in the San Fernando Valley.

Kim Olsen


UC Institute of Geophysics & Planetary Physics, IGPP-99-048

Fault Rupture, Nonlinear Ground Motion
and Whole Basin Attenuation Modeling

Dr. Chris Bradley, Los Alamos Natl. Laboratories and Prof. Steve Day, San Diego state University has implemented the coarse-grained visco-elasticity approach in the 3D staggered-grid finite-difference code developed by Olsen (1994). We are in the process of estimating the improvement of the ground motion computation with a minimum S-wave velocity of 500 m/s compared to previous results using a minimum S-wave velocity of 1000 m/s in the SCEC Phase 3 report. The coarse-grained visco-elastic finite-difference code is being used for this analysis. Preliminary results show that synthetics with the lower minimum velocity provide an improved fit to the strong motion waveforms. The near-surface sediments (Vs of 0.5-1.0 km/s) clearly affect the ground motion in the Los Angeles region and should be included in the simulations. However, the relation Qs = 0.1 Vs (m/s) and Qp = 1.5 Qs seems to break down for velocities less than about 1 km/s. In order to estimate accurate ground motions, an improved estimate of the distribution of Qs and Qp must be obtained.

Kim Olsen

8/15/96 - 7/31/00

National Science Foundation, EAR -9628682

Simulation of Three-Dimensional Ground Motion
in Los Angeles from Large Earthquakes in Southern California

Dr. Paul Spudich of the USGS in Menlo Park and I have studied the effects of a 1-2 km wide low velocity zone in which the Calaveras fault is embedded. Strong motion accelerographs G06, located in the low velocity zone 1.2 km from the Calaveras fault, and G07, 4 km from G06, recorded both the M 6.2 Morgan Hill and the M 6.9 1989 Loma Prieta earthquakes. Comparison of the ground motions show that a large 0.6-1.0 Hz velocity pulse observed at G06 during the Morgan Hill event my be amplified by focusing caused by the low velocity zone. Such amplified waves might be a mappable seismic hazard, and the zone of increased hazard can extend as much as 1.2 km from the surface trace of the fault.

In collaboration with Dr. Nicolas Shapiro, now at University of Colorado at Boulder, who visited ICS in the spring of 2000, I have simulated seismic wave propagation in models of the subduction zone off the Pacific coast of Mexico using a three-dimensional finite difference method. The results show that a significant part of the energy generated by earthquakes in the subduction zone can be trapped inside a wave-guide generated by the low-velocity material of the accretionary prism. The trapped waves in the simulations appear as long-duration wave-trains that follow the onset of the surface waves, similar to the signature of broadband seismograms recorded by stations along the Pacific coast of Mexico for subduction earthquakes.

Dr. Shapiro and I have also used finite-difference simulations in simple two-dimensional (2D) models of the lithosphere to model the duration of long-period (<0.5 Hz) ground motion incident onto the Valley of Mexico for subduction zone earthquakes. Our simulations suggest that two heterogeneous structures extend the duration of the ground motion between the subduction zone and Mexico City by more than one minute: (1) the Mexican Volcanic Belt and (2) two low-velocity layers in the coastal region: the accretionary prism and the water layer. The duration generated by a crustal model including these structures is similar to that for earthquake records observed in between the coast and Mexico City. In the Valley of Mexico, our models including only regional-scale heterogeneity reproduce approximately one half of the observed duration. The results suggest that both the regional- and the local-scale low-velocity structures must be taken into account in order to explain the observed extended signal duration in the Valley of Mexico.

Kim Olsen


National Science Foundation Award EAR-0003275

Forward and Inverse Modeling of Rupture Dynamics in Three Dimesions

The proposed research is an integrated approach to significantly advance our knowledge in the field of earthquake rupture dynamics. The work consists of three parts. The first part is the development of a hybrid method for flexible and efficient modeling of dynamic rupture propagation on curved fault geometry and its radiation in a heterogeneous three-dimensional medium. The hybrid method enables dynamic modeling of rupture propagation on curved or multi-segmented faults in realistic earth models, which is necessary to understand the evolution of fault systems. The first tests of the method show promise while extensive code development still remains. The second part of the proposed research is a continuation of ongoing efforts on defining the parameters and conditions under which dynamic rupture starts, propagates, and stops. Preliminary results indicate that there is a simple non-dimensional parameter that controls rupture propagation. The implications of this simple number are wide ranging for rupture dynamics but require extensive tests, for example of the transitions of rupture speed, in earthquake models. The third part is to examine the feasibility of inverting for the friction or the initial stress field of a large earthquake. Preliminary studies show that the radiated waves are highly sensitive to the values of friction parameters as well as the initial stress field on the fault. If indeed dynamic friction is feasible, the proposed research may open up the most promising approach toward understanding the dynamics of earthquakes.

Toshiro Tanimoto


National Science Foundation, EAR-9972855

Nature & Cause of Long Period Background Oscillations

In 1997, we co-discovered the background oscillations along with some Japanese researchers. The whole Earth seems to be constantly shaking for frequencies between 2 and 7 millihertz, irrespective of earthquake occurrence; such a phenomenon was predicted by Benioff et al. (1959) before the earthquake-generated free oscillations were found in the records of the 1960 Chilean earthquake. But somehow it escaped detection by seismologists until 1997, almost 40 years later. This NSF grant supported the study of this interesting phenomenon. Through a series of papers, we have shown that these oscillations are caused by the atmosphere-solid Earth interaction; turbulent atmosphere seems to be hitting the surface of the Earth constantly, which results in ringing of the whole Earth body at its resonant frequencies. The critical evidence was seasonal variations, found in amplitudes of these excited oscillations (Tanimoto and Um, 1999). We developed a theory of atmospheric excitation (Tanimoto 1999, GJI paper), and analyzed data from global seismological network data (Tanimoto et al., 1999 in GRL and Tanimoto and Um 1999 in JGR). This new phenomenon was apparently intuitively appealing to the general public. It led to some extra-curricular activities such as quotes by a few general scientific journals, interview by a radio station and invitations to write articles for encyclopedia and journals for general public: they were (1) Science News (vol.154, July 1998) (2) NewScientist (September, 1999) (3) GEO (basically German version of Scientific American). (4) GRL paper (my first report on observation) was one of the four featured articles by Physics Today right after its publication (Physics Today. July, 1998) (4) Interviewed by National Public Radio in September, 1999. (5) An article for Encyclopedia (2001 McGraw-Hill). Toshiro Tanimoto

Toshiro Tanimoto


National Aeronautic & Science Administration, NAG5-7981

Atmospheric Excitation of Planetary Normal Modes
New Approach in Planetary Normal Mode Observation

This project expands the idea pursued in the Caltech grant and asks how the background oscillations are excited in other planets, in addition to Mars. If there is any other planet that has a large effect, a potentially new seismological approach to the planet may be planned in order to understand its interior structure.

Venus has been studied so far. Venus and Mars make an interesting contrast; in Mars, the atmosphere is thin (low density), but the wind is very strong. So the strong wind may play an important role in the excitation of background oscillations. In Venus, wind may be milder, but the atmosphere has high density and thus its force on the solid body has some inertial effects. In other words, the density may somewhat make up for the lack of winds. Our preliminary results indicate that background oscillations are excited as well as those in Mars.

We are extending this analysis further to other planets. We plan to come up with a comprehensive analysis of this effect for all planets.

Hazardous Waste

Lorne Everett


US Navy, N47408-96-A-7023

Support of the DOD National Environmental Technology Test Site for
Remediation Technologies

Lorne G. Everett provides commentary on technological approach and provides suggestions for increasing the technical quality of the proposed work in research, development and implementation of innovative remediation technologies at DOD sites.

Next Page || Previous Page || Table of Contents

Home Page || E-mail