Interplay of phosphogenesis and hydrothermalism in the latest Ediacaran rift of the High Atlas,Morocco

Two dissolution surfaces, capped by phosphatic crusts that contain up to 40 wt.% of P₂O₅, are well exposed in the latest Ediacaran Tabia Member (Adoudou Formation) of the Ouneïn region, High Atlas. Both the Tabia dolostones and capping phosphorites are crosscut by several generations of synsedimentary faults, fissures, and hydrothermal dikes hosting Cu ore bodies, in a scenario of down-faulted blocks. Stratigraphic and facies relationships suggest the preservation of a major pulse of intra-cratonic rifting, responsible for the development of an unstable carbonate platform subjected to sharp uplift and tilting perturbations and dissolution processes, associated with hydrothermal dike-swarm injection. The whole framework is sealed with the Tamjout Bed (Tifnout Member), which marks the end of abundant regional tilting (syn-rift phase) and the beginning of a thermal subsidence-dominated regime (post-rift phase) in the Moroccan margin of West Gondwana.     

Modelling suitable estuarine habitats for Zostera noltii, using Ecological Niche Factor Analysis and Bathymetric LiDAR

Predicting species distribution and habitat suitability is of considerable use in supporting the implementation of environmental legislation, protection and conservation of marine waters and ecosystem-based management. As other seagrasses, Zostera noltii has declined worldwide, mainly due to human pressures, such as eutrophication and habitat loss. In the case of the Basque Country (northern Spain), the species is present only in 3 out of 12 estuaries. From the literature, it is known that at least 6 of these estuaries were formerly vegetated by this seagrass. Consequently, efforts to monitor and restore (potential) habitats have been enhanced. Therefore, we aim: (i) to determine the main environmental variables explaining Zostera noltii distribution, within the Basque estuaries based upon the Oka estuary; (ii) to model habitat suitability for this species, as a wider applicable management-decision tool for seagrass restoration; and (iii) to assess the applicability and predicted accuracy of the model by using internal and external validation methods. For this purpose, Ecological Niche Factor Analysis (ENFA) has been used to model habitat suitability, based upon topographical variables, obtained from bathymetric Light Detection And Ranging (LiDAR); sediment characteristics variables; and hydrodynamic variables. The results obtained from the ecological factors of the ENFA (Marginality: 1.00; Specialization: 2.59) indicate that the species habitat differs considerably from the mean environmental conditions over the study area; likewise, that the species is restrictive in the selection of the range of conditions within which it dwells. The main environmental variables relating to the species distribution, in order of importance, are: mean grain size; redox potential; intertidal height; sediment sorting; slope of intertidal flat; percentage of gravels; and percentage of organic matter content. The model has a high predicted accuracy (Boyce index: 0.92). Model validation using an independent dataset in the Bidasoa estuary has shown the applicability but also the limitations in extrapolating the habitat suitability model to select suitable transplantation areas in other estuaries with similar morphological and biogeographical characteristics. ENFA-technique, applied with an accurate selection of environmental predictors, could be a promising tool for predicting seagrass habitat suitability which could assist on seagrass conservation and restoration programs worldwide.     

Potential collapse of the upper slope and tsunami generation on the Great Barrier Reef margin, north-eastern Australia

Analysis of high-resolution multibeam bathymetry and seismic profiles in the Noggin Passage region, north-eastern Australia, has identified a small area (Noggin block) in the upper-slope offshore Cairns that may potentially collapse and generate a tsunami wave. The Noggin block extends from 340 to 470 m depth covering a roughly circular (2.4 km long and 3.7 km wide) area of about 5.3 km². The well-defined margins of the block correspond to different bounding seabed features. These features include steep headscarps, small landslides and a group of aligned circular pockforms up to 500 m wide and 20 m deep. Slope stability simulations indicate that the Noggin block is stable under normal present-day gravitational conditions on the upper slope. However, block failure may result under external loads, such as those produced by earthquakes. Failure modelling shows that critical peak horizontal accelerations of 0.2–0.4 g could lead to the collapse of the Noggin block. In north-eastern Australia, these acceleration values would involve earthquakes generated at short hypocentral distances and short periods. The collapse of the potential sediment slide mass of about 0.86 km³ (162 m average thickness) may lead to the formation of a landslide-generated tsunami wave. Semi-empirical equations indicate the collapse of this mass would yield a 7–11-m high three-dimensional tsunami wave. These waves could reach an estimated run-up height at the coast of 5–7 m. Our first-order approach highlights the potential consequences for nearby coastal communities, the need for better sediment characterisation in the study area, and the systematic identification of other areas prone to slope failures along the Great Barrier Reef margin.     

Combined U/Pb and (U–Th)/He geochronometry of basalt maars in Western Carpathians: implications for age of intraplate volcanism and origin of zircon metasomatism

The age of intraplate volcanism in northern Pannonian Basin of Carpathians is revisited using a combination of zircon U/Pb, zircon (U–Th)/He and apatite (U–Th)/He dating techniques, complemented by electron microprobe (EMP) characterisation of dated minerals. A total of six maar structures and diatremes in the South-Slovakian Volcanic Field (SSVF) were dated and the obtained new ages yielded the following key findings: Two isolated maars in SE part indirectly dated by geomorphologic constraints to Late Pleistocene are actually of Pliocene (2.8 ± 0.2 Ma) and Late Miocene (5.5 ± 0.6 Ma) ages. In contrast, two maars in NW part of the study area are of Late Pliocene age (4.1 ± 0.4 and 5.2–5.4 Ma), younger than the Late Miocene age (~6.5 Ma) inferred previously from K/Ar data on the proximal basaltic lava flows. These maars therefore belong to the second volcanic phase that was previously identified only in SE part of the SSVF. In the light of the new geochronologic data, it seems likely that the Pliocene phreatomagmatic eruptions may have occurred along extension-related, NW- and NE-trending orthogonal faults. EMP analyses and imaging revealed an extensive syn- and post-growth metasomatic replacement by dissolution-reprecipitation in the majority of zircons. Abundant silicate melt inclusions in porous metasomatised parts of the zircons are diagnostic of magmatic rather than hydrothermal metasomatism. Consistent ages of the metasomatised and non-metasomatised zones do not indicate disturbance of the U–Pb system during the metasomatism. Enrichment in U and Th loss in the metasomatised zircons are diagnostic of an increasing oxygen fugacity triggered by degassing of the volatile residual melt during the final stages of alkali basalt fractionation. Rare zircon-to-baddeleyite transformation was probably connected with lowered silica activity in carbonated basaltic magmas in south-eastern part of the study area.     

Are burnt sediments reliable recorders of geomagnetic field strength?

This study tests if burnt soils and sediments can provide reliable records of geomagnetic field strength at the time of burning by carrying out an experiment to reproduce the prehistoric use of fire on a clayish soil substratum. Rock magnetic experiments showed that in the upper 0–1 cm of the central part of the burnt surface, remanence is a thermoremanent magnetization carried by single-domain magnetite and that samples are thermally stable. Fourteen specimens from that area were subjected to paleointensity experiments with the Coe method (1967). An intensity of 42.9 ± 5.7 μT was estimated below 440°C, whereas at higher temperatures magneto-mineralogical alterations were observed. Corresponding successful microwave intensity determinations from two specimens gave a mean value of 47.6 μT. Both results are in reasonable agreement with the expected field value of 45.2 μT. Burnt soils of archeological fires thus have the potential to record accurately the paleofield strength and may be useful targets for archeointensity investigations. Coincident results obtained from two different paleointensity determination methods support this conclusion.     

Capturing and modelling metre‐scale spatial facies heterogeneity in a Jurassic ramp setting (Central High Atlas,Morocco)

Each simulation algorithm, including Truncated Gaussian Simulation, Sequential Indicator Simulation and Indicator Kriging is characterized by different operating modes, which variably influence the facies proportion, distribution and association of digital outcrop models, as shown in clastic sediments. A detailed study of carbonate heterogeneity is then crucial to understanding these differences and providing rules for carbonate modelling. Through a continuous exposure of Bajocian carbonate strata, a study window (320 m long, 190 m wide and 30 m thick) was investigated and metre‐scale lithofacies heterogeneity was captured and modelled using closely‐spaced sections. Ten lithofacies, deposited in a shallow‐water carbonate‐dominated ramp, were recognized and their dimensions and associations were documented. Field data, including height sections, were georeferenced and input into the model. Four models were built in the present study. Model A used all sections and Truncated Gaussian Simulation during the stochastic simulation. For the three other models, Model B was generated using Truncated Gaussian Simulation as for Model A, Model C was generated using Sequential Indicator Simulation and Model D was generated using Indicator Kriging. These three additional models were built by removing two out of eight sections from data input. The removal of sections allows direct insights on geological uncertainties at inter‐well spacings by comparing modelled and described sections. Other quantitative and qualitative comparisons were carried out between models to understand the advantages/disadvantages of each algorithm. Model A is used as the base case. Indicator Kriging (Model D) simplifies the facies distribution by assigning continuous geological bodies of the most abundant lithofacies to each zone. Sequential Indicator Simulation (Model C) is confident to conserve facies proportion when geological heterogeneity is complex. The use of trend with Truncated Gaussian Simulation is a powerful tool for modelling well‐defined spatial facies relationships. However, in shallow‐water carbonate, facies can coexist and their association can change through time and space. The present study shows that the scale of modelling (depositional environment or lithofacies) involves specific simulation constraints on shallow‐water carbonate modelling methods.     

Palaeoclimate of the Late Jurassic of Portugal: comparison with the Western United States

Investigation of the palaeoclimatic conditions associated with Upper Jurassic strata in Portugal and comparison with published palaeoclimate reconstructions of the Upper Jurassic Morrison Formation in western North America provide important insights into the conditions that allowed two of the richest terrestrial faunas of this period to flourish. Geochemical analyses and observations of palaeosol morphology in the informally named Upper Jurassic Lourinhã formation of western Portugal indicate warm and wet palaeoclimatic conditions with strongly seasonal precipitation patterns. Palaeosol profiles are dominated by carbonate accumulations and abundant shrink‐swell (vertic) features that are both indicative of seasonal variation in moisture availability. The δ¹⁸O_SMOW and δD_SMOW values of phyllosilicates sampled from palaeosol profiles range from +22·4‰ to +22·7‰ and ?53·0‰ to ?37·3‰, respectively. These isotope values correspond to temperatures of formation between 32°C and 39°C ± 3°, with an average of 36°C, which suggest surface temperatures between 27°C and 34°C (average 31°C). On average, these surface temperature estimates are 1°C higher than the highest summer temperatures modelled for Late Jurassic Iberia using general circulation models. Elemental analysis of matrix material from palaeosol B‐horizons provides proxy (chemical index of alteration minus potassium) estimates of mean annual precipitation ranging from 766 to 1394 mm/year, with an average of approximately 1100 mm/year. Palaeoclimatic conditions during deposition of the Lourinhã formation are broadly similar to those inferred for the Morrison Formation, except somewhat wetter. Seasonal variation in moisture availability does not seem to have negatively impacted the ability of these environments to support rich and relatively abundant faunas. The similar climate between these two Late Jurassic terrestrial ecosystems is probably one of the factors which explains the similarity of their vertebrate faunas.     

Approaches for tsunami risk assessment and application to the city of Cádiz, Spain

Tsunamis can represent a significant risk to the population and cause huge economic damage in many costal regions. In order to be able to identify risk hot spots and implement targeted risk reduction measures, decision makers need to have a clear picture of the risk situation in their countries or regions. This work reviews existing approaches for tsunami risk assessment and recommends a five-step process for assessing tsunami risk. As a case study, a qualitative risk assessment for a worst-case tsunami scenario was carried out to understand the tsunami risk to the population in Cádiz. Moreover, a sensitivity analysis of the tsunami hazard input parameters was performed as a strong influence of the variability of the input parameters on the resultant tsunami hazard and risk zonation maps was observed. The study shows that regardless of the assumptions made a non-negligible tsunami risk to Cádiz exists.     

A review of Markov Chain Monte Carlo and information theory tools for inverse problems in subsurface flow

Parameter identification is one of the key elements in the construction of models in geosciences. However, inherent difficulties such as the instability of ill-posed problems or the presence of multiple local optima may impede the execution of this task. Regularization methods and Bayesian formulations, such as the maximum a posteriori estimation approach, have been used to overcome those complications. Nevertheless, in some instances, a more in-depth analysis of the inverse problem is advisable before obtaining estimates of the optimal parameters. The Markov Chain Monte Carlo (MCMC) methods used in Bayesian inference have been applied in the last 10 years in several fields of geosciences such as hydrology, geophysics or reservoir engineering. In the present paper, a compilation of basic tools for inference and a case study illustrating the practical application of them are given. Firstly, an introduction to the Bayesian approach to the inverse problem is provided together with the most common sampling algorithms with MCMC chains. Secondly, a series of estimators for quantities of interest, such as the marginal densities or the normalization constant of the posterior distribution of the parameters, are reviewed. Those reduce the computational cost significantly, using only the time needed to obtain a sample of the posterior probability density function. The use of the information theory principles for the experimental design and for the ill-posedness diagnosis is also introduced. Finally, a case study based on a highly instrumented well test found in the literature is presented. The results obtained are compared with the ones computed by the maximum likelihood estimation approach.     

Tracing groundwater flow using chemical data

Abstract

In the southern part of the Great Hungarian Plain there are two different groundwater flow systems located in the study area. Between these two systems there is a narrow groundwater divide where groundwater flow cannot be detected. In the western part, groundwater of greater hardness moves from the northwest, west and southwest. In the eastern part, groundwater flow is from southeast to northwest. The directions of groundwater flow have been established on the basis of dissolved mono- and di-valent cation concentrations. The major direction of groundwater flow detected by a statistical evaluation of water chemical data agrees with previous geological investigations.