ALBERT

Notes: Abstract Regionally distributed pelitic granulites in the Wilson Lake region contain the assemblage sapphirine + hypersthene + sillimanite + quartz. Geochronology and geobarometry suggest it developed in early Proterozoic rocks at temperatures approaching 900°C and pressures above 10 kbar. Vein-like metasomatized rocks around a suite of mafic to ultramafic intrusions, emplaced near the peak of metamorphism about 1700 Ma ago, contain sapphirine, but these assemblages developed at temperatures near 750°C and pressures of 4.5 kbar. Both types of assemblage occur as relics in amphibolite-grade (biotite–sillimanite) migmatites. P–T determinations indicate rapid isothermal uplift of 20 km accompanied by mafic intrusion and hydration. The metamorphic history and tectonic setting suggest exposure of deep continental crust by thrusting during continental collision, followed by essentially isothermal decompression.

Notes: Abstract In the contact metamorphic aureole of the Tinaroo Batholith (north Queensland, Australia), mylonitic rocks were metamorphosed during a regional folding/crenulation event (D2) synchronous with the emplacement of muscovite-bearing granitoids. Prismatic and skeletal andalusite porphyoblasts grew in carbonaceous schists, mainly from the dissolution of staurolite. Muscovite, quartz and biotite played a dual role in this reaction, acting in a catalytic capacity as well as reactants or products. Staurolite was replaced by coarse-grained muscovite ± biotite, whereas andalusite locally replaced quartz ± muscovite ± biotite, with diffusion of H, Al, Si, Mg, Fe and K ionic species linking sites of dissolution and growth.Graphite contributed to the reaction mechanism in a number of ways. Accumulations of graphite in front of advancing andalusite crystal faces led to skeletal growth and the formation of chiastolite structure, where incremental growth occurred on adjacent {110} faces, with subsequent filling in and inclusion of graphite along the diagonal zones. The presence of graphite in some layers in the schist matrix prevented recrystallization of strained muscovite grains. The muscovite grains in these layers, in contrast to adjacent thin non-graphitic layers, were preferentially replaced by quartz. This resulted in muscovite-depletion haloes in graphitic layers around andalusite porphyroblasts. Somewhat arcuate zones of graphite, concentrated during dissolution of quartz along a crenulation cleavage, occur on some andalusite faces. Reactivation of the mylonitic foliation during the formation of D2 crenulations led to a preferential dissolution of quartz in zones of progressive shearing localized near andalusite porphyroblasts and hence the accumulation of graphite.Lack of deflection of the pre-existing mylonitic foliation and anastomosing of the axial planes of D2 crenulations around andalusite porphyroblasts demonstrate not only the timing of growth, but also that growing porphyroblasts do not push aside existing foliations.

Notes: THE YOUNG EARTH: AN INTRODUCTION TO ARCHAEAN GEOLOGY. By E.G. Nisbet. Allen and Unwin, Boston, 1987. pp. 402.

Notes: Abstract Mica porphyroblasts in schists from several regions show nearly planar inclusion trails that are parallel over areas much larger than the wavelengths of later folds. This indicates that the porphyroblasts have not rotated, with respect to geographical co-ordinates, during deformation. Instead, the matrix has rotated, as suggested by Ramsay (1962). Even in zones of marked shortening in the matrix adjacent to large rigid porphyroblasts (e.g. of cordierite or staurolite), small biotite porphyroblasts have not rotated, but have become thinned by solution, as indicated by parallelism of inclusion trails in separate biotite grains and by evidence of truncation of inclusion trails by the matrix foliation. Less common are biotite porphyroblasts that have single asymmetrical microfolds in the matrix adjacent to the porphyroblasts and so appear to have rotated; these porphyroblasts are characterized by kinking.

Notes: Abstract Two major problems which exist in the use of illite crystallinity to define low-grade metamorphic zones are the variety of values chosen for the zone boundaries and the persistent use of three different indices of crystallinity. Although measurement techniques, which cause much of the interlaboratory variation, can be standardized, it is shown that there is, nevertheless, significant additional variation which demands calibration on standards. The greatest variations are due to choices of different absolute values of crystallinity to define zone boundaries. The problem of relating measurements between different indices is approached by fitting mathematical relationships to pairs of measurements from the same sample. A power–law relationship is a satisfactory fit to the Kubler–Weaver and Weaver–Weber pairs, while the Kubler–Weber indices are linearly related. These relationships are used to transform definitions of the diagenetic zone, anchizone and epizone from one index to the others, although they apply strictly only to the data set from which they are derived. This results in compatibility between the three zones and shows that previous definitions to the anchizone in different indices have been chosen at incompatible values. The boundaries of Kubler's anchizone (0.42 and 0.25 Δ2θ) are 0.4 and 0.215 Δ2θ in this study, which become 5.1 and 14.6 in the Weaver index and 278 and 149 in the Weber index. An error analysis shows that percentage errors in both Kubler and Weaver indices increase with crystallinity; the Kubler measurements are marginally preferred at all grades.

Notes: Abstract Migmatites in the Quetico Metasedimentary Belt contain two types of leucosome: (1) Layer-parallel leucosomes that grew during deformation and prograde metamorphism. These are enriched in SiO2, Sr, and Eu, but depleted in TiO2, Fe2O3, MgO, Cs, Rb, REE, Sc, Th, Zr, and Hf relative to the Quetico metasediments. (2) Discordant leucosomes that formed after the regional folding events when metamorphic temperatures were at their peak. These are enriched in Rb, Ba, Sr and Eu, but display a wide range of LREE, Th, Zr, and Hf contents relative to the Quetico metasediments.Layer-parallel leucosomes formed by a subsolidus process termed tectonic segregation. This stress-induced mass transfer process began when the Quetico sediments were deformed during burial, and continued whilst the rocks were both stressed and heterogeneous. Subsolidus leucosome compositions are consistent with the mobilization of quartz and feldspar from the host rocks by pressure solution. The discordant leucosomes formed by partial melting of the Quetico metasediments, possibly during uplift of the belt. The range of composition displayed by the anatectic leucosomes arises from crystal fractionation during leucosome emplacement. Some anatectic leucosomes preserve primary melt compositions and have smooth REE patterns, but those with negative Eu anomalies represent fractionated melts, and others with positive Eu anomalies represent accumulations of feldspar plus trapped melt.

Notes: Abstract The effects of varying amounts of partial melt on the deformation of granitic aggregates have been tested experimentally at conditions (900°C, 1500 MPa, 10-4 to 10-6/s) where melt-free samples deform by dislocation creep, with microstructures approximately equivalent to those of upper greenschist facies. Experiments were performed on samples of various grain sizes, including an aplite (150 μm) and sintered aggregates of quartz-albitemicrocline (10–50 and 2–10 μm). Water was added to the samples to obtain various amounts of melt (1–15% in the aplite, 1–5% in the sintered aggregates). Optical and TEM observations of the melt distribution in hydrostatically annealed samples show that the melt in the sintered aggregates is homogeneously distributed along an interconnected network of triple junction channels, while the melt in the aplites is inhomogeneously distributed.The effect of partial melt on deformation depends an melt amount and distribution, grain size and strain rate. For samples deformed with ˜ 1% melt, all grain sizes exhibit microstructures indicative of dislocation creep. For samples deformed with 3–5% melt, the 150 μm and 10–50 μm grain size samples also exhibit dislocation creep microstructures, but the 2–10 μm grain size samples exhibit abundant TEM-scale evidence of dissolution-precipitation and little evidence of dislocation activity, suggesting a switch in deformation mechanism to predominantly melt-enhanced diffusion creep. At natural strain rates melt-enhanced diffusion creep would predominate at larger grain sizes, although probably not for most coarse-grained granites.The effects of melt percentage and strain rate have been studied for the 150 μm aplites. For samples with ˜ 5 and 10% melt, deformation at 10–6/s squeezes excess melt out of the central compressed region allowing predominantly dislocation creep. Conversely, deformation at 10-5/s produces considerable cataclasis presumably because the excess melt cannot flow laterally fast enough and a high pore fluid pressure results. For samples with 15% melt, deformation at both strain rates produces cataclasis, presumably because the inhomogeneous melt distribution resulted in regions of decoupled grains, which would produce high stress concentrations at point contacts. At natural strain rates there should be little or no cataclasis if an equilibrium melt texture exists and if the melt can flow as fast as the imposed strain rate. However, if the melt is confined and cannot migrate, a high pore fluid pressure should promote brittle deformation.

Notes: Abstract Crystal-chemical relationships between coexisting sodic and calcic amphiboles have been studied in eclogitic metagabbros from the Aosta Valley, Western Alps. Textural analysis gives evidence of three successive high-pressure parageneses:1. Pre-kinematic high-grade blueschist assemblages, preserved as polymineralic inclusions in garnet cores and made of glaucophane and actinolite (stage A).2. Synkinematic eclogite assemblages, composed of garnet + omphacite + glaucophane ± actinolite ± white mica ° Clinozoisite + quartz + rutile (stage B).3. Post-kinematic epitactic overgrowths of barroisitic amphibole on glaucophane and actinolite (stage C). P–T conditions of the eclogitic metamorphism have been estimated at around 500–550°C, 16 kbar.Glaucophane and actinolite coexist as discrete grains in stage A and B assemblages. This texture and the chemistry of the amphiboles unambiguously denotes the existence of a miscibility gap between sodic and calcic amphiboles (from NaM4= 0.80 in actinolite to NaM4= 1.70 in glaucophane at T= 500–550°C). A comparison with published analyses allows a new solvus along the glaucophane–actinolite join to be drawn.The later barroisitic amphibole (stage C) exhibits strong chemical zonation indicating disequilibrium growth. This amphibole cannot either be used to define a miscibility gap with glaucophane or actinolite or be considered as an intermediate stage between these two end-members.

Notes: Abstract Garnet granulites from Sri Lanka preserve textural and chemical evidence for prograde equilibration at temperatures of at least 700–750°C and pressures in the vicinity of 6–8 kbar. Associated strain patterns suggest prograde metamorphism occurred during and immediately following an episode of crustal thickening, with the prograde P–T conditions probably reflecting a combination of the conductive and advective transport of heat at the mid-levels of tectonically thickened crust. The occurrence of prograde wollastonite provides evidence for internally buffered fluid compositions, or fluid absent conditions, during peak metamorphism and precludes pervasive advection of a CO2-rich fluid. The advective heat component is therefore likely to have been provided by the transport of silicate melt. Intricate symplectitic textures record partial re-equilibration of the garnet granulites to lower pressures (˜ 4–6 kbar) at high temperatures (600–750°C), and testify either to the erosional denudation of the overthick crust prior to significant cooling (i.e. quasi-isothermal decompression) or to a subsequent static heating possibly of early Palaeozoic age (Pan-African). The metamorphic history of the Sri Lankan granulites is compared with high grade terrains in the neighbouring fragments of Gondwana, with the emphasis on similarities with Proterozoic granulites of the East Antarctic craton.

Notes: TECTONIC SETTINGS OF REGIONAL METAMORPHISM. Edited by E.R. Oxburgh, B.W.D. Yardley and P.C. England. Royal Society of London. 1987.

Notes: Abstract This paper provides methods and a description of a Pascal computer program, thermocalc, for various thermodynamic calculations using the thermodynamic dataset presented in earlier papers in this series (Holland & Powell, 1985; Powell & Holland, 1985). The dataset involves uncertainties on the thermodynamic parameters and therefore allows uncertainties to be calculated on results, for example in geothermometry and geobarometry. Recommendations are made for the uncertainties on activities to be used in calculations on rocks, particular emphasis being placed on preventing underestimates of these uncertainties at small mole fractions. Apposite examples of phase diagram and rock calculations are presented with ouput from thermocalc, demonstrating the utility of the program. Of the rock calculations, the most valuable are considered to be those involving simultaneous combination ‘least squares’of calculated conditions for a set of reactions applicable to a rock. This set of reactions involves the independent reactions which can be written between the end-members in the minerals in a rock and in the thermodynamic dataset. In contrast to an approach based on specific geothermometers and geobarometers, this approach maximizes the benefit of having an internally consistent thermodynamic dataset. thermocalc is available in IBM PC and Mac versions, from Roger Powell for A$25 or Tim Holland for £10 per version.

Notes: Abstract An inverted metamorphic gradient associated with the northern mylonite zone of the Cheyenne belt, a deeply eroded Precambrian suture in southern Wyoming, has been documented within metasedimentary rocks of the Early Proterozoic Snowy Pass Supergroup. Metamorphic grade in the steeply dipping supracrustal sequence increases from the chlorite through the biotite, garnet, and staurolite zones both stratigraphically and structurally upward toward the northern mylonite zone. A minimum temperature increase of approximately 100° C over a km-wide zone is required for this transition. Parallelism of inverted isograds with the trace of the northern mylonite zone implies a genetic relationship between deformation associated with that zone and the inverted metamorphic gradient within the Snowy Pass Supergroup.Field evidence together with microstructural and petrofabric analysis indicate northward thrusting of amphibolite-grade rocks over rocks of the Snowy Pass Supergroup along the northern mylonite zone. Mineral equilibria and garnet-biotite geothermometry on synkinematic mineral assemblages within the Snowy Pass metasedimentary rocks indicate deformation at minimum temperatures of 480° C and pressures of 350–400 MPa (3°5–4°0 kbar). This implies tectonic burial or upper plate thickness of 13–15 km.The narrow character of metamorphic zonation and microtextures within the Snowy Pass Supergroup which indicate late synkine-matic growth of garnet and staurolite, preclude rotation of pre-existing isograds by folding as a mechanism for development of the inverted gradient. Conductive transport of heat from the upper into the lower plate across the originally low-angle thrust is insufficient to produce the necessary temperatures in the lower plate. Shear heating is considered insufficient to produce the observed metamorphic transition unless high shear stresses are postulated. Up-dip advection of metamorphic fluids is a feasible, but unproven, mechanism for heat transport. The possibility that rapid uplift due to stacking of several thrust sheets may have played a role in preserving the inverted metamorphic gradient cannot be evaluated at present.

Notes: Abstract Partial melting of tonalitic gneisses in the 2.7 Ga Badcallian granulite facies metamorphic episode in the Scourian complex of north-west Scotland produced a suite of granitic to trondhjemitic liquids. On cooling and excavation of the complex, these melts underwent fractional crystallization and the residual liquids eventually became water saturated. Comparison with experimental data suggests that water saturation would have occurred in these melts at around 620–700°C. From the retrograde P–T-time path followed by the complex it is estimated that H2O-dominated fluids were exsolved from these melts at c. 2.5 Ga. It is proposed that these fluids were the cause of the 2.5 Ga Inverian retrogression of the Scourian complex and that water-saturated melts formed during the crystallization of the leucogneisses were intruded as a suite of pegmatites. The timing of pegmatite intrusion is consistent with this proposition as are the temperature estimates, timing, distribution and nature of the Inverian phase of metamorphism. It is likely that the crystallization of melts is an important process in bringing about hydrous retrogressive metamorphic episodes in a number of other basement terrains, such as West Greenland and Australia.

Notes: Abstract Partitioning of Fe and Mg between garnet and phengitic muscovite was calibrated as a geothermometer by Green & Hellman (1982) using experimental data at 25–30 kbar. When the thermometer is applied to pelites regionally metamorphosed at pressures of between 3 and 7 kbar it yields temperatures much higher than those from the garnet–biotite thermometer. A new empirical calibration is proposed for use with such rocks, with particular application where garnet occurs at lower grades than biotite. The new calibration is where K is given by: In K= In Kd and Xii are mole fractions in the garnets.The calibration was derived from comparison with the garnet–biotite thermometer of Ferry & Spear (1978), assuming no pressure-dependence for the partitioning between garnet and muscovite, no ferric iron partitioning, ideal mixing in muscovite, and the garnet mixing model of Ganguly & Saxena (1984) modified for a non-linear Ca effect. This latter garnet mixing model was selected because it gave the geologically most reasonable results. It has not proved possible to distinguish a pressure effect from a ferric-iron effect.Despite the simplifying assumptions used to derive the calibration, it yields temperatures generally within 15°C of those given by the garnet–biotite thermometer, and has been used to supply thermometric data in a low-grade region of the Canadian Rockies.

Notes: Abstract Mineral assemblages in different samples of amphibolite facies pelitic schists collected from two separate outcrops in the Moosilauke area, NH, record differences in the chemical potential of water during metamorphism. Mineralogical, petrological, and field relations indicate that mineral assemblages at both outcrops equilibrated at 520°C and 3.5–4.0 kbar. Thermodynamic analysis of the mineral assemblages demonstrates that maximum chemical potential differences at each outcrop were of the order of 150 calories, over distances of 10–20 m.The differences in the chemical potential of water recorded in both bed-to-bed and outcrop-to-outcrop relations are consistent with the following conclusions: (1) mineral assemblages on a specific outcrop did not equilibrate with an external reservoir of fluid of fixed composition, (2) the relatively small magnitude of the chemical potential differences suggests little or no infiltration of externally derived fluid, (3) these differences on the outcrop scale are probably related to initial compositional variations and the buffer capacity of the mineral assemblage, and (4) the different values of the chemical potential of water exhibited by the various mineral assemblages permits an understanding of the effects of variable μH2O for amphibolite facies pelitic schists.

Notes: Abstract A detailed study of garnet–chloritoid micaschists fom the Sesia zone (Western Alps) is used to constrain phase relations in high pressure (HP) metapelitic rocks. In addition to quartz, phengite, paragonite and rutile, the micaschists display two distinct parageneses, namely garnet + chloritoid + chlorite and garnet + chloritoid + kyanite. Talc has never been observed. Garnet and chloritoid are more magnesian when chlorite is present instead of kyanite. The distinction of the two equilibria results from different bulk rock chemistries, not from P–T conditions or redox state. Estimated P–T conditions for the eclogitic metamorphism are 550–600°C, 15–18 kbar.The presence of primary chlorite in association with garnet and chloritoid leads us to construct two possible AFM topologies for the Sesia metapelites. The paper describes a KFMASH multisystem for HP pelitic rocks, which extends the grid of Harte & Hudson (1979) towards higher pressures and adds the phase talc. Observed parageneses in HP metapelites are consistent with predicted phase relations. Critical associations are Gt–Ctd–Chl and Gt–Ctd–Ky at relatively low temperatures and Gl–Chl–Ky and Gt–Tc–Ky at relatively high temperatures.

Notes: Abstract Low-pressure prograde metamorphism of pelitic rocks in the Cooma Complex, south-east Australia, has produced cordierite-andalusite schists at intermediate grades. The first foliation (S1) is preserved largely as inclusion trails in cordierite porphyroblasts. Microstructural evidence indicates that the cordierite porphyroblasts grew during the early stages of development of a crenulation-foliation (S2) and that andalusite porphyroblasts grew during the development of a later crenulation-foliation (S3). Microstructural evidence also indicates that the andalusite was a product of the prograde reaction: cordierite + muscovite ± andalusite + biotite + quartz. The occurrence of the products of this reaction in ‘beard’structures between cordierite microboudins formed by extension in S3 confirms that the andalusite grew during the development of S3. The investigation shows that porphyroblast-matrix relationships can preserve the orientation of an early S-surface that has been largely obliterated from the matrix, as well as providing relatively direct evidence of sequential mineral growth and metamorphic reactions.

Notes: Abstract In the northeastern part of the Grenville Province, along the gulf of St Lawrence, cordierite is widespread in the migmatites of Baie Jacques Cartier (BJC) and Baie des Ha! Ha! (BHH). In the BJC area, rafts of mesosome occur in a pervasive network of leucosome consisting of cordierite-bearing pegmatite. In BHH, however, the mesosome and leucosome are well segregated and locally separated by thin biotite –hornblende melanosomes.Leucosomes in the BJC area record the highest temperatures (oxide thermometry = 900°C), whereas leucosomes of BHH and mesosomes of both areas indicate peak temperatures around 800°C (oxide thermometry; biotite–garnet thermometry with fluorine-rich biotite). Peak pressures were constrained at 720 MPa using the Ilm-Sil–Qtz–Grt–Rt (GRAIL) equilibrium.The area is thought to have undergone extensive melting under relatively modest pressures. The highest temperatures recorded in the BJC area are probably related to a pervasive impregnation of this terrane by aluminous granitic melts.Most post-peak P–T estimates for the mesosomes fall on a nearly isobaric, clockwise, P–T path (0.6 MPa/°C) with the exception of the high-temperature leucosomes of BJC, which fall about 100°C away from this path; this is additional evidence for the external origin of these leucosomes. The ultimate source of heat that generated the migmatites is thus though to be an underlying plutonic complex (anorthosite?).

Notes: Abstract The introduction of externally derived fluids into rocks of the Zermatt–Saas zone of the Swiss Alps gave rise to the simultaneous formation of shear and hydraulic fractures. These fractures are now filled with albite-rich assemblages and surrounded by alteration halos up to c. 2 m wide. The alteration assemblages are zoned and an examination of reactions in P–T–aH2O space implies that the parageneses developed by the hydration of fluid-absent eclogites. A mechanical analysis of the veins (after Sibson, 1981) shows that Pfluid/Pload must have been at least 0.96. Fluid migration into the country rocks must have been driven by excess hydraulic head either derived from the vertical extent of the veins or due to their connection to a deeper, external reservoir, possibly tapped along thrust surface(s). Diffusive and capillary transport were insignificant. The fluids may have been derived from underlying metasediments that were dehydrating during the quasi-isothermal uplift of this part of the Alps, or they may have originated during the prograde mesoalpine metamorphism documented in the area.

Notes: Abstract The Catalina Schist of southern California is a subduction zone metamorphic terrane. It consists of three tectonic units of amphibolite-, high-P greenschist- and blueschist-facies rocks that are structurally juxtaposed across faults, forming an apparent inverted metamorphic gradient. Migmatitic and non-migmatitic metabasite blocks surrounded by a meta-ultramafic matrix comprise the upper part of the Catalina amphibolite unit. Fluid-rock interaction at high-P, high-T conditions caused partial melting of migmatitic blocks, metasomatic exchange between metabasite blocks and ultramafic rocks, infiltration of silica into ultramafic rocks, and loss of an albitic component from nonmigmatitic, clinopyroxene-bearing metabasite blocks.Partial melting took place at an estimated P=˜8–11 kbar and T=˜640–750°C at high H2O activity. The melting reaction probably involved plagioclase + quartz. Trondhjemitic melts were produced and are preserved as leucocratic regions in migmatitic blocks and as pegmatitic dikes that cut ultramafic rocks.The metasomatic and melting processes reflected in these rocks could be analogous to those proposed for fluid and melt transfer of components from a subducting slab to the mantle wedge. Aqueous fluids rather than melts seem to have accomplished the bulk of mass transfer within the mafic and ultramafic complex.

Notes: Abstract Migmatite structures in the Coast Plutonic-Metamorphic Complex are well exposed in the inlet of Boca de Quadra, southeast Alaska. Two types of anatectic migmatites are present. Patch migmatites formed by in situ melting and subsequent crystallization of melt. Diktyonitic migmatites comprise a discontinuous veined network of leucocratic material, in which leucosomes enclose boudins of host rock. The margins of these boudins show the development of both melanosomes and shear band fabrics.Strain analysis of diktyonitic melanosomes indicates that these regions have undergone volume decreases of 20-27%. This volume decrease is attributed to melt extraction into the adjacent fracture-filling leucosomes. Thus, diktyonitic migmatites formed by shear-induced segregation of partial melt, whereas in patch migmatites the lack of shear stresses inhibited melt segregation. The variable structural style of anatectic migmatites in Boca de Quadra is not related to host-rock composition, but may be due to differences in the amount of differential stress during migmatization. These in turn may be controlled by host-rock strength and/or diachroneity of migmatization and deformation.Determination of volume changes during migmatization using strain analysis is potentially capable of discriminating intrusive and anatectic migmatites and consequently of documenting melt segregation and subsequent migration across crustal levels.

Notes: Abstract Chloritoid and pyrophyllite occur together in all major goldfields of the Witwatersrand Basin and are widespread in virtually all rock types of the upper Witwatersrand Supergroup, including metaconglomeratic reefs and altered mafic rocks. Both minerals are particularly characteristic of the pelitic horizons intimately associated with reef packages, but they are also developed locally in the regionally persistent metapelites that have basin-wide extent. Pyrophyllite is particularly common in foliated zones, adjacent to quartz veins, and near unconformably overlying auriferous conglomerates.The wide distribution of chloritoid and pyrophyllite in metapelites of the Witwatersrand Basin is attributed to alteration of chlorite-rich shales, rather than to unusual premetamorphic starting materials. This alteration event involved the redistribution of many elements, with up to 40% volume loss, mainly due to removal of silica. Removal of most of the Mg and some Fe accounts for the stabilization of chloritoid and pyrophyllite. Relatively immobile elements included Al, Ti, Nb, Cr, V, P, La and Ce, whereas Si, Fe, Mn, Zn, Co, Ni, Cu, Mg and Ca were lost, and K, Rb and Ba were introduced by an infiltrating fluid.The alteration event is inferred to have been within the chloritoid and pyrophyllite stability field (and thus syn-metamorphic) as bulk chemical changes in metapelites are from chlorite directly towards chloritoid and then pyrophyllite, rather than to lower grade minerals such as kaolinite. Muscovite–chlorite–chloritoid and muscovite–chloritoid–pyrophyllite assemblages are attributed to fluid buffering along appropriate curves, as their production by metamorphism of lower grade mineral mixes is considered unlikely, based on the present bulk rock compositional data. A metamorphic timing for the alteration accounts for the correlation of strongly foliated areas with greater degrees of inferred alteration. The transitions from chlorite to chloritoid to pyrophyllite define zones of increasing alteration.Widespread infiltration as part of peak metamorphism is suggested by the distribution of chloritoid and pyrophyllite, quartz veining and textures. Fluid:rock ratios calculated from a silica budget in one metapelitic horizon exceed 100:1 over many square kilometres. These values need not imply multi-pass fluid flow, as much of the silica migration may be redistribution on a scale of a few metres, from source rocks into veins. Although infiltration during metamorphism may have affected much of the upper Witwatersrand succession, channelized fluid flow within reef packages, along faults and unconformities and in certain metaconglomerates and metapelites is inferred.

Notes: Abstract Biotite and cordierite occur in a 1-km wide zone of pelitic hornfelses around the McGerrigle pluton. These phases display systematic changes in XFe that can be attributed to continuous reactions involving chlorite or andalusite in the system KFMASH. Through much of the zone biotite and cordierite were products of the ‘breakdown’of chlorite. Close to the pluton this continuous reaction was terminated by a discontinuous reaction that introduced andalusite. Pelites which interdigitate with apophyses of the intrusive at the pluton margin contain assemblages that record a continuous reaction between biotite, cordierite, andalusite, muscovite, and quartz or, alternatively, the discontinuous breakdown of muscovite and quartz to K-feldspar and andalusite.The mole fraction of Fe in biotite and cordierite increased significantly with the progress of the first continuous reaction and apparently decreased during the second continuous reaction. The KD of Fe-Mg between the minerals decreased and apparently increased, respectively, during the two reactions.Biotite-cordierite-chlorite assemblages are interpreted to have been stable at temperatures between 525° C and 615° C and biotite-cordierite-andalusite assemblages stable at temperatures between 615° C and 635° C. The confining pressure was estimated to have been 〈 2 kbar.The results of this study suggest that the KD of Fe-Mg between biotite and cordierite is a function of temperature, the Fe-Mg exchange characteristics of the controlling continuous reaction and non-ideal mixing of Fe and Mg.

Notes: Abstract The textural and compositional features of phlogopites in a contact-metamorphic dolomite marble inclusion in the Bergell intrusion (central Alps) and in a metasomatic reaction vein cutting through this marble suggest different origins for vein phlogopites:(a) High-Al vein phlogopite represents former marble phlogopite which has been compositionally modified by reaction with the vein forming fluid.(b) Low-Al vein phlogopite represents phlogopite precipitated from the vein forming fluid.As both types of vein phlogopite were in contact with the same vein forming fluid at the same time, low-Al phlogopite most likely represents an equilibrium phlogopite composition, whereas high-Al phlogopite does not. High-Al vein phlogopite retained its Al-content from the contact-metamorphic marble parent phlogopite and only underwent Fe-Mg exchange with the metasomatic fluid.All the vein phlogopites studied are strongly enriched in Fe relative to marble phlogopite. The data may suggest in general that phlogopite Al/Si ratios may be retained from the conditions under which the phlogopites first formed, whereas the Mg/Fe-ratios may be substantially modified by exchange with other ferromagnesian solid phases and/or a metamorphic fluid at later stages in their metamorphic history. This may have significant effects on calculated pressures and temperatures from thermobarometers involving biotite.

Notes: Abstract The Saint-Martin des Noyers Formation is interpreted as a slice of an island-arc system of Lower to Middle Palaeozoic age, located in the internal part of the Variscan orogen in Vendée (Armorican Massif, France). Metamorphosed igneous rocks range in composition from ultramafic to rhyolitic. The regular increase in the FeO/(FeO+MgO) ratio, from mafic to silicic samples, results in a systematic variability in the nature and composition of the metamorphic phases. In basaltic samples, the occurrence of relict garnet-barroisite assemblages suggests relatively high-pressure conditions for the peak of metamorphism. During a subsequent retrograde evolution, the primary barroisitic hornblendes recrystallized to texturally complex mixtures of actinolite and hornblende. Despite this complication, it is possible to decipher a P–T-t path based on amphibole chemistry. The P–T trajectory deduced is dominated by the effect of pressure and consistent with early underthrusting and subsequent tectonic uplift of the ancient arc of Saint-Martin des Noyers.

Notes: Abstract In order to study the thermal structure of active thrust belts, we have developed a numerical model of conductive heat transfer between thrust sheets during deformation. Our finite difference approach alternates small, instantaneous increments of displacement and isotherm translation with conductive relaxation of perturbed isotherms. In each step, conduction occurs for a length of time equal to the displacement increment divided by the thrust velocity. Computer simulations demonstrate that conductive heat transfer is significant during deformation and that temperatures in hanging-wall rocks decrease while temperatures in foot-wall rocks increase over distances of up to 10 km from the thrust surface. When the effects of internal heat production are also calculated, heating of foot-wall rocks exceeds cooling of hanging-wall rocks. Rocks located between two thrusts may experience a complicated temperature–time path of early heating followed by cooling. These models help to explain the rapid metamorphism of rocks in the Taconian thrust belt in the northern Appalachians of New England soon after deposition of the youngest sediments.

Notes: Abstract 40Ar/39Ar age spectrum analysis of phengite separates from Naxos, part of the Attic Cycladic Metamorphic Belt in Greece, indicates that cooling following high-pressure, low- to medium-temperature metamorphism, M1, occurred about 50 Ma ago. Phengite has 40Ar* gradients that suggest that part of the scatter observed in conventional K–Ar ages was caused by diffusion of radiogenic argon from the minerals during a younger metamorphism, M2. In central Naxos, this metamorphism (M2) has overprinted the original mineral assemblages completely, and is associated with development of a thermal dome. Excellent 40Ar/39Ar plateaus at 15.0 ± 0.1 Ma, 11.8 ± 0.1 Ma, and 11.4 ± 0.1 Ma, obtained on hornblende, muscovite and biotite, respectively, from the migmatite zone, indicate that relatively rapid cooling followed the M2 event, and that no significant thermal overprinting occurred subsequent to M2. Toward lower M2 metamorphic grade, 40Ar/39Ar plateau ages of hornblendes increase to 19.8 ± 0.1 Ma; concomitantly the proportion of excess 40Ar in the spectra increases as well. We propose that the peak of M2 metamorphism occurred beween 15.0 and 19.8 Ma ago. K–Ar ages of biotites from a granodiorite on the west coast are indistinguishable from those found in the metamorphic complex, and hornblende K–Ar ages from the same samples are in the range 12.1–13.6 Ma. As the latter ages are somewhat younger than most ages obtained from the metamorphic complex, intrusion of the granodiorite most likely followed the peak of the M2 metamorphism.The metamorphic evolution of Naxos is consistent with rapid crustal thickening during the Cretaceous or early Tertiary, causing conditions at which supracrustal rocks experienced pressures in the range 900–1500 MPa. Transition to normal crustal thicknesses ended the M1 metamorphism about 50 Ma ago. The M2 metamorphism and granodiorite intrusion occurred during a period of heat input into the crust, possibly related to the migration of the Hellenic volcanic ar°C in a southerly direction through the area.

Notes: Abstract The layers of six stromatic migmatites from Northern, Western, and Central Europe display small but systematic chemical and mineralogical differences. At least five of these migmatites do not show any signs of largescale metamorphic differentiation, metasomatism, or segregation of melts. It is concluded, therefore, that the compositional layering observed in most of the investigated migmatites is due to compositional differences inherited from the parent rocks. Almost isochemical partial melting seems to be the most probable process transforming layered paragneisses, metavolcanics, or schists into migmatites.The formation of neosomes is believed to be caused by higher amounts of partial melts formed due to higher amounts of water moving into these layers. The neosomes have less biotite and more K-feldspar, if K-feldspar is present at all, than the adjacent mesosomes. These differences are small but systematic and seem to control the access of different amounts of water to the various rock portions. Petrographical observations, chemical data, and theoretical considerations indicate a close relationship between rock composition, rock deformation, transport of water, partial melting, and formation of layered migmatites.

Notes: Abstract Anatectic migmatites of contrasting structural style are found adjacent to the contacts of the Ballachulish Igneous Complex, Argyllshire, Scotland. On the east flank, evidence for migmatization is largely restricted to the local development of millimetre-centimetre scale Kfs + Qtz-rich leucocratic segregations, which accompany fragmentation of brittle hornfels layers and ductile deformation of mm-cm scale semipelitic layers. Large volumes of semipelitic rock rich in feldspar and quartz on the east flank show no migmatitic features, and bedding is usually preserved undisturbed right up to the contact. On the west flank, in contrast, similar semipelitic rocks show widespread migmatitic features and disruption of layering is substantial and widespread over a 400 m wide zone. Within the west-flank migmatites, 1–100 cm scale rigid bedding fragments (schollen) may be suspended and disoriented in a semipelitic matrix that underwent ductile deformation. The P-T conditions on both flanks are in the same range: 3 kbar and 650–700°C.The contrast in gross structural style is believed to result from differences in the volumes of melt produced and differences in the proportion of rock in which the critical melt fraction of the rocks was exceeded. On the east flank, only on a mm-cm scale was enough melt locally accumulated to cause disruption of some layers and segregation of melt. On the west flank, melting proceeded substantially in a broad tract of semipelitic rocks, resulting in larger scale contrasts in rheology that led to the present chaotic structures in this zone.Because migmatization occurred at a pressure too low for muscovite dehydration melting, and at temperatures too low for substantial biotite dehydration melting, the different amounts of melting on the east and west flanks most probably resulted from the introduction of differing amounts of externally derived water. On the east flank, and throughout most of the aureole, the absence of melting even in quartzofeldspathic protoliths indicates that there was no substantial movement of fluid towards or away from the igneous complex during migmatization. The contrasting situation on the west flank may have resulted from devolatilization of underlying quartz diorite magma (˜ 690–710°C), which released heat and fluids into the overlying quartz- and feldspar-rich semipelites (solidus temperature ˜ 650–680°C).

Notes: Abstract Edenite/tremolite and edenite/magnesio-hornblende in equilibrium with plagioclase, chlorite, epidote, quartz and vapour involve several types of reactions for which KD can be related to T and P. Thermodynamic calculation of these equilibria leads to isopleth systems. Given knowledge of the progressive changes of end-member activities in zoned Ca–Mg amphiboles (based on microprobe analyses), it is possible to construct precise pressure–temperature–time paths (P–T–t paths) which have been followed by metabasites during polyphase metamorphism. When applied to basic rocks from the River Vilaine area, this method allows us to construct a P–T–t path that can be compared directly to the P–T–t path constructed from interbedded acid rocks (aluminous micaschists) in the same structural unit. Through time, both basic and acid rocks underwent the same complex deformation history that can be described conveniently in the L–S fabric system of Flinn. This allows us to construct a P–T–t deformation path for this structural unit.These paths are interpreted in terms of an under/overthrusting continental collision belt (the Hercynian belt), and represent an illustration of the time delay caused by stacking of more than two crustal units.

Notes: Abstract The Western Baja terrane (WBt) of west-central Baja California is an uplifted subduction complex that is divided into smaller ‘subterranes’on the basis of bounding faults and petrological differences. Each subterrane contains coherent Early Jurassic to Early Cretaceous sedimentary and mafic volcanic rocks (not melange) that have been metamorphosed under blueschist facies conditions. Key phases in metabasites and metaturbidites include jadeitic to acmitic clinopyroxene, sodic amphibole, lawsonite, aragonite, chlorite, titanite and white mica. Pressure indicators include the jadeite content of clinopyroxene and the presence of aragonite. Temperature indicators include the presence of lawsonite, the absence of greenschist facies minerals and results from vitrinite reflectance studies. Conditions at the peak of metamorphism were 〉8 kbar, 225–325°C for subterrane 1, 7–8 kbar, 170–220°C for subterrane 2, and 5–6 kbar, 175–200°C for subterrane 3; these correspond to cold geothermal gradients (6–9/km). Vein assemblages that include aegerine–jadeite and aegerine, albite, aragonite, lawsonite and sodic amphibole indicate uplift during continued cold conditions, probably during steady-state subduction.

Notes: LOW TEMPERATURE METAMORPHISM. Edited by M. Frey. Blackie & Son Limited, Glasgow and London. 1987. pp. 364.

Notes: Abstract Endogreisens which replace K-feld-spar-quartz dykes in a Devonian (360 Ma) tin deposit at Mt Bischoff, north-west Tasmania, formed from the interaction of unusual solutions, probably derived from an underlying leucogranite pluton, porphyry dykes and limited quantities of local dolomitic country rock components. The intensity of greisenization and pH of the solutions increase inward to the greisenized dykes’cores and downward. The following types of greisen assemblages indicate increasing degrees of greisenization: ‘sericite’muscovite + quartz ± tourmaline ± fluorite, topaz + quartz ± tourmaline ± fluorite, weberite, prosopite, ralstonite, Ca-ralstonite; and quartz ± topaz ± fluorite. Where the solutions interacted with dolomite, exogreisens consisting of topaz- or tourmaline-bearing assemblages were formed. The greisens were subsequently overprinted to varying degrees by siderite, sulphides and hydrous silicates (talc, serpentine, chlorite, micas).The temperature during greisenization ranged from 180 to 414°C, based on fluid inclusions in topaz, quartz, fluorite, sellaite and cassiterite. The main greisen-forming event occurred at temperatures of 360±20°C. The fluids boiled intermittently. Their salinities ranged from 31.5 to 38.9 wt% total dissolved salts, consisting of Ca–K–Na–Fe–Cl±hydrocarbon species. Fluid inclusion data indicate that only 0.5–1.5 km of cover were present above this deposit at the time of formation.The greisenized dykes were intruded by and intrude different stages of breccias. The breccias consist mainly of country rock and greisenized dyke fragments, with rock-flour and later tourmaline alteration. The Mt Bischoff greisen system is possibly part of a ‘porphyry tin’style deposit formed at near-surface conditions (0.5–1.0 km).

Notes: Abstract Dehydration-melting reactions, in which water from a hydrous phase enters the melt, leaving an anhydrous solid assemblage, are the dominant mechanism of partial melting of high-grade rocks in the absence of externally derived vapour. Equilibria involving melt and solid phases are effective buffers of aH2,o. The element-partitioning observed in natural rocks suggests that dehydration melting occurs over a temperature interval during which, for most cases, aH2o is driven to lower values. The mass balance of dehydration melting in typical biotite gneiss and metapelite shows that the proportion of melt in the product assemblage at T± 850°C is relatively small (10–20%), and probably insufficient to mobilize a partially melted rock body.Granulite facies metapelite, biotite gneiss and metabasic gneiss in Namaqualand contain coarse-grained, discordant, unfoliated, anhydrous segregations, surrounded by a finer grained, foliated matrix that commonly includes hydrous minerals. The segregations have modes consistent with the hypothesis that they are the solid and liquid products of the dehydration-melting reactions: Bt + Sil + Qtz + PI = Grt ° Crd + Kfs + L (metapelite), Bt + Qtz + Pl = Opx + Kfs + L (biotite gneiss), and Hbl + Qtz = Opx + Cpx + Pl + L (metabasic gneiss). The size, shape, distribution and modes of segregations suggest only limited migration and extraction of melt. Growth of anhydrous poikiloblasts in matrix regions, development of anhydrous haloes around segregations and formation of dehydrated margins on metabasic layers enclosed in migmatitic metapelites all imply local gradients in water activity. Also, they suggest that individual segregations and bodies of partially melted rock acted as sinks for soluble volatiles. The preservation of anhydrous assemblages and the restricted distribution of late hydrous minerals suggest that retrograde reaction between hydrous melt and solids did not occur and that H2O in the melt was released as vapour on crystallization.This model, combined with the natural observations, suggests that it is possible to form granulite facies assemblages without participation of external fluid and without major extraction of silicate melt.

Notes: Abstract Nearly pure CO2 fluid inclusions are abundant in migmatites although H2O-rich fluids are predicted from the phase equilibria. Processes which may play a role in this observation include (1) the effects of decompression on melt, (2) generation of a CO2-bearing volatile phase by the reaction graphite + quartz + biotite + plagioclase = melt + orthopyroxene + CO2-rich vapour, (3) selective leakage of H2O from CO2+ H2O inclusions when the pressure in the inclusion exceeds the confining pressure during decompression, and (4) enrichment of grain-boundary vapour in CO2 by subsolidus retrograde hydration reactions.

Notes: Abstract The Hidaka metamorphic terrane in the Meguro-Shoya area, Hokkaido, Japan is divided into four progressive metamorphic zones: A—biotite zone; B—cordierite zone; C—cordierite–K-feldspar zone; and, D—sillimanite–K-feldspar zone of the andalusite–sillimanite facies series type of metamorphism. The metamorphic grade ranges from the higher temperature part of the greenschist facies (zone A) through the amphibolite facies (zones B and C) to the lower temperature part of the granulite facies (zone D). The zone boundaries intersect the bedding planes at high angles. P–T conditions estimated are 450–550°C and 2 kbar for zone A, 550–600°C and 2–2.5 kbar for zone B, 600–650°C and 2.5–3 kbar for zone C and 650–750°C and 3–4 kbar for zone D. The metapelites of zone D were partially melted.At the later stage of the regional metamorphism which is early Oligocene to early Miocene in age, cordierite tonalite and biotite tonalite intrusives associated with segments of the highest grade rocks (zone D) were emplaced into the lower temperature part of the regional metamorphic rocks, giving rise to a contact metamorphic aureole. The thermally metamorphosed terrain (zone C') belongs to the amphibolite facies and its P–T conditions are estimated to have been 550–700°C and 2 kbar.The P–T–t paths of the Hidaka metamorphism show a thickening–heating–uplifting process. The metamorphism is inferred to have taken place beneath an active island arc accompanied by partial melting of the crust.

Notes: Abstract Textural relationships between porphyroblasts of biotite and garnet in metasediments in the Nordkinn Peninsula area of the Finnmarkian Caledonides of North Norway are apparently complex. There is evidence for two textural zones in both mineral phases and superficially the development of these appears to have overlapped, at least in part, in time and space. This apparently complex porphyroblast growth history can be considerably simplified if only one period of garnet growth occurred and if different inclusion fabrics developed where garnet replaced biotite porphyroblasts and where it overgrew the matrix foliation. The possibility that porphyroblasts with textural evidence for multiphase growth histories actually grew during a single crystallization event is of importance in the interpretation and elucidation of tectonometamorphic relationships.

Notes: Abstract Retrograde and prograde mineral assemblages from metapelitic and metabasic rocks of the Iforas Granulitic Unit (Mali) were generated by the superimposition of two granulite facies metamorphic events. They clearly result from a polycyclic evolution and can be related to a late Eburnean unroofing followed by a Pan-African burial.Thermobarometry on Pan-African garnet-bearing assemblages yields (P, T) estimates of 620±50°C and 5± Ikbar. The nearly anhydrous conditions produced in the Eburnean appear to be the direct cause of the unusually lowtemperature granulite-facies metamorphism in the Pan-African. These P, T estimates are compared with those obtained on the underlying unit (Kidal Assemblage) upon which the Iforas Granulitic Unit was thrust. A P-T-t path, during the Pan-African orogeny, is proposed and discussed for both the Iforas Granulites and Kidal Assemblage.

Notes: Abstract Fluids, some of which are CO2-rich (up to 40 mol.% CO2) and some of which are highly saline (up to 18 wt% NaCl equivalent), are trapped as fluid inclusions in quartz-calcite (∼ metallic minerals) veins which cross-cut the pumpellyite-actinolite to amphibolite facies rocks of the Alpine Schist. Fluids were commonly trapped as immiscible liquid-vapour mixes in quartz and calcite showing open-space growth textures. Fluid entrapment occurred at fluid pressures near 500 bars (possibly as low as 150 bars) at temperatures ranging from 260 to 330° C. Saline fluids may have formed by partitioning of dissolved salts into an aqueous phase on segregation of immiscible fluids from a low-density CO2-rich fluid. Calcite deposited by these fluids has δ13C ranging from – 8.4 to – 11.5 and δ18O from + 4 to + 13. Isotopic data, fluid compositions and mode of occurrence suggest that the fluids are derived from high-grade metamorphic rocks. Fluid interaction with wall-rock has caused biotite crystallization and/or recrystallization in some rocks and retrogression of biotite to chlorite in other rocks.Fluid penetration through the rock is almost pervasive in many areas where permeability, probably related to Alpine Fault activity, has focussed fluids on a regional scale into fractured rocks. The fluid flow process is made possible by high uplift-rates (in excess of 10 mm/year) bringing hot rocks near to the surface.

Notes: Abstract The Tal y Fan Intrusion is a 110 m thick sub-concordant metabasite sheet intruded into volcaniclastic and pyroclastic rocks of Ordovician age in North Wales. Despite low grade metamorphism, primary textural zones resulting from initial cooling of the sheet are preserved and retain primary mineralogical and chemical variations which influenced the nature and extent of metamorphic recrystallization. This has resulted in a vertical sequence of secondary mineral assemblages through the intrusion. During early hydrothermal alteration K-feldspar replaced plagioclase micropheno-crysts in the marginal and contact zones, and olivine in the central zone was replaced by saponite. Subsequent regional metamorphism resulted in the development of (metastable) prehnite-pumpellyite-epidote assemblages in two sub-zones characterized by high Fe2O3. Elsewhere the assemblage prehnite-actinolite-epidote developed except in the contact and marginal zones where activity of CO2 suppressed both prehnite and pumpellyite. Both assemblages contain excess albite, quartz and chlorite and, on the basis of uniform mineral compositions over the area of an individual thin section, are considered to represent buffered equilibrium assemblages indicative of prehnite-pumpellyite and prehnite-actinolite facies conditions. A metamorphic temperature of 310° C at 1.85 kbar is obtained using the P-T-X grid of Liou, Maruyama & Cho (1985), which implies a field gradient of ∼ 44° C km-1. Assuming that metamorphism relates to burial, an overburden thickness of ∼ 7 km is indicated. Total maximum thicknesses, however, of Ordovician, Silurian and Lower Devonian strata, in the area, do not exceed 6 km indicating a field gradient of 52° C km-1. These relatively high gradients may possibly be related to concealed late Caledonian intrusions, or alternatively may result from high heat flow as a consequence of crustal thinning, rapid sedimentation and intense magmatic activity in a marginal basin setting.

Notes: Abstract Mineralogical and geochemical evidence indicates that partial melting and desulphidation have occurred in the Big Bell gold deposit. Decarbonation may also have occurred, to account for the lack of a carbonate alteration halo; this is compatible with the present data, but difficult to test.The Big Bell deposit consists of auriferous sulphide-bearing (‘lode’;) schists with muscovite and K-feldspar, and surrounding biotite schists, all derived by intense premetamorphic alteration of rocks of mafic composition. Assemblages which include cordierite-sillimanite-K-feldspar-garnet-biotite-quartz suggest peak metamorphic conditions of 4–5 kbar, and 650–700° C, based on phase relations, geobaro-meters and garnet-biotite Fe-Mg exchange partitioning. Partial melting occurred at peak metamorphism, particularly in the altered mafic rocks in and around the deposit, and its occurrence may have been essential to the preservation of the deposit. Melting greatly limited the importance of devolatilization reactions, resulting in negligible aqueous fluids and no means of removing appreciable gold. Minor gold loss may have accompanied desulphidation. A diversity of complex metamorphic assemblages occurs around the mine, compared to the assemblages developed regionally; variable bulk rock composition influences this contrast, but there is no evidence of higher metamorphic grades at the mine, nor that this might have been the prime control on the different assemblages in this narrow belt.It is suggested that the Big Bell and Hemlo deposits are the higher metamorphic grade equivalents of the more abundant greenschist facies gold deposits within Archaean greenstone belts. This interpretation is favoured by the host rock setting and geochemical characteristics of Big Bell. Alternative models that suggest that this class of deposit is a new type must account for the absence of high-grade equivalents of the greenschist facies deposits and also the lack of low-grade equivalents of the Big Bell/Hemlo type.Archaean gold deposits in high-grade metamorphic terrains have undergone a series of processes that are not recorded in the more typical gold deposits of the greenschist facies.

Notes: Abstract. Simple predictions of deep drainage in the soil profile are often required for preliminary planning of land management where the cost of direct measurement is not warranted. Soil hydraulic conductivity and drainage of water below the root zone can be related to the salt content at the bottom of the root zone, assuming steady-state balances of water and salt. A physically based empirical model uses readily measured soil properties to predict the quantity of drainage below the root zone under varying regimes of water management and shows a good relationship with ponded infiltration rate.

Notes: Abstract. The applicability of the ‘threshold concentration’ concept in formulating guidelines for irrigating with saline water was tested under field conditions on red-brown earths from different field experiments in south eastern Australia. Infiltration of water in the field and the effect of rainfall impact were studied using ring infiltrometers and a rotating-disc rainfall simulator.Three threshold concentration lines relating sodium adsorption ratio (SAR) and total cation concentration (TCC) in soil extracts were tested. These lines were based on laboratory tests of soil permeability, spontaneous dispersion and mechanical dispersion. They were found to predict the infiltration problems due to rainfall impact under three different surface soil conditions - bare soil without cultivation, soil with no tillage and complete crop cover, and cultivated soils without any crop cover.Infiltration rates in continuous pasture plots were predicted by the threshold concentration line of spontaneous dispersion except in a high salt treated soil where reduced plant growth affected evapotranspiration and water intake during infiltration.

Notes: Abstract. The soil solution may contain both plant nutrients and toxic ions. The total salt concentration affects both osmotic pressure and plant water stress. This review describes the main methods of evaluating soil salinity. They are listed as extraction methods (saturation and other soil extracts, suction cups), displacement methods (pressure membrane, centrifugation) and electrical methods of total salinity measurement (salinity sensors, four-electrode methods and time-domain reflectometry). The methods are compared so that the reader may choose the one most suitable for his purpose, based on cost, on the inherent advantages or drawbacks of the methods themselves, on his need for single or repeated measurements and either estimates of total salinity or the concentration of selected ions.

Notes: Abstract. The interaction of caesium isotopes with soil has been widely investigated and the influence of important soil properties studied. From the results of such work and a detailed knowledge of the physico-chemical properties of soils it is possible to classify Cumbrian soil according to its ability to immobilize caesium. The ‘immobilization capacity’ is a reflection principally of the clay mineral content and type, organic content, pH, ammonium content and potassium status. Although it is not quantifiable, the immobilization capacity permits ranking of the soils and indicates which areas may give rise to persistent caesium problems. Combination of the soil sensitivity classification with deposition data for Cumbria indicates that the mountainous region in the south-west of the county is the most vulnerable. This conclusion is supported by field evidence, since the area identified coincides closely with that where sheep movement and slaughter are restricted and where caesium remains persistently available to the plant-animal chain.

Notes: Abstract. Mathematical models describing radionuclide transport in soil developed for radiological assessment have increased in complexity over the last decade. In particular fairly simple ‘black box’, equilibrium approaches have given way to more complex, time-dependent, process-orientated methods. The increase in complexity of these models has outstripped the available data to specify, test and validate them. Current issues in model development include those that are associated with times up to a million years. Further development requires new laboratory and field research to provide adequate data to justify the inclusion or omission of known soil processes.

Notes: Abstract. The factors influencing uptake of radionuclides from soils into vegetation are discussed with reference to soil type, radionuclide, plant species and organ, and time since initial contamination. Gaps in knowledge are identified, particularly as highlighted by the unexpected behaviour of radiocaesium in many upland areas of Britain, following deposition after the Chernobyl accident. The importance of resuspended soil for contamination of aerial plant parts is also considered in relation to radionuclide type, vegetation height, and meteorological conditions.The development of an international database by the International Union of Radioecologists for soil to crop transfer factors of radionuclides derived from European experiments is briefly described. This database is now being used for statistical analyses aimed at quantifying the importance of environmental and biological factors in influencing uptake of radionuclides from the soil.

Notes: Abstract. The Rothamsted Drain Gauges built in 1870 comprise blocks of soil, 0.5, 1.0 and 1.5 m deep, isolated laterally by brickwork and undermined for the collection of drainage water but otherwise undisturbed. The soil has not been cropped, manured or cultivated subsequently. The annual nitrate leakages from these blocks were recorded for the 38 years from 1877/8 to 1914/5. The soil in the 0.5 and 1.5-m gauges lost on average about 45 kg ha−1 of nitrate-N per year during the first seven years of this period; that in the 1.0-m gauge lost slightly less. The overall decline in leakage was masked by large annual fluctuations attributable mainly to variation in rainfall. Fitting a simple function that assumed an exponential decline and took account of rainfall fluctuations gave a rate constant for each gauge from which the half-life could be estimated for the organic nitrogen feeding the leakage. The half-life for the 1.5-m gauge was 41 yr. The average nitrate leakage during the first seven years of the record differs little from estimates of the current leakage from soil carrying fully fertilized crops of winter wheat. This and the long half-life of the leakage show that pollution of drainage water by nitrate will not be controlled by limiting the use of fertilizer in the short term.

Notes: Abstract. Fifteen soil profiles in the Alltcailleach Forest in NE Scotland have been resampled after almost 40 years. The pH, in 0.01 M CaCl2, of the soil has decreased by 0.07 to 1.28 units in 80% of the surface organic horizons and by 0.16 to 0.54 units in 73% of the mineral horizons below 40 cm. The key factors governing increases and decreases in soil pH are changes in ground vegetation and tree canopy, although some effects of acid deposition cannot be ruled out.

Notes: Abstract. Erosion and excessive runoff from a crusting and hard-setting red-brown earth may he ameliorated with suitable management. A field trial, near Cowra, New South Wales, to assess the long-term effect of different tillage systems was used to compare the effect of direct drilling with conventional district cultivation practices under continuous wheat. The soil was sampled in the eighth year for assessment of the soil macropore structure, measurement of bulk density and hydraulic conductivity under tension. Vertical faces were prepared from resin impregnated blocks and the macropore structure described mathematically and visually using digital images and data generated from these images. Infiltration, bulk density and image analysis data all lead to the same conclusions about changes in pore structure. Under direct drilling no crust was evident, and there was greater macroporosity (〉 0.175 mm diameter in section). The treatment effects appeared to be significant to about 30 to 35 mm depth at the time of sampling. Greater root and faunal activity were observed under direct drilling.

Notes: Abstract. People in rural areas have disposed of their wastes on land for centuries, relying on the soil's ability to degrade and render harmless any toxic elements the wastes may contain. Leather tanneries produce a large amount of sludge and liquid wastes. The liquids contain much sodium and in most circumstances adversely affect both the soil and groundwater. However, in some countries they are used for irrigation. The sludges contain nitrogen, calcium, magnesium, phosphorus, trivalent chromium and some sodium. Given careful management these sludges can be used as soil amendments, either directly or after composting. Application rates of about 200 tonnes ha−1 have proved toxic to crops in pot trials, though larger applications than this have not adversely affected crops in the field. Much smaller rates, of less than 20 tonnes ha−1, have been used in the field to minimize nitrate contamination of groundwater. The effects of CrIII depend on complex interactions between the sludge, the soil to which it is added and the plant species grown. Safe limits for the disposal of the sludges and their long term effects are not known.

Notes: Abstract. Research on soil erosion under forest in Britain is reviewed. Erosion can increase as a result of afforestation in the uplands, sometimes with undesirable consequences for surface water quality. Published rates of erosion are usually close to ‘natural’ ones, at around 500 kg ha−1 yr−1. Of the forest operations that can lead to erosion ploughing is the most important. Clearfelling may also increase erosion, but little is known of its long-term effects. New developments in forestry may do much to reduce the risk of soil erosion, particularly the replacement of ploughing by subsoiling and the control of drain gradients.

Notes: Abstract. A simple model for droughtiness, when linked to the distribution of soil types in England and Wales, is potentially valuable for estimating drought-induced losses of yield in winter wheat at either particular places or in larger areas of the country. The model defines droughtiness, D, in terms of the soil water extractable by the crop, AP, and the adjusted potential moisture deficit, MD: D=AP−MD.The model should represent well the growth of actual crops of winter wheat if AP, which is based on laboratory measurements, accurately simulates the extraction of soil water by roots, if MD represents the cumulative transpiration of water by wheat crops by mid-July, and if the latter is an appropriate date for testing the effect of drought on grain growth. These three assumptions have been investigated using measurements of artificially draughted crops of winter wheat.The results indicate that mid-July is a good choice for a single date and MD a good representation of the water requirement of a wheat crop that has been draughted to the point where yield is beginning to be affected. For the deep-rooting crops studied, AP underestimates the soil water extracted by the crop, and therefore overestimates the susceptibility of the soil to droughtIf average MD values are replaced by means and a standard deviation the resulting normal distribution of D-values can be used to assess the probability that drought will limit yields. When applied to a droughtiness map of England and Wales with AP-values used for the soils the model predicts that the soils growing wheat will be susceptible to drought in 16–84 years out of every 100. Our results suggest, however, that this probably applies to shallow rooting or diseased crops and that for deep-rooting, healthy crops the drought risk is much less serious.

Notes: Abstract. Ploughing and tine cultivation to depths of 5, 15 and 25 cm were used to prepare land before sowing winter wheat. The methods were compared for two years with or without the presence of straw residues from the previous crop and the effects on crop growth and yield were assessed.Some combinations produced large differences in yield. In 1985–86 volunteer cereals were a problem where straw was not burnt, but deeper ploughing controlled them. In both years the concentrations and uptake of nitrogen, phosphorus and potassium differed during the early period of growth but not at maturity.

Notes: Abstract. The experiment conducted jointly by the Forestry Commission (FC) and the Institute of Terrestrial Ecology (ITE) at Beddgelert Forest, North Wales, studies the effect of conventional clear-felling (CF) and whole-tree harvesting (WTH) of Sitka spruce (Picea sitchensis) on the local ecosystem and on the future productivity of forestry at the site. Bulked soil samples were taken from Block 2 of the experiment just before felling in 1984 and two years later in 1986. Sub-samples from the horizons Ah, E and B+C were analysed for exchangeable and short- and long-term reserves of K using Ca-resin and strong acid extraction procedures. The flux of K through the soil profile after both CF and WTH resulted in a small increase in exchangeable K throughout the profile after both treatments, but in a loss of short-term reserves from the surface Ah horizons of both and an overall loss of these after WTH. The nutrient flux down-slope through the Ah horizon could result in differential nutrient deficiency in future. The data suggest that exchangeable and short-term reserves of K will support about two further cycles of conifers, with either CF or WTH, but that long-term reserves are likely to be released quickly enough to meet the needs of such a slow-growing crop; these would support about 30 cycles. Other nutrients, such as Ca or P, may prove to be more limiting than K.

Notes: Abstract. The effect of eight years of applications of five rates (0, 134, 269, 538 and 1075 m3 ha−1 a−1) of pig slurry on the soil strength two years later were studied in a field experiment. Soil strength in the 0–150 mm depth was measured on five occasions in winter using a hand-held recording cone penetrometer. On one occasion the penetration resistance at some depths greater than 100 mm was significantly (P 〈 0.001) decreased by adding more than 269 m3 of slurry ha−1 a−1. On three occasions different amounts of slurry caused significant differences in the rate of increase of penetration resistance with depth. Large applications of slurry may decrease penetration resistance because they increase organic matter, thereby increasing the water retention of the soil.

Notes: Abstract. In a field experiment over two years, broad beans (Vicia faba), cabbage (Brassica oleracea var. capitata), leeks (Allium porrum) and red beet (Beta vulgaris var. esculenta) were grown on a sandy clay loam soil in which a range of bulk densities and penetration resistances had been established by (1) thorough loosening to 0.9 m by trenching, (2) artificially compacting with tractor wheelings or (3) leaving unloosened.Loosening the soil substantially increased, and compacting it decreased, yields of all four crops. The mean penetration resistance of the subsoil at field capacity correlated negatively with dry matter production. The relationship was broadly similar for all crops and years, showing a decrease in dry matter production of about 1 t ha-1 per 0.5 MPa increase in resistance over the range examined.

Notes: Abstract. Soil water samples from five horizons in a stagnopodzol were collected regularly over a five-year period in a Sitka spruce plantation at Beddgelert Forest, North Wales. Samples were analysed for nitrate-N and ammonium-N. After felling, inorganic-N concentrations increased markedly in the C horizon, generally decreased in the surface horizons and showed little change in the E and Bs horizons. Fluxes through the C horizon increased after felling from 10 to 70 kg N ha-1 a-1, the latter being equivalent to leaching losses in intensive lowland agricultural systems. Trends in concentration and flux were attributed to seasonal temperature and rainfall variations.Nitrate-N dominated the dissolved inorganic-N, especially in the lower horizons. Nitrification was obviously active, despite the acid soil. Nitrate leaching losses occurred, even beneath the standing crop. On felling, cessation of nitrogen uptake allowed substantially more nitrate to be leached as no alternative sink was immediately available.

Notes: Abstract. Plots of spring barley, winter wheat and winter barley were sheltered with netting in a coastal area of North-East England. Crop height, yield components and grain yield were measured and compared with unsheltered crops receiving the same agronomic treatments. Plant height, ears m-2, 1000 grain weight and grain yield were increased by shelter in seasons with windy, dry weather during the tillering and stem extension phases. Negligible response was found in a wet season. One experiment suggests that greater applications of nitrogen fertilizer can counteract the effects of exposure.

Notes: Abstract . A combination of crop photography and crop radiance (near-infrared and red) values are used to measure ground cover establishment (leaf area). The high correlation between the green leaf area and near-infrared/red ratio (NIR/R) can provide a non-destructive method for monitoring crop growth. Using this relationship the effects of pre- and post-emergence herbicide treatments (Metazachlor) on winter oilseed rape were studied. It is shown that pre-emergence spraying reduces germination and autumn and spring growth rates. One management aspect of pre-emergence spraying is the delay in achieving ground cover. For marginal crops on soils with a high risk of erosion critical ground cover may be delayed by as much as two weeks.

Notes: Abstract. Simple equations are presented to predict the movement of a solute front in response to a given net cumulative flux of water. These equations do not require the assumption that the water content is at field capacity following water inflow, but may be simplified to that case. Equations are presented both for inert solutes and for anions and cations that react with soils through which they pass. The use of the equations is illustrated with simple hypothetical examples and, in one case, with field data.

Notes: Abstract. Cracks in dry saline montmorillonitic clay allow the soil to wet rapidly when flooded with negligible redistribution of salts. Once closed the only effective pathways remaining for the movement of leaching water are old root channels and faunal burrows. However, their effectiveness in conducting water and for leaching is severely restricted because of the lack of horizontal connections between them. Restructuring of clay can introduce sufficient permeability to a depth of about 0.7 to 0.8 m to allow salt to be leached provided that the soil's initial moisture content is sufficiently large to prevent disintegration upon wetting (about 24%) and that the clay is not allowed to become unsaturated during the leaching.

Notes: Abstract. The leaching of salt from a restructured saline heavy montmorillonite clay soil was studied experimentally. Restructured clay was placed in a 25-m long ± 0.75-m deep polythene-lined trench and leaching studied under flood irrigation. The hydraulic conductivity of the soil always exceeded 25 m d-1 during the experiment, and 85 % of the salt was leached within 16 days, most being leached in the first 8 days. The results indicate that restructuring clay could make the drainage of heavy clay soils economically feasible.

Notes: Abstract. Radionuclides in the environment come from a wide range of sources, some natural and some artificial. Their biogeochemical behaviour is influenced both by their own physico-chemical properties and by those of the soil with which they interact. The source of the radionuclides is important, as are any changes in physico-chemical characteristics that occur during transport and deposition. Once in the soil further reactions can occur as the radionuclides equilibrate. These include immobilization by ion exchange and precipitation and remobilization by complexing and dissolution. Land management and environmental changes that alter conditions in the soil result in changes in the interactions of radionuclides with the soil and thereby to changes in mobility in the soil.

Notes: Abstract. The soil of a restored opencast site in Northumberland is described and compared with the surrounding undisturbed soil. The most serious change in the soil was increased bulk density at the surface of the subsoil (0.3 m depth). Only a small overall reduction in pore space occurred, at the expense of continuous macropores. More of the pore space in the restored soil occurred as fine fissures as opposed to biotic channels, and the restored soil was weaker because of subsoiling and a loss of stability in microstructures. These changes are discussed in relation to root exploration, water movement and land management.

Notes: A new time-domain method is introduced for the calculation of theoretical seismograms which include frequency dependent effects like absorption. To incorporate these effects the reflection and transmission coefficients become convolutionary operators. The method is based on the communication theory approach and is applicable to non-normal incidence plane waves in flat layered elastic media. Wave propagation is simulated by tracking the wave amplitudes through a storage vector inside the computer memory representing a Goupillaud earth model discretized by equal vertical transit times. Arbitrary numbers of sources and receivers can be placed at arbitrary depth positions, while the computational effort is independent of that number. Therefore, the computation of a whole plane-wave vertical seismic profile is possible with no extra effort compared to the computation of the surface seismogram. The new method can be used as an aid to the interpretation of plane-wave decomposed reflection data where the whole synthetic vertical seismic profile readily gives the interpreter the correct depth position of reflection events.

Notes: A parameter estimation or inversion procedure is incomplete without an analysis of uncertainties in the results. In the fundamental approach of Bayesian parameter estimation, discussed in Part I of this paper, the a posteriori probability density function (pdf) is the solution to the inverse problem. It is the product of the a priori pdf, containing a priori information on the parameters, and the likelihood function, which represents the information from the data. The maximum of the a posteriori pdf is usually taken as a point estimate of the parameters. The shape of this pdf, however, gives the full picture of uncertainty in the parameters. Uncertainty analysis is strictly a problem of information reduction. This can be achieved in several stages. Standard deviations can be computed as overall uncertainty measures of the parameters, when the shape of the a posteriori pdf is not too far from Gaussian. Covariance and related matrices give more detailed information. An eigenvalue or principle component analysis allows the inspection of essential linear combinations of the parameters.The relative contributions of a priori information and data to the solution can be elegantly studied. Results in this paper are especially worked out for the non-linear Gaussian case. Comparisons with other approaches are given. The procedures are illustrated with a simple two-parameter inverse problem.

Notes: For converted waves stacking requires a true common reflection point gather which, in this case, is also a common conversion point (CCP) gather. We consider converted waves of the PS- and SP-type in a stack of horizontal layers.The coordinates of the conversion points for waves of PS- or SP-type, respectively, in a single homogeneous layer are calculated as a function of the offset, the reflector depth and the velocity ratio vp/vs. Knowledge of the conversion points enables us to gather the seismic traces in a common conversion point (CCP) record. Numerical tests show that the CCP coordinates in a multilayered medium can be approximated by the equations given for a single layer. In practical applications, an a priori estimate of vp/vs is required to obtain the CCP for a given reflector depth.A series expansion for the traveltime of converted waves as a function of the offset is presented. Numerical examples have been calculated for several truncations. For small offsets, a hyperbolic approximation can be used. For this, the rms velocity of converted waves is defined. A Dix-type formula, relating the product of the interval velocities of compressional and shear waves to the rms velocity of the converted waves, is presented.

Notes: An increase in the grain/pore size can cause the received echo signal to be distorted due to the preferential attenuation of the high-frequency components of the transmitted acoustic signal. Therefore, a rock acts as a low-pass filter and the property of the filter has a linear relation to the grain/pore size. The results have shown that P- and S-wave attenuations can be used to characterize selected internal characteristics of a rock. For example, an increase in the differential stress can cause an increase in the attenuation, due to the sensitivity of the S-wave to microcracks.

Notes: Three-component recordings permit the construction of particle trajectories. These three-dimensional pictures of particle motion show successive predominant directions of polarization and allow wave modes with distinct polarization directions to be recognized. A polarization selection, called ‘spatial directional filtering’, can be accomplished by several methods; four techniques are described. The application of these polarization filters to a noise shot, offset VSPs and CDP stack are also presented. This type of filtering is shown to cancel waves with undue polarization and to enhance the signal-to-noise ratio.

Notes: Gold, as a trace element, is not directly detectable by geophysical borehole logging techniques. Geophysical logging methods have therefore been used to infer the presence of structural features and alteration processes associated with gold. Since these features and processes differ with the style of gold mineralization, a unique relation between geophysical anomalies and gold cannot be established. However, in a particular environment, such an association can be established. We examined geophysical log data (self potential, induced polarization, resistivity, gamma, temperature, and temperature gradient) and drill core from the Barber-Larder property in NE Ontario. The geophysically detectable alteration processes associated with gold were identified as sericitization and pyritization, by means of factor analysis. A linear discriminant function was constructed which allowed the zones of economic mineralization (nominally 0.05 oz/ton or 1.5 ppm) to be identified with a 75% success rate based on geophysical log data alone. This rate was achieved without sophisticated non-linear data transformations.

Notes: A model of parallel slip interfaces simulates the behaviour of a fracture system composed of large, closely spaced, aligned joints. The model admits any fracture system anisotropy: triclinic (the most general), monoclinic, orthorhombic or transversely isotropic, and this is specified by the form of the 3 × 3 fracture system compliance matrix. The fracture system may be embedded in an anisotropic elastic background with no restrictions on the type of anisotropy. To compute the long wavelength equivalent moduli of the fractured medium requires at most the inversion of two 3 × 3 matrices. When the fractures are assumed on average to have rotational symmetry (transversely isotropic fracture system behaviour) and the background is assumed isotropic, the resulting equivalent medium is transversely isotropic and the effect of the additional compliance of the fracture system may be specified by two parameters (in addition to the two isotropic parameters of the isotropic background). Dilute systems of flat aligned microcracks in an isotropic background yield an equivalent medium of the same form as that of the isotropic medium with large joints, i.e. there are two additional parameters due to the presence of the microcracks which play roles in the stress-strain relations of the equivalent medium identical to those played by the parameters due to the presence of large joints. Thus, knowledge of the total of four parameters describing the anisotropy of such a fractured medium tells nothing of the size or concentration of the aligned fractures but does contain information as to the overall excess compliance due to the fracture system and its orientation. As the aligned microcracks, which were assumed to be ellipsoidal, with very small aspect ratio are allowed to become non-fiat, i.e. have a growing aspect ratio, the moduli of the equivalent medium begin to diverge from the standard form of the moduli for flat cracks. The divergence is faster for higher crack densities but only becomes significant for microcracks of aspect ratios approaching 0.3.

Notes: Ray theories are a class of methods often chosen to compute synthetic seismograms due to their efficiency and ability to deal with complex, three-dimensional inhomogeneous media. To deal with the large number of rays needed to compute synthetic seismograms, a ray generation algorithm is given which is capable of generating a numerical code describing each ray. The code describes a subset of all possible rays by considering only pre-critical reflections. In a horizontally plane-layered medium the generation of rays and computation of amplitudes and traveltimes can be efficiently accomplished by grouping the rays into reflection order and dynamic analogue groups. Expressions summing all unconverted rays and rays with a single mode conversion are given for source and receiver located at arbitrary positions within the medium. Examples of zero-offset synthetic VSPs obtained by this method are given.

Notes: The cross-correlation technique for increasing the anomaly-to-noise ratio is applied to the interpretation of resistivity profiles.To verify the method theoretically, resistivity profiles above a body having the shape of a parallelepiped were simulated with correlated and uncorrelated noise. Seven different electrode geometries are discussed. As a practical test, we considered profiles of a geoelectrical survey to locate tombs at a site of archaeological interest.When the shape and dimensions of the anomaly can be foreseen, the cross-correlation method can be applied and it gives a substantial improvement in the signal-to-noise ratio.

Notes: In the quantitative data interpretation for HLEM induction prospecting, a vertical half-plane model in an insulating medium is widely employed. For this assumption to be valid, the steeply dipping massive sulphide dykes must have large strike lengths and depth extents, but small thickness.We report investigations, using the laboratory scale-modelling method, on the response variation of large vertical conductors as the thickness is varied. We conclude that a steeply dipping large dyke can be approximated by a half-plane model only if its thickness is less than half the skin depth. An inductively thick conductor produces larger amplitudes and relatively higher quadrature compared to a thin conductor, even if both have the same induction number.

Notes: A model has been developed to relate the velocities of acoustic waves Vp and Vs in unconsolidated permafrost to the porosity and extent of freezing of the interstitial water. The permafrost is idealized as an assemblage of spherical quartz grains embedded in a matrix composed of spherical inclusions of water in ice. The wave-scattering theory of Kuster and Toksoz is used to determine the effective elastic moduli, and hence the acoustic velocities. The model predicts Vp and Vs to be decreasing functions of both the porosity and the water-to-ice ratio. The theory has been applied to laboratory measurements of Vp and Vs in 31 permafrost samples from the North American Arctic. Although no direct measurements were made of the extent of freezing in these samples, the data are consistent with the predictions of the model. Electrical resistivity measurements on the permafrost samples have demonstrated their essentially resistive behaviour. The ratio of resistivity of permafrost in its frozen state to that in its unfrozen state has been related to the extent of freezing in the samples.Electromagnetic and seismic reflection surveys can be used together in areas of permafrost: firstly an EM survey to determine the extent of freezing and then the acoustic velocity model to predict the velocities in the permafrost. The necessary transit time corrections can thus be made on seismic reflection records to compensate for the presence of permafrost.

Notes: The classical aim of non-linear inversion of seismograms is to obtain the earth model which, for null initial conditions and given sources, best predicts the observed seismograms. This problem is currently solved by an iterative method: each iteration involves the resolution of the wave equation with the actual sources in the current medium, the resolution of the wave equation, backwards in time, with the current residuals as sources; and the correlation, at each point of space, of the two wavefields thus obtained.Our view of inversion is more general: we want to obtain a whole set of earth model, initial conditions, source functions, and predicted seismograms, which are the closest to some a priori values, and which are related through the wave equation. It allows us to justify the previous method, but it also allows us to set the same inverse problem in a different way: what is now searched for is the best fit between calculated and a priori initial conditions, for given sources and observed surface displacements. This leads to a completely different iterative method, in which each iteration involves the downward extrapolation of given surface displacements and tractions, down to a given depth (the‘bottom’), the upward extrapolation of null displacements and tractions at the bottom, using as sources the initial time conditions of the previous field, and a correlation, at each point of the space, of the two wavefields thus obtained. Besides the theoretical interest of the result, it opens the way to alternative numerical methods of resolution of the inverse problem. If the non-linear inversion using forward-backward time propagations now works, this non-linear inversion using downward-upward extrapolations will give the same results but more economically, because of some tricks which may be used in depth extrapolation (calculation frequency by frequency, inversion of the top layers before the bottom layers, etc.).

Notes: Two particular sources of distortion, which may be encountered when applying tomographic imaging techniques to crosshole seismic data, have been investigated.Errors in survey locations of the shots and receivers can produce significant distortions in the images obtained. A simple method for solving simultaneously for the velocity field and shot and receiver location errors is presented and applied to synthetic and real data.Reflection and refraction of rays at velocity interfaces may produce poor density and angular coverage of the rays within the region of interest. It is shown that the effect of the velocity field on the ray coverage can significantly affect the resolution in the velocity image, even if ray bending is taken into account. One consequence of this effect is that, in some cases, little improvement in image quality is achieved by using curvi-ray rather than straight-ray inversion techniques, despite the occurrence of pronounced ray bending.

Notes: The WKB-method is used for the derivation of both the complex dispersion relation and displacement functions for Love channel-waves that propagate in a coal seam of varying thickness. The constant Q-model is used to describe the anelastic friction. With numerical solutions of the absorption-dispersion relation, the influence of thickness changes on the phase velocity and absorption coefficient of Love seam-waves is analysed at various frequencies. It is shown that the changes in the seam thickness can be optimally detected around the average Airy-phase frequency. An equivalence is pointed out between the wave guide structures: homogeneous with varying seam thickness and horizontally inhomogeneous with constant seam thickness.

Notes: Methods of minimum entropy deconvolution (MED) try to take advantage of the non-Gaussian distribution of primary reflectivities in the design of deconvolution operators. Of these, Wiggins’(1978) original method performs as well as any in practice. However, we present examples to show that it does not provide a reliable means of deconvolving seismic data: its operators are not stable and, instead of whitening the data, they often band-pass filter it severely. The method could be more appropriately called maximum kurtosis deconvolution since the varimax norm it employs is really an estimate of kurtosis. Its poor performance is explained in terms of the relation between the kurtosis of a noisy band-limited seismic trace and the kurtosis of the underlying reflectivity sequence, and between the estimation errors in a maximum kurtosis operator and the data and design parameters.The scheme put forward by Fourmann in 1984, whereby the data are corrected by the phase rotation that maximizes their kurtosis, is a more practical method. This preserves the main attraction of MED, its potential for phase control, and leaves trace whitening and noise control to proven conventional methods. The correction can be determined without actually applying a whole series of phase shifts to the data. The application of the method is illustrated by means of practical and synthetic examples, and summarized by rules derived from theory. In particular, the signal-dominated bandwidth must exceed a threshold for the method to work at all and estimation of the phase correction requires a considerable amount of data.Kurtosis can estimate phase better than other norms that are misleadingly declared to be more efficient by theory based on full-band, noise-free data.

Notes: The Offset Wenner resistivity sounding system provides for the extrapolation of the Wenner resistivity curve. The extrapolation technique was applied to data measured in the Solomon Islands and it is shown to be unreliable. An accurate method of predicting the reliability of extrapolation using measured resistances could not be found.

Notes: The concept of multifrequency induction logging simulations in the frequency range of 10 kHz to 1 MHz, applied to two-dimensional, axial symmetric model geometries, is presented. The scalar Helmholtz equation has been solved by a finite-element procedure. The model domain has been discretisized under the condition that the discontinuities in conductivity are represented by the nodes of the adjacent triangular elements. The modification of the signal distribution by the skin effect is illustrated for several models. Several sets of induction logs have been calculated with particular consideration of the frequency-dependent conductivities and permittivities.The improvement of a multifrequency inversion technique, based on standard least-squares methods, is shown for a two-layer model including borehole and invasion zones. Using this improved inversion technique we can state, as an additional inversion parameter, the frequency dependence of each inverted rock conductivity.

Notes: The wavefield in, and at the surface of, a homogeneous, isotropic, perfectly elastic half-space, excited by a traction distribution at the surface of the medium is investigated. The emitted wavefield is a spatial convolution of the surface tractions and the spatial impulse response. The properties of the wavefield in the far-field of the medium are derived and it is shown that the far-field particle velocity is essentially equal to a weighted sum of the time derivative of the integrated surface tractions, that is, of the components of the ‘ground force’. The theory is valid for an arbitrary geometry and orientation of the surface tractions, and is independent of the boundary conditions at the surface of the medium.The surface tractions are related to a source that consists of a mass distribution with an arbitrary force distribution imposed upon it. A boundary condition is introduced that accounts for the mass load and the forces applied to it but neglects vibrations within the mass. The boundary condition follows from the equation of motion of the surface mass load.The theory is applied to the Vibroseis configuration, using a P-wave vibrator model with a uniformly distributed force imposed on top of the baseplate, and assuming that horizontal surface traction components are absent. The distribution of displacement and stress directly underneath the baseplate of a single vibrator and an array of vibrators is investigated. Three different boundary conditions are used: (1) assuming uniform pressure, (2) assuming uniform displacement, (3) using the equation of motion of the baseplate as a boundary condition. The calculations of the distribution of stress and displacement over the plate for different elastic media and several frequencies of operation show that only the results obtained with the mixed boundary condition agree with measurements made in the field.The accuracy of three different phase-feedback signals is compared using synthetic data. Baseplate velocity phase-feedback leads to huge deviations in the determination of the far-field wavelet; reaction mass acceleration phase-feedback looks stable but neglects the differentiating earth filter; and phase-feedback to a weighted sum of baseplate and reaction mass accelerations becomes unstable with increasing frequency. The instability can be overcome using measurements over the whole baseplate.The model can be extended to a lossy layered earth.

Notes: This paper describes a method of generating pseudovelocity logs using measurements of electrical resistivity. A theoretical relation between electrical resistivity and transit time, which is applicable to a wide range of lithologies, has been developed. The application of this relation using a method which defines lithoresistivity zones as lithological intervals related to the same formation and showing small resistivity variations, has been tested in the Recôncavo sedimentary basin in Bahia, Brazil. A comparison of derived pseudovelocity logs with actual sonic logs for five wells whows the validity of the present approach.

Notes: Determination of petrography and pore fluid content is an ultimate goal of an integrated seismic-petrophysical study. For lack of a general inversion technique, forward modelling is useful in studying the relations between lithology, stratigraphy, pore fluid content and the seismic response. This report describes a study of two clastic sequences in Utah, from which 32 rock samples were analysed. A detailed petrographic study was done. Laboratory measurements were made of ultrasonic compressional- and shear-wave velocity as a function of pressure. We computed the velocities at seismic frequencies for the samples when dry, over-pressured, brine saturated, and oil saturated. The velocities were sensitive to the porosity, carbonate cementation and the depositional facies. We generated velocity profiles for hypothetical reservoirs for a range of saturation states. The velocity profiles were used to generate synthetic seismic shot gathers to study the seismic response of these clastic reservoirs. The fluid-saturation strongly affects the seismic respone, as does the presence of a coal seam. An amplitude change with offset is often observed. However, stratigraphy appears to have a stronger effect on the seismic response.

Notes: Two-dimensional VSP surveys are often conducted to provide structural illumination of the subsurface away from the borehole. The illumination is achieved through offsetting the source with respect to the downhole geophone. This inevitably gives rise to mode-conversions in both downgoing and upgoing wavefields.Migration of mixed-mode wavefields is complex because the velocity profile used for wavefield extrapolation is valid only for a particular propagation mode; the other mode always propagates at a different velocity. It is therefore advisable to separate the wave-types (P-wave and SV-wave) prior to migration. This may be achieved through wavemode filtering, a multichannel process which exploits the relation between propagation velocity, slowness of events at the recording array and particle motion. The necessary information about particle motion is available only if the VSP data are acquired with a three-component downhole geophone assembly.The wavemode filter partitions wave-types at the recording array; it provides no information about the various changes of propagation mode experienced by the energy as it travels from source to geophone. For the purpose of migration, the intermediate modes of propagation must be deduced.Much of the energy arriving at the receivers is P-wave which has followed the P-wave velocity profile from the source. It can therefore be imaged by conventional (Kirchhoff) migration. As an example of SV-wave imaging, a common mode-code is P-wave from source to reflector and SV-wave from reflector to geophone. Migration of such data calls for back-propagation of the geophone array wavefield, at SV-wave velocity, to the point in the subsurface where it is time-coincident with the forward propagated downwave, at P-wave velocity.

Notes: This paper gives a review of Bayesian parameter estimation. The Bayesian approach is fundamental and applicable to all kinds of inverse problems. Its basic formulation is probabilistic. Information from data is combined with a priori information on model parameters. The result is called the a posteriori probability density function and it is the solution to the inverse problem. In practice an estimate of the parameters is obtained by taking its maximum. Well-known estimation procedures like least-squares inversion or l1 norm inversion result, depending on the type of noise and a priori information given. Due to the a priori information the maximum will be unique and the estimation procedures will be stable except (in theory) for the most pathological problems which are very unlikely to occur in practice. The approach of Tarantola and Valette can be derived within classical probability theory.The Bayesian approach allows a full resolution and uncertainty analysis which is discussed in Part II of the paper.

Notes: A triple axis borehole magnetometer is described that consists of a Förster-probe (fluxgate) triplet (sensitivity 1 n T), a Förster-probe gradiometer (sensitivity 2 nT/40 cm), a gyro unit (mean angular drift approx. 0.5°/h) which is equipped with accelerometers (sensitivity 1/100°), and a data transmission unit (with multiplexer and 16-bit AD converter). The sensitive fluxgate-magnetometer can detect weakly magnetic or small source bodies. Data from the gyro and the accelerometers allow the 3-component magnetic field values to be transformed to north, east and vertical components. Since they do not rely on magnetically-determined directional data, the results are not disturbed by local anomalies of the magnetic declination. Furthermore, the magnetometer can also be used in vertical boreholes. 3-component measurements are carried out at discrete points in the neighbourhood of a source body to locate its position, and within the source body to determine the direction of magnetization. The strength of magnetization and information on magnetic classification are obtained by continuous measurement of one or more components within the source body. Calculation algorithms and computer programs are available to simplify data processing and interpretation. Survey examples are discussed.

Notes: Resistivity measurements were carried out in a survey area in the south of Germany. This area is characterized by complicated subsurface geology. Schlumberger full-arrays and their respective half-arrays were recorded simultaneously. The results obtained by the one-dimensional (1D) interpretation of the full-array measurements were incorrect because of a resistivity discontinuity. This discontinuity, under a relatively thick overburden, could only be located by the half-array soundings. Its exact location and the resistivity distribution in the subsurface were ascertained by comparing the sounding curves with 2D model curves, which are calculated by a finite-difference method.

Notes: A novel process is used to coat glass spheres with a polymerizing thermo-setting resin, only microns thick. Synthetic rocks of known grain size distribution and pore space characteristics are then made by heating the resin-coated glass spheres under compression in a special mould. The dynamic Young's moduli of these rocks are found to be affected by the percentage of resin content (cement) and the synthetic diagenesis rather than the grain size and permeability.