Igneous Petrology, 2nd edn, by Anthony Hall. Addison Wesley, Longman, 1996. 551 pp. {pound}26.99 (paperback).

The second edition of Dr Anthony Hall's undergraduate textbookhas been published some 9 years after the first edition. Thisedition is a significant revision although the approach hasnot changed, and many of the comments of this journal's reviewerof the first edition are still appropriate (Macdonald, 1988).The structure of this edition is largely the same as     

Compressibility and crystal structure of sillimanite, Al2SiO5, at high pressure

The unit-cell dimensions and crystal structure of sillimanite at various pressures up to 5.29 GPa have been refined from single-crystal X-ray diffraction data. As pressure increases, a and b decrease linearly, whereas c decreases nonlinearly with a slightly positive curvature. The axial compression ratios at room pressure are β_a:β_b:β_c=1.22:1.63:1.00. Sillimanite exhibits the least compressibility along c, but the least thermal expansivity along a (Skinner et al. 1961; Winter and Ghose 1979). The bulk modulus of sillimanite is 171(1) GPa with K′=4 (3), larger than that of andalusite (151 GPa), but smaller than that of kyanite (193 GPa). The bulk moduli of the [Al1O₆], [Al2O₄], and [SiO₄] polyhedra are 162(8), 269(33), and 367(89) GPa, respectively. Comparison of high-pressure data for Al₂SiO₅ polymorphs reveals that the [SiO₄] tetrahedra are the most rigid units in all these polymorphic structures, whereas the [AlO₆] octahedra are most compressible. Furthermore, [AlO₆] octahedral compressibilities decrease from kyanite to sillimanite, to andalusite, the same order as their bulk moduli, suggesting that [AlO₆] octahedra control the compression of the Al₂SiO₅ polymorphs. The compression of the [Al1O₆] octahedron in sillimanite is anisotropic with the longest Al1-OD bond shortening by ～1.9% between room pressure and 5.29 GPa and the shortest Al1-OB bond by only 0.3%. The compression anisotropy of sillimanite is primarily a consequence of its topological anisotropy, coupled with the compression anisotropy of the Al-O bonds within the [Al1O₆] octahedron.     

A fast response,open path,infrared hygrometer,using a semiconductor source

Investigations were carried out as to the feasibility of using a semiconductor source in the design of a new rapid response, open-path hygrometer. A single-beam instrument was constructed employing an infrared light emitting diode (LED) as a source instead of the usual high energy, wideband filament. The spectral emission envelope encompassed the 1.87 m water absorption band. Electronic modulation and thermoelectrical cooling of the diode eliminated the conventional chopper wheel and stabilized the peak wavelength emission. Path length was 200 mm. Over a water vapour concentration range of 0–16 g m^–3, absorption varied by 2% in a linear fashion. At 10 Hz, the noise level was 0.1 g m^–3 rms. Hygrometer resolution and stability are constrained by the detector noise level, the small source emission in the absorption band and low frequency drift in the optical filter. Despite these problems, the new instrument showed comparable performance characteristics to a commercial Lyman- hygrometer. Latent heat fluxes measured with both instruments and a Kaijo-Denki, 3-D sonic anemometer agreed to within 4% over a range 0–350 W m^–2. Further improvements in performance can be anticipated with advances in detector and LED technology.     

A multifractal analysis of lidar measured water vapour

Accurate and fast-response measurements of space-time observations of specific humidity were made above a drying land surface at the University of California at Davis, using the Los Alamos water Raman-lidar. In an attempt to quantify the space-time intermittency features of turbulent flows in the lower atmosphere, a multifractal analysis of these water vapour measurements was performed. The structure of the specific humidity, (x, t), was analyzed quantifying a scalar gradient measure both in time and space, for all possible one-dimensional cuts, i.e. and . The results confirm the multifractal nature of this scalar gradient measure (a type of scalar dissipation rate) and show that humidity measurements at fixed times (x) are more intermittent (e.g. have less entropy dimension) than those at fixed locations in space (t). Similar multifractal behaviour of the spatial data, with and without a transformation from the observed wind velocities, supports the validity of Taylor's hypothesis for the studied fields.     

Geology of the volcanic-hosted Brockman rare-metals deposit,Halls Creek Mobile Zone,northwest Australia. l. Volcanic environment,geochronology and petrography of the Brockman volcanics

Summary Rare-metals mineralization of the Brockman deposit (Halls Creek Mobile Zone, NW Australia) is hosted in a fluorite-bearing, rhyolitic volcaniclastic unit informally termed the Niobium Tuff. The Tuff, more correctly described as a tuffaceous volcaniclastic deposit, is the lowermost unit of a sequence of trachyte-to-rhyolite lavas, trachyandesite subvolcanic rocks, and volcaniclastic units of the Brockman volcanics located within the Halls Creek Group, a thick, early Proterozoic volcano-sedimentary sequence. High precision SHRIMP ion-microprobe zircon dating of the Niobium Tuff gives an eruption age of 1870 ± 4 Ma. Regional geochronological constraints indicate the Niobium Tuff was deposited about 15 Ma before major orogenic activity affected the area. Despite folding, faulting and low-grade metamorphism, the Brockman volcanics show excellent preservation of primary volcanological features, including pillow-lavas and pillow-breccias, that suggest a dominantly subaqueous, below-wave-base emplacement environment. The style of eruption products and magma volume constraints suggest the trachyte-dominated volcanics were erupted from a small shield volcanic complex probably in a rift-related basin in a shallow-marine setting. The tectonic setting was intraplate but differs from most Tertiary to recent trachyte volcanic complexes which are largely subaerial, are built on relatively thick continental crust, and show no post-eruptive orogenic history. Brockman-style rare-metal deposits are characterized by preservation of subaqueous volcanics beneath a thick sedimentary sequence, eruption of early incompatible-element enriched products followed by less differentiated magmas, and fine-grained mineralogy influenced by alteration processes. Prospects exist for discovery of analogous deposits, particularly in early Proterozoic mobile belts and Tertiary intraplate shield volcanic provinces.

---

Geologie der an Vulkanite gebundenen Seltene-Metalle-Lagerstätte Brockman, Halls Creek Mobile Zone, Nordwest-Australien. I. Vulkanologie, Geochronologie und Petrographie der Brockman-Vulkanite

Zusammenfassung Selten-Metall-Vererzung der Brockman-Lagerstätten (Halls Creek Mobile Zone, NW Australien) sitzt in einer Fluorit-führenden, rhyolitischen vulkanoklastischen Einheit auf, die informell als der Niob-Tuff bezeichnet wird. Es handelt sich hier um ein vulkanoklastisches Tuff-Sediment; dieses ist die unterste Einheit einer Abfolge vor Trachyt bis Rhyolit-Laven, trachyandesitischen Subvulkaniten and vulkanoklastischen Einheiten der Brockman Vulkanite innerhalb der Halls Creek Gruppe, einer mächtigen frühproterozoischen vulkanosedimentdren Abfolge. Prüzisionsdatierungen des Niob Tuffs mit der SHRIMP Ionen-Mikrosonde ergeben ein Eruptionsalter von 1870 +-4 Millionen Jahren. Regionale geochronologische Zusammenhänge zeigen, daß der Niob Tuff etwa 15 Millionen Jahre vor einer größeren Orogenese, die das Gebiet betroffen hat, abgelagert wurde. Trotz Faltung, Bruchtektonik and niedriggradiger Metamorphose zeigen die Brockman-Vulkanite einen hervorragenden Erhaltungszustand primärer vulkanologischer Erscheinungen. Diese urnfassen auch Kissenlaven und KissenBreckzien, die eine vorwiegend subaquatische Ablagerung in ruhigem Wasser erkennen lassen. Die Art der Eruptionsprodukte and das Magmavolumen zeigen, daß die Trachytdominierten Vulkanite von einem kleinen Schildvulkan stammer, wahrscheinlich in einem Becken in einer Rift-Situation im seichten marinen Milieu. Die tektonische Situation war intraplate, aber unterscheidet sich von den meisten tertiären bis rezenten trachytischen Vulkan-Komplexen, die hauptsächlich subaerisch sind, auf einer relativ mächtigen kontinentalen Kruste aufsitzen, and keine post-eruptive orogene Entwicklung zeigen. Seltene-Metalle-Lagerstatten des Brockman-Typs rind durch die Erhaltung subaquatischer-Vulkanite unterhalb einer machtigen sedimentdren Abfolge gekennzeichnet; welters durch frühe Eruptionsprodukte, die an inkompatiblen Elementen angereichert sind, auf die dann weniger differenzierte Magmen folgten, und schließlich durch einen feinkörnigen Mineralbestand, der vor Umwandlungsprozessen betroffen war. Es besteht die Möglichkeit der Entdeckung analoger Lagerstätten, besonders in frühproterozoischen mobilen Gürteln und in tertidren intraplate Schildvulkan Provinzen.

---

---

With 5 Figures

---

Deceased     

Textures exhibited by feldspars in the Giant Mottled Anorthosite (GMA) of the bastard unit in the Upper Critical Zone,Western Bushveld Complex

Summary The Bastard Unit, overlying the Merensky Unit, is the uppermost and most complete of the cyclic units in the Upper Critical Zone and its upper contact is taken as the boundary between the Critical and Main Zones (as recommended by SACS, 1980). The Unit is different from underlying units for three reasons: it does not have a well-defined top contact (with an overlying unit); it is much thicker than the Merensky and Footwall Units; and it has an unusually thick sequence of mottled anorthosite at its top. This distinctive composite anorthositic sequence is termed the Giant Mottled Anorthosite (GMA) which is predominantly made up of poikilitic anorthosite. The GMA can broadly be subdivided into three distinct parts; the Lower Giant Mottled Anorthosite (LGMA), Giant Mottled Middling (GMM) and the Upper Giant Mottled Anorthosite (UGMA).Petrographic examination of 82 samples taken in the GMA from 11 borehole profiles around the western limb of the Complex revealed that the LGMA and GMM are essentially adcumulates while the UGMA is orthocumulate in character. By its very nature the dominant phase throughout the GMA is cumulus plagioclase feldspar ( 85%) with the balance being made up of intercumulus orthopyroxene and clinopyroxene and minor biofite — all of which constitute the darker mottles. Intercumulus inverted pigeonite occurs in the upper part of the GMM, throughout the UGMA and for c. 30 m into the Main Zone. K-feldspar also joins the paragenesis within the UGMA and occurs both as intercumulus orthoclase and antiperthite. Cumulus plagioclase grains, showing varying degrees of complex oscillatory-zonal growth are common in the upper part of the GMA. These grains are interpreted as records of growth episodes during earlier influxes of magma which were responsible for the development of the underlying cyclic units, but escaped incorporation into the cumulate pile.

---

Feldspat-Texturen im Giant Mottled Anorthosite der Bastard-Einheit in der oberen kritischen Zone des westlichen Bushveld-Komplexes

Zusammenfassung Die Bastard-Einheit, die die Merensky-Einheit überlagert, ist die oberste und am besten entwickelte der zyklischen Einheiten in der oberen kritischen Zone; ihr oberer Kontakt gilt als die Grenze zwischen der Kritischen und der Main-Zone (entsprechend den Empfehlungen von SACS, 1980). Die Einheit unterscheidet sich aus drei Gründen von den Einheiten im liegenden: sie hat keinen gut definierten oberen Kontakt (mit einer darüberliegenden Einheit), sie ist wesentlich mächtiger als die Merensky und die Footwall-Einheiten, und sie hat eine ungewöhnliche mächtige Abfolge von Mottled Anorthosite als obersten Teil der Abfolge. Diese wohldefinierte anorthositische Abfolge wird als Giant Mottled Anorthosite (GMA) bezeichnet, und besteht hauptsächlich aus poikilitischem Anorthosit. Die GMA kann in 3 Teile, den unteren (LGMA), den mittleren (GMM) un den oberen (UGMA) unterteilt werden.Die petrographische Untersuchung von 82 Proben in der GMA aus 11 Bohrlöchern im Westteil des Bushveld-Komplexes zeigt, daß LGMA und GMM im wesentlichen Adkumulate sind, während die UGMA ein Orthokumulat darstellt. Die Hauptphase in der gesamten GMA ist Kumulus-Plagioklas ( 85 %); der Rest besteht aus Intercumulus-Orthopyroxen, Klinopyroxen, und geringen Mengen von Biotit. Diese sind die Komponenten der dunkleren mottles. Intercumulus invertierter Pigeonit kommt im oberen Teil der GMM, in der gesamten UGMA und in den unteren 30 Metern der Main-Zone vor. Innerhalb der UGMA tritt K-Feldspat hinzu, dieser kommt als Intercumulus Orthoklas und Antiperthit vor. Körner von Cumulusplagioklas mit verschiedenen Intensitäten komplexer Zonierung sind im oberen Teil der GMA verbreitet. Diese Körner werden als Hinweise auf Wachstumsepisoden während früherer Magmenzufuhr-Phasen interpretiert, die für die Entwicklung der zyklischen Einheiten im Liegenden verantwortlich waren, aber selbst nicht Teil der Cumulatabfolge wurden.

---

---

With 5 Figures     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

Crystallization processes in an artificial magma: variations in crystal shape,growth rate and composition with melt cooling history

A large (4.8 m³, 1.3x10⁷ g) artificial mafic melt with a bulk composition similar to that of a basalt (but with a high CaO content of 17 wt%) was generated during a demonstration of in situ vitrification and was allowed to cool naturally. During the melting process, convection was vigorous, resulting in a chemically and thermally homogeneous melt body. Once heating was complete, the cooling rate was rapid with the temperature dropping from 1500°C to 500°C in 6 days within the interior of the 3 m diameter, 1.5 m thick body. A 20 h period of constant temperature (1140°C) observed during colling was the result of latent heat released by widespread crystallization. The final crystalline assemblage consists of diopsidic to hedenbergitic pyroxene and anorthitic feldspar, with a subordinate amount of potassic feldspar, plus a small amount of evolved glass. The compositions and proportions of phases agree well with those predicted by the MELTS thermodynamic model. Thermal and textural evidence suggest that convection within the melt ceased coincident with formation of the first crystals. Textural investigation of core samples reveals large (up to 1 cm in length) acicular diopsidic pyroxenes in a matrix of smaller feldspar and zoned pyroxene crystals (500 m in length). Crystal shape and pyroxene composition vary as a function of position within the solidified body, as a function of cooling rate. Both crystal size and degree of crystallinity are highest in the central, most slowly-cooled parts of the rock. Crystal shape is characterized by tabular, equilibrium-growth forms in the slowly-cooled areas, grading to highly skeletal, dendritic forms at the rapidly-cooled edges of the body. The pyroxene crystals are dominantly homogeneous diopside, but crystals are characterized by thin Fe-rich hedenbergitic rims. These rims were deposited when Mg in the melt was depleted by diopside growth, and melt temperature had cooled sufficiently to allow Fe-rich pyroxene growth. Crystal growth rates can be calculated based on thermal behavior of the melt, reinforced by thermodynamic modelling, and are determined to be between 10^-7 and 10^-8 cm/s in the central part of the melt. These estimates agree well with growth rates in natural systems with similar cooling rates. Pyroxene crystals that formed at a higher cooling rates are characterized by higher Al and lower Mg contents relative to tabular equilibrium crystalline forms, presumably as a result of disequilibrium melt compositions at the crystal-melt interface.New Mexico Bureau of Mines and Mineral Resources, New Mexico Institute of Mining and Technology, Socorro, NM 87801, USA     

Geochemical Evidence for Subduction Fluxes, Mantle Melting and Fractional Crystallization Beneath the South Sandwich Island Arc

The volcanoes of the South Sandwich island arc follow threedistinct series: low-K tholeiitic (followed by Zavodovski, Candlemas,Vindication, Montagu and Bristol), tholeiitic (followed by Visokoi,Saunders and Bellinghausen) and calcalkaline (followed by Leskov,Freehand and part of Cook and Thule). Flux calculations indicatethat the percentage contribution of the subduction componentto the mantle source of all three series varies from undetectable(e.g. Zr) through small (e.g. Nd=20%) and moderate (e.g. La,Ce, Sr=50–80%) to dominant (e.g. Pb, K, Ba, Rb, Cs >90%)with little change along the arc. Isotope systematics (Pb, Nd,Sr) show that this subduction component obtains a greater contributionfrom altered oceanic crust than from pelagic sediment. Elementsfor which the subduction contribution is small show that themantle is already depleted relative to N-MORB mantle (equivalentto loss of an 2•5% melt fraction) before melting beneaththe arc. After addition of the subduction component, dynamicmelting of this depleted mantle then causes the variations inK that distinguish the three series. The estimated degree ofpartial melting (20%) is slightly greater than that beneathocean ridges, though geothermometry suggests that the primarymagma temperature (1225C) is similar to that of primary MORB.About half of the melting may be attributed to volatile addition,and half to decompression. Dynamic melting involving three-dimensional,two-phase flow may be needed to explain fully the inter-islandvariations. KEY WORDS: geochemistry; petrology; fluxes; melting; subduction ^*Corresponding author