A semi-implicit method for incompressible three-phase flow in porous media

In this paper, we present a semi-implicit method for the incompressible three-phase flow equations in two dimensions. In particular, a high-order discontinuous Galerkin spatial discretization is coupled with a backward Euler discretization in time. We consider a pressure-saturation formulation, decouple the pressure and saturation equations, and solve them sequentially while still keeping each equation implicit in its respective unknown. We present several numerical examples on both homogeneous and heterogeneous media, with varying permeability and porosity. Our results demonstrate the robustness of the scheme. In particular, no slope limiters are required and a relatively large time step may be taken.     

Spatial analysis approaches for the evaluation and protection of groundwater resources in large watersheds of the Canadian Shield

This study focuses on the development of two GIS-based approaches that are used jointly to evaluate the groundwater resources associated with granular aquifers in shield environments. The first approach is a multi-criteria analysis (MCA) using an analytical hierarchic process (AHP) based on geological and hydrogeological data for ranking the probability of finding readily available groundwater resources in a specific territory. The second approach relies on GIS-based geometric calculations that were developed for evaluating the extent and volume of aquifers. The approaches are applied on a 100?×?100 m grid in a 185,000-km² area corresponding to watersheds of the James Bay area in Quebec, Canada. The MCA-AHP approach revealed that the unconfined granular aquifers that present the highest aquifer potential (AP) are sparsely distributed and mostly associated with glaciofluvial formations such as the Harricana and Sakami moraines. The geometric calculations approach allowed for estimating that the total volume of groundwater stored in the unconfined granular aquifers reaches approximately 40 km³. When used jointly, the two approaches reveal that the shallow unconfined aquifers that require increased groundwater protection account for approximately 5% of the territory. In areas of confined granular aquifers, the highest APs are located in river valleys and lowlands. A sensitivity analysis conducted on the MCA-AHP approach revealed that the grid size does not significantly affect the results. Therefore, the approach was expanded northward, to a 490,000-km² territory reaching the Ungava Bay area. The proposed method could be adapted and applied in other shield areas.

     

Tectono-sedimentary evolution of Jurassic–Cretaceous diapiric structures: Miravete anticline,Maestrat Basin,Spain

Integration of extensive fieldwork, remote sensing mapping and 3D models from high-quality drone photographs relates tectonics and sedimentation to define the Jurassic–early Albian diapiric evolution of the N–S Miravete anticline, the NW-SE Castel de Cabra anticline and the NW-SE Cañada Vellida ridge in the Maestrat Basin (Iberian Ranges, Spain). The pre shortening diapiric structures are defined by well-exposed and unambiguous halokinetic geometries such as hooks and flaps, salt walls and collapse normal faults. These were developed on Triassic salt-bearing deposits, previously misinterpreted because they were hidden and overprinted by the Alpine shortening. The Miravete anticline grew during the Jurassic and Early Cretaceous and was rejuvenated during Cenozoic shortening. Its evolution is separated into four halokinetic stages, including the latest Alpine compression. Regionally, the well-exposed Castel de Cabra salt anticline and Cañada Vellida salt wall confirm the widespread Jurassic and Early Cretaceous diapiric evolution of the Maestrat Basin. The NE flank of the Cañada Vellida salt wall is characterized by hook patterns and by a 500-m-long thin Upper Jurassic carbonates defining an upturned flap, inferred as the roof of the salt wall before NE-directed salt extrusion. A regional E-W cross section through the Ababuj, Miravete and Cañada-Benatanduz anticlines shows typical geometries of salt-related rift basins, partly decoupled from basement faults. These structures could form a broader diapiric region still to be investigated. In this section, the Camarillas and Fortanete minibasins displayed well-developed bowl geometries at the onset of shortening. The most active period of diapiric growth in the Maestrat Basin occurred during the Early Cretaceous, which is also recorded in the Eastern Betics, Asturias and Basque-Cantabrian basins. This period coincides with the peak of eastward drift of the Iberian microplate, with speeds of 20 mm/year. The transtensional regime is interpreted to have played a role in diapiric development.     

Estimating average shock pressures recorded by impactite samples based on universal stage investigations of planar deformation features in quartz—Sources of error and recommendations

Planar deformation features (PDFs) in quartz are the most widely used indicator of shock metamorphism in terrestrial rocks. They can also be used for estimating average shock pressures that quartz-bearing rocks have been subjected to. Here we report on a number of observations and problems that we have encountered when performing universal stage measurements and crystallographically indexing of PDF orientations in quartz. These include a comparison between manual and automated methods of indexing PDFs, an evaluation of the new stereographic projection template, and observations regarding the PDF statistics related to the c-axis position and rhombohedral plane symmetry. We further discuss the implications that our findings have for shock barometry studies. Our study shows that the currently used stereographic projection template for indexing PDFs in quartz might induce an overestimation of rhombohedral planes with low Miller–Bravais indices. We suggest, based on a comparison of different shock barometry methods, that a unified method of assigning shock pressures to samples based on PDFs in quartz is necessary to allow comparison of data sets. This method needs to take into account not only the average number of PDF sets/grain but also the number of high Miller–Bravais index planes, both of which are important factors according to our study. Finally, we present a suggestion for such a method (which is valid for nonporous quartz-bearing rock types), which consists of assigning quartz grains into types (A–E) based on the PDF orientation pattern, and then calculation of a mean shock pressure for each sample.     

A recent history of science cases for optical interferometry

Optical long-baseline interferometry is a unique and powerful technique for astronomical research. Since the 1980’s (with I2T, GI2T, Mark I to III, SUSI, ...), optical interferometers have produced an increasing number of scientific papers covering various fields of astrophysics. As current interferometric facilities are reaching their maturity, we take the opportunity in this paper to summarize the conclusions of a few key meetings, workshops, and conferences dedicated to interferometry. We present the most persistent recommendations related to science cases and discuss some key technological developments required to address them. In the era of extremely large telescopes, optical long-baseline interferometers will remain crucial to probe the smallest spatial scales and make breakthrough discoveries.     

Why chromatic imaging matters

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u ? v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction.     

The variability of ruthenium in chromite from chassignite and olivine‐phyric shergottite meteorites: New insights into the behavior of PGE and sulfur in Martian magmatic systems

The Martian meteorites comprise mantle‐derived mafic to ultramafic rocks that formed in shallow intrusions and/or lava flows. This study reports the first in situ platinum‐group element data on chromite and ulvöspinel from a series of dunitic chassignites and olivine‐phyric shergottites, determined using laser‐ablation ICP‐MS. As recent studies have shown that Ru has strongly contrasting affinities for coexisting sulfide and spinel phases, the precise in situ analysis of this element in spinel can provide important insights into the sulfide saturation history of Martian mantle‐derived melts. The new data reveal distinctive differences between the two meteorite groups. Chromite from the chassignites Northwest Africa 2737 (NWA 2737) and Chassigny contained detectable concentrations of Ru (up to ~160 ppb Ru) in solid solution, whereas chromite and ulvöspinel from the olivine‐phyric shergottites Yamato‐980459 (Y‐980459), Tissint, and Dhofar 019 displayed Ru concentrations consistently below detection limit (<42 ppb). The relatively elevated Ru signatures of chromite from the chassignites suggest a Ru‐rich (~1–4 ppb) parental melt for this meteorite group, which presumably did not experience segregation of immiscible sulfide liquids over the interval of mantle melting, melt ascent, and chromite crystallization. The relatively Ru‐depleted signature of chromite and ulvöspinel from the olivine‐phyric shergottites may be the consequence of relatively lower Ru contents (<1 ppb) in the parental melts, and/or the presence of sulfides during the crystallization of the spinel phases. The results of this study illustrate the significance of platinum‐group element in situ analysis on spinel phases to decipher the sulfide saturation history of magmatic systems.     

On the occurrence and origin of anthropogenic radionuclides found in a fragment of the Chelyabinsk (LL5) meteorite

A piece of the 2013 Chelyabinsk meteorite was investigated for its content of anthropogenic radionuclides. In addition to traces of cesium‐137 that had been previously reported for this particular fragment, we found an unusually high amount of strontium‐90, which indicates that the source of this contamination was the Kyshtym accident (1957). A high Sr‐90/Cs‐137 activity ratio is characteristic for Kyshtym‐derived contaminations. Based on the cesium‐137 content in the soil from the finding site, it is estimated that the fragment was contaminated with soil particles in the milligram range upon impact. Investigation of the soil revealed very unusual ferromagnetic characteristics and an iron‐rich chemical composition. Mössbauer spectroscopy indicated the presence of steel components in this soil, suggesting that the investigated meteorite fragment was found in an industrial dumping site rather than natural soil.