Eruption episodes and magma recharge events in andesitic systems: Mt Taranaki,New Zealand

Eruption episodes, where a series of eruption events are generically related, can include the eruption of a wide spectrum of volcanic activity over decadal periods. This paper concentrates on the opening phases of an eruption episode which occurred approximately 1800 yrs BP from Mt Taranaki, New Zealand. These events spanned the eruption of differing bulk compositions and styles from two distinct vent locations; an andesitic sub-plinian eruption from the summit vent and a scoria cone-building eruption of basaltic magma from a satellite vent. Compositional profiles and zoning textures of plagioclase, amphibole and clinopyroxene phenocrysts from the opening andesitic event show evidence of magma mixing and subsequent crystallisation just prior to the initiation of the eruption episode. Titanomagnetite grain morphology and Ti variation suggest that the magma mixing event occurred within a few days to weeks before the eruption acting as a trigger for it. We present a magmatic model which is constrained by the petrological observations and eruptions of the episode. In this model magma differentiation at depth causes its rise and recharging of a mid-crustal magma storage area at 5–7 km. Although the recharging magma differed slightly in oxygen fugacity and temperature, it was compositionally and physically similar enough to the residing andesitic magma to allow efficient mixing. The petrological characteristics described here can be readily observed and enable identification of mixing events in other recent eruption episodes.     

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.     

Spatial Variation of Sea Level Trend Along the Bangladesh Coast

The Bangladesh coast is threatened by rising sea level due to various factors. The results based on the analysis of past 22 years of tidal data of the Bangladesh coast reveal that the annual mean tidal level in the eastern Bangladesh coast is rising at an alarmingly high rate of 7.8 mm/year, which is almost twice the observed rate in the western region. This type of sea level trend seems to be the result of changing local conditions like increased precipitation and land subsidence during the recent decades. It seems that the higher rate of land subsidence in the eastern Bangladesh coast is the main causative factor for the steeper sea level trends there. The differential sea level trends show that the subsidence component in the sea level rise may be as high as 4 mm/year in the eastern Bangladesh coast. However, this needs to be verified with actual geological observations.     

Recycling in the Puerto Rican mantle wedge, Greater Antilles Island Arc

Abstract The initial volcanic phase of Cretaceous island arc strata in central Puerto Rico, at the eastern end of the extinct Greater Antilles Arc, comprises a 6‐km thick pile of lava and volcanic breccia (Río Majada Group). Preserved within the sequence is a conspicuous shift in absolute abundances of the more incompatible elements, including Th, Nb, and the light rare earth elements (LREE: La, Ce, Pr and Nd). The compositional shift is marked by a decrease in La/Sm from averages of 2.11 in the lowest third of the pile (Formation A) to 1.48 at the top (Formation C), and by a distinctive flattening of LREE segments of chondrite‐normalized REE patterns. i⁸⁷Sr/⁸⁶Sr and ?_Nd average about 0.7035 and 8.2, respectively, in early Formation A basalts. These ranges normally overlap samples from later Formations B and C. Isotope compositions of the latter group are more variable, however, and several samples are considerably more radiogenic than Formation A basalts, such that i⁸⁷Sr/⁸⁶Sr averages almost 0.7042 while ?_Nd‐values decrease to 7.5 in Formation B and C basalts. Theoretical models of non‐modal melting processes in both amphibole peridotite and spinel lherzolite sources provide insight into the origin of depleted Th, Nb, and LREE abundances in Puerto Rican basalts. Low Nb concentrations less than normal mid‐oceanic ridge basalts in Formation A basalts indicate the wedge was slightly depleted by low‐volume decompression fusion due to induced convection in the back‐arc region prior to entry of the source into the arc melting zone. However, depleted patterns in Formation C basalts cannot be generated by relatively greater degrees of decompression fusion in the back‐arc, because addition of the La‐enriched slab‐derived component to more depleted source material invariably produces elevated rather than decreased La/Sm. Refluxing of Formation A harzburgitic residua is similarly precluded. In contrast, the observed patterns are readily reproduced by multistage melting models involving hybridized sources containing normal Formation A lherzolite source material blended with recycled, unrefluxed harzburgite residua. Successful models require hybrid sources containing large volumes of recycled harzburgite (up to 50%) during generation of Formation C basalts. Slightly elevated radiometric Sr and Nd isotopes in a few flows from Formation C are attributed to partial refluxing of the hybrid sources within the wedge.     

Investigation of Failure Mechanisms and Migration of Organic Chemicals at Wilsonville, Illinois

Ground water contamination was discovered in 1981 in a monitoring well at the Earthline disposal facility near Wilsonville, Illinois. Organic chemicals had migrated at a rate 100 to 1000 times greater than predicted when the site received its permit to operate in 1978. Postulated failure mechanisms included migration through previously unmapped permeable zones, subsidence of an underground mine, organic-chemical and clay-mineral interactions, acid-mine drainage and clay interactions, trench-cover settlement, and erosion.
In this investigation, the Illinois State Geological Survey found the primary reason for the rapid migration: the presence of previously undetermined fractures and joints in glacial till. The inaccurate predictions of hydraulic conductivity were based on laboratory-determined values that did not adequately measure the effects of fractures and joints on the transit time calculations. Field-measured hydraulic conductivity values were generally 10 to 1000 times greater than their laboratory-measured counterparts, thus largely accounting for the discrepancy between predicted and actual migration rates in the transit time calculations. The problem was compounded, however, by the burial of liquid wastes and by trench covers that allowed excess surface runoff to enter the trenches. Organic-chemical and clay-mineral interactions may also have exacerbated the problem in areas where liquid organic wastes were buried.     

Oxygen isotope speedometry in granulite facies garnet recording fluid/melt–rock interaction (Sør Rondane Mountains,East Antarctica)

In situ analysis of a garnet porphyroblast from a granulite facies gneiss from Sør Rondane Mountains, East Antarctica, reveals discontinuous step‐wise zoning in phosphorus and large δ¹⁸O variations from the phosphorus‐rich core to the phosphorus‐poor rim. The gradually decreasing profile of oxygen isotope from the core (δ¹⁸O = ~15‰) to the rim (δ¹⁸O = ~11‰) suggests that the ¹⁸O/¹⁶O zoning was originally step‐wise, and modified by diffusion after the garnet rim formation at ~800°C and 0.8 GPa. Fitting of the ¹⁸O/¹⁶O data to the diffusion equation constrains a duration of the high‐T event (~800°C) to c. 0.5–40 Ma after the garnet rim formation. The low δ¹⁸O value of the garnet rim, together with the previously reported low δ¹⁸O values in metacarbonates, indicates regional infiltration, probably along a detachment fault, of low δ¹⁸O fluid/melt possibly derived from meta‐mafic to ultramafic rocks.     

Structure and metamorphism of the Calliope Volcanic Assemblage: Implications for middle to Late Devonian orogeny in the northern New England Fold Belt

The Late Silurian to Middle Devonian Calliope Volcanic Assemblage in the Rockhampton region is deformed into a set of northwest‐trending gently plunging folds with steep axial plane cleavage. Folds become tighter and cleavage intensifies towards the bounding Yarrol Fault to the east. These folds and associated cleavage also deformed Carboniferous and Permian rocks, and the age of this deformation is Middle to Late Permian (Hunter‐Bowen Orogeny). In the Stanage Bay area, both the Calliope Volcanic Assemblage and younger strata generally have one cleavage, although here it strikes north to northeast. This cleavage is also considered to be of Hunter‐Bowen age. Metamorphic grade in the Calliope Volcanic Assemblage ranges from prehnite‐pumpellyite to greenschist facies, with higher grades in the more strongly cleaved rocks. In the Rockhampton region the Calliope Volcanic Assemblage is part of a west‐vergent fold and thrust belt, the Yarrol Fault representing a major thrust within this system.

A Late Devonian unconformity followed minor folding of the Calliope Volcanic Assemblage, but no cleavage was formed. The unconformity does not represent a collision between an exotic island arc and continental Australia as previously suggested.     

Effects of nonlinear model response on allocation of streamflow depletion: exemplified by the case of Beaver Creek,USA

Anomalies found when apportioning responsibility for streamflow depletion are examined. The anomalies arise when responsibility is assigned to the two states that contribute to depletion of Beaver Creek in the Republican River Basin in the United States. The apportioning procedure for this basin presumes that the sum of streamflow depletions, computed by comparing simulation model runs with and without groundwater pumping from individual states, approximates the streamflow depletion when both states are pumping. In the case study presented here, this presumed superposition fails dramatically. The stream drying and aquifer-storage depletion, as represented in the simulation model used for allocation, are examined in detail to understand the hydrologic and numerical basis for the severe nonlinear response. Users of apportioning procedures that rely on superposition should be aware of the presence and likely magnitude of nonlinear responses in modeling tools.