PARAXIAL RAY TRACING IN 3D INHOMOGENEOUS,ANISOTROPIC MEDIA1

Many investigations of the propagation elastic waves within the earth require a technique for producing synthetic seismograms which is capable of modelling 3D propagation effects. Ray methods are an excellent option for these problems, because they can be made fully 3D and allow a relatively quick and flexible computation of synthetic seismograms. However, the two point problem of finding the ray which connects exactly a specific source and receiver, may still be difficult and time consuming. Therefore, application of the paraxial method, which allows extrapolation of the information on a given ray to nearby receiver locations, is very valuable. With this approach, great savings in computation time and significant simplification of computer codes are possible. We investigate the application of the paraxial ray method to two problems in which the effects of 3D seismic wave propagation are important. The first is a model of a reef structure. In this case, we consider synthetic seismograms for a VSP experimental configuration. When the SV source and well are located along the axis of the reef, only 2D propagation effects are observed. If the source-receiver plane is located to the side of the reef, however, the 3D shape of the reef causes significant amplitudes to be predicted for shear-wave arrivals on the transverse component of the synthetic seismograms. The second example is a ID, layered earth model, but it includes two layers which are azimuthally anisotropic due to the presence of aligned, vertical fractures. This anisotropy leads to 3D raypaths. Synthetic seismographs are presented for a cross-hole geometry both for an equivalent isotropic model and for the direction parallel to the fractures and at an angle of 45° to the fractures in the anisotropic case. These synthetics show that the differences between the isotropic case and the case for source and receivers aligned with the plane of the fractures are small and subtle. On the other hand, the predictions for the direction at 45° to the cracks show shear-wave splitting and significant transverse component signal. These results have important implications for both modelling and for applications such as tomography. It is clear that for some of the cases considered, a 2D algorithm will lead to errors in interpretation of data. In addition, ray-based tomographic techniques will have great difficulty in obtaining a well-defined 2D planar image when the signals are propagating in regions outside the image plane.     

Rank exposure index on a Silurian carbonate tidal flat

Adequate discrimination between ancient intertidal and supratidal limestones and dolomites is often difficult because of vagueness of terms, problems of resolution, and the fact that most ‘diagnostic’sedimentary features are not restricted to narrow environmental zones. Furthermore, the rise and fall of sea-level in restricted, epeiric seas is usually due more to storm surges than to true tides, and hence the occurrence and extent of tidal-flat inundations are quite irregular. ‘Rank exposure index’is a concept introduced (1) to quantify the relative percentage of time that sediments are subaerially exposed and (2) to produce a more refined subdivision of shoreline carbonate rocks. The scheme employed to determine this index is a combination of numerical methods: cluster analysis separates 46 Tonoloway Limestone samples (Silurian of Maryland) into groups based on 16 sedimentary, mineralogic, and organic features, and minimum spanning tree produces a connected series of the five resultant groups. These groups are then ranked in sequence from 1 (greatest subaerial exposure) to 5 (least). Rocks of rank 1 are thought to have been deposited at the highest position within the tidal zone, as on a beach ridge. Rocks of decreasing rank reflect greater frequency of flooding. Those of rank 5 represent sediments deposited closest to low-water level, as in a tidal pond. Rock type, mineralogy, fossil content, algal stromatolites, and sedimentary structures are closely related to exposure rank as interpreted for the Tonoloway carbonates. Rank exposure index, therefore, is a semi-quantitative measure to aid in the detailed reconstruction of ancient shoreline environments.     

Modern sediments and sedimentary processes in Lake Rudolf (Lake Turkana) eastern Rift Valley, Kenya

Tertiary sediments around Lake Rudolf (now Lake Turkana) in the East African Rift Valley have yielded abundant palaeontological and palaeoanthropological remains. The present study provides a basis for interpreting the ancient lake environment and furthering our knowledge of rift valley lacustrine deposits. Bottom sediments in Lake Rudolf are fine-grained (average 71% clay) well laminated and have montmorillonite, kaolinite and illite as the principal clay minerals. The sediments are relatively poor in silica (40–45%) but rich in Fe₂O₃ (10%). Both mineral proportions and chemical composition change systematically over the area of the lake and delineate four sedimentological provinces: (1) iron-rich, silty kaolinitic muds (Omo Delta); (2) iron-rich, fine-grained montmorillonite muds (North Basin); (3) silty montmorillonite muds rich in Na₂O and K₂O (Central Delta); and (4) argillaceous calcite silts (South Basin). Omo Delta and North Basin sediments are derived from the volcanics of the Ethiopian plateau; the source of Central Delta sediments is the Precambrian metamorphic terrain of the rift valley margin; the South Basin has a restricted detrital input. The water in the lake is alkaline (pH 9.2) and moderately saline (TDS = 2500 p.p.m.). Comparisons with influent water from the Omo River indicate a 200-fold concentration for the lake water. Models based on equilibrium between sediments and water column account for most of the non-conservative chemical components in the lake water. Sedimentation rates are high (about 1 m per 1000 years) and the dominance of detrital sediments makes Lake Rudolf unusual in comparison with other closed-basin lakes in the African Rift Valley although some similarities with ancient rift valley deposits are suggested.     

Calcareous crust (caliche) profiles and early subaerial exposure of Carboniferous carbonates, northeastern Kentucky

Sequences of laminated limestones found within thin Carboniferous carbonate strata of northeastern Kentucky were studied to determine their origin and palaeo-environmental significance. These laminated zones are strikingly similar to Holocene and Pleistocene surficial calcareous crusts (caliche) profiles that occur in various parts of the world. Carboniferous laminated carbonates are associated with shallow marine carbonate units, palaeokarst, and overlying palaeosol zones. A typical laminated profile ranges in thickness from 1 to 2 m and contains brecciated, light olive-grey to brown micrite that lacks distinctive bedding. Structures and textures common in most profiles include: (1) calcareous and silicious laminae (laminae form diffuse, alternating light and dark bands that generally parallel bedding but often fill fractures and vugs within the rocks); (2) particles (allochems, and micrite and microspar fragments) coated by brown microcrystalline calcite; (3) brecciated texture; (4) circular to elliptical fossil moulds (occur in sinuous patterns and fill fractures within the rocks); (5) large and small scale fracture patterns. Subaerial weathering and vadose diagenesis of carbonate mud banks or islands is suggested as a mechanism for the formation of these Carboniferous calcareous crust profiles. These ‘crusts’ formed by a combination of solution (karsting), brecciation, and soil development that transformed an exposed marine biomicrite (‘host’ rock) into a porous subsoil rubble. Laminated ‘crusts’ and coated particles developed as the result of dissolution and reprecipitation of CaCO₃ and SiO₃ from the soil and carbonate bedrock. Carboniferous laminated carbonates in northeastern Kentucky are often referred to as ‘algal limestones’ because of their superficial similarity to some modern and ancient algal structures. This study, however, reveals numerous characteristics that can only be explained by diagenesis in a subaerial environment.     

Phonolitic Diatremes within the Dunedin Volcano, South Island, New Zealand

The Port Chalmers Breccia is a vent-filling, clastic volcanicunit exposed within the Miocene Dunedin Volcano of South Island,New Zealand. Clasts (up to in excess of 1 m but generally <20cm) are supported in ash and fine lapilli of phonolitic (ne-benmoreiteor tephro-phonolite) composition and the dominant clast type(55 to almost 100%) is also phonolitic. Less abundant lithologiesinclude ne-normative basalt (basanite), hawaiite, mugeariteand trachyandesite, syenites and microsyenites, coarse-grainedmafic (gabbros) and ultramafic rocks (pyroxenites, hornblendites),schists and sediments. The breccias were emplaced as diatremesassociated with localized, but highly explosive, eruptive eventsin which mantle-derived CO₂ was an important component. Thesyenitic and ultramafic clasts could represent intrusive suitesproduced by crystal fractionation acting on parental ne-benmoreitemagmas that may themselves have been derived by crystal fractionationfrom basanitic precursors. An alternative variation on thismodel is that the parental ne-benmoreites were generated throughpartial melting of an alkalic igneous underplate. Sr, Nd andPb isotopic compositions are strikingly similar to those ofintraplate igneous rocks, ranging in age from 100 to less than10 Ma, from elsewhere in the South Island, and New Zealand'ssub-Antarctic islands, the south Tasman Sea and the Ross Searegion. This regional, HIMU-influenced, isotopic signature isbelieved to be derived from within the lithospheric mantle. KEY WORDS: phonolite; diatreme; nepheline syenite; Dunedin Volcano; alkalic rocks; fractional crystallization     

The Role of Open-System Processes in the Development of Silicic Magma Chambers: a Chemical and Isotopic Investigation of the Fogo A Trachyte Deposit, Sao Miguel, Azores

The processes operating in the development of chemical zonationin silicic magma chambers have been addressed with a Sr–Nd–Pb–Hf–Thisotope study of the chemically zoned trachyte pumice depositof the Fogo A eruption, Fogo volcano, Azores. Sr isotopic variationis observed in whole rocks, glass separates and sanidine phenocrysts(whole-rock ⁸⁷Sr/⁸⁶Sr: 0·7049–0·7061; glass⁸⁷Sr/⁸⁶Sr: 0·7048–0·7052; sanidine ⁸⁷Sr/⁸⁶Sr:0·7048–0·7062). Thorium isotopic variationis observed in glass separates, with (²³⁰Th/²³²Th)_o rangingfrom 0·8737 to 0·8841, and exhibiting a negativecorrelation with Sr isotopes. The Nd, Pb and Hf isotopic compositionsof the whole-rock trachytic pumices are invariant and indistinguishablefrom basalts flanking the volcano. The Sr isotope variationsin the whole rocks are proposed to be the result of three distinctprocesses: contamination of the Fogo A magma by assimilationof radiogenic seawater-altered syenite wall rock, to explainthe Sr and Th isotopic compositions of the glass separates;incorporation of xenocrysts into the trachytic magma, requiredto explain the range in feldspar Sr isotopic compositions; andpost-eruptive surface alteration. This study emphasizes theimportance of determining the isotopic composition of glassand mineral separates rather than whole rocks when pre-eruptivemagmatic processes are being investigated. KEY WORDS: Azores; open-system processes; Sr isotopes; trachytic pumices; zoned magma chambers