Crustal contamination and mantle source characteristics in continental intra-plate volcanic rocks: Pb, Hf and Os isotopes from central European volcanic province basalts

We report new Os-Pb-Hf isotope data for a suite of alkaline to basaltic (nephelinites, basanites, olivine tholeiites to quartz-tholeiites) lavas from the Miocene Vogelsberg (Germany), the largest of the rift-related continental volcanic complexes of the Central European Volcanic Province (CEVP). ¹⁸⁷Os/¹⁸⁸Os in primitive (high-MgO) alkaline lavas show a much wider range than has been observed in alkaline basalts and peridotite xenoliths from elsewhere in the CEVP, from ratios similar to those in modern MORB and OIB (0.1260-0.1451; 58.9-168 ppt Os) to more radiogenic ratios (0.1908 and 0.2197; 27.6-15.1 ppt Os). Radiogenic Os is associated with high ε_Hf and ε_Nd, low ⁸⁷Sr/⁸⁶Sr and does not correlate with Mg^∗ or incompatible trace elements (e.g. Ce/Pb), suggesting the presence of a radiogenic endmember in the mantle rather than crustal contamination as the source of radiogenic Os. This contrasts with another high-Mg alkaline lava characterized by highly radiogenic ¹⁸⁷Os/¹⁸⁸Os (0.4344, 10.3 ppt Os), lower ε_Hf and ε_Nd, higher ⁸⁷Sr/⁸⁶Sr, and Pb isotope signatures than the other alkaline lavas with similar trace element composition suggestive of contamination with crustal material. Hafnium (ε_Hf: +8.9 to +5.0) and Pb isotope compositions (²⁰⁶Pb/²⁰⁴Pb: 19.10-19.61; ²⁰⁷Pb/²⁰⁴Pb: 15.56-15.60) of the alkaline rocks fall within the range of enriched MORB and some OIB. The Vogelsberg tholeiites show even more diverse ¹⁸⁷Os/¹⁸⁸Os, ranging from 0.1487 in Os-rich olivine tholeiite (31.7 ppt) to ratios as high as 0.7526 in other olivine-tholeiites and in quartz-tholeiites with lower Os concentrations (10.3-2.0 ppt). Low-¹⁸⁷Os/¹⁸⁸Os tholeiites show Pb-Hf isotope ratios (²⁰⁶Pb/²⁰⁴Pb:18.81; ²⁰⁷Pb/²⁰⁴Pb: 15.61; ε_Hf: +2.7) that are distinct from those in alkaline lavas with similar ¹⁸⁷Os/¹⁸⁸Os and originate from a different mantle source. By contrast, the combination of radiogenic Os and low ²⁰⁶Pb/²⁰⁴Pb and ε_Hf in the other tholeiites probably reflects crustal contamination.The association at Vogelsberg of primitive alkaline and tholeiitic lavas with a range of MORB- to OIB-like Os-Pb-Hf-Nd-Sr isotopic characteristics requires at least two asthenospheric magma sources. This is consistent with trace element modelling which suggests that the alkaline and tholeiitic parent magmas represent mixtures of melts from garnet and spinel peridotite sources (both with amphibole), implying an origin of the magmas in the garnet peridotite-spinel peridotite transition zone, probably at the asthenosphere-lithosphere interface. We propose that uncontaminated Vogelsberg lavas originated in ‘metasomatized’ mantle, involving a 3-stage model: (1) early carbonatite metasomatism several 10-100 Ma before the melting event (2) deposition of low-degree asthenospheric melts from carbonated peridotite at the lithosphere-asthenosphere thermal boundary produces hydrous amphibole-bearing veins or patches, and (3) remobilization of this modified lithospheric mantle into other asthenospheric melts passing through the same area later. In keeping with ‘metasomatized’ mantle models for other continental basalt provinces, we envisage that stage (2) is short-lived (few Ma), thus producing a prominent lithospheric trace element signature without changing the asthenospheric isotopic signatures. Models of this type can explain the peculiar mix of lithospheric (prominent depletions of Rb and K) and asthenospheric (OIB-like high ¹⁸⁷Os/¹⁸⁸Os, ¹⁴³Nd/¹⁴⁴Nd and ¹⁷⁶Hf/¹⁷⁷Hf) signatures observed in the Vogelsberg and many other continental basalt suites.     

Quantifying two-way influences between the Arctic and mid-latitudes through regionally increased CO2 concentrations in coupled climate simulations

Climate Dynamics - In which direction is the influence larger: from the Arctic to the mid-latitudes or vice versa? To answer this question, CO2 concentrations have been regionally increased in...     

The role of mantle-hybridization and crustal contamination in the petrogenesis of lithospheric mantle-derived alkaline rocks: constraints from Os and Hf isotopes

The Rhön area as part of the Central European Volcanic Province (CEVP) hosts an unusual suite of Tertiary 24-Ma old hornblende-bearing alkaline basalts that provide insights into melting and fractionation processes within the lithospheric mantle. These chemically primitive to slightly evolved and isotopically (Sr, Nd, Pb) depleted basalts have slightly lower Hf isotopic compositions than respective other CEVP basalts and Os isotope compositions more radiogenic than commonly observed for continental intraplate alkaline basalts. These highly radiogenic initial ¹⁸⁷Os/¹⁸⁸Os ratios (0.268–0.892) together with their respective Sr–Nd–Pb isotopic compositions are unlikely to result from crustal contamination alone, although a lack of Os data for lower crustal rocks from the area and limited data for CEVP basalts or mantle xenoliths preclude a detailed evaluation. Similarly, melting of the same metasomatized subcontinental lithospheric mantle as inferred for other CEVP basalts alone is also unlikely, based on only moderately radiogenic Os isotope compositions obtained for upper mantle xenoliths from elsewhere in the province. Another explanation for the combined Nd, Sr and Os isotope data is that the lavas gained their highly radiogenic Os isotope composition through a mantle “hybridization”, metasomatism process. This model involves a mafic lithospheric component, such as an intrusion of a sublithospheric primary alkaline melt or a melt derived from subducted oceanic material, sometime in the past into the lithospheric mantle where it metasomatized the ambient mantle. Later at 24 Ma, thermal perturbations during rifting forced the isotopically evolved parts of the mantle together with the peridotitic ambient mantle to melt. This yielded a package of melts with highly correlated Re/Os ratios and radiogenic Os isotope compositions. Subsequent movement through the crust may have further altered the Os isotope composition although this effect is probably minor for the majority of the samples based on radiogenic Nd and unradiogenic Sr isotope composition of the lavas. If the radiogenic Os isotope composition can be explained by a mantle-hybridization and metasomatism model, the isotopic compositions of the hornblende basalts can be satisfied by ca. 5–25% addition of the mafic lithospheric component to an asthenospheric alkaline magma. Although a lack of isotope data for all required endmembers make this model somewhat speculative, the results show that the Re–Os isotope system in continental basalts is able to distinguish between crustal contamination and derivation of continental alkaline lavas from isotopically evolved peridotitic lithosphere that was contaminated by mafic material in the past and later remelted during rifting. The Hf isotopic compositions are slightly less radiogenic than in other alkaline basalts from the province and indicate the derivation of the lavas from low Lu–Hf parts of the lithospheric mantle. The new Os and Hf isotope data constrain a new light of the nature of such metasomatizing agents, at least for these particular rocks, which represent within the particular volcanic complex the first product of the volcanism.     

Assessment of co-seismic landslide susceptibility using LR and ANCOVA in Barpak region,Nepal

Nepal was affected by a catastrophic earthquake with Mw 7.8 on 25th April, 2015 with its epicenter in the central part of Barpak village. A number of co-seismic landslides were triggered by the main shock of the event and associated aftershocks. Due to the rugged topography and vicinity of the main shock, the village was extremely affected by co-seismic landslides. In total, 59 landslides were identified using Google Earth and were verified during the field survey in Barpak village. Furthermore, 11 conditioning factors, including Peak ground acceleration (PGA), epicenter proximity, fault proximity, geology, slope, elevation, plan curvature, profile curvature, topographic wetness index, drainage proximity and the sediment transport index were selected as independent variables for analysis. In this study, logistic regression (LR) and analysis of covariance (ANCOVA) models were used and their performance was assessed. Finally, the landslide susceptibility classes were produced and an evaluation of models was done by using receiver operating characteristic curves. The area under the curve for LR and ANCOVA showed 85.38 and 78.4% accuracy, respectively. Based on the overall assessments, the LR model was more accurate than the ANCOVA model for co-seismic landslide prediction in the study area. The result of this study can be used to mitigate landslide-induced hazards and for land-use planning.     

On the soil water characteristic curves of poorly graded granular materials in aqueous polymer solutions

Tempe cell with a self-developed horizontal hanging column attachment was used to measure the soil water characteristic curves (SWCCs) of granular materials initially saturated with either water or polymer solutions. For SWCCs of six poorly graded granular materials with d ₅₀ ranging from 0.04 to 0.7 mm in water, it was found that (1) as grain size decreases, air entry value increases, and matric suction (ψ) for the funicular and pendular regimes increases and that (2) steep desaturation curves over narrow ψ range in the funicular regime were observed. Air entry values obtained from the fitting parameter 1/α in van Genuchten SWCC equation fall in the boundaries calculated from the pore throat sizes in both simple cubic and face-centered cubic packings. A toroidal meniscus water model, which incorporates the measured surface tension and contact angle values between aqueous solutions and solid surfaces, was proposed for the SWCC in pendular regime and was compared to well-received numerical methods. This toroid model successfully depicts SWCC of poorly graded granular materials in water. However, SWCCs predicted by this toroid model underestimate the degree of saturation in the pendular regime for Ottawa 20–30 sands in polymer solutions. Herschel–Bulkley fluid, which is a type of non-Newtonian fluids, is postulated to increase the ψ needed to drain the polymer solution due to the nonzero shear stress intercept. In addition, it is also postulated by scanning electron microscopy and optical confocal imaging results that the rough surfaces of Ottawa 20–30 sand, which have many micron-sized “kinks”, together with the possible chemical attractions, help retain the polymer solutions on the solid surfaces, or water film.     

Large-scale rainfall–runoff–inundation modeling for upper Citarum River watershed,Indonesia

The Citarum River is one of the strategic rivers in West Java, Indonesia. Its total watershed area is approximately 1800 km². Almost every year, the overflow from the Citarum River causes the inundation of most of the upper Citarum River watershed. To prevent and mitigate flood damage, it is necessary to understand the flooding characteristics. The region, however, suffers from a lack of observational data. Therefore, to analyze the inundation caused by flooding in the upper Citarum River watershed, a rainfall–runoff–inundation (RRI) model was employed. It used the following multiple satellite-derived datasets as input data as well as for model verification: Global Satellite Mapping of Precipitation, Hydrological data and maps based on Shuttle elevation Derivatives at multiple scales, Global Mosaics of the standard MODIS land cover type data product, and Landsat 7 satellite images. Parameter calibration was performed using a Monte Carlo simulation. The simulation was performed for February 2010. The results of this study show that the RRI model identifies inundation areas in large-scale river watersheds more effectively when using multiple satellite-derived datasets compared with the observed inundation map obtained from JICA in 2010 and Landsat 7 images. The model results can be improved if high-quality observed rainfall data, topographic data, and river cross-sectional data are available.     

Landslide susceptibility mapping using random forest and boosted tree models in Pyeong-Chang,Korea

Landslides susceptibility maps were constructed in the Pyeong-Chang area, Korea, using the Random Forest and Boosted Tree models. Landslide locations were randomly selected in a 50/50 ratio for training and validation of the models. Seventeen landslide-related factors were extracted and constructed in a spatial database. The relationships between the observed landslide locations and these factors were identified by using the two models. The models were used to generate a landslide susceptibility map and the importance of the factors was calculated. Finally, the landslide susceptibility maps were validated. Finally, landslide susceptibility maps were generated. For the Random Forest model, the validation accuracy in regression and classification algorithms showed 79.34 and 79.18%, respectively, and for the Boosted Tree model, these were 84.87 and 85.98%, respectively. The two models showed satisfactory accuracies, and the Boosted Tree model showed better results than the Random Forest model.     

Non-parametric simulations-based conditional stochastic predictions of geologic heterogeneities and leakage potentials for hypothetical CO2 sequestration sites

The present study focuses on understanding the leakage potentials of the stored supercritical CO₂ plume through caprocks generated in geostatistically created heterogeneous media. For this purpose, two hypothetical cases with different geostatistical features were developed, and two conditional geostatistical simulation models (i.e., sequential indicator simulation or SISIM and generalized coupled Markov chain or GCMC) were applied for the stochastic characterizations of the heterogeneities. Then, predictive CO₂ plume migration simulations based on stochastic realizations were performed and summarized. In the geostatistical simulations, the results from the GCMC model showed better performance than those of the SISIM model for the strongly non-stationary case, while SISIM models showed reasonable performance for the weakly non-stationary case in terms of low-permeability lenses characterization. In the subsequent predictive simulations of CO₂ plume migration, the observations in the geostatistical simulations were confirmed and the GCMC-based predictions showed underestimations in CO₂ leakage in the stationary case, while the SISIM-based predictions showed considerable overestimations in the non-stationary case. The overall results suggest that: (1) proper characterization of low-permeability layering is significantly important in the prediction of CO₂ plume behavior, especially for the leakage potential of CO₂ and (2) appropriate geostatistical techniques must be selectively employed considering the degree of stationarity of the targeting fields to minimize the uncertainties in the predictions.     

Sensitivity analysis of the Ångstrom exponent for multimodal aerosol size distributions

This study characterizes the Ångstrom exponent for polydispersed aerosol size distributions. Under the assumption of a lognormal size distribution, the dependence of Ångstrom exponent on the size distribution and the refractive index with varying real and imaginary parts are determined. Further, the influence of coarse mode particles on the Ångstrom exponent is investigated quantitatively. The results show that the nuclei mode has less influence under the simulation conditions considered in this study. It is also shown that the refractive index is an important factor influencing the Ångstrom exponent. The effect of the coarse mode on the Ångstrom exponent computed with different aerosol number concentrations and as a function of a geometric standard deviation and a geometric mean diameter is tested. It is shown that the coarse mode is crucial for determining the Ångstrom exponent.