P‐wave stacking‐velocity tomography for VTI media

A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity V_P0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces.     

High‐temperature gold‐copper extraction with chloride flux in lava tubes of Tolbachik volcano (Kamchatka)

Subvolcanic environments in supra‐subduction zones are renowned for hosting epithermal deposits that often contain electrum and native gold, including bonanza examples. This study examined mineral assemblages and processes occurring in shallow‐crust volcanic settings using recent eruption (2012–2013) of the basaltic Tolbachik volcano in the Kamchatka arc. The Tolbachik eruptive system is characterized by an extensive system of lava tubes. After cessation of magma input, the tubes maintained the flow of hot oxidized gases that episodically interacted with the lava surfaces and sulphate‐chloride precipitates from volcanic gases on these surfaces. The gas‐rock interaction had strong pyrometamorphic effects that resulted in the formation of molten salt, oxidized (tenorite, hematite, Cu‐rich magnesioferrite) and skarn‐like silicate mineral assemblages. By analogy with experimental studies, we propose that a combination of these processes was responsible for extraction of metals from the basaltic wall rocks and deposition of Cu‐, Fe‐ and Cu‐Fe‐oxides and native gold.     

An Anisotropic Cosmic-Ray Enhancement Event on 07-June-2015: A Possible Origin

A usual event, called anisotropic cosmic-ray enhancement (ACRE), was observed as a small increase (\({\leq}\,5\%\)) in the count rates of polar neutron monitors during 12?–?19 UT on 07 June 2015. The enhancement was highly anisotropic, as detected only by neutron monitors with asymptotic directions in the southwest quadrant in geocentric solar ecliptic (GSE) coordinates. The estimated rigidity of the corresponding particles is \({\leq}\,1\) GV. No associated detectable increase was found in the space-borne data from the Geostationary Operational Environmental Satellite (GOES), the Energetic and Relativistic Nuclei and Electron (ERNE) on board the Solar and Heliospheric Observatory (SOHO), or the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) instruments, whose sensitivity was not sufficient to detect the event. No solar energetic particles were present during that time interval. The heliospheric conditions were slightly disturbed, so that the interplanetary magnetic field strength gradually increased during the event, followed by an increase of the solar wind speed after the event. It is proposed that the event was related to a crossing of the boundary layer between two regions with different heliospheric parameters, with a strong gradient of low-rigidity (\({<}\,1\) GV) particles. It was apparently similar to another cosmic-ray enhancement (e.g., on 22 June 2015) that is thought to have been caused by the local anisotropy of Forbush decreases, with the difference that in our case, the interplanetary disturbance was not observed at Earth, but passed by southward for this event.     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

Origins and evolution of rhyolitic magmas in the central Snake River Plain: insights from coupled high-precision geochronology,oxygen isotope,and hafnium isotope analyses of zircon

We present new high-precision CA-ID-TIMS and in situ U–Pb ages together with Hf and O isotopic analyses (analyses performed all on the same grains) from four tuffs from the 15?10 Ma Bruneau–Jarbidge center of the Snake River Plain and from three rhyolitic units from the Kimberly borehole in the neighboring 10?6 Ma Twin Falls volcanic center. We find significant intrasample diversity in zircon ages (ranges of up to 3 Myr) and in δ¹⁸O (ranges of up to 6‰) and ε_Hf (ranges of up to 24 ε units) values. Zircon rims are also more homogeneous than the associated cores, and we show that zircon rim growth occurs faster than the resolution of in situ dating techniques. CA-ID-TIMS dating of a subset of zircon grains from the Twin Falls samples reveals complex crystallization histories spanning 10⁴–10⁶ years prior to some eruptions, suggesting that magma genesis was characterized by the cyclic remelting of buried volcanic rocks and intrusions associated with previous magmatic episodes. Age-dependent trends in zircon isotopic compositions show that rhyolite production in the Yellowstone hotspot track is driven by the mixing of mantle-derived melts (normal δ¹⁸O and ε_Hf) and a combination of Precambrian basement rock (normal δ¹⁸O and ε_Hf down to ??60) and shallow Mesozoic and Cenozoic age rocks, some of which are hydrothermally altered (to low δ¹⁸O values) by earlier stages of Snake River Plain magmatism. These crustal melts hybridize with juvenile basalts and rhyolites to produce the erupted rhyolites. We also observe that the Precambrian basement rock is only an important component in the erupted magmas in the first eruption at each caldera center, suggesting that the accumulation of new intrusions quickly builds an upper crustal intrusive body which is isolated from the Precambrian basement and evolves towards more isotopically juvenile and lower-δ¹⁸O compositions over time.     

Oxygen isotope study of the Long Valley magma system, California: isotope thermometry and convection in large silicic magma bodies

Products of voluminous pyroclastic eruptions with eruptive draw-down of several kilometers provide a snap-shot view of batholith-scale magma chambers, and quench pre-eruptive isotopic fractionations (i.e., temperatures) between minerals. We report analyses of oxygen isotope ratio in individual quartz phenocrysts and concentrates of magnetite, pyroxene, and zircon from individual pumice clasts of ignimbrite and fall units of caldera-forming 0.76 Ma Bishop Tuff (BT), pre-caldera Glass Mountain (2.1-0.78 Ma), and post-caldera rhyolites (0.65-0.04 Ma) to characterize the long-lived, batholith-scale magma chamber beneath Long Valley Caldera in California. Values of '¹⁸O show a subtle 1‰ decrease from the oldest Glass Mountain lavas to the youngest post-caldera rhyolites. Older Glass Mountain lavas exhibit larger (~1‰) variability of '¹⁸O(quartz). The youngest domes of Glass Mountain are similar to BT in '¹⁸O(quartz) values and reflect convective homogenization during formation of BT magma chamber surrounded by extremely heterogeneous country rocks (ranging from 2 to +29‰). Oxygen isotope thermometry of BT confirms a temperature gradient between "Late" (815 °C) and "Early" (715 °C) BT. The '¹⁸O(quartz) values of "Early" and "Late" BT are +8.33 and 8.21‰, consistent with a constant '¹⁸O(melt)=7.8ǂ.1‰ and 100 °C temperature difference. Zircon-melt saturation equilibria gives a similar temperature range. Values of '¹⁸O(quartz) for different stratigraphic units of BT, and in pumice clasts ranging in pre-eruptive depths from 6 to 11 km (based on melt inclusions), and document vertical and lateral homogeneity of '¹⁸O(melt). Worldwide, five other large-volume rhyolites, Lava Creek, Lower Bandelier, Fish Canyon, Cerro Galan, and Toba, exhibit equal '¹⁸O(melt) values of earlier and later erupted portions in each of the these climactic caldera-forming eruptions. We interpret the large-scale '¹⁸O homogeneity of BT and other large magma chambers as evidence of their longevity (>10⁵ years) and convection. However, remaining isotopic zoning in some quartz phenocrysts, trace element gradients in feldspars, and quartz and zircon crystal size distributions are more consistent with far shorter timescales (10²-10⁴ years). We propose a sidewall-crystallization model that promotes convective homogenization, roofward accumulation of more evolved and stagnant, volatile-rich liquid, and develops compositional and temperature gradients in pre-climactic magma chamber. Crystal + melt + gas bubbles mush near chamber walls of variable '¹⁸O gets periodically remobilized in response to chamber refill by new hotter magmas. One such episode of chamber refill by high-Ti, Sr, Ba, Zr, and volatile-richer magma happened 10³-10⁴ years prior to the 0.76-Ma caldera collapse that caused magma mixing at the base, mush thawing near the roof and walls, and downward settling of phenocrysts into this hybrid melt.     

Monocotyledon-like Leaves from the Middle Jurassic of East Siberia (Russia)

Leaf fragments of Herbifolia antiqua A. Frolov et Enushchenko gen. et sp. nov. from the Middle Jurassic (Aalenian) of the Irkutsk Coal Basin (Eastern Siberia, Russia) are reported. The following features are characteristic of H. antiqua gen. et sp. nov.: the presence of leaf sheaths, linear smooth-edged leaves with parallel venation, anastomoses between the veins, anomocytic stomata, rhomb-shaped ordinary epidermal cells. Such a combination of characters is widespread in modern monocotyledonous plants and is absent in fossil and modern cryptogam and gymnosperms. Due to the lack of evidence of an enclosed ovule in H. antiqua gen. et sp. nov.we attribute it to a typological angiosperm, based on its unique leaf structure characteristic of monocotyledons. The leaf epidermal structure of Herbifolia gen. nov. is most similar to those of modern Asparagales and Liliales.     

Experimental calibration of oxygen isotope fractionation between quartz and zircon

We report the results of an experimental calibration of oxygen isotope fractionation between quartz and zircon. Data were collected from 700 to 1000 °C, 10–20 kbar, and in some experiments the oxygen fugacity was buffered at the fayalite–magnetite–quartz equilibrium. Oxygen isotope fractionation shows no clear dependence on oxygen fugacity or pressure. Unexpectedly, some high-temperature data (900–1000 °C) show evidence for disequilibrium oxygen isotope partitioning. This is based in part on ion microprobe data from these samples that indicate some high-temperature quartz grains may be isotopically zoned. Excluding data that probably represent non-equilibrium conditions, our preferred calibration for oxygen isotope fractionation between quartz and zircon can be described by:

This relationship can be used to calculate fractionation factors between zircon and other minerals. In addition, results have been used to calculate WR/melt–zircon fractionations during magma differentiation. Modeling demonstrates that silicic magmas show relatively small changes in δ¹⁸O values during differentiation, though late-stage mafic residuals capable of zircon saturation contain elevated δ¹⁸O values. However, residuals also have larger predicted melt–zircon fractionations meaning zircons will not record enriched δ¹⁸O values generally attributed to a granitic protolith. These results agree with data from natural samples if the zircon fractionation factor presented here or from natural studies is applied.     

Inversion of multicomponent seismic time shifts for reservoir pressure and length: a feasibility study

Pressure drops associated with reservoir production generate excess stress and strain that cause travel‐time shifts of reflected waves. Here, we invert time shifts of P‐, S‐, and PS‐waves measured between baseline and monitor surveys for pressure reduction and reservoir length. The inversion results can be used to estimate compaction‐induced stress and strain changes around the reservoir. We implement a hybrid inversion algorithm that incorporates elements of gradient, global/genetic, and nearest neighbour methods and permits exploration of the parameter space while simultaneously following local misfit gradients. Our synthetic examples indicate that optimal estimates of reservoir pressure from P‐wave data can be obtained using the reflections from the reservoir top. For S‐waves, time shifts from the top of the reservoir can be accurately inverted for pressure if the noise level is low. However, if noise contamination is significant, it is preferable to use S‐wave data (or combined shifts of all three modes) from reflectors beneath the reservoir. Joint wave type inversions demonstrate improvements over any single pure mode. Reservoir length can be estimated using the time shifts of any mode from the reservoir top or deeper reflectors. We also evaluate the differences between the actual strain field and those corresponding to the best‐case inversion results obtained using P‐ and S‐wave data. Another series of tests addresses the inversion of the time shifts for the pressure drops in two‐compartment reservoirs, as well as for the associated strain field. Numerical testing shows that a potentially serious source of error in the inversion is a distortion in the strain‐sensitivity coefficients, which govern the magnitude of stiffness changes. This feasibility study suggests which wave types and reflector locations may provide the most accurate estimates of reservoir parameters from compaction‐induced time shifts.     

Solar modulation in surface atmospheric electricity

The solar wind modulates the flux of galactic cosmic rays impinging on Earth inversely with solar activity. Cosmic ray ionisation is the major source of air's electrical conductivity over the oceans and well above the continents. Differential solar modulation of the cosmic ray energy spectrum modifies the cosmic ray ionisation at different latitudes, varying the total atmospheric columnar conductance. This redistributes current flow in the global atmospheric electrical circuit, including the local vertical current density and the related surface potential gradient. Surface vertical current density and potential gradient measurements made independently at Lerwick Observatory, Shetland, from 1978 to 1985 are compared with modelled changes in cosmic ray ionisation arising from solar activity changes. Both the lower troposphere atmospheric electricity quantities are significantly increased at cosmic ray maximum (solar minimum), with a proportional change greater than that of the cosmic ray change.