ALBERT

Description: The inversion of local earthquake data (LED) for three-dimensional velocity structure requires the simultaneous solution of the coupled hypocenter-model problem. The Aki-Christoffersson-Husebye method (ACH) involves the inversion of large matrices, a task that is often performed by approximative solutions when the matrices become too big, as is the case for most LED, considering the coupled inverse problem. Such an approximate method (herein referred to as approximate geotomographic method) is used to perform tests with LED to obtain the best suited inversion parameters, such as velocity damping and number of iteration steps. The ACH method has been proposed for use of teleseismic data. Several adjustments to the original ACH method, which are necessary for use of LED, have been developed and are discussed. Such adjustments are the separation of the unknown hypocentral from the velocity model parameters for the inversion, the use of geometric weighting and step length weighting, the calculation of a minimum one-dimensional (1D) model as the starting three-dimensional (3D) model for the model inversion, and the display of an approximate resolution matrix (ray density tensors) before the inversion is performed. The ray density tensors allow the block cutting, e.g., the definition of the 3D velocity grid, to better correspond with the resolution capability of the specific data set. The adjustments to the method are tested by inversion of realistic LED of known variance. Synthetic LED are also used to demonstrate the effects of systematic errors, such as mislocations of seismic stations, on the resulting velocity field. Using the data sets from Long Valley, California, Yellowstone National Park, Wyoming, and Borah Peak, Idaho, the effects of improvements to the ACH method and of the data filtering process are shown. The use of the minimum 1D models for routine earthquake location improves this location procedure, as shown with the relocation of shots for the Long Valley and Yellowstone areas. The three-dimensional velocity fields obtained by the ACH method for the Long Valley and Yellowstone areas show local anomalies in the p velocity that can be correlated with tectonic and volcanic features. A pronounced anomaly of low p velocity below the Yellowstone caldera can be interpreted as a large magma chamber. However, the bulk of the paper addresses problems of the inversion method. The LED from the areas mentioned above are used to numerically and theoretically tune the inversion method for the defects that all real data contain. It is shown that one of the most important steps for any inversion of LED is the selection of the data for quality and for geometrical distribution.

Description: High-resolution seismic reflection and Sea Beam bathymetric data provide insights into the processes of sediment offscraping and accretion in the Middle America Trench off southern Mexico. Thick terrigenous sediments that are transported down Ometepec Canyon and accumulate along the trench floor are scraped off the oceanic plate and accreted in thrust packets to the lower trench slope. The packets offscraped represent most of the trench strata. Underlying hemipelagic deposits that accumulate on the seafloor seaward of the trench are subducted landward of the toe of the slope. Horizontal displacement on the thrust is less than 1 km. Leading edge folds are the surface expressions of the thrusts and strike subparallel to the base of the trench slope. The folds are continuous for as much as 10 km and have amplitudes as high as 200 m and wavelengths of 0.5 to 2 km. Folds are best developed along sections of the trench with interbedded silty turbidite and mud deposits. Fold are absent where thick coarse-grained fan deposits occur. Thickening of the thrust packets occurs by large-scale thrust duplication, by layer-parallel shortening, and by deposition of material that slumps off the leading edge of older upslope thrust blocks.

Description: We have determined the centroid depths and source mechanisms of 12 large earthquakes on transform faults of the northern Mid-Atlantic Ridge from an inversion of long-period body waveforms. The earthquakes occurred on the Gibbs, Oceanographer, Hayes, Kane, 15°20′, and Vema transforms. We have also estimated the depth extent of faulting during each earthquake from the centroid depth and the fault width. For five of the transforms, earthquake centroid depths lie in the range 7–10 km beneath the seafloor, and the maximum depth of seismic faulting is 14–20 km. On the basis of a comparison with a simple thermal model for transform faults, this maximum depth of seismic behavior corresponds to a nominal temperature of 900° ± 100°C. In contrast, the nominal temperature limiting the maximum depth of faulting during oceanic intraplate earthquakes with strike-slip mechanisms is 700° ± 100°C. The difference in these limiting temperatures may be attributed to the different strain rates characterizing intraplate and transform fault environments. Three large earthquakes on the 15°20′ transform have shallower centroid depths of 4–5 km and a maximum depth of seismic faulting of 10 km, corresponding to a limiting temperature of 600°C. The shallower extent of seismic behavior along the 15°20′ transform may be related to a recent episode of extension across the transform associated with the northward migration of the triple junction among North American, South American, and African plates to its present position near the transform. The source mechanisms for all events in this study display the strike-slip motion expected for transform fault earthquakes; slip vector azimuths agree to within 2°–3° of the local strike of the zone of active faulting. The only anomalies in mechanism were for two earthquakes near the western end of the Vema transform which occurred on significantly nonvertical fault planes. Secondary faulting, occurring either precursory to or near the end of the main episode of strike-slip rupture, was observed for five of the 12 earthquakes. For three events the secondary faulting was characterized by reverse motion on fault planes striking oblique to the trend of the transform. In all three cases the site of secondary reverse faulting is near a compressional jog in the current trace of the active transform fault zone. We find no evidence to support the conclusions of Engeln, Wiens, and Stein that oceanic transform faults in general are either hotter than expected from simple thermal models or weaker than normal oceanic lithosphere.

Description: A seismic array consisting of nine Hawaii Institute of Geophysics (HIG) ocean bottom seismometers (OBSs) was deployed at the eastern intersection of the Oceanographer Fracture Zone (OFZ) and the Mid‐Atlantic Ridge (MAR). The 12‐day experiment was designed to relate low‐magnitude earthquakes to the structure and tectonics of the MAR‐OFZ intersection. An average of 10 locatable events with duration‐based magnitudes between −1.0 and 2.0 were recorded per day. Excellent hypocentral locations of 112 events were obtained. Earthquake locations based on more than eight observations generally show 50% confidence volume constraints within OBS location errors. The earthquake locations cover a broad swath across the corner of the intersection zone. Magnitude‐weighted earthquake location likelihood maps suggest a decline in magnitudes near the intersection bathymetry low. Composite focal plane solutions suggest source mechanisms which indicate that the region is dominated by extensional tectonics. Alternative source solutions indicating translational movement are presented but are inconsistent with apparent bathymetric trends. The transition from the diverging (MAR) to translational (OFZ) plate margin occurs in the context of reduced magma genesis and crustal thinning due to the influence of the adjacent older lithosphere. The region may be described in terms of semirigid plate tectonics accompanying transform valley genesis.

Description: Two buoy types have been tested with respect to their drift performance under drogued and undrogued conditions. Additionally, forces acting on the buoys were measured directly. Quadratic drag laws have been confirmed for the drag in water and the combined drag of wind and waves. Stokes drift contributes about one half to the wind factor of 0.023, which is obtained for undrogued buoys in the Atlantic. The forces on a windowshade drogue are given by a linear relation between force and water velocity for speeds exceeding 10 cm/s. They have been extrapolated to speeds of less than 10 cm/s by both a linear and a quadratic relationship. Correlations between drift and wind speed in the Atlantic suggest that the linear law is a better approximation under realistic conditions. According to these measurements in the Atlantic the described buoy-drogue system with a windowshade drogue in 100-m depth is a good current-measuring device. Slippage is negligible for wind speeds of less than 15 m/s and is less than 2 cm/s under gale conditions. Undrogued buoys are strongly affected by wind and cannot be used for the analysis of currents without correction, even under light winds.

Description: The constitutive relations describing the fluid pressure response of a porous medium to changes in stress and temperature must reflect the microscopic processes that are operative over the time scale allowed for the deformation. Short‐duration deformations are readily described by undrained moduli, and intermediate duration deformations by drained moduli, both of which are formulated through linear elastic theory. Long‐term deformations that operate over geologic time are normally dominated by irreversible processes and result in considerably larger deformations, for the same applied stress conditions, than would be expected from their elastic counterparts. Model constitutive equations are developed for both elastic and irreversible processes and the magnitude and interpretation of the relevant material properties examined. Although the theory is presented in general terms, a sample calculation shows that for sandstone the inelastic deformation is one and one half orders of magnitude greater than the elastic deformation at the same applied stress. This difference in magnitude has a significant effect on the effective hydraulic diffusivity, various pore pressure coefficients, and the prospective fluid pressure development of the sediment.

Description: A parametric model is developed to describe relative permeability‐saturation‐fluid pressure functional relationships in two‐ or three‐fluid phase porous media systems subject to monotonic saturation paths. All functions are obtained as simple closed‐form expressions convenient for implementation in numerical multiphase flow models. Model calibration requires only relatively simple determinations of saturation‐pressure relations in two‐phase systems. A scaling procedure is employed to simplify the description of two‐phase saturation‐capillary head relations for arbitrary fluid pairs and experimental results for two porous media are presented to demonstrate its applicability. Extension of two‐phase relations to three‐ phase systems is obtained under the assumption that fluid wettability follows the sequence water 〉 nonaqueous phase liquid 〉 air. Expressions for fluid relative permeabilities are derived from the scaled saturation‐capillary head function using a flow channel distribution model to estimate effective mean fluid‐conducting pore dimensions. Constraints on model application are discussed.

Description: A worldwide investigation of continental erosion is carried out by the study of large drainage basins, on the basis of hydrological data, environmental factors, and basin relief distribution. Inside each basin, mean geochemical and mechanical denudation rates are defined. A multicorrelation analysis shows that the mechanical denudation rates Ds are uncorrelated with environmental factors and correlated with mean basin elevation H, while chemical denudation rates Dd are insensitive to relief but correlated with mean annual precipitation. Furthermore, two linear relationships between H and Ds are detected: (1) Ds (m/10³ yr) = 419×10−6 H (m) ‐ 0.245, with V (explained variance) = 95.1%; this law concerns basins related to orogenies younger than 250 Ma. The negative intercept is interpreted as a continental sedimentation rate of 245 m/m.y. An alternative model in which one invokes a critical elevation, separating erosion from sedimentation, is equally successful and leads to lower sedimentation rates (60–110 m/m.y.). For both models, one derives from the slope of the adjustments, erosion time constants on the order of 2.5 m.y. (2) Ds (m/10³ yr) = 61×10−6 H (m), with V = 86.5%; this law concerns basins related to older orogenies. The null intercept suggests the lack of continental storage. Because of the more important dispersion of the data, the erosion time constant is calculated separately for each basin; it ranges from 15 to 360 m.y. The tectonic implications of these results are discussed. In particular, the short time constant 2.5 m.y. agrees with orogenic uplift rates on the order of 1 mm/yr, observed in active mountain chains.

Description: An ocean carbon pump is defined as a process that depletes the ocean surface of σCO2 relative to the deep‐water σCO2. Three pumps are recognized: a carbonate pump, a soft‐tissue pump, and a solubility pump. The first two result from the biological flux of organic and CaCO3 detritus from the ocean's surface. The third results from the increased CO2 solubility in downwelling cold water and is demonstrated by a one‐dimensional upwelling‐diffusion model of an abiotic ocean. In the soft‐tissue and solubility pumps, working strengths are defined in terms of the ΔσCO2 each creates between surface and deep‐water. Efficiencies of each pump are quantified as a ratio of working strength to potential maximum strength. Using alkalinity, nitrate, and σCO2 to remove the carbonate pump signal from ocean or model data, the individual working strengths of the soft‐tissue and solubility pumps can be calculated by scaling the soft‐tissue's ΔσCO2 to the surface‐to‐deep ΔPO4. This technique is applied to a three‐box ocean model known to demonstrate high‐latitude control of atmospheric CO2 through a variety of circulation and biological changes. Considering each pump separately reveals that the various changes which lower pCO2atm in the model are caused primarily by an increased solubility pump. Analysis of global ocean data indicates a positive solubility pump signal, subject to uncertainties in the C:P Redfield ratio and in the preindustrial pCO2atm. If C:P = 105 and pCO2atm = 270 μatm, the efficiency of the solubility pump is about 0.5. We suggest that this type of analysis of relative carbon pump strengths will be an effective method for inter‐model and intra‐model comparison and diagnosis of underlying oceanic mechanisms for pCO2atm changes.

Description: The structure and evolution of the northern New Guinea collision zone is deduced from International Seismological Center (ISC) seismicity (1964–1985), new and previously published focal mechanisms and a reexamination of pertinent geological data. A tectonic model for the New Guinea margin is derived which illustrates the sequential stages in the collision and suturing of the Bewani‐Toricelli‐Adelbert‐Finisterre‐Huon‐New Britain arc to central New Guinea followed by subduction polarity reversal in the west. East of 149°E, the Solomon plate is being subducted both to the north and south; bringing the New Britain and Trobriand forearcs toward collision. West of 149°E the forearcs have collided, and together they override a fold in the doubly subducted Solomon plate lithosphere, which has an axis that is parallel to the strike of the Ramu‐Markham suture and that plunges westward at an angle of 5° beneath the coast ranges of northern New Guinea. Active volcanism off the north coast of New Guinea is related to subduction of the Solomon plate beneath the Bismarck plate. Active volcanism of the Papuan peninsula and Quaternary volcanism of the New Guinea highlands are related to slow subduction of the Solomon plate beneath the Indo‐Australian plate along the Trobriand Trough and the trough’s former extension to the west, respectively. From 144°–148°E, seismicity and focal mechanisms reveal that convergence between the sutured Bismarck and Indo‐Australian plates is accommodated by thrusting within the Finisterre and Adelbert ranges and compression of the New Guinea orogenic belt, together with basement‐involved foreland folding and thrusting to the south. The Finisterre block overthrusts the New Guinea orogenic belt, whereas the Adelbert block is sutured to New Guinea and overthrusts the oceanic lithosphere of the Bismarck Sea. Along the New Guinea Trench, west of 144°E, seismicity defines a southward dipping Wadati‐Benioif zone, and focal mechanisms indicate oblique subduction. Only this oldest, westernmost portion of the collision has progressed past suturing to a full reversal in subduction polarity.