Engineering topics
Publications and source records attributed to Stein, Seth.
The main objective of this project was to better understand the geodynamics controlling the mountain building and topographic evolution in the central Andes using an integrated approach that combines GPS, seismological, and numerical studies.
Although seafloor depth and heat flow for young oceanic lithosphere can be descibed by modeling the lithosphere as the boundary layer of a cooling halfspace, a long standing question has been why data at older ages deviate from those expected for a halfspace. Two classes of models have been proposed for these deviations. In one, heat added from below 'flattens' depth and heat flow. In the other, asthenospheric flow beneath the lithosphere perturbs the depths. We compare recent versions of the model classes: the GDH1 thin-lithosphere plate model (Stein and Stein, 1992) and an asthenospehric flow model (Phipps Morgan and Smith, 1992). The plate model fits heat flow data better than the flow model for all cases considered, and topographic data in all but one case. The flow model significantly overpredicts depths for the North Atlantic, because the assumed asthenospheric flow in the plate motion direction would yield deepening for old ages rather than the observed flattening. Overall, the GDH1 global average model does better than this flow model, whose parameters were fit to specific plates. Moreover, the plate models fit to specific plates do better than the flow model. Plate models thus appear more useful than this flow model, suggesting that deviations from a cooling halfspace are largely thermal in origin.
Spreading center segments that have experienced a complex tectonic history including rift propagation may have a complicated signature in bathymetric and magnetic anomaly data. To gain insight into the history of such regions, we have developed techniques in which both the magnetic anomaly patterns and seafloor fabric trends are predicted theoretically, and the combined predictions are compared numerically with the data to estimate best fitting parameters for the propagation history. Fitting functions are constructed to help determine which model best matches the digitized fabric and magnetic anomaly data. Such functions offer statistical criteria for choosing the best fit model. We use this approach to resolve the propagation history of the Cobb Offset along the Juan de Fuca ridge. In this example, the magnetic anomaly data prove more useful in defining the geometry of the propagation events, while the fabric, with its greater temporal resolution, is more useful for constraining the rate of propagation. It thus appears that joint inversion of magnetic and seafloor fabric data can be valuable in tectonic analyses.
A significant discrepancy exists between the heat flow measured at the seafloor and the higher values predicted by thermal models of the cooling lithosphere. This discrepancy is generally interpreted as indicating that the upper oceanic crust is cooled significantly by hydrothermal circulation. The magnitude of this heat flow discrepancy is the primary datum used to estimate the volume of hydrothermal flow, and the variation in the discrepancy with lithospheric age is the primary constraint on how the hydrothermal flux is divided between near-ridge and off-ridge environments. The resulting estimates are important for investigation of both the thermal structure of the lithosphere and the chemistry of the oceans. We reevaluate the magnitude and age variation of the discrepancy using a global heat flow data set substantially larger than in earlier studies, and the GDHI (Global Depth and Heat Flow) model that better predicts the heat flow. We estimate that of the predicted global oceanic heat flux of 32 x 10(exp 12) W, 34% (11 x 10(exp 12) W) occurs by hydrothermal flow. Approximately 30% of the hydrothermal heat flux occurs in crust younger than 1 Ma, so the majority of this flux is off-ridge. These hydrothermal heat flux estimates are upper bounds, because heat flow measurements require sediment at the site and so are made preferentially at topographic lows, where heat flow may be depressed. Because the water temperature for the near-ridge flow exceeds that for the off-ridge flow, the near-ridge water flow will be even a smaller fraction of the total water flow. As a result, in estimating fluxes from geochemical data, use of the high water temperatures appropriate for the ridge axis may significantly overestimate the heat flux for an assumed water flux or underestimate the water flux for an assumed heat flux. Our data also permit improved estimates of the 'sealing' age, defined as the age where the observed heat flow approximately equals that predicted, suggesting that hydrothermal heat transfer has largely ceased. Although earlier studies suggested major differences in sealing ages for different ocean basins, we find that the sealing ages for the Atlantic, Pacific, and Indian oceans are similar and consistent with the sealing age for the entire data set, 65 +/- 10 Ma. The previous inference of a young (approximately 20 Ma) sealing age for the Pacific appears to have biased downward several previous estimates of the global hydrothermal flux. The heat flow data also provide indirect evidence for the mechanism by which the hydrothermal heat flux becomes small, which has often been ascribed to isolation of the igneous crust from seawater due to the hydraulic conductivity of the intervening sediment. We find, however, that even the least sedimented sites show the systematic increase of the ratio of observed to predicted heat flow with age, although the more sedimented sites have a younger sealing age. Moreover, the heat flow discrepancy persists at heavily sedimented sites until approximately 50 Ma. It thus appears that approximately 100-200 m of sediment is neither necessary nor sufficient to stop hydrothermal heat transfer. We therefore conclude that the age of the crust is the primary control on the fraction of heat transported by hydrothermal flow and that sediment thickness has a lesser effect. This inference is consistent with models in which hydrothermal flow decreases with age due to reduced crustal porosity and hence permeability.
Recent revisions to the geomagnetic time scale indicate that global plate motion model NUVEL-1 should be modified for comparison with rates of motion estimated from space geodetic measurements. The optimal recalibration, which is a compromise among slightly different calibrations appropriate for slow, medium, and fast spreading rates is to multiply NUVEL-1 angular velocities by a multiplicative constant, alpha, of 0.9562. We refer to this simply recalibrated plate motion model as NUVEL-1A, and give correspondingly revised tables of angular velocities and uncertainties. Prior work has shown that space geodetic rates are slower on average than those calculated from NUVEL-1 by 6±1%. This average discrepancy is reduced to less that 2% when space geodetic rates are instead compared to NUVEL-1A.
We investigate angular velocity vectors of the Philippine Sea (PH) plate relative to the adjacent major plates, Eurasia (EU) and Pacific (PA), and the smaller Caroline (CR) plate. Earthquake slip vector data along the Philippine Sea plate are inverted, subject to the constraint that EU-PA motion equals that predicted by the global relative plate model NUVEL-1. The resulting solution fails to satisfy geological constraints along the Caroline-Pacific boundary: convergence along the Mussau Trench and divergence along the Sorol Trough. We then seek solutions satisfying both the CR-PA boundary conditions and the Philippine Sea slip vector data, by adjusting the PA-PH and EU-PH best fitting poles within their error ellipses. We also consider northern Honshu to be part of the North American plate and impose the constraint that the Philippine Sea plate subducts beneath northern Honshu along the Sagmi Trough in a NNW-NW direction. Of the solutions satisfying these conditions, we select the best EU-PH as 48.2 deg N, 157.0 deg E, 1.09 deg/my, corresponding to a pole far from Japan and south of Kamchatka, and PA-PH, 1.2 deg N, 134.2 deg E, 1.00 deg/my. Predicted NA-PH and EU-PH convergence rates in central Honshu are consistent with estimated seismic slip rates. Previous estimates of the EU-PH pole close to central Honshu are inconsistent with extension within the Bonin backarc implied by earthquake slip vectors and NNW-NW convergence of the Bonin forearc at the Sagami Trough.
Models for the thermal evolution of oceanic lithosphere are primarily constrained by variations in seafloor depth and heat flow with age. These models have been largely based on data from the Pacific and Atlantic Ocean basins. We construct seafloor age relations for the Indian Ocean which we combine with bathymetric, sediment isopach and heat flow data to derive curves for depth and heat flow versus age. Comparison of these curves with predictions from three thermal models shows that they are better fit by the shallower depths and higher heat flow for the GDH1 model, which is characterized by a thinner and hotter lithosphere than previous models.
Variations in sea-floor depth and heat flow with age provide the main constraints on the thermal structure and evolution of the oceanic lithosphere. Joint fitting of heat flow and bathymetry yields a model with a hotter, thinner lithosphere than in previous models. The new model provides a significantly better fit to the data, including those from older lithosphere previously treated as anomalous. This will facilitate the analysis of lithospheric processes, including the effects of mid-plate volcanism and swells, regional subsidence, and hydrothermal circulation near spreading centers.
Much of the success of plate tectonics can be attributed to the near rigidity of tectonic plates and the availability of data that describe the rates and directions of motion across narrow plate boundaries of about 1 to 60 kilometers. Nonetheless, many plate boundaries in both continental and oceanic lithosphere are not narrow but are hundreds to thousands of kilometers wide. Wide plate boundary zones cover approximately 15 percent of earth's surface area. Space geodesy, which includes very long baseline radio interferometry, satellite laser ranging, and the global positioning system, provides the accurate long-distance measurements needed to estimate the present motion across and within wide plate boundary zones. Space geodetic data show that plate velocities averaged over years are remarkably similar to velocities avaraged over millions of years.
A model for the present-day motion of the Rivera plate relative to the North America, Cocos, and Pacific plates is derived using new data from the Pacific-Rivera rise and Rivera transform fault, together with new estimates of Pacific-Rivera motions. The results are combined with the closure-consistent NUVEL-1 global plate motion model of DeMets et al. (1990) to examine present-day deformation in southwestern Mexico. The analysis addresses several questions raised in previous studies of the Rivera plate. Namely, do plate motion data from the northern East Pacific rise require a distinct Rivera plate? Do plate kinematic data require the subduction of the Rivera plate along the seismically quiescent Acapulco trench? If so, what does the predicted subduction rate imply about the earthquake recurrence interval in the Jalisco region of southwestern Mexico?
Results are reported from an investigation of the resolving power of ISC/NEIC P travel-time data in tomographic inversions for the geometry of the subduction zones in the NW Pacific. From thermal models for the Kurile, Janan, Izu-Bonin, Mariana, and Ryukyu slabs, three-dimensional synthetic velocity anomalies for subducting slabs are generated and projected onto a cell model for the uppermost 1400 km of the mantle. These synthetic models are used to compute synthetic delay times for ray paths corresponding to the source and receiver locations used for the actual data, add Gaussian noise, invert the synthetic data, and compare the resulting velocity structure to the initial synthetic models. This comparison is illustrated for sections through the Kuriles and the Mariana arcs. A variety of resolution artifacts are observed, which in many cases resemble features visible in the tomographic results obtained from inverting the actual ISC/NEIC data.
The current motions of the African, Eurasian, and North American plates are examined. The problems addressed include whether there is resolvable motion of a Spitsbergen microplate, the direction of motion between the African and North American plates, whether the Gloria fault is an active transform fault, and the implications of plate circuit closures for rates of intraplate deformation. Marine geophysical data and magnetic profiles are used to construct a model which predicts about 4 mm/yr slip across the Azores-Gibraltar Ridge, and west-northwest convergence near Gibraltar. The analyzed data are consistent with a rigid plate model with the Gloria fault being a transform fault.
Possible rigid plate models and shear models for the formation of curved seafloor lineaments by rift propagation are investigated and are applied to the Galapagos propagation rift system at 95 deg W. It is shown that the geometry noted at the Galapagos can result from either rift propagation acceleration or a spreading rate decrease during the last few hundred thousand years. It is postulated that the reverse curvature could result from either deceleration of rift propagation or an increase in spreading rate. The data interpreted as requiring a shear zone are found to be equally consistent with two distinct models based on rigid plate tectonics.
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.
The kinematics of rift propagation and the resulting goemetries of various tectonic elements for two plates is reviewed with no overlap zone. The formation and evolution of overlap regions using schematic models is discussed. The models are scaled in space and time to approximate the Easter plate, but are simplified to emphasize key elements. The tectonic evolution of overlap regions which act as rigid microplates and shear zones is discussed, and the use of relative motion and structural data to discriminate between the two types of models is investigated. The effect of propagation rate and rise time on the size, shape, and deformation of the overlap region is demonstrated.
The new NUVEL-1 data set for global relative plate motions is used here to discriminate between the two prevailing models for Caribbean plate motion. One model, by Jordan (1975), assumes that North America-Caribbean motion is reflected by the spreading rate inferred from magnetic anomalies at the Cayman Spreading Center and the azimuths of nearby transforms. The other model, by Sykes et al. (1982), uses rates and azimuths inferred from the geometry of the Lesser Antilles Wadati-Benioff zone. Overall, it is found that the data fit the Jordan geometry better, that the data used in global plate motion models are more suitable than rates and azimuths inferred from the geometry of the Wadati-Benioff zone for determining relative motions, and that incorporation of all relevant plate boundaries is essential.
A simple kinematic model is used to study the effects of various parameters on the evolution of zero-offset transforms and very-long-offset transforms. Consideration is given to the effects of initial configuration, degree of asymmetry, and degree of bias in asymmetry on the generation of these ridge transform geometries and on the possible steady-state nature of the transform length spectra. Of the parameters tests, only lack of 'memory' of zero-offset transforms affects transform length distribution.
A new bathymetric chart of part of the Chile transform system is constructed, based mainly on an R/V Endeavor survey from 100 deg W to its intersection with the East Ridge of the Juan Fernandez microplate. A generally continuous lineated trend can be followed through the entire region, with the transform valley being relatively narrow and well-defined from 109 deg W to approximately 104 deg 30 min W. The fracture zone then widens to the east, with at least two probable en echelon offsets to the south at 104 deg and 102 deg W. Six new strike-slip mechanisms along the Chile Transform and one normal fault mechanism near the northern end of the Chile Rise, inverted together with other plate-motion data from the eastern portion of the boundary, produce a new best-fit Euler pole for the Nazca-Antarctic plate pair, providing tighter constraints on the relative plate motions.