Geosciences, Volume 11, Issue 4 (April 2021) – 34 articles

Leak-off pressure (LOP) is a key parameter to determine the allowable weight of drilling mud in a well and the in situ horizontal stress. The LOP test is run in situ and is frequently used by the petroleum industry. If the well pressure exceeds the LOP, wellbore instability may occur, with hydraulic fracturing and large mud losses in the formation. A reliable prediction of LOP is required to ensure safe and economical drilling operations. The prediction of LOP is challenging because it is affected by the usually complex earlier geological loading history, and the values of LOP and their measurements can vary significantly geospatially. This paper investigates the ability of machine learning algorithms to predict leak-off pressure on the basis of geospatial information of LOP measurements. About 3000 LOP test data were collected from 1800 exploration wells offshore Norway. Three machine learning algorithms (the deep neural network (DNN), random forest (RF), and support vector machine (SVM) algorithms) optimized by three hyperparameter search methods (the grid search, randomized search and Bayesian search) were compared with multivariate regression analysis. The Bayesian search algorithm needed fewer iterations than the grid search algorithms to find an optimal combination of hyperparameters. The three machine learning algorithms showed better performance than the multivariate linear regression when the features of the geospatial inputs were properly scaled. The RF algorithm gave the most promising results regardless of data scaling. If the data were not scaled, the DNN and SVM algorithms, even with optimized parameters, did not provide significantly improved test scores compared to the multivariate regression analysis. The analyses also showed that when the number of data points in a geographical setting is much smaller than that of other geographical areas, the prediction accuracy reduces significantly. Full article

We present a generalized Geertsma solution that can consider any number of finite-thickness layers in the subsurface whose mechanical properties are different from layer to layer. In addition, each layer can be assumed either isotropic or anisotropic. The accuracy of the generalized solution is validated against a numerical reference solution. The generalized Geertsma solution is further extended by a linear superposition framework that enables a response simulation due to an arbitrarily-distributed non-uniform pressure anomaly. The linear superposition approach is tested and validated by solving a realistic synthetic model based on the In Salah CO₂ storage model and compared with a full 3D finite element solution. Finally, by means of a simple inversion exercise (based on the linear superposition approach), we learn that the stiffnesses of cap rock and reservoir are the most influencing parameter on the inversion result for a given layering geometry, suggesting that it is very important to estimate high-confidence mechanical properties of both cap rock and reservoir. Full article

The seabed of the Pontine Archipelago (Tyrrhenian Sea) insular shelf is peculiar as it is characterized by a mixed siliciclastic–carbonate sedimentation. In order to reconstruct the Late Quaternary paleoenvironmental evolution of the Pontine Archipelago, this study investigates the succession of facies recorded by two sediment cores. For this purpose, benthic foraminifera and rhodoliths assemblages were considered. The two cores (post-Last Glacial Maximum in age) were collected at 60 (CS1) and 122 m (Caro1) depth on the insular shelf off Ponza Island. The paleontological data were compared with seismo-stratigraphic and lithological evidence. The cores show a deepening succession, with a transition from a basal rhodolith-rich biodetritic coarse sand to the surface coralline-barren silty sand. This transition is more evident along core Caro1 (from the bottom to the top), collected at a deeper water depth than CS1. In support of this evidence, along Caro1 was recorded a fairly constant increase in the amount of planktonic foraminiferal and a marked change in benthic foraminiferal assemblages (from Asterigerinata mamilla and Lobatula lobatula assemblage to Cassidulina carinata assemblage). Interestingly, the dating of the Caro1 bottom allowed us to extend to more than 13,000 years BP the rhodolith record in the Pontine Archipelago, indicating the possible presence of an active carbonate factory at that time. Full article

This study presents a numerical investigation of the source aspect ratio (AR) influence on tsunami decay characteristics with an emphasis in near and far-field differences for two initial wave shapes Pure Positive Wave and N-wave. It is shown that, when initial total energy for both tsunami types is kept the same, short-rupture tsunami with more concentrated energy are likely to be more destructive in the near-field, whereas long rupture tsunami are more dangerous in the far-field. The more elongated the source is, the stronger the directivity and the slower the amplitude decays in the intermediate- and far-fields. We present evidence of this behavior by comparing amplitude decay rates from idealized sources and showing their correlation with that observed in recent historical events of similar AR. Full article

We describe the opening of back-arc basins and the associated formation of accretionary wedges through the application of techniques of deformable plate kinematics. These methods have proven to be suitable to describe complex tectonic processes, such as those that are observed along the Africa–Europe collision belt. In the central Mediterranean area, these processes result from the passive subduction of the lithosphere belonging to the Alpine Tethys and Ionian Ocean. In particular, we focus on the opening of the Tyrrhenian basin and the contemporary formation of the Apennine chain. We divide the area of the Apennine Chain and the Tyrrhenian basin into deformable polygons that are identified on the basis of sets of extensional structures that are coherent with unique Euler pole grids. The boundaries between these polygons coincide with large tectonic lineaments that characterize the Tyrrhenian–Apennine area. The tectonic style along these structures reflects the variability of relative velocity vectors between two adjacent blocks. The deformation of tectonic elements is accomplished, allowing different rotation velocities of lines that compose these blocks about the same stable stage poles. The angular velocities of extension are determined on the basis of the stratigraphic records of syn-rift sequences, while the rotation angles are obtained by crustal balancing. Full article

The west Florida shelf (WFS; Gulf of Mexico, USA) is an important area for commercial and recreational fishing, yet much of it remains unmapped and unexplored, hindering effective monitoring of fish stocks. The goals of this study were to map the habitat at an intensively fished area on the WFS known as “The Elbow”, assess the differences in fish communities among different habitat types, and estimate the abundance of each fish taxa within the study area. High-resolution multibeam bathymetric and backscatter data were combined with high-definition (HD) video data collected from a near-bottom towed vehicle to characterize benthic habitat as well as identify and enumerate fishes. Two semi-automated statistical classifiers were implemented for obtaining substrate maps. The supervised classification (random forest) performed significantly better (p = 0.001; α = 0.05) than the unsupervised classification (k-means clustering). Additionally, we found it was important to include predictors at a range of spatial scales. Significant differences were found in the fish community composition among the different habitat types, with both substrate and vertical relief found to be important with rock substrate and higher relief areas generally associated with greater fish density. Our results are consistent with the idea that offshore hard-bottom habitats, particularly those of higher vertical relief, serve as “essential fish habitat”, as these rocky habitats account for just 4% of the study area but 65% of the estimated total fish abundance. However, sand contributes 35% to total fish abundance despite comparably low densities due to its large area, indicating the importance of including these habitats in estimates of abundance as well. This work demonstrates the utility of combining towed underwater video sampling and multibeam echosounder maps for habitat mapping and estimation of fish abundance. Full article

In the frame of a Belgo-Haitian cooperation project (PIC 2012–2016), a study of the local seismic hazard was performed in Fond Parisien, an area located on the foothills of the “Massif de la Selle”, along the easternmost portion of the Enriquillo Plantain Garden Fault (EPGF). The H/V Spectral Ratio (HVSR) technique was applied to study the resonance frequency of the target areas and the azimuth of the wave field. The amplification factors were estimated using Standard Spectral Ratios obtained from earthquakes recorded by a temporary seismic network. Using the Multichannel Analysis of Surface Waves method, the seismic properties of the shallow layers were investigated. Then, the results were compared to local Electrical Resistivity Tomography data. These results highlight, in the central part of Fond Parisien, an E-W zone of low velocities ranging from 200 m/s to 450 m/s and low resistivities between 1 Ωm and 150 Ωm, due both to tectonic folding of the rocks and to the presence of sediment filling in the eastern part. The latter is marked, in most of its sites, by resonances at one or more frequencies ranging from 0.7 Hz to 20 Hz. Infiltration and storage of brackish water in the underground layers also contribute to the low resistivity values. With the noise HVSR data, we also evidenced a significant influence of the EPGF on the main orientation of the seismic wavefield as in the vicinity of this fault, the azimuths are parallel to the orientation of the fault. Overall, the results also show greater potential for site effects in the block formed by the sedimentary basin and strong amplification of the seismic ground motion for the sites bordering the basin to the north and west. We interpret the amplification in the north and south-west as probably originating from topographic irregularities locally coupled with sediment deposits, while in the center of the western part, the site effects could be explained by the presence of folds and related weakened and softened rocks. By the integration of several geophysical methods, we could distinguish areas where it is possible to build more safely. These zones are located in the northern part and encompass Quisqueya Park and neighboring areas as well as the village “La Source” in the southern part. In the rest of Fond Parisien, i.e., in the more central and eastern parts, buildings should be erected with caution, taking into account the nearby presence of the EPGF and the influence of fine sediments on the amplification of the seismic motion. Full article

The availability of high-quality surface data acquired by recent Mars missions and the development of increasingly accurate methods for analysis have made it possible to identify, describe, and analyze many geological and geomorphological processes previously unknown or unstudied on Mars. Among these, the slow and large-scale slope deformational phenomena, generally known as Deep-Seated Gravitational Slope Deformations (DSGSDs), are of particular interest. Since the early 2000s, several studies were conducted in order to identify and analyze Martian large-scale gravitational processes. Similar to what happens on Earth, these phenomena apparently occur in diverse morpho-structural conditions on Mars. Nevertheless, the difficulty of directly studying geological, structural, and geomorphological characteristics of the planet makes the analysis of these phenomena particularly complex, leaving numerous questions to be answered. This paper reports a synthesis of all the known studies conducted on large-scale deformational processes on Mars to date, in order to provide a complete and exhaustive picture of the phenomena. After the synthesis of the literature studies, the specific characteristics of the phenomena are analyzed, and the remaining main open issued are described. Full article

We aimed to contribute to a shift in higher education teaching and learning methods by considering problem-based learning (PBL) as an approach capable of positively affecting students from a geology and environment (GE) curricular unit. In a convenience sample from a Portuguese public university, two groups of students were defined: (1) an experimental group (n = 16), to which an intervention program (IP) based on PBL was applied, and (2) a comparison group (n = 17), subjected to the traditional teaching approach. For nine weeks, students subject to the IP faced four problem scenarios about different themes. A triangulation of methods was chosen. The study involved two phases: (1) qualitative (sustained on content analysis of driving questions raised by students, registered in a monitoring sheet) and (2) quantitative (quasi-experimental study, based on data from a prior and post-test knowledge assessment). The qualitative results point to the development of more complex cognitive-level questioning skills after increasing familiarity with PBL. The data obtained in the quantitative study, which included both a “within-subjects” and a “between-subjects” design, show higher benefits in the experimental group, documenting gains in terms of scientific knowledge when using the PBL methodology. Full article

Southern Ozark Mississippi Valley-type ores are enriched in radiogenic Pb, with isotopic signatures suggesting that metals were supplied by two end-member components. While the less radiogenic component appears to be derived from various shale and sandstone units, the source of the more radiogenic component has not yet been identified. Analyses of cherts from the Early Ordovician Cotter Dolomite and tripolitic chert from the Early Mississippian Boone Formation contain highly radiogenic Pb, with isotopic ratios comparable to those of ores. However, most samples have lower ²⁰⁸Pb/²⁰⁴Pb and ²⁰⁷Pb/²⁰⁴Pb for a given ²⁰⁶Pb/²⁰⁴Pb compared to ores. These relationships demonstrate that the enriched Pb isotopic values of the ore array cannot be related to the host and regional lithologies sampled, suggesting that the source of high ratios may lay further afield. The slope of the linear trend defined by the Pb isotope ratios of ores corresponds to an age of about 1.19 Ga. Therefore, an alternative for the linear array is the involvement of Precambrian basement in supplying ore Pb. Rare earth element patterns show that diagenetic processes involving the action of groundwater and hydrothermal fluids affected the sampled lithologies to various degrees, with Cotter Dolomite having experienced the highest degree of alteration. Full article

Cobbles and boulders on the seafloor are of high ecological value in their function as habitats for a variety of benthic species, contributing to biodiversity and productivity in marine environments. We investigate the origin, physical shape, and structure of habitat-forming cobbles and boulders and reflect on their dynamics in coastal environments of the southwestern Baltic Sea. Stone habitats are not limited to lag deposits and cannot be sufficiently described as static environments, as different dynamic processes lead to changes within the physical habitat structure and create new habitats in spatially disparate areas. Dynamic processes such as (a) ongoing exposure of cobbles and boulders from glacial till, (b) continuous overturning of cobbles, and (c) the migration of cobbles need to be considered. A distinction between allochthonous and autochthonous habitats is suggested. The genesis of sediment types indicates that stone habitats are restricted to their source (glacial till), but hydrodynamic processes induce a redistribution of individual cobbles, leading to the development of new coastal habitats. Thus, coastal stone habitats need to be regarded as dynamic and are changing on a large bandwidth of timescales. In general, wave-induced processes changing the physical structure of these habitats do not occur separately but rather act simultaneously, leading to a dynamic type of habitat. Full article

At the very end of the year 2020, on 29 December, a hazardous earthquake of M_w = 6.2 hit the area of Petrinja and its surroundings, in the NW of Croatia. The earthquake was felt across the area of 400 km, leaving an inconceivable damage in the vicinity of the epicenter, devastated towns and ruined lives. In order to map the spreading of earthquake waves and to determine the coseismic ground displacement after the mainshock, we have analyzed open satellite radar images of Sentinel-1 and the GNSS data from the nearest CORS station related to the epicenter, along with the seismic faults. In this paper, we addressed and mapped the displacement linear surface ruptures detected by the SAR interferometry. The results show the vertical ground displacement to the extent of −12 cm in the southern area and up to 22 cm in the north-western part of a wide area struck by the earthquake impact, related to the epicenter. Subsidence and uplift in a range of ±5 cm over a wider affected area indicate a spatial extent and hazardous impact made by the earthquake. The ground displacement of 30 cm to the West and 40 cm to the East has been identified considering the intersection of Pokupsko and Petrinja strike-slip fault system in the seismic zone of Pannonian basin. Accordingly, we obtained matching results of 5 cm south-easting shift and −3 cm subsidence on Sisak GNSS CROPOS station, addressing the tectonic blocks movement along the activated complex fault system. The results compared with the geology data confirm the existence of two main faults; the Pokupsko and the Petrinja strike-slip faults and interpret the occurrence of secondary post-seismic events over the observed area. Full article

Since 2016, Stromboli volcano has shown an increase of both frequency and energy of the volcanic activity; two strong paroxysms occurred on 3 July and 28 August 2019. The paroxysms were followed by a series of major explosions, which culminated on January 2021 with magma overflows and lava flows along the Sciara del Fuoco. This activity was monitored by the soil CO₂ flux network of Istituto Nazionale di Geofisica e Vulcanologia (INGV), which highlighted significant changes before the paroxysmal activity. The CO₂ flux started to increase in 2006, following a long-lasting positive trend, interrupted by short-lived high amplitude transients in 2016–2018 and 2018–2019. This increasing trend was recorded both in the summit and peripheral degassing areas of Stromboli, indicating that the magmatic gas release affected the whole volcanic edifice. These results suggest that Stromboli volcano is in a new critical phase, characterized by a great amount of volatiles exsolved by the shallow plumbing system, which could generate other energetic paroxysms in the future. Full article

Recently, a new strain rate map of Italy and the surrounding areas has been obtained by processing data acquired by the persistent scatterers (PS) of the synthetic aperture radar interferometry (InSAR) satellites—ERS and ENVISAT—between 1990 and 2012. This map clearly shows that there is a link between the strain rate and all the shallow earthquakes (less than 15 km deep) that occurred from 1990 to today, with their epicenters being placed only in high strain rate areas (e.g., Emilia plain, NW Tuscany, Central Apennines). However, the map also presents various regions with high strain rates but in which no damaging earthquakes have occurred since 1990. One of these regions is the Apennine sector, formed by Sannio and Irpinia. This area represents one of the most important seismic districts with a well-known and recorded seismicity from Roman times up to the present day. In our study, we merged historical records with new satellite techniques that allow for the precise determination of ground movements, and then derived physical dimensions, such as strain rate. In this way, we verified that in Irpinia, the occurrence of new strong shocks—forty years after one of the strongest known seismic events in the district that occurred on the 23 November 1980, measuring Mw 6.8—is still a realistic possibility. The reason for this is that, from 1990, only areas characterized by high strain rates have hosted significant earthquakes. This picture has been also confirmed by analyzing the historical catalog of events with seismic completeness for magnitude M ≥ 6 over the last four centuries. It is easy to see that strong seismic events with magnitude M ≥ 6 generally occurred at a relatively short time distance between one another, with a period of 200 years without strong earthquakes between the years 1732 and 1930. This aspect must be considered as very important from various points of view, particularly for civil protection plans, as well as civil engineering and urban planning development. Full article

The abandoned mining district of Linares (South Spain) is marked with waste from the mining and the processing of metal ores that pose an environmental hazard to watercourses. A combined analysis of waste, sediments and water was carried out to analyse the impact of a smelter on Baños Creek. The composition of the facility waste was determined using X-ray diffractometry and scanning electron microscopy. The total contents of the metal(loid)s in the waters and sediments of the watercourse were analysed, and sequential metal(loid) extraction of solid samples was carried out. The facility wastes consisted mainly of secondary minerals, such as natropharmacosiderite and spertiniite, as well as rare metal salts, such as mopungite and NaPb₂(CO₃)₂(OH). The leachates generated by these wastes were highly alkaline, with a pH of 10 and a total dissolved solids concentration of approximately 9 g L⁻¹. This Na-bicarbonate-type water had an As concentration above 200 mg L⁻¹ and elevated levels of Pb, Sb and Zn (5029 µg L⁻¹, 841 µg L⁻¹ and 525 µg L⁻¹, respectively). This highly contaminated lixiviate had a significant effect on the chemical quality of the waters and the bioavailability of metal(loid)s in the creek sediments, especially in the headwaters. In this zone, the As, Pb, Sb and Zn concentrations in the most mobile fraction of the sediments reached 1035 mg kg⁻¹, 261 mg kg⁻¹, 45 mg kg⁻¹ and 30 mg kg⁻¹, respectively. By comparison, smelter slag and mining waste have a much lower impact on the waters and the mobile fraction of the sediments, while significantly increasing the total concentration of these potentially toxic elements in creek sediments. Full article

The Late Holocene was characterized by several centennial-scale climate oscillations including the Roman Warm Period, the Dark Ages Cold Period, the Medieval Warm Period and the Little Ice Age. The detection and investigation of such climate anomalies requires paleoclimate archives with an accurate chronology as well as a high temporal resolution. Here, we present ²³⁰Th/U-dated high-resolution multi-proxy records (δ¹³C, δ¹⁸O and trace elements) for the last 2500 years of four speleothems from Bunker Cave and the Herbstlabyrinth cave system in Germany. The multi-proxy data of all four speleothems show evidence of two warm and two cold phases during the last 2500 years, which coincide with the Roman Warm Period and the Medieval Warm Period, as well as the Dark Ages Cold Period and the Little Ice Age, respectively. During these four cold and warm periods, the δ¹⁸O and δ¹³C records of all four speleothems and the Mg concentration of the speleothems Bu4 (Bunker Cave) and TV1 (Herbstlabyrinth cave system) show common features and are thus interpreted to be related to past climate variability. Comparison with other paleoclimate records suggests a strong influence of the North Atlantic Oscillation at the two caves sites, which is reflected by warm and humid conditions during the Roman Warm Period and the Medieval Warm Period, and cold and dry climate during the Dark Ages Cold period and the Little Ice Age. The Mg records of speleothems Bu1 (Bunker Cave) and NG01 (Herbstlabyrinth) as well as the inconsistent patterns of Sr, Ba and P suggests that the processes controlling the abundance of these trace elements are dominated by site-specific effects rather than being related to supra-regional climate variability. Full article

Charred materials (anthracomass) stored within a soil constitute a major part of its pyrogenic archive and could provide evidence of past fire events, both natural and anthropogenic. However, the dynamics of man-made contributions to the total anthracomass of soil at different time scales are insufficiently understood. In this study, we determined the anthracomass concentrations, stocks, and particle-size distribution in anthropogenically transformed soils of different genesis and ages. Materials were collected from the following archaeological sites within Central Russia—3 Upper Paleolithic sites (Avdeevo, Khotylevo-2 and Yudinovo-1), 2 Early Iron Age settlements (Khotylevo-2 and Yaroslavl), and 1 Medieval site (Yaroslavl). Samples from different cultural layers (CLs), plough layers, and native soils (control) were studied. We identified anthracomass accumulation over a wide chronological scale starting from the Upper Paleolithic Period. The high degree of preservation of anthracomass in ancient anthropogenically transformed soils was explained by the presence of large fragments of charred bones, which are more durable in comparison to wood charcoal. The anthracomass concentrations and stocks in the Early Iron Age plough layer were lower than those in the Medieval plough layer. CLs were generally more enriched in the anthracomass than plough layers, due to their sedimentational genesis, which is more favorable for anthracomass preservation than the turbational genesis of plough layers. However, the differences between charred particle sizes in synlithogenic CLs and turbational plough layers were less clear than expected, due to the specific conditions of formation of each particular layer, e.g., burial rate, duration of ploughing, and type of agricultural land use. Full article

This work addresses the contamination of the sediments of an alluvial plain and riverbed of a tributary of the San Francisco River, in the Brazilian state of Minas Gerais, by potentially toxic elements from an industrial unit of metallic alloys production. This area was subdivided into four areas (A1, A2, A3, and A0 (background area)) where sediment samples have been collected followed by geochemical characterization and spatial distribution of the contaminants. This characterization was based on the (1) analysis of dissolved elements in the interstitial waters, (2) identification of exchangeable and carbonates bounded fractions, and (3) leaching tests using deionized water adjusted to the local pH. This analysis revealed high levels mainly in Cd, Pb, and Zn, in the interstitial waters and in the more soluble phases of sediments. The comparison between the levels of these elements in the leached extracts and the more soluble fractions corroborates the high capacity of these elements to be leached from the alluvium following precipitation episodes. The geochemical characterization and spatial distribution of the contaminants will allow, in the near future, a choice of the most appropriate environmental remediation technique(s) for the environmental requalification of this area. Full article

The existence of iron carbide in the upper mantle allows an assumption to be made about its possible involvement in the abyssal abiogenic synthesis of hydrocarbons as a carbon donor. Interacting with hydrogen donors of the mantle, iron carbide can form hydrocarbon fluid. In order to investigate the role of iron carbide in the abiogenic synthesis of hydrocarbons, the chemical reaction between cementite Fe₃C and water was modeled under thermobaric conditions, corresponding to the upper mantle. A series of experiments were conducted using a high-pressure high-temperature Toroid-type large reactive volume unit with further analysis by means of gas chromatography. The results demonstrated the formation of hydrocarbon fluid in a wide range of thermobaric conditions (873–1223 K, 2.5–6.0 GPa) corresponding to the upper mantle. A strong correlation between the composition of the fluid and the pT conditions of the synthesis was illustrated in the investigation. The higher temperature of the synthesis resulted in the formation of a “poor” hydrocarbon mixture, primarily comprising methane, while a higher pressure yielded the opposite effect, converting iron carbide into a complex hydrocarbon system, containing normal and iso-alkanes up to C₇ and benzene. This correlation explains the diversity of hydrocarbon systems produced experimentally, thus expanding the thermobaric range of the possible existence of complex hydrocarbon systems in the upper mantle. The results support the suggestion that the carbide—water reaction can be a source of both the carbon and hydrogen required for the abyssal abiogenic synthesis of hydrocarbons. Full article

The South-Pyrenean Basin (northeastern Spain) has yielded a rich and diverse record of Upper Cretaceous (uppermost Campanian−uppermost Maastrichtian) vertebrate fossils, including the remains of some of the last European dinosaurs prior to the Cretaceous-Paleogene (K-Pg) extinction event. In this work, we update and characterize the vertebrate fossil record of the Arén Sandstone and Tremp formations in the Western Tremp Syncline, which is located in the Aragonese area of the Southern Pyrenees. The transitional and continental successions of these sedimentary units are dated to the late Maastrichtian, and exploration of their outcrops has led to the discovery of numerous fossil remains (bones, eggshells, and tracks) of dinosaurs, including hadrosauroids, sauropods, and theropods, along with other tetrapods such as crocodylomorphs, testudines, pterosaurs, squamates, and amphibians. In particular, this fossil record contains some of the youngest lambeosaurine hadrosaurids (Arenysaurus and Blasisaurus) and Mesozoic crocodylomorphs (Arenysuchus and Agaresuchus subjuniperus) in Europe, complementing the lower Maastrichtian fossil sites of the Eastern Tremp Syncline. In addition, faunal comparison with the fossil record of Hațeg island reveals the great change in the dinosaur assemblages resulting from the arrival of lambeosaurine hadrosaurids on the Ibero-Armorican island, whereas those on Haţeg remained stable. In the light of its paleontological richness, its stratigraphic continuity, and its calibration within the last few hundred thousand years of the Cretaceous, the Western Tremp Syncline is one of the best places in Europe to study the latest vertebrate assemblages of the European Archipelago before the end-Cretaceous mass extinction. Full article

This project aims to re-assess our understanding of the shape of fluvial bedload gravels by drawing together existing information on fluvial gravel shape. At its crux, however, is the interpretation of a large, high-quality set of measurements made on bedload gravels from the Sabeto River of western Viti Levu, Fiji. This work reveals that the apparent simplicity displayed by most studies of downstream rounding disguises a complex pattern of stepwise reversals to an angular state, the result of the splitting of cobble- and boulder-sized particles. Particle sphericity changes rapidly during the initial stages of transport. Along the Sabeto, this seems to be the result of attrition, with breakage generating the low and continuing presence of low sphericity particles in the system. Elsewhere, however, sphericity is a consequence of shape sorting and we speculate that rivers globally exist along a sorting–attrition continuum. The form of fluvial gravels is not what would be expected were sorting the dominant control on gravel form. Instead measurements of form display a complex relationship with roundness (and thus with breakage and abrasion). Fluvial gravels appear to evolve to a distinctive shape that may offer a means of distinguishing the products of riverine deposition. Full article

The polyphase structural evolution of a sector of the internal Central Apennines, where the significance of pelagic deposits atop neritic carbonate platform and active margin sediments has been long debated, is here documented. The results of a new geological survey in the Volsci Range, supported by new stratigraphic constraints from the syn-orogenic deposits, are integrated with the analysis of 2D seismic reflection lines and available wells in the adjacent Latin Valley. Late Cretaceous syn-sedimentary faults are documented and interpreted as steps linking a carbonate platform to the adjacent pelagic basin, located to the west. During Tortonian time, the pelagic deposits were squeezed off and juxtaposed as mélange units on top of the carbonate platform. Subsurface data highlighted stacked thrust sheets that were first involved into an initial in-sequence propagation with top-to-the-ENE, synchronous to late Tortonian foredeep to wedge-top sedimentation. We distinguish up to four groups of thrust faults that occurred during in-sequence shortening (thrusts 1–3; about 55–60 km) and backthrusting (thrust 4). During Pliocene to recent times, the area has been uplifted and subsequently extended by normal faults cross-cutting the accretionary wedge. Beside regional interest, our findings bear implications on the kinematic evolution of an orogenic wedge affected by far-traveled units. Full article

Constraining the pre-eruptive processes that modulate the chemical evolution of erupted magmas is a challenge. An opportunity to investigate this issue is offered by the interrogation of the crystals carried in lavas. Here, we employ clinopyroxene crystals from back-arc lavas in order to identify the processes driving basalt to andesite magma evolution within a transcrustal plumbing system. The assembled clinopyroxene archive reveals that mantle melts injected at the crust-mantle transition cool and crystalize, generating a clinopyroxene-dominated mush capped by a melt-rich domain. Magma extracted from this deep storage zone fed the eruption of basalt to basaltic andesite lavas. In addition, chemically evolved melts rapidly rising from this zone briefly stalled at shallow crustal levels, sourcing crystal-poor andesite lavas. Over time, hot ascending primitive magmas intercepted and mixed with shallower cooling magma bodies forming hybrid basic lavas. The blended clinopyroxene cargoes of these lavas provide evidence for the hybridization, which is undetectable from a whole-rock chemical perspective, as mixing involved chemically similar basic magmas. The heterogeneity we found within the clinopyroxene archive is unusual since it provides, for the first time, a complete set of mush-related scenarios by which mantle melts evolve from basalt to andesite compositions. Neither the whole-rock chemistry alone nor the record of the mineral phases crystallizing subsequent to clinopyroxene can provide insights on such early magmatic processes. The obtained clinopyroxene archive can be used as a template for interpretation of the record preserved in the clinopyroxene cargoes of basalt to andesite lavas elsewhere, giving insights into the magma dynamics of the feeding plumbing system that are lost when using whole-rock chemistry. Full article

There is an increasing need for quantitative rockfall hazard and risk assessment that requires a precise definition of the terms and concepts used for this particular type of landslide. This paper suggests using terms that appear to be the most logic and explicit as possible and describes methods to derive some of the main hazards and risk descriptors. The terms and concepts presented concern the rockfall process (failure, propagation, fragmentation, modelling) and the hazard and risk descriptors, distinguishing the cases of localized and diffuse hazards. For a localized hazard, the failure probability of the considered rock compartment in a given period of time has to be assessed, and the probability for a given element at risk to be impacted with a given energy must be derived combining the failure probability, the reach probability, and the exposure of the element. For a diffuse hazard that is characterized by a failure frequency, the number of rockfalls reaching the element at risk per unit of time and with a given energy (passage frequency) can be derived. This frequency is relevant for risk assessment when the element at risk can be damaged several times. If it is not replaced, the probability that it is impacted by at least one rockfall is more relevant. Full article

The body of literature seeking to evaluate particle interception in vegetated, aquatic environments is growing; however, comparing the results of these studies is difficult due to large variation in flow regime, particle size, vegetation canopy density, and stem configuration. In this work, we synthesize data from these studies and develop a functional form of particle interception efficiency ($\eta$) as a function of stem Reynolds number (${\mathrm{Re}}_c$), stem diameter, vegetation frontal area, particle–collector diameter ratio, flow velocity, and kinematic viscosity. We develop this functional relationship based on a dimensional analysis and hypothesize that the coefficients would exhibit regimes within different ${\mathrm{Re}}_c$ ranges. We test this hypothesis by synthesizing data from 80 flume experiments reported in the literature and in-house flume experiments. Contrary to our hypothesis, data from different ${\mathrm{Re}}_c$ ranges follow a single functional form for particle interception. In this form, $\eta$ varies strongly with collector density and particle–collector diameter ratio, and weakly with ${\mathrm{Re}}_c$ and particle–fluid density ratio. This work enables more accurate modeling of the flux terms in sedimentation budgets, which can inform ongoing modeling and management efforts in marsh environments. For example, we show that by integrating the new functional form of particle interception into established models of marsh elevation change, interception may account for up to 60% of total sedimentation in a typical silt-dominated marsh ecosystem with emergent vegetation. Full article

In spite of picturesque landscapes, natural beauties and authentic traditional lifestyles to be seen in East Kazakhstan, tourism is far from being developed. The Kazakh Altai (called the Kazakh Switzerland) is one the most colourful parts of the country and, indeed, all Central Asia. The attractiveness of this geographically isolated region (formerly a part of the Imperial Russia), consisting of rocky semi-deserts, vast parkland-steppes, and rugged mountain terrains, is reflected in its distinctive geological and geomorphological character, its pristine nature, and its extraordinary geodiversity and biodiversity. This study presents a roster of geotourism and ecotourism loci for the broader Altai area within a framework of sustainable development. The modelled assessment of the tourism and recreation potential is based upon multi-proxy analyses of GIS, DEM, and cartographic data. It integrates the most appealing natural (biotic and abiotic) site-specific natural features across all physiographic zones within a broad region. The most significant and representative geosites fall within three geographic sectors suitable for geo- and ecotourism. Prospects for travel to these places are enhanced by the presence of numerous prehistoric archaeological sites and historical monuments, which document the rich, multi-ethnic background of Kazakhstan and the ancient Silk Road that traverses it. These geological, environmental and cultural resources, and the regional geoheritage and environmental conservation concepts have been figured into strategies for economic growth of rural Kazakhstan. Visitors travelling to this most appealing region are constrained by climate of pronounced continentality, seasonality, geographic accessibility, the international border-zone regulations and a lack of services of an international standard. Full article

The continental margin is a key area for studying the sedimentary processes related to the advance and retreat of the Ross Ice Shelf (Antarctica); nevertheless, much remains to be investigated. The aim of this study is to increase the knowledge of the last glacial/deglacial dynamics in the Central Basin slope–basin system using a multidisciplinary approach, including integrated sedimentological, micropaleontological and tephrochronological information. The analyses carried out on three box cores highlighted sedimentary sequences characterised by tree stratigraphic units. Collected sediments represent a time interval from 24 ka Before Present (BP) to the present time. Grain size clustering and data on the sortable silt component, together with diatom, silicoflagellate and foraminifera assemblages indicate the influence of the ice shelf calving zone (Unit 1, 24–17 ka BP), progressive receding due to Circumpolar Deep Water inflow (Unit 2, 17–10.2 ka BP) and (Unit 3, 10.2 ka BP–present) the establishment of seasonal sea ice with a strengthening of bottom currents. The dominant and persistent process is a sedimentation controlled by contour currents, which tend to modulate intensity in time and space. A primary volcanic ash layer dated back at around 22 ka BP is correlated with the explosive activity of Mount Rittmann. Full article

Being a country exposed to strong seismicity, the estimation of seismic hazard in Tajikistan is essential for urbanized areas, such as the rapidly growing capital city Dushanbe. To ensure people’s safety and adequate construction work, a detailed seismic microzonation is the key to proper hazard planning. Existing estimations of seismic hazard date back to 1978; they are based on engineering geological investigations and observed macroseismic data. Thereupon relies the Tajik Building Code, which considers seismic intensities according to the Medvedev–Sponheuer–Karnik Scale, MSK-64. However, this code does not accurately account for soil types, which vary considerably in Dushanbe—not only by their nature, but also due to increasing anthropogenic influences. In this study, we performed a series of analyses based on microtremor array measurements, seismic refraction tomography, and instrumental data recording from permanent stations for standard spectral ration and from mobile seismic stations for the horizontal to vertical spectral ratio in order to provide a comprehensive full-cover microzonation of Dushanbe accounting for soil types. Our results identify several critical areas where major damage is likely to occur during strong earthquakes. Full article

Most rock masses contain natural fractures. In many engineering applications, a detailed understanding of the characteristics of fluid flow through a fractured rock mass is critically important for design, performance analysis, and uncertainty/risk assessment. In this context, rock fractures and fracture networks play a decisive role in conducting fluid through the rock mass as the permeability of fractures is in general orders of magnitudes greater than that of intact rock matrices, particularly in hard rock settings. This paper reviews the modelling methods developed over the past four decades for the generation of representative fracture networks in rock masses. It then reviews some of the authors’ recent developments in numerical modelling and experimental studies of linear and non-linear fluid flow through fractures and fracture networks, including challenging issues such as fracture wall roughness, aperture variations, flow tortuosity, fracture intersection geometry, fracture connectivity, and inertia effects at high Reynolds numbers. Finally, it provides a brief review of two applications of methods developed by the authors: the Habanero coupled hydro-thermal heat extraction model for fractured reservoirs and the Kapunda in-situ recovery of copper minerals from fractures, which is based on a coupled hydro-chemical model. Full article

The quality of water is crucial for the qualification of river beaches. The Cávado River watershed (Northern Portugal) contains five river beaches with a regular and specific mandatory monitorization. The main subject of this research is the evaluation of spatial and temporal water microbiological and physicochemical parameters to assess the water quality improvement and consequently watershed management. The results of monitoring surface water, considering microbiological parameters from the five river beaches (2015/19), and physicochemical parameters from three water points along the Cávado River (2018/19) were considered. The river beaches located upstream of the town of Braga has an “excellent” and “good” quality, while the river beach located downstream shows a lower water quality. The physicochemical water results indicated that there is a progressive degradation of water quality from upstream to downstream of the river, which is associated with the influence of domestic and industrial activities. To improve water quality, continuous monitoring will be necessary, with the implementation of adequate awareness-raising programs and strategic water quality management by the population and local agents. Full article