ALBERT

Description: What is the relationship between volcanic eruptions and climate change? More than 200 years after the connection was first proposed, it remains a thorny question. This article provides a brief historical overview of the problem and a review of the various data bases used in evaluating volcanic events and associated climatic change. We use the term “climate” to describe changes in the atmosphere over wide regions for periods of several months and longer. We use “weather” to describe shorter-term, variable atmospheric fluctuations experienced over more restricted areas. We appraise the present state of knowledge and highlight some pitfalls involved in using available information. Cautiously, we suggest future avenues for study, including the possibility of “volcanic winters,” or severe eruption-induced coolings.

Description: Magnetic lineation mapping in the western central Pacific has revealed a pair of opposite-sensed, fanned lineation patterns that define the accretionary boundaries of the fossil Magellan microplate. This tectonic synthesis results from extensive magnetic mapping of two new lineation patterns over a large area and extended mapping of previously identified lineations. The entire evolutionary history of the Magellan microplate is well constrained to a 9-m.y. period in the Early Cretaceous by synchronous spreading patterns and associated geologic data. During this period the microplate grew and evolved as a generally rectangular structure to a final size of 700 km×600 km with spreading centers on two opposing sides and transform faults on the other two sides. The lifetime and size of the Magellan microplate are somewhat longer and larger, respectively, than presently active microplates on the East Pacific Rise. However, these modern structures are still evolving and growing, and the tectonic behavior of the modern and Cretaceous systems appears to be similar. Study of both active and fossilized microplates should provide additional insights on their common tectonic histories. In particular, we show that the Magellan Trough spreading center behaved as an asymmetric accretionary plate boundary that can be described with two separate poles of motion very close to this spreading center during much of its history. The Magellan Trough spreading center then failed as a result of a larger ridge reorganization at the triple junction of the Pacific, Farallon, and Phoenix plates at Ml0N time. Microplate activity ceased when the microplate became welded to the Pacific plate at M9 time.

Description: The Pacific seafloor is littered with small fragments of lithosphere captured from adjacent plates by past plate boundary reorganizations. One of the clearest examples of such a reorganization is documented in the Mathematician Seamounts region, where a distinctive geomorphology and well-developed magnetic anomalies are present. This reorganization involved a short-lived microplate between the failing Mathematician Ridge and a new propagating spreading center: the East Pacific Rise. It produced a transfer of a fragment of lithosphere from the Farallon to the Pacific plate, and also created a number of landforms and magnetic patterns, within and on the margins of the captured fragment; these make up the Mathematician paleoplate. In many cases, two sides of a microplate are active spreading ridges. A microplate evolves into a paleoplate when dual spreading ceases and full spreading resumes at the prevailing spreading ridge. We define a paleoplate as the area of the seafloor, from the axis of a failed rift to the boundary of resumed, full spreading. It includes a fragment of captured lithosphere and the lithosphere slowly accreted to it during the period of dual spreading, prior to complete abandonment of the failed rift. The Mathematician paleoplate has the following boundaries and components from west to east: the axis of the Mathematician failed rift, the fragment of captured Farallon plate, a complex rift initiation site at the Moctezuma Trough, a zone of slow spreading, and an as yet ill-defined eastern boundary where dual spreading stopped and full spreading resumed. The northern boundary of the paleoplate is the Rivera fracture zone; its southeastern boundary a now-inactive transform fault, the West O'Gorman fracture zone. In this case, as well as in other more poorly documented ones, relict landforms and magnetic patterns are carried on the aging lithosphere, away from the spreading ridge, recording a former geometry.

Description: Methane carbon isotopic composition ranged from −76.9 to −62.6‰ in a tidal freshwater estuary (the White Oak River, North Carolina, United States) with site specific seasonal variations ranging from 6 to 10‰. During warmer months, tidally induced bubble ebullition actively transported this methane to the atmosphere. At two sites, these seasonally varying fluxes ranged from 1.2 ± 0.3 to 1.3 ± 0.3 mol CH4 m−2yr−1 (19.2 to 20.8 g CH4m−2yr−1), with flux-weighted average isotopic compositions at two sites of −66.3 ± 0.4 and −69.5 ± 0.6‰. The carbon isotopic composition of naturally released bubbles was shown to be indistinguishable from the sedimentary methane bubble reservoir at three sites, leading to the conclusion that isotopic fractionation did not occur during the ebullition of methane. The hypothesis was developed that ebullitive methane fluxes are depleted in 13CH4 relative to fluxes transported via molecular diffusion or through plants, as zones of 13C enriching microbial methane oxidation are bypassed.

Description: We examine the diffusive behavior of the flow field in an eddy-resolving, primitive equation circulation model. Analysis of fluid particle trajectories illustrates the transport mechanisms, which are leading to uniform tracer and potential vorticity distributions in the interior of the subtropical thermocline. In contrast to the assumption of weak mixing in recent analytical theories, the numerical model indicates the alternative of tracer and potential vorticity homogenization on isopycnal surfaces taking place in a nonideal fluid with strong, along-isopycnal eddy mixing. The eastern, ventilated portion of the gyre appears to be sufficiently homogeneous to allow the concept of an eddy diffusivity to apply. A break in a random walk behavior of particle statistics occurs after about 100 days when along-flow dispersion sharply increases, indicative of mean shear effects. During the first months of particle spreading, eddy dispersal and mean advection are of similar magnitude. Eddy kinetic energy is of O(60–80 cm2 s−2) in the model thermocline, comparable to the pool of weak eddy intensity found in the eastern parts of the subtropical oceans. Eddy diffusivity in the model thermocline (Kxx = 8 × 107, Kyy = 3 × 107 cm2 s−1) seems to be higher by a factor of about 3 than oceanic values estimated for these area. Below the thermocline, model diffusivity decreases substantially and becomes much more anisotropic, with particle dispersal preferentially in the zonal direction. The strong nonisotropic behavior, prominent also in all other areas of water eddy intensity, appears as the major discrepancy when compared with the observed behavior of SOFAR floats and surface drifters in the ocean.

Description: Replacement dolomitization by seawater has been modeled in order to quantify the Sr-isotope signature in Cenozoic dolomites as a function of precursor mineralogy and 87Sr/86Sr ratio, reaction stoichiometry and 87Sr/86Sr ratio of the dolomitizing fluids. High Sr carbonates, such as aragonite, may introduce a significant precursor memory into an otherwise seawater dominated Sr-isotope signature if small quantities of seawater per unit volume of precursor carbonate are involved. Dolomitization of low Sr carbonates (i.e. low-Mg calcite) are shown to create an isotopic signature indistinguishable from that of the seawater involved in the reaction. Therefore, by comparison with the Sr-isotope evolution curve of seawater, the- 87Sr/86Sr ratios of the dolomites can be used to record the oldest possible age of dolomitization and the youngest age of deposition. The implications for this approach have been applied to data obtained from a dolomitized core from Little Bahama Bank, Bahamas. Two periods of dolomitization are recognized, one in the early Late Miocene involving Middle Miocene or older rocks, and a second one around 2.4 Ma ago affecting early Pliocene carbonates.

Description: A radiocarbon-calibrated box model for today's ocean suggests that a lag of about 1750 years should exist between the arrival of the midpoint of the deglaciation 18O signal in the deep Atlantic Ocean and its arrival in the deep Pacific Ocean. In order to assess the actual lag, we have carried out accelerator radiocarbon measurements on two cores from the Atlantic Ocean and one core from the Pacific Ocean. Although the results are not definitive, there is a suggestion that the actual time lag was about 1000 years.

Description: Based on detailed reconstructions of global distribution patterns, both paleoproductivity and the benthic δ13C record of CO2, which is dissolved in the deep ocean, strongly differed between the Last Glacial Maximum and the Holocene. With the onset of Termination I about 15,000 years ago, the new (export) production of low- and mid-latitude upwelling cells started to decline by more than 2-4 Gt carbon/year. This reduction is regarded as a main factor leading to both the simultaneous rise in atmospheric CO2 as recorded in ice cores and, with a slight delay of more than 1000 years, to a large-scale gradual CO2 depletion of the deep ocean by about 650 Gt C. This estimate is based on an average increase in benthic δ13C by 0.4–0.5‰. The decrease in new production also matches a clear 13C depletion of organic matter, possibly recording an end of extreme nutrient utilization in upwelling cells. As shown by Sarnthein et al., [1987], the productivity reversal appears to be triggered by a rapid reduction in the strength of meridional trades, which in turn was linked via a shrinking extent of sea ice to a massive increase in high-latitude insolation, i.e., to orbital forcing as primary cause.

Description: The eastern part of the North Atlantic subtropical gyre is found in the region between the Azores and the Cape Verde Islands. A study of the gyre structure in the area east of 35°W between 8°N and 41°N is presented. The geostrophic flow field determined from historical temperature-salinity data sets by objective analysis indicates seasonal variations in shape but no significant changes in the magnitude of volume transports. The eastern part of the gyre has a larger east-west and smaller north-south extension in summer compared with the winter season. The center shifts by about 2° latitude to the south from winter to summer. Long-term temperature time series (6.5 years) from a mooring near the Azores are consistent with these results, showing always a consistent temperature increase at the beginning of the year which is apparently due to the displacement of the northeastern part of the gyre. A comparison between the mean flow fields and fields obtained from individual zonal sections indicates large deviations north and south of the gyre but small deviations within the gyre.

Description: It has long been recognized that the transition from the last glacial to the present interglacial was punctuated by a brief and intense return to cold conditions. This extraordinary event, referred to by European palynologists as the Younger Dryas, was centered in the northern Atlantic basin. Evidence is accumulating that it may have been initiated and terminated by changes in the mode of operation of the northern Atlantic Ocean. Further, it appears that these mode changes may have been triggered by diversions of glacial meltwater between the Mississippi River and the St. Lawrence River drainage systems. We report here Accelerator Mass Spectrometry (AMS) radiocarbon results on two strategically located deep-sea cores. One provides a chronology for surface water temperatures in the northern Atlantic and the other for the meltwater discharge from the Mississippi River. Our objective in obtaining these results was to strengthen our ability to correlate the air temperature history for the northern Atlantic basin with the meltwater history for the Laurentian ice sheet.

Description: Radiocarbon ages for benthic and planktonic foraminifera from the late glacial sections of two Atlantic and two Pacific cores are reported. The differences for benthic-planktonic pairs suggest that the radiocarbon age for deep Atlantic water was somewhat larger than today's (i.e., 600±250, as opposed to 400 years) and that the radiocarbon age for deep Pacific water was also slightly larger than today's (2100±400, as opposed to 1600, years). Our results suggest that during glacial time, the deep Pacific was, as it is today, significantly depleted in radiocarbon relative to the deep Atlantic. As many questions remain unanswered regarding the reliability of this approach, these conclusions must be considered to be preliminary.

Description: As a test of the reliability of paleocean ventilation rates reconstructed from radiocarbon age differences between planktonic and benthic foraminifera, measurements have been made on coexisting species of planktonic foraminifera. While ideally no differences should exist, we do find them. In this paper we discuss the possible causes for these differences and attempt to evaluate their impact on the interpretation of benthic-planktonic age differences.

Description: The total mass of sediments on the ocean floor is estimated to be 262 × 1021 g. The overall mass/age distribution is approximated by an exponential decay curve: (11.02 × 1021 g)e−0.0355t Ma. The mass/age distribution is a function of the area/age distribution of ocean crust, the supply of sediment to the deep sea, and submarine erosion and redeposition. About 140 × 1021 g of the sediment on the ocean floor is pelagic sediment, consisting of about 74% CaCO3, with the remainder opaline silica and red clay. Of the sediment on the ocean floor, 122 × 1021 g is detritus, mostly terrigenous, but a small portion (about 6 × 1021 g) is volcanic. Because very little pelagic sediment is obducted, virtually all of the pelagic sediment mass and some fraction of the terrigenous sediment is being subducted at a rate estimated to be about 1 × 1021 g per million years. The composition of sediment on the ocean floor differs significantly from that of average passive margin and continental sediment, so that the loss of ocean floor sediment through subduction may drive the composition of global sediment toward enrichment in silica, alumina, and potash and toward depletion in calcium.

Description: Sedimentological, isotopic and magnetostratigraphic investigations of Ocean Drilling Program and Deep Sea Drilling Project sites 642, 643, 644 and 610 document the oceanographic and climatic evolution of the Norwegian Sea and the northeastern Atlantic over the last 2.8 m.y.. The results show that a major expansion of the Scandinavian Ice Sheet to the coastal areas took place at about 2.56 Ma. Relatively severe glacials appeared until about 2 Ma. The period 2.6 ‐ 1.2 Ma experienced in general cold surface water conditions with only a weak influx of temperate Atlantic water as compared with late Quaternary interglacials. The Norwegian Sea was a sink of deep water through this period but deepwater ventilation was reduced and calcite dissolution was high compared with the Holocene. Deep water formed by other mechanisms than it does today. Between 2 and 1.2 Ma the glaciations in Scandinavia were small. A transition toward larger glacials took place during the period 1.2 to 0.6 Ma, corresponding to warmer interglacials and reduced calcite dissolution. Only during the last 0.6 m.y. has the oceanographic and climatic system of the Norwegian Sea varied in the manner described in previous studies of the late Quaternary. A strong thermal gradient was present between the Norwegian Sea and the northeastern Atlantic during the Matuyama (2.5–0.7 Ma). This is interpreted as a sign of a more zonal and less meridional climatic system over the region compared with the present situation.

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: 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: The Gulf of California is a long, narrow marginal sea lying between the Baja California peninsula and the mainland of Mexico. Air‐sea fluxes of heat and moisture in the gulf are enhanced because of geographical isolation from the Pacific provided by the mountainous Baja California peninsula. In the northern gulf, annual evaporation rates are about 1 m y−1. Unlike most evaporative basins, however, the gulf gains heat from the atmosphere at an annual average rate of 20 to 80 W m−2 (Bray, 1988). Given the unusual air‐sea forcing of the gulf, what form or forms should water mass formation take? The annual moisture loss and heat gain require that high‐salinity surface water be transported downward to an intermediate depth and that cold, fresh inflow be transported upward to an intermediate depth. This is accomplished through several mechanisms. (1) Winter convection: this occurs only in a limited geographical region, the Wagner Basin of the far northern gulf, except in El Niño–Southern Oscillation years, when convection appears to be more widespread. (2) Dispersion of convected water in small eddylike features: this occurs within a large‐scale southward transport possibly driven by the large‐scale density gradient associated with atmospheric fluxes. (3) An anticyclonic circulation in the northern gulf: this is found south of the convection region and transports high‐salinity water off the shallow shelves and to substantial depths, where it may mix with water of central gulf origin. (4) Tidal mixing: most of the energy available for mixing in the northern gulf derives from the tides. In particular, tidal mixing over the sill in the island region is responsible for the substantial reduction in salinity of northern gulf waters as they enter the central gulf.

Description: A coupled ocean-atmosphere general circulation model has been developed for TOGA related problems. The coupled model consists of an ocean model of the tropical Pacific and a global low-order spectral atmosphere model. The two models interact via wind stress and sea surface temperature. In order to avoid a climate drift within the coupled model, a flux correction method is applied.Experiments were performed by introducing a westerly wind stress burst over the western equatorial Pacific for one month. Thereafter, the wind burst is turned off and the response of the coupled model to the initial disturbance is investigated. The results are compared with the response of the ocean model run with the same disturbance in an uncoupled mode.It is shown that the coupling leads to a significant increase of the duration of anomalous conditions in the ocean. SST anomalies persist for about 12 months in the coupled run, while they have already disappeared after 4 months in the uncoupled case. The increase in persistence is due to the feedback of the atmosphere, which responds with an eastward shift of the ascending branch of the Walker Circulation.In a second experiment with the coupled model the initial disturbance was introduced within another season. The results show no basic differences to the results of the first experiment.An interesting result of the coupled model runs is the occurrence of spontaneous westerly wind bursts over the western Pacific, which developed by internal dynamics. Location and duration of these spontaneous wind bursts show some correspondence with the time-space structure of observed westerly wind stress episodes over the western Pacific.

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: 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.