Spatiotemporal variability in hydrochemistry of shallow groundwater in a small pre‐alpine catchment: The importance of landscape elements

Topography and landscape characteristics affect the storage and release of water and, thus, groundwater dynamics and chemistry. Quantification of catchment scale variability in groundwater chemistry and groundwater dynamics may therefore help to delineate different groundwater types and improve our understanding of which parts of the catchment contribute to streamflow. We sampled shallow groundwater from 34 to 47 wells and streamflow at seven locations in a 20‐ha steep mountainous catchment in the Swiss pre‐Alps, during nine baseflow snapshot campaigns. The spatial variability in electrical conductivity, stable water isotopic composition, and major and trace ion concentrations was large and for almost all parameters larger than the temporal variability. Concentrations of copper, zinc, and lead were highest at sites that were relatively dry, whereas concentrations of manganese and iron were highest at sites that had persistent shallow groundwater levels. The major cation and anion concentrations were only weakly correlated to individual topographic or hydrodynamic characteristics. However, we could distinguish four shallow groundwater types based on differences from the catchment average concentrations: riparian zone‐like groundwater, hillslopes and areas with small upslope contributing areas, deeper groundwater, and sites characterized by high magnesium and sulfate concentrations that likely reflect different bedrock material. Baseflow was not an equal mixture of the different groundwater types. For the majority of the campaigns, baseflow chemistry most strongly resembled riparian‐like groundwater for all but one subcatchment. However, the similarity to the hillslope‐type groundwater was larger shortly after snowmelt, reflecting differences in hydrologic connectivity. We expect that similar groundwater types can be found in other catchments with steep hillslopes and wet areas with shallow groundwater levels and recommend sampling of groundwater from all landscape elements to understand groundwater chemistry and groundwater contributions to streamflow.     

Geochronological and isotopic records of crustal storage and assimilation in the Wolverine Creek–Conant Creek system,Heise eruptive centre,Snake River Plain

Understanding the processes of differentiation of the Yellowstone–Snake River Plain (YSRP) rhyolites is typically impeded by the apparent lack of erupted intermediate compositions as well as the complex nature of their shallow interaction with the surrounding crust responsible for their typically low O isotopic ratios. A pair of normal-δ¹⁸O rhyolitic eruptions from the Heise eruptive centre in eastern Idaho, the Wolverine Creek Tuff and the Conant Creek Tuff, represent unique magmatic products of the Yellowstone hotspot preserving abundant vestiges of the intermediate differentiation steps leading to rhyolite generation. We address both shallow and deep processes of magma generation and storage in the two units by combining high-precision ID–TIMS U–Pb zircon geochronology, trace element, O and Hf isotopic studies of zircon, and Sr isotopic analyses of individual high-Mg# pyroxenes inherited from lower- to mid-crustal differentiation stages. The zircon geochronology confirms the derivation of both tuffs from the same rhyolitic magma reservoir erupted at 5.5941 ± 0.0097 Ma, preceded by at least 92 ± 14 ky of continuous or intermittent zircon saturation approximating the length of pre-eruptive magma accumulation in the upper crust. Some low-Mg# pyroxenes enclosing zircons predate the eruption by at least 45 ± 27 ky, illustrating the co-crystallisation of major and accessory phases in the near-liquidus rhyolitic melts of the YSRP over a significant period of time. Coeval zircon crystals are isotopically heterogeneous (two populations at εHf ~?5 and ?13), requiring the assembly of isotopically distinct melt pockets directly prior to, or during, the eruption. The primitive Mg# 60–90 pyroxenes are out of isotopic equilibrium with the host rhyolitic melt (⁸⁷Sr/⁸⁶Sr_i = 0.70889), covering a range of ⁸⁷Sr/⁸⁶Sr_i = 0.70705–0.70883 corresponding to ratios typical of the most radiogenic YSRP basalts to the least radiogenic YSRP rhyolites. Together with the low εHf in zircon, the Sr isotopic ratios illustrate limited assimilation dominated by radiogenic Archean crustal source materials incorporated into variably evolved YSRP melts as they progress towards rhyolitic compositions by assimilation–fractional crystallisation.     

Tectonic significance of Triassic mafic rocks in the June Complex,Sanandaj–Sirjan zone,Iran

This study concentrates on the petrological and geochemical investigation of mafic rocks embedded within the voluminous Triassic June Complex of the central Sanandaj–Sirjan zone (Iran), which are crucial to reconstruct the geodynamics of the Neotethyan passive margin. The Triassic mafic rocks are alkaline to sub-alkaline basalts, containing 43.36–49.09 wt% SiO₂, 5.19–20.61 wt% MgO and 0.66–4.59 wt% total alkalis. Based on MgO concentrations, the mafic rocks fall into two groups: cumulates (Mg# = 51.61–58.94) and isotropic basaltic liquids (Mg# = 24.54–42.66). In all samples, the chondrite-normalized REE patterns show enrichment of light REEs with variable (La/Yb)_N ratios ranging from 2.48 to 9.00, which confirm their amalgamated OIB-like and E-MORB-like signatures. Enrichment in large-ion lithophile elements and depletion in high field strength elements (HFSE) relative to the primitive mantle further support this interpretation. No samples point to crustal contamination, all having undergone fractionation of olivine + clinopyroxene + plagioclase. Nevertheless, elemental data suggest that the substantial variations in (La/Sm)_PM and Zr/Nb ratios can be explained by variable degrees of partial melting rather than fractional crystallization from a common parental magma. The high (Nb/Yb)_PM ratio in the alkaline mafic rocks points to the mixing of magmas from enriched and depleted mantle sources. Abundant OIB alkaline basalts and rare E-MORB appear to be linked to the drifting stage on the northern passive margin of the Neotethys Ocean.     

Evaluating igneous sources of the Taveyannaz formation in the Central Alps by detrital zircon U–Pb age dating and geochemistry

Late Palaeogene syn-tectonic volcanic products have been found in the Northern Alpine foreland basin and in the South Alpine hemipelagic basin. The source of abundant volcanic fragments is still in debate. We analyzed the geochronology and geochemistry of detrital zircons, and evaluated their temporal and genetic relationships with potential volcanic sources. The study shows that the detrital zircon U–Pb age patterns have two major age groups: a dominance (ca. 90%) of pre-Alpine zircons was found, as commonly observed in other Alpine flysch formations. These zircons apparently derived from erosion of the early Alpine nappe stack in South Alpine and Austroalpine units. Furthermore, a few Neo-Alpine zircons (ca. 10%) have ages ranging from Late Eocene to Early Oligocene (~ 41–29 Ma). Both source materials were mixed during long riverine transport to the basin margins before being re-deposited by gravity flows. These Palaeogene ages match with the activity of Peri-Adriatic magmatism, including the Biella volcanic suite as well as the Northern Adamello and Bergell intrusions. The values of REE and ¹⁷⁶Hf/¹⁷⁷Hf_(t) ratios of the Alpine detrital zircons are in line with the magmatic signatures. We observe an in time and space variable supply of syn-sedimentary zircons. From late Middle Eocene to Late Eocene, basin influx into the South Alpine and Glarus (A) basins from the Northern Adamello source is documented. At about 34 Ma, a complete reorganisation is recorded by (1) input of Bergell sources into the later Glarus (B) basin, and (2) the coeval volcaniclastic supply of the Haute-Savoie basin from the Biella magmatic system. The Adamello source vanished in the foreland basin. The marked modification of the basin sources at ~ 34 Ma is interpreted to be initiated by a northwestern shift of the early Alpine drainage divide into the position of the modern Insubric Line.     

An experimental study of Fe–Ni exchange between sulfide melt and olivine at upper mantle conditions: implications for mantle sulfide compositions and phase equilibria

The behavior of nickel in the Earth’s mantle is controlled by sulfide melt–olivine reaction. Prior to this study, experiments were carried out at low pressures with narrow range of Ni/Fe in sulfide melt. As the mantle becomes more reduced with depth, experiments at comparable conditions provide an assessment of the effect of pressure at low-oxygen fugacity conditions. In this study, we constrain the Fe–Ni composition of molten sulfide in the Earth’s upper mantle via sulfide melt–olivine reaction experiments at 2 GPa, 1200 and 1400 °C, with sulfide melt \(X_{{{\text{Ni}}}}^{{{\text{Sulfide}}}}=\frac{{{\text{Ni}}}}{{{\text{Ni}}+{\text{Fe}}}}\) (atomic ratio) ranging from 0 to 0.94. To verify the approach to equilibrium and to explore the effect of \({f_{{{\text{O}}_{\text{2}}}}}\) on Fe–Ni exchange between phases, four different suites of experiments were conducted, varying in their experimental geometry and initial composition. Effects of Ni secondary fluorescence on olivine analyses were corrected using the PENELOPE algorithm (Baró et al., Nucl Instrum Methods Phys Res B 100:31–46, 1995), “zero time” experiments, and measurements before and after dissolution of surrounding sulfides. Oxygen fugacities in the experiments, estimated from the measured O contents of sulfide melts and from the compositions of coexisting olivines, were 3.0?±?1.0 log units more reduced than the fayalite–magnetite-quartz (FMQ) buffer (suite 1, 2 and 3), and FMQ ??1 or more oxidized (suite 4). For the reduced (suites 1–3) experiments, Fe–Ni distribution coefficients \(K_{{\text{D}}}^{{}}=\frac{{(X_{{{\text{Ni}}}}^{{{\text{sulfide}}}}/X_{{{\text{Fe}}}}^{{{\text{sulfide}}}})}}{{(X_{{{\text{Ni}}}}^{{{\text{olivine}}}}/X_{{{\text{Fe}}}}^{{{\text{olivine}}}})}}\) are small, averaging 10.0?±?5.7, with little variation as a function of total Ni content. More oxidized experiments (suite 4) give larger values of K_D (21.1–25.2). Compared to previous determinations at 100 kPa, values of K_D from this study are chiefly lower, in large part owing to the more reduced conditions of the experiments. The observed difference does not seem attributable to differences in temperature and pressure between experimental studies. It may be related in part to the effects of metal/sulfur ratio in sulfide melt. Application of these results to the composition of molten sulfide in peridotite indicates that compositions are intermediate in composition (\(X_{{{\text{Ni}}}}^{{{\text{sulfide}}}}\)?~?0.4–0.6) in the shallow mantle at 50 km, becomes more Ni rich with depth as the O content of the melt diminishes, reaching a maximum (0.6–0.7) at depths near 80–120 km, and then becomes more Fe rich in the deeper mantle where conditions are more reduced, approaching (\(X_{{{\text{Ni}}}}^{{{\text{sulfide}}}}\)?~?0.28)?>?140 km depth. Because Ni-rich sulfide in the shallow upper mantle melts at lower temperature than more Fe-rich compositions, mantle sulfide is likely molten in much of the deep continental lithosphere, including regions of diamond formation.     

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.