SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Corrective opposing vibrations reduce Earthquake damage. Algorithm for controlling reactions of buildings and other structures to Earthquake initiates counterpulses when incoming shock pulses exceed threshold level. By using "go/no-go" criterion rather than continuously adjusting response to every input excitation no matter how small, control system needs less computation effort than comparable closed-loop systems.
Some of the correlations concerning the properties of subduction zones are reviewed. A quantitative global comparison of many subduction zones reveals that the largest earthquakes occur in zones with young lithosphere and fast convergence rates. Maximum earthquake size is directly related to the asperity distribution on the fault plane. This observation can be translated into a simple model of seismic coupling where the horizontal compressive stress between two plates is proportional to the ratio of the summed asperity area to the total area of the contact surface. Plate age and rate can control asperity distribution directly through the horizontal compressive stress associated with the vertical and horizontal velocities of subducting slabs. The basalt to eclogite phase change in the down-going oceanic crust may be largely responsible for the uncoupling of subduction zones below a depth of about 40 km.
Repeated surveys of short level lines in the Shumagin Islands, Alaska, reveal coherent tilt signals associated with subduction of the Pacific plate beneath the North American plate in the Shumagin seismic gap. Ten years of steady tilt down toward the trench is interrupted during 1978-1980 by a rapid episode of reverse tilt. The 'normal' tilt represents surface deformation as subduction occurs, with the plate boundary locked to at least 60 km depth. Using all available tilt, sea level, and seismic data, the tilt reversal is interpreted as due to an episodic reverse slip of about 80 cm magnitude on the plate boundary between about 70 km and 20 km depth, downdip from the main seismogenic zone, which remains locked. This event causes an increase of stress on the locked main thrust zone. It is speculated that such events may be a regular process at subduction zones, that great plate boundary earthquakes may be more common during their occurrence, and that their onset may be detectable early enough to give warning of an increase in probability for the occurrence of a great earthquake.
The progress of the NASA Crustal Dynamics Project is summarized. The Project is performed to measure regional deformation and strain accumulation due to large earthquakes in the plate boundary region of western N. America, to measure relative plate motions around the globe, to measure internal deformations of continental and oceanic lithospheric plates, to relate the rotational dynamics of the earth to earthquakes, plate motions and other geophysical phenomena and to measure regional fault motions and strain accumulation globally. VLBI and satellite laser ranging (SLR) have been employed using several base and mobile stations in N. America, Europe, S. America and the S. Pacific. Satellite data are being gathered with the LAGEOS spacecraft. Measurement campaigns in S. America and the Carribbean and Mediterranean regions will be carried out with new instrumentation in 1985.
Because the constitutive laws for soils are governed mainly by interparticle friction, all aspects of their mechanical behavior depend strongly on gravitational body forces. This fact poses serious limitations on the formulation of a materially objective soil constitutive theory, based on experimentation performed on Earth. In particular, the presence of the Earth's gravity prohibits the design of controlled experiments to properly simulate a variety of critical phenomena associated with the dynamic response of soils to seismic excitation in a very low effective confining stress field. For these reasons, laboratory-controlled experiments in the Space Shuttle, under essentially zero-gravity conditions, could offer invaluable opportunities for developing a quantitative understanding of fundamental aspects of soil behavior during or after an earthquake, which, in turn, could result in significant technological advances in geotechnical earthquake engineering.
Stressing and rupture of a locked zone adjacent to a creeping fault segment was studied with special reference to strength heterogeneity depthwise and along-strike. The resulting precursory temporal and spatial variations of surface strain rate profiles were compared to geodetic measurements on the San Andreas fault in central California. Crustal deformation in great California earthquake cycles was also studied with special reference to the temporal decay of strain rate observed since the 1957 and 1906 great earthquakes, and comtemporary surface strain rate and velocity profiles at several locations along the San Andreas. The effect of viscoelastic response in the deep aseismic shear zone on the surface deformation behavior was examined. Work was begun on a fundamental reformulation of the crustal deformation problem focusing on the crustal deformation process affected by deep aseismic slip as the slip zone progresses toward an instability and as deep seismic slip continues postseismically, the 3-D nature of the problem due to geometry and material heterogeneity, and the time-dependent source coming from the lithosphere/astenospheric coupling process.
The assessment of earthquake hazards and mineral and oil potential of a given region requires a detailed knowledge of geological structure, including the configuration of faults. Delineation of faults is traditionally based on three types of data: (1) seismicity data, which shows the location and magnitude of earthquake activity; (2) field mapping, which in remote areas is typically incomplete and of insufficient accuracy; and (3) remote sensing, including LANDSAT images and high altitude photography. Recently, high resolution radar images of tectonically active regions have been obtained by SEASAT and Shuttle Imaging Radar (SIR-A and SIR-B) systems. These radar images are sensitive to terrain slope variations and emphasize the topographic signatures of fault zones. Techniques were developed for using the radar data in conjunction with the traditional types of data to delineate major faults in well-known test sites, and to extend interpretation techniques to remote areas.
The hypothesis that impact-induced faults have been preferred sites for later deformation in response to lithospheric stresses has been suggested for several planets and satellites. This hypothesis is investigated on earth by examining whether terrestrial impact structures show higher rates of nearby earthquake activity than do surrounding intraplate regions. For 28 of 30 probable impact structures having an original crater 20 km or more in diameter, the rates of nearby seismicity have been no higher than the regional background rates. For two large probable impact structures, Vredefort and Charlevoix, with higher than normal rates of nearby seismicity, factors other than slip on impact-induced faults appear to control the occurrence of earthquakes. It is concluded that impact-induced faults, at least on earth, do not persist as lithospheric 'weak zones' for periods in excess of several million years after the impact event.
It is hypothesized that crustal deformation occurs on a scale-invariant matrix of faults. For simplicity, a two-dimensional pattern of hexagons on which strike-slip faulting occurs is considered. The behavior of the system is controlled by a single parameter, the fractal dimension. Deformation occurs on all scales of faults. The fractal dimension determines the fraction of the total displacement that occurs on the first-order or primary faults. The value of the fractal dimension can be obtained from the frequency-magnitude relation for earthquakes. The results are applied to the San Andreas fault system in central California. Earthquake studies give D = 1.90. The main strand of the San Andreas fault is associated with the primary faults of the fractal system. It is predicted that the relative velocity across the main strand is 2.93 cm/yr. The remainder of the relative velocity of 5.5 cm/yr between the Pacific and North American plates occurs on higher-order faults. The predicted value is in reasonably good agreement with the value 3.39 + or - 0.29 cm/yr obtained from geological studies.
The average rates of line-length change and the rates of the shallow fault slip observed near Parkfield, California, were used to estimate the slip rate at depth on the San Andreas fault (SAF), using the Matsu'ura et al. (1986) procedure to model the aseismic deformation as a perturbation to steady block motion across the SAF. Five different models were tested against the data, with one model, featuring a nearly monotonic transition from creeping behavior on the northwest to locked behavior on the southeast, found to fit the data acceptably well. It is shown that stored elastic energy southeast of Parkfield could produce an earthquake of magnitude up to 8.0, that would have displacement and rupture length similar to historic earthquakes of the same size.
Intermediate and deep focus earthquakes in Wadati-Benioff zones are considered, and it is suggested that although seismicity in downgoing slabs is restricted to their cold and strong interior, the seismicity distribution is not easily described as being controlled by a single lithospheric strength value. A numerical thermal model is used to compute strength contours for flow law constants reported from laboratory experiments, and it is shown that the expected pressure strengthening is large enough that the slab should have considerable strength well below the deepest seismicity. The results indicate that either a strongly depth-dependent limiting strength in downgoing slabs exists, or factors in addition to strength control the distribution of subduction zone earthquakes.
Comparison of angles from historical triangulation observations dating as far back as 1932 with Global Positions System (GPS) measurements taken in 1987 indicates that rapid convergence may be taking place on decade timescales in the central and eastern part of the Ventura basin, an east-west trending trough bounded by thrust faults. Changes in angles over this time were analyzed using Prescott's modified Frank's method and in terms of a model which assumes that the regions to the north and south of the basin are rigid blocks undergoing relative motion. For the two block model, inversion of the observed angle changes over the last 28 years for the relative motion vector leads to north-south convergence across the basin of 30 + or - 5 mm/yr, with a left lateral component of 10 + or - 1 mm/yr in the Fillmore-Santa Paula area in the central part of the basin. The modified Frank's method yields strain rates of approximately 2 microrad/yr in both the east and central parts of the basin for measurements spanning the 1971 San Fernando earthquake. Assuming no east-west strain yeilds north-south compression of approximately 3.5 + or - .2 cm/yr. Comparison of triangulation data prior to the earthquake shows no strain outside the margin of error. The convergence rates determined by geodetic techniques are consistent with geologic observations in the area. Such large geodetic deformation rates, with no apparent near-surface creep on the major thrust, can be understood if these faults become subhorizontal at relatively shallow depths and if the subhorizontal portions of the faults are creeping. An alternative explanation of the large displacement rates might be that the pumping of oil in the vicinity of the benchmarks caused large horizontal motions, although it is unlikely that meter scale horizontal motions are due to oil withdrawal. These and other hypotheses are evaluated to better constrain the tectonics of this active region.
The analysis of earthquake fault lineament patterns within the alluvial valley of west Tennessee, which is often made difficult by the presence of unconsolidated sediments, is presently undertaken through a synergistic use of Landsat satellite images in conjunction with gravity anomaly data, which were quantitatively analyzed and compared by means of two-dimensional histograms and rose diagrams. The northeastern trend revealed for the lineaments corresponds to faults and is in keeping with reactivation of the Reelfoot rift near the Mississippi River; this suggests that deeper features, perhaps at earthquake focal depth, may extend to the land surface as Landsat-detectable lineaments.
Differential radar interferometry uses data from synthetic-aperture radar (SAR). Three passes of SAR yield three amplitude-and-phase images, from which two interferograms (phase-difference images) made. Interferograms used to make third, "double-difference" interferogram indicating vertical motion of terrain between passes. Vertical earthquake motions as small as 1 cm detectable. Used to make extensive, accurate maps of such geophysical phenomena as heaving and buckling in fault zones, motions to tectonic plates, residual displacements from earthquakes, motions from prevolcanic swelling, motions of glaciers, tides, and thermal expansion of mountains from diurnal heating.
A recent reanalysis of the International Latitude Serivce (ILS) polar motion data-day has been processed using Kalman filtering techniques to generate the polar motion excitation function over the time-span from 1960 to 1965. The resulting excitation function has been examined for the effects of 1960 Chile in an attempt to determine experimentally how large earthquake affect polar motion. The resulting upper bound of about 75 x 10 to the 22nd N-m for a 10-deg dip (about 36 x 10 to the 22nd N-m for 20-deg dip) is consistent with results obtained from previous seismic studies, including a recent normal mode excitation result. Following future great earthquakes, monitoring of polar motion by space-based techniques such as VLBI should continue at high temporal resolution for several weeks in order to directly measure the rheological parameters of the upper mantle.
Shallow dip angles (not greater than 45 deg) suggested by field observations of continental extensional faults are not predicted by classical isotropic Mohr-Coulomb-Anderson theory. Earthquake data indicate that normal faults exist in the upper crust with dip angles commonly as shallow as 30 deg. One explanation of the apparent conflict between theory and structural/seismic observations is that intrinsically weak, shallow-dipping pre-existing faults are preferentially reactivated. An analytical treatment clarifies the roles of geometry, intact/preexisting fault strengths and fluid pore pressures. Frictional strength ratios of 3 or greater could account for extremely shallow normal faults (dips 10-20 deg) without consideration of pore pressures in excess of the least principal stress or of principal stress systems rotated away from the gravity vector. Moderate reduction in friction (3/4) with respect to wall rock can reduce the dip to 30 deg and can account for shallow normal-slip earthquakes.
It is shown here that seismicity around the margins of deglaciated areas provides a constraint on the viscosity of the lower mantle. Calculations using a spherical, viscoelastic earth model show that the present-day magnitude of the stress fields induced in the lithosphere beneath the Laurentide and Fennoscandian ice sheets is very sensitive to the value of the lower-mantle viscosity. Stress of about 100 bar, sufficient to cause seismicity, can still remain in the lithosphere for lower-mantle viscosities greater than about 10 to the 22nd Pa-s; for lower-mantle viscosities of about 10 to the 21st Pa-s, only a few tens of bars of stress persist in the lithosphere today. This influence of lower-mantle viscosity on the state of stress in the lithosphere also has implications for the migration of stress from earthquakes, and hence for earthquake recurrence times.
A novel iterative technique, based on the concept of fault segmentation and computed using 2D static dislocation theory, for building models of seismicity and fault interaction which are physically acceptable and geometrically and kinematically correct, is presented. The technique is applied in two steps to seismicity observed at the Middle America Trench. The first constructs generic models which randomly draw segment strengths and lengths from a 2D probability distribution. The second constructs predictive models in which segment lengths and strengths are adjusted to mimic the actual geography and timing of large historical earthquakes. Both types of models reproduce the statistics of seismicity over five units of magnitude and duplicate other aspects including foreshock and aftershock sequences, migration of foci, and the capacity to produce both characteristic and noncharacteristic earthquakes. Over a period of about 150 yr the complex interaction of fault segments and the nonlinear failure conditions conspire to transform an apparently deterministic model into a chaotic one.