Richard J. Blakely

Bio: Richard J. Blakely is an academic researcher from United States Geological Survey. The author has contributed to research in topics: Fault (geology) & Magnetic anomaly. The author has an hindex of 39, co-authored 141 publications receiving 7810 citations. Previous affiliations of Richard J. Blakely include Oregon State University.

Potential theory in gravity and magnetic applications

Abstract: Introduction 1. The potential 2. Consequences of the potential 3. Newtonian potential 4. Magnetic potential 5. Magnetization 6. Spherical harmonic analysis 7. Regional gravity fields 8. The geomagnetic field 9. Forward method 10. Inverse method 11. Fourier-domain modeling 12. Transformations A. Review of vector calculus B. Subroutines C. Review of sampling theory D. Conversion of units.

Approximating edges of source bodies from magnetic or gravity anomalies

Abstract: Cordell and Grauch (1982, 1985) discussed a technique to estimate the location of abrupt lateral changes in magnetization or mass density of upper crustal rocks. The final step of their procedure is to identify maxima on a contoured map of horizontal gradient magnitudes. We attempt to automate their final step. Our method begins with gridded magnetic or gravity anomaly data and produces a plan view of inferred boundaries of magnetic or gravity sources. The method applies to both local surveys and to continent-wide compilations of magnetic and gravity data (e.g., Zietz, 1982; Simpson et al., 1983a; Kane et al., 1982).

Fore-arc migration in Cascadia and its neotectonic significance

Abstract: Neogene deformation, paleomagnetic rotations, and sparse geodetic data suggest the Cascadia fore arc is migrating northward along the coast and breaking up into large rotating blocks. Deformation occurs mostly around the margins of a large, relatively aseismic Oregon coastal block composed of thick, accreted seamount crust. This 400-km-long block is moving slowly clockwise with respect to North America about a Euler pole in eastern Washington, thus increasing convergence rates along its leading edge near Cape Blanco, and creating an extensional volcanic arc on its trailing edge. Northward movement of the block breaks western Washington into smaller, seismically active blocks and compresses them against the Canadian Coast Mountains restraining bend. Arc-parallel transport of fore-arc blocks is calculated to be up to 9 mm/yr, sufficient to produce damaging earthquakes in a broad deformation zone along block margins.

Basin‐centered asperities in great subduction zone earthquakes: A link between slip, subsidence, and subduction erosion?

Abstract: [1] Published areas of high coseismic slip, or asperities, for 29 of the largest Circum-Pacific megathrust earthquakes are compared to forearc structure revealed by satellite free-air gravity, bathymetry, and seismic profiling. On average, 71% of an earthquake's seismic moment and 79% of its asperity area occur beneath the prominent gravity low outlining the deep-sea terrace; 57% of an earthquake's asperity area, on average, occurs beneath the forearc basins that lie within the deep-sea terrace. In SW Japan, slip in the 1923, 1944, 1946, and 1968 earthquakes was largely centered beneath five forearc basins whose landward edge overlies the 350°C isotherm on the plate boundary, the inferred downdip limit of the locked zone. Basin-centered coseismic slip also occurred along the Aleutian, Mexico, Peru, and Chile subduction zones but was ambiguous for the great 1964 Alaska earthquake. Beneath intrabasin structural highs, seismic slip tends to be lower, possibly due to higher temperatures and fluid pressures. Kilometers of late Cenozoic subsidence and crustal thinning above some of the source zones are indicated by seismic profiling and drilling and are thought to be caused by basal subduction erosion. The deep-sea terraces and basins may evolve not just by growth of the outer arc high but also by interseismic subsidence not recovered during earthquakes. Basin-centered asperities could indicate a link between subsidence, subduction erosion, and seismogenesis. Whatever the cause, forearc basins may be useful indicators of long-term seismic moment release. The source zone for Cascadia's 1700 A.D. earthquake contains five large, basin-centered gravity lows that may indicate potential asperities at depth. The gravity gradient marking the inferred downdip limit to large coseismic slip lies offshore, except in northwestern Washington, where the low extends landward beneath the coast. Transverse gravity highs between the basins suggest that the margin is seismically segmented and could produce a variety of large earthquakes.

A new isostatic residual gravity map of the conterminous United States with a discussion on the significance of isostatic residual anomalies

Abstract: To display more clearly the gravity anomalies caused by geologic bodies in the upper parts of the crust, a new colored isostatic residual gravity map of the conterminous United States has been prepared using the gravity data set compiled for the Gravity Anomaly Map of the United States (Society of Exploration Geophysicists, 1982). The new isostatic residual gravity map is based on an Airy-Heiskanen model of local compensation, in which the surface load requiring compensation is defined by 5-min topographic and bathymetric data sets. A colored first-vertical-derivative map of isostatic residual gravity further enhances the short-wavelength anomalies produced by bodies at or near the surface and emphasizes the regional fabrics and trends in the gravity field. For the purpose of displaying gravity anomalies caused by shallow bodies of geologic significance, the nature of the isostatic model and the values of its parameters are of lesser importance than the application of an isostatic correction of some sort. Most isostatic models result in residual gravity maps that appear nearly identical in their main patterns and features. Anomalies on isostatic residual gravity maps should not be casually interpreted in terms of “undercompensation” or “overcompensation” because large-amplitude anomalies can be produced by crustal bodies in complete local isostatic equilibrium. Many isostatic residual gravity anomalies less than several hundred kilometers wide can be related to known geologic bodies. We present here a classification scheme that attempts to categorize such anomalies on the basis of tectonic environment. In general, highs correlate with intruded or accreted mafic material or with upthrusted crustal sections, whereas lows occur over low-density sedimentary or volcanic sections, felsic intrusive bodies, or down-warped crustal sections. Although some longer-wavelength anomalies, such as the broad gravity high centered over Montana, could be manifestations of density contrasts deep in the mantle, many such anomalies can also be modeled by geologically reasonable density contrasts that are isostatically compensated and confined to depths of less than several hundred kilometers, so that their source bodies need not be deep. The fact that certain of these broader anomalies have well-defined boundaries which correlate with near-surface geologic features increases the likelihood that their sources lie entirely within the lithosphere. If so, then the density contrasts required to explain the gravity data imply fundamental anomalies in chemical composition or thermal state for the crustal and upper mantle columns under these regions. We have investigated spectral analysis as a method to quantitatively characterize regional anomaly patterns. Contoured plots of normalized amplitude spectra were prepared for various areas of the isostatic residual gravity field of the United States. These Fourier domain representations show characteristic patterns that can be interpreted in terms of the trends and wavelengths of anomalies and may help to more objectively distinguish geologic basements with different origins or tectonic histories.

Current plate motions

Abstract: A global plate motion model, named NUVEL-1, which describes current plate motions between 12 rigid plates is described, with special attention given to the method, data, and assumptions used Tectonic implications of the patterns that emerged from the results are discussed It is shown that wide plate boundary zones can form not only within the continental lithosphere but also within the oceanic lithosphere; eg, between the Indian and Australian plates and between the North American and South American plates Results of the model also suggest small but significant diffuse deformation of the oceanic lithosphere, which may be confined to small awkwardly shaped salients of major plates

Current plate motions

Abstract: SUMMARY We determine best-fitting Euler vectors, closure-fitting Euler vectors, and a new global model (NUVEL-1) describing the geologically current motion between 12 assumed-rigid plates by inverting plate motion data we have compiled, critically analysed, and tested for self-consistency. We treat Arabia, India and Australia, and North America and South America as distinct plates, but combine Nubia and Somalia into a single African plate because motion between them could not be reliably resolved. The 1122 data from 22 plate boundaries inverted to obtain NUVEL-1 consist of 277 spreading rates, 121 transform fault azimuths, and 724 earthquake slip vectors. We determined all rates over a uniform time interval of 3.0m.y., corresponding to the centre of the anomaly 2A sequence, by comparing synthetic magnetic anomalies with observed profiles. The model fits the data well. Unlike prior global plate motion models, which systematically misfit some spreading rates in the Indian Ocean by 8–12mm yr−1, the systematic misfits by NUVEL-1 nowhere exceed ∼3 mm yr−1. The model differs significantly from prior global plate motion models. For the 30 pairs of plates sharing a common boundary, 29 of 30 P071, and 25 of 30 RM2 Euler vectors lie outside the 99 per cent confidence limits of NUVEL-1. Differences are large in the Indian Ocean where NUVEL-1 plate motion data and plate geometry differ from those used in prior studies and in the Pacific Ocean where NUVEL-1 rates are systematically 5–20 mm yr−1 slower than those of prior models. The strikes of transform faults mapped with GLORIA and Seabeam along the Mid-Atlantic Ridge greatly improve the accuracy of estimates of the direction of plate motion. These data give Euler vectors differing significantly from those of prior studies, show that motion about the Azores triple junction is consistent with plate circuit closure, and better resolve motion between North America and South America. Motion of the Caribbean plate relative to North or South America is about 7 mm yr−1 slower than in prior global models. Trench slip vectors tend to be systematically misfit wherever convergence is oblique, and best-fitting poles determined only from trench slip vectors differ significantly from their corresponding closure-fitting Euler vectors. The direction of slip in trench earthquakes tends to be between the direction of plate motion and the normal to the trench strike. Part of this bias may be due to the neglect of lateral heterogeneities of seismic velocities caused by cold subducting slabs, but the larger part is likely caused by independent motion of fore-arc crust and lithosphere relative to the overriding plate.

Present-day plate motions

Abstract: A data set comprising 110 spreading rates, 78 transform fault azimuths, and 142 earthquake slip vectors has been inverted to yield a new instantaneous plate motion model, designated Relative Motion 2 (RM2). The model represents a considerable improvement over our previous estimate, RM1 [Minster et al., 1974]. The mean averaging interval for the spreading rate data has been reduced to less than 3 m.y. A detailed comparison of RM2 with angular velocity vectors which best fit the data along individual plate boundaries indicates that RM2 performs close to optimally in most regions, with several notable exceptions. The model systematically misfits data along the India-Antarctica and Pacific-India plate boundaries. We hypothesize that these discrepancies are manifestations of internal deformation within the Indian plate; the data are compatible with northwest-southeast compression across the Ninetyeast Ridge at a rate of about 1 cm/yr. RM2 also fails to satisfy the east-west trending transform fault azimuths observed in the French-American Mid-Ocean Undersea Study area, which is shown to be a consequence of closure constraints about the Azores triple junction. Slow movement between North and South America is required by the data set, although the angular velocity vector describing this motion remains poorly constrained. The existence of a Bering plate, postulated in our previous study, is not necessary if we accept the proposal of Engdahl and others that the Aleutian slip vector data are biased by slab effects. Absolute motion models are derived from several kinematical hypotheses and compared with the data from hot spot traces younger than 10 m.y. Although some of the models are inconsistent with the Wilson-Morgan hypothesis, the overall resolving power of the hot spot data is poor, and the directions of absolute motion for the several slower-moving plates are not usefully constrained.

Present‐day plate motions

Abstract: A data set comprising 110 spreading rates, 78 transform fault azimuths and 142 earthquake slip vectors was inverted to yield a new instantaneous plate motion model, designated RM2. The mean averaging interval for the relative motion data was reduced to less than 3 My. A detailed comparison of RM2 with angular velocity vectors which best fit the data along individual plate boundaries indicates that RM2 performs close to optimally in most regions, with several notable exceptions. On the other hand, a previous estimate (RM1) failed to satisfy an extensive set of new data collected in the South Atlantic Ocean. It is shown that RM1 incorrectly predicts the plate kinematics in the South Atlantic because the presently available data are inconsistent with the plate geometry assumed in deriving RM1. It is demonstrated that this inconsistency can be remedied by postulating the existence of internal deformation with the Indian plate, although alternate explanations are possible.

A new geomagnetic polarity time scale for the Late Cretaceous and Cenozoic

Abstract: We have constructed a magnetic polarity time scale for the Late Cretaceous and Cenozoic based on an analysis of marine magnetic profiles from the world's ocean basins. This is the first time, since Heirtzler et al. (1968) published their time scale, that the relative widths of the magnetic polarity intervals for the entire Late Cretaceous and Cenozoic have been systematically determined from magnetic profiles. A composite geomagnetic polarity sequence was derived based primarily on data from the South Atlantic. Anomaly spacings in the South Atlantic were constrained by a combination of finite rotation poles and averages of stacked profiles. Fine-scale information was derived from magnetic profiles on faster spreading ridges in the Pacific and Indian Oceans and inserted into the South Ariantic sequence. Based on the assumption that spreading rates in the South Atlantic were smoothly varying but not necessarily constant, a time scale was generated by using a spline function to fit a set of nine age calibration points