ALBERT

Notes: The equilibrium thermodynamics of the reaction:And the equilibrium constant is composed of activities formulated using ideal mixing on sites. Consideration is given to the evaluation of uncertainties in pressures calculated using the geobarometer. Preliminary testing suggests that the geobarometer has considerable potential. Much wider testing is now required.

Notes: Abstract. Pink piemontite-spessartine-bearing and grey-green spessartine-bearing manganiferous quartzose schists derived from siliceous pelagites, and green quartzofeldspathic schists, are described from the greenschist facies of the Haast Schist terrane, near Arrow Junction, western Otago. Electron microprobe data are reported for sphene, spessartine-rich garnet, manganoan epidote, piemontite, tourmaline, phengitic muscovite, chlorite, albite, haematite, rutile, manganoan calcite and chalcopyrite.Metamorphism occurred at about 6.4kbar, 400°C. Xco2 was above the quartz-rutile-calcite-sphene buffer (Xco2± 0.02) throughout the recorded metamorphic history of the piemontite schists. It dropped from above to below this critical buffering value in a spessartine-rich schist and it was close to or below the buffering value in the quartzofeldspathic schists. Production of piemontite required high fO2, believed to be inherited from MnOx in the parent pelagite. Substantial loss of O2 (e.g. minimum of 0.19% by weight in one rock) during diagenesis and/or metamorphism is inferred. In the grey-green schists this inhibited piemontite formation. Slight loss of O2 and Ca2+ accompanied minor late-stage replacement of piemontite by second generation spessartine. Observed zoning and mineral replacements indicate rise of temperature, drop in pressure, or invasion by solutions of lower fO2 and XCO2 equilibrated with surrounding schists.The detailed chemistry of the minerals studied correlates with available Mn and with bulk-rock (Fe3+ x 100)/(Fe2++ Fe3+). The oxidation ratio ranges from 24 in average green quartzofeldspathic schist, through 78 in average grey-green manganiferous quartzose schist, to almost 100 in some piemontite-bearing schists. As Fe2+ gives way to Fe3+, Mg/Fe ratios tend to rise in chlorite, phengite, tourmaline, spessartine, and calcite, Mn increases and Ti decreases in haematite, Mn increases in spessartine and calcite, and Fe increases in rutile. Available divalent cations are depleted relative to Al; chlorite is more aluminous, and phengite more paragonitic than in typical Haast schists.

Notes: Field, petrographic and microprobe investigations of metaclastic rocks, calcareous schists, marbles, chloritic calcareous meta-volcanic units and schists/paragneisses which crop out along the eastern portion of the Central East-West Cross Island Highway in Taiwan demonstrate that metamorphic intensity gradually increases eastward. The lower greenschist facies Slate Formation on the W contains completely recrystallized, pure albitic plagioclase, but at least some of the white micas (± chlorites) probably represent relict detrital flakes. Neo-blastic biotite and epidote occur sporadically in the Pihou(?) Formation, and increase dramatically eastward; concomitantly the abundance of carbonaceous matter decreases to zero in the eastern Tailuko zone, and the amount of chlorite + white mica diminishes somewhat. Epidote becomes more aluminous at higher metamorphic grade. Eastward, phengites change progressively to more muscovitic compositions as the proportion of biotite increases.A close approach to chemical equilibrium for the pre-Cenozoic, complexly deformed metamorphic basement assemblages is suggested by regular, systematic, major and minor element partitioning between analysed coexisting phases. Fractionation is less pronounced on the E, reflecting higher temperatures. Estimated physical conditions of recrystallization with αH2O and αCO2 moderate, are: T 〉 325 ± 75°C, P 〉 3 kbar (W); T 〉 425 ± 75°C, P 〉 4kbar(E).The gradual eastward increase in metamorphic intensity from the Slate Formation through the Pihou(?) Formation and the three Tailuko zones, as well as the relict precursor textures in the pre-Cenozoic layered basement rocks indicate that the observed paragenetic sequence could represent a synchronous Neogene recrystallization event, probably accompanying the Plio-Pleistocene collision of the Asiatic continental margin and the Luzon (Coastal Range) andesitic arc.

Notes: The structure, microstructure and petrology of a small area close to the village of Bard in Val d'Aosta (Italy) has been studied in detail. The area lies across the contact between the Gneiss Minuti (GM) and the Eclogitic Micaschist (EMS) Complexes of the Lower element of the Sesia portion of the Sesia-Lanzo Zone (Western Alps). Both complexes have undergone high-pressure metamorphism, but the metamorphic assemblages indicate a sudden increase in pressure in going across the contact from the GM to the EMS. Therefore, we interpret the contact as a thrust dividing the lower element of the Sesia into two sub-elements. This interpretation is supported by structural evidence.The early Alpine (90-70 Ma) metamorphic history is best preserved in the EMS and is one of increasing pressure associated with thrusting. The maximum P/T recorded in the EMS is 〉1500 MPa (〉15kbar) and 550°C and in the GM is 〈 1500-1300 MPa (〈 15-13 kbar) and 500-550°C. We suggest that the rocks were probably in an active Benioff zone during this time.From then on the histories of the GM and EMS are the same. Deformation continued and the thrust and thrust slices were folded during decreasing pressure. We interpret the first postthrusting deformation in terms of uplift associated with continued shortening of the crust and underplating after the Benioff zone had become inactive and a new Benioff zone had developed further to the north-west.A still later deformation and the Lepontine metamorphism (38 Ma) are related to continued uplift. Much of this deformation is characterized by structures indicative of vertical shortening and lateral spreading as the mountains rose above the general level of the surface.

Notes: In the Boi Massif of Western Timor the Mutis Complex, which is equivalent to the Lolotoi Complex of East Timor, is composed of two lithostratigraphical components: various basement schists and gneisses; and the dismembered remnants of an ophiolite. Cordierite-bearing pelitic schists and gneisses carry an early mineral assemblage of biotite + garnet + plagioclase + Al-silicate, but contain no prograde muscovite; sillimanite occurs in a textural mode which suggests that it replaced and pseudomorphed kyanite at an early stage and some specimens of pelitic schist contain tiny kyanite relics in plagioclase. Textural relations between, and mineral chemistries of, ferro-magnesian phases in these pelitic chists and gneisses suggest that two discontinuous reactions and additional continuous compositional changes have been overstepped, possibly with concomitant anatexis, as a result of decrease in Pload during high temperature metamorphism. The simplified reactions are: garnet and/or biotite + sillimanite + quartz + cordierite + hercynite + ilmenite + excess components. P-T conditions during the development of the early mineral assemblage in the pelitic gneisses are estimated to have been P + 10 kbar and T 〉 750°C, based upon the plagioclase-garnet-Al-silicate-quartz geobarometer and the garnet-biotite geothermometer. P-T conditions during the subsequent development of cordierite-bearing mineral assemblages in the pelitic gneisses are estimated to have been P + 5 kbar and T + 700°C with XH2O 〈 0.5, based upon the Fe content of cordierite occurring in the assemblage quartz + plagioclase + sillimanite + biotite + garnet + cordierite coexisting with melt.Final equilibration between some of the phases suggests that conditions dropped to P 〉 2.3 kbar and T 〉 600°C. A similar exhumation P-T path is suggested for the pelitic schists with early metamorphic conditions of P 〉 6.2 kbar and T 〉 745°C and subsequent development of cordierite under conditions in the range P = 3-4 kbar and T = 600-700°C. The tectonic implications of these P-T estimates are discussed and it is concluded that the P-T path followed by these rocks was caused by decompression during rifting and synmetamorphic ophiolite emplacement resulting from processes during the initiation and development of a convergent plate junction located in Southeast Asia during late Jurassic to Cretaceous time.

Notes: Plagioclase compositions vary from An0.1–2.5 to An32 with increasing grade in chlorite zone to oligoclase zone quartzofeldspathic schists, Franz Josef-Fox Glacier area, Southern Alps, New Zealand. This change is interrupted by the peristerite composition gap in rocks transitional between greenschist and amphibolite facies grade. Oligoclase (An20-24) and albite (An0.1–0.5) are found in biotite zone schists below the garnet isograd. With increasing grade, the plagioclase compositions outline the peristerite gap, which is asymmetric and narrows to compositions of An12 and An6 near the top of the garnet zone. In any one sample, oligoclase is the stable mineral in mica-rich layers above the garnet isograd, whereas albite and oligoclase exist in apparent textural equilibrium in adjacent quartz-plagioclase layers. The initial appearance of oligoclase in both layers results from the breakdown of epidote and possibly sphene. Carbonate is restricted to the quartz-plagioclase rich layers and probably accounts for the more sodic composition of oligoclase in these layers. The formation of more Ca-rich albite and more Na-rich oligoclase near the upper limit of the garnet zone coincides with the disappearance of carbonate and closure of the peristerite gap. Garnet appears to have only a localized effect on Ca-enrichment of plagioclase in mica-rich layers within the garnet zone. The Na-content of white mica increases sympathetically with increasing Ca-content of oligoclase and metamorphic grade.Comparison of the peristerite gap in the Franz Josef-Fox Glacier schists and schists of the same bulk composition in the Haast River area, 80 km to the S, indicates that oligoclase appears and epidote disappears at lower temperatures, and that the composition gap between coexisting albite and oligoclase is narrower in the Franz Josef-Fox Glacier area. It is suggested that a higher thermal gradient (38-40°C/km) and variations in Si/Al ordering during growth of the plagioclases between the two areas may account for these differences. In the Alpine schists the peristerite gap exists over a temperature and pressure interval of about 370-515°C and 5.5-7 kbar (550-700 MPa) PH2O.

Notes: In north-central Wopmay Orogen, syntectonic low-P(Buchan-type) suites of mineral isograds outline regional metamorphic temperature culminations that are associated, at the higher structural levels, with emplacement of early Proterozoic plutons in the west part of a deformed and eastward transported continental margin prism. The mapped isograds mark the first occurrence of biotite, staurolite, andalusite, sillimanite, sillimanite-K feldspar and K feldspar-plagioclase-quartz ± muscovite (granitic) pods in metapelites, with increasing proximity to the plutons.Microprobe analyses and field observations have resulted in the formulation of reactions for the ‘ideal’pelitic system K2O-Na2O-FeO-MgO-Al2O3-SiO2-H2O-Al2O3-SiO2-H2O, to account for the various mineral assemblages of each metamorphic zone. A P-T petrogenetic grid showing erosion surface P-T curves for the northern Wopmay Orogen pelites, compiled on the basis of the mapped isograds and the inferred reaction(s) for each metamorphic zone, documents a variation in exposed metamorphic pressure ranging between 2 and 4 kbar.The configuration of a new bathograd, based on the invariant model reaction sillimanite + K feldspar + plagioclase + biotite + quartz + vapor ± muscovite + liquid and interpolated across three metamorphic suites, is consistent with a major regional structure culmination and with independently determined pressures obtained from anorthite-grossular-quartz-Al2SiO5 geobarometry. The positive correlation between the configuration of the bathograd and the structural and pressure culmination points to the pressure-dependence of anatectic-granitic-pod mineral associations.

Notes: Variations in assemblage and composition of the constituent minerals in basic and intermediate metavolcanics encountered in the Zarouchla Group of the Phyllite-Quartzite Series are consistent with a progressive sequence, corresponding to temperature conditions estimated at 290-380°C (minimum values) under a total pressure greater than 3°5kbar and possibly as high as 5 kbar. In the absence of more critical evidence, the parageneses recorded in the metavolcanic rocks are interpreted as belonging to a prograde facies series from the lawsonite-albitechlorite facies through the pumpellyite-actinolite facies to the greenschist facies. The present distribution of mineral assemblages does not show a simple increase of metamorphic grade in a given direction but is apparently related to the tectonic evolution of the metamorphic sequence.

Notes: Abstract Mylonites from shear zones cutting Hercynian gneisses in the central Pyrenees have been studied in thin section and using the electron microprobe. The shear zones contain retrogressive greenschist facies assemblages implying introduction of an aqueous fluid during deformation in the zones. Textural evidence suggests that fluid-rock interaction occurred throughout the active life of the shear zones.Whole-rock chemical changes during deformation are documented in a variety of mylonitic lithologies and retrogressed country rocks. The overall effect was to reduce chemical differences between lithologies. Activity diagrams show that this would be expected if a hydrous fluid was circulating between different lithologies during deformation. In most cases fluid/rock ratios were relatively small resulting in gradual chemical changes and repeated recrystallization. ‘Open-system’behaviour with reduction in the number of phases is seen in some granite mylonites, suggesting focusing of fluid movement in parts of the shear zones. Continual fluid-rock interaction may have led to reaction-enhanced ductility in the shear zones over a long period of time. The source of fluid is uncertain, but may be related to underthrusting of material beneath the area investigated.

Notes: Abstract The hornblende-bearing basic gneisses in the Uvete area, central Kenya, were metamorphosed under a narrow range of P and T (6.5 ± 0.5kbar and 530 ± 40°C) of the staurolitekyanite zone in the Mozambique metamorphic belt. They show a wide variety of divariant and trivariant mineral assemblages consisting of hornblende, cumminatonite, gedrite, anthophyllite, chlorite, garnet, epidote, clinopyroxene, plagio-clase and quartz. The bulk and mineral chemistries and the graphical representation of phase relations show that each mineral assemblage approaches chemical equilibrium and defines a unique composition volume in the A′(Al + Fe3+− (13/7)Na)-F(Fe2+)-M′(Mg)-C′(Ca-(3/7)Na) tetrahedron. The composition volumes are distributed quite regularly and do not overlap each other.The phase relations in the Uvete area are in contrast with those in the staurolite-kyanite zone amphibolites in the Mt. Cube quadrangle, Vermont. The amphibolites there contain low-variance mineral assemblages formed under different values of μH2O and μCO2. These assemblages define overlapping composition volumes in the A′-F′-M′-C’tetrahedron.The mineral assemblages in the Uvete area are interpreted as having formed in equilibrium with fluid at a high and nearly constant μH2O value. Such a fluid composition was externally controlled by the supply of H2O-rich fluid expelled from the surrounding pelitic and psammitic rocks. The body size of the basic gneisses in the Uvete area (less than 400m in thickness) was small enough for the fluid to migrate completely.

Notes: Abstract Blueschist-facies rocks on the Seward Peninsula constitute a structurally coherent terrane measuring at least 100 × 150 km. Radiometric age data indicate that high-pressure metamorphism probably occurred in Jurassic rather than in Palaeozoic or Precambrian time, as previously suggested. Protolith sediments (Nome Group) are of intracontinental basin or continental margin type, and of lower Palaeozoic and possibly late Precambrian age, thus predating the high pressure metamorphism by more than 200 m.y.Blueschist-facies mineral assemblages were developed in almost all lithologies of the Nome Group, and are best preserved in FeTi-rich metabasites (glaucophane + almandine + epidote) and pelites (glaucophane + chloritoid + phengite). A lawsonite–crossite subfacies was developed in possible Nome Group rocks on the east flank of the Darby Mountains. Albite–epidote–amphibolite facies assemblages characterize Nome Group rocks in the southwestern part of the Peninsula. Metamorphism in the central zone of the terrane passed from early lawsonitic to subsequent epidote–almandine–glaucophane schist subfacies with the local development (east of the Nome River) of eclogitic assemblages.The high pressure metamorphic minerals were synkinematic with the development of mesoscopic-scale intrafolial isoclinal folds and a flattening foliation of consistent orientation. Initiation of uplift probably corresponded to the growth of barroisite rims on earlier sodic and actinolitic amphiboles, and partial post-kinematic greenschist facies replacements record later stages of decompression. Ophiolites and melange are not associated with the Seward Peninsula blueschists. The high-pressure metamorphism was caused by tectonic loading of a continental plate by an allochthon of indeterminate origin. The PT conditions of high pressure metamorphism were approximately 9–11 kbar, 400–450°C, thus falling between the PT paths of the Shuksan and Franciscan terranes.

Notes: Abstract Two periods of garnet growth (Gt1 and Gt2) have been found in the Finnmarkian nappes of north Norway. In the Kolvik Nappe (the lowest nappe) Gt1 has preserved an S2 syntectonic spiral inclusion fabric; in the Olderfjord Nappe an earlier S1 fabric and an interkinematic inter-D1–D2 fabric have been preserved in Gt1 whilst only the S1 fabric has been found in Gt1 in the Brennsvik Nappe (the highest nappe). In each nappe Gt2 overgrew a penetrative fabric (S2) wrapped around Gt1. In the Kolvik Nappe inclusion fabrics may be continuous from Gt1 into Gt2 but in the higher nappes there is a distinct break. Gt2 may have been partially syntectonic with D3 in the Brennsvik Nappe.Chemically Gt1 in the Kolvik Nappe and in parts of the Olderfjord and Brennsvik Nappes has antithetic Fe-Mn zoning. In all nappes XCa and XMg are weakly zoned in Gt1; XMg increases outwards and is greater in the higher nappes in Gt1 suggesting higher nucleation temperatures. In the Olderfjord and Brennsvik Nappes Gt2 is marked by increasing XCa, probably due to changing garnet-plagioclase equilibria, although the Fe/Mg ratio remains constant. XMg is higher in Gt2 than Gt1.Basement rocks within the nappe pile have an early pre-Finnmarkian growth (Gt1) and a later Finnmarkian growth (GtH) correlated with Gt2 on the basis of chemical zoning patterns.The diachroneity of Gt1 is ascribed to progressively earlier (compared to the structural development) cessation of overstepping of garnet-forming reactions before peak metamorphism in the higher nappes, resulting in earlier structural events being preserved.

Notes: Abstract Fluid inclusion studies of rocks from the late Archaean amphibolite-facies to granulite-facies transition zone of southern India provide support for the hypothesis that CO2,-rich H2O-poor fluids were a major factor in the origin of the high-grade terrain. Charnockites, closely associated leucogranites and quartzo-feldspathic veins contain vast numbers of large CO2-rich inclusions in planar arrays in quartz and feldspar, whereas amphibole-bearing gray gneisses of essentially the same compositions as adjacent charnockites in mixed-facies quarries contain no large fluid inclusions. Inclusions in the northernmost incipient charnockites, as at Kabbal, Karnataka, occasionally contain about 25 mol. % of immiscible H2O lining cavity walls, whereas inclusions from the charnockite massif terrane farther south do not have visibile H2OMicrothermometry of CO2 inclusions shows that miscible CH4 and N2 must be small, probably less than 10mol.%combined. Densities of CO2 increase steadily from north to south across the transitional terrane. Entrapment pressures calculated from the CO2 equation of state range from 5 kbar in the north to 7.5 kbar in the south at the mineralogically inferred average metamorphic temperature of 750°C, in quantitative agreement with mineralogic geobarometry. This agreement leads to the inference that the fluid inclusions were trapped at or near peak metamorphic conditions.Calculations on the stability of the charnockite assemblage biotite-orthopyroxene-K-feldspar-quartz show that an associated fluid phase must have less than 0.35 H2O activity at the inferred P and T conditions, which agrees with the petrographic observations. High TiO2 content of biotite stabilizes it to lower H2O activities, and the steady increase of biotite TiO2 southward in the area suggests progressive decrease of aH2O with increasing grade. Oxygen fugacities calculated from orthopyroxene-magnetite-quartz are considerably higher than the graphite CO2-O2 buffer, which explains the absence of graphite in the charnockites.The present study quantifies the nature of the vapours in the southern India granulite metamorphism. It remains to be determined whether CO2-flushing of the crust can, by itself, create large terranes of largeion lithophile-depleted granulites, or whether removal of H2O-bearing anatectic melts is essential.

Notes: Abstract There are discrete masses of un-deformed metabasite within the blueschist series of the island of Syros. Greece. Around the margins of these masses are zonal sequences through rocks showing intracrystalline deformation but without a geometric fabric, to rocks with discrete and anastomosing shear zones, and finally to penetratively foliated rocks with isolated relics of the original undeformed texture. Textural relics suggest that this spatial sequence is at least qualitatively also a temporal sequence.This progressive shear zone deformation took place concurrently with a glaucophane-epidote to eclogite reaction. The reaction pathways in the rocks that underwent the shear zone deformation can be compared with those in rocks of a similar composition that suffered a longer deformation history and show no relics of an undeformed parent. Although the final assemblages are in both cases the same, the pathways are different. These differences are in part related to reactions promoted by the change from local to bulk equilibrium on the onset of deformation in the rocks. They are also related to the crystallization and later breakdown during the sequence of progressive equilibration of a metastable phase, in this case an impure glaucophane.

Notes: Abstract The Rockley Volcanics from near Oberon, New South Wales occur within the aureole of the Carboniferous Bathurst Batholith and have been contact metamorphosed at P ∼ 100 ± 50MPa (10.5kbar) and a maximum T ∼ 565°C in the presence of a C–O–H fluid. Prior to contact metamorphism the volcanics were regionally metamorphosed and altered with the extensive development of actinolite, chlorite, plagioclase, quartz and calcite. The contact metamorphosed equivalents of these rocks have been subdivided into: Ca-poor (cordierite + gedrite), Mg-rich (amphibole + olivine + spinel), mafic (amphibole + plagioclase) and Ca-rich (amphibole + garnet + diopside; diopside + plagioclase; garnet + diopside + wollastonite) rocks.The chemistry of the minerals in the hornfelses was controlled by the bulk rock chemistry and fluid composition. Pargasites and hastingsites as well as an unusual phlogopite with blue green pleochroism, are found in Ca-rich hornfelses. A comparison of the assemblages with experimentally derived equilibria suggests that the fluid phase associated with the Ca-rich hornfelses was water-rich (Xco2= 0.1 to 0.3) while that associated with the Mg-rich hornfelses was enriched in CO2 (Xco2 〉 0.7). The different hornfels types have reacted to contact metamorphism independently in both their solid and fluid chemistries.

Notes: Abstract A garnet–hornblende Fe–Mg exchange geothermometer has been calibrated against the garnet–clinopyroxene geothermometer of Ellis & Green (1979) using data on coexisting garnet + hornblende + clinopyroxene in amphibolite and granulite facies metamorphic assemblages. Data for the Fe–Mg exchange reaction between garnet and hornblende have been fitted to the equation. In KD=Δ (XCa,g) where KD is the Fe–Mg distribution coefficient, using a robust regression approach, giving a thermometer of the form: with very satisfactory agreement between garnet–hornblende and garnet–clinopyroxene temperatures. The thermometer is applicable below about 850°C to rocks with Mn-poor garnet and common hornblende of widely varying chemistry metamorphosed at low aO2.Application of the garnet–hornblende geothermometer to Dalradian garnet amphibolites gives temperatures in good agreement with those predicted by pelite petrogenetic grids, ranging from 520°C for the lower garnet zone to 565–610°C for the staurolite to kyanite zones. These results suggest that systematic errors introduced by closure temperature problems in the application of the garnet–clinopyroxene geothermometer to the ‘calibration’data set are not serious. Application to ‘eclogitic’garnet amphibolites suggests that garnet and hornblende seldom attain Fe–Mg exchange equilibrium in these rocks.Quartzo-feldspathic and mafic schists of the Pelona Schist on Sierra Pelona, Southern California, were metamorphosed under high pressure greenschist, epidote–amphibolite and (oligoclase) amphibolite facies beneath the Vincent Thrust at pressures deduced to be 10±1 kbar using the phengite geobarometer, and 8–9kbar using the jadeite content of clinopyroxene in equilibrium with oligoclase and quartz. Application of the garnet–hornblende thermometer gives temperatures ranging from about 480°C at the garnet isograd through 570°C at the oligoclase isograd to a maximum of 620–650°C near the thrust. Inverted thermal gradients beneath the Vincent Thrust were in the range 170 to 250°C per km close to the thrust.

Notes: Sapphirine occurs in a 3-5 m wide zone between amphibole-lherzolite and garnetiferous metagabbro at Finero in the Ivrea Zone, NW Italian Alps. Layers consisting of plag + hb + sa + cpx + opx + sp + gt are interbanded with spinel pyroxenites, which may contain sapphirine replacing spinel. All minerals are very magnesian, with XMg between 0.78 and 0.92. Bulk rock analyses suggest that precursors to the sapphirine-bearing rocks were igneous cumulates of plagioclase + olivine + hornblende + spinel. Up to 16wt% CaO does not inhibit sapphirine formation and it is the unusually Mg-rich nature of the host rocks which allows sapphirine development. The early igneous assemblage was replaced by one of cpx + sa + hb +± plag at a pressure of 9 ± 1 kbar and temperatures of 900 ± 50°C. Subsequent rapid uplift caused the instability of gt, gt + hb, hb and sa + cpx to form opx + plag ± sp ± sa symplectites.

Notes: The Gran Paradiso basement complex of the French and Italian Alps is composed of metasediments, termed the gneiss minuti, and metabasic rocks, both of which are intruded by a late Hercynian granite. The Bonneval gneiss, which crops out at the western edge of the complex, is composed of highly deformed metasediments, volcanics and volcaniclastic rocks. Eclogites, now highly altered, occur in the metabasic rocks. Kyanite and blue-green amphibole are locally present in the gneiss minuti and aegirine plus riebeckite occur in the Bonneval gneiss. A moderately high pressure - low temperature metamorphic event of probable Alpine age occurred in the basement complex. This metamorphic event differs from that in the overlying Sesia unit and ophiolites of the Schistes lustrés nappe in being at lower pressures (below the ab = jd100+ qz transition) and post-dating the major (D2A) deformation. The origin of the metamorphism is discussed and interpreted as a probable consequence of the overlying nappe pile which was emplaced during the D2A event. Subsequent greenschist facies metamorphism in the basement complex is a consequence of thermal relaxation during uplift.

Notes: Abstract. A method for the quantitative analysis of the spatial relations of minerals is described. Dispersed distributions are formed by annealing and destroyed in post-tectonic migmatization. Aggregate distributions characterize solid-state differentiation, whereas leucosomes formed in systems of high fluid:rock ratio (in the examples studied, anatectic melts) show random distributions.Quantitative textural analysis can be used to indicate whether migmatization was post-tectonic or earlier, though caution is necessary if post-migmatite cooling is slow or if there is some minor deformation. More importantly, it can be used to discriminate melt-present from melt-absent leucosomes; this is exemplified by a suite of metamorphic and anatectic migmatites from the Scottish Caledonides.The textural evolution of anatexites with increasing melt percentage is traced. Initial feldspar porphyroblastesis occurs by Ostwald ripening via grain boundary melts; subsequently ophthalmites develop with fabrics and chemistry inherited from the palaeosome. At greater than 30% melt these inherited fabrics are wholly destroyed. Deformation prompts segregation into melanosome and leucosome; resultant leucosomes contain no inherited crystals. The scale of anatectic systems is fixed at the point at which segregation begins; ophthalmites provide evidence for melt and crystal transfer beyond original palaeosome boundaries. In contrast, metamorphic migmatites are necessarily small-scale systems because of diffusive constraints, and melanosomes are invariably produced.

Notes: The oligoclase-biotite zone of the Bessi area, central Shikoku is characterized by sodic plagioclase (XCa= 0.10–0.28)-bearing assemblages in pelitic schists, and represents the highest-grade zone of the Sanbagawa metamorphic terrain. Mineral assemblages in pelitic schists of this zone, all with quartz, sodic plagioclase, muscovite and clinozoisite (or zoisite), are garnet + biotite + chlorite + paragonite, garnet + biotite + hornblende + chlorite, and partial assemblages of these two types. Correlations between mineral compositions, mineral assemblages and mineral stability data assuming PH2O = Psolid suggests that metamorphic conditions of this zone are about 610 ± 25°C and 10 ± 1 kbar.Based upon a comparative study of mineralogy and chemistry of pelitic schists in the oligoclase-biotite zone of the Sanbagawa terrain with those in the New Caledonia omphacite zone as an example of a typical high-pressure type of metamorphic belt and with those in a generalized‘upper staurolite zone’as an example of a medium-pressure type of metamorphic belt, progressive assemblages within these three zones can be related by reactions such as:

Notes: An assemblage consisting of corundum, sapphirine, spinel, cordierite, garnet, biotite and bronzite is described from the Messina area of the Limpopo Mobile Belt, and consideration given to its petrogenesis. Various geothermometers and geobarometers have been applied in an attempt to determine the temperatures and pressures of metamorphism.A former coexistence of garnet and corundum is suggested to have developed during the earliest high pressure phase of the metamorphism, where temperatures exceeded 800°C and pressures as high as 10kbar may have been experienced. Subsequently, continuous retrograding reactions from medium pressure granulite facies at about 800°C and 8kbar towards amphibolite facies generated spinel, cordierite, sapphirine and possibly also bronzite. The most notable reaction was probably of the form: garnet + corundum = cordierite + sapphirine + spinel.

Notes: A review of currently available information relevant to the Basal Gneiss Complex (BGC) of Western South Norway, combined with the authors’own observations, leads to the following conclusions.1. Most of the BGC consists of Proterozoic crystalline rocks and probably subordinate Lower Palaeozoic cover.2. The last major deformation of these rocks was during the Caledonian orogeny and involved large-scale thrusting, recumbent folding and doming. The structural development of the BGC is closely tied in with that of the Caledonian allochthon.3. The whole eclogite-bearing part of the BGC has suffered a high pressure metamorphism with conditions of between 550°C, 12.5 kbar (Sunnfjord) and about 750°C, 20 kbar (Møre og Romsdal) at the metamorphic climax.4. This metamorphism was of Caledonian age, probably rather early in the Caledonian tectonic history of the BGC and is considered to have been a rather transient event.By setting these conclusions in a framework provided by geophysical evidence for the deep structure of the crust in southern Norway we have constructed a geotectonic model to explain the recorded metamorphic history of the BGC. It is suggested that considerable crustal thickening was caused by imbrication of the Baltic plate margin during continental collision with the Greenland plate. This resulted in high pressure metamorphism in the resulting nappe stack. Progradation of the suture caused underthrusting of the Baltic foreland below the eclogite-bearing terrain causing it to emerge at the Earth's surface, aided by tectonic stripping and erosion.Application of isostacy equations to the model shows that eclogites can be formed by in-situ metamorphism in crustal rocks and reappear at the land surface above a normal thickness of crust in a single orogenic episode of approximately 65-70 Ma duration.

Notes: Detailed geochronological, structural and petrological studies reveal that the geological evolution of the Field Islands area, East Antarctica, was substantially similar to that of the adjacent Archaean Napier Complex, though with notable differences in late and post Archaean times. These differences reflect the area's proximity to the Proterozoic Rayner Complex and consequent vulnerability to tectonic process involved in the formation of the latter. Distinctive structural features of the Field Islands are (1) consistent development of a discordant, pervasive S3 axial-plane foliation; (2) re-orientation of S3 axial planes to approximate to the subsequent E-W tectonic trend of the nearby Rayner Complex; (3) selective retrogression by a post-D3 static thermal overprint; and (4) relatively common development of retrogressive, E-W-trending, mylonitic shear zones.Peak metamorphic conditions in excess of 800°C at 900 ± 100 M Pa (9 kbar) were attained at one locality following, but probably close to the time of D2 folding. D3 took place in late Archaean times when metamorphic temperatures were about 650°C and pressures were about 600 MPa (6 kbar). Later, temperatures of 600 ± 50°C and pressures of 700 MPa (7kbar) were attained in an amphibolite-facies event, presumably associated with the widespread granulite to amphibolite-facies metamorphism and intense deformation involved in the formation of the Rayner Complex at about 1100 Ma. The area was subsequently subjected to near-isothermal uplift.Rb-Sr isotopic data indicate that the pervasive D3 fabric developed at about 2400–2500 Ma, and this age can be further refined to 2456+8-5 Ma by concordant zircon analyses from a syn-D3 pegmatite. All zircons were affected by only minor (〈7–10%) Pb loss and/or new zircon growth during the Rayner event at about 1100Ma. Thus the 450–850 μg/gU concentrations of these zircons were too low to cause sufficient lattice damage over the 1350 Ma (from 2450 Ma) for excessive Pb to be lost during the 1100 Ma event. The emplacement of pegmatite at 522 ± 10 Ma substantially changed the Rb-Sr systematics of the only analysed rock that developed a penetrative fabric during the 1100 Ma event. Monazite in this pegmatite contains an inherited Pb component, which probably resides in small opaque inclusions.A good correlation is found between Rb-Sr total-rock ages and rock fabric. U-Pb zircon intercepts with concordia also mostly correspond to known events. However, in one example a near perfect alignment of zircon analyses, probably developed by mixing of unrelated components, produced concordia intercepts that appear to have no direct geochronological significance.

Notes: Coexisting Ca-poor and Ca-rich pyroxenes in granulites at Cape Riche, in the Precambrian Albany-Fraser Province, Western Australia, are dominantly chemically homogeneous within individual samples, suggesting a major episode of equilibration. However, occasional grains in a few samples contain exsolved domains interpreted as relics of an earlier, higher-T assemblage. Pyroxene pairs in ten, presumably isothermal, samples from a restricted area are used to (i) assess the suitability of several versions of the two-pyroxene thermometer for application to metamorphic rocks, and (ii) determine the thermal history of the Cape Riche pyroxenes.The various versions of the two-pyroxene thermometer applied to the well-equilibrated homogeneous pyroxene grains show poor to good precision and yield mean temperatures varying widely from 683° to 893°C, in the following order of increasing T: Lindsley (1983; opx version), 683°± 11°C; Kretz (1982; KD version), 705°± 19°C; Ross & Huebner (1975), 709°± 30°C; Kretz (1982; solvus version), 735°± 24°C; Fonarev & Graphchikov (1982; opx version), 〈750°C; Lindsley (1983; cpx version), 784°± 40°C; Fonarev & Graphchikov (1982; cpx version), ~820°± 30°C; Wood & Banno (1973), 849°± 16°C; Powell (1978), 854°± 23°C; Wells (1977), 893°± 10°C. Independent T estimates, based on mafic assemblages and garnet-biotite thermometry, suggest that the major episode of metamorphism occurred at 700-800°C (P ~ 5 kbar). Therefore the Wells, Powell, Wood & Banno and Fonarev & Graphchikov (cpx) temperatures are almost certainly too high. In the absence of a more precise independent T estimate it is difficult to assess the relative merits of the results obtained from the remaining versions of the two-pyroxene thermometer, none of which can be unequivocally demonstrated to be seriously in error, though the Lindsley (opx) T is probably too low. Other significant shortcomings evident in the results include the relatively poor precision obtained from the three methods based on purely graphical representation of the augite limb of the solvus (i.e., the Ross & Huebner, Fonarev & Graphchikov (cpx) and Lindsley (cpx) versions), and the apparent dependence of derived T on Mg/Fe2+ ratio for the Powell, Wood & Banno and Lindsley (cpx) methods.For the bulk compositions of exsolved domains, the different versions of the two-pyroxene thermometer yield mean temperatures 23° to 82°C (overall mean, 65°C) higher than for homogeneous grains in the same samples. These exsolved domains are interpreted as relics of a higher-T (peak?) metamorphic assemblage, rather than an igneous precursor.

Notes: This work presents the results of a fluid inclusion study of an amphibolite-granulite facies transition in West Uusimaa, S.W. Finland. Early fluid-inclusions in the granulite facies area are characteristically carbonic (CO2), in contrast to predominantly aqueous early inclusions in the amphibolite facies area. These early inclusions can be related to peak metamorphic conditions (750-820°C and 3-5 kbar for peak granulite facies metamorphism). Relatively young CO2 inclusions with low densities (〈0.8g/cm3) indicate that the first part of the cooling history of the rocks was characterized by a near isothermal uplift.N2-CH4 inclusions, with compositions ranging between pure CH4 and pure N2 (Raman spectral analysis), were found in the whole area. They are probably syn- or even pre-early inclusions. Only nearly critical homogenizing inclusions have been found (low density). Pressure estimates, based on densities of early fluid inclusions, show that the rapid transition of amphibolite towards granulite facies metamorphism is virtually isobaric. Granulite facies metamorphism in West Uusimaa is a thermal event, probably induced by the influx of hot, CO2-bearing fluids.

Notes: Multisystems of n+k (k 〉 3) phases are very complicated and knowledge of them has suffered as a result. The successful solution of the topological relationships in n+ 3 phase multisystems by Zen (1966, 1967) and Zen & Roseboom (1972) has aroused much interest regarding what will happen in a multisystem of more than n+ 3 phases. Since 1979, some important research results on this topic have been published. These results have expounded the substantial rules governing the appearance of phase relations in phase diagrams of n - k (k 〉 3) phase multisystems. The most significant conclusions include: (1) It is impossible to incorporate all the possible phase relations in an n+k (k 〉 3) phase multisystem in a single closed net. Therefore, it is no longer enough to use only a single closed net to depict the topological relations involved in these types of multisystems. Instead, one or more groups of closed nets, namely the complete system(s) of closed nets are necessary for this purpose. (2) A principle called the Combination Principle has been proposed and proved. It states: Any closed net of one n+k (k 〉 3) phase multisystem must be a combination of two or more distinct n+ 3 order submultisystem closed nets belonging to the given n+k phase multisystem, if it is not one of the n+ 3 order submultisystem closed nets itself. The combination principle provides both a theoretical basis and a practical method for the construction of closed nets and, hence, for the derivation of the real phase diagrams for any n+k (k 〉 3) phase multisystem. (3) A theorem on divariant-assemblage-characteristic-stability-polygons is also important to our understanding of the n+k (k± 3) phase multisystem closed nets. This theorem can be stated as follows: A divariant assemblage of an n+k (k± 3) phase multisystem will be stable in an l-polygon lacking diagonals in an appropriate set of closed-net-diagrams, and this l-polygon may be at least a triangle, and at most a k-polygon. In addition, the closed-net-diagrams of unary and binary n+ 4 phase multisystems derived respectively by Guo (1980b, 1980c, 1981a) and by Roseboom & Zen (1982) have also been summarized. The combination principle is applied to a practical petrological problem in this paper, dealing with 7 phases in the system FeO-Fe2O3-SiO2.

Notes: Abstract The stability of quartz-chloritoid-staurolite-almandine-cordierite and aluminium silicates is used to constrain both metamorphic conditions and pressure-temperature trajectories for two localities within the 2700 Ma Archaean Yilgarn Block in Western Australia. Available experimental data are used to calculate thermodynamic data for a self-consistent set of equilibria between these minerals. A lower amphibolite facies locality from the margin of a lower strain area contains assemblages including quartz-chloritoid-staurolite-garnet-biotite with altered cordierite replacing chloritoid, quartz-staurolite-andalusite, and quartz-cordierite-andalusite-biotite. This locality was heated to 530–560°C in the andalusite field, at 4.2 kbar. A sample from a mid- to upper-amphibolite facies, highly strained locality contains relict staurolite enclosed by andalusite, in turn replaced by cordierite and muscovite with biotite and sillimanite in the matrix. The assemblage was heated isobarically from conditions near the maximum experienced by the lower grade locality of 560°C at 4.2 kbar to temperatures in excess of the andalusite-sillimanite transition but within the quartz plus muscovite stability field (600–650°C). The higher grade locality is close to a granitoid dome and sections based on gravity profiles reveal that this locality is underlain by granitoid at shallow depths. The higher grade metamorphism apparently reflects superposition of the thermal aureole on regional metamorphic conditions similar to those in the lower grade areas.

Notes: Abstract A major system of steep Caledonian shear zones, of regional extent, has been identified in NE Scotland. The shear zones affect a wide range of lithologies, including Argyll and Southern Highland Group Dalradian, ‘Younger Basic’intrusives and their hornfelses, and also the earlier of the more acid intrusions. The observed fabrics and parageneses are consistent with low-pressure amphibolite facies metamorphism. These shear zones represent a phase of movement which occurred in the 490-465 Ma interval when ambient temperatures were still high, and it is concluded that this is the principal control on the metamorphic grade achieved within the shear zones, although local anomalies may exist.

Notes: Abstract An experimental study of the system CaCO3–MgCO3–FeCO3 was undertaken in order to calibrate the iron correction to the calcite–dolomite geothermometer, which is based on the solubility of magnesium in calcite in the assemblage calcite + dolomite. The experiments, at 450°C and lower temperatures, resulted in products with a very small grain size and incomplete equilibration. However, application of a carefully-devised automatic data processing algorithm to analyses of the phases in experimental charges, combined with a thermodynamic analysis, results in geothermometer diagrams which should be preferred to previous theoretical predictions.

Notes: Abstract Biotite, garnet, staurolite and kyanite isograds in pelitic metasedimentary rocks are developed as a result of thermal metamorphism around syntectonic granitoids in Eastern Rouergue (France). Temperature estimates range between 400°C and 650°C at about 6.5 kbar. Geothermobarometry shows a steep isobaric T gradient which is consistent with the interpretation that the metamorphic highs are thermal aureoles. High grade rocks show evidence of two staurolite forming reactions in the presence of plagioclase and the absence of chlorite that have not been described previously in the literature. The reaction that occurs in the middle staurolite zone, alm-rich ga + Ca-rich pla + Na-rich mu gro-rich ga + Na-rich pla + st + Na-poor mu, is considered to be prograde, whereas the reaction that occurs in the kyanite zone, alm-rich ga + Ca-rich pla + w st + Ca-rich ga + Na-rich pla + qz, is retrograde. The topology of these reactions is illustrated in terms of end member compositions for the systems KNaFASH and KCaFASH, respectively.

Notes: An occurrence of quartz-eclogite is described from the Inner Schieferhülle unit of the Pennine Basement Complex in the SE Tauern Window, Austria.Field relations strongly suggest a pre-Alpine age for the primary eclogitic mineral assemblage (garnet + omphacite + quartz + rutile). This implies that there was no connection between the formation of these eclogites and the late Cretaceous and Tertiary tectonic evolution of the Eastern Alps. The quartz-eclogite mineral assemblage crystallized under conditions of 620 ± 100°C and at pressures in excess of 12 kbar, and suffered amphibolitic overprinting of Alpine and possibly Hercynian age.A four-stage polymetamorphic history is proposed for the Inner Schieferhülle:

Notes: Mafic and ultramafic xenoliths in a basaltic cone at The Anakies in south-eastern Australia are geochemically equivalent to continental basaltic magmas and cumulates. The xenolith microstructures range from recognizably meta-igneous for intrusive rocks to granoblastic for garnet pyroxenites. Contact relationships between different rock types within some xenoliths suggest a complex petrogenesis of multiple intrusive, metamorphic and metasomatic events at the crust/mantle boundary during the evolution of south-eastern Australia. Unaltered spinel lher-zolite, typical of the uppermost eastern Australian mantle, is interleaved with or veined by the metamorphosed intrusive rocks of basaltic composition.Geothermobarometry calculations by a variety of methods show a concordance of equilibration temperatures ranging from 880°C to 980°C and pressures of 12 to 18 kbar (1200-1800 mPa). These physical conditions span the gabbro to granulite to eclogite transition boundaries. The water-vapour pressure during equilibration is estimated to be about 0.5% of the load pressure, using amphibole breakdown data. Large fluid inclusions of pure CO2 are abundant in the mineral phases in the xenoliths, and it is suggested that flux of CO2 from the mantle has been an important heat source and fluid medium during metamorphism of the mafic and ultramafic protoliths at the lower crust/upper mantle boundary.The calculated pressures and temperatures suggest that the south-eastern Australian crust has sustained a high geothermal gradient. In addition, the nature of the mineral assemblages and the contact relationships of granulitic rock with spinel lherzolite, characteristic of mantle material, suggest that the Moho is not a discrete feature in this region, but is represented by a transition zone approximately 20 km thick. These inferences are in agreement with geophysical data (including seismic, heat-flow and electrical resistivity data) determined for south-eastern Australia.Underplating at the crust/mantle boundary by continental basaltic magmas may be an important alternative or additional mechanism to the conventional andesite model for crustal accretion.

Notes: The occurrence of lawsonite is described from pelitic schists of the lower-grade part of the pumpellyite-bearing subzone of the chlorite zone in the Asemi River area of central Shikoku. The lawsonite-bearing parageneses are consistent with the generally accepted view that the Sanbagawa facies series represents higher pressures than the lawsonite-bearing facies series in New Zealand.

Notes: Abstract Late Archaean orthogneisses and aluminous and iron-rich metasedimentary rocks intruded by anorthosite and a ferrodiorite-granite suite were completely recrystallized during Proterozoic granulite facies metamorphism. Geobarometry and geothermometry indicate P-T conditions of around 7.5kbar. 700°C, with a CO2-rich fluid phase and logfO2 at or below -16. A two-stage high-grade history of near isochemical corona growth is preserved in metasediments with the reaction cycle opx + plag + H2O → hbl+gar+SiO2→ opx+plag+H2O. End product compositions resemble those of the initial phases, and the only mobile components were SiO2 and/or H2O. The coronas reflect shortlived fluctuations in chemical activity at essentially constant P and T, contrary to simple progressive change in equilibrium parameters recorded by most corona-bearing textures.

Notes: In a polymetamorphic, felsic, biotite-bearing gneiss, biotite has reacted to form magnetite and microcline. The resulting structure is a magnetite core surrounded by a mantle of feldspar and quartz normally not exceeding 20mm in diameter. Measurements of oxygen isotope ratios disclose disequilibrium between mantle microcline and mantle quartz and also between mantle and matrix minerals of the same species. A clustering of temperature estimates from the oxygen isotope distribution between magnetite and quartz and between magnetite and microcline in the interval 550 to 600°C suggests an approach to oxygen isotope equilibrium. No signs of a re-equilibriation of the reacting biotite can be found.

Notes: Abstract Clinopyroxenes and garnets from 11 blueschist-facies Fe-rich eclogite samples from the Voltri Group show a wide range of chemical compositions. Detailed analyses of single pyroxene and garnet grains show wide and scattered chemical inhomogeneity, the KD(KD= (Fe2+/Mg)Gt/(Fe2+/Mg)Cpx) ranges from 20 to 87 based on rim analyses only. The data obtained indicate that the mineral pairs never attained equilibrium under uniform P-T conditions and that the compositions of the metamorphic minerals were influenced mainly by the composition of the pre-metamorphic minerals and by topotactical reactions.

Notes: Abstract At Sulitjelma, Norway, there is a major inversion of metamorphic isograds beneath an inverted but undisrupted ophiolite. The flysch-like Furulund schist in which the inverted isograds occur is also inverted and the early folds in it are downward facing. The isograds cut across the axial surfaces of early folds and across the schistosity. These relationships are explained as the consequence of metamorphism during the progressive development of a large overfold. The inverted limb of the overfold is regarded as a major, thick, gently-dipping shear zone, separating the lower-grade, lower part of the Caledonian allochthon below from the higher-grade upper part of the allochthon above. The association between stratigraphical inversion, downward-facing of syn-schistosity folds and metamorphic inversion is explained by the progressive development of the shear zone. It is suggested that the presence of such shear zones is a common feature of orogenic belts formed by continental collision.

Notes: Abstract The Hercynian granitic basement which forms the Tenda Massif in NE Corsica represents part of the leading edge of the European Plate during middle-to-late Cretaceous (Eoalpine) high P metamorphism. The metamorphism of this basement, induced by the overthrusting of a blueschist facies (schistes lustrés) nappe, was confined to a major ductile shear zone (c. 1000m thick) within which deformation increases upwards towards the overlying nappe. Metamorphism within the basement mostly records lower blueschist facies conditions (crossite + epidote) except near the base of the shear zone where the greenschist facies assemblage albite + actinolitic amphibole has developed instead of crossite. Study of the primary mafic phase breakdown reactions within hornblende granodiorite reveals the following metamorphic zonation. Zone 1: biotite to chlorite. Towards zone 2: biotite to phengite. Zone 2: Hornblende to actinolitic Ca-amphibole + albite + sphene, and biotite to actinolitic Ca-amphibole + albite + phengite + Ti-ore + epidote. Zone 3: Hornblende to crossite + low Ti-biotite + phengite + sphene, and biotite to crossite + low Ti-biotite + phengite + Ti-ore + sphene ± epidote. P-T conditions at the base of the shear zone are estimated to have been 390-490°C at 600-900 M Pa (6-9kbar) and the Corsican basement is therefore deduced to have been buried to 20-30 km during metamorphism. This relatively shallow metamorphism contrasts with some other areas in the Western Alps where the Eoalpine event apparently buried the European continental crust to depths of 80 km or more. As there is no evidence for a long history of blueschist facies metamorphism prior to the involvement of the European continent, it is deduced that the Eoalpine blueschists were produced during the collision of the Insubric plate with Europe, rather than during Tethyan intraoceanic subduction. Coherent blueschist terrains such as the schistes lustres probably record buovant feature collision and obduction tectonics rather than any preceding oceanic subduction.

Notes: Abstract Three major blastomylonitic synmetamorphic (epidote amphibolite to mid amphibolite facies) shear zones are seen on the NW coast of the Mullet Peninsula in NW Mayo. These shear zones occur at the contacts of major structural units and in an imbricated slice where rocks of the Erris Complex are deformed and chemically modified. Chemical changes associated with individual shear zones have been deduced by comparing the compositions of various gneisses both within and adjacent to the shear zones. Compositional changes are different in the constituent rock-types within each unit and many elements normally considered immobile have been selectively mobilized within the shear zones. Little evidence of wholesale metasomatic introduction of components into these shear zones was found to accompany the selective mobilization.

Notes: Abstract. In the Kamuikotan zone, central Hokkaido, Japan, two distinct types of metamorphic rocks are tectonically mixed up, along with a great quantity of ultramafic rocks; one type consists of high-pressure metamorphic rocks, and the other of low-pressure ones. The high-pressure metamorphic rocks are divided into two categories. (1) Prograde greenschist to glaucophaneschist facies rocks derived from mudstone, sandstone, limestone, a variety of basic rocks such as pillow and massive lavas, hyaloclastite and tuff, and radiolarian (Valanginian to Hauterivian) chert, among which the basic rocks and the chert, and occasionally the sandstone, occur as incoherent blocks (or inclusions) enveloped by mudstone. (2) Retrograde amphibolites with minor metachert and glaucophane-calcite rock, which are tectonic (or exotic) blocks enclosed within prograde mudstone or serpentinite, or separated from these prograde rocks by faults. The K-Ar ages of the prograde metamorphic rocks (72, 107 and 116 Ma on phengitic muscovites) are younger than those of the retrograde rocks (109, 132, 135 and 145 Ma on muscovites, and 120 Ma on hornblende). The low-pressure metamorphic rocks consist of the mafic members of an ophiolite sequence with a capping of radiolarian (Tithonian) chert with the metamorphic grade ranging from the zeolite facies, through the greenschist (partly, actinolite-calcic plagioclase) facies to the amphibolite (partly, hornblende-granulite) facies. The low-pressure metamorphism has a number of similarities with that described for‘ocean-floor’metamorphism. The tectonic evolution of such a mixed-up zone is discussed in relation to Mesozoic plate motion.

Notes: The rocks of the Scourian Complex have been intensively studied, but there is still no consensus as to the conditions of the granulite-facies metamorphism preserved in these rocks. Recent estimates of these conditions fall into two groups, one at 820-920°C and ca. 11 kbar and the second at ca. 1000°C and 〉12 kbar. Investigation of a variety of rocks shows that the recorded conditions vary with grain-size, with higher-grade conditions recorded by the cores of coarser (ca. 10 mm) crystals, and lower-grade conditions recorded by the rims of coarser grains and by finer grains. This observation suggests that re-equilibration during recovery of these rocks to the surface has been important which may account for the discrepancy in estimated P-T conditions. Revised estimates of the equilibration conditions of the Scourian Complex of T 〉 1000°C and P 〉 8.5 kbar are presented. The conditions suggested for the peak of metamorphism mean that the role of anatexis in the genesis of these rocks must be considered and the nature of the fluid phase thoroughly investigated.

Notes: Abstract. Plagioclase porphyroblasts from silvergrey schists belonging to the Nevado Filabride Complex in the Sierra Alhamilla (Betic Zone, SE Spain) are interpreted as having been formed preand synkinematically with respect to the second phase of deformation. Different types of inclusion patterns represent 'snap-shots’(high growth-rate/strain-rate ratio features) of the formation of a diffentiated crenulation cleavage during this second phase of deformation, by the processes of kinking, crenulation and associated differentiation.Regional considerations indicate an Alpine age for this tectono-metamorphic event, which can be explained by the‘hot emplacement’of the higher Nevado Filabride units. The observed structural evolution is not consistent with a pre-Alpine polyphase deformation history.

Notes: Evidence from rock microstructures, mass transfer and isotopic exchange indicates that substantial quantities of aqueous fluids are involved in low- and medium-grade regional metamorphism. Similar conclusions are drawn from many retrograde environments, whereas high-grade metamorphic fluids may be melt dominated. The mobile fluids play essential roles in metamorphic reactions, mass transport and deformation processes. These processes are linked by the mechanical consequences of metamorphic fluid pressures (Pf) generally being greater than or equal to the minimum principal compressive stress. Under such conditions metamorphic porosity comprises grain boundary tubules and bubbles together with continuously generated (and healed) microfractures. Deformation results in significant interconnected porosity and hence enhanced permeability. Lithologically and structurally controlled permeability variations may cause effective fluid channelling.Simple Rayleigh-Darcy modelling of a uniformly permeable, crustal slab shows that convective instability of metamorphic fluid is expected at the permeabilities suggested for the high Pf metamorphic conditions. Complex, large-scale convective cells operating in overpressured, but capped systems may provide a satisfactory explanation for the large fluid/rock ratios and extensive mass transport demonstrated for many low- and medium-grade metamorphic environments. Such large-scale fluid circulation may have important consequences for heat transfer in and the thermal evolution of metamorphic belts.

Notes: Hydrothermal synthesis of Fe-pum-pellyites was conducted using high pressure cold-seal apparatus and solid oxygen buffering techniques at temperatures between 250°C and 550°C and 2.0–9.1 kbar Pfluid. Fe-pumpellyites were synthesized from partially crystalline gel mixtures of compositions: 4CaO - 2.1Al2O3_1.5FeO - 0.3MgO - 6SiO2 (II) and 3CaO - 1.5 Al2O3 - 2.7FeO - 0.3MgO - 6SiO2 (III) in the presence of excess H2O at Pfluid of 5–9.1 Kbar, temperatures between 275°C and 325°C, and fO2 defined by the QFM and HM buffers; for both of these compositions (II and III), the condensed synthetic run products included minor 7Å chlorite ± garnet ± Fe-oxide. The cell dimensions and aggregate refractive index (a= 19.13(2)Å, b= 5.940(4)Å, c= 8.847(5)Å, ±= 97.37(6)±, and n= 1.702(2)) of the pum-pellyite synthesized from the bulk composition II mix are compatible with those of natural pumpellyites containing similar total Fe contents. Attempts at synthesizing Fe-pumpellyites from a Mg-free bulk composition were not successful; these results are consistent with the total absence of natural Mg-free pumpellyites.The higher temperature, higher oxygen fugacity assemblages of the equivalent bulk compositions (II and III) consist of epidote ± minor amounts of chlorite, garnet, quartz, hematite, and magnetite. The results of these synthesis experiments accord with the mineral parageneses observed in low-grade metabasites which imply that Fe-pumpellyites are replaced by epidote with increasing temperature and/or fO2 and that Fe3+ is preferentially partitioned into epidote with respect to coexisting pum-pellyite. In addition, these synthesis experiments indicate that Fe-bearing pumpellyites crystallize at and are stable to lower temperatures than more aluminous pumpellyites—a result also consistent with natural systems.

Notes: This is the first of a series of papers giving the solution of the inverse problem in seismic exploration. The acoustic approximation is used together with the assumption that the velocity field has the form 〈displayedItem type="mathematics" xml:id="mu1" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR998:GPR_998_mu1"/〉 . The forward problem is then linearized (thus neglecting multiple reflected waves) and the inverse problem of estimating δ is set up. Its rigorous solution can be obtained using an iterative algorithm, each step consisting of a classical Kirchhoff migration (hyperbola summation) plus a classical forward modeling step (circle summation).

Notes: The response of a seismic model to excitation by a source can be represented in terms of the action of reflection and transmission operators for portions of the structure. This approach provides a flexible framework for both modeling and processing problems.The operator development provides a physical description of the wave propagation process and, via the expansion of reverberation operators, gives a mechanism for assessing the accuracy of approximate developments. The representation suggests new ways of developing modeling algorithms by balancing the computational effort expended on minor and major features of the model.For processing problems, the operator representation shows the relation of processing stages to the seismic wave field and thereby indicates effective sequences of operations. For migration it is possible to specify an ideal pre-stack migration procedure in terms of the inverse of the propagation operators and to examine the problems which need to be overcome by practical algorithms.

Notes: A numerical method is presented for direct interpretation of resistivity sounding measurements. The early part of the resistivity transform curve derived from field observations by standard methods is approximated by a two-layer curve. The resistivity of the first layer is determined from the arithmetic mean of the successive computations which are carried on each of three successive discrete values of the resistivity transform curve. Using this mean value of the resistivity, the thickness of the first layer is computed from the sample values in pairs of the resistivity transform curve. After these determinations, the top layer is removed by Pekeris's reduction equation. The parameters of the second layer are obtained from the discrete values of the reduced transform curve (which corresponds to the second part of the resistivity transform curve) by the same procedure as described for the first layer.The same computational scheme is repeated until the parameters of all intermediate layers are obtained. The resistivity of the substratum is determined from the reduction equation.

Notes: In odd-depth structure the two-way traveltime to each boundary is constrained to be an odd integer. The odd-depth property of a model is exposed to possible refutation under a seismogram test. Test function is a simple transformation of a synthetic seismogram. For an odd-depth model the test function has identically the value 1.The testability of a synthetic seismogram over an odd-depth structure provides a method of deterministic deconvolution. There is no need of specialized assumptions, like the minimum-phase property, about the source wavelet. The deconvolution may be performed in the absence of the early segment of a seismogram.

Notes: Zero-offset-source VSP surveys provide information about the subsurface only within the Fresnel zone centered at the well. Offsetting the source location moves the reflection zones away from the well thus providing lateral cover.Conventional processing of this type of data gives rise to a distorted image of the subsurface. Using a simple ray-tracing scheme, this image may be reconstructed into the more familiar coordinate system of the surface seismic section. This simple data-independent mapping is based on the assumption of horizontal layering and requires a vertical velocity profile.The technique of placing the source away from the borehole was first applied to the single-offset-source VSP survey. However, data from any survey geometry (such as deviated well with rig source, walkaway VSP, etc.) can be mapped to the coordinate system defined by the appropriate seismic section.To obtain the best results from this type of survey the target area must be defined and simple modeling techniques used to optimize the source location(s). These pre-survey modeling methods may also be used to anticipate—and hence avoid a number of problem areas which experience has highlighted.The data from any VSP survey is the result of a realizable experiment and as such obeys the wave equation. This implies that the wave equation may be used to migrate the data to its true subsurface location. Theoretically, such a process is more secure than ray-tracing techniques, although its practice presents many difficulties.

Notes: A mathematical expression for potential of a direct current point source in an inhomoge-neous anisotropic earth is derived. The coefficient of anisotropy is given by f= (σr/σz), where σr and σz are the conductivities parallel and perpendicular to the bedding plane. It is assumed that σz varies with depth, whereas σr varies transversely. This potential may be useful in interpretation of geoelectrical data in specified geological situations. Master curves for Wenner and Schlumberger configurations are presented

Notes: Total magnetic intensity contour maps for the study region (between 2°E to 10°E and 56°N to 60°N) were digitized and converted to a regular grid of 285 × 285 points. The study area measures approximately 444 km × 444 km and the grid spacing is thus 1. 56 km. The International Geomagnetic Reference Field for 1975 was gridded for the above-used net, and from the two data sets a further grid of the ▵T field was generated. A large number of profiles were constructed which were suitable for depth determinations. The regular grid ▵T data is also convenient for the computation of the second vertical derivative. Using the method of vertical prisms of Vacquier et al. (1963), a large suite of curvature-depth indices was measured to complement the depths obtained from the intensity slopes and from boreholes which reach the crystalline basement. The depth to the magnetic basement has been contoured, and the resulting map is shown to be in good agreement with what is known about the deeper geology of the study area.The work reported here is part of a research project supported by Amoco Norway, BP Petroleum Development Ltd, Elf Aquitaine, Esso Exploration and Production, Norwegian Gulf, Norsk Hydro, Mobil Exploration Norway, Norwegian Petroleum Directorate, Royal Norwegian Council for Scientific and Industrial Research (NTNF), Norske Shell, and Statoil.

Notes: Geophysical inversion involves the estimation of the parameters of a postulated earth model from a set of observations. Since the associated model responses can be nonlinear functions of the model parameters, nonlinear least-squares techniques prove to be useful for performing the inversion. A common type of inversion applies iterative damped linear least squares through use of the Marquardt-Levenberg method. Traditionally, this method has been implemented by solving the associated normal equations in conventional ways. However, Singular Value Decomposition (SVD) produces significant improvements in computational precision when applied to the same system of normal equations. Iterative least-squares modeling finds application in a wide variety of geophysical problems. Two examples illustrate the approach: (1) seismic wavelet deconvolution, and (2) the location of a buried wedge from surface gravity data. More generally, nonlinear least-squares inversion can be used to estimate earth models for any set of geophysical observations for which an appropriate mathematical description is available.

Notes: In this study we derive expressions for particle displacement or particle velocity anywhere inside a stratified earth and at its surface due to horizontal torque source located in the top layer. Equivalently, invoking Green's function reciprocity theorem, the solution applies also to the case of a surface or subsurface source when the resulting displacement or velocity is measured within the top layer.In order to evaluate the closed-form analytical solution economically and accurately it is advisable to introduce inelastic attenuation. Causal inelastic attenuation also lends the necessary realism to the computed seismic trace. To provide proof that the analytical solution is indeed correct and applicable to the multilayer case, a thick uniform overburden was assumed to consist of many thin layers. The correctness of the computed particle velocity response can be very simply verified by inspection. The computed response can also serve as a check on other less accurate methods of producing synthetic seismograms, such as the techniques of finite differences, finite elements, and various sophisticated ray-tracing techniques.It is not difficult to construct horizontal surface torque source. It appears that such source is well suited for seismic exploration in areas with a high-velocity surface layer. A realistic source function is analyzed in detail and normalized displacement response evaluated at different incidence angles in the near and the far fields.In an effort to distinguish the features of an SH torque seismogram from a pressure seismogram two models with identical layerings and layer parameters have been set up. As expected the torque seismogram is very different from the compressional seismogram. One desirable feature of a torque seismogram is the fast decay of multiples.Exact synthetic seismograms have many uses; some of them, such as the study of complex interference phenomena, phase change at wide angle reflection, channeling effects, dispersion (geometrical and material), absolute gain, and inelastic attenuation, can be carried out accurately and effortlessly. They can also be used to improve basic processing techniques such as deconvolution and velocity analysis.The numerical evaluation of the analytical solution of the wave equation as described in this paper has a long history. Most of the work leading to this paper was carried out by one of us (M. J. K.) in the years 1957 to 1968 at the Geophysical Research Corporation. However, the full testing of the various computer codes was carried out only very recently at the Phillips Petroleum Company.

Notes: In South and Central Goa iron ore occurs in two parallel belts with the general NW-SE Dharawar trend. The ore occurrence, however, is not continuous. There are barren zones as well as zones of very high concentration in some of which there are mining activities.Landsat MSS data have been interpreted over a zone covering both mineral belts in order to delineate the ore occurrences. As a guide line a known ore-bearing area has been considered along with the unknown zones.On the basis of two-dimensional plotting of gray level values it has been found that the MSS bands 4 and 7 are most suitable for the studies over iron-rich areas in Goa. Two techniques are described here for the processing of the MSS data; the separation of residual from the regional and MSS band-ratioing. It is observed that (i) the gray level residual maps of MSS bands 4 and 7 are of use in demarcating the iron-ore-bearing zones, and (ii) an existing mine, an abandoned mine, and a proved iron ore zone could be delineated by MSS band-ratioing. On the basis of the latter technique, a few areas with ore occurrence potential have been indicated.

Notes: The most common source of seismic energy is an explosion at some depth in a borehole. The radiated waves are reflected not only at the subsurface layers but also at the free surface. The earth's surface acts as a generator of both P- and S-waves.If the source depth is much less than the dominant wavelength the reflected waves resemble closely the waves generated by a single force. Theoretical seismograms were computed with different methods to look for the relevance of the surface-reflected waves. The numerical experiments show reflected shear waves even for small shotpoint—receiver distances. Due to their polarization these waves can be detected most easily on in-line horizontal geophones. The existence of these waves was examined during a conventional survey in Northern Germany. Conventional data analysis shows a large variability in the νp/νs ratio. The method used here produced a shear-wave section with a rather good signal-to-noise ratio down to 4 s S-wave reflection time.

Notes: We describe the concept of physico-geological models (PGMs) in geophysical exploration. They represent a “general model”, a spatial combination of a set of particular models (disturbing bodies). The modeling is called complete, incomplete or approximate, depending on the degree of characterization of the PGM by parameters such as dimension, shape and petrophysical property. Each of the three modeling types can be realized as a conceptual, and analytical, or a material PGM. Both deterministic and stochastic PGMs exist; deterministic models are mainly used to investigate the possibilities of a geophysical method, while stochastic models serve to substantiate complex geophysical interpretations.Depending on the geological problem, PGMs are subdivided into multi-alternative models (geological mapping, prediction, general prospecting) and double alternative models (specialized prospecting).An exploration-oriented classification of the PGMs of mineral deposits is discussed. According to this classification the variety of known genetic deposit types is reduced to a limited number of generalized PGM types. The development of typical PGMs is illustrated with examples of magnetitic deposits of Siberia.

Notes: An inverse problem is one in which the parameters of a model are determined from measured seismic data. Important to the solution of inverse problems is the issue of whether or not a solution exists. In this paper we show, in a constructive manner, that a solution does exist to the specific inverse problem of determining the parameters of a horizontally stratified, lossless, isotropic and homogeneous layered system that is excited by a non-normal incidence (NNI) plane wave. Mode conversion between P- and S-waves is included.We develop a seven-step layer-recursive procedure for determining all of the parameters for layer j. These parameters are P-wave and S-wave velocities and angles of incidence, density, thickness, traveltimes, and reflection- and transmission-coefficient matrices. Downward continuation of data from the top of one layer to the top of the next lower layer is an important step in our procedure, just as it is in normal incidence (NI) inversion. We show that, in order to compute all parameters of layer j, we need to (and can) compute some parameters for layer j+ 1. This is a non-causal phenomenon that seems to be necessary in NNI inversion but is not present in NI inversion.

Notes: The suitability of VLF ground surveys in the investigation of shallow two-dimensional structures is analyzed. For such structures the polar formalism is derived, necessary in practice since the transmitters are generally not in the structural strike or profile. A simple vertical dike is considered to demonstrate the striking anisotropy which can be expected over such a structure, in particular the high apparent resistivity along the direction of a well-conducting dike and the low resistivity across it. The theory is then confronted with the practical example of an asymmetrical vertical dike resulting from a strike-slip fault. Modelling of the survey results is very successful and yields good confirmation of the polar behaviour. VLF ground surveys thus provide a quick and powerful tool for the study of geological accidents within about 100 m of the surface.

Notes: A new mode of operation for the Turam electromagnetic exploration system is proposed in which the transmitter loop is placed across the expected trend of a conductor and the receiver is operated along lines parallel to one side of the transmitter. The concept appears to offer several benefits which include greatly extended traverse length, the use of large coil spacing, rejection of the effects of conductive environments, and consistency in the indication of target dip.

Notes: Determination of thickness of sediments (usually of high conductivity) overlying a high-resistivity basement is one of the basic problems of electrical exploration methods. This paper proposes to determine horizontal electrical conductance on the basis of impedance calculated from electrical and magnetic fields of distant quasi-static (low-frequency) point sources. Using the proposed method, horizontal conductance of the sediments can be determined also from artificial quasi-static noise-impulses coming from sources of unknown position and intensity. The results of analogue modeling and field examples prove the potential of the proposed technique.

Notes: Since the early years of seismic surveying, field engineers and observers have been faced with the appearance of interference at power-line frequency (and harmonics) on seismic records: 60 Hz in the Americas, 50 Hz in Eurasia; and 16 2/3 Hz (the train supply frequency in Germany, Austria, and Switzerland) in part of Europe. Such interference is picked up from the ground surface through electric leakage between geophone leads and ground.To minimize power-line interference, the autobalancer automatically performs a balancing procedure by adjusting two potentiometers on each channel—essentially the same operation as manually performed by seismic operators for decades. The first 24-channel system became operational in 1977.The first-generation design comprised two 256-position electronic switching potentiometers. Over 4 years of experience with this system in all parts of the world has shown its value under a variety of field conditions. Improvements up to 30 dB were often found, sometimes even up to 40 dB. Balancing time was about 30 s.A second-generation design was field-tested in 1981. Potentiometer stepping is now performed at three levels: coarse, medium, and fine. This leads to faster operation and better resolution: balancing time is about 10 s; improvements up to 60 dB are occasionally found.A later modification gives a further reduction of the balancing time to 5 s in 60 Hz areas or 6 s in 50 Hz areas.

Notes: The technique of digital linear filtering is used for transformation of apparent resistivity data from one electrode configuration into another. Usually filter spectra are determined via the discrete Fourier transforms of input and output functions: the filter characteristic is the quotient of the spectra of the output function and input function.In this paper, the transformation of the apparent resistivities is presented for four electrode configurations (Wenner, the two-electrode, Schlumberger, and dipole configurations). In our method, there is no need to use the discrete Fourier transform of the input and output functions in order to determine the filter spectrum for converting apparent resistivity in one electrode configuration to any other configuration. Sine responses for determination of the derivative of apparent resistivities are given in analytical form.If the filter spectrum for converting the apparent resistivity to the resistivity transform for one electrode configuration is known, the filter spectra for transforming the apparent resistivity to the resistivity transform for any electrode configurations can be calculated by using newly derived expressions.

Notes: In the western coal-mining area of Ruhrkohle AG, reflection seismic prospecting for the Carboniferous coal measures is severely impaired by structures with halokinetic features. These structures make the interface between Mesozoic and Paleozoic layers, i.e., the top of Zechstein in general, very rugged. Unfortunately the velocity contrast at this interface is very high in that area, the ratio of velocities being 1.5 to 2.0. Therefore, migration and stacking become a problem.Three types of migration are presented:〈list xml:id="l1" style="custom"〉1(f, x)-time-migration with vertical time-to-depth conversion as a second step.2Kirchhoff migration down to a level determined approximately by the highest points of the top of Paleozoics, i.e., 0.35 s, and Kirchhoff-downward continuation for all times exceeding 0.35 s. Intermediate static corrections for these latter times with subsequent (f, k)-time-migration and final vertical time-to-depth conversion.3Direct depth migration in the (f, x)-domain using three interval velocities.In all cases an intermediate picking of the velocity interfaces is necessary. In case 2 this occurs at an earlier stage of the process than in case 1, and in case 3 at a still earlier stage.The results of the second and third migration procedures are superior to those of the first. Possibilities for misinterpretation of faults are reduced considerably when the second or third migration procedure is applied.

Notes: Airborne electromagnetic methods are most commonly used in mineral exploration. However, new developments, such as multifrequency capability and digital on-board field recording, as well as improvements in instrumentation resulting in high signal-to-noise ratios in recorded data, have made their application in geological mapping possible.A three-frequency airborne EM survey carried out over an area northwest of Timmins, Ontario, was interpreted in terms of thickness and resistivity of the layers of a two-layer earth section. Since both in-phase and quadrature components are measured, this provides six independent parameters at each point in space. Based on prior geological information and a preliminary interpretation of the field records, two two-layer models of the subsurface seemed to be appropriate for most of the survey area. An automatic computerized interpretation procedure was devised to interpret the field data at each point in terms of thickness and resistivity parameters of those two models. When the geology is more complex, the data do not fit the models and no interpretations are made. Two maps illustrating the variation of the resistivity and the thicknesses of the layers were constructed from the interpreted data. These maps agree with the known geological information about the distribution of glacial clay in the area. Areas where the layered models do not fit are known to be areas where the geology is complex with a large number of dykes and other lateral inhomogeneities. The study shows that multifrequency airborne EM surveys can be very useful in geological mapping over inaccessible terrain and can significantly help the mapping geologist where outcrops are scarce.

Notes: Shear (= S) wave studies in addition to compressional (= P) wave surveys have revealed that S-wave velocities are much more variable than P-wave velocities. This strongly affects Poisson's ratio σ, especially in young sediments. It is shown that σ has a great influence on the directivity pattern, i.e., on the radiation and receiving characteristic of horizontal sources and receivers. For their calculation, well-known surface boundary conditions and White's statement of the reciprocity relation (White 1965) are used. They also form the basis for the combination of source and receiver patterns essential for practical field work. The various combinations of horizontally and vertically arranged sources and receivers are investigated with regard to their directivity pattern, their amplitude ratio and, in general, to their usefulness in seismic prospecting.

Notes: Previously ignored characteristics of the seismic recording instrument are presently experienced as limitations as more sophisticated interpretive methods using wider frequency ranges are developed to extract stratigraphic information from seismic land data for hydrocarbon and mineral exploration. Most of these limitations arise from inadequate characteristics of the first element of the seismic instrument: the geophone. A geophone does not faithfully follow the motion of the earth for higher frequencies due to poor geophone-earth coupling. This filtering effect brings about time shifts that are dependent on the frequency and the soil type. A geophone can also produce spurious outputs, brought about by the motion of the suspended part of the geophone, with a magnitude comparable to that of the desired output. The suspension is made very compliant to obtain the required sensitivity. A compliant suspension, however, gives a large sag. The geophone can therefore only be used in one position, tolerating little tilt. A compliant suspension also widens the traveling range of the movable part. Minor sensitivity changes with travel are then noticeable as nonlinearity, since the surface wave is large with respect to the reflected wave. A compliant suspension is usually realized in the form of thin, spirally shaped spring-spiders. Such suspensions exhibit transverse or rotational resonances that are in or close to the seismic frequency band. Excited by ground roll, they can produce considerable undesirable output.The novel geophone we describe is a light-weight (17 g) acceleration-sensitive transducer which gives good ground coupling and partial correction for the increasing damping in the earth with increasing frequencies. It employs internal hybrid electronics for a magnetodynamic velocity-nulling feedback system. Velocity nulling makes the movable part of the geophone virtually rigid with respect to the housing. This makes the geophone characteristics independent of the suspension. The springs used are stiff in a transverse and rotational direction so that the suspension resonances are well outside the useful frequency band. This suspension also allows the geophone to be used in any orientation while being only sensitive to the vibration component along the main axis. The feedback system makes the sensitivity flat within 1 dB from 2 Hz to 500 Hz, with a phase tolerance smaller than 5°. The geophone is robust, has no moving internal wires, employs a current output [sensitivity 1 mA/(m s−2)] and internal gain so that the signal-to-cable-noise ratio is improved. This type of output allows parallel connection without any interaction between the geophones.

Notes: There are two types of masked layers in seismic refraction work: the velocity reversal (low-velocity layer) and the hidden layer (insufficient velocity contrast or layer thickness). On the basis of an analytical formulation of the general case of a masked layer under an overburden of plane and parallel multiple refractors the two limiting cases are discussed: the solution resulting from an uncritical interpretation of the measured time-distance curve and the blind zone solution. Between these two limiting cases there is a variety of possible masked layer solutions. These no-blind zone solutions—as well as the blind zone solution itself—are formulated separately for the velocity inversion and the hidden layer case.For the evaluation of some no-blind zone solution a diagram is presented which can be used for any case of multiple refractors in the overburden of the masked layer. However, it is only for the three- and the four-layer case that a blind zone interpretation by use of diagrams is advisable. Such diagrams are presented together with the basic sets of formulae which contain as parameters only ratios of velocities and layer thicknesses. As the velocity of the masked layer is usually unknown the diagrams are principally constructed to show the dependence on the masked layer velocity. This is useful for estimation of the largest possible error.

Notes: Wave equation migration techniques have shown the limits of traditional stacking methods with data from tectonically complicated areas. An improved stack can be obtained utilizing the dip-moveout correction technique based on offset continuation. The properties and the limits of the algorithms used are summarized briefly.Several synthetic and real data examples are shown and compared with the results obtained using conventional processing in order to show the focusing effects and the strong improvement in signal-to-noise ratios, both at the stacked and migrated section level. The possibility of exploiting this technique to transform multiple coverage into increased spatial resolution is illustrated with examples.

Notes: The far-field signatures from a comprehensive and systematic airgun pulse test have been analyzed. Empirical relations between the characteristic signature parameters and depth (5–12 m), pressure (100–137 bar = 10–13.7 MPa) and total chamber volume (0.65–9.5 l) have been derived. Also, the influence of using waveshape kits in different positions within the chamber has been tested.The results indicate that:〈list xml:id="l1" style="custom"〉1 The amplitude is proportional to chamber pressure to the power 3/4.2 The bubble period is nearly independent of the position of the waveshape plate.3 The increase in primary/bubble amplitude ratio is inversely proportional to the chamber volume above the waveshape plate.4 The amplitude is independent of airgun depth.Suggestions and comments regarding this work from Dr B. Ursin and Dr A. Ziolkowski are appreciated. The field work was supported by the Norwegian Petroleum Directorate through the Continental Shelf Project at the Seismological Observatory, University of Bergen.An airgun allowing for continuous variation of the chamber volumes was supplied by GECO (Geophysical Company of Norway). The purchase of two airguns was financed by Norske Getty Exploration A/S.

Notes: The accuracy of short length digital linear filter operators can be substantially increased if the sampling interval as well as the abscissa shift are properly adjusted. This may be done by a trial and error process of adjustment of these parameters until the error made by the filter operator, applied to a suitably chosen test function, is smallest.As an illustration of the application of this method, 7-, 11- and 19-point filters for the calculation of Schlumberger apparent resistivity from a known resistivity transform are designed. Errors with the new 7-point filter are seen to be less than those with a 19-point filter of conventional design. The errors with the new 19-point filter are two to three orders of magnitude smaller than those made by the conventional 19-point filter.The new method should provide digital linear operators that allow significant improvements in accuracy for comparable computation efforts, or substantial reduction in computation for comparable accuracy of results, or something of both.

Notes: A well-known technique for the migration of normal-incidence two-way travel-time maps is extended to common-source-point travel-time data. The travel time and the travel-time gradient are used to compute the parameters defining the tangent plane of the reflecting interface. It is also shown how the curvature matrix of the received wavefront can be used to compute the curvature of the reflecting interface. The method is initially derived for common-source-point data and then extended to common-midpoint data.In a three-dimensional medium the wavefront curvature matrix is computed by solving a 2 × 2 symmetric matrix Riccati equation. In a two-dimensional medium and in a medium with constant velocity gradient, the wavefront curvature matrix is computed by solving a scalar Riccati equation and two linear equations. The migration procedures are also simplified.When the velocity function is unknown, the migration procedures cannot be used. An inverse modeling algorithm which simultaneously performs the migration and estimates the velocity function must then be applied. Two different inversion schemes are discussed briefly.

Notes: Time-domain-induced polarization (IP) laboratory measurements were performed on about 200 fine sediment samples with varying water content. The results permitted an analysis of IP properties of clays, loams, silts, and sands.Particular emphasis has been given to the analysis of the chargeability m as a function of lithotype and the water content.By analyzing decay curves, a new parameter was identified. It is a statistically specific characteristic of the lithotype and is independent of the water content. Therefore, it provides a diagnostic parameter for lithotype identification. In association with the values of chargeability and electrical resistivity, this parameter permits a reliable evaluation of water content and yields useful information about the porosity and permeability of the lithotype.

Notes: Approaches to the reduction of bias in the computation of the elements of the magnetotelluric impedance tensor have been proposed in the past by several authors. In this paper a clear distinction is made between random errors and bias errors. No effort is made to reduce either, but the emphasis is on their estimation. Both types of errors depend critically upon the polarization of the magnetic field. The random error increases with increasing noise-to-signal ratio in the electrical field, and it is rather insensitive to noise in the magnetic field. The bias error increases with increasing noise-to-signal ratio in the magnetic field. Expressions for random errors and maximum bias errors are developed and discussed using a single station set-up. Random errors with a reference station set-up are also calculated.

Notes: Dipole-dipole induced polarization (IP) data are displayed typically as multi-level profiles, or as contours on vertical sectional plots referred to as pseudo-sections. The dipole-dipole array tends to yield IP anomalies in which the most anomalous values are displaced laterally from the source body. The data patterns are fairly interpretable on pseudo-sections or on multi-level profiles but are sufficiently complex to discourage the contouring of the data in plan.A method was developed for the presentation of dipole-dipole IP data on a contour map. The method consists of a simple averaging of data which can be performed manually if desired. It yields a single output value per station which reflects all levels of the pseudo-section, and is suitable for contouring in plan. The advantage of the technique is that it provides a quantitative picture of IP anomalies in their background or regional setting.

Notes: Constant offset sections can be mapped to a fixed offset and compared in order to provide a method of velocity analysis. The direct mapping to zero offset prior to stack might provide an alternative processing procedure to NMO and stack. The main advantage of such a procedure would lie in the correct treatment of cross-dips, but interpretational advantages might also follow from the performance of partial stacks biased to either high or low offset information.

Notes: Induced polarization and resistivity model studies over thin dykes for varying resistivity contrasts, depth of burial, and dip angles show striking parallelism of the surface apparent resistivity contours with the boundary of the body. This effect may be utilized for the estimation of strike length of the body. Results show that intermediate values of the electrode spacing is satisfactory for detection under widely varying conditions of resistivity contrast and depth. The percentage frequency effect (P.F.E.) and the metal factor (M.F.) responses are found to be more sensitive to the variation in the depth of burial than the resistivity responses. Pseudosections for P.F.E. and M.F. are concentrated much closer to the body than the resistivity pseudo-sections.

Notes: A direct interpretation scheme is developed which is capable of determining most of the geological features of a ground which can be assumed to be two dimensional in structure. This scheme extends the earlier work of Pekeris (1940) and Koefoed (1968) to the case where the basal layer of a ground is undulating. It also has a limited use for finding the parameters of a dipping dyke in the lower medium. Though the top and dip of the dyke can be determined, this is not true for the thickness.

Notes: Laboratory studies of absorption-frequency behavior in rocks often use spectral ratios of digitally recorded ultrasonic signals which have been transmitted through a rock sample and a reference sample of very low absorption, respectively. It is proposed to treat the digitally recorded signals as an autoregressive-moving average (ARMA) process which, using recursive filter concepts, can be represented as a ratio of two polynomials in the z-transform variable z. The numerator polynomial contains only that part of the signal that is modified by anelastic effects, whereas the denominator contains the elastic effects of the physical apparatus such as reverberations. Examples are given which show that this separation of the recorded signal greatly facilitates the laboratory investigation of loss mechanisms and absorption-frequency behavior based on spectral ratios.

Notes: The generalized integral equation for the electric potential governed by a quasi-harmonic equation can be derived via a variational formulation. For surface current distributions it is not always a Fredholm integral equation of the second kind. Numerical solutions of the general heterogeneous problem can be obtained with the “reciprocal averaging technique”, where the solution is obtained a second time after exchange of source and field points.

Notes: A method to calculate the resistivity transform of Schlumberger VES curves has been developed. It consists in approximating the field apparent resistivity data by utilizing a linear combination of simple functions, which must satisfy the following requirements: (i) they must be suitable for fitting the resistivity data; (ii) once the fitting function has been obtained they allow the kernel to be determined in an analytic way.The fitting operation is carried out by the least mean squares method, which also accomplishes a useful smoothing of the field curve (and therefore a partial noise filtering). It gives the possibility of assigning different weights to the apparent resistivity values to be approximated according to their different reliability.For several examples (theoretical resistivity curves in order to estimate the precision of the method and with field data to verify the practicality) yield good results with short execution time independent of shape the apparent resistivity curve.

Notes: We present a new method for the extraction and removal of the source wavelet from the reflection seismogram. In contrast to all other methods currently in use, this one does not demand that there be any mathematically convenient relationship between the phase spectrum of the source wavelet and the phase spectrum of the earth impulse response. Instead, it requires a fundamental change in the field technique such that two different seismograms are now generated from each source-receiver pair: the source and receiver locations stay the same, but the source used to generate one seismogram is a scaled version of the source used to generate the other. A scaling law provides the relationship between the two source signatures and permits the earth impulse response to be extracted from the seismograms without any of the usual assumptions about phase.We derive the scaling law for point sources in an homogeneous isotropic medium. Next, we describe a method for the solution of the set of three simultaneous equations and test it rigorously using a variety of synthetic data and two types of synthetic source waveform: damped sine waves and non-minimum-phase air gun waveforms. Finally we demonstrate that this method is stable in the presence of noise.

Notes: In this paper a theorem is demonstrated which allows—after the introduction of a suitable dipole kernel function or dipole resistivity transform function—to write the apparent resistivity function as an Hankel transformable integral expression.As a practical application of the theorem a procedure of quantitative interpretation of dipole soundings is suggested in which the dipole resistivity transform function obtained after inversion of the original dipole apparent resistivity data is used to control the goodness of the set of layering parameters which have been derived with our previous method of transformation of dipole sounding curves into equivalent Schlumberger diagrams.

Notes: The Hankel transform theorem can be applied to the inversion of gravity data for the buried sphere, the horizontal cylinder, and the vertical rod. This new approach leads to exact solutions of the transforms for the assumed bodies. A comparison with the classical procedure by Fourier transform reveals that for the sphere and the vertical rod, the Hankel transform is preferable.

Notes: Air guns have been used in various applications for a number of years. They were first used in coal-mining operations and were operated at up to 16000 psi charge pressures. Later, single air guns, operated at 2000 psi, found application as an oceanographic survey tool. Air gun arrays were first used in offshore seismic exploration in the mid-1960's. These early arrays were several hundred cubic inches in total volume and were operated at 2000 psi; they were either tuned arrays or several large guns of the same size with wave-shape kits. Today's arrays have total volumes greater than 5000 cu in. and are typically operated at 2000 psi. Recently, higher-pressure, lower-volume arrays operated at 4000–5000 psi have been introduced; guns used in these arrays are descendants of the coal-mining gun.On first thought one would equate increased gun pressure linearly with the amplitude of the initial pulse. This is approximately true for the signature radiated by a “free-bubble” (no confining vessel) and recorded broadband. The exact relation depends on the depth at which the gun is operated; from solution of the free-bubble oscillation equation, the relation is 〈displayedItem type="mathematics" xml:id="mu1" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu1"/〉 If Pc,1= 6014.7 psia, Pc,2= 2014.7 psia and PO, 1=PO, 2= 25.8 psia (corresponding to absolute pressure at 25 ft water depth), then 〈displayedItem type="mathematics" xml:id="mu2" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu2"/〉 Experiments were conducted offshore California in deep water to determine the performance of several models of air guns at pressures ranging from 2000 to 6000 psi and gun volumes ranging from 5 to 300 cu in. At a given gun pressure, the initial acoustic pulse Pa correlated with gun volume Vc according to the classical relation 〈displayedItem type="mathematics" xml:id="mu3" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu3"/〉 For 1 ms sampled data the ratio 〈displayedItem type="mathematics" xml:id="mu4" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu4"/〉 varied between 4.5 and 5.5 dB depending on gun model. Pulse width of the 2000 psi signatures indicated they are compatible with 2 ms sample-rate recording while pulse width of the 6000 psi signatures was greater, indicating they are less compatible with 2 ms sample-rate recording.Conclusions reached were that 2000 psi air guns are more efficient than higher pressure guns and are more compatible with 2 ms sample-rate requirements.