Feasibility of geoelectrical monitoring and multiphase modeling for process understanding of gaseous CO2 injection into a shallow aquifer

Potential pathways in the subsurface may allow upwardly migrating gaseous CO₂ from deep geological storage formations to be released into near surface aquifers. Consequently, the availability of adequate methods for monitoring potential CO₂ releases in both deep geological formations and the shallow subsurface is a prerequisite for the deployment of Carbon Capture and Storage technology. Geoelectrical surveys are carried out for monitoring a small-scale and temporally limited CO₂ injection experiment in a pristine shallow aquifer system. Additionally, the feasibility of multiphase modeling was tested in order to describe both complex non-linear multiphase flow processes and the electrical behavior of partially saturated heterogeneous porous media. The suitability of geoelectrical methods for monitoring injected CO₂ and geochemically altered groundwater was proven. At the test site, geoelectrical measurements reveal significant variations in electrical conductivity in the order of 15?C30?%. However, site-specific conditions (e.g., geological settings, groundwater composition) significantly influence variations in subsurface electrical conductivity and consequently, the feasibility of geoelectrical monitoring. The monitoring results provided initial information concerning gaseous CO₂ migration and accumulation processes. Geoelectrical monitoring, in combination with multiphase modeling, was identified as a useful tool for understanding gas phase migration and mass transfer processes that occur due to CO₂ intrusions in shallow aquifer systems.     

Simulation of DNAPL infiltration and spreading behaviour in the saturated zone at varying flow velocities and alternating subsurface geometries

The influence of varying groundwater flow velocities on DNAPL infiltration and spreading behaviour was investigated by multiphase modelling using TMVOC and PetraSim. The multiphase models were calibrated by results of previously conducted laboratory experiments for the complete spatio-temporal range of the experiments. The small scale 2D scenario modelling was applied to qualify and quantify changes in position, architecture, geometry and dissolution of a TCE body in a fully saturated homogeneous sandy medium. The applied flow velocities ranging from 0.05 up to 40.00 m/day exhibited that the DNAPL TCE is affected even at the lowest flow velocity in its position, its size and its architecture. Additionally, several impermeable lenses with simple geometry were assumed in the model, to investigate the influence of stratified subsoil. In the experimental set-ups, the DNAPL body reacts more sensitive to the applied groundwater flow velocities than to the geometrical set-up of the scenarios. A possible consequence can be the transportation and displacement of a DNAPL pool due to natural or anthropogenic induced high groundwater flow velocities, as by Pump and Treat facilities, complicating site investigation process and planning of remediation activities.     

Mendigite,Mn2Mn2MnCa(Si3O9)2, a new mineral species of the bustamite group from the Eifel volcanic region,Germany

A new mineral, mendigite (IMA no. 2014-007), isostructural with bustamite, has been found in the In den Dellen pumice quarry near Mendig, Laacher Lake area, Eifel Mountains, Rhineland-Palatinate (Rheinland-Pfalz), Germany. Associated minerals are sanidine, nosean, rhodonite, tephroite, magnetite, and a pyrochlore-group mineral. Mendigite occurs as clusters of long-prismatic crystals (up to 0.1 × 0.2 × 2.5 mm in size) in cavities within sanidinite. The color is dark brown with a brown streak. Perfect cleavage is parallel to (001). D _calc = 3.56 g/cm³. The IR spectrum shows the absence of H₂O and OH groups. Mendigite is biaxial (–), α = 1.722 (calc), β = 1.782(5), γ = 1.796(5), 2V _meas = 50(10)°. The chemical composition (electron microprobe, mean of 4 point analyses, the Mn²⁺/Mn³⁺ ratio determined from structural data and charge-balance constraints) is as follows (wt %): 0.36 MgO, 10.78 CaO, 37.47 MnO, 2.91 Mn₂O₃, 4.42 Fe₂O₃, 1.08 Al₂O₃, 43.80 SiO₂, total 100.82. The empirical formula is Mn_2.00(Mn_1.33Ca_0.67) (Mn_0.50 ²⁺ Mn_0.28 ³⁺ Fe_0.15 ³⁺ Mg_0.07)(Ca_0.80 (Mn_0.20 ²⁺)(Si_5.57 Fe_0.27 ³⁺ Al_0.16O₁₈). The idealized formula is Mn₂Mn₂MnCa(Si₃O₉)₂. The crystal structure has been refined for a single crystal. Mendigite is triclinic, space group \(P\bar 1\); the unit-cell parameters are a = 7.0993(4), b = 7.6370(5), c = 7.7037(4) Å, α = 79.58(1)°, β = 62.62(1)°, γ = 76.47(1)°; V = 359.29(4) ų, Z = 1. The strongest reflections on the X-ray powder diffraction pattern [d, Å (I, %) (hkl)] are: 3.72 (32) (020), 3.40 (20) (002, 021), 3.199 (25) (012), 3.000 (26), (\(01\bar 2\), \(1\bar 20\)), 2.885 (100) (221, \(2\bar 11\), \(1\bar 21\)), 2.691 (21) (222, \(2\bar 10\)), 2.397 (21) (\(02\bar 2\), \(21\bar 1\), 203, 031), 1.774 (37) (412, \(3\bar 21\)). The type specimen is deposited in the Fersman Mineralogical Museum, Russian Academy of Sciences, Moscow, registration number 4420/1.     

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.     

Magmatic-to-hydrothermal crystallization in the W–Sn mineralized Mole Granite (NSW, Australia): Part I: Crystallization of zircon and REE-phosphates over three million years—a geochemical and U–Pb geochronological study

Zircon, monazite and xenotime crystallized over a temperature interval of several hundred degrees at the magmatic to hydrothermal transition of the Sn and W mineralized Mole Granite. Magmatic zircon and monazite, thought to have crystallized from hydrous silicate melt, were dated by conventional U–Pb techniques at an age of 247.6 ± 0.4 and 247.7 ± 0.5 Ma, respectively. Xenotime occurring in hydrothermal quartz is found to be significantly younger at 246.2 ± 0.5 Ma and is interpreted to represent hydrothermal growth. From associated fluid inclusions it is concluded that it precipitated from a hydrothermal brine ≤ 600 °C, which is below the accepted closure temperature for U–Pb in this mineral. These data are compatible with a two-stage crystallization process: precipitation of zircon and monazite as magmatic liquidus phases in deep crustal magma followed by complete crystallization and intimately associated Sn–W mineralization after intrusion of the shallow, sill-like body of the Mole Granite. Later hydrothermal formation of monazite in a biotite–fluorite–topaz reaction rim around a mineralized vein was dated at 244.4 ± 1.4 Ma, which distinctly postdates the Mole Granite and is possibly related to a younger hidden intrusion and its hydrothermal fluid system.

Obtaining precise age data for magmatic and hydrothermal minerals of the Mole Granite is hampered by uncertainties introduced by different corrections required for multiple highly radiogenic minerals crystallising from evolved hydrous granites, including ²³⁰Th disequilibrium due to Th/U fractionation during monazite and possibly xenotime crystallization, variable Th/U ratios of the fluids from which xenotime was precipitating, elevated contents of common lead, and post-crystallization lead loss in zircon, enhanced by the fluid-saturated environment. The data imply that monazite can also survive as a liquidus phase in protracted magmatic systems over periods of 10⁶ years. The outlined model is in agreement with prominent chemical core-rim variation of the zircon.     

Climatic variability in the Late Copper Age: stable isotope fluctuation of prehistoric Unio pictorum (Unionidae) shells from Lake Balaton (Hungary)

Geochemical records of bivalve shells have been increasingly studied in the last decade to obtain information on climate conditions. In this paper we present stable isotope compositions of living and prehistoric shells of freshwater mussels (Unionidae) and their relationships with climate conditions in a shallow lake environment of Lake Balaton, West-Central Hungary. Physical conditions and stable oxygen isotope compositions of lake water samples were monitored where living bivalves were collected. Comparisons between seasonal variations in ambient temperature, water composition and within-shell isotopic variations indicate that the shells of Unio pictorum do reflect local changes at high resolution and thus can be used to study past conditions. Additionally, shells covering the last two decades were gathered at several locations along the lake in order to determine spatial and temporal variations in the shells’ isotopic compositions as a function of weather conditions. As an application, prehistoric shells collected in archaeological excavations were analysed in order to study past environmental variations. Climate variations during the Late Copper Age (5460–4870 cal. yr BP) have been assumed on the basis of geomorphological and archaeozoological observations at the site Balatonkeresztúr-Réti-dűlő (south of Lake Balaton), that suggested increasing humidity as a cause of changes in settlement location and domestic livestock husbandry. Stable carbon and oxygen isotope compositions of prehistoric bivalve shells were analyzed from excavations representing five archaeological subphases (Boleráz subphase, 5460–5310 cal. yr BP; two transitional subphases around 5310 cal. yr BP; Early Classic subphase, 5310–5060 cal. yr BP; Late Classic subphase, 5040–4870 cal. yr BP). The analyses revealed significant negative C and O isotope shifts in the transitional subphases relative to the earlier and later subphases. The isotopic variations indicate that the local climate became relatively wet and possibly cold around 5310 cal. yr BP, then it returned to drier (and likely warmer) conditions during the Classic subphases. This interpretation is in agreement with previous studies on climate changes related to the “5.3 ky event” in the European continental area and the North Atlantic Region, indicating an Atlantic influence in the Carpathian Basin.     

Uncertainty in GPS Networks due to Remaining Systematic Errors: The Interval Approach

A realistic assessment of the total uncertainty budget of Global Positioning System (GPS) observations and its adequate mathematical treatment is a basic requirement for all analysis and interpretation of GPS-derived point positions, in particular GPS heights, and their respective changes. This implies not only the random variability but also the remaining systematic errors. At present in geodesy, the main focus is on stochastic approaches in which errors are modeled by means of random variables. Here, an alternative approach based on interval mathematics is presented. It allows us to model and to quantify the impact of remaining systematic errors in GPS carrier-phase observations on the final results using deterministic error bands. In this paper, emphasis is given to the derivation of the observation intervals based on influence parameters and to the study of the complex linear transfer of this type of uncertainty to estimated point positions yielding zonotopes. From the presented simulation studies of GPS baselines, it turns out that the uncertainty due to remaining systematic effects dominates the total uncertainty budget for baselines longer than 200 km.     

Storm-induced turbidity currents on a sediment-starved shelf: Insight from direct monitoring and repeat seabed mapping of upslope migrating bedforms

The monitoring of turbidity currents enables accurate internal structure and timing of these flows to be understood. Without monitoring, triggers of turbidity currents often remain hypothetical and are inferred from sedimentary structures of deposits and their age. In this study, the bottom currents within 20 m of the seabed in one of the Pointe-des-Monts (Gulf of St. Lawrence, eastern Canada) submarine canyons were monitored for two consecutive years using Acoustic Doppler Current Profilers. In addition, multibeam bathymetric surveys were carried out during deployment of the Acoustic Doppler Current Profilers and recovery operations. These new surveys, along with previous multibeam surveys carried out over the last decade, revealed that crescentic bedforms have migrated upslope by about 20 to 40 m since 2007, despite the limited supply of sediment on the shelf or river inflow in the region. During the winter of 2017, two turbidity currents with velocities reaching 0·5 m sec⁻¹ and 2·0 m sec⁻¹, respectively, were recorded and were responsible for the rapid (<1 min) upstream migration of crescentic bedforms measured between the autumn surveys of 2016 and 2017. The 200 kg (in water) mooring was also displaced 10 m down-canyon, up the stoss side of a bedform, suggesting that a dense basal layer could be driving the flow during the first minute of the event. Two other weaker turbidity currents with speeds <0·5 m sec⁻¹ occurred, but did not lead to any significant change on the seabed. These four turbidity currents coincided with strong and sustained wind speed >60 km h⁻¹ and higher than normal wave heights. Repeat seabed mapping suggests that the turbidity currents cannot be attributed to a canyon-wall slope failure. Rather, sustained windstorms triggered turbidity currents either by remobilizing limited volumes of sediment on the shelf or by resuspending sediment in the canyon head. Turbidity currents can thus be triggered when the sediment volume available is limited, likely by eroding and incorporating canyon thalweg sediment in the flow, thereby igniting the flow. This process appears to be particularly important for the generation of turbidity currents capable of eroding the lee side of upslope migrating bedforms in sediment-starved environments and might have wider implications for the activity of submarine canyons worldwide. In addition, this study suggests that a large external trigger (in this case storms) is required to initiate turbidity currents in sediment-starved environments, which contrasts with supply-dominated environments where turbidity currents are sometimes recorded without a clear triggering mechanism.     

Determination of vertical temperature profiles for the atmospheric boundary layer by ground-based microwave radiometry

A useful method for remote sensing of vertical temperature profiles in the atmospheric boundary layer is described. From angular measurements of brightness temperature at 58 GHz, profiles have been inferred up to an altitude of 700 m. Calculations were done with an iterative inversion procedure (Smith et al., 1972) using Twomey-type smoothing (Twomey, 1963). It is shown how an initial-guess profile can be directly derived from the radiation measurements using, a nomogram.     

The Ordovician chondrite from Brunflo, central Sweden: III. Geochemistry of terrestrial alteration

The fossil H chondrite Brunflo, found in a slab of Ordovician limestone from central Sweden, is pervasively altered to an assemblage dominated by calcite and barite. The meteorite is surrounded by a 15–20 cm wide zone of lighter colors than the unaffected limestone due to dissolution of hematite. Here we present detailed geochemical analyses of two meteorite samples, 14 limestone samples at distances from 0 to 29 cm along two profiles from the meteorite, and a reference sample of Brunflo limestone. Element concentrations in Brunflo and surrounding bleached limestone have been strongly disturbed during two stages of alteration (early oxygenated and deep burial). In the meteorite, the Ni/Co ratio has changed from an initial value of 20 to 0.8 and redox sensitive elements like V, As, Mo, Re and U are strongly enriched. The sulfur isotope composition of barite from Brunflo (δ³⁴S=+35‰) indicates initial loss of meteoritic sulfide, followed by later accumulation of sea water sulfate as barite. During deep burial under more reducing conditions, reduction processes supported by an externally derived reductant possibly derived from alum shale underlying the limestone, were largely responsible for the observed redox phenomena. In spite of massive redistribution of many elements, concentrations of Pt, Ir and Au remain at chondritic levels. The geochemistry and mineralogy of alteration determined for Brunflo are similar to those in “reduction spots” in red beds, where accumulation of a similar suite of elements (except Mo, Re) occurred as a result of isolated reduction activity.