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: 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: 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: 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: 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: 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: 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 interpretation of Turam data is treated by means of a development of the conventional approach to this type of data. The development consists of a curve matching technique based upon a line current approximation to the true induced current pattern.The technique requires computation facilities which are readily available in modern programmable calculators.Practical examples of the use of the method illustrate that it can reveal information concerning the geometry of a conducting body which would not be gained by the conventional approach.

Notes: “Remote sensing” techniques have been used for many years as an exploration tool; broadly speaking, any method using sensors not directly connected to the earth's surface falls into this group.It is preferable, however, to use the term solely with reference to a number of new sensors which have recently become available to earth scientists for the detection of various radiations of the electromagnetic spectrum in the wave length range from 0.1 micron (UV) to some centimeters including the visible and infrared regions.Radiations artificially induced in the earth (active sensing) or spontaneous radiations (passive sensing) can be detected. Some of the bands investigated seem promising for application to geological surveys, especially using sensors on aircraft and satellites, which allow synoptic and large scale investigations.A brief account is given of the results already achieved and of future developments such as the ERTS and manned satellite (post Apollo) “earth observation” programmes.The potential role of Europe in the field is discussed.Details are given of some applications of thermal I.R. imagery to geological and vulcanological problems, with special emphasis on the research carried out in Italy.

Notes: In mathematical statistical filtering the deconvolution problem can be solved by two different methods:〈list xml:id="l1" style="custom"〉1by inverse filtering2by calculating the prediction error.Both methods are well known in the theory of Wiener filters.If, however, the generating process of the signal is known and can be described by a set of linear first order differential equations, then the Kalman filter can also be used to solve the deconvolution problem. In the case of the inverse filtering method this was shown by Bayless and Brigham (1970). But, while their method can only be used if the original signal is a colored random process, this paper shows that in the case of a white process the prediction error filtering method is a more appropriate approach. The method is extremely efficient and simple. This can be demonstrated by an example which maybe of special interest for seismic exploration.

Notes: The use of the electrical image method in calculating the potential of a source in the presence of a buried conductive sphere proves to be difficult because of the great number of terms making up the potential formula. The disadvantages contained by this method are avoided here by resorting to the theory of graphs. In view of this aim the “propagation” of the electrical images due to the sphere has been represented by an infinite graph to which an infinite matrix has been attached.Taking this as a basis several FORTRAN IV programs for the computation of resistivity pseudo-sections for various electrode configurations have been coded, and a great deal of such pseudo-sections have already been calculated with an IBM-360 computer. Using these results, recommendations regarding both the methodology of the geoelectrical prospecting of spherical shaped conductors and the most suitable way of building resistivity pseudo-sections are developed.

Notes: An account is given of the use of the square array technique in investigating the surface effects of rotational anisotropy when the axes of anisotropy are inclined to the surface. It is shown that, as with other arrays, two anisotropy parameters 〈inlineGraphic alt="inline image" href="urn:x-wiley:00168025:GPR249:GPR_249_mu1" location="equation/GPR_249_mu1.gif"/〉 and n can be derived by varying the array orientation.On the basis of these considerations, the effects of such anisotropy on the values of the mean apparent resistivity and azimuthal inhomogeneity ratio normally obtained in square array measurements is reviewed.Particular attention is paid to the variation of resistivity with orientation and it is noted that, in areas of moderate anisotropy, this variation is lower for the square than for the Schlumberger array.In addition to this advantage, the azimuthal inhomogeneity ratio obtained from square array measurements may be used to indicate the severity of anisotropy in an area and two field examples of this use are given.Where anisotropy is severe, gross variations of apparent resistivity with orientation are obtained with either square or collinear arrays. In these circumstances, the use of crossed measurements is considered and the particular stability of the crossed square array demonstrated.

Notes: Electrical resistivity sounding was attempted with success in 1960 on the surface of the Nile water in Khor Kundi El-Bahari, the starting place for the construction of the offset channel of the High Dam. The method was applied to determine the depth to the upper surface of the granite below the alluvium and Nile water. The success of the method conducted at two stations in the Khor is attributed to the favourable geological conditions of the section which consists of nearly homogeneous layers besides the good earthing conditions of electrodes dipped in water. The high resistivity contrast between water and granite facilitates interpretation of sounding curves. Factors interfering with the results of interpretation are mentioned.

Notes: This paper presents a dot chart for calculating gravity effects from two-dimensional bodies with the addition of a supplementary set of curves which can be used to determine the “end corrections” to modify the calculated effects when the length of the body is not relatively long with respect to the horizontal dimensions of the cross-section. The use of the chart is demonstrated by different examples. The results of theoretical cases amenable to mathematical treatment and computer results both compare favorably to those computed using the present chart.A fundamental incentive of constructing this chart, in a final form ready to use, is to make it available to gravity interpreters for utilization. Future effort in chart construction is thereby avoided.

Notes: By applying dynamic corrections a seismic trace recorded at a distance x from the energy source should be varied in such a way as to obtain a trace which would be recorded at zero-distance, i.e. at the source itself. Only such a zero-offset-trace contains the correct sequence of reflection coefficients (reflectivity function), whilst all other traces contain a distorted reflectivity function. In the simplest case, the reflectivity function is compressed over a shorter time whereas in more complicated cases a partial inversion of the reflectivity function results. This happens when some of the reflection hyperbolae intersect one another.The reconstruction of the true zero-offset reflectivity function by the application of dynamic corrections can only be an approximative process. In the first case mentioned we must expect a decrease in accuracy of the corrected trace in comparison with a zero-offset-trace. In the second case, where intersections of the hyperbolae occur, accurate reconstruction is clearly impossible.The problems are discussed with the help of theoretical and practical examples.

Notes: Fourier transform techniques have been used to calculate the theoretical filter (amplitude) response function of Nth order vertical derivative continuation operation. The amplitude response functions of the vertical gradient and its continuation follow from the same. These response functions are subsequently used to calculate the weighting coefficients suitable for two dimensional equispaced data. A shortening operator has been incorporated to limit the extent of the operator. For comparative study, some of the developed coefficient sets and the one presented in this paper are analysed in the frequency domain and their merits and demerits are discussed.

Notes: In this paper it is shown how one may obtain a generalized Ohm's law which relates the induced polarization electric field to the steady-state current density through the introduction of a fictitious resistivity defined as the product of the chargeability and the resistivity of a given medium.The potential generated by the induced polarization is calculated at any point in a layered earth by the same procedure as used for calculating the potential due to a point source of direct current.On the basis of the definition of the apparent chargeability ma, the expressions of ma for different stratigraphie situations are obtained, provided the IP measurements are carried out on surface with an appropriate AMNB array. These expressions may be used to plot master curves for IP vertical soundings.Finally some field experiments over sedimentary formations and the quantitative interpretation procedure are reported.

Notes: In the linear filter method of interpreting resistivity sounding data, as developed by Ghosh (1971), it appears that the filter function in the x-domain approaches an oscillating function for both large positive and large negative abscissa values. In the present note the reason for this oscillating behaviour is derived, and a possible practical application is indicated.

Notes: Amplitude spectra of input FM signals used in the vibratory source method of seismic exploration often show undesirable oscillations near the initial and terminal frequencies. These oscillations have an effect on the correlation background and distort the output signal. Considerable improvement in reducing the amplitude of these oscillations is obtained using a proper taper fuction. Attention is given to the relation between the tapering time and bandwidth of the spectrum.Analyses of the spectra of the received data from vibratory sources show considerable attenuation in comparison with the original field sweep. Since the matched filtering process will result in a series of waveforms which have the shape of the autocorrelation of the input signal, consideration is given to the autocorrelation function and its zero-lag coefficient of the FM signal in the presence of attenuation. A method has been developed which compensates for the attenuation and recovers the distortion of waveforms when the received data is correlated.The design of a waveform shaping filter for vibratory source data is given to reduce the influence of phase distortion on the received waveforms as well as to increase S/N ratio resolution. Parameters used for this filter are based on the properties of the FM signal and its autocorrelation function.Several examples from field data are presented to illustrate the methods. The results indicate that the use of the above techniques yields sections with good frequency resolution and improved S/N ratio.

Notes: The suggested interpretation technique is based on discrete linear filtering of VLF data. The output of the described filtering results is expressed in terms of an equivalent current density at a specified depth that would cause the measured magnetic field. The most practical six-point filter gives an accuracy of 8%. The filter is an extension of the commonly used Fraser filter to process VLF dip-angle data.Filtering the same data set for various depths gives an idea about the change of current density with depth. Areas with high current-density correspond to good conductors. The conductor dip can also be determined.The use of the technique is illustrated on theoretical and field examples. In all cases a good correlation with original models and other types of geophysical measurements was obtained. As shown in the last example, the filtering technique is also applicable in interpretation of other electromagnetic methods.

Notes: Using approximate boundary conditions, expressions for electromagnetic fields have been derived for a thin, highly resistive layer lying between two homogeneous layers excited by an electric dipole grounded on the surface of the earth. The variations of the fields with the parameter T/T1 (ratio of the transverse resistance of the thin layer to the transverse resistance of the first layer) were studied in relation to frequency, time, the normalized separation source—receiver, and the angle between the source and the radius to the observation point. For a value of h2/h1 (ratio of thickness of second layer to the thickness of the first layer) approximately equal to 0.2, the general three-layer medium case gives the same results as this approach. It was found that the electric fields have a very strong dependence on the parameter T (transverse resistance) which characterizes the thin, highly resistive layer. However, the magnetic fields depend only very weakly on this parameter.

Notes: Microgravity can be used for predicting rock bursts. For the first time gravity anomalies related to rock bursts have been recorded. The methodology developed has lead to the first successful predictions. The depth of the rock-burst focus might be determined on the basis of analytical downward continuation of related gravity anomalies: the focus is treated as a singular point of the gravity potential and its derivatives.The rock-burst gravity anomalies might be explained on the basis of the assumed dilatancy process that causes the rock burst. The mean density change of the rock mass threatened with rock bursts can be estimated from the corresponding rock-burst gravity anomalies.

Notes: A model has been set up for the interpretation of geoelectrical sounding data for certain kinds of clay deposits containing gradually varying amounts of sand. The model assumes that in the so-called gradient layer, resistivity varies linearly with depth. Model calculations show how such a layer can be replaced by two homogeneous layers.An inversion procedure using the Marquardt algorithm has been developed for the interpretation of sounding data obtained with the Schlumberger array; it assumes a gradient layer beneath several covering layers. The procedure is demonstrated on two clay deposits. A comparison is made between the newly developed interpretation, the traditional approach using model curves, and computerized inversion for a homogeneous layer model.

Notes: Gravity data have been transformed into a three-layer, three-dimensional model by using an automatic procedure based on linear filtering. The Bouguer anomaly is first transformed by linear filtering into density variations located between two planes 1100 and 2500 m deep. These densities are then transformed into thicknesses with a constant density contrast of 0.4 g/cm3 with two geological constraints for the second and third interface:—minimum at 2500 m depth;—maximum below a variable limit given by geology.This gives the contact between the second and third layer. Differences between measured and computed gravity are then applied by a similar procedure to a layer located between depths of 0 and 500 m, giving the contact between the first and second layer. Interesting secondary anticlines and transverse faults are shown by various structural maps.

Notes: Migration of an observed zero-offset wavefield can be performed as the solution of a boundary value problem in which the data are extrapolated backward in time. This concept is implemented through a finite-difference solution of the two-dimensional acoustic wave equation. All depths are imaged simultaneously at time 0 (the imaging condition), and all dips (right up to vertical) are correctly migrated. Numerical examples illustrate this technique in both constant and variable velocity media.

Notes: Synthetic seismograms can be very useful in aiding understanding of wave propagation through models of real media, verification of geologic models derived from interpretation of field seismic data, and understanding the nature and complexity of wave phenomena. If meaningful results are to be obtained from synthetic seismograms, the method of their computation must, in general, include three-dimensional geometrical spreading of wavefronts associated with highly concentrated (i.e., point) sources. The method should also adequately represent the seismic response of solid-layered media by including enough primaries, multiples, and converted phases to accurately approximate the total wavefield. In addition to these features, it is also very helpful, although not always essential, if the method of seismogram computation provides for explicit identification of wave type and ray path for each arrival. Various seismograms, computed via asymptotic ray theory and an automatic ray generation scheme, are presented for a highly simplified North Sea velocity structure. This is done to illustrate the importance of the above features and to demonstrate the inadequacy of the plane-wave synthesis method of seismogram computation for point sources and the limitations of acoustic models of solid-layered media.

Notes: A periodically stratified elastic medium can be replaced by an equivalent homogeneous transverse isotropic medium in the long wavelength limit. The case of a homogeneous medium with equally spaced parallel interfaces along which there is imperfect bonding is a special instance of such a medium. Slowness surfaces are derived for all plane wave modes through the equivalent medium and reflection coefficients for a half-space of such a medium are found. The slowness surface for the SH mode is an ellipsoid. The exact solution for the reflection of SH-waves from a half-space with parallel slip interfaces is found following the matrix method of K. Gilbert applied to elastic waves. Explicit results are derived and in the long wavelength limit, shown to approach the results for waves in the equivalent homogeneous medium. Under certain conditions, a half-space of a medium with parallel slip interfaces has a reflection coefficient independent of the angle of incidence and thus acts like an acoustic reducing mirror. The method for the reflection of P- and SV-waves is fully outlined, and reflection coefficients are shown for a particular example. The solution requires finding the eigenvalues of a 4 × 4 transfer matrix, each eigenvalue being associated with a particular wave. At higher frequencies, unexpected eigenvalues are found corresponding to refracted waves for which shear and compressional parameters are completely coupled. The two eigenvalues corresponding to the transmitted wavefield give amplitude decay perpendicular to the stratification along with up- and downgoing phase propagation in some other direction.Much of this work was performed while the author was at the Department of Geophysics and Planetary Sciences, Tel-Aviv University, Ramat-Aviv, Israel. The author is grateful for illuminating discussions with K. Helbig and K. Gilbert.

Notes: Observations made from space have become almost a standard method of geologic investigation. However, interpretation of space photographs is not always unambiguous. Particularly, the nature and spatial location of various objects photographed from space remain obscure. It is especially difficult to investigate territories where basement rocks are overlain by sedimentary cover.From the example of investigations carried out in the central part of the Russian Platform it is shown how separate lineaments and ovoids reflect the peculiar features of geologic, hydrogeologic and geophysical structure of the region. It has been established that—within the sedimentary cover—they correspond usually to zones of increased fissuredness and are associated with low-amplitude tectonic dislocations.The authors concentrate on the elucidation of the transfer of deep geologic information to the surface. The sedimentary cover is regarded as a communication channel which makes it possible to employ the concepts of information theory for assessing its “transparency”. The contrast due to heterogeneities of the sedimentary cover in space photographs and, consequently, the degree of their manifestation in geologic and physical fields, depends on the water content of the rocks composing the section. Along with hydrogeologic characteristics, geo-electric parameters, such as specific electric resistivity, longitudinal conductivity and transverse resistivity, can serve as criteria of the degree of rock moistening.When the location of lineaments does not coincide with the plane axes of geophysical anomalies, an inclination of the fault plane can be supposed. The angle of inclination is characterized by the deviation of the lineament from the axis of the anomaly.The data resulting from the investigation of selected standard structures were used to develop a classification based on the degree of their expression in hydrodynamic, geoelectric, gravitational, and magnetic fields and depending on the spatial position and peculiarities of geologic arrangement. For geophysical investigations of unknown structures detected by space observations, this classification can be used as the basis for solving inverse problems by the methods of mathematical statistics.

Notes: How do papers on electrical and electromagnetic methods published in Geophysics and Geophysical Prospecting reflect the actual field application of the various techniques? Our study based on 10 volumes of the journals (1973–1982) suggests a rather good agreement. EM techniques are the most important both in the field and in print. However, an imbalance was detected in the subject matter. Too many published papers deal with theory, and too few with new instrumentation and techniques, case histories, and physical properties of rocks. Papers presented at annual meetings of the EAEG and SEG give a more even-handed and accurate view of geophysical activities. Unfortunately, many of our colleagues around the world are not able to attend the annual meetings regularly and depend on journals for their information. Practicing and research geophysicists should submit more papers on subjects that at present are not adequately covered.

Notes: Gravity and magnetic anomalies may have the same source but it is always difficult to achieve correlated solutions if interpretations are carried out separately. Therefore it is useful to invert both anomalous gravity and magnetic profiles at the same time, so that the solution mav be more constrained.Existing inversion techniques do not lend themselves equally to this kind of extension, since the parameters—such as density and susceptibility contrasts—should not be related. The algorithms more easily adaptable to simultaneous inversion are those which look for the shape and the position of anomalous bodies with constant unknown density or susceptibility contrasts.In this study, we use the generalized inverse method in the 2 1/2 D case. Parameters are the coordinates of the vertices of the polygonal cross-sections of the anomalous bodies and both density and susceptibility contrasts of each body. The two types of profile to be processed must be superposable but may have different lengths, spacings, and coordinate origins. Both synthetic examples and field data from geologically known areas have been processed, and it appears that the simultaneous inversion technique may provide an important additional degree of control in the interactive interpretation process.

Notes: A model of the seismic trace is generally given as a convolution between the propagating wavelet and the reflectivity series of the earth and normally it is assumed that a white noise is added to the trace. The knowledge of the propagating wavelet is the basic point to estimate the reflectivity series from the seismic trace.In this paper a statistical method of wavelet extraction from several seismic traces, assuming the wavelet to be unique, is discussed. This method allows one to obtain the propagating wavelet without any classical limitative assumptions on the phase spectrum. Furthermore, a phase unwrapping method is suggested and some statistical properties of the phase spectrum of the reflectivity traces are examined.

Notes: An interactive least-squares method for the interpretation of VES curves was proposed by Johansen (1977). The method permits one to select some parameters (thicknesses and/or resistivities of individual layers) and to change the rest in such a way that the interpreted model approaches the measured data. This note suggests a modification of Johansen's method, in which not only the individual parameters can be selected but also linear combinations of parameters—in particular, the sum of thicknesses of several layers.

Notes: A borehole magnetometric resistivity (MMR) experiment is described in which an attempt is made to determine the extent and orientation of zones of fissuring within an otherwise massive sequence of carboniferous limestone. The region under investigation lies beneath a landfill site and the main parameter of interest is the direction of flow of ground water, which will be influenced by the orientation of the faults or fissures within the rock. The MMR method possesses some extremely advantageous features for application to hydro-geological problems of this kind. The method is sensitive to electrical current channelling within an otherwise relatively resistive medium, and the detection of the magnetic field within the borehole does not depend upon electrical contact of the receiver with the ground. Consequently, the method can be used in dry or plastic-cased boreholes. A direction of maximum electrical current flow is deduced from the MMR data which coincides with the predominant direction of jointing within the region.We are most grateful to the Water Research Centre, Medmenham, Berkshire, for supporting this research and particularly to Dr M. Fleet and Mr K. Beesley for their help and advice.

Notes: The geophone spurious frequency is modeled as the resonance of the planar motion of a spider spring carrying a moving mass. An analytic solution is found using the Castigliano method by assuming that the spring arm is a single-mode vibrating cantilever beam. The spring shape is found from this analysis. When the typical spring has a circumferentially varying cross-section, the Castigliano method is no longer applicable. A dimensional analysis is used as an approximate method for general design. Based on the theoretical result, a rotational fixture and a translational fixture were designed for experimental purposes. A low-noise fixture and a phase-averaging technique provide the amplitude and spurious resonance in the frequency domain. Finally, a test is run by using a rotational fixture to compare with the approximate method of spurious frequency prediction. Very good agreement between prediction and experiment is found.

Notes: A finite-difference solution and an integral algorithm are developed for computing time-domain electromagnetic fields generated by an arbitrary source located in horizontally stratified earth.The finite-difference problem is first solved for the kernel of an integral Bessel transform of the desired field and then an inverse transformation is performed either by numerical integration or by linear filtering. A description of the appropriate algorithm as well as numerical results are presented for two cases: an axially symmetric source (a vertical magnetic dipole or a circular loop) placed on the surface of a horizontally layered section overlying either an insulating or an ideally conducting basement.In comparison with the conventional methods used for computing transient electromagnetic fields, the integral-finite-difference method is more universally applicable.

Notes: Temperature modeling around the Vejrum salt structure in Denmark shows that the heat flow near the top is twice the regional heat flow. The temperature in the vicinity of the top of the diapir is thus anomalously high. Overlying permeable formations could therefore be attractive for geothermal energy exploitation. A negative temperature anomaly of −20°C is calculated around the root of the salt body. Comparisons between calculated and measured temperatures are used to test the model. It is shown that measurements of heat flow (or temperature gradients) in wells penetrating the region of groundwater circulation could be used in identification of salt diapirs.

Notes: An electromagnetic sounding experiment with a large square loop as source was carried out on the dried-up bed of a water reservoir near the town of Dharmavaram in Andhra Pradesh. The sounding was performed in both geometric and parametric modes, and involved measurements of phase as well as amplitude of the vertical magnetic field inside and outside the loop. The six-frequency EM system used for the experiment was found to be a workable system for electromagnetic sounding. The results of the experiment more or less confirm the conclusions from theoretical modeling. They also show that even though the earth is not always electrically horizontally layered over a sufficiently large lateral extent for the one-dimensional model to be strictly valid, it is still possible to apply such modeling to sounding curves taken one part at a time and obtain layer parameters which check qualitatively with the layer parameters obtained from direct current resistivity sounding.

Notes: Finite-offset seismic reflection modeling of acoustic waves, propagating in a two-dimensional depth section of arbitrary complexity, is discussed. The procedure developed employs the principles of simplified (far-field) diffractor theory and ray tracing. Each reflector is represented by a set of discrete secondary sources or diffractors and the wavefield associated with each diffractor is calculated directly in the time domain by ray tracing. Reflections and diffractions are subsequently built up by the numerical superposition of these wavefields. This superposition is nondispersive for all frequencies for which the Fresnel zones are large compared with the diffractor separation.All primary travel paths connecting the shot to diffractor and diffractor to geophone are accounted for together with phase changes induced by focal events. The method allows the modeling of arbitrary trace gathers for energy originating from selected reflectors. The nonsequential nature of the algorithm makes it suited to machines capable of carrying out many similar operations in parallel or concurrently. Diffractor theory also provides physical insight into wave scattering and focusing. In particular, the half-differential waveform associated with a line diffractor leads to an explanation of the 90° phase lead induced by a cylindrical focus and, similarly, the full differential waveform of a point diffractor can be used to explain the 180° phase shift induced by a point focus.

Notes: The addition of a third dimension to seismic data acquisition offers significant improvement of our reflection images of geologic structures, but potentially at a substantial increase in data processing cost. Here we study an approach to migration of three-dimensional (3-D) stacked data that can considerably reduce the burden of computation and data manipulation. Cascading two-dimensional (2-D) migrations in orthogonal lateral directions, we obtain an efficient 3-D scheme that is exact for homogeneous media. As applied to examples of synthetic data from inhomogeneous media, the scheme introduces errors well below those attributable, in practice, to uncertainties in migration velocity. Application of this two-step approach, using both Kirchhoff-summation and finite-difference algorithms, to a 3-D seismic land survey over an area of known hydrocarbon traps additionally proves the importance of 3-D migration in stratigraphic imaging.

Notes: Downhole periodic seismic sources have been proposed as a means for studying fluid zones in hydrothermal/magma systems encountered in deep scientific holes drilled as part of the Continental Scientific Drilling Program. Scientific measurements in these holes will require the use of downhole, nondestructive, high resolution, seismic tools. Downhole periodic seismic sources tested in a shallow geothermal zone showed that a pneumatic oscillator could effectively generate downhole low-frequency (10–100 Hz) seismic waves and that these waves could effectively penetrate a highly fractured formation. An improved version of these downhole sources, capable of swept-frequency operation, is currently under development. This type of downhole seismic source has applications in crustal studies and in the evaluation of petroleum, groundwater and geothermal reservoirs.

Notes: Numerical investigations have been carried out on suitable interpretation schemes for gravity, based upon Backus-Gilbert inverse theory. A two-dimensional linear model consisting of horizontal prisms (grid shape) has been adopted for interpretation of gravity profiles. An inversion technique, which uses a linear approach suitable for rapid computation, and studies with synthetic bodies have led to improvements in the method, mainly by applying weighting factors, which improve the solution in the final stages of computational work. Gravity data from the well-known Cyprus island positive anomaly were inverted to demonstrate the feasibility of the method.

Notes: The problem of a plane wave incident on a non-isotropic dipping layer lying over an isotropic conducting substratum has been studied and some numerical results are presented to show the effects of〈list xml:id="l1" style="custom"〉1) degree of anisotropy m,2) conductivity contrast between the substratum to the upper layer b,3) angle of inclination of the axis of anisotropy α,on the variation of amplitude and phase difference of the wave impedance with the normalized skin-depth (d/δ1).

Notes: Consideration of the spectral representation of gravity and magnetic fields shows the field to be the result of the convolution of factors depending on the parameters of the field. Wiener filters, calculated using model transformations, provide an optimum method for altering these factors and hence effecting field transformations.

Notes: A discussion is given of the requirements, the advantages and the methods to be considered in attempting the quantitative interpretation of gravity or magnetic fields from computed maps of the vertically derived field.The transform which is used here is the first vertical derivative (or vertical gradient) with or without downward continuation, but the computed maps are in fact obtained by a controllable Fourier method in which two kinds of operations can be simultaneously performed in complete independance: on one hand the separation of any part of the data by frequency cut-off, and on the other hand the transformation by vertical derivation or continuation of the part which is retained.Taking as raw data either actual surveys or artificially constructed maps, it is first shown how separated and transformed maps of this type can be efficiently obtained under quite flexible conditions, using a special computer program. It is further seen that for correctly controlled filterings the accuracy of the computed maps actually permits to take them as the basis for quantitative interpretation.To effect this, any one of the conventional methods which make use of equivalent model computations may in the first place be adapted to the interpretation of gradients, with the benefits, however, of an enhanced lateral separation of the anomalies and of a large attenuation of the regional effects.Particularly, the delineation of horizontal contours for even fairly complex models can often be made directly in a sufficiently safe way on the anomalies as they show on the vertical gradient maps. This greatly accelerates the process of determining equivalent model bodies.More special methods of interpretation can also be designed by taking into account first the fact that the vertical derivation of the field amounts to an operation of separation on the field's sources themselves, and in addition the availability of the frequency form of the information as a result of using a Fourier method of transform computation.Trial utilisations of various such interpretation processes of either the conventional or the less conventional type, are presented especially in connection with an interpretation study on transformed maps of parts of an offshore aeromagnetic survey (English Channel).

Notes: Generation of electrical power through the use of natural steam obtained from the Earth has the potential to be a significant source of pollution-free power. Successful future exploitation of geothermal power depends upon the development of cost-effective exploration techniques. A field experiment in the Imperial Valley of California indicates that there is a close empirical relationship between hot water deposits at depth and an anomalously high seismic background level at the surface. If this relationship proves to be a general one, it will serve as the basis of a simple and inexpensive method of geothermal exploration.

Notes: In massive rocks ground waters mainly flow in fracture zones whose permeability greatly changes depending on their filling. When studying ground-water flow in fissures, the results of observations of the electric fields of filtrational origin—which, in this case, considerably differ from those in porous media—can be used. Therefore the authors have made experiments on fissured media models. The measurements have been made in a special filtration tube with the fissured media simulated by a system of quartz glass plates. The spacings between plates were regarded as fissure widths. The observations have been made in fully open “fissures” and in those partially filled with sand or sandy-clay material.These experiments have resulted in establishing a dependence between the values of streaming potentials and pressure drops. The SP values have particularly been found to decrease with the opening of “fissures”. The most intensive electro-filtrational fields were observed at 40 per cent filling with medium grained sand. Additions of argillaceous material to the sand filler brought about sharp reductions in the intensity of electro-filtrational potentials.

Notes: During the last few years the airpulser, or air gun, has become very common as an energy source for marine seismic surveys. This paper describes the physical processes which take place during the operation of the pulser and develops theoretical results concerning the energy and frequency of the radiated signal and the amplitude decay of the secondary bubble pulses. The theory takes into account the presence of the airpulser itself which is assumed to be a rigid sphere within the bubble of released air. The theoretical results are combined and compared with measurements made of the pressure within the airpulser, the acceleration of the body of the pulser, and the amplitude and frequency of the signal radiated into the surrounding water. A formula for calculating the bubble frequency is given and a diagram made of the energy partition between mechanical losses, radiated energy, etc. Finally, a comparison is made of the energy release from the airpulser with that from TNT.

Notes: An equivalent stratum model is used to provide an explanation for the power spectrum characteristics of potential fields.The power spectrum of observed magnetic fields is found to consist of effects due to sources which can be represented by apparent monopoles and dipoles. Similarly, the power spectrum of observed gravity fields is found to consist of two groups of effects. A matched filter is proposed whereby the effects due to the two components in the potential field can be separated. As a consequence of the established theoretical expressions a scheme is suggested whereby the aliasing power of sampled data can be estimated. Also, the concepts of downward continuation, reduction to the pole, and reduction to pseudo-gravity of magnetic fields are re-examined in light of the theoretical expressions due to the equivalent stratum model.

Notes: The oscillatory character of the solution obtained by Qureshi and Mula's method for the direct determination of two-dimensional mass distributions from gravity anomalies is examined.

Notes: A number of resistivity profiles over dolerite and metadolerite dykes have been taken using Wenner and Schlumberger configurations around Dhanbad in order to find the potential of these for water accumulation. These dykes give rise to both resistivity highs and lows as compared to the background. The results are discussed in the light of known theories. Some suitable sites for water accumulation have been delineated.

Notes: Differential and difference soundings of resistivity for real, imaginary and complex values of the newly introduced parameters ν and μ are considered. The transverse and longitudinal differential soundings, known up to the present, correspond with the values ν= 1 and ν=— 1.The point of departure is the generalization of the concept of Dar Zarrouk parameters and the dimensional analogy between D.Z. resistivity and transverse and longitudinal differential resistivities. Properties of generalized D.Z. curves as well as the method of their construction are given. The examples of differential and difference curves are presented and the properties of these curves are discussed; particulary the effective spacings are determined. The problems of the realization of differential and difference soundings are considered; the bisymmetrical and asymmetrical arrays are proposed as a practical ones.

Notes: One of the methods for studying the near-surface low-velocity zone and for subsequent determination of static corrections is the technique of employing first arrivals. During the past few years several computer algorithms, based on some simplifying assumptions, have been suggested for automatic determination of first arrivals.This paper suggests a new method for automatic picking of first arrivals, which works under quite general assumptions concerning the character of the data. The method is based on the use of correlation properties of signals and the application of a statistical criterion for the estimation of first arrivals time. A number of dimensionless parameters is used in the algorithm making it possible to regulate the level of reliability and the resolution of the picking procedure.The second stage of the algorithm is the parameterization of the traveltime curve, that is a division of the previously obtained t—x curve into separate rectilinear segments. The suggested algorithm of parameterization is based on an heuristic use of some properties of maximum likelihood estimates. This permits location of the breakpoints of the t—x curve and the estimation of the parameters of each rectilinear segment.A computer program has been written based on the picking and parameterization algorithm. This program has been tested on a large amount of field data and the results show that it works at least as well as the hand procedure.

Notes: A theoretical study is presented of the effect of source/receiver patterns with fixed phase difference on the radiation/receiving characteristics of the single elements of the pattern in linear elastic media. The radiation characteristic of two vertical vibrators in phase and in counterphase is investigated in detail.By introducing source/receiver patterns the characteristics become dependent on dimension of the pattern, frequency and angle of observation in the horizontal plane. The types of waves generated by source patterns are, independent of phase differences, equal to those generated by single sources. With two vertical vibrators in phase radiation of P-waves can be improved relative to that of SV-waves. With two vibrators in counterphase radiation of SV-waves in vertical directions is enhanced relative to that of P-waves; this may be a valuable method of SV-shear wave generation in exploration seismics. SV-wave amplitudes larger than P-wave amplitudes can be expected at angles larger than 5°–10° with the vertical.

Notes: Three current migration methods derived from the wave-equation model and based on approximations by Fourier series, finite-differences, and a combination of both (mixed method) are reviewed. Each method is examined in the context of horizontally layered media, where the velocity of wave propagation varies with depth alone. It is shown how, under such conditions, the special properties of the linear equations obtained from finite-difference and mixed approximations can be exploited for computational efficiency. A numerical algorithm is described which substantially reduces the complexity and the computational cost of the solution. This algorithm is directly applicable to vector processing.The description of each method includes derivation, dispersion relation, stability properties, and applicable numerical algorithms. Comparative analyses show that, for the same computational cost, the finite-difference approximation has poor accuracy characteristics and should be used only when the boundary conditions are difficult to express in the frequency domain. The results are illustrated with numerical experiments on synthetic data.

Notes: Recently the author developed and demonstrated (Safar 1980) an efficient method for operating the airgun. The method involves the generation of a short seismic pulse from the pressure bubble pulses radiated by an airgun when fired several times at the same optimum depth but with different chamber pressures.The purpose of this paper is to present and discuss the test results obtained when implementing the same method using a two-dimensional airgun array. The array consists of seven 0.65 liter airguns fired simultaneously at the same depth but with different chamber pressures.It is shown that the far-field pressure pulse radiated by the seven 0.65 liter airgun array is similar to that radiated by the Flexichoc seismic source. It is concluded that the proposed airgun array can be used as a subarray to form an extremely powerful super-long array suitable for deep seismic exploration.The author would like to thank the Chairman and Board of Directors of the British Petroleum Co. Ltd for permission to publish this paper. Thanks are also due to Mike Symes and Lovell Cox for carrying out the field tests and Seismograph Service (England) Ltd for providing the airguns.

Notes: Recently, a heuristic and empirical analysis of a method of three-dimensional (3-D) migration has been presented by Gibson, Larner, and Levin. It consists of sets of two-dimensional (2-D) migrations in orthogonal directions. In this paper, we examine this two-pass method using scalar wave theory and show its action to be identical to that of so-called full 3-D migration for constant-velocity media. The demonstration of equivalence is carried out in both the frequency-space and frequency-wavenumber domains, the latter providing a particularly transparent view of the mechanism of the two-pass approach. We then extend our analysis to consider how various approximate algorithms used in practical implementations affect the validity of our conclusion.

Notes: Variations of seismic interval velocities within the cable length cause anomalies in the stacking velocity analyses. Utilizing the approximation of rectilinear ray propagation, i.e. supposing that the velocity changes cause time delays only, it is shown that the stacking velocity anomalies are linearly related to the interval velocity variations. In particular, the stacking velocity anomaly is calculated when the interval velocity of an intermediate layer undergoes a stepwise variation. The amplitude of the anomaly increases with the ratio between horizon depth and cable length.From the forward model, a program for the inversion is derived in order to identify lateral changes of interval velocities from unsmoothed stacking velocity analyses.Some examples of the application of this technique to synthetic and real data are presented.

Notes: A theoretical solution is presented to the problem where a VLF anomaly is generated by a conducting half-plane or a perfectly conducting wedge below a stratified overburden. The solution is obtained by the use of a scattering matrix for plane-wave eigenfunctions. VLF anomalies have been computed for different values of the conductance and dip of the half-plane. The phase of the VLF anomaly due to a conducting half-plane depends on the conductance and the distance to the half-plane. Close to the half-plane the tilt angle and ellipticity are of opposite sign for a perfect conductor, but the ellipticity will change sign for a poor conductor. The VLF anomaly for a perfectly conducting wedge is essentially determined by the position of the upper surface of the wedge, i.e. the anomaly will closely resemble the anomaly of a perfectly conducting half-plane in the same position as the upper surface of the wedge.