Indirect unstructured hex-dominant mesh generation using tetrahedra recombination

Corner-point gridding is widely used in reservoir and basin modeling but generally yields approximations in the representation of geological interfaces. This paper introduces an indirect method to generate a hex-dominant mesh conformal to 3D geological surfaces and well paths suitable for finite-element and control-volume finite-element simulations. By indirect, we mean that the method first generates an unstructured tetrahedral mesh whose tetrahedra are then merged into primitives (hexahedra, prisms, and pyramids). More specifically, we focus on determining the optimal set of primitives that can be recombined from a given tetrahedral mesh. First, we detect in the tetrahedral mesh all the feasible volumetric primitives using a pattern-matching algorithm (Meshkat and Talmor Int. J. Numer. Meth. Eng. 49(1-2), 17–30 2000) that we re-visit and extend with configurations that account for degenerated tetrahedra (slivers). Then, we observe that selecting the optimal set of primitives among the feasible ones can be formalized as a maximum weighted independent set problem (Bomze et al. 1999), known to be \(\mathcal {N}\mathcal {P}\)-Complete. We propose several heuristic optimizations to find a reasonable set of primitives in a practical time. All the tetrahedra of each selected primitive are then merged to build the final unstructured hex-dominant mesh. This method is demonstrated on 3D geological models including a faulted and folded model and a discrete fracture network.     

Geology and petroleum potential of the Fairway Basin in the Tasman Sea

The Fairway Basin is a large, generally north – south-trending, sediment-filled structure in water 1500 – 3000 m deep, on the eastern slope of the Lord Howe Rise in the Tasman Sea, and is partly within Australian jurisdiction. It was poorly known until a few years ago, when seismic profiling and piston coring cruises were carried out. The basin, about 1100 km long and 120 – 200 km wide, can be divided into three segments—north, central and south—that trend northwest, north and north-northwest, respectively. All three segments probably formed by thinning of continental crust during breakup of Lord Howe Rise and surrounding aseismic continental ridges in the Late Cretaceous and Paleocene. Normal faulting, large inputs of terrigenous sediment and subsidence to bathyal marine depths occurred during that time. A period of compression, perhaps related to overthrusting on New Caledonia, occurred in the Eocene, leading to uplift (and in parts, erosion) of northern Lord Howe Rise, and reversal of faulting in the basin. By the Oligocene, the area was again in bathyal depths, and pelagic ooze and some turbidites accumulated. The basinal sequence is generally 2000 – 4000 m thick, with 1200 – 3200 m of Cretaceous to Eocene sediment concentrated in depocentres, capped by 500 – 800 m of Oligocene and younger sediment. In the depocentres, numerous sedimentary diapirs pierce sedimentary sequences. The sedimentary diapirs appear to be fed by Cretaceous muds deposited during rifting. Often, these diapirs are overlain by faults extending to the seafloor, and hummocky bathymetry is possibly caused by fluid escape. The overall geology suggests that the Fairway Basin may be a large frontier hydrocarbon province. Seismic profiles display a bottom-simulating reflector above many depocentres, 500 – 700 m below the seafloor. The bottom-simulating reflector has positive polarity, which could result from a diagenetic phase transformation, a thin gas hydrate layer with a sharp top, or from the sharp base of a gas layer (probably beneath gas hydrates). Standard piston cores taken above diapirs and apparent fluid-escape features have recovered little gas. Other than drilling, the next steps in assessing petroleum potential are to clearly document fluid-escape structures, and to sample any fluids emitted for hydrocarbons.     

Alpine tectonics in the Calabrian–Peloritan belt (southern Italy): New 40Ar/39Ar data in the Aspromonte Massif area

This study provides new ⁴⁰Ar/³⁹Ar geochronological constraints on the age of the Alpine tectonics in the Aspromonte Massif (southern part of the Calabrian–Peloritan belt). This massif exposes the upper units of the Calabride Complex which originated from the European continental margin. The Calabride Complex was incorporated in the Alpine orogenic wedge and then integrated into the Apennines and Maghrebides fold-and-thrust belts. Throughout the Calabride Complex there is evidence for a two stage tectonic history, which remains however rather poorly dated: Alpine nappe stacking is followed by extensional reworking along the former thrust contacts or along new detachment surfaces. Our new ages suggest that exhumation of the uppermost units, which accompanied nappe stacking, probably started at 45 Ma and that the deepest units were almost completely exhumed at 33 Ma. This kinematics probably corresponds to syn-orogenic extension while the end of exhumation is clearly related to the extensional tectonics dated at 28.6 Ma along detachment structures.Our geochronological data reveal a very short lag time between accretional and extensional processes in this part of the Mediterranean Alpine orogenic belt. The direction of extension, when the units are restored to their initial position (i.e. before the opening of the Western Mediterranean basins and the bending of the arc) is NNE–SSW. Such a direction does not fit with the eastward slab-retreat model generally put forward to explain extension in the Western Mediterranean. In contrast, we provide evidence for roughly N–S middle Oligocene extension in the accretionary prism, not previously described in this part of the Mediterranean domain.     

Tableau de répartition stratigraphique des Grands Foraminifères caractéristiques du Crétacé moyen de la région Méditerranéenne

Le Groupe de Travail Européen des Grands Foraminifères présente un tableau de répartition stratigraphique de 42 espèces bien définies du Crétacé moyen de la région méditerranéenne. La répartition stratigraphique proposée pour chaque espèce est fondée soit sur des propres observations soit sur des données de la littérature supprimer et est contrôlée par la présence d'Ammonites ou de Foraminifères planctoniques.A summarizing account is presented of the deliberations of the research group for Large Foraminifers of the IGCP Project “Mid-Cretaceous Events”. Large Foraminifers are of incontestable value for dating carbonate platform sequences owing to the absence of many other diagnostic groups of organísms. A table of stratigraphical distributions for 42 species is presented.     

External controls on the distribution,fabrics and mineralization of modern microbial mats in a coastal hypersaline lagoon,Cayo Coco (Cuba)

Active, carbonate‐mineralizing microbial mats flourish in a tropical, highly evaporative, marine‐fed lagoonal network to the south of Cayo Coco Island (Cuba). Hypersaline conditions support the development of a complex sedimentary microbial ecosystem with diverse morphologies, a variable intensity of mineralization and a potential for preservation. In this study, the role of intrinsic (i.e. microbial) and extrinsic (i.e. physicochemical) controls on microbial mat development, mineralization and preservation was investigated. The network consists of lagoons, forming in the interdune depressions of a Pleistocene aeolian substratum; they developed due to a progressive increase in sea‐level since the Holocene. The hydrological budget in the Cayo Coco lagoonal network changes from west to east, increasing the salinity. This change progressively excludes grazers and increases the saturation index of carbonate minerals, favouring the development and mineralization of microbial mats in the easternmost lagoons. Detailed mapping of the easternmost lagoon shows four zones with different flooding regimes. The microbial activity in the mats was recorded using light–dark shifts in conjunction with microelectrode O₂ and HS^? profiles. High rates of O₂ production and consumption, in addition to substantial amounts of exopolymeric substances, are indicative of a potentially strong intrinsic control on mineralization. Seasonal, climate‐driven water fluctuations are key for mat development, mineralization, morphology and distribution. Microbial mats show no mineralization in the permanently submersed zone, and moderate mineralization in zones with alternating immersion and exposure. It is suggested that mineralization is also driven by water‐level fluctuations and evaporation. Mineralized mats are laminated and consist of alternating trapping and binding of grains and microbially induced magnesium calcite and dolomite precipitation. The macrofabrics of the mats evolve from early colonizing Flat mats to complex Cerebroid or Terrace structures. The macrofabrics are influenced by the hydrodynamic regime: wind‐driven waves inducing relief terraces in windward areas and flat morphologies on the leeward side of the lagoon. Other external drivers include: (i) storm events that either promote (for example, by bioclasts covering) or prevent (for example, by causing erosion) microbial mat preservation; and (ii) subsurface degassing, through mangrove roots and desiccation cracks covered by Flat mats (i.e. forming Hemispheroids and Cerebroidal structures). These findings provide in‐depth insights into understanding fossil microbialite morphologies that formed in lagoonal settings.     

The lacustrine microbial carbonate factory of the successive Lake Bonneville and Great Salt Lake,Utah, USA

The Bonneville Basin is a continental lacustrine system accommodating extensive microbial carbonate deposits corresponding to two distinct phases: the deep Lake Bonneville (30 000 to 11 500 ¹⁴C bp ) and the shallow Great Salt Lake (since 11 500 ¹⁴C bp ). A characterization of these microbial deposits and their associated sediments provides insights into their spatio‐temporal distribution patterns. The Bonneville phase preferentially displays vertical distribution of the microbial deposits resulting from high‐amplitude lake level variations. Due to the basin physiography, the microbial deposits were restricted to a narrow shoreline belt following Bonneville lake level variations. Carbonate production was more efficient during intervals of relative lake level stability as recorded by the formation of successive terraces. In contrast, the Great Salt Lake microbial deposits showed a great lateral distribution, linked to the modern flat bottom configuration. A low vertical distribution of the microbial deposits was the result of the shallow water depth combined with a low amplitude of lake level fluctuations. These younger microbial deposits display a higher diversity of fabrics and sizes. They are distributed along an extensive ‘shore to lake’ transect on a flat platform in relation to local and progressive accommodation space changes. Microbial deposits are temporally discontinuous throughout the lake history showing longer hiatuses during the Bonneville phase. The main parameters controlling the rate of carbonate production are related to the interaction between physical (kinetics of the mineral precipitation, lake water temperature and runoff), chemical (Ca²⁺, Mg²⁺ and HCO₃^? concentrations, Mg/Ca ratio, dilution and depletion) and/or biological (trophic) factors. The contrast in evolution of Lake Bonneville and Great Salt Lake microbial deposits during their lacustrine history leads to discussions on major chemical and climatic changes during this interval as well as the role of physiography. Furthermore, it provides novel insights into the composition, structure and formation of microbialite‐rich carbonate deposits under freshwater and hypersaline conditions.     

How do granitoid magmas mix with each other? Insights from textures,trace element and Sr–Nd isotopic composition of apatite and titanite from the Matok pluton (South Africa)

In plutonic systems, magma mixing is often modelled by mass balance based on whole-rock geochemistry. However, magma mixing is a chaotic process and chemical equilibration is controlled by non-linear diffusive–advective processes unresolved by the study of bulk samples. Here we present textural observations, LA-(MC-)ICP-MS trace element and Sr–Nd isotopic data of accessory apatites and titanites from a hybrid granodiorite of the Neoarchean Matok pluton (South Africa), collected in a zone of conspicuous mixing between mafic and felsic magmas. Apatite grains mostly show a pronounced zoning in CL images, corresponding to abrupt changes in REE and HFSE concentrations recording their transfer through compositionally different melt domains during mixing. These grains crystallized early, at temperatures of 950–1000 °C. Titanite grains crystallized at temperatures of 820–900 °C (Zr-in-sphene thermometry). They show limited intra-grain chemical variations but huge inter-grain compositional scatter in REE and HFSE, pinpointing crystallization within a crystal mush, from isolated melt pockets having different composition from one another owing to incomplete chemical homogenization and variable Rayleigh fractionation. These chemical–textural characteristics, in combination with partitioning models and Polytopic Vector Analysis, point to “self-mixing” between co-genetic dioritic and granodioritic/granitic magmas. Both resulted from differentiation of mantle-derived mafic melts, showing that mixing does not necessarily involve magmas from contrasted (crust vs. mantle) sources. Systematic variations in εNd _t (?4.5 to ?2.5) and ⁸⁷Sr/⁸⁶Sr_(i) (0.703–0.707) of titanite and apatite grains/domains crystallized from the two magmas point to an isotopically inhomogeneous mantle source, which is not resolved by bulk-rock isotopic data. Interaction between the two magmas must have occurred at relatively high temperatures (ca. 900°C) so that their viscosity contrast remained low, allowing efficient mechanical mixing. Despite this, chemical homogenization was incomplete, as recorded by diffusive fractionation between REE–HFSE and Sr. Modelling thereof reveals that chemical exchange between the liquid phases of the two mixed magmas did not last more than a few tens to hundreds of years. The chemical equilibration between mixed magmas thus strongly depends on the considered elements and observational length scales.     

An Eocene magmatic belt across central Tibet: mantle subduction triggered by the Indian collision?

Mountain domes rising to ≈ 6600 m along Tibet's Tanggula range herald the Eocene intrusion of calcalkaline granites into terranes accreted much earlier. Together with coeval, cogenetic volcanics, such intrusives, which have similar crystallization and cooling ages, may be part of an ESE-trending belt cutting central Tibet in half. This magmatic belt may have marked a former northern boundary of the plateau, testifying to S-directed subduction of the Asian mantle. Such subduction would have developed soon after India's collision but long before the rise of the plateau's present rims, along one of central Tibet's Mesozoic sutures.     

Evaluating flash-flood warnings at ungauged locations using post-event surveys: a case study with the AIGA warning system

Abstract

This article presents a comparison between real-time discharges calculated by a flash-flood warning system and post-event flood peak estimates. The studied event occurred on 15 and 16 June 2010 at the Argens catchment located in the south of France. Real-time flood warnings were provided by the AIGA (Adaptation d’Information Géographique pour l’Alerte en Crue) warning system, which is based on a simple distributed hydrological model run at a 1-km² resolution using radar rainfall information. The timing of the warnings (updated every 15 min) was compared to the observed flood impacts. Furthermore, “consolidated” flood peaks estimated by an intensive post-event survey were used to evaluate the AIGA-estimated peak discharges. The results indicated that the AIGA warnings clearly identified the most affected areas. However, the effective lead-time of the event detection was short, especially for fast-response catchments, because the current method does not take into account any rainfall forecast. The flood peak analysis showed a relatively good correspondence between AIGA- and field-estimated peak values, although some differences were due to the rainfall underestimation by the radar and rainfall–runoff model limitations.

Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Javelle, P., Demargne, J., Defrance, D., Pansu, J. and Arnaud, P., 2014. Evaluating flash-flood warnings at ungauged locations using post-event surveys: a case study with the AIGA warning system. Hydrological Sciences Journal, 59 (7), 1390–1402. http://dx.doi.org/10.1080/02626667.2014.923970