ALBERT

Notes: A petrogenetic grid in the model system CaO–FeO–MgO–Al2O3–SiO2–H2O is presented, illustrating the phase relationships among the minerals grunerite, hornblende, garnet, clinopyroxene, chlorite, olivine, anorthite, zoisite and aluminosilicates, with quartz and H2O in excess. The grid was calculated with the computer software thermocalc, using an upgraded version of the internally consistent thermodynamic dataset HP98 and non-ideal mixing activity models for all solid solutions. From this grid, quantitative phase diagrams (P–T pseudosections) are derived and employed to infer a P–T path for grunerite–garnet-bearing amphibolites from the Endora Klippe, part of the Venetia Klippen Complex within the Central Zone of the Limpopo Belt. Agreement between calculated and observed mineral assemblages and garnet zonation indicates that this part of the Central Zone underwent a prograde temperature and pressure increase from c. 540 °C/4.5 kbar to 650 °C/6.5 kbar, followed by a post-peak metamorphic pressure decrease. The inferred P–T path supports a geotectonic model suggesting that the area surrounding the Venetia kimberlite pipes represents the amphibolite-facies roof zone of migmatitic gneisses and granulites that occur widely within the Central Zone. In addition, the P–T path conforms to an interpretation that the Proterozoic evolution of the Central Zone was controlled by horizontal tectonics, causing stacking and differential heating at c. 2.0 Ga.

Notes: A review and reinterpretation of previous experimental data on the deformation of partially melted crustal rocks reveals that the relationship of aggregate strength to melt fraction is non-linear, even if plotted on a linear ordinate and abscissa. At melt fractions, Φ 〈 0.07, the dependence of aggregate strength on Φ is significantly greater than at Φ 〉 0.07. This melt fraction (Φ = 0.07) marks the transition from a significant increase in the proportion of melt-bearing grain boundaries up to this point to a minor increase thereafter. Therefore, we suggest that it is the increase of melt-interconnectivity that causes the dramatic strength drop between the solidus and a melt fraction of 0.07. We term this drop the ‘melt connectivity transition’ (MCT). A second, less-pronounced strength drop occurs at higher melt fractions and corresponds to the breakdown of the solid (crystal) framework. This is the ‘solid-to-liquid transition’ (SLT), corresponding to the well known ‘rheologically critical melt percentage’. Although the strength drop at the SLT is about four orders of magnitude, the absolute value of this drop is small compared with the absolute strength of the unmelted aggregate, rendering the SLT invisible in a linear aggregate strength v. melt-fraction diagram. On the other hand, the more important MCT has been overlooked in previous work because experimental data usually are plotted in logarithmic strength v. melt-fraction diagrams, obscuring large strength drops at high absolute strength values. We propose that crustal-scale localization of deformation effectively coincides with the onset of melting, pre-empting attainment of the SLT in most geological settings. The SLT may be restricted to controlling flow localization within magmatic bodies, especially where melt accumulates.

Notes: Contact metamorphism caused by the Glenmore plug in Ardnamurchan, a magma conduit active for 1 month, resulted in partial melting, with melt now preserved as glass. The pristine nature of much of the aureole provides a natural laboratory in which to investigate the distribution of melt. A simple thermal model, based on the first appearance of melt on quartz–feldspar grain boundaries, the first appearance of quartz paramorphs after tridymite and a plausible magma intrusion temperature, provides a time-scale for melting. The onset of melting on quartz–feldspar grain boundaries was initially rapid, with an almost constant further increase in melt rim thickness at an average rate of 0.5–1.0 × 10−9 cm s−1. This rate was most probably controlled by the distribution of limited amounts of H2O on the grain boundaries and in the melt rims.The melt in the inner parts of the aureole formed an interconnected grain-boundary scale network, and there is evidence for only limited melt movement and segregation. Layer-parallel segregations and cross-cutting veins occur within 0.6 m of the contact, where the melt volume exceeded 40%. The coincidence of the first appearance of these signs of the segregation of melt in parts of the aureole that attained the temperature at which melting in the Qtz–Ab–Or system could occur, suggests that internally generated overpressure consequent to fluid-absent melting was instrumental in the onset of melt movement.

Notes: New eclogite localities and new 40Ar/39Ar ages within the Western Gneiss Region of Norway define three discrete ultrahigh-pressure (UHP) domains that are separated by distinctly lower pressure, eclogite facies rocks. The sizes of the UHP domains range from c. 2500 to 100 km2; if the UHP culminations are part of a continuous sheet at depth, the Western Gneiss Region UHP terrane has minimum dimensions of c. 165 × 50 × 5 km. 40Ar/39Ar mica and K-feldspar ages show that this outcrop pattern is the result of gentle regional-scale folding younger than 380 Ma, and possibly 335 Ma. The UHP and intervening high-pressure (HP) domains are composed of eclogite-bearing orthogneiss basement overlain by eclogite-bearing allochthons. The allochthons are dominated by garnet amphibolite and pelitic schist with minor quartzite, carbonate, calc-silicate, peridotite, and eclogite. Sm/Nd core and rim ages of 992 and 894 Ma from a 15-cm garnet indicate local preservation of Precambrian metamorphism within the allochthons. Metapelites within the allochthons indicate near-isothermal decompression following (U)HP metamorphism: they record upper amphibolite facies recrystallization at 12–17 kbar and c. 750 °C during exhumation from mantle depths, followed by a low-pressure sillimanite + cordierite overprint at c. 5 kbar and c. 750 °C. New 40Ar/39Ar hornblende ages of 402 Ma document that this decompression from eclogite-facies conditions at 410–405 Ma to mid-crustal depths occurred in a few million years. The short timescale and consistently high temperatures imply adiabatic exhumation of a UHP body with minimum dimensions of 20–30 km. 40Ar/39Ar muscovite ages of 397–380 Ma show that this extreme heat advection was followed by rapid cooling (c. 30 °C Myr−1), perhaps because of continued tectonic unroofing.

Notes: High-grade gneisses (amphibolite–granulite facies) of the Namche Barwa and Gyala Peri massifs, in the eastern Himalayan syntaxis, have been unroofed from metamorphic depths in the late Tertiary–Recent. Rapid exhumation (2–5 mm year−1) has resulted in a pronounced shallow conductive thermal anomaly beneath the massifs and the intervening Tsangpo gorge. The position of the 300 °C isotherm has been estimated from fluid inclusions using CO2–H2O immiscibility phase equilibria to be between 2.5 and 6.2 km depth below surface. Hence, the near-surface average thermal gradient exceeds 50 °C km−1 beneath valleys, although the thermal gradient is relatively lower beneath the high mountains. The original metamorphic fluid in the gneisses was 〉90% CO2. This fluid was displaced by incursion of brines from overlying marine sedimentary rocks that have since been largely removed by erosion. Brines can exceed 60 wt% dissolved salts, and include Ca, Na, K and Fe chlorides. These brines were remobilized during the earliest stages of uplift at 〉500 °C. During exhumation, incursion of abundant topography-driven surface waters resulted in widespread fracture-controlled hydrothermal activity and brine dilution down to the brittle–ductile transition. Boiling water was particularly common at shallow levels (〈2.5 km) beneath the Yarlung Tsangpo valley, and numerous hot springs occur at the surface in this valley. Dry steam is not a major feature of the hydrothermal system in the eastern syntaxis (in contrast to the western syntaxis at Nanga Parbat), but some dry steam fluids may have developed locally.

Notes: The Limousin ophiolite is located at the suture zone between two major thrust sheets in the western French Massif Central. This ophiolitic section comprises mantle-harzburgite, mantle-dunite, wehrlites, troctolites and layered gabbros. It has recorded a static metamorphic event transforming the gabbros into undeformed amphibolites and the magmatic ultramafites into serpentinites and/or pargasite-bearing chloritites. With various thermobarometric methods, it is possible to show that the different varieties of amphibole have registered low-P (c. 0.2 GPa) conditions with temperature ranging from high-T, late-magmatic conditions to greenschist–zeolite metamorphic facies. The abundance of undeformed metamorphic rocks (which is typical of the lower oceanic crust), the occurrence of Ca–Al (–Mg) metasomatism illustrated by the growth of Ca–Al silicates in veins or replacing the primary magmatic minerals, the P–T conditions of the metamorphism and the numerous similarities with oceanic crustal rocks from Ocean Drilling Program and worldwide ophiolites are the main arguments for an ocean-floor hydrothermal metamorphism in the vicinity of a palaeo-ridge. Among the West-European Variscan ophiolites, the Limousin ophiolites constitute an extremely rare occurrence that has not been involved in any HP (subduction-related) or MP (orogenic) metamorphism as observed in other ophiolite occurrences (i.e. France, Spain and Germany).

Notes: A structural, metamorphic and geochronological study of the Staré Město belt implies the existence of two distinct metamorphic events of similar peak P–T conditions (700–800 °C, 8–10 kbar) during the Cambro-Ordovician and the Carboniferous tectonometamorphic events. The hypothesis of two distinct periods of metamorphism was suggested on the basis of structural discordance between an undoubtedly Carboniferous granodiorite sill intrusion and earlier Cambro-Ordovician fabrics of a banded amphibolite complex. The analysis of crystal size distribution (CSD) shows high nucleation density (N0) and low average growth rate (Gt) for Carboniferous mylonitic metagabbros and mylonitic granodiorites. The parameter N0 decreases whereas the quantity Gt increases towards higher temperatures progressively approaching the values obtained from the Cambro-Ordovician banded amphibolite complex. The spatial distribution of amphibole and plagioclase shows intense mechanical mixing for lower-temperature mylonitic metagabbros. In high-temperature mylonites a strong aggregate distribution is developed. Cambro-Ordovician amphibolites unaffected by Carboniferous deformation show a regular to anticlustered spatial distribution resulting from heterogeneous nucleation of individual phases. This pattern, together with CSD, was subsequently modified by the grain growth and textural equilibration controlled by diffusive mass transfer during Carboniferous metamorphism. The differences between the observed textures of the amphibolites are interpreted to be a consequence of the different durations of the Carboniferous and Cambro-Ordovician thermal events.

Notes: Serpentinite mylonites from the Happo ultramafic complex show evidence of two stages of mylonitization at different temperature conditions. Peridotite mylonites exhibit two types of olivine – porphyroclasts and neoblasts – produced at the earlier stage. The olivine neoblasts have a stretching lineation with a fabric suggesting plastic deformation along (0 1 0) [0 0 1]. In addition to the olivine fabric, the stable association of olivine, orthopyroxene and tremolite in the peridotites that survived later serpentinization, and the Si and Na contents of tremolite, suggest that the earlier mylonitization took place at temperatures between 700 and 800 °C. Later mylonitization was associated with high-temperature serpentinization to form serpentinite mylonites. In contrast to a common type of serpentinite in orogenic belts, the serpentinite mylonites are cohesively foliated, rich in olivine and diopside, and poor in antigorite. The diopside has low Al, Cr and Na contents typical of a retrograde origin, and the olivine has a homogeneous composition except in areas subjected to contact metamorphism at a later stage. Modal composition and mineral chemistry suggest that the serpentinite mylonites were formed by a hydration reaction of tremolite and olivine to produce diopside and antigorite under stable conditions of olivine, at temperatures between 400 and 600 °C. Later-stage mylonitization has preferentially been superimposed on the earlier-stage mylonite zone with a common direction of foliation. The difference in temperature between the two mylonitization stages suggests that the shear zone was episodically active during the emplacement of the Happo complex. Conditions of relatively high temperature for serpentinization at a convergent plate boundary and high permeability caused by the early mylonitization favoured the formation of the serpentinite mylonites.

Notes: High-pressure kyanite-bearing felsic granulites in the Bashiwake area of the south Altyn Tagh (SAT) subduction–collision complex enclose mafic granulites and garnet peridotite-hosted sapphirine-bearing metabasites. The predominant felsic granulites are garnet + quartz + ternary feldspar (now perthite) rocks containing kyanite, plagioclase, biotite, rutile, spinel, corundum, and minor zircon and apatite. The quartz-bearing mafic granulites contain a peak pressure assemblage of garnet + clinopyroxene + ternary feldspar (now mesoperthite) + quartz + rutile. The sapphirine-bearing metabasites occur as mafic layers in garnet peridotite. Petrographical data suggest a peak assemblage of garnet + clinopyroxene + kyanite + rutile. Early kyanite is inferred from a symplectite of sapphirine + corundum + plagioclase ± spinel, interpreted to have formed during decompression. Garnet peridotite contains an assemblage of garnet + olivine + orthopyroxene + clinopyroxene. Thermobarometry indicates that all rock types experienced peak P–T conditions of 18.5–27.3 kbar and 870–1050 °C. A medium–high pressure granulite facies overprint (780–820 °C, 9.5–12 kbar) is defined by the formation of secondary clinopyroxene ± orthopyroxene + plagioclase at the expense of garnet and early clinopyroxene in the mafic granulites, as well as by growth of spinel and plagioclase at the expense of garnet and kyanite in the felsic granulite. SHRIMP II zircon U-Pb geochronology yields ages of 493 ± 7 Ma (mean of 11) from the felsic granulite, 497 ± 11 Ma (mean of 11) from sapphirine-bearing metabasite and 501 ± 16 Ma (mean of 10) from garnet peridotite. Rounded zircon morphology, cathodoluminescence (CL) sector zoning, and inclusions of peak metamorphic minerals indicate these ages reflect HP/HT metamorphism. Similar ages determined for eclogites from the western segment of the SAT suggest that the same continental subduction/collision event may be responsible for HP metamorphism in both areas.

Notes: A new quantitative approach to constraining mineral equilibria in sapphirine-bearing ultrahigh-temperature (UHT) granulites through the use of pseudosections and compatibility diagrams is presented, using a recently published thermodynamic model for sapphirine. The approach is illustrated with an example from an UHT locality in the Anápolis–Itauçu Complex, central Brazil, where modelling of mineral equilibria indicates peak metamorphic conditions of about 9 kbar and 1000 °C. The early formed, coarse-grained assemblage is garnet–orthopyroxene–sillimanite–quartz, which was subsequently modified following peak conditions. The retrograde pressure–temperature (P–T) path of this locality involves decompression across the FeO–MgO–Al2O3–SiO2 (FMAS) univariant reaction orthopyroxene + sillimanite = garnet + sapphirine + quartz, resulting in the growth of sapphirine–quartz, followed by cooling and recrossing of this reaction. The resulting microstructures are modelled using compatibility diagrams, and pseudosections calculated for specific grain boundaries considered as chemical domains. The sequence of microstructures preserved in the rocks constrains a two-stage isothermal decompression–isobaric cooling path. The stability of cordierite along the retrograde path is examined using a domainal approach and pseudosections for orthopyroxene–quartz and garnet–quartz grain boundaries. This analysis indicates that the presence or absence of cordierite may be explained by local variation in aH2O. This study has important implications for thermobarometric studies of UHT granulites, mainly through showing that traditional FMAS petrogenetic grids based on experiments alone may overestimate P–T conditions. Such grids are effectively constant aH2O sections in FMAS-H2O (FMASH), for which the corresponding aH2O is commonly higher than that experienced by UHT granulites. A corollary of this dependence of mineral equilibria on aH2O is that local variations in aH2O may explain the formation of cordierite without significant changes in P–T conditions, particularly without marked decompression.

Notes: The gabbroic/dioritic Pembroke Hornblende Granulite (PHG) of Milford Sound displays a geometrically simple mesoscopic network of sub-planar garnet reaction zones (GRZ) in which the meta-igneous hornblende granulite has been depleted of Na, Si, and H2O, and c. 25 vol.% almandine-rich garnet has formed. Some studies postulate the initial presence of melt along the centres of all GRZ, explaining the frequent absence of feldspathic veins by selective melt loss. A more parsimonious model is necessitated by structural evidence and, together with chemical data, suggests a relationship between mid-range metasomatic transport and anatexis. The Pembroke outcrops show a process of incipient melting of gabbro/diorite in an environment of relatively low aH2O in lithologies that have limited free quartz. A non-equilibrium steady state is proposed, in which a sodic dehydration fluid moves some distance via the GRZ network towards areas of partial melting. Only in these areas are Na and Si reconstituted as albite, with more garnet as byproduct, having avoided the need for melt percolation. The combined structural and chemical evidence directs a focus on mass transport in low-aH2O gabbroic environments. In subsequent events of shearing and complete transposition, both sets of garnet – the atypical GRZ residue and partial melt melanosomes – were inherited by the Milford Gneiss ‘facies’ of the PHG.

Notes: Hydration reactions are direct evidence of fluid–rock interaction during regional metamorphism. In this study, hydration reactions to produce retrograde actinolite in mafic schists are investigated to evaluate the controlling factors on the reaction progress. Mafic schists in the Sanbagawa belt contain amphibole coexisting with epidote, chlorite, plagioclase and quartz. Amphibole typically shows two types of compositional zoning from core to rim: barroisite → hornblende → actinolite in the high-grade zone, and winchite → actinolite in the low-grade zone. Both types indicate that amphibole grew during the exhumation stage of the metamorphic belt. Microstructures of amphibole zoning and mass-balance relations suggest that: (1) the actinolite-forming reactions proceeded at the expense of the preexisting amphibole; and (2) the breakdown reaction of hornblende consumed more H2O fluid than that of winchite, when one mole of preexisting amphibole was reacted. Reaction progress is indicated by the volume fraction of actinolite to total amphibole, Yact, with the following details: (1) reaction proceeded homogeneously in each mafic layer; (2) the extent of the hornblende breakdown reaction is commonly low (Yact 〈 0.5), but it increases drastically in the high-grade part of the garnet zone (Yact 〉 0.7); and (3) the extent of the winchite breakdown reaction is commonly high (Yact 〉 0.7). Many microcracks are observed within hornblende, and the extent of hornblende breakdown reaction is correlated with the size reduction of the hornblende core. Brittle fracturing of hornblende may have enhanced retrograde reaction progress by increasing of influx of H2O and the surface area of hornblende. In contrast to high-grade rocks, the winchite breakdown reaction is well advanced in the low-grade rocks, where reaction progress is not associated with brittle fracturing of winchite. The high extent of the reaction in the low-grade rocks may be due to small size of winchite before the reaction.

Notes: The combination of metamorphic petrology tools and in situ laser 40Ar/39Ar dating on phengite (linking time of growth, compositions and P–T conditions) enables us to identify a detailed P–T–d–t path for the still debated tectonometamorphic evolution of the Nevado-Filabride complex and infer new geodynamic-scale constraints. Our data show an isothermal decompression (at 550 °C) from 20 kbar for the Bédar-Macael unit and 14 kbar for the Calar Alto unit down to c. 3–4 kbar for both units at 2.8 mm year−1. At 22–18 Ma, this first part of the exhumation is followed by a final exhumation at 0.6 mm year−1 along a high-temperature low-pressure (HTLP) gradient of c. 60 °C km−1. The age of the peak of pressure is not precisely known but it is shown that it is around 30 Ma and possibly older, which is at variance with recent models suggesting a younger age for high-pressure (HP) metamorphism. Most of the exhumation is related to late-orogenic extension from c. 30 to 22–18 Ma. Thus the formation of the main ductile extensional shear zone, the Filabres Shear Zone (FSZ), occurred at 22–18 Ma and is clearly associated with a top-to-the-west shear sense once the FSZ is well localized. The transition from ductile to brittle then occurred at c. 14 Ma. The final exhumation, accommodated by brittle deformation, occurred from c. 14 to 9 Ma and was accompanied, from 12 to 8 Ma, by the formation of nearby extensional basins. The duration of the extensional process is c. 20 Myr which argues in favour of a progressive slab retreat from c. 30 to 9 Ma. The change in the shape of the P–T path at 22–18 Ma together with strain localization along the main top-to-the-west shear zone suggests that this date corresponds to a change in the direction of slab retreat from southwards to westwards.

Notes: Thermal zoning of the Highland Complex, Sri Lanka has been delineated using the Fe2+–Mg distribution coefficient between garnet and biotite from garnet–biotite gneiss samples collected with wide geographical distribution. In order to minimize the potential for retrograde Fe–Mg exchange and maximize the potential for retaining peak equilibrium KD (garnet–biotite) and temperature, garnet and biotite included within feldspar and quartz without other mineral inclusions have been selected. The calculated results indicate four distinct temperature contours with KD values varying from 1.84 to 6.38 and temperature varying from 996 to 591 °C. From the present results, it is possible to divide the Highland Complex into two major metamorphic zones: a high-temperature area in the central region and a low-temperature area in the south-western and north-eastern region. In conjunction with the metamorphic pressure variations estimated from the granulites of the Highland Complex in previous studies, it is shown that the high- and low-temperature areas are complemented by a high-pressure region towards the eastern side and a low-pressure region towards the western side of this complex. This thermal dome is interpreted to be an artifact of the different crustal levels exhumed following Pan-African metamorphism.

Notes: Monazite grains from Greater Himalayan Sequence gneisses, Langtang valley, Nepal, were chemically mapped and then dated in situ via Th–Pb ion-microprobe analysis. Correlation of ages and chemistry reveals at least five different generations of monazite, ranging from c. 9 to 〉300 Ma. Petrological models of monazite chemistry provide a link between these generations and the thermal evolution of these rocks, yielding an age for the melting of Greater Himalayan rocks within the Main Central Thrust sheet (c. 16 Ma), and for the timing of thrust sheet emplacement that are younger than commonly viewed. Chemical characterization of monazite is vital prior to chronological microanalysis, and many ages previously reported for monazite from the Greater Himalayan Sequence are interpretationally ambiguous.

Notes: This study presents calcite–graphite carbon isotope fractionations for 32 samples from marble in the northern Elzevir terrane of the Central Metasedimentary Belt, Grenville Province, southern Ontario, Canada. These results are compared with temperatures calculated by calcite–dolomite thermometry (15 samples), garnet–biotite thermometry (four samples) and garnet–hornblende thermometry (three samples). Δcal-gr values vary regularly across the area from 〉6.5‰ in the south to 4.0‰ in the north, which corresponds to temperatures of 525 °C in the south to 650 °C in the north. Previous empirical calibration of the calcite–graphite thermometer agrees very well with calcite–dolomite, garnet–biotite and garnet–hornblende thermometry, whereas, theoretical calibrations compare less well with the independent thermometry. Isograds in marble based on the reactions rutile + calcite + quartz =titanite and tremolite + calcite + quartz = diopside, span temperatures of 525–600 °C and are consistent with calculated temperature–X(CO2) relations. Results of this study compare favourably with large-scale regional isotherms, however, local variation is greater than that revealed by large-scale sampling strategies. It remains unclear whether the temperature–Δcal-gr relationship observed in natural materials below 650 °C represents equilibrium fractionations or not, but the regularity and consistency apparent in this study demonstrate its utility for thermometry in amphibolite facies marble.

Notes: Oxygen-isotope compositions of kyanite, andalusite, prismatic sillimanite and fibrolite from the Proterozoic terrane in the Truchas Mountains, New Mexico differ from one another, suggesting that these minerals did not grow in equilibrium at the Al2SiO5 (AS) polymorph-invariant point as previously suggested. Instead, oxygen-isotope temperature estimates indicate that growth of kyanite, andalusite and prismatic sillimanite occurred at c. 575, 615 and 640 °C respectively. Temperature estimates reported in this paper are interpreted as those of growth for the different AS polymorphs, which are not necessarily the same as peak metamorphic temperatures for this terrane. Two distinct temperature estimates of c. 580 °C and c. 700 °C are calculated for most fibrolite samples, with two samples yielding clear evidence of quartz-fibrolite oxygen-isotope disequilibrium. These data indicate that locally, and potentially regionally, oxygen-isotope disequilibrium between quartz and fibrolite may have resulted from rapid fibrolite nucleation. Pressures of mineral growth that were extrapolated from oxygen-isotope thermometry results and calculated using petrological constraints suggest that kyanite and one generation of fibrolite grew during M1 at 5 kbar, and that andalusite, prismatic sillimanite and a second generation of fibrolite grew during M2 at 3.5 kbar. M1 and M2 therefore represent two distinct metamorphic events that occurred at different crustal levels. The ability of the AS polymorphs to retain δ18O values of crystallization make these minerals ideal to model prograde-growth histories of mineral assemblages in metamorphic terranes and to understand more clearly the pressure–temperature histories of multiple metamorphic events.

Notes: The integration of information which can be gained from accessory [i.e. age (t)] and rock-forming minerals [i.e. temperature (T) and pressure (P)] requires a more profound understanding of the equilibration kinetics during metamorphic processes. This paper presents an approach comparing conventional P–T estimate from equilibrated assemblages of rock-forming minerals with temperature data derived from yttrium-garnet-monazite (YGM) and yttrium-garnet-xenotime (YGX) geothermometry. Such a comparison provides an initial indication on differences between equilibration of major and trace elements. Regarding this purpose, two migmatites, two polycyclic and one monocyclic gneiss from the Central Alps (Switzerland, northern Italy) were investigated. While the polycyclic samples exhibit trace-element equilibration between monazite and garnet grains assigned to the same metamorphic event, there are relics of monazite and garnet obviously surviving independent of their textural position. These observations suggest that surface processes dominate transport processes during equilibration of those samples. The monocyclic gneiss, on the contrary, displays rare isolated monazite with equilibration of all elements, despite comparably large transport distances. With a nearly linear crystal-size distribution of the garnet grain population, growth kinetics, related to the major elements, were likely surface-controlled in this sample. In contrast to these completely equilibrated examples, the migmatites indicate disequilibrium between garnet and monazite with a change in REE patterns on garnet transects. The cause for this disequilibrium may be related to a potential disequilibrium initiated by a changing bulk chemistry during melt segregation. While migmatite environments are expected to support high transport rates (i.e. high temperatures and melt presence), the evolution of equilibration in migmatites is additionaly related to change in chemistry. As a key finding, surface-controlled equilibration kinetics seem to dominate transport-controlled processes in the investigated samples. This may be decisive information towards the understanding of age data derived from monazite.

Notes: The presence of ternary feldspar in high-grade meta-igneous rocks, and the recognition of the thermometric significance of this mineral, has led recent researchers to postulate peak metamorphic temperatures in excess of 1000 °C. However, it needs to be established that such ternary feldspar is not in fact a survivor of the original high-temperature crystallization of the igneous protolith. After exsolution, the host and lamellae in the ternary feldspar grains may be stable throughout subsequent history as long as recrystallization does not occur. Such a history may involve rehydration and metamorphism, including H2O-saturated conditions, with the compositions and proportions of the host and lamellae being modified to reflect the P–T conditions experienced. In the case of the high-grade meta-igneous rocks from the Moldanubian of the Bohemian Massif, some samples that contain ternary feldspar preserve a substantial measure of their igneous heritage. Orthopyroxene-bearing granulites not only include types that are barely affected by the metamorphism, but also others that have undergone hydration of the igneous protolith prior to the development of a metamorphic overprint. A key to establishing the igneous origin of the ternary feldspar grains is their preservation in garnet that is either itself igneous, or of a relatively low-temperature metamorphic origin. Applying the logic to the other ternary feldspar-bearing meta-igneous rocks deprives the Moldanubian of its ultrahigh temperature (UHT) metamorphic status.

Notes: The effect of fluids on recrystallization behaviour is well known; however, the detailed microscale distribution of fluid in grain boundaries and the influence of fluid on grain boundary migration are still unresolved. In this study, in-situ deformation experiments in transmitted light microscopy were undertaken, as this allows continuous and direct observation of the whole range of processes involved in fluid-assisted grain boundary migration. A new see-through deformation apparatus was developed to enable the control of fluid pressure. Bischofite containing small amounts of aqueous fluid was deformed at temperatures between 50 and 90 °C, over a range of fluid pressure from 0.5 to 1 MPa, and strain rates of 5 × 10−6 to 1 × 10−4 s−1. The rates of grain boundary migration were measured at different temperatures and strain rates. Detailed observations during and after the deformation illustrate the evolution of migrating fluid-filled grain boundaries and show that the incorporation of fluids from inclusions as well as their pinch-off is dependent on the grain boundary velocity, the thickness of the grain boundary and the size and shape of the inclusions. Direct evidence is presented for the contraction of the grain boundary fluids into isolated inclusions after equilibrium conditions are attained.

Notes: A mid-ocean ridge basalt (MORB)-type eclogite from the Moldanubian domain in the Bohemian Massif retains evidence of its prograde path in the form of inclusions of hornblende, plagioclase, clinopyroxene, titanite, ilmenite and rutile preserved in zoned garnet. Prograde zoning involves a flat grossular core followed by a grossular spike and decrease at the rim, whereas Fe/(Fe + Mg) is also flat in the core and then decreases at the rim. In a pseudosection for H2O-saturated conditions, garnet with such a zoning grows along an isothermal burial path at c. 750 °C from 10 kbar in the assemblage plagioclase-hornblende-diopsidic clinopyroxene-quartz, then in hornblende-diopsidic clinopyroxene-quartz, and ends its growth at 17–18 kbar. From this point, there is no pseudosection-based information on further increase in pressure or temperature. Then, with garnet-clinopyroxene thermometry, the focus is on the dependence on, and the uncertainties stemming from the unknown Fe3+ content in clinopyroxene. Assuming no Fe3+ in the clinopyroxene gives a serious and unwarranted upward bias to calculated temperatures. A Fe3+-contributed uncertainty of ±40 °C combined with a calibration and other uncertainties gives a peak temperature of 760 ± 90 °C at 18 kbar, consistent with no further heating following burial to eclogite facies conditions. Further pseudosection modelling suggests that decompression to c. 12 kbar occurred essentially isothermally from the metamorphic peak under H2O-undersaturated conditions (c. 1.3 mol.% H2O) that allowed the preservation of the majority of garnet with symplectitic as well as relict clinopyroxene. The modelling also shows that a MORB-type eclogite decompressed to c. 8 kbar ends as an amphibolite if it is H2O saturated, but if it is H2O-undersaturated it contains assemblages with orthopyroxene. Increasing H2O undersaturation causes an earlier transition to SiO2 undersaturation on decompression, leading to the appearance of spinel-bearing assemblages. Granulite facies-looking overprints of eclogites may develop at amphibolite facies conditions.

Notes: A series of striking migmatitic structures occur in rectilinear networks through western Fiordland, New Zealand, involving, for the most part, narrow anorthositic dykes that cut hornblende-bearing orthogneiss. Adjacent to the dykes, host rocks show patchy, spatially restricted recrystallization and dehydration on a decimetre-scale to garnet granulite. Although there is general agreement that the migration of silicate melt has formed at least parts of the structures, there is disagreement on the role of silicate melt in dehydrating the host rock. A variety of causal processes have been inferred, including metasomatism due to the ingress of a carbonic, mantle-derived fluid; hornblende-breakdown leading to water release and limited partial melting of host rocks; and dehydration induced by volatile scavenging by a migrating silicate melt. Variability in dyke assemblage, together with the correlation between dehydration structures and host rock silica content, are inconsistent with macroscopic metasomatism, and best match open system behaviour involving volatile scavenging by a migrating trondhjemitic liquid.

Notes: Upper amphibolite facies felsic gneiss from Broken Hill records the metatexite to schlieren diatexite to massive diatexite transition in a single rock type over a scale of tens to hundreds of metres. The metatexites are characterized by centimetre-scale segregation of melt into leucosomes to form stromatic migmatite. The schlieren diatexites are characterized by the disaggregation of the rocks and the development of schlieren migmatite. The massive diatexites represent a higher degree of disaggregation, lack schlieren and contain plagioclase and K-feldspar phenocrysts. The transition from metatexite to schlieren diatexite and massive diatexite was heterogeneous with both disaggregation of the rock on a grain scale and disaggregation of the rock into centimetre- to metre-scale rafts. As melt contents increased, the proportion of material disaggregated on a grain scale increased. The high proportion of melt needed to form diatexites at upper amphibolite facies conditions was the result of an influx of hydrous fluid at temperatures just above the solidus of the diatexites. Nearby metapelitic rocks, with a slightly higher solidus temperature, undergoing subsolidus muscovite breakdown are the likely source of the fluid. Continued heating during and after the influx of fluid led to melt contents of up to c. 60 mol.% in the massive diatexite. The metatexite zone probably involved little added fluid. Continued deformation during cooling and melt crystallization resulted in the extensive development of schlieren and late-stage melt segregations and melt-rich shear bands in the schlieren diatexite zone. The rocks of the massive diatexite zone lack these late-stage segregations, consistent with the cessation of D2 deformation prior to them developing a crystal framework.

Notes: Cambrian orogenesis (550–490 Ma) in the Lambert Province of the southern Prince Charles Mountains resulted in three successive stages of deformation. The earliest of these deformations resulted in the development of a layer-parallel foliation (S1) that was folded into macro-scale recumbent folds (F2). Subsequent deformation buckled the rocks into long-wavelength (c. 20 km), SW- to NW-trending antiformal closures (F3) mostly separated from each other by west to SW trending, steeply dipping, high-strain zones. Metapelitic rocks from the region are divisible into two compositional types: a high-Al, -Fe and -K type and a high-Mg, -Ca and -Na type. In rocks of both composition, relic staurolite preceded the formation of upper amphibolite facies garnet + biotite + sillimanite ± muscovite mineral assemblages that record peak pressures and temperatures of c. 650–700 °C and 6–7 kbar. Subsequent decompression of c. 3 kbar is implied from texturally late plagioclase and a reduction in the modal abundance of garnet in the high-Al, -Fe and -K metapelites, and from texturally late cordierite in the more magnesium rocks. This clockwise P–T–t path, with prograde heating followed by rapid decompression, is: (i) equivalent to that recorded in the same-aged rocks at Prydz Bay located 600 km to the north, and (ii) similar to the modelled response of the crust to thickening following continent–continent collision. These results indicate that large areas of East Antarctica were thickened and rapidly exhumed, probably in response to collisional orogenesis during the Early Cambrian. This supports the inference that Early Cambrian orogenesis in the Prydz Bay–Prince Charles Mountains region of East Antarctica marks one of the fundamental lithospheric boundaries within Gondwana.

Notes: Interlayered graphitic and non-graphitic schists from the Tauern Window, Eastern Alps, record contrasting mechanical behaviour during extensional exhumation. Graphitic schists contain mesoscale extension fractures, pervasive microcracks in garnet, and abundant secondary fluid inclusion planes; all three types of structures are oriented perpendicular to the stretching lineation. Crack spacings in garnet from graphitic samples are tightly clustered around a mean of 180 μm. Non-graphitic schists have fewer and more randomly oriented microcracks and fluid inclusion planes and maintained strain compatibility via crystal plasticity. The presence or absence of graphite appears to have exerted a fundamental control on rheology during unroofing. Calculations for a model graphitic rock at 500 °C and fO2 = 10−24 MPa show that the equilibrium metamorphic fluid evolves from XCO2 = 0.07 to 0.38 during decompression from 700 to 400 MPa, in agreement with microcrack fluid inclusion data that show a change from XCO2 〈 0.1 to 0.45 in graphitic samples over the same pressure interval. This compositional shift results in 〉60% expansion of the pore fluid during decompression. H2O-rich fluid in non-graphitic rocks expands 〈15% over the same pressure interval. The greater pore fluid expansion in low-permeability graphitic horizons likely promoted tensile failure during unroofing. These results suggest that microcracking should be an inevitable consequence of decompression in many graphitic schists, whereas rocks that lack graphite are less likely to undergo microcracking. Microseismicity is predicted to be more common in graphitic than non-graphitic rocks during unroofing of mountain belts.

Notes: The 3D shape, size and orientation data for white mica grains sampled along two transects of increasing metamorphic grade in the Otago Schist, New Zealand, reveal that metamorphic foliation, as defined by mica shape-preferred orientation (SPO), developed rapidly at sub-greenschist facies conditions early in the deformation history. The onset of penetrative strain metamorphism is marked by the rapid elimination of poorly oriented large clastic mica in favour of numerous new smaller grains of contrasting composition, higher aspect ratios and a strong preferred orientation. The metamorphic mica is blade shaped with long axes defining the linear aspect of the foliation and intermediate axes a partial girdle about the lineation. Once initiated, foliation progressively intensified by an increase in the aspect ratio, size and alignment of grains, although highest grade samples within the chlorite zone record a decrease in aspect ratio and reduction in SPO strength despite continued increase in grain size. These trends are interpreted in terms of progressive competitive anisotropic growth of blade-shaped grains so that the fastest growth directions and blade lengths tend to parallel the extension direction during deformation. The competitive nature of mica growth is indicated by the progressive increase in size and resultant decrease in number of metamorphic mica with increasing grade, from c. 1000 relatively small mica grains per square millimetre of thin section at lower grades, to c. 100 relatively large grains per square millimetre in higher grade samples. Reversal of SPO intensity and grain aspect ratio trends in higher grade samples may reflect a reduction in the strain rate or reduction in the deviatoric component of the stress field.

Notes: The assumption of oxygen isotope and major element equilibrium during prograde metamorphism was tested using staurolite-grade pelitic schists that have undergone sequential porphyroblast growth and multiple episodes of recrystallization of matrix minerals and foliation development. Textural relationships are used to infer a metamorphic history that involves garnet growth followed by staurolite growth, with each porphyroblast growth event followed by at least one period of recrystallization of matrix minerals. Conventional geothermobarometry using Qtz–Grt–Pl–Ms–Bt ± St equilibria yields peak P–T conditions of c. 625 °C at 9–11 kbar, consistent with KMnFMASH petrogenetic grid predictions for stability of the assemblage Grt + St + Bt. Qtz–Grt oxygen isotope fractionations yield apparent temperatures of c. 590 °C and Qtz–St fractionations yield an apparent temperature of c. 595 °C. Diffusional modelling indicates that quartz isotopic compositions were reset by c. 30 °C via retrograde isotopic diffusional exchange with micas. The isotopic temperatures appear to be in excellent agreement with one another, and suggest oxygen isotope equilibrium was attained between garnet and staurolite at c. 625 °C. However, the agreement of Qtz–Grt and Qtz–Str isotopic temperatures is not consistent with petrographic observations (garnet grew before staurolite) and petrogenetic grid constraints that predict that garnet grows over a temperature interval of c. 525–550 °C. Given that: (i) oxygen diffusion rates in staurolite and garnet are slow enough to render an individual porphyroblast effectively closed to exchange after it forms; and (ii) matrix minerals are able to exchange isotopes via recrystallization during each period of deformation; garnet and staurolite could not have simultaneously achieved oxygen isotope equilibrium with each other or with minerals in the recrystallized matrix. Thus, the Qtz–Grt fractionations, which yield apparent temperatures that are in apparent agreement with peak metamorphic temperature and apparent temperatures for Qtz–St fractionations, cannot be fractionations resulting from equilibrium isotopic exchange. Instead, they are apparent fractionations between porphyroblasts formed at different temperature and times in the prograde P–T–D path, and quartz that recrystallized and exchanged with micas and plagioclase during several phases of deformation.

Notes: Shape, size and orientation measurements of quartz grains sampled along two transects that cross zones of increasing metamorphic grade in the Otago Schist, New Zealand, reveal the role of quartz in the progressive development of metamorphic foliation. Sedimentary compaction and diagenesis contributed little to the formation of a shape-preferred orientation (SPO) within the analysed samples. Metamorphic foliation was initiated at sub-greenschist facies conditions as part of a composite S1-bedding structure parallel to the axial planes of tight to isoclinal F1 folds. An important component of this foliation is a pronounced quartz SPO that formed dominantly by the effect of dissolution–precipitation creep on detrital grains in association with F1 strain. With increasing grade, the following trends are evident from the SPO data: (i) a progressive increase in the aspect ratio of grains in sections parallel to lineation, and the development of blade-shaped grains; (ii) the early development of a strong shape preferred orientation so that blade lengths define the linear aspect of the foliation (lineation) and the intermediate axes of the blades define a partial girdle about the lineation; (iii) a slight thinning and reduction in volume of grains in the one transect; and (iv) an actual increase in thickness and volume in the survivor grains of the second transect. The highest-grade samples, within the chlorite zone of the greenschist facies, record segregation into quartz- and mica-rich layers. This segregation resulted largely from F2 crenulation and marks a key change in the distribution, deformation and SPO of the quartz grains. The contribution of quartz SPO to defining the foliation lessens as the previously discrete and aligned detrital quartz grains are replaced by aggregates and layers of dynamically recrystallized quartz grains of reduced aspect ratio and reduced alignment. Pressure solution now affects the margins of quartz-rich layers rather than individual grains. In higher-grade samples, therefore, the rock structure is characterized increasingly by segregation layering parallel to a foliation defined predominantly by mica SPO.

Notes: The Tres Arboles ductile fault zone in the Eastern Sierras Pampeanas, central Argentina, experienced multiple ductile deformation and faulting events that involved a variety of textural and reaction hardening and softening processes. Much of the fault zone is characterized by a (D2) ultramylonite, composed of fine-grained biotite + plagioclase, that lacks a well-defined preferred orientation. The D2 fabric consists of a strong network of intergrown and interlocking grains that show little textural evidence for dislocation or dissolution creep. These ultramylonites contain gneissic rock fragments and porphyroclasts of plagioclase, sillimanite and garnet inherited from the gneissic and migmatitic protolith (D1) of the hangingwall. The assemblage of garnet + sillimanite + biotite suggests that D1-related fabrics developed under upper amphibolite facies conditions, and the persistence of biotite + garnet + sillimanite + plagioclase suggests that the ultramylonite of D2 developed under middle amphibolite facies conditions. Greenschist facies, mylonitic shear bands (D3) locally overprint D2 ultramylonites. Fine-grained folia of muscovite + chlorite ± biotite truncate earlier biotite + plagioclase textures, and coarser-grained muscovite partially replaces relic sillimanite grains. Anorthite content of shear band (D3) plagioclase is c. An30, distinct from D1 and D2 plagioclase (c. An35). The anorthite content of D3 plagioclase is consistent with a pervasive grain boundary fluid that facilitated partial replacement of plagioclase by muscovite. Biotite is partially replaced by muscovite and/or chlorite, particularly in areas of inferred high strain. Quartz precipitated in porphyroclast pressure shadows and ribbons that help define the mylonitic fabric. All D3 reactions require the introduction of H+ and/or H2O, indicating an open system, and typically result in a volume decrease. Syntectonic D3 muscovite + quartz + chlorite preferentially grew in an orientation favourable for strain localization, which produced a strong textural softening. Strain localization occurred only where reactions progressed with the infiltration of aqueous fluids, on a scale of hundreds of micrometre. Local fracturing and microseismicity may have induced reactivation of the fault zone and the initial introduction of fluids. However, the predominant greenschist facies deformation (D3) along discrete shear bands was primarily a consequence of the localization of replacement reactions in a partially open system.

Notes: Metapelites, migmatites and granites from the c. 2 Ga Mahalapye Complex have been studied for determining the P–T–fluid influence on mineral assemblages and local equilibrium compositions in the rocks from the extreme southwestern part of the Central Zone of the Limpopo high-grade terrane in Botswana. It was found that fluid infiltration played a leading role in the formation of the rocks. This conclusion is based on both well-developed textures inferred to record metasomatic reactions, such as Bt ⇒ And + Qtz + (K2O) and Bt ± Qtz ⇒ Sil + Kfs + Ms ± Pl, and zonation of Ms | Bt + Qtz | And + Qtz and Grt | Crd | Pl | Kfs + Qtz reflecting a perfect mobility (Korzhinskii terminology) of some chemical components. The conclusion is also supported by the results of a fluid inclusion study. CO2 and H2O (〈inlineGraphic alt="inline image" href="urn:x-wiley:02634929:JMG579:JMG_579_mu1" location="equation/JMG_579_mu1.gif"/〉 = 0.6) are the major components of the fluid. The fluid has been trapped synchronously along the retrograde P–T path. The P–T path was derived using mineral thermobarometry and a combination of mineral thermometry and fluid inclusion density data. The Mahalapye Complex experienced low-pressure granulite facies metamorphism with a retrograde evolution from 770 °C and 5.5 kbar to 560 °C and 2 kbar, presumably at c. 2 Ga.

Notes: Quartzo-feldspathic veins emplaced within a migmatite terrane near Wilson Lake in the Grenville Province of central Labrador record a metamorphic event not evident in the host rocks. The discordant veins are undeformed and have undisturbed primary igneous/hydrothermal textures. Most of the veins contain euhedral kyanite, as well as aggregates of kyanite, K-feldspar, phlogopite and minor dumortierite which are likely pseudomorphs after primary phengite. The reconstructed phengite compositions range from 3.1 to 3.2 Si per 11 oxygen formula unit. The pseudomorph assemblage is interpreted as the product of phengite + quartz melting under H2O-undersaturated conditions, which brackets P–T conditions of formation to about 9–16 kbar and 775–875 °C. A parallel vein that is likely of the same generation contains the borosilicate phases, dumortierite, prismatine and grandidierite, but no kyanite. The borosilicate assemblages constrain the P–T conditions of vein crystallization to ≥10 kbar and c. 750–850 °C. Vein emplacement is constrained to T ≤ 875 °C at the same pressures, which is well within the kyanite zone.Because the host rocks and veins must have experienced the same P–T history following vein emplacement, the presence of unreacted sillimanite in the host migmatites implies insufficient time for host rock equilibration. Slow reaction rates because of anhydrous conditions are not a likely explanation given the abundance of biotite and hornblende in the host rocks. The ductility implied by the breakdown of a hydrous phase (phengite) and the production of an H2O-undersaturated melt in the veins contrasts with the apparently brittle behaviour of the host rocks. The absence of deformation since the time of vein emplacement, even at temperatures above 750 °C, suggests that the deep crust in this part of Labrador had a very short residence time under conditions of the kyanite zone. Rapid decompression from those conditions is consistent with quartz + phengite melting and accounts for the relatively brittle behaviour of the terrane as it was uplifted.

Notes: Regionally metamorphosed metapelites from Rogaland, SW Norway, contain zircon formed during the decompression reaction garnet + sillimanite + quartz → cordierite. The zircon, which occurs as inclusions in cordierite coronas around garnet, is texturally, chemically and isotopically distinct from older zircon in other textural settings in the matrix. A SHRIMP U–Pb age of 955 ± 8 Ma based on analyses in thin section on the decompression zircon from the cordierite coronas, therefore dates a point on the retrograde path, estimated from garnet–cordierite equilibria to be 5.6 kbar, 710 °C. This population was under-represented in conventional SHRIMP analyses of individual zircon in a mono-mineralic grain mount and, in the absence of a textural context, its significance unknown. The dominant age identified from SHRIMP analyses of the grain mount, in combination with analyses from matrix zircon in thin section, was 1035 ± 9 Ma. Based on the lack of consistent textural relationships with any specific minerals in thin section, as well as rare earth element chemistry, the 1035-Ma population is interpreted to represent zircon growth during incipient migmatization of the rocks at 6–8 kbar and c. 700 °C. This is consistent with previous estimates for the age of regional M1 metamorphism during the Sveconorwegian Orogeny. The most important outcome of this study is the successful analysis of zircon grains in a specific, well-constrained reaction texture. Not only does this allow a precise point on the regional P–T path to be dated, but it also emphasizes the possibility of zircon formation during the retrograde component of a typical metamorphic cycle.

Notes: The subduction and exhumation of accretionary prism metasedimentary rocks are accompanied by large-strain ductile deformations which may be recorded in microstructures. Porphyroblast microstructures have been a key to unravel the kinematics in such deformed belts. Shape-preferred orientation (SPO) of epidote and amphibole inclusions that define S-shaped trails in prograde cores of plagioclase porphyroblasts were analysed from the high-P/T Sambagawa metamorphic rocks. Inclusions are found to be elongate parallel to the [010] and [001] directions, respectively, and their long-axis orientations define an internal foliation Si (best-fit great circle) and lineation Li (maximum on the Si). S-shaped inclusion trails in the orthogonal sections do not exhibit the same geometries, but rather are grouped into two types, where the foliation intersection axes (FIAs) are nearly perpendicular and parallel to Li, respectively. These two types of S-shaped inclusion trails are seen in the sections inclined at low and high angles to the Li, respectively. However, the latter type commonly consists of composite trails, where the Si is first rotated about an FIA perpendicular to the Li (i.e. unique axis), and then about an FIA parallel to the Li. The S-shaped inclusion trails are interpreted to have formed by the successive overgrowth of matrix minerals and rotation of the plagioclase porphyroblast cores about a unique axis in non-coaxial deformation. The rotation of Si about an FIA nearly parallel to the Li is perhaps an apparent rotation, caused by the deflection of foliation around the growing prismatic plagioclase grain prior to inclusion into the porphyroblast. This study has for the first time documented the 3-D geometry of S-shaped inclusion trails in porphyroblasts from accretionary prism metasedimentary rocks and identified their origin, which helps to understand the flow kinematics in the deeper part of a subduction channel.

Notes: In the Orlica-Śnieżnik complex at the NE margin of the Bohemian Massif, high-pressure granulites occur as isolated lenses within partially migmatized orthogneisses. Sm–Nd (different grain-size fractions of garnet, clinopyroxene and/or whole rock) and U–Pb [isotope dilution-thermal ionization mass spectrometry (ID-TIMS) single grain and sensitive high-resolution ion microprobe (SHRIMP)] ages for granulites, collected in the surroundings of Červený Důl (Czech Republic) and at Stary Gierałtów (Poland), constrain the temporal evolution of these rocks during the Variscan orogeny. Most of the new ages cluster at c. 350–340 Ma and are consistent with results previously reported for similar occurrences throughout the Bohemian Massif. This interval is generally interpreted to constrain the time of high-pressure metamorphism. A more complex evolution is recorded for a mafic granulite from Stary Gierałtów and concerns the unknown duration of metamorphism (single, short-lived metamorphic cycle or different episodes that are significantly separated in time?). The central grain parts of zircon from this sample yielded a large spread in apparent 206Pb/238U SHRIMP ages (c. 462–322 Ma) with a distinct cluster at c. 365 Ma. This spread is interpreted to be indicative for variable Pb-loss that affected magmatic protolith zircon during high-grade metamorphism. The initiating mechanism and the time of Pb-loss has yet to be resolved. A connection to high-pressure metamorphism at c. 350–340 Ma is a reasonable explanation, but this relationship is far from straightforward. An alternative interpretation suggests that resetting is related to a high-temperature event (not necessarily in the granulite facies and/or at high pressures) around 370–360 Ma, that has previously gone unnoticed. This study indicates that caution is warranted in interpreting U–Pb zircon data of HT rocks, because isotopic rejuvenation may lead to erroneous conclusions.

Notes: The development of shear zones at mid-crustal levels in the Proterozoic Willyama Supergroup was synchronous with widespread fluid flow resulting in albitization and calcsilicate alteration. Monazite dating of shear zone fabrics reveal that they formed at 1582 ± 22 Ma, at the end of the Olarian D3 deformational event and immediately prior to the emplacement of regional S-type granites. Two stages of fluid flow are identified in the area: first an albitizing event which involved the addition of Na and loss of Si, K and Fe; and a second phase of calcsilicate alteration with additions of Ca, Fe, Mg and Si and removal of Na. Fluid fluxes calculated for albitization and calcsilicate alteration were 5.56 × 109 to 1.02 × 1010 mol m−2 and 2.57 × 108–5.20 × 109 mol m−2 respectively. These fluxes are consistent with estimates for fluid flow through mid-crustal shear zones in other terranes. The fluids associated with shearing and alteration are calculated to have δ18O and δD values ranging between +8 and +11‰, and −33 and −42‰, respectively, and ɛNd values between −2.24 and −8.11. Our results indicate that fluids were derived from metamorphic dehydration of the Willyama Supergroup metasediments. Fluid generation occurred during prograde metamorphism of deeper crustal rocks at or near peak pressure conditions. Shear zones acted as conduits for major crustal fluid flow to shallow levels where peak metamorphic conditions had been attained earlier leading to the apparent ‘retrograde’ fluid-flow event. Thus, the peak metamorphism conditions at upper and lower crustal levels were achieved at differing times, prior to regional granite formation, during the same orogenic cycle leading to the formation of retrograde mineral assemblages during shearing.

Notes: The thermodynamic properties of silicate minerals can be described as a linear combination of the fractional properties of their constituent polyhedra. In contrast, given the thermodynamic properties of these polyhedra, the thermodynamic properties of minerals can be estimated, where only the crystallography of the mineral needs to be known. Such estimates are especially powerful for hypothetical mineral end-members or for minerals where experimental determination of their thermodynamic properties is difficult. In this contribution the fractional enthalpy, entropy and molar volume for 35 polyhedra have been determined using weighted multiple linear regression analysis on a data set of published mineral thermodynamic properties. The large number of polyhedra determined, allows calculation of a much larger variety of phases than was previously possible and the larger set of minerals used provides more confident fractional properties. The OH-bearing minerals have been described by partial and total hydroxide coordinated components, which gives better results than previous models and precludes the need of a S–V term to improve estimates of entropy. However, the fractional thermodynamic properties only give adequate results for silicate minerals and double oxides, and should therefore not be used to estimate the properties of other minerals. The thermodynamic properties of ‘new’ minerals are calculated from a linear stoichiometric combination of their constituent polyhedra, resulting in estimates generally with associated uncertainty of 〈5%. The quality of such data appears to be of sufficient accuracy for thermodynamic modelling as shown for meta-bauxites from the Alps and the Aegean, where the effect of Zn on the P–T stability of staurolite can be both qualitatively and quantitatively reproduced.

Notes: Reaction progress exhibited by multivariant assemblages in micaceous limestones can provide an excellent record of metamorphic fluid flow. However, it is necessary to understand the sensitivity of these assemblages to bulk-composition parameters. Here, analysis of bulk composition on different scales and pseudosection construction are used to draw conclusions on relationships between bulk composition, fluid flow and reaction progress. Issues addressed include the effects of bulk composition on the mineralogical evolution of micaceous carbonates, the sensitivity of bulk composition to bulk-composition sampling methods, the magnitude of cross-layer fluid-composition gradients, the potential for metasomatism to drive reaction progress, and the relative timing of reaction in adjacent layers. Pseudosections successfully represent observed mineral assemblages, constrain the position of reactions in T–X(CO2) space, and allow assessment of the sensitivity of reaction position, inferred reaction progress and calculated fluid fluxes to uncertainties in bulk composition. The scale of bulk-composition sampling affects bulk compositions, calculated modes, predicted mineral assemblages and calculated fluid compositions. Larger samples record an average of different lithological subdomains, while point-count-derived bulk compositions are subject to uncertainties related to the small number of sample points. The optimum bulk composition for pseudosection purposes probably lies between measured bulk compositions. Results suggest that reaction progress in some extensively reacted layers was driven by infiltration of H2O-rich fluid which flowed or diffused parallel to layering, perpendicular to layering in response to fluid-composition gradients, and out of veins. Small variations in fluid composition across layering (ΔX(CO2) 〈 0.02) were maintained by internal buffering by the mineral assemblages. Internal buffering must also have driven samples up a sequence of narrow low-variance fields in T–X(CO2) space, and so reaction in adjacent layers must have close to simultaneous. Metasomatic effects on reaction progress are likely to have been small, so long as the porosity was low.

Notes: Metapelitic residual enclaves in the Neogene Volcanic Province of SE Spain are residues left after melt extraction. Glass (quenched melt) of granitic composition occurs as inclusions in most minerals and as intergranular pockets. The most common enclave types show one stage of garnet growth that is interpreted to have occurred at the same time as glass production. Some of these show a well-developed foliation outlined by fibrolite, biotite, graphite and glass, which wraps around elongate garnet crystals that have aspect ratios up to 10:1. Based on microstructures and chemistry, the garnet within these rocks shows clear core and mantle structure. The core has an average composition of Alm76–Prp08–Sps14–Grs03 and contains primary inclusions of biotite and melt, trapped during garnet growth. A thin (c. 100 μm), irregular mantle overgrows the garnet core, enclosing oriented fibrolite inclusions in strain caps, and biotite in strain shadows. In places, the overgrowths form skeletal elongated arms, which extend parallel to the foliation. The garnet mantle contains less Mn and higher XMg, but both core and mantle display flat Mn profiles, the contact being a sharp break. Ternary feldspar and Grt–Bt thermometry yield temperatures in the range 800–900 °C, with no systematic differences among the different microstructural domains of elliptical garnet. Based on the observed intracrystalline microstructures, the high amount of melt extraction in the rock by flattening component strain and the chemical zoning of garnet, the formation of elliptical garnet is modelled by a multistage sequence. This involves pressure solution and reprecipitation of the core, followed by post-kinematic, partly mimetic growth of the garnet mantle.

Notes: Thermal models for Barrovian metamorphism driven by doubling the thickness of the radiogenic crust typically meet difficulty in accounting for the observed peak metamorphic temperature conditions. This difficulty suggests that there is an additional component in the thermal budget of many collisional orogens. Theoretical and geological considerations suggest that viscous heating is a cumulative process that may explain the heat deficit in collision orogens. The results of 2D numerical modelling of continental collision involving subduction of the lithospheric mantle demonstrate that geologically plausible stresses and strain rates may result in orogen-scale viscous heat production of 0.1 to 〉1 μW m−3, which is comparable to or even exceeds bulk radiogenic heat production within the crust. Thermally induced buoyancy is responsible for crustal upwelling in large domes with metamorphic temperatures up to 200 °C higher than regional background temperatures. Heat is mostly generated within the uppermost mantle, because of large stresses in the highly viscous rocks deforming there. This thermal energy may be transferred to the overlying crust either in the form of enhanced heat flow, or through magmatism that brings heat into the crust advectively. The amplitude of orogenic heating varies with time, with both the amplitude and time-span depending strongly on the coupling between heat production, viscosity and collision strain rate. It is argued that geologically relevant figures are applicable to metamorphic domes such as the Lepontine Dome in the Central Alps. We conclude that deformation-generated viscous dissipation is an important heat source during collisional orogeny and that high metamorphic temperatures as in Barrovian type metamorphism are inherent to deforming crustal regions.

Notes: Metanorites from two eclogitized metagabbros of the Hercynian French Massif Central preserve coronitic textures of hornblende, garnet, quartz and/or kyanite produced at the expense of the primary magmatic assemblage orthopyroxene and plagioclase. Using a petrogenetic grid in the CFMASH system, two possible P–T evolutions for the origin of the coronas are evaluated. The sequence of reactions involving the formation of Hbl (–Ky) ± Grt and Qtz coronitic assemblages is consistent with an isobaric cooling at high pressure (c. 1–2 GPa) under hydrated conditions. However, this P–T path, inferred by using only petrographical observations, is inconsistent with the geochronological constraints: emplacement of the gabbro at 490 Ma and high-pressure metamorphism at 410 Ma. In order to reconcile petrographical observations with geochronological constraints, we propose a discontinuous two-stage evolution involving a change in water activity with time. (1) Emplacement and cooling of the norite at low pressure under anhydrous conditions, at 490 Ma. (2) During the Hercynian orogeny, the norite experienced an increase in pressure and temperature under fluid-present conditions. Adding water to the system implies a dramatic change in the petrogenetic grid topology, restricting the orthopyroxene–plagioclase assemblage only to high temperatures. Therefore, the breakdown of the unstable magmatic assemblage, through apparent retrograde reactions, occurred along the prograde P–T path which never crossed the equilibrium boundaries of these reactions.

Notes: Detailed petrographic analysis was performed on samples from five localities within the southern Adirondacks. Textures and zoning patterns in garnet from all samples provide evidence for dehydration melting of biotite. Zoning of grossular in garnet – providing a record of prograde growth – shows both increasing and decreasing trends in garnet from the same sample. However, Ca concentrations at the garnet rims of most samples are identical (grossular = 3.4%). These observations have been interpreted as evidence for the differential timing of garnet nucleation and growth. All Fe/(Fe + Mg) and some spessartine distributions are consistent between samples, displaying diffusive profiles established largely upon cooling. Only one sample, in which retrogression was minimal, contains garnet with flat Fe/(Fe + Mg) profiles. A general pelitic pseudosection constructed in the system MnNCKFMASH reveals a maximum for Ca in garnet where the plagioclase-out isopleth intersects the solidus (muscovite = 0). The pseudosection predicts bell-shaped core-to-rim profiles of grossular during anatexis, similar to those observed in the rocks. Garnet–biotite thermometry and GASP barometry indicate peak temperatures of at least 790 °C at about 7–9 kbar, similar to conditions determined for the central Adirondacks. Cooling rates determined from finite difference modelling of spessartine and Fe/(Fe + Mg) diffusional profiles indicate a multi-stage cooling history in which some period of rapid cooling (〉200 °C Myr−1) is required.

Notes: Oxide–sulphide–Fe–Mg–silicate and titanite–ilmenite textures as well as their mineral compositions have been studied in felsic and intermediate orthogneisses across an amphibolite (north) to granulite facies (south) traverse of lower Archean crust, Tamil Nadu, south India. Titanite is limited to the amphibolite facies terrane where it rims ilmenite or occurs as independent grains. Pyrite is widespread throughout the traverse increasing in abundance with increasing metamorphic grade. Pyrrhotite is confined to the high-grade granulites. Ilmenite is widespread throughout the traverse increasing in abundance with increasing metamorphic grade and occurring primarily as hemo-ilmenite in the high-grade granulite facies rocks. Magnetite is widespread throughout the traverse and is commonly associated with ilmenite. It decreases in abundance with increasing metamorphic grade. In the granulite facies zone, reaction rims of magnetite + quartz occur along Fe–Mg silicate grain boundaries. Magnetite also commonly rims or is associated with pyrite. Both types of reaction rims represent an oxidation effect resulting from the partial subsolidus reduction of the hematite component in ilmenite to magnetite. This is confirmed by the presence of composite three oxide grains consisting of hematite, magnetite and ilmenite. Magnetite and magnetite–pyrite micro-veins along silicate grain boundaries formed over a wide range of post-peak metamorphic temperatures and pressures ranging from high-grade SO2 to low-grade H2S-dominated conditions. Oxygen fugacities estimated from the orthopyroxene–magnetite–quartz, orthopyroxene–hematite–quartz, and magnetite–hematite buffers average 2.5 log units above QFM. It is proposed that the trends in mineral assemblages, textures and composition are the result of an external, infiltrating concentrated brine containing an oxidizing component such as CaSO4 during high-grade metamorphism later acted upon by prograde and retrograde mineral reactions that do not involve an externally derived fluid phase.

Notes: Silica-undersaturated, sapphirine-bearing granulites occur in a large number of localities worldwide. Such rocks have historically been under-utilized for estimating P–T evolution histories because of limited experimental work, and a consequent poor understanding of the topology and P–T location of silica-undersaturated mineral equilibria. Here, a calculated P–T projection for sapphirine-bearing, silica-undersaturated metapelitic rock compositions is constructed using THERMOCALC for the FeO-MgO-Al2O3-SiO2 (FMAS) and KFMASH (+K2O + H2O) chemical systems, allowing quantitative analysis of silica-undersaturated mineral assemblages. This study builds on that for KFMASH sapphirine + quartz equilibria [Kelsey et al. (2004) Journal of Metamorphic Geology, vol. 22, pp. 559–578]. FMAS equilibria are significantly displaced in P–T space from silicate melt-bearing KFMASH equilibria. The large number of univariant silica-undersaturated KFMASH equilibria result in a P–T projection that is topologically more complex than could be established on the basis of experiments and/or natural assemblages. Coexisting sapphirine and silicate melt (with or without corundum) occur down to c. 900 °C in KFMASH, some 100 °C lower than in silica-saturated compositions, and from pressures of c.≤1 to ≥12 kbar. Mineral compositions and composition ranges for the calculated phases are consistent with natural examples. Bulk silica has a significant effect on the stability of sapphirine-bearing assemblages at a given P–T, resulting in a wide variety of possible granulite facies assemblages in silica-undersaturated metapelites. Calculated pseudosections are able to reproduce many naturally occurring silica-undersaturated assemblages, either within a single assemblage field or as the product of a P–T trajectory crossing several fields. With an understanding of the importance of bulk composition on sapphirine stability and textural development, silica-undersaturated assemblages may be utilized in a quantitative manner in the detailed metamorphic investigation of high-grade terranes.

Notes: Several petrographic studies have linked accessory monazite growth in pelitic schist to metamorphic reactions involving major rock-forming minerals, but little attention has been paid to the control that bulk composition might have on these reactions. In this study we use chemographic projections and pseudosections to argue that discrepant monazite ages from the Mount Barren Group of the Albany–Fraser Orogen, Western Australia, reflect differing bulk compositions. A new Sensitive High-mass Resolution Ion Microprobe (SHRIMP) U–Pb monazite age of 1027 ± 8 Ma for pelitic schist from the Mount Barren Group contrasts markedly with previously published SHRIMP U–Pb monazite and xenotime ages of c. 1200 Ma for the same area. All dated samples experienced identical metamorphic conditions, but preserve different mineral assemblages due to variable bulk composition. Monazite grains dated at c. 1200 Ma are from relatively magnesian rocks dominated by biotite, kyanite and/or staurolite, whilst c. 1027 Ma grains are from a ferroan rock dominated by garnet and staurolite. The latter monazite population is likely to have grown when staurolite was produced at the expense of garnet and chlorite, but this reaction was not intersected by more magnesian compositions, which are instead dominated by monazite that grew during an earlier, greenschist facies metamorphic event. These results imply that monazite ages from pelitic schist can vary depending on the bulk composition of the host rock. Samples containing both garnet and staurolite are the most likely to yield monazite ages that approximate the timing of peak metamorphism in amphibolite facies terranes. Samples too magnesian to ever grow garnet, or too iron-rich to undergo garnet breakdown, are likely to yield older monazite, and the age difference can be significant in terranes with a polymetamorphic history.

Notes: Fan-shaped polycrystalline staurolite porphyroblasts, 3–4 cm in length and 0.5 cm in width, occur together with centimetre-sized euhedral prismatic staurolite porphyroblasts in pelitic schists of the Littleton Formation on the western overturned limb of the Bolton syncline in eastern Connecticut. The fans consist of intergrown planar splays of [001] elongated prisms, which are crudely radial from a single apex. The apical angles of the radial groupings range up to 70°. The orientations of the individual staurolite prisms are related by a rigid rotation about an axis perpendicular to the fan plane. The zone axes [001] always lie in the plane of the fan. Although the angle between the [100] zone axes of the individual prisms is uniform in each fan, it ranges from +30° to −30° in different fans. Internally, the fans display: (i) remnants of a passively captured Si foliation defined by disc-shaped quartz blebs (type 1 inclusions) and layers of very fine carbonaceous material and tabular ilmenite platelets; (ii) bent staurolite blades and undulose extinction along low-angle (010) subgrain boundaries near the apex of the fans; (iii) wedge-shaped dilatational zones containing equigranular inclusion-free quartz, mica and staurolite, and (iv) growth-related quartz inclusion trails roughly perpendicular to a crystal face (type 2 inclusions). The Si inclusion trails are typically perpendicular to the fan surface, radiate parallel to the blades, and show little to no curvature except at the very edge of the fans where they abruptly curve through nearly 90° into parallelism with an external Se foliation. Careful examination of the three-dimensional geometry of fans based on U-stage measurements, serial sections and two-circle optical goniometric measurements permits a detailed reconstruction of their sequential development. The origin of a fan involves limited intracrystalline deformation and brittle crack dilation, spalling, rotation, and growth of small marginal fragments and of new staurolite along wedge-shaped zones along the Si inclusion surfaces. Fans preferentially develop in porphyroblasts in which Si is subparallel to the 010 cleavage. These internal features reflect the rotation and deformation of a brittle porphyroblast relative to syn-growth shear stresses.

Notes: Detailed microtextural observations and bulk chemical analysis were undertaken on a garnet-pyroxenite nodule within retrograde eclogites from the NE Sulu ultrahigh-pressure metamorphic (UHPM) terrane. The results suggest that the protolith was a cumulate from a gabbroic body. The nodule consists primarily of coarse clinopyroxene grains with a very high content of the Ca-Tschermakite molecule. Microscopic observations and back-scattered electron images (BSE) demonstrate a complicated intergrowth of clinopyroxene, garnet and ilmenite, which represents the peak metamorphic assemblage. The primary clinopyroxene grains are armoured with a thin garnet corona up to 0.5 mm wide that forms an interconnected network. Within the clinopyroxene grains, four sets of garnet lamellae are distributed along crystallographic planes; locally, a vermicular intergrowth of garnet and diopside is developed. Besides the garnet, parallel arrays of ilmenite blebs are common within the clinopyroxene. Hydrous minerals such as amphibole, zoisite and titanite formed at later stages, and replaced diopside, garnet and ilmenite respectively. The P–T conditions determined for the formation of the garnet lamellae indicate that the garnet pyroxenite experienced UHP metamorphism at the same peak P–T condition as its host eclogite. The very high Ca-Tschermakite content (31–34 mol.%) of the primary clinopyroxene indicates crystallization at about 9–17 kbar and 1250–1450 °C, and together with the microtextural observations, suggests that the protolith of the garnet pyroxenite was a cumulate from a former gabbroic body, in which case, the host eclogite might represent the gabbroic body.

Notes: Thermodynamic calculations in petrology are generally performed at pressures and temperatures beyond the standard state conditions. Accurate prediction of mineral equilibria therefore requires knowledge of the heat capacity, thermal expansion and compressibility for the minerals involved. Unfortunately, such data are not always available. In this contribution we present a data set to estimate the heat capacity, thermal expansion and compressibility of mineral end-members from their constituent polyhedra, based on the premise that the thermodynamic properties of minerals can be described by a linear combination of the fractional properties of their constituents. As such, only the crystallography of the phase of interest needs to be known. This approach is especially powerful for hypothetical mineral end-members and for minerals, for which the experimental determination of their thermodynamic properties is difficult. The data set consists of the properties for 35 polyhedra in the system K–Na–Ca–Li–Be–Mg–Mn–Fe–Co–Ni–Zn–Al–Ti–Si–H, determined by multiple linear regression analysis on a data set of 111 published end-member thermodynamic properties. The large number of polyhedra determined allows calculation of a much larger variety of phases than was previously possible, and the choice of constituents together with the large number of thermodynamic input data results in estimates with associated uncertainty of generally 〈5%. The quality of the data appears to be sufficiently accurate for thermodynamic modelling as demonstrated by modelling the stability of margarite in the CASH system and the position of the talc–staurolite–chloritoid–pyrope absent invariant point in the KMASH system. In both cases, our results overlap within error with published equivalents.

Notes: The prograde pressure–temperature (P–T) path for the complexly polydeformed Proterozoic Broken Hill Block (Australia) has been reconstructed through detailed structural analysis in conjunction with calculation of compositionally specific P–T pseudosections of pelitic rock units within a high-temperature shear zone that formed early in the tectonic evolution of the terrane. Whilst the overall P–T path for the Broken Hill Block has been interpreted to be anticlockwise, the prograde portion of this path has been unresolved. Our results have constrained part of this prograde path, showing an early heating event (M1) at P–T conditions of at least c. 600 °C and c. 2.8–4.2 kbar, associated with an elevated geothermal gradient (c. 41–61 °C km−1). This event is interpreted to be the result of rifting at c. 1.69–1.67 Ga, or at c. 1.64–1.61 in the Broken Hill Block. Early rifting was followed by an episode of lithospheric thermal relaxation and burial, during which time sag-phase sediments of the upper Broken Hill stratigraphy (Paragon Group) were deposited. Following sedimentation, a second tectonothermal event (M2/D2) occurred. This event is associated with peak low-pressure granulite facies metamorphism (c. 1.6 Ga) and attained conditions of at least 740 °C at c. 5 kbar. A regionally pervasive, high-temperature fabric (S2) developed during the M2/D2 event, and deformation was accommodated along lithology-parallel, high-temperature shear zones. The larger-scale deformation regime (extensional or shortening) of this event remains unresolved. The M2/D2 event was terminated by intense crustal shortening during the Olarian Orogeny, during which time the first mappable folds within the Broken Hill Block developed.

Notes: A calibration is presented for an activity–composition model for amphiboles in the system Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O–O (NCFMASHO), formulated in terms of an independent set of six end-members: tremolite, tschermakite, pargasite, glaucophane, ferroactinolite and ferritschermakite. The model uses mixing-on-sites for the ideal-mixing activities, and for the activity coefficients, a macroscopic multicomponent van Laar model. This formulation involves 15 pairwise interaction energies and six asymmetry parameters. Calibration of the model is based on the geometrical constraints imposed by the size and shape of amphibole solvi inherent in a data set of 71 coexisting amphibole pairs from rocks, formed over 400–600 °C and 2–18 kbar. The model parameters are calibrated by combining these geometric constraints with qualitative consideration of parameter relationships, given that the data are insufficient to allow all the model parameters to be determined from a regression of the data. Use of coexisting amphiboles means that amphibole activity–composition relationships are calibrated independently of the thermodynamic properties of the end-members. For practical applications, in geothermobarometry and the calculation of phase diagrams, the amphibole activity–composition relationships are placed in the context of the stability of other minerals by evaluating the properties of the end-members in the independent set that are in internally consistent data sets. This has been performed using an extended natural data set for hornblende–garnet–plagioclase–quartz, giving the small adjustments necessary to the enthalpies of formation of tschermakite, pargasite and glaucophane for working with the Holland and Powell data set.

Notes: Integrated petrological and structural investigations of eclogites from the eclogite zone of the Voltri Massif (Ligurian Alps) have been used to reconstruct a complete Alpine P–T deformation path from burial by subduction to subsequent exhumation. The early metamorphic evolution of the eclogites has been unravelled by correlating garnet zonation trends with the chemical variations in inclusions found in the different garnet domains. Garnet in massive eclogites displays typical growth zoning, whereas garnet in foliated eclogites shows rim-ward resorption, likely related to re-equilibration during retrogressive evolution. Garnet inclusions are distinctly different from core to rim, consisting primarily of Ca-, Na/Ca-amphibole, epidote, paragonite and talc in garnet cores and of clinopyroxene ± talc in the outer garnet domains. Quantitative thermobarometry on the inclusion assemblages in the garnet cores defines an initial greenschist-to-amphibolite facies metamorphic stage (M1 stage) at c. 450–500 °C and 5–8 kbar. Coexistence of omphacite + talc + katophorite inclusion assemblage in the outer garnet domains indicate c. 550 °C and 20 kbar, conditions which were considered as minimum P–T estimates for the M2 eclogitic stage. The early phase of retrograde reactions is polyphase and equilibrated under epidote–blueschist facies (M3 stage), characterized by the development of composite reaction textures (garnet necklaces and fluid-assisted Na-amphibole-bearing symplectites) produced at the expense of the primary M2 garnet-clinopyroxene assemblage. The blueschist retrogression is contemporaneous with the development of a penetrative deformation (D3) that resulted in a non-coaxial fabric, with dominant top-to-the-N sense of shear during rock exhumation. All of that is overprinted by a texturally late amphibolite/greenschist facies assemblages (M4 & M5 stages), which are not associated with a penetrative structural fabric. The combined P–T deformation data are consistent with an overall counter-clockwise path, from the greenschist/amphibolite, through the eclogite, the blueschist to the greenschist facies. These new results provide insights into the dynamic evolution of the Tertiary oceanic subduction processes leading to the building up of the Alpine orogen and the mechanisms involved in the exhumation of its high-pressure roots.

Notes: A convolution integral is developed to evaluate transient, drilling speed–dependent discharge rates into a tunnel gradually excavated in a homogeneous infinite aquifer. Comparison with the classical instantaneous-drilling evaluation commonly used in practice reveals drastically reduced early-time discharge rates, higher maximum rates, and similar long-term rates. Dimensionless-type curves are provided to help assess total discharge sensitivity to drilling time and predict safer maximum flow rates.

Notes: Conceptual geological models based on geophysical data can elucidate aquifer architecture and heterogeneity at meter and smaller scales, which can lead to better predictions of preferential flow pathways. The macrodispersion experiment (MADE) site, with 〉2000 measurements of hydraulic conductivity obtained and three tracer tests conducted, serves as an ideal natural laboratory for examining relationships between subsurface flow characteristics and geophysical attributes in fluvial aquifers. The spatial variation of hydraulic conductivity measurements indicates a large degree of site heterogeneity. To evaluate the usefulness of geophysical methods for better delineating fluvial aquifer heterogeneities and distribution of preferential flow paths, a surface grid of two-dimensional ground penetrating radar (GPR) and direct current (DC) resistivity data were collected. A geological model was developed from these data that delineate four stratigraphic units with distinct electrical and radar properties including (from top to bottom) (1) a meandering fluvial system (MFS); (2) a braided fluvial system (BFS); (3) fine-grained sands; and (4) a clay-rich interval. A paleochannel, inferred by other authors to affect flow, was mapped in the MFS with both DC resistivity and GPR data. The channel is 2 to 4 m deep and, based on resistivity values, is predominantly filled with clay and silt. Comparing previously collected hydraulic conductivity measurements and tracer-plume migration patterns to the geological model indicates that flow primarily occurs in the BFS and that the channel mapped in the MFS has no influence on plume migration patterns.

Notes: Discrete-fracture and dual-porosity models are infrequently used to simulate solute transport through fractured unconsolidated deposits, despite their more common application in fractured rock where distinct flow regimes are hypothesized. In this study, we apply four fracture transport models—the mobile-immobile model (MIM), parallel-plate discrete-fracture model (PDFM), and stochastic and deterministic discrete-fracture models (DFMs)—to demonstrate their utility for simulating solute transport through fractured till. Model results were compared to breakthrough curves (BTCs) for the conservative tracers potassium bromide (KBr), pentafluorobenzoic acid (PFBA), and 1,4-piperazinediethanesulfonic acid (PIPES) in a large-diameter column of fractured till. Input parameters were determined from independent field and laboratory methods. Predictions of Br BTCs were not significantly different among models; however, the stochastic and deterministic DFMs were more accurate than the MIM or PDFM when predicting PFBA and PIPES BTCs. DFMs may be more applicable than the MIM for tracers with small effective diffusion coefficients (De) or for short timescales due to differences in how these models simulate diffusion or incorporate heterogeneities by their fracture networks. At large scales of investigation, the more computationally efficient MIM and PDFM may be more practical to implement than the three-dimensional DFMs, or a combination of model approaches could be employed. Regardless of the modeling approach used, fractures should be incorporated routinely into solute transport models in glaciated terrain.

Notes: System dynamics is a computer-aided approach to evaluating the interrelationships of different components and activities within complex systems. Recently, system dynamics models have been developed in areas such as policy design, biological and medical modeling, energy and the environmental analysis, and in various other areas in the natural and social sciences. The Idaho National Engineering and Environmental Laboratory, a multipurpose national laboratory managed by the Department of Energy, has developed a system dynamics model in order to evaluate its utility for modeling large complex hydrological systems. We modeled the Bear River basin, a transboundary basin that includes portions of Idaho, Utah, and Wyoming. We found that system dynamics modeling is very useful for integrating surface water and ground water data and for simulating the interactions between these sources within a given basin. In addition, we also found that system dynamics modeling is useful for integrating complex hydrologic data with other information (e.g., policy, regulatory, and management criteria) to produce a decision support system. Such decision support systems can allow managers and stakeholders to better visualize the key hydrologic elements and management constraints in the basin, which enables them to better understand the system via the simulation of multiple “what-if” scenarios. Although system dynamics models can be developed to conduct traditional hydraulic/hydrologic surface water or ground water modeling, we believe that their strength lies in their ability to quickly evaluate trends and cause-effect relationships in large-scale hydrological systems, for integrating disparate data, for incorporating output from traditional hydraulic/hydrologic models, and for integration of interdisciplinary data, information, and criteria to support better management decisions.

Notes: This paper presents a general model of optimal water management for a transboundary aquifer under three different management approaches: cooperative, noncooperative, and myopic. Comparing the results from the approaches, we find the cooperative solution, where a single management plan is executed for all parts of the aquifer, results in the highest level of net social welfare, followed by the noncooperative and then the myopic. The trade-offs we find for the higher levels of welfare are lower use levels in the earlier periods. We present a short discussion of factors that can increase the inefficiencies of the models and suggest directions for future research.

Notes: The goal for any property rights system is to achieve equity, efficiency, and certainty. Trying to achieve these goals for ground water is difficult because a ground water right is not exclusive. To make matters more complicated, ground water is often under the jurisdiction of more than one political unit. The result is transboundary conflicts. Two critical elements must be included in any system of ground water rights. The system must define how the ground water can be used and define the relationships that each user and each use has with the other users and uses in the system. Unfortunately, these relationships are seldom completely defined and are made more complex by the transboundary scales at which they operate. As ground water moves horizontally across boundaries, different users or different jurisdictions have sequential control, creating conflicts between the first users and subsequent ones. Other problems occur because of vertical relationships, with more than one person or entity having control over ground water at the same time. This simultaneous exercise of authority can create conflicts between an individual who possesses a right to use ground water and a state or federal agency that regulates the same water. Transboundary conflicts occur at different scales and include conflicts between neighboring property owners as well as conflicts between countries. Scale, the property rights structure, and the nature of the relationship between users influence the way transboundary ground water conflicts are resolved.

Notes: In Europe, a long history of cooperation over transboundary rivers—most notably the Rhine and Danube rivers—exists. To help foster cooperation and communication vis-à-vis transboundary ground water, the United Nations Economic Commission for Europe (UNECE), as part of its ground water program, conducted a survey on transboundary aquifers in Europe. The survey produced 25 responses from 37 countries and identified 89 transboundary aquifers. Respondents reported on the degree of ground water use within their own boundaries, transboundary aspects (agreements, joint commissions, etc.) of ground water, and transboundary aquifers themselves. The inventory proved useful, but a number of problems were identified: different map scales and symbols, difficulty in identifying transboundary aquifers, inconsistent labeling of aquifers, and data discrepancies. The UNECE ground water program also drafted guidelines for monitoring and assessment of transboundary ground water. These guidelines are not legally binding but have been adopted by 25 countries, deal mainly with monitoring and assessment, and are being implemented through a number of pilot projects. Other organizations—the United Nations Scientific, Educational and Cultural Organization, the Food and Agriculture Organization, the International Association of Hydrogeologists, and the European Union—are all supporting the investigation of transboundary aquifers in an effort to facilitate data sharing and coordinated management of these valuable resources.

Notes: More than one-half of the world's population is dependent on ground water for everyday uses such as drinking, cooking, and hygiene. In fact, it is the most extracted natural resource in the world. As a result of growing populations and expanding economies, many aquifers today are being depleted while others are being contaminated. Notwithstanding the world's considerable reliance on this resource, ground water resources have long received only secondary attention as compared to surface water, especially among legislatures and policymakers. Today, while there are hundreds of treaties governing transboundary rivers and lakes, there is only one international agreement that directly addresses a transboundary aquifer. Given that many of the aquifers on which humanity so heavily relies cross international borders, there is a considerable gap in the sound management, allocation, and protection of such resources. In order to prevent future disputes over transboundary aquifers and to maximize the beneficial use of this resource, international law must be clarified as it applies to transboundary ground water resources. Moreover, it must be defined with a firm basis in sound scientific understanding. In this paper we offer six conceptual models is which ground water resources can have transboudary consequences. The models are intended to help in assessing the applicability and scientific soundness of existing and proposed rules governing transboundary ground water resources. In addition, we consider the development of international law as it applies to ground water resources and make recommendations based on the models and principles of hydrogeology. The objective is the development of clear, logical, and science-based norms of state conducts as they relate to aquifers that traverse political boundaries.

Notes: A tracer test in a carbonate aquifer is analyzed using the method of moments and two analytical advection-dispersion models (ADMs) as well as a numerical model. The numerical model is a coupled continuum-pipe flow and transport model that accounts for two different flow components in karstified carbonate aquifers, i.e., rapid and often turbulent conduit flow and Darcian flow in the fissured porous rock. All techniques employed provide reasonable fits to the tracer breakthrough curve (TBC) measured at a spring. The resulting parameter estimates are compared to investigate how each conceptual model of flow and transport processes that forms the basis of the analyses affects the interpretation of the tracer test. Numerical modeling results suggest that the method of moments and the analytical ADMs tend to overestimate the conduit volume because part of the water discharged at the spring is wrongly attributed to the conduit system if flow in the fissured porous rock is ignored. In addition, numerical modeling suggests that mixing of the two flow components accounts for part of the dispersion apparent in the measured TBC, while the remaining part can be attributed to Taylor dispersion. These processes, however, cannot reasonably explain the tail of the TBC. Instead, retention in immobile-fluid regions as included in a nonequilibrium ADM provides a possible explanation.

Notes: Previous studies have used flowmeters in environments that are within the expectations of their published ranges. Electromagnetic flowmeters have a published range from 0.1 to 79.0 m/min, and impeller flowmeters have a published range from 1.2 to 61.0 m/min. Velocity-log data collected in five long-screened production wells in the Pleasant Valley area of southern California showed that (1) electromagnetic flowmeter results were comparable within ±2% to results obtained using an impeller flowmeter for comparable depths; (2) the measured velocities from the electromagnetic flowmeter were up to 36% greater than the published maximum range; and (3) both data sets, collected without the use of centralizers or flow diverters, produced comparable and interpretable results. Although either method is acceptable for measuring wellbore velocities and the distribution of flow, the electromagnetic flowmeter enables collection of data over a now greater range of flows. In addition, changes in fluid temperature and fluid resistivity, collected as part of the electromagnetic flowmeter log, are useful in the identification of flow and hydrogeologic interpretation.

Notes: While the discussion of model uncertainty has centered on spatial heterogeneity, it is possible that ground water models have not enjoyed much success as predictive tools often because the sources that were eventually imposed in the field differed from those represented in the simulations. This is because deterministic prediction of future conditions is often inaccurate due to the random nature of contaminant sources, in terms of their timing, location, and magnitude. This paper presents a stochastic framework for accommodating random contaminant sources in conventional, deterministic advection-dispersion transport models. The contaminant sources are first classified into two types: those occurring continuously with a deterministic component and random variations and those occurring randomly at instantaneous discrete-time intervals. For the first type, the governing partial differential equation (PDE) is replaced by a stochastic PDE. The random variations are modeled by Gaussian noise or Brownian motion, and the solution is obtained by using Ito's integration technique. For the second type, Markovian analysis is used for discrete-time contamination events. Both approaches use a deterministic transport model to generate response functions at any observation location and time. The response functions are then integrated to yield probabilistic description of contaminant transport, from which key statistical properties such as mean, standard deviation, and confidence interval can be drawn.

Notes: A model of porosity development has been developed to investigate general relationships between simple fracture aperture growth laws and fracture porosity in evolved fracture arrays in aquifers. The growth of fracture apertures in two-dimensional orthogonal arrays with initially spatially uncorrelated lognormal aperture distributions has been studied, where aperture growth rate is proportional to an exponent of the flow rate through each fracture. The evolved arrays show geometrical phase changes as a function of the aperture growth rate exponent, e, and the standard deviation of the initial aperture distribution, σz. Low values of e and σz lead to bimodal aperture distributions, where apertures parallel to flow are preferentially enlarged. At moderate values of e and σz, there is a transition to a regime of more complex geometries consisting of networks of channel-like structures of preferentially enlarged apertures. At larger values of e, array-spanning channel-like paths of preferentially enlarged apertures develop, where the tortuosity of the channel-like paths is a linear function of σz. Following an initial growth phase, during which dynamically stable aperture configurations develop, arrays undergo simple amplification. The geometry of the evolved aperture fields is diverse and they can be highly complex; consequently, parameterization and prediction of their evolution in terms of the initial aperture distributions and growth rate laws is not trivial.

Notes: A field-scale demonstration project was conducted to evaluate the capability of eastern cottonwood trees (Populus deltoides) to attenuate trichloroethene (TCE) contamination of ground water. By the middle of the sixth growing season, trees planted where depth to water was 〈3 m delivered enough dissolved organic carbon to the underlying aquifer to lower dissolved oxygen concentrations, to create iron-reducing conditions along the plume centerline and sulfate-reducing or methanogenic conditions in localized areas, and to initiate in situ reductive dechlorination of TCE. Apparent biodegradation rate constants for TCE along the centerline of the plume beneath the phytoremediation system increased from 0.0002/d to 0.02/d during the first six growing seasons. The corresponding increase in natural attenuation capacity of the aquifer along the plume centerline, from 0.0004/m to 0.024/m, is associated with a potential decrease in plume-stabilization distance from 9680 to 160 m. Demonstration results provide insight into the amount of vegetation and time that may be needed to achieve cleanup objectives at the field scale.

Notes: The hydraulic conductivity (K) variation has important ramifications for ground water flow and the transport of contaminants in ground water. The delineation of the nature of that variation can be critical to complete characterization of a site and the planning of effective and efficient remedial measures. Site-specific features (such as high-conductivity zones) need to be quantified. Our alluvial field site in the Kansas River valley exhibits spatial variability, very high conductivities, and nonlinear behavior for slug tests in the sand and gravel aquifer. High-resolution, multilevel slug tests have been performed in a number of wells that are fully screened. A general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the wellbore, and a Hvorslev model for the aquifer has been used to analyze the data, employing an automated processing system that runs within the Excel spreadsheet program. It is concluded that slug tests can provide the necessary data to identify the nature of both horizontal and vertical K variation in an aquifer and that improved delineation or higher resolution of K structure is possible with shorter test intervals. The gradation into zones of higher conductivity is sharper than seen previously, and the maximum conductivity observed is greater than previously measured. However, data from this project indicate that well development, the presence of fines, and the antecedent history of the well are important interrelated factors in regard to slug-test response and can prevent obtaining consistent results in some cases.

Notes: To improve understanding of property measurements in heterogeneous media, an energy-based weighting function concept is developed. In (assumed) homogeneous media, the instrument spatial weighting function (ISWF) depends only on the energy dissipation distribution set up by the measurement procedure and it reduces to simply inverse sample volume (uniform weighting) for 1-D parallel flow case (ideal permeameter). For 1-D transient flow in homogeneous media, such as with slug tests, the ISWF varies with position and time, with 95% of the total weighting contained within 115 well radii, even late in the test. In the heterogeneous case, the determination of the ISWF is connected to the problem of determining an equivalent hydraulic conductivity (K), where the criterion for equivalence is based on equal energy dissipation rate rather than equal volume discharge. The discharge-based equivalent K (KE) and the energy-based equivalent K in heterogeneous media (Keh) are not equal in general, with Keh typically above the nodal arithmetic mean K. The possibly more fundamental problem is that as one makes K measurements in heterogeneous media at different locations or on different cores of heterogeneous materials, the ISWF will be heterogeneity dependent, implying that the averaging process resulting in the equivalent K value also varies with position. If the testing procedure is transient, then the averaging process varies with time. This suggests a fundamental ambiguity in the interpretation of hydraulic conductivity measurements in heterogeneous media that may impact how we approach modeling and prediction in a practical sense (Molz 2003). Further research is suggested.

Notes: At a number of sites, a plume of methyl tert-butyl ether (MTBE) in ground water has dived below the screen of conventional monitoring wells and escaped detection. Techniques are needed to predict the vertical extent of MTBE in an aquifer. Two techniques that are emerging in the site characterization market are electrical conductivity logging and pneumatic slug testing performed in temporary push wells. These techniques were evaluated at a diving plume of MTBE in the aquifer that supplies water to the village of East Alton, Illinois. The plume stayed near the water table for the first 100 m from the potential sources and then dived below conventional monitoring over the next 100 m. At the location where the MTBE plume dived, the depth to water was 9.1 m below land surface. The first 10 m of material below the water table had an electrical conductivity near 100 mS/m, indicating silts and clays. An electrical conductivity near 25 mS/m, indicating sands or gravels, was encountered at a depth of 10.6 m below the water table, and the sands and gravel extended to a depth of at least 15.2 m below the water table. Pneumatic slug tests measured low hydraulic conductivity in the interval of silt and clay (0.34 and 0.012 m/d) and higher hydraulic conductivity in the interval with sands and gravels (12.5, 11.6, and 11.3 m/d). Ground water with the highest concentration of MTBE was produced just below the contact between the silt and clay and the sands and gravel.

Notes: Steam injected below the water table tends to move upward because of buoyancy. This limits the horizontal steam zone development, which determines the optimal spacing between injection wells. In this study, observations indicating steam override at a full-scale remediation of an unconfined aquifer are analyzed by numerical modeling using the code T2VOC. A simplified three-dimensional (3D) numerical model is set up, which qualitatively shows the same mechanisms as observed at the site. By means of the model, it is found that it will be possible to achieve a larger horizontal extent of the steam zone in a layered geology compared to the homogeneous case. In the homogeneous case, the steam injection rate increases dramatically when the injection pressure is increased, which is necessary to achieve a larger horizontal development. The development of the steam zone under unconfined conditions is found to be a complex function of the geologic layering, the water table at steady-state extraction, and the injection/extraction system. Because of this complexity, it will be difficult to predict steam behavior without site-specific 3D numerical modeling.

Notes: This work describes the integration of data acquisition (DAQ) hardware and software for the purpose of acquiring not only data but real-time transport model parameter estimates in the context of subsurface flow and transport problems. Integrated DAQ parameter estimation systems can be used to reduce data storage requirements, trigger event recognition and more detailed sampling actions, and otherwise enhance remote monitoring capabilities. The contaminant transport problem is posed here as the analogous heat transfer problem in a three-dimensional, intermediate-scale physical aquifer model. A constant source of warm water is fed into a sandy aquifer undergoing steady, unidirectional flow. The spatial distribution of temperature in the medium is monitored over time using 17 thermocouples embedded in the medium. These sensors log temperatures via conventional analog-to-digital conversion hardware driven by commercially available DAQ software (LabVIEW™). Parameter estimation routines programmed in MATLAB™-based M-files are embedded in the LabVIEW DAQ routine and access parameter estimation libraries, such as the descent method employed here, via the Internet. The integrated DAQ parameter estimation system is demonstrated for the estimation of (1) the thermal dispersion coefficients (analogous to mass dispersion coefficients), given a known heat source; and (2) the location of a heat source, given known thermal dispersion coefficients. In both cases, the parameter estimation procedure is executed repeatedly as the data are acquired. For the case of source location, the effect of the number of sensors on the parameter estimation procedure is also demonstrated. Reasonable parameter estimates are provided rapidly during both the transient and steady-state phases of the experiments, with accuracy increasing with time and with the number of observations employed.

Notes: Herbicide and herbicide degradate compounds were measured in tile drains, monitoring wells, domestic wells, and in the receiving stream at a farm in south-central Michigan, in which corn, soybeans, and livestock are the major products. Immunoassay screens of two herbicide groups (triazines and alachlor) in 2000 showed a decrease in the number of detections from June to September, after the well-known spring flush, of 48% to 26% for alachlor and 62% to 22% for the triazines. Gas chromatography/mass spectrometry and liquid chromatography/mass spectrometry analyses of 2001 samples documented a fall flush as well, mobilizing relatively high levels of both parent and degradate compounds to the tile drain system. The compounds included both metolachlor and atrazine, which were applied on the farm earlier that year, as well as the degradates of alachlor, which had not been applied on the farm for several years. In four tile drain samples, metolachlor was detected at concentrations ranging from 18 to 25 μg/L, whereas atrazine ranged from nondetect to 2.25 μg/L. Contrary to previous studies, metolachlor concentrations were higher than metolachlor ethane sulfonic acid, by a factor of 3 to 7. The high soil permeability, high water table, and shallow screened intervals of monitoring wells directly below the point of application may explain the lesser degree of breakdown at the point of sampling. Although most detections were below maximum contaminant level or health advisory levels, the unknown effects of degradates and the aggregate levels of compounds measured in many samples may pose higher health risks than for individual compounds.

Notes: Monitoring data collected over a 6-year period show that a plume of chlorinated ethene–contaminated ground water has contracted significantly following treatment of the contaminant source area using in situ oxidation. Prior to treatment (1998), concentrations of perchloroethene (PCE) exceeded 4500 μg/L in a contaminant source area associated with a municipal landfill in Kings Bay, Georgia. The plume emanating from this source area was characterized by vinyl chloride (VC) concentrations exceeding 800 μg/L. In situ oxidation using Fenton's reagent lowered PCE concentrations in the source area below 100 μg/L, and PCE concentrations have not rebounded above this level since treatment. In the 6 years following treatment, VC concentrations in the plume have decreased significantly. These concentration declines can be attributed to the movement of Fenton's reagent–treated water downgradient through the system, the cessation of a previously installed pump-and-treat system, and the significant natural attenuation capacity of this anoxic aquifer. While in situ oxidation briefly decreased the abundance and activity of microorganisms in the source area, this activity rebounded in 〈6 months. Nevertheless, the shift from sulfate-reducing to Fe(III)-reducing conditions induced by Fenton's treatment may have decreased the efficiency of reductive dechlorination in the injection zone. The results of this study indicate that source-area removal actions, particularly when applied to ground water systems that have significant natural attenuation capacity, can be effective in decreasing the areal extent and contaminant concentrations of chlorinated ethene plumes.

Notes: Two small-diameter (≤0.5-inch/13-mm O.D.) bladder pumps were evaluated under the U.S. EPA Environmental Technology Verification Program to assess their ability to provide representative water quality samples for several volatile organic compounds (VOCs) and major element cations. Both a pneumatic bladder pump (PnBP) and a newly designed mechanical bladder pump (MBP) were evaluated. Since the MBP does not require compressed air to operate, the field logistics are simplified and cost of conducting low-flow sampling can be reduced. Replicate samples were collected from an above-ground standpipe spiked with selected VOCs and cations. Low-flow sampling protocol was followed during collection of replicate samples from six direct push (DP)–installed monitoring wells at Tyndall Air Force Base, Florida. Simple statistical methods are used to evaluate the sample quality provided by the pumps as compared to reference samples. F-tests are used to evaluate pump sample precision, and t-tests are applied to evaluate comparability of the pump sample results to the colocated reference or port sample results. Statistical analysis reveals that the samples collected with the pumps provide precision comparable to that observed in the reference samples. Significantly greater variability was observed between the pump and reference samples collected from the DP wells than was observed during the standpipe trials. The results and statistics from the trials at the standpipe indicate that both the PnBP and MBP provide water quality samples that are comparable to the reference samples, for the analytes investigated.

Notes: Previous site-specific studies designed to assess the impacts of unsewered subdivisions on ground water quality have relied on upgradient monitoring wells or very limited background data to characterize conditions prior to development. In this study, an extensive monitoring program was designed to document ground water conditions prior to construction of a rural subdivision in south-central Wisconsin. Previous agricultural land use has impacted ground water quality; concentrations of chloride, nitrate-nitrogen, and atrazine ranged from below the level of detection to 296 mg/L, 36 mg/L, and 0.8 μg/L, respectively, and were highly variable from well to well and through time. Seasonal variations in recharge, surface topography, aquifer heterogeneities, surficial loading patterns, and well casing depth explain observed variations in ground water chemistry. This variability would not have been detected if background conditions were determined from only a few monitoring wells or inferred from wells located upgradient of the subdivision site. This project demonstrates the importance of characterizing both ground water quality and chemical variability prior to land-use change to detect any changes once homes are constructed.

Notes: From the mid-1940s through the 1980s, large volumes of waste water were discharged at the Hanford Site in southeastern Washington State, causing a large-scale rise (〉20 m) in the water table. When waste water discharges ceased in 1988, ground water mounds began to dissipate. This caused a large number of wells to go dry and has made it difficult to monitor contaminant plume migration. To identify monitoring wells that will need replacement, a methodology has been developed using a first-order uncertainty analysis with UCODE, a nonlinear parameter estimation code. Using a three-dimensional, finite-element ground water flow code, key parameters were identified by calibrating to historical hydraulic head data. Results from the calibration period were then used to check model predictions by comparing monitoring wells’ wet/dry status with field data. This status was analyzed using a methodology that incorporated the 0.3 cumulative probability derived from the confidence and prediction intervals. For comparison, a nonphysically based trend model was also used as a predictor of wells’ wet/dry status. Although the numerical model outperformed the trend model, for both models, the central value of the intervals was a better predictor of a wet well status. The prediction interval, however, was more successful at identifying dry wells. Predictions made through the year 2048 indicated that 46% of the wells in the monitoring well network are likely to go dry in areas near the river and where the ground water mound is dissipating.

Notes: A method is proposed to estimate the uncertainty of the location of pathlines in two-dimensional, steady-state confined or unconfined flow in aquifers due to the uncertainty of the spatially variable unconditional hydraulic conductivity or transmissivity field. The method is based on concepts of the semianalytical first-order theory given in Stauffer et al. (2002, 2004), which allows estimates of the lateral second moment (variance) of the location of a moving particle. However, this method is reformulated in order to account for nonuniform recharge and nonuniform aquifer thickness. One prominent application is the uncertainty estimation of the catchment of a pumping well by considering the boundary pathlines starting at a stagnation point. In this method, the advective transport of particles is considered, based on the velocity field. In the case of a well catchment, backtracking is applied by using the reversed velocity field. Spatial variability of hydraulic conductivity or transmissivity is considered by taking into account an isotropic exponential covariance function of log-transformed values with parameters describing the variance and correlation length. The method allows postprocessing of results from ground water models with respect to uncertainty estimation. The code PPPath, which was developed for this purpose, provides a postprocessing of pathline computations under PMWIN, which is based on MODFLOW. In order to test the methodology, it was applied to results from Monte Carlo simulations for catchments of pumping wells. The results correspond well. Practical applications illustrate the use of the method in aquifers.

Notes: Extremely alkaline ground water has been found underneath many shuttered steel mills and slag dumps and has been an impediment to the cleanup and economic redevelopment of these sites because little is known about the geochemistry. A large number of these sites occur in the Lake Calumet region of Chicago, Illinois, where large-scale infilling of the wetlands with steel slag has created an aquifer with pH values as high as 12.8. To understand the geochemistry of the alkaline ground water system, we analyzed samples of ground water and the associated slag and weathering products from four sites. We also considered several potential remediation schemes to lower the pH and toxicity of the water. The principal cause of the alkaline conditions is the weathering of calcium silicates within the slag. The resulting ground water at most of the sites is dominated by Ca2+ and OH− in equilibrium with Ca(OH)2. Where the alkaline ground water discharges in springs, atmospheric CO2 dissolves into the water and thick layers of calcite form. Iron, manganese, and other metals in the metallic portion of the slag have corroded to form more stable low-temperature oxides and sulfides and have not accumulated in large concentrations in the ground water. Calcite precipitated at the springs is rich in a number of heavy metals, suggesting that metals can move through the system as particulate matter. Air sparging appears to be an effective remediation strategy for reducing the toxicity of discharging alkaline water.

Notes: Artificial neural networks (ANNs) were developed to accurately predict highly time-variable specific conductance values in an unconfined coastal aquifer. Conductance values in the fresh water lens aquifer change in response to vertical displacements of the brackish zone and fresh water–salt water interface, which are caused by variable pumping and climate conditions. Unlike physical-based models, which require hydrologic parameter inputs, such as horizontal and vertical hydraulic conductivities, porosity, and fluid densities, ANNs can “learn” system behavior from easily measurable variables. In this study, the ANN input predictor variables were initial conductance, total precipitation, mean daily temperature, and total pumping extraction. The ANNs were used to predict salinity (specific conductance) at a single monitoring well located near a high-capacity municipal-supply well over time periods ranging from 30 d to several years. Model accuracy was compared against both measured/interpolated values and predictions were made with linear regression, and in general, excellent prediction accuracy was achieved. For example, although the average percent change of conductance over 90-d periods was 39%, the absolute mean prediction error achieved with the ANN was only 1.1%. The ANNs were also used to conduct a sensitivity analysis that quantified the importance of each of the four predictor variables on final conductance values, providing valuable insights into the dynamics of the system. The results demonstrate that the ANN technology can serve as a powerful and accurate prediction and management tool, minimizing degradation of ground water quality to the extent possible by identifying appropriate pumping policies under variable and/or changing climate conditions.

Notes: A new distance-drawdown method for aquifers with anisotropy on the horizontal plane is presented. The method uses scalar transformation to convert to an equivalent, isotropic medium, thus permitting application of the Cooper-Jacob Method. The method is applicable to cases where at least one ellipse of equal drawdown can be delineated but can also be applied where no ellipse can be discerned from the data. In the latter case, a least-squares regression approach can be employed to estimate the orientation and magnitude of the anisotropy. The regression R2 value provides a quantitative assessment of the degree to which the drawdown data are indicative of a systematic areal anisotropy in the aquifer or whether the data simply reflect natural aquifer heterogeneity. In addition to confined aquifers, this methodology, like the Cooper-Jacob Method, is also applicable to unconfined aquifers either before the onset of delayed drainage or following the completion of delayed drainage provided that the u value meets the recommended criterion.

Notes: A new multilayer approach is presented for the modeling of ground water flow to radial collector wells. The approach allows for the inclusion of all aspects of the unique boundary condition along the lateral arms of a collector well, including skin effect and internal friction losses due to flow in the arms. The hydraulic conductivity may differ between horizontal layers within the aquifer, and vertical anisotropy can be taken into account. The approach is based on the multilayer analytic element method, such that regional flow and local three-dimensional detail may be simulated simultaneously and accurately within one regional model. Horizontal flow inside a layer is computed analytically, while vertical flow is approximated with a standard finite-difference scheme. Results obtained with the proposed approach compare well to results obtained with three-dimensional analytic element solutions for flow in unconfined aquifers. The presented approach may be applied to predict the yield of a collector well in a regional setting and to compute the origin and residence time, and thus the quality, of water pumped by the collector well. As an example, the addition of three lateral arms to a collector well that already has three laterals is investigated. The new arms are added at an elevation of 2 m above the existing laterals. The yield increase of the collector well is computed as a function of the lengths of the three new arms.