Crustal xenoliths from Calbuco Volcano,Andean Southern Volcanic Zone: implications for crustal composition and magma-crust interaction

Crustal xenoliths in the 1961 andesite flow of Calbuco Volcano, in the southern Southern Volcanic Zone (SSVZ) of the Andes, consist predominantly of pyroxene granulites and hornblende gabbronorites. The granulites contain plagioclase+pyroxene+magnetite±amphibole, and have pronounced granoblastic textures. Small amounts of relict amphibole surrounded by pyroxene-plagioclase-magnetite-glass symplectites are found in some specimens. These and similar textures in the gabbronorites are interpreted as evidence of dehydration melting. Mineral and bulk rock geochemical data indicate that the granulites are derived from an incompatible trace element depleted basaltic protolith that underwent two stages of metamorphism: a moderate pressure, high temperature stage accompanied by melting and melt extraction from some samples, followed by thermal metamorphism after entrainment in the Calbuco andesite lavas. High _Nd ^T values (+4.0 to +8.6), Nd-isotope model ages of 1.7–2.0 Ga, and trace element characteristics like chondrite normalized La/Yb< and La/Nb1 indicate that the protoliths were oceanic basalts. Similar oceanic metabasalts of greenschist to amphibolite facies are found in the Paleozoic metamorphic belt that underlies the Chilean coastal ranges. Mineral and bulk rock compositions of the gabbronorite xenoliths indicate that they are cognate, crystallizing from the basaltic andesite magma at Calbuco. Crystallization pressures for the gabbros based on total Al contents in amphibole are 6–8 kbar. These pressures point to middle to lower crustal storage of the Calbuco magma. Neither granulite nor gabbro xenoliths have the appropriate geochemical characteristics to be contaminants of Calbuco andesites, although an ancient sedimentary contaminant is indicated by the lava compositions. The presence of oceanic metabasaltic xenoliths, together with the sedimentary isotopic imprint, suggests that the lower crust beneath the volcano is analogous to the coastal metamorphic belt, which is an accretionary complex of intercalated basalts and sediments that formed along the Paleozoic Gondwanan margin. If this is the case, the geochemical composition of the lower and middle crust beneath the SSVZ is significantly different from that of most recent SSVZ volcanic rocks.     

Properties of the clouds of Venus,as inferred from airborne observations of its near-infrared reflectivity spectrum

We have measured the shape and absolute value of Venus' reflectivity spectrum in the 1.2-to 4.0-μm spectral region with a circular variable filter wheel spectrometer having a spectral resolution of 1.5%. The instrument package was mounted on the 91-cm telescope of NASA Ames Kuiper Airborne Observatory, and the measurements were obtained at an altitude of about 41,000 feet, when Venus had a phase angle of 86°. Comparing these spectra with synthetic spectra generated with a multiple-scattering computer code, we infer a number of properties of the Venus clouds. We obtain strong confirmatory evidence that the clouds are made of a water solution of sulfuric acid in their top unit optical depth and find that the clouds are made of this material down to an optical depth of at least 25. In addition, we determine that the acid concentration is 84 ± 2% H_2SO4 by weight in the top unit optical depth, that the total optical depth of the clouds is 37.5 ± 12.5, and that the cross-sectional weighted mean particle radius lies between 0.5 and 1.4 μm in the top unit optical depth of the clouds. These results have been combined with a recent determination of the location of the clouds' bottom boundary [Marov et al., Cosmic Res.14, 637–642 (1976)] to infer additional properties about Venus' atmosphere. We find that the average volume mixing ratio of H₂SO₄ and H₂O contained in the cloud material both equal approximately 2× 10^?6. Employing vapor pressure arguments, we show that the acid concentration equals 84 ± 6% at the cloud bottom and that the water vapor mixing ratio beneath the clouds lies between 6 × 10^?4 and 10^?2.     

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

The   M _w γ 9.0  2004 December 26 Sumatra-Andaman and   M _w = 8.7  2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well.     

A highlight of environmental and engineering geology in Fargo, North Dakota, USA

Several unfavorable environmental and engineering geologic conditions exist in Fargo, North Dakota. Dominantly, the behavior of smectitic clays within the proglacial Lake Agassiz sediments of the Sherack and Brenna Formations creates subsoil instability beneath engineered structures in the Fargo area and slope instability within cutbank meanders of the Red River of the North. Unfavorable engineering geologic conditions encountered include: the elastic deformation of clayey glaciolacustrine soils, shrink-swell properties, inadequate bearing capacities, and mass movements. These conditions are responsible for structural failures including the Fargo Grain Elevator in 1955 and the Northern Pacific railroad grade. Bank failures along the Red River are common due to the inherent instability of Brenna Formation smectitic clays which are subject to plastic deformation in the subsurface, with resultant block failure of overlying Sherack Formation. Recent alluvial sediments due to typical fluvial action and the continued seasonal saturation of cutbank meanders within the floodplain also add to soil instability.     

Implications of the relationship between catchment vegetation type and the variability of annual runoff

The impact of changing catchment vegetation type on mean annual runoff has been known for some time, however, the impact on the variability of annual runoff has been established only recently. Differences in annual actual evapotranspiration between vegetation types and the potential effect of changing vegetation type on mean annual runoff and the variability of annual runoff are briefly reviewed. The magnitude of any change in the variability of annual runoff owing to a change in catchment vegetation type is related to the pre‐ and post‐change vegetation types and the seasonality of precipitation, assuming that the variability of annual precipitation remains constant throughout. Significant implications of the relationship between vegetation type and the variability of annual runoff are presented and discussed for water resource management, stream ecology and fluvial geomorphology. Copyright © 2002 John Wiley & Sons, Ltd.     

A Corrected and Generalized Successive Random Additions Algorithm for Simulating Fractional Levy Motions

Simulation of subsurface heterogeneity is important for modeling subsurface flow and transport processes. Previous studies have indicated that subsurface property variations can often be characterized by fractional Brownian motion (fBm) or (truncated) fractional Levy motion (fLm). Because Levy-stable distributions have many novel and often unfamiliar properties, studies on generating fLm distributions are rare in the literature. In this study, we generalize a relatively simple and computationally efficient successive random additions (SRA) algorithm, originally developed for generating Gaussian fractals, to simulate fLm distributions. We also propose an additional important step in response to continued observations that the traditional SRA algorithm often generates fractal distributions having poor scaling and correlation properties. Finally, the generalized and modified SRA algorithm is validated through numerical tests.     

Crustal and upper mantle seismic structure of the Australian Plate, South Island, New Zealand

Seismic reflection and refraction data were collected west of New Zealand's South Island parallel to the Pacific–Australian Plate boundary. The obliquely convergent plate boundary is marked at the surface by the Alpine Fault, which juxtaposes continental crust of each plate. The data are used to study the crustal and uppermost mantle structure and provide a link between other seismic transects which cross the plate boundary. Arrival times of wide-angle reflected and refracted events from 13 recording stations are used to construct a 380-km long crustal velocity model. The model shows that, beneath a 2–4-km thick sedimentary veneer, the crust consists of two layers. The upper layer velocities increase from 5.4–5.9 km/s at the top of the layer to 6.3 km/s at the base of the layer. The base of the layer is mainly about 20 km deep but deepens to 25 km at its southern end. The lower layer velocities range from 6.3 to 7.1 km/s, and are commonly around 6.5 km/s at the top of the layer and 6.7 km/s at the base. Beneath the lower layer, the model has velocities of 8.2–8.5 km/s, typical of mantle material. The Mohorovicic discontinuity (Moho) therefore lies at the base of the second layer. It is at a depth of around 30 km but shallows over the south–central third of the profile to about 26 km, possibly associated with a southwest dipping detachment fault. The high, variable sub-Moho velocities of 8.2 km/s to 8.5 km/s are inferred to result from strong upper mantle anisotropy. Multichannel seismic reflection data cover about 220 km of the southern part of the modelled section. Beneath the well-layered Oligocene to recent sedimentary section, the crustal section is broadly divided into two zones, which correspond to the two layers of the velocity model. The upper layer (down to about 7–9 s two-way travel time) has few reflections. The lower layer (down to about 11 s two-way time) contains many strong, subparallel reflections. The base of this reflective zone is the Moho. Bi-vergent dipping reflective zones within this lower crustal layer are interpreted as interwedging structures common in areas of crustal shortening. These structures and the strong northeast dipping reflections beneath the Moho towards the north end of the (MCS) line are interpreted to be caused by Paleozoic north-dipping subduction and terrane collision at the margin of Gondwana. Deeper mantle reflections with variable dip are observed on the wide-angle gathers. Travel-time modelling of these events by ray-tracing through the established velocity model indicates depths of 50–110 km for these events. They show little coherence in dip and may be caused side-swipe from the adjacent crustal root under the Southern Alps or from the upper mantle density anomalies inferred from teleseismic data under the crustal root.     

The Okavango giant mafic dyke swarm (NE Botswana): its structural significance within the Karoo Large Igneous Province

The structural organization of a giant mafic dyke swarm, the Okavango complex, in the northern Karoo Large Igneous Province (LIP) of NE Botswana is detailed. This N110°E-oriented dyke swarm extends for 1500 km with a maximum width of 100 km through Archaean basement terranes and Permo-Jurassic sedimentary sequences. The cornerstone of the study is the quantitative analysis of N>170 (exposed) and N>420 (detected by ground magnetics) dykes evidenced on a ca. 80-km-long section lying in crystalline host-rocks, at high-angle to the densest zone of the swarm (Shashe area). Individual dykes are generally sub-vertical and parallel to the entire swarm. Statistical analysis of width data indicates anomalous dyke frequency (few data <5.0 m) and mean dyke thickness (high value of 17 m) with respect to values classically obtained from other giant swarms. Variations of mean dyke thicknesses from 17 (N110°E swarm) to 27 m (adjoining and coeval N70°E giant swarm) are assigned to the conditions hosting fracture networks dilated as either shear or pure extensional structures, respectively, in response to an inferred NNW–SSE extension. Both fracture patterns are regarded as inherited brittle basement fabrics associated with a previous (Proterozoic) dyking event. The Okavango N110°E dyke swarm is thus a polyphase intrusive system in which total dilation caused by Karoo dykes (estimated frequency of 87%) is 12.2% (6315 m of cumulative dyke width) throughout the 52-km-long projected Shashe section. Assuming that Karoo mafic dyke swarms in NE Botswana follow inherited Proterozoic fractures, as similarly applied for most of the nearly synchronous giant dyke complexes converging towards the Nuanetsi area, leads us to consider that the resulting triple junction-like dyke/fracture pattern is not a definitive proof for a deep mantle plume in the Karoo LIP.