ALBERT

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: 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: 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: The paper presents research and development at the Thematic Information Service of the Natural Environment Research Council (NERC) on the topic of integrated geographic information systems. The term integrated is used to describe the incorporation of vector digitised map data and raster scanned remotely sensed data within a single computational environment and with a single interface to the user. Point sampled and statistical data can also be handled by the system. The general structure and working of the Thematic Information Service (TIS) is described only to the extent that it is relevant to providing the background to the research or to highlighting why the specific approach described has been taken.

Notes: A procedure for digital image correlation is described which is based on least squares window matching. The immediate aim is high precision parallax assessment, point transfer, and point measurement. Experiments and theory have confirmed the high accuracy potential of the method. By implementation of charge coupled device (CCD) video cameras in an analytical plotter, an experimental hardware and software configuration has been established with which the operational on line application of digital image correlation for conventional photogrammetric measuring tasks can be tested. First results of calibration and performance of the system are presented. They allow optimistic conclusions as to the further development and practical application of digital image processing in photogrammetry.

Notes: Book review in this aticle:DESERTS AND ARID LANDS. Edited by F. El-Baz.IMAGES OF THE WORLD AN ATLAS OF SATELLITE IMAGERY AND MAPS. By C. Mueller-Wille with an interpretive supplement by R. M. Smith.ELEMENTS OF PHOTOGRAMMETRY. By P. R. Wolf.

Notes: The essential features of equipment for mobile processing are outlined with reference to the development of small photographic laboratory vehicles for use in the Jield. The resulting air film processor, the Wainco Products Limited 392. is described and details are given of trials carried out to establish its performance for a range of tasks in both reconnaissance and survey. The tests cover the sensitometric and drying behaviour of several air films in current use and include an analysis of the dimensional stability o f j l m in the processor. Applications using direrent containers and eflectiveness in improvised situations are also described and possible future variants of the processor are suggested.

Notes: The author reviews the historical background of the analytical orthoprojector and examines the design requirements for this instrument. The construction and performance of the Zeiss (Oberkochen)Z2 orthoprojector are described and an assessment is made of its practical advantages.

Notes: The new APY (Analytical Photogrammetry Yzerman) system is described and a detailed explanation is given of the principles involved. The APY depends on a unique modular design. A list of the straightforward operations required for its application in data compilation and revision is included.

Notes: The application of phorogrammetry to the measurement of deformation is discussed and four different approaches are identified and compared. References to and summaries of recent work are made.

Notes: The development and operation of an inexpensive retractable vertical camera rig is described. Unlike several commercially available camera packages, the cameras on this rig are accesible in flight thus permitting the changing of magazines and other operations. Since the camera carriage is only deployed when over the site to be photographed, the aircraft performance is not degraded during positioning flights.

Notes: An individual view is provided of the function of the art and science of photogrammetry, over the whole of the history of attempts to make measurements from optical observations. The predominating influence, in past years, of the need for mapping from air photographs is compared with the prospect that increasing demands for similar facilities in close range measurement, for application in fields such as mechanical engineering, medicine and architecture, may soon be seen to give further encouragement to the development of the subject. The author makes use of illustrations from his experience to point out simple cases where the potential of photogrammetric techniques, in the wider sense, has not been appreciated.

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: 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: Experiments were carried out on seven sites in Fermanagh to investigate the yield loss caused by leatherjackets in grassland. Populations on the sites ranged from 65,000 to 865,000 ha−1 and significant damage was recorded from all sites except the least populated. Greater yield increases were obtained by controlling leatherjackets in September when compared with control in March. Regression models were fitted and it was found that 125,000 leatherjackets ha−1 in March caused a yield loss of 50 kg herbage dry matter (DM) ha−1 by mid-May. Based on these figures the average yield loss to leatherjackets in Northern Ireland, from 1965 to 1982, was 208 kg DM ha−1. Controlling leatherjackets in September, rather than March, increased the potential avoidable yield loss by a factor of 2.72. The average yield loss at first silage cut is therefore 566 kg DM ha−1. On average, 100,000 ha grassland in Northern Ireland may be suffering an annual loss of 1t herbage DM ha−1.

Notes: The utilization of sown and indigenous plant species was studied in three experiments in which plots were stocked with similar live weights of sheep and goats. In the first experiment the animals grazed plots containing 0–5 ha of rush (Juncus effusus)- infested reseeded pasture and 0–5 ha of unimproved blanket bog. The second and third experiments took place on old rush-infested improved pasture; in one experiment two levels of herbage mass of grass were provided while in the other the rushes were cut in spring or remained uncut.The goats grazed the rushes readily in all three experiments. Reduction in herbage mass of grass increased utilization of rushes by goats although these animals still grazed rushes when grass supply was plentiful. Sheep scarcely grazed J. effusus even when hard-pressed by shortage of grass. Both sheep and goats grazed J. acutiflorus.In Experiment I sheep utilized reseeded pasture more heavily than did goats. Sheep grazed similar proportions of the leaves of grass and clover while goats grazed a lower proportion of clover leaves as compared with grass. Among the grasses sheep discriminated in favour of Lolium perenne whereas goats did not. On the blanket bog vegetation, Eriophorum vaginatum and Calluna vulgaris were grazed both earlier in the season and more heavily by goats than by sheep.The relationships between sward structure and grazing height in accounting for differences in species selection by sheep and goats are briefly discussed. The possibility of using goals to control coarse weeds in hill pasture and for strategic grazing to manipulate floristic composition is outlined.

Notes: The performances of a rising-plate meter and a single-probe electronic capacitance meter were compared for estimating the herbage mass of irrigated pastures both before and after grazing by dairy cows. Mean CV for the rising plate meter and the single-probe meter were 12.7% and 13.3% respectively when herbage mass was measured before grazing, and 21.8% and 15.4% respectively immediately after grazing. These coefficients of variation indicate that while the performance of the two meters was similar when herbage mass was measured before grazing, the single-probe meter was a more precise instrument for estimating herbage mass after grazing.Trampling of herbage prior to measurement had a major effect on the regressions of both meters; the regression intercept was increased by 2800 kg dry matter (DM) ha−1 for the rising-plate meter and by 2580 kg DM ha−1 for the single-probe meter. The slopes of the graphs, however, remained unchanged.It was concluded that while the single-probe meter performed better than the rising-plate meter for estimating the mass of post-grazing herbage when trampling was negligible, the large effect that trampling had on the post-grazing calibration regressions precluded both instruments from general use in dairy cattle research. This indicates that another technique is required in circumstances in which pastures are heavily trampled.

Notes: Only a very limited number of species, including red fescue (Festuca rubra), perennial ryegrass (Lolium perenne), timothy (Phleum protense) and white clover (Trifolium repens) are considered as suitable for wet hill land. Detailed comparison of red fescue, perennial ryegrass and timothy on contrasting acid brown earth and stagnogley soils from 1966 to 1969 at Pant-y-dŵr Hill Centre (305 m OD) highlighted the advantages of red fescue in annual DM production and its more even distribution over a longer growing season but also indicated drawbacks of slow initial establishment and lower feed quality. In vitro DOMD of red fescue averaged 645 g kg−1 over the year compared with 690 g kg−1 for perennial ryegrass and timothy. Assessment of a narrow range of species under surface seeding conditions on a stagnogley soil from 1971 to 1975 confirmed the value of red fescue in DM production and persistence. The results are discussed in relation to the characteristics required in herbage species for hill land improvement.

Notes: High-yielding grass-legume mixtures play an important role in forage-animal systems but finding compatible and adapted species can sometimes be difficult. The objective of this study was to examine productivity of perennial ryegrass and orchardgrass in pure stands with N and in mixtures with legumes.Broadcast sowings were made on conventionally prepared seedbeds in August 1979 and May 1980 on a Hagerstown silt loam soil (fine, mixed mesic Typic Hapludalf). When sown alone orchardgrass cv. Pennlate and perennial ryegrass cv. Reveille received rates of N ranging from 0 to 448 kg ha−1 a−1. The grasses were also grown in mixtures with alfalfa cv. Arc, red clover cv. Arlington, or birds foot trefoil cv. Viking. The same legumes were sown alone. Four legume seed rates were used in both the pure legume and mixed stands.Orchardgrass-N swards were more productive than ryegrass-N swards over 3 years. More dry matter (DM) was harvested from pure stands of

Notes: Fungicide applied to plots of perennial ryegrass cv. S24 in the autumn of 1981 significantly increased tiller number and total dry matter accumulation in February and March of 1982. However, fungicide application in the autumn of 1982 failed to increase tiller number or dry matter accumulation in the spring of 1983. It is suggested that fungicide application in autumn 1981 prevented the death of tillers caused by fungal attack after snow covered the plots in December 1981.

Notes: A mathematical model of grass production and utilization on a dairy farm is described. Using the model an assessment has been made of the extent to which financial uncertainty arising from year-to-year variability in grass yields, coupled with a preference among farmers for minimizing risks, may explain the relatively low stocking rates and observed nitrogen usages on many dairy farms in England and Wales. The degree of risk has been equated with the probability of profits in a particular year being less than those required to cover the consumption needs and short-term borrowing requirements of the farmer. The results of the analysis indicate that a strategy of minimizing risks may lead to a significantly lower stocking rate than one of maximizing profits. Thus, considerations of risk may lead to stocking rates which are suboptimal from the viewpoint of economic and biological efficiency. However, there is no evidence to suggest that the comparatively low average usage of nitrogen on dairy farms in the UK is determined by risk considerations. On the contrary, increasing nitrogen usage lowers the apparent financial risk at a given stocking rate.

Notes: The herbage dry matter (DM) harvested, N response and herbage quality of monocultures and mixtures of Holcus lanatus (local ecotype) and Lolium perenne (cv. S24) were compared in a field experiment. The species were seeded in the proportions 100:0, 75:25, 50:50, 25:75 and 0: 100; three rates of N were applied (150, 300 and 450 kg ha−1) and five or six cuts were taken in each of four harvest years. H. lanatus established poorly in the seeding year and in the first full harvest year L. perenne contributed 1, 54, 76, 89 and 99% of the total herbage DM harvested from the five sward types (mean of N treatments). Only a small decline in the L. perenne content of the swards occurred during the experiment.Total amounts of herbage harvested were significantly lower for H. lanatus monoculture than for the other sward types in three of the four harvest years However, L perenne monoculture yielded significantly more total DM than mixtures with H lanatus in only one of the four harvest years. The five sward types did not differ overall in response to N fertilization. Nitrogen, phosphate and potassium content of the herbage was greatest for H. lanatus monoculture and decreased with increasing proportions ofL perenne present in the mixtures.Thus the presence of a high proportion of H. lanatus in mixture with L. perenne cv. S24 had

Notes: Two experimental populations of Italian ryegrass selected for high and low magnesium content were established as swards, harvested fresh daily and offered to sheep. Chemical analysis of the crops showed a 44% higher Mg and a 22% higher Ca concentration in the high-Mg selection but no significant differences between the selections for the other minerals.Blood and rumen contents were sampled throughout the feeding period and voluntary intake, digestibility and mineral balances for Na, K, Mg, Ca and P were determined. Magnesium intake, apparent availability and retention were significantly greater with the high-Mg selection. There were no significant differences in the plasma-mineral concentrations and all were within the normal range. Rumen fluid from animals fed the high-Mg selection had a higher proportion of soluble Mg and a higher Na:K ratio in the supernatant liquor.

Notes: A perennial ryegrass sward was managed by continuous stocking with sheep (April–September) for 4 successive years after sowing. The sward was grazed to maintain a leaf (lamina) area index (LAI) close to 1.0. Areas of the sward were released from grazing on three occasions: once during summer in the third year after sowing, and twice during spring and summer in the fourth year after sowing. There were marked changes in the structure and physiology of the continuously stocked sward following release from grazing. After several successive years of continuous stocking, the sward comprised a large population of small tillers and the small LAI resulted in consistently low rates of photosynthesis. Following release from grazing, photosynthesis increased markedly as the LAI increased but this change was associated with the loss of a large proportion of the population of tillers. There were seasonal differences in the pattern of changes in photosynthesis and tiller numbers following release from grazing which were not apparent under continuous stocking. The changes in the structure and physiology of the sward following release from grazing suggest that the net accumulation of herbage in areas of sward from which the animals are excluded, for instance using cages, may be an unreliable estimate of production under continuous stocking.

Notes: A randomized block experiment was conducted to compare unwilted and wilted grass silages and the effects of the feed additive monensin sodium on the silage intake and performance of finishing beef cattle. Two regrowths from a predominantly perennial ryegrass (cv. S24) sward were ensiled either without wilting or after field wilting for 3 d (dry matter (DM) concentrations 161 and 266 g kg−1 respectively). Both silages were treated with formic acid (2·6 and 30 litre t−1 respectively) and were well preserved. The silages were offered ad libitum to forty-eight Charolais-cross cattle (thirty-two steers and sixteen heifers, mean initial live weight 351 kg) for 145 d. All animals received 2·2 kg concentrates per head daily and half of those on each silage treatment received in addition 200 mg monensin sodium per head daily. Silage DM intake was 5.04, 504. 5·48, 5·63 ± 0.134 kg d−l; fasted liveweight gain was 0·69, 0·77. 0·64 and 0·73 ± 0.033 kg d−l and carcass gain was 0·47, 0·50, 0·40 and 0·45 ± 0·020 kg d−1 for the unwilted silage without and with monensin and the wilted silage without and with monensin respectively. It is concluded that wilting grass of low DM concentration for 3 d prior to ensiling reduced the performance of finishing beef cattle below that obtained from well-preserved unwilted silage in spite of a higher DM intake being achieved with the wilted silage. The inclusion of monensin sodium in a silage-based diet increased performance without significantly affecting feed intake.

Notes: Formic acid, sodium diacetate, propionic acid and two of its ammonium salts, ammonium isobutyrate and tributyl phosphate, were examined for their antimicrobial effects with cultures of micro-organisms typical of the microflora of moist hay. The assays, which were conducted at pH 5, 6, 7 and 8, revealed that all potential preservatives discriminated against actinomycetes in particular and were least effective against fungi. In addition, it was apparent that antimicrobial properties were enhanced under slightly acid conditions and for this reason the use of mixtures of organic acids and salts in commercial hay preservatives is justified. Tributyl phosphate not only had the strongest antimicrobial properties of the compounds tested but was also the least affected by pH.

Notes: Yields of grass were measured on three dairy farms in 1981. Farm 1, on free-draining loamy soil over chalk, produced 10.3 t DM ha−1 using 336 kg N ha−1. Farm 2, on poorly-drained loam over Oxford clay, produced 12.8 t DM ha−1 using 354 kg N ha−1. Farm 3, on loam over clay with some free-draining sandy loam, produced 9.5 t DM ha−1 with 169 kg N ha−1. Allowing for differences in soil nitrogen and summer rainfall these yields were quite close to expectations from experimental results, despite considerable poaching damage to some swards.Animal production records were kept and all livestock were weighed every 6 months, so as to calculate the utilized metabolizable energy (UME) output from grassland on each farm. These were 73, 62 and 59 GJ ha−1 on farms 1, 2 and 3 respectively.The UME output, expressed as a proportion of the yield of grass (converted to ME), was taken as the ‘efficiency of utilization’ of grass. This was 70% on farm 1,44% on farm 2 and 58% on farm 3. Results are presented separately for the summer and winter 6 months. Utilization by grazing appeared to be markedly reduced under very wet soil conditions. Major losses occurred in the conservation of grass.The results begin to establish a valuable base of data from commercial farms.

Notes: The establishment and growth of white clover (cv. Grasslands Huia) was investigated in two experiments with seeds sown in slots cut in old permanent pasture in early May 1982. In the first experiment, fully protected from pests and diseases, either 50 (low), 100 (medium) or 200 (high) seeds were sown evenly per metre length of slot; 100 seeds were also sown in a clumped pattern unevenly. In the second experiment, the effects of various combinations of molluscicide, insecticide and fungicide were investigated.Nine weeks after sowing there were 27, 37, 34 and 62 seedlings m−1 on the low, medium even, medium uneven and high seed-rate treatments, representing percentage establishments of 54, 37, 34 and 31 respectively. By September, however, there were no significant differences between the dry weights m−l of any of the treatments. In the second experiment, damage by slugs occurred early in treatments without molluscicide and caused a mean reduction of 70% in seedling establishment. As a result significantly smaller dry weights m−1 of clover were recorded in September in these than in other treatments. Failure to protect against insects resulted in smaller seedlings.The results indicated that variations in seed rate within the range investigated were not critical but that failure to protect against pests decreased percentage establishment and subsequent growth.

Notes: Masham ewes were stocked at 12(L), 16(M) or 20 (H) ewes ha−1, with 1·8 lambs per ewe on average, on two blocks (A and B) in both 1976 and 1977. Silage was offered at pasture to the ewes in early lactation and as stocking rate increased from 12 to 16 and then 20 ewes ha−1, 47, 73 and 100 kg DM per ewe was consumed respectively. Silage was cut mainly in the autumn and 179, 100 and 9 kg per ewe was made for L, M and H respectively. More variation between treatments was measured in net herbage accumulation in summer than spring and 11·4, 10·0 and 9·7t DM ha −1 was grown over the whole season for L, M and H respectively as a mean of both years.High intakes of herbage OM were measured with values up to 40 g per kg live weight for BL lambs in 1976. Rapid lamb growth was achieved. There were effects of stocking rate on lamb performance and for L, M and H the respective growth rates to 4 weeks were 266, 248 and 247 g d−1; growth rates to sale were 274, 263 and 252 g d−1; days to sale were 124, 126 and 129; mean carcass weights were 17·7, 17·1 and 16·7 kg. Most ewes lost weight in early lactation, especially AH and BH in 1976, but there were no stocking rate effects on ewe weight change in 1977.

Notes: Classical methods of interpretation of reflection seismic data are such that interpretation and processing usually occur in the “collected” frame of reference. However, in recent times other data planes have gained increasing acceptance in seismology as a viable alternative. Through linear transformations applied to a record section, both the t—p- and p—x-planes can be produced. The r—p-domain may be obtained from the t—x-plane by a transformation known as slant stacking. Normal practice has been to do most of the data processing in the t—x-plane and then transforming to the r—p-plane. However, many of the procedures used in the t—x-domain can be modified for use in the t—p-plane to increase the coherence.Velocity inversion may be carried out either in the r—p-domain or further transformed to the p—x-plane where the modified Herglotz-Wiechert inversion may be applied. To perform the inversion, the t—p-wavefield is converted to a p—x-representation by the use of a new linear transformation technique, the cross-stack. By a simple sampling process along a particular p—x-trajectory, the Herglotz-Wiechert method can be used to reconstruct an acceptable velocity model of the subsurface. A comparison of derived velocity structures is made between that produced by the Herglotz-Wiechert technique and that of the Dix method.

Notes: The computation of static corrections requires information about subsurface velocities. This information can be obtained by different methods: surface wave analysis, short refraction lines, downhole times, uphole times and first arrivals from seismograms.For pure shear waves generated by SH sources the analysis of first arrivals from seismograms combined, if necessary, with short refraction lines has proved to be most accurate and economic.A comparison of first-arrival plots from P- and S-wave surveys of the same line measured in areas of unconsolidated sediments in northern Germany illustrates the characteristic differences between the two velocity models. P-waves show a marked velocity increase at the water table from about 600 to 1800 m/s. S-wave velocities of the same strata increase gradually from about 100 to 400 m/s. As a consequence, S-wave models are vertically and laterally more complex and, in general, show no significant velocity increase at a defined boundary as P-wave models do. Therefore, other suitable correction levels with specific velocities must be chosen.A comparison of “tgd-corrections” (correction time between geophone position and datum level) for P- and S-waves in areas of unconsolidated sediments shows that their ratio is different from the P-/S-velocity ratio for the respective correction level because of the greater depth of the S-wave refractor. Therefore, P- and S-waves are influenced by different near-surface anomalies, and time corrections calculated for both wave types are largely independent.

Notes: Forward filters to transform the apparent resistivity function over a layered half-space into the resistivity transform have been derived for a number of sample intervals. The filters have no apparent Gibbs' oscillations and hence require no phase shift. In addition, the end points of the filter were modified to compensate for truncation. The filters were tested on simulated ascending and descending two-layer cases. As expected, “dense” filters with sample spacing of In (10)/6 or smaller performed very well. However, even “sparse” filters with spacing of In (10)/2 and a total of nine coefficients have peak errors of less than 5% for p1:p2 ratios of 10–6 to 106. If a peak error of 5.5% is acceptable, then an even sparser filter with only seven coefficients at a spacing of 3 In (10)/5 may be used.

Notes: The concept of minimum phase is clarified for geophysicists by collecting in one place the properties of minimum phase functions. The “earliest energy arrival” property in the time domain, the “minimum phase-slope property” in the frequency domain, and of the role of causal all-pass filters are demonstrated. The emphasis is placed on keeping the mathematics within the realm familiar to geophysicists and on making clear the somewhat imperfect match between physical continuous time functions and their associated discrete time representations.

Notes: Volcanic maar structures sometimes contain important groundwater reserves. The volcanic formations of the Beaunit maar have a much higher magnetization than their granite substrata. The structure of the maar and its volcanoclastic infill can therefore be defined by magnetic surveying and resistivity soundings. Resistivities of the geological formations present in the Beaunit maar are comparable to those recorded in other volcanic systems. The geophysical data, confirmed by drilling tests, contain useful information on the hydrogeological regime of the volcanic system under investigation.

Notes: In migration procedures, the velocity profile of the subsurface is the most important input information. Since, in general, this information is only approximately known, errors in the migration output due to errors in the velocity input occur in all practical applications. In migration, velocity errors and depth errors can be interchanged. This interchange property is perfect in the paraxial approximation. From this result it follows that migration with incorrect velocities may still yield correctly migrated data if the imaging principle is modified. This attractive property can be used in the stripping version of migration (recursive migration).

Notes: The presence of the water layer in marine seismic prospecting provides an effective waveguide for acoustic energy trapped between the sea-bed and the sea-surface. This energy persists to large ranges and can be the dominant early feature on far-offset traces. On airgun records, there is commonly a lower frequency set of arrivals following the water-trapped waves. These arrivals are not as obvious with higher frequency watergun sources. By using a combination of intercept-time/slowness (τ—p) mapping on observational data and theoretical modelling, we are able to identify the origin of the events. If a very rapid increase in a seismic wavespeed occurs beneath the sea-bed sediments, a new waveguide is formed bounded by the sea surface and this transition zone. The low frequency waves are principally guided within this thicker waveguide. Numerical filtering in the τ—p domain followed by trace reconstruction is very effective in removing the low frequency noise.

Notes: The efficacy of the magnetotelluric and audiomagnetotelluric (MT/AMT) methods for detailing the structure of a hypothetical geological section is investigated by using the singular value decomposition (SVD) technique. The section is representative of southeastern Turkey, which is mostly covered by basalt and is a prime area for oil exploration. One of the geological units, the Germav shale at a depth of 600 m, is a problem layer for electromagnetic surveys because of its very low resistivity (on average 3 Ωm) and highly variable thickness across the area (200–900 m). In the MT frequency range (0.0004–40 Hz) its total conductance—or, since its resistivity is known from resistivity log information, its thickness—is the best resolved model parameter. The total depth to the Germav shale and the resistivity of the Cambrian/Precambrian basement are the marginally resolved parameters. In the AMT frequency range (4–10000 Hz) the resistivity of the surface basalt layer strongly affects the resolution of the other, less important, model parameters which are the total depth to the Germav shale and the total conductance of the Germav shale. The errors in the measurements determine the number of model parameters resolvable, and are also important for interpretation of the geological model parameters to within a desired accuracy.It is shown that statistical evaluation of the MT and/or AMT interpretations by using an SVD factorization of the sensitivity matrix can be helpful to define the importance of some particular stage of the interpretation, and also provides a priori knowledge to plan a proposed survey. Arrangements of MT and AMT observations, together with some Schlumberger resistivity soundings, on a large grid will certainly provide three-dimensional detailed information of the deep geoelectric structure of the area.

Notes: A method for interval velocity analysis is formulated on the basis of wavefield extrapolation, i.e., on the basis of wave-equation migration. When this scheme is applied to multioffset seismic sections or to an ensemble of CMP gathers, it allows for the proper treatment of dipping events. The underlying assumptions are that local velocities should be derived from data associated with events within the interval under consideration. To minimize the effect of the region above the layer of interest, the data are first extrapolated to the top of the analysis interval. Subsequent analysis of these data then pertains to the events within this interval. Velocity estimation consists of repeated wavefield extrapolations through the analysis interval using a set of trial velocities. The optimal velocity is chosen on the basis of coherency measures designed to express the collective phase agreement among a set of offset Fourier modes. The reliability of this approach to interval velocity estimation is demonstrated on synthetic multi-offset data.

Notes: A statistical treatment of MT data from the Pannonian Basin, Hungary, based on the distortion theory of the S-effect suggests that the local increase in crustal conductivity is connected with deep fractures.Field model measurements have recently been carried out in a shallow, quasi two-dimensional basin of well known tectonics using MT station distances of 1–3 km. The results of the investigation and numerical modelling of near surface distortions, support the idea that deep fractures (faults) contain the conducting formations. According to 2-D model computations, several conducting dykes at considerable depths can cause anomalies which are indiscernible from an anomaly due to a conducting layer. The significance of these results is discussed from the point of view of geothermal energy exploitation.

Notes: A principal use of the vertical seismic profile (VSP) is to determine the variation of seismic velocity with depth. Presented here is a discussion of the errors involved with the time picks of a VSP survey and several methods currently used to calculate a velocity section from these time picks. Another technique is proposed, based on the least-squares inversion of the traveltimes, to arrive at a better estimate and statistical description of the velocity section. This technique uses the Levenberg-Marquardt damped least-squares formulation and ray tracing through a horizontally-layered medium to iteratively refine the velocity section. The accuracy and robustness of the procedure are investigated by inverting noisy traveltime curves and comparing these results to the original model velocity section. Agreement is found to be good. One interesting feature of the inverse procedure is that, for certain geometries, it can resolve a few velocity layers, even though there are no measurements made inside those layers. Three actual VSP surveys are analyzed and compared to their corresponding sonic logs. In two of the surveys, it appears that there is some velocity dispersion. Velocity changes associated with gas saturation are evident on one of the surveys.

Notes: It is well recognized that in order to realize the full potential of the Vibroseis technique, one needs to ensure accurate phase locking and a meaningful cross-correlation. To achieve these two important objectives we require an accurate estimate of the compressional stress wave radiated by the vibrator into the ground.In this paper a simple method (subject of a patent application) is developed for predicting the compressional stress waves radiated by a vertical vibrator. The main feature of the proposed method is that it involves the field measurement of the acceleration of the reaction mass and the baseplate, respectively.The method is illustrated by computing the compressional stress waves generated by a typical vertical vibrator radiating into ice, chalk, sand, and mud. It is shown that for a seismic vibrator radiating into hard ground the pressure of the downgoing P-wave is 180° out of phase with the baseplate velocity. It is also shown that when the driving force of the seismic vibrator has a flat amplitude spectrum, the amplitude spectrum of the downgoing P-wave falls off by 6 dB/octave towards low frequencies.

Notes: Bedding on a scale small in comparison to the wavelength results in transverse isotropy. On the other hand, anisotropy has been observed in the field, particularly with horizontally polarized shear waves. In this article velocity data from 23 wells are used to estimate the resulting transverse isotropy and to compare these data with anisotropy observations reported in the literature. Since the well data were compressional velocity (or rather transit time) only, the corresponding shear velocity data were estimated on the basis of a reasonable—though arbitrary—assumption of Poisson's ratio. Actual densities were used wherever they were available; for the remainder of the data, density was assumed to be constant throughout the averaging interval.The anisotropy data estimated on this basis were generally considerably smaller than the observational data that have been reported in the literature. Although both the above assumptions—functional dependence of Poisson's ratio on transit time and constant density—tend to underestimate anisotropy, and although 23 wells is but a small subset of all available data, it appears unlikely that bedding alone could cause anisotropy of the magnitude reported.

Notes: In land seismic surveys spectrum equalization can increase the quality of seismic data in a selected frequency band. The power of lower frequencies in the spectrum of input traces is generally greater than that of higher frequencies, particularly in land seismic surveys because of ground roll. In order to improve the quality of seismic data it is necessary to raise the energy of higher frequencies to the same level as that of lower frequencies, without alteration of the phases.The first step of the method is to compute the amplitude spectrum of each input trace to determine a weighting function which is then applied to the amplitude spectrum in order to balance it. The function is the inverse of the short wavelength variation of the amplitude spectrum. The short wavelength variation can be obtained by interpolation between average values of the modulus of the amplitude spectrum computed in narrow bands within a selected band of frequencies. Another way of obtaining the short wavelength variation is to apply a low-pass filter to the amplitude spectrum. The calculations are readily performed in the frequency domain by the Fourier transform.Spectrum equalization is automatically adjusted to each trace and does not modify the average amplitude in the time domain. However, as the frequency band and energy of the ground roll both vary according to the distance from the shot, spectrum equalization tends to make the spectrum of output traces independent of the offset distance.The use of spectrum equalization before any two-dimensional filtering improves ground roll elimination. Continuity and resolution of horizons are also increased by spectrum equalization before CDP stack.Several examples of applications of spectrum equalization to seismic land and marine surveys are shown.