Mixing and differentiation in the Oruanui rhyolitic magma, Taupo, New Zealand: evidence from volatiles and trace elements in melt inclusions

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km³ (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H₂O, CO₂ and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H₂O and CO₂ in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H₂O and CO₂ in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

DMS and SO2 at Baring Head, New Zealand: Implications for the Yield of SO2 from DMS

Atmospheric dimethyl sulfide (DMS) and sulfur dioxide (SO₂) concentrations were measured at Baring Head, New Zealandduring February and March 2000. Anti-correlated DMS and SO₂ diurnalcycles, consistent with the photochemical production of SO₂ from DMS, were observed in clean southerly air off the ocean. The data is used to infer a yield of SO₂ from DMS oxidation. The estimated yields are highly dependent on assumptions about the DMS oxidation rate. Fitting the measured data in a photochemical box model using model-generated OH levels and the Hynes et al. (1986) DMS + OH rate constant suggests that theSO₂ yield is 50–100%, similar to current estimates for the tropical Pacific.However, the observed amplitude of the DMS diurnal cycle suggests that the oxidation rate is higher than that used by the model, and therefore, that theSO₂ yield is lower in the range of 20–40%.     

An integrated hydrological model for rain‐induced landslide prediction

This paper describes an extension to the Combined Hydrology And Stability Model (CHASM) to fully include the effects of vegetation and slope plan topography on slope stability. The resultant physically based numerical model is designed to be applied to site‐specific slopes in which a detailed assessment of unsaturated and saturated hydrology is required in relation to vegetation, topography and slope stability. Applications are made to the Hawke's Bay region in New Zealand where shallow‐seated instability is strongly associated with spatial and temporal trends in vegetation cover types, and the Mid‐Levels region in Hong Kong, an area subject to a variety of landslide mechanisms, some of which may be subject to strong topographic control. An improved understanding of process mechanism, afforded by the model, is critical for reliable and appropriate design of slope stabilization and remedial measures. Copyright © 2002 John Wiley & Sons, Ltd.     

Paleoecology of a Middle Wisconsin Deposit from Southern California

Analysis of a buried deposit in the Diamond Valley of southern California has revealed well-preserved pollen, wood, and diatom remains. Accelerator mass spectrometry dates of 41,200±2100 and 41,490±1380 ¹⁴C yr B.P. place this deposit in marine isotope stage 3. Diatoms suggest a shallow lacustrine environment. Pollen data suggest that several plant communities were present near the site, with grassland, scrub, chaparral, forest, and riparian communities represented. Comparison with modern pollen suggests similarities with montane forests in the nearby San Bernardino and San Jacinto ranges, indicating vegetation lowering by at least 900 m elevation and temperatures 4°–5°C cooler than today. An increase in high-elevation conifer pollen documents climatic cooling near the profile top. Early-profile diatoms are typical of warm water with high alkalinity and conductivity, whereas later diatoms suggest a higher flow regime and input of cooler water into the system. We suggest that the sequence is part of the cooling phase of an interstadial Dansgaard–Oeschger cycle. Records of the middle Wisconsin period are rare in southern California, but the Diamond Valley site is similar to records from Tulare Lake in the San Joaquin Valley and the ODP Site 893A record from Santa Barbara Basin. It is probable that the Diamond Valley assemblage is a local expression of a vegetation type widespread in the ranges and basins of southwestern California during the middle Wisconsin.     

Effect of adsorbed metal ions on the transport of Zn- and Ni-EDTA complexes in a sand and gravel aquifer

Adsorption, complexation, and dissolution reactions strongly influenced the transport of metal ions complexed with ethylenediaminetetraacetic acid (EDTA) in a predominantly quartz-sand aquifer during two tracer tests conducted under mildly reducing conditions at pH 5.8 to 6.1. In tracer test M89, EDTA complexes of zinc (Zn) and nickel (Ni), along with excess free EDTA, were injected such that the lower portion of the tracer cloud traveled through a region with adsorbed manganese (Mn) and the upper portion of the tracer cloud traveled through a region with adsorbed Zn. In tracer test S89, Ni- and Zn-EDTA complexes, along with excess EDTA complexed with calcium (Ca), were injected into a region with adsorbed Mn. The only discernable chemical reaction between Ni-EDTA and the sediments was a small degree of reversible adsorption leading to minor retardation. In the absence of adsorbed Zn, the injected Zn was displaced from EDTA complexes by iron(III) [Fe(III)] dissolved from the sediments. Displacement of Zn by Fe(III) on EDTA became increasingly thermodynamically favorable with decreasing total EDTA concentration. The reaction was slow compared to the time-scale of transport. Free EDTA rapidly dissolved aluminum (Al) from the sediments, which was subsequently displaced slowly by Fe. In the portion of tracer cloud M89 that traveled through the region contaminated with adsorbed Zn, little displacement of Zn complexed with EDTA was observed, and Al was rapidly displaced from EDTA by Zn desorbed from the sediments, in agreement with equilibrium calculations. In tracer test S89, desorption of Mn dominated over the more thermodynamically favorable dissolution of Al oxyhydroxides. Comparison with results from M89 suggests that dissolution of Al oxyhydroxides in coatings on these sediment grains by Ca-EDTA was rate-limited whereas that by free EDTA reached equilibrium on the time-scale of transport. Rates of desorption are much faster than rates of dissolution of Fe oxyhydroxides from sediment-grain surfaces and, therefore, adsorbed metal ions can strongly influence the speciation of ligands like EDTA in soils and sediments, especially over small temporal and spatial scales.     

Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.     

A 1500-year record of climatic and environmental change in ElkLake, Minnesota I: Varve thickness and gray-scale density

The deepest part (29.5 m) of Elk Lake, ClearwaterCounty, northwestern Minnesota, contains a complete Holocene section that iscontinuously varved. The varve components are predominantly autochthonous(CaCO₃, organic matter, biogenic silica, and several iron andmanganese minerals), but the varves do contain a minor detrital-clastic(aluminosilicate) component that is predominantly wind-borne (eolian) andprovides an important record of atmospheric conditions. Singular spectrumanalysis (SSA) and wavelet analysis of varve thickness recognized significantperiodicities in the multicentennial and multidecadal bands that varied inpower (i.e., variable significance) and position (i.e., variable period) withinthe periodic bands. Persistent periodicities of about 10, 22, 40, and 90 years,and, in particular, multicentennial periodicities in varve thickness and otherproxy variables are similar to those in spectra of radiocarbon production, aproxy for past solar activity. This suggests that there may be a solar control,perhaps through geomagnetic effects on atmospheric circulation. Multicentennialand multidecadal periodicities also occur in wavelet spectra of relativegray-scale density. However, gray-scale density does not appear to correlatewith any of the measured proxy variables, and at this point we do not know whatcontrolled gray scale.     

Using Spatial Statistics to Examine the Pattern of Urban Land Development in Halifax–Dartmouth

This article examines the locational patterns of population and employment and the resultant changes in urban form in the Halifax–Dartmouth region between 1970 and 1996. The article employs the univariate and bivariate K function to measure spatial dependence or clustering within and between the classes of residential and commercial land parcels. The results of the K function estimates suggest that residential land parcels cluster together, commercial land parcels cluster together, and residential and commercial land parcels have become more clustered over time. Evidence of clustering provides insight into the changing urban form of the region and possible multinucleation.     

Wind erosion from military training lands in the Mojave Desert, California, U.S.A.

Military training activities reduce vegetation cover, disturb crusts, and degrade soil aggregates, making the land more vulnerable to wind erosion. The objective of this study was to quantify wind erosion rates for typical conditions at the Marine Corps Air Ground Combat Center, Twentynine Palms, CA, U.S.A. Five Big Spring Number Eight (BSNE) sampler stations were installed at each of five sites. Each BSNE station consisted of five BSNE samplers with the lowest sampler at 0·05 m and the highest sampler at 1·0 m above the soil surface. Once a month, sediment was collected from the samplers for analysis. Occurrence of saltating soil aggregates was recorded every hour using Sensits, one at each site. The site with the most erosion had a sediment discharge of 311 kg m⁻¹ over a period of 17 months. Other sites eroded much less because of significant rock cover or the presence of a crust. Hourly sediment discharge was estimated combining hourly Sensit count and monthly sediment discharge measured using BSNE samplers. More simultaneously measured data are needed to better characterize the relationship between these two and reconstruct a detailed time-series of wind erosion. This measured time-series can then be used for comparison with simulation results from process-based wind erosion models such as the Wind Erosion Prediction System (WEPS), once it has been adapted to the unique aspects of military lands.