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.     

Instruments, Detectors and the Future of Astronomy with Large Ground Based Telescopes

Results of a survey of instrumentation and detector systems, either currently deployed or planned for use at telescopes larger than 3.5 m, in ground based observatories world-wide, are presented. This survey revealed a number of instrumentation design trends at optical, near, and mid-infrared wavelengths. Some of the most prominent trends include the development of vastly larger optical detector systems (> 10⁹ pixels) than anything built to date, and the frequent use of mosaics of near-infrared detectors – something that was quite rare only a decade ago in astronomy. Some future science applications for detectors are then explored, in an attempt to build a bridge between current detectors and what will be needed to support the research ambitions of astronomers in the future.     

Robust dual‐sensor data sea‐surface ghost attenuation for quality control purposes

Three‐dimensional receiver ghost attenuation (deghosting) of dual‐sensor towed‐streamer data is straightforward, in principle. In its simplest form, it requires applying a three‐dimensional frequency–wavenumber filter to the vertical component of the particle motion data to correct for the amplitude reduction on the vertical component of non‐normal incidence plane waves before combining with the pressure data. More elaborate techniques use three‐dimensional filters to both components before summation, for example, for ghost wavelet dephasing and mitigation of noise of different strengths on the individual components in optimum deghosting. The problem with all these techniques is, of course, that it is usually impossible to transform the data into the crossline wavenumber domain because of aliasing. Hence, usually, a two‐dimensional version of deghosting is applied to the data in the frequency–inline wavenumber domain. We investigate going down the “dimensionality ladder” one more step to a one‐dimensional weighted summation of the records of the collocated sensors to create an approximate deghosting procedure. We specifically consider amplitude‐balancing weights computed via a standard automatic gain control before summation, reminiscent of a diversity stack of the dual‐sensor recordings. This technique is independent of the actual streamer depth and insensitive to variations in the sea‐surface reflection coefficient. The automatic gain control weights serve two purposes: (i) to approximately correct for the geometric amplitude loss of the Z data and (ii) to mitigate noise strength variations on the two components. Here, Z denotes the vertical component of the velocity of particle motion scaled by the seismic impedance of the near‐sensor water volume. The weights are time‐varying and can also be made frequency‐band dependent, adapting better to frequency variations of the noise. The investigated process is a very robust, almost fully hands‐off, approximate three‐dimensional deghosting step for dual‐sensor data, requiring no spatial filtering and no explicit estimates of noise power. We argue that this technique performs well in terms of ghost attenuation (albeit, not exact ghost removal) and balancing the signal‐to‐noise ratio in the output data. For instances where full three‐dimensional receiver deghosting is the final product, the proposed technique is appropriate for efficient quality control of the data acquired and in aiding the parameterisation of the subsequent deghosting processing.     

Tracer methods used to verify the hypothesis of Cvijić about the underground connection between Prespa and Ohrid Lake

Prespa Lake and Ohrid Lake constitute a hydraulic system shared between Albania, FYR of Macedonia and Greece. Karst rocks separate both lakes. The elevation of Prespa Lake is about 150 m higher than that of Ohrid Lake. Considering these facts, Cvijić formulated in 1906 the hypothesis that Prespa Lake recharges the St. Naum and Tushemisht springs at Ohrid lakeside. Environmental isotopes demonstrated that Prespa Lake recharges about 37–42 and 52–54% of water emerging in St. Naum, and Tushemisht springs, respectively. An artificial tracer experiment carried out in 2002 physically demonstrated the underground connection between both lakes. This experiment confirmed the supposed underground connection and brought important information about the groundwater velocity, transit time, and karst water conduits development.     

The Bayesian approach to an internally consistent thermodynamic database: theory, database, and generation of phase diagrams

An internally consistent thermodynamic dataset has been derived for 148 endmember phases (145 solids and 3 fluids) comprising the elements Li, Na, K, Be, Mg, Ca, Ti, Cr, Mn, Fe, Zn, Al, Si, C, H, and O. This has been achieved by simultaneous treatment of phase property (like standard enthalpy of formation, standard entropy, molar heat capacity, molar volume, thermal expansivity, bulk modulus etc.) and reaction reversal data by the Bayesian method. The theory underlying the approach, and the computational methods involved, are briefly outlined. (For the benefit of readers unfamiliar with inference statistics, the basic concepts of the Bayes method are also presented in such a way that they can be grasped intuitively.) Although not yet addressed, this method can be extended to refine the thermodynamic mixing properties of crystalline solutions. The sources of the input data, culled from the literature, are summarized in the Appendix. The resulting database is succinctly documented in this paper. It includes the enthalpies of formation and entropies, their uncertainties, and the correlation among them. The database allows calculation of P-T, T-X _CO2, P-X _CO2, and T-f _O2 sections, with error propagation into the computed phase diagrams on a routine basis. A user-friendly computer program has been written to generate such phase diagrams. It is public domain software. The software and the thermodynamic database (which includes a complete documentation of the thermodynamic data above and beyond those listed (Table 2, here) may be downloaded from the web site http://homepage.ruhr-uni-bochum.de/niranjan.chatterjee/Index.htm. Examples of computed phase diagrams are given to illustrate the quality of the data and the capabilities of the software. Received: 11 March 1998 / Accepted: 11 June 1998     

Assimilation and diffusion during xenolith-magma interaction: a case study of the Variscan Karkonosze Granite, Bohemian Massif

Magmatic enclaves from the Rudolfov quarry near Liberec (Czech Republic) are interpreted to represent remnants of lamprophyric melt that intruded the Karkonosze granite at a stage at which the granite was not fully solidified. Based on the observation that larger enclaves are generally more circular than the smaller ones, we conclude that bigger blobs of mafic magma became more spherical during flow in the gravity field (sink or float). This flow is also interpreted to be responsible for the incorporation of mineral grains into the enclaves and may have facilitated the assimilation of granitic melt. Linear mixing trends on Harker diagrams for most network-forming and mainly slow-diffusing or fluid-immobile elements indicate such an assimilation process between granite and lamprophyre. In contrast, all fast-diffusing or fluid-mobile elements display scattered element distributions, implying that chemical diffusion also played a role. Pressure and temperature for this late magmatic stage are estimated at around 1 kbar and 500°C. These results suggest that two processes modified the composition of the enclaves in the Karkonosze granite: (1) assimilation (mechanical mixing) of granitic melt during the injection of the lamprophyric melt and the subsequent flow of the forming enclaves in the gravity field (responsible for the linear mixing trends) and (2) diffusion-controlled redistribution of elements between both solidifying rock types at the magmatic stage and/or due to late-stage magmatic fluids (responsible for the scattering and deviation from the linear mixing trends).     

A new method for the determination of the specific kinetic energy (SKE) released to pyroclastic particles at magmatic fragmentation: theory and first experimental results

Brittle magmatic fragmentation plays a crucial role in explosive eruptions. It represents the starting point of hazardous explosive events that can affect large areas surrounding erupting volcanoes. Knowing the initial energy released during this fragmentation process is fundamental for the understanding of the subsequent dynamics of the eruptive gas-particle mixture and consequently for the forecasting of the erupting column’s behavior. The specific kinetic energy (SKE) of the particles quantifies the initial velocity shortly after the fragmentation and is therefore a necessary variable to model the gas-particle conduit flow and eruptive column regime. In this paper, we present a new method for its determination based on fragmentation experiments and identification of the timings of energy release. The results obtained on compositions representative for basaltic and phonolitic melts show a direct dependence on magma material properties: poorly vesiculated basaltic melts from Stromboli show the highest SKE values ranging from 7.3 to 11.8 kJ/kg, while experiments with highly vesiculated samples from Stromboli and Vesuvius result in lower SKE values (3.1 to 3.8 kJ/kg). The described methodology presents a useful tool for quantitative estimation of the kinetic energy release of magmatic fragmentation processes, which can contribute to the improvement of hazard assessment.     

Node-centered finite volume discretizations for the numerical simulation of multiphase flow in heterogeneous porous media

Three node-centered finite volume discretizations for multiphase porous media flow are presented and compared. By combination of these methods two additional discretization methods are generated. The ability of these schemes to describe flows at textural interfaces of different geologic formations is investigated. It was found that models with nonzero-entry pressures for the capillary pressure-saturation relationship in conjunction with the Box discretization may give rise to spurious oscillations for flows around low permeable lenses. Furthermore, the applicability and sensitivity of the discretization methods with regard to the used computational grids is discussed. The schemes are used for the numerical study of two-phase flow in porous media with zones of different material properties. This revised version was published online in July 2006 with corrections to the Cover Date.     

Methane exchange between coal-bearing basins and the atmosphere: the Ruhr Basin and the Lower Rhine Embayment, Germany

A precise knowledge of methane exchange processes is required to fully understand the recent rise of atmospheric methane concentration. Three of these processes take place at the lithosphere/atmosphere boundary: bacterial consumption of methane and emission of bacterial or thermogenic methane. This study was initiated to quantify these processes on a regional scale in the Ruhr Basin and the Lower Rhine Embayment. Since these areas are subject to bituminous coal and lignite mining, natural and anthropogenically-induced methane exchange processes could be studied. The methane emission and consumption rates and their carbon isotope signal were measured at the lithosphere/atmosphere boundary using flux chambers. On most of the soils studied, methane consumption by bacteria was identified. Thermogenic methane was released only at some of the natural faults examined. In active and abandoned bituminous coal mining areas methane emissions were restricted to small areas, where high emission rates were measured. The carbon isotope composition of methane at natural faults and in mining subsidence troughs was typical of thermogenic methane (−45 to −32 ‰ δ¹³C). Methane exchange balancing revealed that natural methane emissions from these two basins represent no source of atmospheric importance. However, methane release by upcast mining shafts dominates the methane exchange processes and is by about two orders of magnitude greater than methane consumption by bacterial oxidation in the soils.