Receiver function study in northern Sumatra and the Malaysian peninsula

In this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities (V _S) near the surface in the range of 3.4–3.6 km s^− 1 attributed to crystalline rocks and 3.6–4.0 km s^− 1 in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30–34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low V _S of 2.6–2.9 km s^− 1 indicate sediment layers. High V _S of 4.2 km s^− 1 are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in V _S between 10 and 20 km depth. Within the subducting slab V _S ≈ 4.7 km s^− 1. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found.     

Fluid-induced mineral composition adjustments during exhumation: the case of Alpine stilbite

Stilbite is locally present as a very late mineral on fractures and fissures of granitic basement in the Central Swiss Alps. Stilbite samples from the Gotthard rail base tunnel provide evidence that they originally formed as a K-absent variety at depth. However, all stilbite samples from surface outcrops above the tunnel display significant potassium concentrations. Interestingly, water from fractures in the tunnel (at 50 °C) is oversaturated with respect to stilbite and essentially potassium-free whereas waters from high-Alpine brooks above the tunnel (and at other high-Alpine areas) have unusually high K/Na ratios. The data suggest that stilbite that may actively form on fissures at tunnel level as a K-absent variety by precipitation from water. Older stilbite that originally formed as coatings on fracture walls was gradually exhumed and uplifted and finally reached the today’s erosion surface about 2,000 m above the tunnel. However, the stilbite reaches the erosion surface as a K-rich variety as a result of interaction of the original low-K stilbite with surface water and near-surface groundwater. This leads to the conclusion that minerals once formed at depth may significantly change their composition once they reach the ground water zone on their way to the erosion surface. In the case of the stilbite, if surface outcrops would have been the only source of samples and data, the K-rich composition could have been mistaken for the composition of the mineral when it formed, which is not the case. Late-stage compositional readjustments may be difficult to discern in samples from surface outcrops. The provided data show that original mineral compositions may be adjusted by late-stage water–rock interaction in a highly selective way.     

Chemical evolution of thermal waters and changes in the hydrothermal system of Papandayan volcano (West Java,Indonesia) after the November 2002 eruption

Papandayan is a stratovolcano situated in West Java, Indonesia. Since the last magmatic eruption in 1772, only few hydrothermal explosions have occurred. An explosive eruption occurred in November 2002 and ejected ash and altered rocks. The altered rocks show that an advanced argillic alteration took place in the hydrothermal system by interaction between acid fluids and rocks. Four zones of alteration have been defined and are limited in extension and shape along faults or across permeable structures at different levels beneath the active crater of the volcano.     

Experimental impact cratering: A summary of the major results of the MEMIN research unit

This paper reviews major findings of the Multidisciplinary Experimental and Modeling Impact Crater Research Network (MEMIN). MEMIN is a consortium, funded from 2009 till 2017 by the German Research Foundation, and is aimed at investigating impact cratering processes by experimental and modeling approaches. The vision of this network has been to comprehensively quantify impact processes by conducting a strictly controlled experimental campaign at the laboratory scale, together with a multidisciplinary analytical approach. Central to MEMIN has been the use of powerful two-stage light-gas accelerators capable of producing impact craters in the decimeter size range in solid rocks that allowed detailed spatial analyses of petrophysical, structural, and geochemical changes in target rocks and ejecta. In addition, explosive setups, membrane-driven diamond anvil cells, as well as laser irradiation and split Hopkinson pressure bar technologies have been used to study the response of minerals and rocks to shock and dynamic loading as well as high-temperature conditions. We used Seeberger sandstone, Taunus quartzite, Carrara marble, and Weibern tuff as major target rock types. In concert with the experiments we conducted mesoscale numerical simulations of shock wave propagation in heterogeneous rocks resolving the complex response of grains and pores to compressive, shear, and tensile loading and macroscale modeling of crater formation and fracturing. Major results comprise (1) projectile–target interaction, (2) various aspects of shock metamorphism with special focus on low shock pressures and effects of target porosity and water saturation, (3) crater morphologies and cratering efficiencies in various nonporous and porous lithologies, (4) in situ target damage, (5) ejecta dynamics, and (6) geophysical survey of experimental craters.     

Uncertainty of ENSO-amplitude projections in CMIP5 and CMIP6 models

There is a long-standing debate on how the El Niño/Southern Oscillation (ENSO) amplitude may change during the twenty-first century in response to global warming. Here we identify the sources of uncertainty in the ENSO amplitude projections in models participating in the Coupled Model Intercomparison Phase 5 (CMIP5) and Phase 6 (CMIP6), and quantify scenario uncertainty, model uncertainty and uncertainty due to internal variability. The model projections exhibit a large spread, ranging from increasing standard deviation of up to 0.6 °C to diminishing standard deviation of up to − 0.4 °C by the end of the twenty-first century. The ensemble-mean ENSO amplitude change is close to zero. Internal variability is the main contributor to the uncertainty during the first three decades; model uncertainty dominates thereafter, while scenario uncertainty is relatively small throughout the twenty-first century. The total uncertainty increases from CMIP5 to CMIP6: while model uncertainty is reduced, scenario uncertainty is considerably increased. The models with “realistic” ENSO dynamics have been analyzed separately and categorized into models with too small, moderate and too large ENSO amplitude in comparison to instrumental observations. The smallest uncertainties are observed in the sub-ensemble exhibiting realistic ENSO dynamics and moderate ENSO amplitude. However, the global warming signal in ENSO-amplitude change is undetectable in all sub-ensembles. The zonal wind-SST feedback is identified as an important factor determining ENSO amplitude change: global warming signal in ENSO amplitude and zonal wind-SST feedback strength are highly correlated across the CMIP5 and CMIP6 models.

     

Field application of calcite Dispersed Alkaline Substrate (calcite-DAS) for passive treatment of acid mine drainage with high Al and metal concentrations

Passive treatment systems are widely used for remediation of acid mine drainage (AMD), but existing designs are prone to clogging or loss of reactivity due to Al- and Fe-precipitates when treating water with high Al and heavy metal concentrations. Dispersed alkaline substrate (DAS) mixed from a fine-grained alkaline reagent (e.g. calcite sand) and a coarse inert matrix (e.g. wood chips) had shown high reactivity and good hydraulic properties in previous laboratory column tests. In the present study, DAS was tested at pilot field scale in the Iberian Pyrite Belt (SW Spain) on metal mine drainage with pH near 3.3, net acidity 1400–1650 mg/L as CaCO₃, and mean concentrations of 317 mg/L Fe (95% Fe(II)), 311 mg/L Zn, 74 mg/L Al, 20 mg/L Mn, and 1.5–0.1 mg/L Cu, Co, Ni, Cd, As and Pb. The DAS-tank removed an average of 870 mg/L net acidity as CaCO₃ (56% of inflow), 25% Fe, 93% Al, 5% Zn, 95% Cu, 99% As, 98% Pb, and 14% Cd, but no Mn, Ni or Co. Average gross drain pipe alkalinity was 181 mg/L as CaCO₃, which increased total Fe removal to 153 mg/L (48%) in subsequent sedimentation ponds. Unfortunately, the tank suffered clogging problems due to the formation of a hardpan of Al-rich precipitates. DAS lifetime could probably be increased by lowering Al-loads.     

Dynamics of snow ablation in a small Alpine catchment observed by repeated terrestrial laser scans

The spatio‐temporal distribution of snow in a catchment during ablation reflects changes in the total amount of snow water equivalent and is thus a key parameter for the estimation of melt water run‐off. This study explores possible rules behind the spatial variability of snow depth during the ablation season in a small Alpine catchment with complex topography. The snow depth observations are based on more than 160 000 terrestrial laser scanner data points with a spatial resolution of 1 m, which were obtained from 11 scanning campaigns of two consecutive ablation seasons. The analysis suggests that for estimating cumulative snow melt dynamics from the catchment investigated, assessing the initial snow distribution prior to the melt season is more important than addressing spatial differences in the melt behaviour. Snow volume and snow‐covered area could be predicted well using a conceptual melt model assuming spatially uniform melt rates. However, accurate results were only obtained if the model was initialized with a pre‐melt snow distribution that reflected measured mean and standard deviation. Using stratified melt rates on the other hand did not improve the model results. At least for sites with similar meteorological and topographical conditions, the model approach presented here comprises an efficient way to estimate snow depletion dynamics, especially if persistent snow accumulation pattern between years facilitate the characterization of the initial snow distribution prior to the melt. Copyright © 2011 John Wiley & Sons, Ltd.