Gravity and magnetic data inversion for 3D topography of the Moho discontinuity in the northern Red Sea area, Egypt

The main goal of our study is to investigate 3D topography of the Moho boundary for the area of the northern Red Sea including Gulf of Suez and Gulf of Aqaba. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. To separate sources of gravity and magnetic field in depth, a method is suggested, based on upward and downward continuation. Both new methods are applied to isolate the contribution of the Moho interface to the total field and to find its 3D topography. At the first stage, we separate near-surface and deeper sources. According to the obtained field of shallow sources a model of the horizontal layer above the depth of 7 km is suggested, which includes a density interface between light sediments and crystalline basement. Its depressions and uplifts correspond to known geological structures. At the next stage, we isolate the effect of very deep sources (below 100 km) and sources outside the area of investigation. After subtracting this field from the total effect of deeper sources, we obtain the contribution of the Moho interface. We make inversion separately for the area of rifts (Red Sea, Gulf of Suez and Gulf of Aqaba) and for the rest of the area. In the rift area we look for the upper boundary of low-density, heated anomalous upper mantle. In the rest of the area the field is satisfied by means of topography for the interface between lower crust and normal upper mantle. Both algorithms are applied also to the magnetic field. The magnetic model of the Moho boundary is in agreement with the gravitational one.     

Ionization of the earth’s atmosphere by solar and galactic cosmic rays

A brief review of the research of atmospheric effects of cosmic rays is presented. Numerical models are discussed, that are capable to compute the cosmic ray induced ionization at a given location and time. Intercomparison of the models, as well as comparison with fragmentary direct measurements of the atmospheric ionization, validates their applicability for the entire atmosphere and the whole range of the solar activity level variations. The effect of sporadic solar energetic particle events is shown to be limited on the global scale, even for the most severe event, but can be very strong locally in polar regions, affecting the physical-chemical properties of the upper atmosphere, especially at high altitudes. Thus, a new methodology is presented to study cosmic ray induced ionization of the atmosphere in full detail using realistic numerical models calibrated to direct observations.     

Separation of diffracted waves in transversely isotropic media

Imaging diffracted waves can provide useful information about complex subsurface geology and fracture networks. Separation of diffractions from typically more intensive reflected events can be done based on specularity, which measures deviation from Snell’s law. Here, we analyze two formulations of specularity and their applicability to diffraction processing in the presence of anisotropy. We show that the most common definition of specularity, originally introduced for pure modes in isotropic media, remains valid for both pure and converted waves in arbitrarily anisotropic models. The other formulation operates directly with the difference between the slowness projections onto the reflector for the incident and reflected waves. Testing on a VTI (transversely isotropic with a vertical symmetry axis) diffraction ramp model demonstrates that both formulations produce satisfactory results for anisotropic media with appropriate tapering of the specularity gathers. Then separation and imaging of diffractions is performed for the structurally complex VTI Marmousi model. We also analyze the sensitivity of diffractions in the specularity gathers to errors in the symmetry-direction velocity and anellipticity parameter η.     

Volcanic arc of Kamchatka: a province with high-δO magma sources and large-scale O/O depletion of the upper crust

We present the results of a regional study of oxygen and Sr-Nd-Pb isotopes of Pleistocene to Recent arc volcanism in the Kamchatka Peninsula and the Kuriles, with emphasis on the largest caldera-forming centers. The δ¹⁸O values of phenocrysts, in combination with numerical crystallization modeling (MELTS) and experimental fractionation factors, are used to derive best estimates of primary values for δ¹⁸O(magma). Magmatic δ¹⁸O values span 3.5‰ and are correlated with whole-rock Sr-Nd-Pb isotopes and major elements. Our data show that Kamchatka is a region of isotopic diversity with high-δ¹⁸O basaltic magmas (sampling mantle to lower crustal high-δ¹⁸O sources), and low-δ¹⁸O silicic volcanism (sampling low-δ¹⁸O upper crust). Among one hundred Holocene and Late Pleistocene eruptive units from 23 volcanic centers, one half represents low-δ¹⁸O magmas (+4 to 5‰). Most low-δ¹⁸O magmas are voluminous silicic ignimbrites related to large >10 km³ caldera-forming eruptions and subsequent intracaldera lavas and domes: Holocene multi-caldera Ksudach volcano, Karymsky and Kurile Lake-Iliinsky calderas, and Late Pleistocene Maly Semyachik, Akademy Nauk, and Uzon calderas. Low-δ¹⁸O magmas are not found among the less voluminous products of stratovolcano eruptions and these volcanoes do not show drastic changes in δ¹⁸O during their evolution. Additionally, high-δ¹⁸O(magma) of +6.0 to 7.5‰ are found among basalts and basaltic andesites of Bezymianny, Shiveluch, Avachinsky, and Koryaksky volcanoes, and dacites and rhyolites of Opala and Khangar volcanoes (7.1-8.0‰). Phenocrysts in volcanic rocks from the adjacent Kurile Islands (ignimbrites and lavas) define normal-δ¹⁸O magmas. The widespread and volumetric abundance of low-δ¹⁸O magmas in the large landmass of Kamchatka is possibly related to a combination of near-surface volcanic processes, the effects of the last glaciation on high-latitude meteoric waters, and extensive geyser and hydrothermal systems that are matched only by Iceland. Sr and Pb isotopic compositions of normal and low-δ¹⁸O, predominantly silicic, volcanic rocks show negative correlation with δ¹⁸O, similar to the trend in Iceland. This indicates that low-δ¹⁸O volcanic rocks are largely produced by remelting of older, more radiogenic, hydrothermally altered crust that suffered δ¹⁸O-depletion during >2 My-long Pleistocene glaciation. The regionally-distributed high-δ¹⁸O values for basic volcanism (ca. + 6 to +7.5‰) in Kamchatka cannot be solely explained by high-δ¹⁸O slab fluid or melt (± sediment) addition in the mantle, or local subduction of hydrated OIB-type crust of the Hawaii-Emperor chain. Overall, Nd-Pb isotope systematics are MORB-like. Voluminous basic volcanism (in the Central Kamchatka Depression in particular) requires regional, though perhaps patchy, remobilization of thick (30-45 km) Mesozoic-Miocene arc roots, possibly resulting from interaction with hot (ca. 1300°C), wedge-derived normal-δ¹⁸O, low-⁸⁷Sr/⁸⁶Sr basalts and from dehydration melting of lower crustal metabasalts, variably high in δ¹⁸O and ⁸⁷Sr/⁸⁶Sr.     

Source-independent waveform inversion for attenuation estimation in anisotropic media

In previous publications, we presented a waveform-inversion algorithm for attenuation analysis in heterogeneous anisotropic media. However, waveform inversion requires an accurate estimate of the source wavelet, which is often difficult to obtain from field data. To address this problem, here we adopt a source-independent waveform-inversion algorithm that obviates the need for joint estimation of the source signal and attenuation coefficients. The key operations in that algorithm are the convolutions (1) of the observed wavefield with a reference trace from the modelled data and (2) of the modelled wavefield with a reference trace from the observed data. The influence of the source signature on attenuation estimation is mitigated by defining the objective function as the ℓ₂-norm of the difference between the two convolved data sets. The inversion gradients for the medium parameters are similar to those for conventional waveform-inversion techniques, with the exception of the adjoint sources computed by convolution and cross-correlation operations. To make the source-independent inversion methodology more stable in the presence of velocity errors, we combine it with the local-similarity technique. The proposed algorithm is validated using transmission tests for a homogeneous transversely isotropic model with a vertical symmetry axis that contains a Gaussian anomaly in the shear-wave vertical attenuation coefficient. Then the method is applied to the inversion of reflection data for a modified transversely isotropic model from Hess. It should be noted that due to the increased nonlinearity of the inverse problem, the source-independent algorithm requires a more accurate initial model to obtain inversion results comparable to those produced by conventional waveform inversion with the actual wavelet.     

On the development of a magnetic susceptibility-based tracer for aeolian sediment transport research

Aeolian processes – the erosion, transport, and deposition of sediment by wind – play important geomorphological and ecological roles in drylands. These processes are known to impact the spatial patterns of soil, nutrients, plant-available water, and vegetation in many dryland ecosystems. Tracers, such as rare earth elements and stable isotopes have been successfully used to quantify the transport and redistribution of sediment by aeolian processes in these ecosystems. However, many of the existing tracer techniques are labor-intensive and cost-prohibitive, and hence simpler alternative approaches are needed to track aeolian redistribution of sediments. To address this methodological gap, we test the applicability of a novel metal tracer-based methodology for estimating post-fire aeolian sediment redistribution, using spatio-temporal measurements of low-field magnetic susceptibility (MS). We applied magnetic metal tracers on soil microsites beneath shrub vegetation in recently burned and in control treatments in a heterogeneous landscape in the Chihuahuan desert (New Mexico, USA). Our results indicate a spatially homogeneous distribution of the magnetic tracers on the landscape after post-burn wind erosion events. MS decreased after wind erosion events on the burned shrub microsites, indicating that these areas functioned as sediment sources following the wildfire, whereas they are known to be sediment sinks in the undisturbed (e.g. not recently burned) ecosystem. This experiment represents the first step toward the development of a cost-effective and non-destructive tracer-based approach to estimate the transport and redistribution of sediment by aeolian processes. © 2018 John Wiley & Sons, Ltd.     

Oxygen isotope heterogeneity and disequilibria of olivine crystals in large volume Holocene basalts from Iceland: Evidence for magmatic digestion and erosion of Pleistocene hyaloclastites

This work considers petrogenesis of the largest Holocene basaltic fissure eruptions of Iceland, which are also the largest in the world: Laki (1783-84 AD, 15 km³), Eldgjá (934 AD, 18 km³), Veidivötn (900, 1480 AD, multiple eruptions, >2 km³), Núpahraun (ca. 4000 BP, >1 km³) and Thjórsárhraun (ca 8000 BP, >20 km³). We present oxygen isotope laser fluorination analyses of 55 individual and bulk olivine crystals, coexisting individual and bulk plagioclase phenocrysts, and their host basaltic glasses with average precision of better than 0.1‰ (1SD). We also report O isotope analyses of cores and rims of 61 olivine crystals by SIMS with average precision on single spots of 0.24‰ (1SD) in 13 samples coupled with electron microprobe data for major and trace elements in these olivines. Within each individual sample, we have found that basaltic glass is relatively homogeneous with respect to oxygen isotopes, plagioclase phenocrysts exhibit crystal to crystal variability, while individual olivines span from the values in equilibrium with the low-δ¹⁸O matrix glass to those being three permil higher in δ¹⁸O than the equilibrium. Olivine cores with maximum value of 5.2‰ are found in many of these basalts and suggest that the initial magma was equilibrated with normal-δ¹⁸O mantle. No olivines or their intracrystalline domains are found with bulk or spot value higher than those found in MORB olivines. The δ¹⁸O variability of 0.3-3‰ exists for olivine grains from different lavas, and variable core-to-rim oxygen isotopic zoning is present in selected olivine grains. Many olivines in the same sample are not zoned, while a few grains are zoned with respect to oxygen isotopes and exhibit small core-to-core variations in Fe-Mg, Ni, Mn, Ca. Grains that are zoned in both Mg# and δ¹⁸O exhibit positive correlation of these two parameters. Electron microprobe analysis shows that most olivines equilibrated with the transporting melt, and thin Fe-richer rim is present around many grains, regardless of the degree of olivine-melt oxygen isotope disequilibrium.The preservation of isotopic and compositional zoning in selected grains, and subtle to severe Δ¹⁸O (melt-olivine) and Δ¹⁸O (plagioclase-olivine) disequilibria suggests rather short crystal residence times of years to centuries. Synglacially-altered upper crustal, tufaceous hyaloclastites of Pleistocene age serve as a viable source for low-δ¹⁸O values in Holocene basalts through assimilation, mechanical and thermal erosion, and devolatilization of stoped blocks. Cumulates formed in response to cooling during assimilation, and xenocrysts derived from hyaloclastites, contribute to the diverse δ¹⁸O crystalline cargo. The magma plumbing systems under each fissure are likely to include a network of interconnected dikes and sills with high magma flow rates that contribute to the efficacy of magmatic erosion of large quantities (10-60% mass) of hyaloclastites required by isotopic mass balance.Olivine diversity and the pervasive lack of phenocryst-melt oxygen isotopic equilibrium suggest that a common approach of analyzing bulk olivine for oxygen isotopes, as a proxy for the basaltic melt or to infer mantle δ¹⁸O value, needs to proceed with caution. The best approach is to analyze olivine crystals individually and demonstrate their equilibrium with matrix.     

Sources of unique rhenium enrichment in fumaroles and sulphides at Kudryavy volcano

Rhenium (Re) is one of the least abundant elements in Earth, averaging 0.28 ppb in the primitive mantle. The unique occurrence of rheniite ReS₂ (74.5 wt% of Re) in Kudryavy volcano precipitates raises questions about recycling of Re-rich reservoirs within the Kurile-Kamchatka volcanic Island arc setting. The sources of this unique Re enrichment have been inferred from studies of Re-Os isotope systematic and trace elements in volcanic gases, sulphide precipitates and host volcanic rocks. The fumarolic gas condensates are enriched in hydrophile trace elements relative to fluid-immobile elements and exhibit high Ba/Nb (133-204), Rb/Y (16-406) and Th/Zr (0.01-0.25) ratios. They are characterised by high Re (7-210 ppb) and Os abundances (0.4-0.9 ppb), with ¹⁸⁷Os/¹⁸⁸Os ratios in a range 0.122-0.152. This Os isotopic compositional range is similar to that of the peridotite xenoliths from the metasomatised mantle wedge above the subducted Pacific plate, the radiogenic isotopic signature of which is probably due to radiogenic addition from a slab-derived fluid.Re- and Os-rich sulphide and oxide minerals precipitate from volcanic gases within fumarolic fields. Molybdenite (MoS₂), powellite (CaMoO₄) and cannizzarite (Pb₄Bi₆S₁₃) contain 1.5-1.7 wt%, 10 ppm, and 65-252 ppb of Re, respectively. Both molybdenite and rheniite contain normal Os concentrations, with total Os abundances in a range from 0.6 to 3.1 ppm for molybdenite, and 2.3-24.3 ppb for the rheniite samples. Repeated analyses of osmium isotope ratios for two rheniite samples form a best-fit line with an initial ¹⁸⁷Os/¹⁸⁸Os ratio of 0.32 ± 0.15 and an age of 79 ± 11 yr, which is the youngest age ever measured in natural samples. The high Re contents in molybdenite and rheniite led to high radiogenic ¹⁸⁷Os values, even in the limited period of time, with ¹⁸⁷Os/¹⁸⁸Os ratios up to 3.3 for molybdenite and up to 4.4 for rheniite.The Os isotopic compositions of andesite-basaltic rocks from the Kudryavy volcano (¹⁸⁷Os/¹⁸⁸Os up to 0.326) are more radiogenic than those of residual peridotites and fumarolic gas condensates that are mainly constituted from magmatic vapor. Such radiogenic values can be attributed either to the addition of a radiogenic Os-rich subduction component to the depleted mantle, or to the assimilation of older dacitic caldera walls (¹⁸⁷Os/¹⁸⁸Os = 0.6) during arc magma ascent and emplacement. The latter hypothesis is supported by the correlation between ¹⁸⁷Os/¹⁸⁸Os ratio and indicators of fractionation such as MgO or Ni, and by low contents of potentially hydrophile trace elements such as Ba, Rb and Th relative to fluid-immobile elements such as Nb, Zr and Y. The high Re flux in the Kudryavy volcano (estimated at ∼46 kg/yr) can be explained by remobilisation of Re by Cl-rich water from an underplated mantle wedge and subducted organic-rich sediments of the Pacific plate.