Migration velocity analysis for tilted transversely isotropic media

Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold‐and‐thrust belts) and in subsalt exploration. Here, we introduce a methodology for P‐wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters ε and δ and linearly varying symmetry‐direction velocity V_P0 controlled by the vertical (k_z) and lateral (k_x) gradients. Since determination of tilt from P‐wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity V_P0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients k_z and k_x and the anisotropy parameters ε and δ in the layer‐stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry‐axis direction is fixed and V_P0 is known, the parameters k_z, k_x, ε and δ can be resolved from reflection data. It should be emphasized that estimation of ε in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas.     

The origin of hydrous,high-δ<Superscript>18</Superscript>O voluminous volcanism: diverse oxygen isotope values and high magmatic water contents within the volcanic record of Klyuchevskoy volcano,Kamchatka, Russia

Klyuchevskoy volcano, in Kamchatka’s subduction zone, is one of the most active arc volcanoes in the world and contains some of the highest δ¹⁸O values for olivines and basalts. We present an oxygen isotope and melt inclusion study of olivine phenocrysts in conjunction with major and trace element analyses of ¹⁴C- and tephrochronologically-dated tephra layers and lavas spanning the eruptive history of Klyuchevskoy. Whole-rock and groundmass analyses of tephra layers and lava samples demonstrate that both high-Mg (7–12.5 wt% MgO) and high-Al (17–19 wt% Al₂O₃, 3–6.5 wt% MgO) basalt and basaltic andesite erupted coevally from the central vent and flank cones. Individual and bulk olivine δ¹⁸O range from normal MORB values of 5.1‰ to values as high as 7.6‰. Likewise, tephra and lava matrix glass have high-δ¹⁸O values of 5.8–8.1‰. High-Al basalts dominate volumetrically in Klyuchevskoy’s volcanic record and are mostly high in δ¹⁸O. High-δ¹⁸O olivines and more normal-δ¹⁸O olivines occur in both high-Mg and high-Al samples. Most olivines in either high-Al or high-Mg basalts are not in oxygen isotopic equilibrium with their host glasses, and Δ¹⁸O_{olivine–glass} values are out of equilibrium by up to 1.5‰. Olivines are also out of Fe–Mg equilibrium with the host glasses, but to a lesser extent. Water concentrations in olivine-hosted melt inclusions from five tephra samples range from 0.4 to 7.1 wt%. Melt inclusion CO₂ concentrations vary from below detection (<50 ppm) to 1,900 ppm. These values indicate depths of crystallization up to ~17 km (5 kbar). The variable H₂O and CO₂ concentrations likely reflect crystallization of olivine and entrapment of inclusions in ascending and degassing magma. Oxygen isotope and Fe–Mg disequilibria together with melt inclusion data indicate that olivine was mixed and recycled between high-Al and high-Mg basaltic melts and cumulates, and Fe–Mg and δ¹⁸O re-equilibration processes were incomplete. Major and trace elements in the variably high-δ¹⁸O olivines suggest a peridotite source for the parental magmas. Voluminous, highest in the world with respect to δ¹⁸O, and hydrous basic volcanism in Klyuchevskoy and other Central Kamchatka depression volcanoes is explained by a model in which the ascending primitive melts that resulted from the hydrous melt fluxing of mantle wedge peridotite, interacted with the shallow high-δ¹⁸O lithospheric mantle that had been extensively hydrated during earlier times when it was part of the Kamchatka forearc. Following accretion of the Eastern Peninsula terrains several million years ago, a trench jump eastward caused the old forearc mantle to be beneath the presently active arc. Variable interaction of ascending flux-melting-derived melts with this older, high-δ¹⁸O lithospheric mantle has produced mafic parental magmas with a spectrum of δ¹⁸O values. Differentiation of the higher δ¹⁸O parental magmas has created the volumetrically dominant high-Al basalt series. Both basalt types incessantly rise and mix between themselves and with variable in δ¹⁸O cumulates within dynamic Klyuchevskoy magma plumbing system, causing biannual eruptions and heterogeneous magma products. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolution of crystal sizes in the series of dissolution and precipitation events in open magma systems

We propose a model that describes the evolution of crystal sizes and crystal size distributions (CSD) of igneous phenocrysts in a sequence of dissolution and crystallization events. This model is based on the assumption that crystal dissolution is rate-limited by diffusion in melt while crystal growth is controlled by the slower kinetic of new nucleation and growth. As a result, the dissolution rate is inversely proportional to crystal size coming into effect through the curvature of the crystal's surface, but the growth rate does not depend on the crystal size. Closed-form analytical solution of equation for CSD is obtained. We apply results of modeling to quartz and zircon, two prime minerals in silicic igneous systems that are widely used in geochemical and isotopic investigations. The time-series of multiple solution–reprecipitation episodes generate concave-downward CSDs and this result fits well with experimental and natural observations on the abundant concave-down CSD in silicic igneous rocks. We suggest that maturation of crystal populations with sizes above several micrometers can not be caused by a size effect on the solubility of the crystals (Ostwald ripening), but is rather driven by thermal oscillations in experiments and in nature. The model predicts that mean crystal size increases with time proportionally to ∼ t^0.20, which is very close to the published experimental results for quartz maturation with the exponent of 0.19–0.22. Our proposed model gives an opportunity to use natural CSDs for interpretation of pre-eruptive magma history, when solubilities and diffusion data are available for constituent elements of the dissolving mineral. In particular, we present time estimates for maturing zircon populations in large volume ignimbrites and estimate that it takes 100–1000 yrs to mature an initially exponential CSD to a lognormal CSD.     

Geochemistry of a new enriched mantle type locality in the northern hemisphere: Implications for the origin of the EM-I source

Late Cretaceous (66.2 ± 0.5 Ma amphibole and 66.7 ± 0.2 Ma phlogopite ⁴⁰Ar/³⁹Ar ages) nephelinitic volcanic rocks from Godzilla Seamount in the eastern North Atlantic (34°N latitude) have trace element and Sr–Nd–Pb–Hf-isotope compositions similar to the Enriched Mantle I (EM-I) endmember, except for their low ²⁰⁷Pb/²⁰⁴Pb relative to ²⁰⁶Pb/²⁰⁴Pb ratios (²⁰⁶Pb/²⁰⁴Pb_in = 17.7, ²⁰⁷Pb/²⁰⁴Pb_in = 15.34) plotting below the Northern Hemisphere Reference Line on the uranogenic Pb isotope diagram. O isotope data on amphibole separates are mantle-like (δ¹⁸O = 5.6–5.8‰). Age and location of the isolated Godzilla Seamount, however, preclude it from being derived from the Madeira or Canary hotspots, making a lower-mantle origin unlikely. Therefore we propose derivation from a shallow (lithospheric/asthenospheric) melting anomaly. As observed in mid-ocean-ridge and ocean-island basalts, there is a systematic decrease of ²⁰⁷Pb/²⁰⁴Pb ratios (and Δ7/4) in the individual EM-I endmember type localities towards northern latitudes with Godzilla lying on the extension of this trend. This trend is mirrored in ultra-potassic volcanic rocks such as lamproites and kimberlites, which reflect the composition of enriched subcontinental lithospheric mantle. Therefore, a global pattern in ²⁰⁷Pb/²⁰⁴Pb ratios and Δ7/4 is suggested. The geochemical composition of EM-I endmember type localities, including Godzilla lavas, and the enriched (DUPAL) anomaly in the southern hemisphere could reflect derivation from ancient, metasomatized subcontinental lithospheric mantle. We propose a two-stage model to explain the trace element and isotopic composition of the EM-I mantle endmember localities worldwide: 1) during the early history of the Earth, subcontinental lithosphere was metasomatized by melts from subducted slabs along convergent margins generating high μ (²³⁸U/²⁰⁴Pb) sources, and 2) as the Earth cooled, hydrous fluids replaced hydrous melts as the main slab component metasomatizing the subcontinental lithospheric mantle (generating EM-I sources with lower μ). In accordance with this model, the global variations in ²⁰⁷Pb/²⁰⁴Pb ratios and Δ7/4 could reflect geographic differences in μ and/or the age at which the transition from stages 1 to 2 took place in the Archaean lithosphere. The model would require a re-definition of the EM-I endmember to low ²⁰⁶Pb/²⁰⁴Pb, high ²⁰⁸Pb/²⁰⁴Pb (positive Δ8/4) but variable ²⁰⁷Pb/²⁰⁴Pb (positive and negative Δ7/4).     

Generalized Dix equation and analytic treatment of normal-moveout velocity for anisotropic media*

Despite the complexity of wave propagation in anisotropic media, reflection moveout on conventional common-midpoint (CMP) spreads is usually well described by the normal-moveout (NMO) velocity defined in the zero-offset limit. In their recent work, Grechka and Tsvankin showed that the azimuthal variation of NMO velocity around a fixed CMP location generally has an elliptical form (i.e. plotting the NMO velocity in each azimuthal direction produces an ellipse) and is determined by the spatial derivatives of the slowness vector evaluated at the CMP location. This formalism is used here to develop exact solutions for the NMO velocity in anisotropic media of arbitrary symmetry. For the model of a single homogeneous layer above a dipping reflector, we obtain an explicit NMO expression valid for all pure modes and any orientation of the CMP line with respect to the reflector strike. The contribution of anisotropy to NMO velocity is contained in the slowness components of the zero-offset ray (along with the derivatives of the vertical slowness with respect to the horizontal slownesses) — quantities that can be found in a straightforward way from the Christoffel equation. If the medium above a dipping reflector is horizontally stratified, the effective NMO velocity is determined through a Dix-type average of the matrices responsible for the ‘interval’ NMO ellipses in the individual layers. This generalized Dix equation provides an analytic basis for moveout inversion in vertically inhomogeneous, arbitrarily anisotropic media. For models with a throughgoing vertical symmetry plane (i.e. if the dip plane of the reflector coincides with a symmetry plane of the overburden), the semi-axes of the NMO ellipse are found by the more conventional rms averaging of the interval NMO velocities in the dip and strike directions. Modelling of normal moveout in general heterogeneous anisotropic media requires dynamic ray tracing of only one (zero-offset) ray. Remarkably, the expressions for geometrical spreading along the zero-offset ray contain all the components necessary to build the NMO ellipse. This method is orders of magnitude faster than multi-azimuth, multi-offset ray tracing and, therefore, can be used efficiently in traveltime inversion and in devising fast dip-moveout (DMO) processing algorithms for anisotropic media. This technique becomes especially efficient if the model consists of homogeneous layers or blocks separated by smooth interfaces. The high accuracy of our NMO expressions is illustrated by comparison with ray-traced reflection traveltimes in piecewise-homogeneous, azimuthally anisotropic models. We also apply the generalized Dix equation to field data collected over a fractured reservoir and show that P-wave moveout can be used to find the depth-dependent fracture orientation and to evaluate the magnitude of azimuthal anisotropy.     

Phase equilibria in the silica-undersaturated part of the KAlSiO4– Mg2SiO4– Ca2SiO4– SiO2– F system at 1 atm and the larnite-normative trend of melt evolution

The evolution of nephelinitic melts in equilibrium with mica-bearing liquidus assemblages and melting relations have been studied on two silica-undersaturated joins of the KAlSiO₄– Mg₂SiO₄– Ca₂SiO₄– SiO₂– F system at atmospheric pressure by quench runs in sealed platinum capsules. Fluorine has been added to the batch compositions by the direct exchange of fluorine for oxygen (2F⁻ = O²⁻). The first join is the pseudo-ternary Forsterite – Diopside – KAlSiO₃F₂ system. Forsterite, diopside, F-phlogopite and leucite crystallisation fields and a fluoride-silicate liquid immiscibility solvus are present on the liquidus surface of the join. Sub-liquidus and sub-solidus phases include akermanite, cuspidine, spinel, fluorite and some other minor fluorine phases. The second system is the pseudo-binary Akermanite – F-phlogopite join that intersects the Forsterite – Diopside – KAlSiO₃F₂ join. Akermanite, forsterite, diopside, F-phlogopite, leucite and cuspidine are found to crystallise on the join. Forsterite (fo) and leucite (lc) are related to F-phlogopite (phl) by a reaction with the fluorine-bearing liquid: fo + lc + l = phl, and the reaction proceeds until forsterite or leucite are completely consumed. The reaction temperature and resulting phase association depend on batch composition. Thus, leucite is not stable in the sub-solidus of the Akermanite – F-phlogopite join, but is preserved in a part of the Forsterite – Diopside – KAlSiO₃F₂ system where forsterite reacts out, or does not crystallise at all. The phlogopite-in reaction has an important effect on the composition of the coexisting liquid. The liquids initially saturated in forsterite evolve to extremely Ca rich, larnite-normative residuals. The experimental data show that larnite-normative melilitolites can crystallise from evolved melilititic melts generated from “normal” melanephelinitic parental magmas with no normative larnite. The evolution towards melilitites requires fractionation of phlogopite-bearing assemblages under volatile pressure. Received: 3 June 1997 / Accepted: 5 January 1998     

Evidence for storm-dominated early progradation of Castle Neck barrier, Massachusetts, USA

Washovers, dune scarps and flattened beach profiles with concentrations of coarse-grained sediment or heavy minerals are the diagnostic geological signatures of large storms on barriers today. It is clear that storms are a major force driving transgressive barriers onshore, but what is not as well understood is the role these powerful erosive events play in the evolution of prograding barriers. Application of ground-penetrating radar (GPR) and a combination of coring techniques have significantly improved our ability to image barrier architecture. Results of these studies reveal a more complex evolution than previously recognized. It is now possible to precisely locate and map storm deposits within prograded barrier lithosomes.

A comprehensive study of northern Castle Neck, Massachusetts was performed using 15 km of GPR surveys, a 120-m-long seismic line, 11 cores, and several radiocarbon dates. Storm-related layers are the most prominent horizons contained in the barrier stratigraphy. The geometry and sedimentology of these layers closely resembles those of a present-day post-storm beach. Twenty closely spaced, curvilinear heavy mineral layers imaged in the landward portion of the barrier suggest that the Castle Neck barrier was heavily influenced by storms during its initial phase of progradation beginning 4000 years BP. Approximately 1800 years BP, two intense storms impacted the coast depositing two extensive coarse-grained units. These layers mimic the flat-lying sand and gravel deposits that occur in front of a nearby eroding till outcrop following major storms. The great number of storm deposits in the early history of Castle Neck is related to either a period of greater storm activity and/or a slow rate of barrier progradation. The occurrence and preservation of these earlier storm layers are likely a product of the exposure of nearby drumlins resulting in greater availability of iron oxide and ferromagnesian sands. The supply of heavy-mineral sands has gradually diminished as the barrier prograded and the proximal drumlin source was buried by beach and dune sands.     

Large-volume silicic volcanism in Kamchatka: Ar–Ar and U–Pb ages,isotopic, and geochemical characteristics of major pre-Holocene caldera-forming eruptions

The Kamchatka Peninsula in far eastern Russia represents the most volcanically active arc in the world in terms of magma production and the number of explosive eruptions. We investigate large-scale silicic volcanism in the past several million years and present new geochronologic results from major ignimbrite sheets exposed in Kamchatka. These ignimbrites are found in the vicinity of morphologically-preserved rims of partially eroded source calderas with diameters from ～ 2 to ～ 30 km and with estimated volumes of eruptions ranging from 10 to several hundred cubic kilometers of magma. We also identify and date two of the largest ignimbrites: Golygin Ignimbrite in southern Kamchatka (0.45 Ma), and Karymshina River Ignimbrites (1.78 Ma) in south-central Kamchatka. We present whole-rock geochemical analyses that can be used to correlate ignimbrites laterally. These large-volume ignimbrites sample a significant proportion of remelted Kamchatkan crust as constrained by the oxygen isotopes. Oxygen isotope analyses of minerals and matrix span a 3‰ range with a significant proportion of moderately low-δ¹⁸O values. This suggests that the source for these ignimbrites involved a hydrothermally-altered shallow crust, while participation of the Cretaceous siliceous basement is also evidenced by moderately elevated δ¹⁸O and Sr isotopes and xenocryst contamination in two volcanoes. The majority of dates obtained for caldera-forming eruptions coincide with glacial stages in accordance with the sediment record in the NW Pacific, suggesting an increase in explosive volcanic activity since the onset of the last glaciation 2.6 Ma. Rapid changes in ice volume during glacial times and the resulting fluctuation of glacial loading/unloading could have caused volatile saturation in shallow magma chambers and, in combination with availability of low-δ¹⁸O glacial meltwaters, increased the proportion of explosive vs effusive eruptions. The presented results provide new constraints on Pliocene–Pleistocene volcanic activity in Kamchatka, and thus constrain an important component of the Pacific Ring of Fire.     

Isotopically diverse rhyolites coeval with the Columbia River Flood Basalts: evidence for mantle plume interaction with the continental crust

The Columbia River Flood Basalts (CRB) of the northwestern USA are coeval with eruptions of several thousand km³ of rhyolite. A broad survey of major phenocryst oxygen isotopes and of O and Hf isotopes in zircons from these rhyolites reveals significant diversity in inferred δ¹⁸O_melt values, ranging from +1.9 to +10.5‰ (SMOW), and in zircon Hf isotope compositions, which range from ε_Hf = ?39 to +9. This newly identified isotopic diversity shows that the syn‐CRB rhyolites were derived from high‐percentage melting of the crust. Low‐δ¹⁸O rhyolites, which fingerprint the melting of hydrothermally altered crust, are concentrated at the edge of the North American craton. This suggests that the conditions of crustal heating, faulting, and hydrothermal alteration required for the production of these rhyolites were concentrated there by the contrasts in crustal thickness and rheology associated with the boundary between the North American craton and younger accreted terranes.     

Cumulate xenoliths from St. Vincent, Lesser Antilles Island Arc: a window into upper crustal differentiation of mantle-derived basalts

In order to shed light on upper crustal differentiation of mantle-derived basaltic magmas in a subduction zone setting, we have determined the mineral chemistry and oxygen and hydrogen isotope composition of individual cumulus minerals in plutonic blocks from St. Vincent, Lesser Antilles. Plutonic rock types display great variation in mineralogy, from olivine–gabbros to troctolites and hornblendites, with a corresponding variety of cumulate textures. Mineral compositions differ from those in erupted basaltic lavas from St. Vincent and in published high-pressure (4–10 kb) experimental run products of a St. Vincent high-Mg basalt in having higher An plagioclase coexisting with lower Fo olivine. The oxygen isotope compositions (δ¹⁸O) of cumulus olivine (4.89–5.18‰), plagioclase (5.84–6.28‰), clinopyroxene (5.17–5.47‰) and hornblende (5.48–5.61‰) and hydrogen isotope composition of hornblende (δD = −35.5 to −49.9‰) are all consistent with closed system magmatic differentiation of a mantle-derived basaltic melt. We employed a number of modelling exercises to constrain the origin of the chemical and isotopic compositions reported. δ¹⁸O_Olivine is up to 0.2‰ higher than modelled values for closed system fractional crystallisation of a primary melt. We attribute this to isotopic disequilibria between cumulus minerals crystallising at different temperatures, with equilibration retarded by slow oxygen diffusion in olivine during prolonged crustal storage. We used melt inclusion and plagioclase compositions to determine parental magmatic water contents (water saturated, 4.6 ± 0.5 wt% H₂O) and crystallisation pressures (173 ± 50 MPa). Applying these values to previously reported basaltic and basaltic andesite lava compositions, we can reproduce the cumulus plagioclase and olivine compositions and their associated trend. We conclude that differentiation of primitive hydrous basalts on St. Vincent involves crystallisation of olivine and Cr-rich spinel at depth within the crust, lowering MgO and Cr₂O₃ and raising Al₂O₃ and CaO of residual melt due to suppression of plagioclase. Low density, hydrous basaltic and basaltic andesite melts then ascend rapidly through the crust, stalling at shallow depth upon water saturation where crystallisation of the chemically distinct cumulus phases observed in this study can occur. Deposited crystals armour the shallow magma chamber where oxygen isotope equilibration between minerals is slowly approached, before remobilisation and entrainment by later injections of magma.