Petrology and geochronology of basalt breccia from the 1996 earthquake swarm of Loihi seamount, Hawaii: magmatic history of its 1996 eruption

 Samples of basalt were collected during the Rapid Response cruise to Loihi seamount from a breccia that was probably created by the July to August 1996 Loihi earthquake swarm, the largest swarm ever recorded from a Hawaiian volcano. ²¹⁰Po–²¹⁰Pb dating of two fresh lava blocks from this breccia indicates that they were erupted during the first half of 1996, making this the first documented historical eruption of Loihi. Sonobuoys deployed during the August 1996 cruise recorded popping noises north of the breccia site, indicating that the eruption may have been continuing during the swarm. All of the breccia lava fragments are tholeiitic, like the vast majority of Loihi's most recent lavas. Reverse zoning at the rim of clinopyroxene phenocrysts, and the presence of two chemically distinct olivine phenocryst populations, indicate that the magma for the lavas was mixed just prior to eruption. The trace element geochemistry of these lavas indicates there has been a reversal in Loihi's temporal geochemical trend. Although the new Loihi lavas are similar isotopically and geochemically to recent Kilauea lavas and the mantle conduits for these two volcanoes appear to converge at depth, distinct trace element ratios for their recent lavas preclude common parental magmas for these two active volcanoes. The mineralogy of Loihi's recent tholeiitic lavas signify that they crystallized at moderate depths (∼8–9 km) within the volcano, which is approximately 1 km below the hypocenters for earthquakes from the 1996 swarm. Taken together, the petrological and seismic evidence indicates that Loihi's current magma chamber is considerably deeper than the shallow magma chamber (∼3–4 km) in the adjoining active shield volcanoes. Received: 21 August 1997 / Accepted: 15 February 1998     

Contrasting origins of Cenozoic silicic volcanic rocks from the western Cordillera of the United States

Two fundamentally different types of silicic volcanic rocks formed during the Cenozoic of the western Cordillera of the United States. Large volumes of dacite and rhyolite, mostly ignimbrites, erupted in the Oligocene in what is now the Great Basin and contrast with rhyolites erupted along the Snake River Plain during the Late Cenozoic. The Great Basin dacites and rhyolites are generally calc-alkaline, magnesian, oxidized, wet, cool (<850°C), Sr-and Al-rich, and Fe-poor. These silicic rocks are interpreted to have been derived from mafic parent magmas generated by dehydration of oceanic lithosphere and melting in the mantle wedge above a subduction zone. Plagioclase fractionation was minimized by the high water fugacity and oxide precipitation was enhanced by high oxygen fugacity. This resulted in the formation of Si-, Al-, and Sr-rich differentiates with low Fe/Mg ratios, relatively low temperatures, and declining densities. Magma mixing, large proportions of crustal assimilation, and polybaric crystal fractionation were all important processes in generating this Oligocene suite. In contrast, most of the rhyolites of the Snake River Plain are alkaline to calc-alkaline, ferroan, reduced, dry, hot (830–1,050°C), Sr-and Al-poor, and Nb-and Fe-rich. They are part of a distinctly bimodal sequence with tholeiitic basalt. These characteristics were largely imposed by their derivation from parental basalt (with low fH₂O and low fO₂) which formed by partial melting in or above a mantle plume. The differences in intensive parameters caused early precipitation of plagioclase and retarded crystallization of Fe–Ti oxides. Fractionation led to higher density magmas and mid-crustal entrapment. Renewed intrusion of mafic magma caused partial melting of the intrusive complex. Varying degrees of partial melting, fractionation, and minor assimilation of older crust led to the array of rhyolite compositions. Only very small volumes of distinctive rhyolite were derived by fractional crystallization of Fe-rich intermediate magmas like those of the Craters of the Moon-Cedar Butte trend. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Precise Relative Location of 25-ton Chemical Explosions at Balapan Using IMS Stations

—?We test how well low-magnitude (m _bLg 1.8 to 2.6), 25-ton chemical explosions at Balapan, Kazakhstan, can be located using IMS stations and standard earth models, relying on precisely determined relative arrival times of nearly similar, regional and teleseismic waveforms. Three 1997 Balapan explosions were recorded by a number of currently reporting and surrogate IMS stations. Three regional stations and two teleseismic arrays yielded consistent waveforms appropriate for relative picking. Master-event locations based on the AK135 model and ground-truth information from the first, shallowest and best-recorded explosion, fell under 1 km from known locations, for depths constrained to that of the master event. The resulting 90% confidence ellipses covered 12–13?km² and contained the true locations; however, results for depth constrained to true depth were slightly less satisf actory. From predictions based on ground truth, we found a P _g -coda phase at Makanchi, Kazakhstan to be misidentified and poorly modeled. After accounting for this, 90% ellipses shrank to 2–3?km² and true-depth mislocation vectors became more consistent with confidence-ellipse orientations. These results suggest that a high level of precision could be provided by a tripartite array of calibration shots in cases where models are poorly known. We hope that the successful relocation of these small Balapan shots will support the role of calibration explosions in verification monitoring and special event studies, including on-site inspection.     

An Efficient Numerical Method for Computing Synthetic Seismograms for a Layered Half-space with Sources and Receivers at Close or Same Depths

— It is difficult to compute synthetic seismograms for a layered half-space with sources and receivers at close to or the same depths using the generalized R/T coefficient method (Kennett, 1983; Luco and Apsel, 1983; Yao and Harkrider, 1983; Chen, 1993), because the wavenumber integration converges very slowly. A semi-analytic method for accelerating the convergence, in which part of the integration is implemented analytically, was adopted by some authors (Apsel and Luco, 1983; Hisada, 1994, 1995). In this study, based on the principle of the Repeated Averaging Method (Dahlquist and Björck, 1974; Chang, 1988), we propose an alternative, efficient, numerical method, the peak-trough averaging method (PTAM), to overcome the difficulty mentioned above. Compared with the semi-analytic method, PTAM is not only much simpler mathematically and easier to implement in practice, but also more efficient. Using numerical examples, we illustrate the validity, accuracy and efficiency of the new method.     

The Influence of Topology on Hydraulic Conductivity in a Sand-and-Gravel Aquifer

A field experiment consisting of geophysical logging and tracer testing was conducted in a single well that penetrated a sand-and-gravel aquifer at the U.S. Geological Survey Toxic Substances Hydrology research site on Cape Cod, Massachusetts. Geophysical logs and flowmeter/pumping measurements were obtained to estimate vertical profiles of porosity ϕ, hydraulic conductivity K, temperature, and bulk electrical conductivity under background, freshwater conditions. Saline-tracer fluid was then injected into the well for 2 h and its radial migration into the surrounding deposits was monitored by recording an electromagnetic-induction log every 10 min. The field data are analyzed and interpreted primarily through the use of Archie's (1942) law to investigate the role of topological factors such as pore geometry and connectivity, and grain size and packing configuration in regulating fluid flow through these coarse-grained materials. The logs reveal no significant correlation between K and ϕ, and imply that groundwater models that link these two properties may not be useful at this site. Rather, it is the distribution and connectivity of the fluid phase as defined by formation factor F, cementation index m, and tortuosity α that primarily control the hydraulic conductivity. Results show that F correlates well with K, thereby indicating that induction logs provide qualitative information on the distribution of hydraulic conductivity. A comparison of α, which incorporates porosity data, with K produces only a slightly better correlation and further emphasizes the weak influence of the bulk value of ϕ on K.     

Lake Winnipeg: geological setting and sediment seismostratigraphy

Lake Winnipeg, the seventh largest lake in North America, is located at the boundary between the Interior Plains and the Canadian Shield in Manitoba, Canada. Seismic profiles were obtained in Lake Winnipeg on two geoscientific cruises in 1994 and 1996. These data indicate the morphology of the bedrock surface. In most cases, a clear distinction between low relief Paleozoic carbonate rock and high relief Precambrian rock can be made. In northern Lake Winnipeg, the eastern limit of Paleozoic rock is clearly demarcated 30 km west of the previous estimate of its position. In southern Lake Winnipeg, all or most of the Paleozoic sequence terminates at a prominent buried escarpment in the centre of the lake. This indicates that Paleozoic rock on the eastern shore, known from drilling and outcrops, is an outlier. Major moraines are apparent as abrupt, large ridges having a chaotic internal reflection pattern. These include the Pearson Reef Moraine, the George Island Moraine and the offshore extension of The Pas Moraine. Little evidence for extensive or thick till was observed. Instead, fine-grained sediments deposited in glacial Lake Agassiz rest directly on bedrock over most of the lake basin. Hence an episode of erosion to bedrock was associated with glaciation and/or deglaciation. The Agassiz Sequence sediments are well-stratified, drape underlying relief and in some areas are over 100 m thick. In places, stratification in these sediments is disrupted, perhaps by dewatering. Evidence of erosion of Agassiz Sequence sediments by recent currents was observed. The contact between the Agassiz Sequence and the overlying Winnipeg Sequence sediments is a marked angular unconformity. The Agassiz Unconformity indicates up to 10 m of erosion in places. The low-relief character of this unconformity precludes subaerial erosion and the lack of till, moraines, or extensive deformation precludes glacial erosion. Waves appear to be the most likely erosional agent, either in waning Lake Agassiz or early Lake Winnipeg time. Winnipeg Sequence sediments, in places very thin, mantle most of the lakefloor. These sediments were deposited in the present Lake Winnipeg and are faintly stratified to massive and reach about 10 m in thickness in deep water. On the surface of the Winnipeg Sequence, vigorous, episodic currents are thought to contribute to the construction of flow-transverse sand waves as much as 6 m high in a deep, narrow constriction in the lake.     

Holocone paleovegetation and paleohydrology of a prairie pothold in southern Saskatchewan, Canada

The Andrews site represents one of countless prairie potholes found in areas of hummocky moraine on the northern Great Plains. Sediments from a depth of 5.8 to 3.1 m at this 'kettle-fill' site in the Missouri Coteau upland of southern Saskatchewan, Canada, provides a record of vegetation, climate, and hydrologic changes within a small, ca 30 m diameter, closed-drainage basin from ca 10.2 to 5.8 ka. Plant macrofossil analyses of 67 samples, 6 ¹⁴C ages, and stratigraphy were used to identify 5 zones, representing the paleohydrological changes that followed deglaciation in southern Saskatchewan.Results of this study indicate that with the melting of residual stagnant ice a pond (>2 m deep) with abundant aquatic, emergent, and shoreline plants developed in the basin at ca 10.2 ka and persisted until at least ca 8.8 ka. During this time there was a shift in upland vegetation from a white spruce forest (Zone II) to a deciduous parkland at ca 10 ka (Zone III). As climate warmed, brackish and alkaline conditions developed coincident with shallowing of the pond at the end of Zone III. The perennial water phase ended at ca 8.8 ka and was followed by a low-water stand lasting ca 1100 years. Prairie fires and slopewash from unstable slopes were dominant sedimentological processes until ca 7.7 ka (Zone IV). Water levels began to rise and between ca 7.7 and 5.8 ka a semi-permanent pond was established in a grassland setting (Zone V). After ca 5.8 ka this prairie pothole wetland became ephemeral, to the point that plant macrofossils could not be preserved, a situation continuing today. Interactions between climate change, variability in local groundwater supply, and sedimentological processes likely account for the paleohydrologic events reconstructed at the Andrews site.     

Holocene lake-level fluctuations and beach-ridge development along the northern shore of Lake Michigan, USA

Geomorphology of a beach-ridge complex and adjacent lake basins along the northern shore of Lake Michigan records fluctuations in the level of Lake Michigan for the last 8000 to 10 000 ¹⁴C yr B.P. (radiocarbon years Before Present). A storm berm at 204.7–206 m (671.6–675.9 ft) exposed in a sandpit provides evidence of a pre-Chippewa Low lake level that is correlated with dropping water levels of Glacial Lake Algonquin (c. 10 300–10 100 ¹⁴C yr B.P.). Radiocarbon dates from organic material exposed in a river cutbank and basal sediments from Elbow Lake, Mackinac Co., Michigan, indicate a maximum age of a highstand of Lake Michigan at 6900 ¹⁴C yr B.P., which reached as high as 196.7 m (645 ft), during the early-Nipissing transgression of Lake Michigan. Basal radiocarbon dates from beach swales and a second lake site (Beaverhouse Lake, Mackinac Co.) provide geomorphic evidence for a subsequent highstand which reached 192.6 m (632 ft) at 5390±70 ¹⁴C yr B.P.Basal radiocarbon dates from a transect of sediment cores, along with tree-ring data, and General Land Office Surveyor notes of a shipwreck, c. A.D. 1846, reveal a late-Holocene rate for isostatic rebound of 22.6 cm/100 radiocarbon years (0.74 ft/100 radiocarbon years) for the northern shore of Lake Michigan, relative to the Lake Michigan-Lake Huron outlet at Port Huron, Michigan. Changes in sediment stratigraphy, inter-ridge distance, and sediment accumulation rates document a mid- to late-Holocene retreat of the shoreline due to isostatic rebound. This regression sequence was punctuated by brief, periodic highstands, resulting in progressive development over the past 5400 ¹⁴C yr of 75 pairs of dune ridges and swales each formed over an interval of approximately 72 years. Times of lake-level fluctuation were identified at 3900, 3200, and 1000 ¹⁴C yr B.P. based on changes in inter-ridge spacing, shifts in the course of Millecoquins River, and reorientation of beach-ridge lineation. Soil type, dune development, and selected pollen data provide supporting evidence for this chronology. Late-Holocene beach-ridge development and lake-level fluctuations are related to a retreat of the dominant Pacific airmass and the convergence of the Arctic and Tropical airmasses resulting in predominantly meridional rather than zonal air flow across the Great Lakes region.This is the 13th in a series of papers published in this special AMQUA issue. These papers were presented at the 1994 meeting of the American Quaternary Association held 19–22 June, 1994, at the University of Minnesota, Minneapolis, Minnesota, USA. Dr Linda C. K. Shane served as guest editor for these papers.