Geosciences, Volume 11, Issue 1 (January 2021) – 37 articles

Passive seismic methods are increasingly used in monitoring unstable rock slopes that are likely to cause rockfalls. Event classification is a basic step in microseismic monitoring. However, the classification of events generated by the propagation of fractures and rockfalls is still uncertain due to their similar features in the time and frequency domains. Hypocenter localization might be a powerful tool to distinguish events generated by fracture propagation from those caused by rockfalls. In this study, a classification procedure based on hypocenter location was validated using a selected subset of high-quality data recorded by a five-geophone network installed on a steep rock slope in Northern Italy. Considering the complexity and heterogeneity of the rock mass, a 3D velocity model that was derived from a tomographic experiment was used. We performed the localization using the equal differential time method. The location results fairly fit our expectations on suspected rockfall events because most signals were located near the rock face. However, only 4 out of 20 suspected fracture events were unquestionably confirmed as fractures being located inside the rock mass and far enough from the rock face. Further improvements in location accuracy are still necessary to distinguish suspected fracture events located close to the rock face from rockfalls. This study demonstrates that hypocenter location is a promising method to improve the final classification of microseismic events. Full article

The present study deals with the characterization of a ceramic assemblage from the Late Mycenaean (Late Helladic III) settlement of Kastrouli, at Desfina near Delphi, Central Greece using various analytical techniques. Kastrouli is located in a strategic position supervising the Mesokampos plateau and the entire peninsula and is related to other nearby coeval settlements. In total 40 ceramic sherds and 8 clay raw materials were analyzed through mineralogical, petrographic and microstructural techniques. Experimental briquettes (DS) made from clayey raw materials collected in the vicinity of Kastrouli, were fired under temperatures (900 and 1050 °C) in oxidizing conditions for comparison with the ancient ceramics. The petrographic analysis performed on thin sections prepared from the sherds has permitted the identification of six main fabric groups and a couple of loners. The aplastic inclusions recognized in all fabric groups but one confirmed the local provenance since they are related to the local geology. Fresh fractures of representative sherds were further examined under a scanning electron microscope (SEM/EDS) helping us to classify them into calcareous (CaO > 6%) and non-calcareous (CaO < 6%) samples (low and high calcium was noted in earlier pXRF data). Here, the ceramic sherds with broad calcium separation are explored on a one-to-one comparison on the basis of detailed mineralogical microstructure. Moreover, their microstructure was studied, aiming to estimate their vitrification stage. The mineralogy of all studied samples was determined by means of X-ray powder diffraction (XRPD), permitting us to test the validity of the firing temperatures revealed by the SEM analysis. The results obtained through the various analytical techniques employed are jointly assessed in order to reveal potters’ technological choices. Full article

In Norway, shallow landslides are generally triggered by intense rainfall and/or snowmelt events. However, the interaction of hydrometeorological processes (e.g., precipitation and snowmelt) acting at different time scales, and the local variations of the terrain conditions (e.g., thickness of the surficial cover) are complex and often unknown. With the aim of better defining the triggering conditions of shallow landslides at a regional scale we used the physically based model TRIGRS (Transient Rainfall Infiltration and Grid-based Regional Slope stability) in an area located in upper Gudbrandsdalen valley in South-Eastern Norway. We performed numerical simulations to reconstruct two scenarios that triggered many landslides in the study area on 10 June 2011 and 22 May 2013. A large part of the work was dedicated to the parameterization of the numerical model. The initial soil-hydraulic conditions and the spatial variation of the surficial cover thickness have been evaluated applying different methods. To fully evaluate the accuracy of the model, ROC (Receiver Operating Characteristic) curves have been obtained comparing the safety factor maps with the source areas in the two periods of analysis. The results of the numerical simulations show the high susceptibility of the study area to the occurrence of shallow landslides and emphasize the importance of a proper model calibration for improving the reliability. Full article

The Cenozoic Malaguide Basin from Sierra Espuña (Internal Betic Zone, S Spain) due to the quality of outcropping, areal representation, and continuity in the sedimentation can be considered a key-basin. In the last 30 years, a large number of studies with very different methodological approaches have been done in the area. Models indicate an evolution from passive margin to wedge-top basin from Late Cretaceous to Early Miocene. Sedimentation changes from limestone platforms with scarce terrigenous inputs, during the Paleocene to Early Oligocene, to the deep basin with huge supplies of turbidite sandstones and conglomerates during the Late Oligocene to Early Miocene. The area now appears structured as an antiformal stack with evidence of synsedimentary tectonics. The Cenozoic tectono-sedimentary basin evolution is related to three phases: (1) flexural tectonics during most of the Paleogene times to create the basin; (2) fault and fold compartmentation of the basin with the creation of structural highs and subsiding areas related to blind-fault-propagation folds, deforming the basin from south to north during Late Oligocene to Early Aquitanian times; (3) thin-skin thrusting tectonics when the basin began to be eroded during the Late Aquitanian-Burdigalian. In recent times some works on the geological heritage of the area have been performed trying to diffuse different geological aspects of the sector to the general public. A review of the studies performed and the revisiting of the area allow proposing different key-outcrops to follow the tectono-sedimentary evolution of the Cenozoic basin from this area. Eight sites of geological interest have been selected (Cretaceous-Cenozoic boundary, Paleocene Mula Fm, Lower Eocene Espuña-Valdelaparra Fms, Middle Eocene Malvariche-Cánovas Fms, Lowermost Oligocene As Fm, Upper Oligocene-Lower Aquitanian Bosque Fm, Upper Oligocene-Aquitanian Río Pliego Fm, Burdigalian El Niño Fm) and an evaluation has been performed to obtain four parameters: the scientific value, the educational and touristic potential, and the degradation risk. The firsts three parameters obtained values above 50 being considered of “high” or “very high” interest (“very high” in most of the cases). The last parameter shows always values below 50 indicating a “moderate” or “low” risk of degradation. The obtained values allow us considering the tectono-sedimentary evolution of this basin worthy of being proposed as a geological heritage. Full article

With evidence suggesting that climate change is resulting in changes within the hydrologic cycle, the ability to robustly model hydroclimatic response is critical. This paper assesses how extreme runoff—1:2- and 1:30-year return period (RP) events—may change at a regional level across the UK by the 2080s (2069–2098). Capturing uncertainty in the hydroclimatic modelling chain, flow projections were extracted from the EDgE (End-to-end Demonstrator for improved decision-making in the water sector in Europe) multi-model ensemble: five Coupled Model Intercomparison Project (CMIP5) General Circulation Models and four hydrological models forced under emissions scenarios Representative Concentration Pathway (RCP) 2.6 and RCP 8.5 (5 × 4 × 2 chains). Uncertainty in extreme value parameterisation was captured through consideration of two methods: generalised extreme value (GEV) and generalised logistic (GL). The method was applied across 192 catchments and aggregated to eight regions. The results suggest that, by the 2080s, many regions could experience large increases in extreme runoff, with a maximum mean change signal of +34% exhibited in East Scotland (1:2-year RP). Combined with increasing urbanisation, these estimates paint a concerning picture for the future UK flood landscape. Model chain uncertainty was found to increase by the 2080s, though extreme value (EV) parameter uncertainty becomes dominant at the 1:30-year RP (exceeding 60% in some regions), highlighting the importance of capturing both the associated EV parameter and ensemble uncertainty. Full article

This paper presents a geospatial methodology for zoning the earthquake-induced soil liquefaction risk at a continental scale and set-up in a Geographic Information System (GIS) environment by coupling data-driven and knowledge-driven approaches. It is worth mentioning that liquefaction is a phenomenon of soil instability occurring at a very local spatial scale; thus, the mega-zonation of liquefaction risk at a continental scale is a hard facing challenge. Since the risk from natural disasters is the convolution of hazard, vulnerability, and exposure, the liquefaction risk mapping is based on the combination of geospatial explanatory variables, available at the continental scale, of the previously listed three assumed independent random variables. First, by applying a prediction model calibrated for Europe, the probability of liquefaction is mapped for the whole continent. Then, the Analytical Hierarchy Process (AHP) is adopted to identify areas that have a high risk of liquefaction, taking into account proxy data for exposure. The maps are computed for different levels of severity of ground shaking specified by three return periods (i.e., 475, 975, and 2475 years). A broad variety of stakeholders would benefit from the outcomes of this study, such as civil protection organizations, insurance and re-insurance companies, and infrastructure operators. Full article

This paper presents a chemical-mineralogical characterization of construction materials from medieval Renaissance buildings of Ferrara (NE Italy) to provide an insight into the nature and provenance of the raw materials used. Biagio Rossetti was an Italian architect and urbanist from the city of Ferrara. From 1483, he was the architect of the Duke of Ferrara Ercole I d’Este who in 1492 assigned him the project of enlarging the city of Ferrara. Biagio Rossetti is still famous because he designed and built many notable palaces and churches in Ferrara, e.g., the Palazzo Roverella, the monastery of Santa Maria delle Grazie and the renovation of the church of San Andrea. To date, only the first two historic buildings are still in use and consequently restored, while the church of San Andrea has been abandoned over the years and the remains have been subject to decay. Different kinds of samples (bricks, cotto, plaster and mortars) were collected from the three sampling sites and analyzed in X-ray fluorescence and X-ray diffractometer to investigate the construction materials through the evaluation of their chemical composition, historic building activity and degradation degree. These investigations should provide knowledge useful for restoration and conservation processes. Full article

From the middle of the 1990s, geotourism was introduced through the first geotrails, their evolution, and the first geopark worldwide in Gerolstein/Vulkaneifel, Germany. The latter is one of the founding members of the European Geoparks Network, which was established in 2000 at the International Tourism Bourse (ITB) in Berlin. The main goal of the first geopark was to link geological heritage with tourism in a rural area that was trying to create new perspectives to inspire young people to stay in their home territory. Geotourism was initiated as part of sustainable tourist development and for future sustainable development at that time in the Gerolstein region. The first steps to implement the Gerolstein/Vulkaneifel Geopark, Germany, were taken in 1992. The core aspects included geological heritage, science transfer, and education as tools for developing geotourism in the broad sense and integrating local people and municipalities in the geopark activities of the rural region. Close collaboration with the local and regional tourism organizations highlighted the need to both define tools and demonstrate their success. Up to now, practice has shown that such success can be demonstrated by the infrastructure created, as well as adjacent measures and activities. A network of factors was determined to play a significant role in ensuring the successful sustainable development in a geopark across the field of geotourism. There are many activities and publications on geological heritage, geosite assessment, significance, and use, but there are fewer which reflect on the network of factors highlighted in this contribution which were first presented in 2002. In many publications and investigations, one factor in particular is emphasized, e.g., infrastructure development, such as panels or other items. Not all of the mentioned factors, however, are being addressed. As a result, a selection of good practice examples of UNESCO Global Geoparks (UGGP) working on the network-oriented conceptual basis has been studied here, in line with the conceptual principle set forth about 25 years ago. The geopark examples in this study include Lesvos Island UGGp (Greece), Naturtejo UGGP (Portugal), Vulkaneifel UGGp (Germany), and Hong Kong UGGP (China), as well as the example of the Messel Pit World Heritage Site (WHS) (Germany). The latter was integrated to present an example which is not a geopark, showing that this concept can also be transferred to a WHS as a tool for sustainable development according the UN 2030 Agenda. The information on the development of the selected examples was obtained by visiting the geopark territory and from the geopark’s websites and published material as a combined methodology. Full article

The transport of dissolved minerals and groundwater flow plays a crucial role in the ecosystem of many wetlands. Nonetheless, installing equipment to monitor groundwater seepage is invasive, harms vegetation, and can impact biodiversity. By remotely mapping surface temperature in late summer, when there is the greatest difference between warm soil and cold groundwater, temperature patterns can expose areas with the greatest upward gradient and flow. The conventional method of using tensiometers to measure hydraulic gradient and estimate flux using Darcy’s law was applied and compared with thermal imaging to characterize groundwater seepage at two contrasting sites within a central North Dakota fen (groundwater discharge wetland). Both sites exhibited variable gradients between the shallow and deep tensiometers. The temperature trend determined from the thermal imaging showed a closer relationship to the measured hydraulic gradients at the herbaceous (Sedge) site than at the wooded (Willow) site. Saturated hydraulic conductivity K ranged from 6 × 10⁻⁵ to 2 × 10⁻⁴ m/s for the Willow site; and 6 × 10⁻⁶ to 1 × 10⁻⁴ m/s for Sedge site. The flux calculated for the Willow site ranged from 1.4 × 10⁻⁵ to 2.7 × 10⁻⁴ m/s and that of the Sedge site ranged from 2.2 × 10⁻⁶ to 6.3 × 10⁻⁵ m/s. The gradients are affected at shallow depth because of heterogeneous soil stratigraphy, which is likely the reason that seepage faces at the sites cannot be mapped solely by thermal imaging. Full article

Seismic sequences are a powerful tool to locally infer geometrical and mechanical properties of faults and fault systems. In this study, we provided detailed location and characterization of events of the 3–7 July 2020 Irpinia sequence (southern Italy) that occurred at the northern tip of the main segment that ruptured during the 1980 Irpinia earthquake. Using an autocorrelation technique, we detected more than 340 events within the sequence, with local magnitude ranging between −0.5 and 3.0. We thus provided double difference locations, source parameter estimation, and focal mechanisms determination for the largest quality events. We found that the sequence ruptured an asperity with a size of about 800 m, along a fault structure having a strike compatible with the one of the main segments of the 1980 Irpinia earthquake, and a dip of 50–55° at depth of 10.5–12 km and 60–65° at shallower depths (7.5–9 km). Low stress drop release (average of 0.64 MPa) indicates a fluid-driven initiation mechanism of the sequence. We also evaluated the performance of the earthquake early warning systems running in real-time during the sequence, retrieving a minimum size for the blind zone in the area of about 15 km. Full article

Amfilochia Bay (Eastern Amvrakikos Gulf, Western Greece), a complex marine area affected by tectonism, was investigated for seabed seepage manifestations and for possible inter-relationships between shallow gas accumulations and hypoxia. For this purpose, an integrated research methodology that combined geophysical, geochemical, and hydrographic surveys was applied. Marine geophysical and bathymetric surveys led to the discovery of a gas-induced pockmark group in the study area. Oceanographic surveying confirmed that the bay is hypoxic/anoxic below a depth of 15 m (dissolved oxygen from ~4 to 60 μM). Very weak CH₄ seepage was detected in correspondence to the pockmark group that seemed to slightly enhance the hypoxic environment locally and close to the seabed, with no effect in shallower waters and the oxycline. Methane isotopic analysis showed variable carbon isotopic composition (from −41‰ to −86‰) which could be either related to differential CH₄ oxidation or mixing between microbial and thermogenic gas. However, the pathway of degassing is clearly related to the fault-controlled pockmark group. A protrusion mound, which has erroneously been reported as a submarine “volcano” since 1876, could be the result of mud volcanism based on the geophysical data of this study. Full article

Airborne LiDAR is a widely accepted tool for archaeological prospection. Over the last decade an archaeology-specific data processing workflow has been evolving, ranging from raw data acquisition and processing, point cloud processing and product derivation to archaeological interpretation, dissemination and archiving. Currently, though, there is no agreement on the specific steps or terminology. This workflow is an interpretative knowledge production process that must be documented as such to ensure the intellectual transparency and accountability required for evidence-based archaeological interpretation. However, this is rarely the case, and there are no accepted schemas, let alone standards, to do so. As a result, there is a risk that the data processing steps of the workflow will be accepted as a black box process and its results as “hard data”. The first step in documenting a scientific process is to define it. Therefore, this paper provides a critical review of existing archaeology-specific workflows for airborne LiDAR-derived topographic data processing, resulting in an 18-step workflow with consistent terminology. Its novelty and significance lies in the fact that the existing comprehensive studies are outdated and the newer ones focus on selected aspects of the workflow. Based on the updated workflow, a good practice example for its documentation is presented. Full article

Preparation of a flood probability map serves as the first step in a flood management program. This research develops a probability flood map for floods resulting from climate change in the future. Two models of Flexible Discrimination Analysis (FDA) and Artificial Neural Network (ANN) were used. Two optimistic (RCP2.6) and pessimistic (RCP8.5) climate change scenarios were considered for mapping future rainfall. Moreover, to produce probability flood occurrence maps, 263 locations of past flood events were used as dependent variables. The number of 13 factors conditioning floods was taken as independent variables in modeling. Of the total 263 flood locations, 80% (210 locations) and 20% (53 locations) were considered model training and validation. The Receiver Operating Characteristic (ROC) curve and other statistical criteria were used to validate the models. Based on assessments of the validated models, FDA, with a ROC-AUC = 0.918, standard error (SE = 0.038), and an accuracy of 0.86% compared to the ANN model with a ROC-AUC = 0.897, has the highest accuracy in preparing the flood probability map in the study area. The modeling results also showed that the factors of distance from the River, altitude, slope, and rainfall have the greatest impact on floods in the study area. Both models’ future flood susceptibility maps showed that the highest area is related to the very low class. The lowest area is related to the high class. Full article

In urban areas, the human influence on the city-ecosystem often results in a Subsurface Urban Heat Island (SUHI), which can be used geothermally. Unfortunately, a model of a SUHI does not consider the geology and hydrogeology of the subsoil. These can vary significantly over short distances, and are of considerable importance for the energy balance. In this work, we calculated the energy and its density stored in the subsoil via a SUHI. For this so-called energy-SUHI (e-SUHI), we evaluated the geology and its physical parameters for the first 20 m below ground level in the German city of Nuremberg and linked them to measured underground temperatures in a GIS application. This approach revealed stored energy of 1.634 × 10¹⁰ MJ within the soil and water for the study area with an area of 163 km² and a volume of 3.26 × 10⁹ m³. It corresponds to an average energy density of 5.0 MJ/m³. The highest energy density of 16.5 MJ/m³ was found in the city center area and correlated well to increases in subsurface temperature. As expected, our model reacts sensitively to thickness changes in the geological layers and the unsaturated zone. Full article

This work critically reviews the Eocene–Oligocene source-to-sink systems accumulating volcanogenic sequences in the basins around the Alps. Through the years, these volcanogenic sequences have been correlated to the plutonic bodies along the Periadriatic Fault System, the main tectonic lineament running from West to East within the axis of the belt. Starting from the large amounts of data present in literature, for the first time we present an integrated 4D model on the evolution of the sediment pathways that once connected the magmatic sources to the basins. The magmatic systems started to develop during the Eocene in the Alps, supplying detritus to the Adriatic Foredeep. The progradation of volcanogenic sequences in the Northern Alpine Foreland Basin is subsequent and probably was favoured by the migration of the magmatic systems to the North and to the West. At around 30 Ma, the Northern Apennine Foredeep also was fed by large volcanogenic inputs, but the palinspastic reconstruction of the Adriatic Foredeep, together with stratigraphic and petrographic data, allows us to safely exclude the Alps as volcanogenic sources. Beyond the regional case, this review underlines the importance of a solid stratigraphic approach in the reconstruction of the source-to-sink system evolution of any basin. Full article

Research into the geometric nests built by white-spotted pufferfish indicated the nest’s potential for flow control and reduction in flow velocity. However, studies to date have only focused on the construction process and behaviour of the male pufferfish. Hence, the form and functions of the unique features of the nest remain unclear. The present study aims to explore the flow features most useful in understanding the habitat conditions of the nest through a combination of photogrammetric reconstructions of the nest features and two-dimensional (2D) computational fluid dynamic simulations. The findings show the role of the nest structure in reducing the flow velocity and shear stress within the nesting site. Analysis of shear stress indicates that male pufferfish build the outer zones of the nest with coarser material that improves the overall shear strength of these areas. The study identified the function of the nest structure in the protection of the eggs through reduction in flow variations and improved aeration. The addition of shell fragments to the nest peaks by the male pufferfish contributes to the resiliency of the nest structure and ensures a stable bed surface at the central zone. Full article

Rapid climate warming at northern high latitudes is driving geomorphic changes across the permafrost zone. In the Yamal and Gydan peninsulas in western Siberia, subterranean accumulation of methane beneath or within ice-rich permafrost can create mounds at the land surface. Once over-pressurized by methane, these mounds can explode and eject frozen ground, forming a gas emission crater (GEC). While GECs pose a hazard to human populations and infrastructure, only a small number have been identified in the Yamal and Gydan peninsulas, where the regional distribution and frequency of GECs and other types of land surface change are relatively unconstrained. To understand the distribution of landscape change within 327,000 km² of the Yamal-Gydan region, we developed a semi-automated multivariate change detection algorithm using satellite-derived surface reflectance, elevation, and water extent in the Google Earth Engine cloud computing platform. We found that 5% of the landscape changed from 1984 to 2017. The algorithm detected all seven GECs reported in the scientific literature and three new GEC-like features, and further revealed that retrogressive thaw slumps were more abundant than GECs. Our methodology can be refined to detect and better understand diverse types of land surface change and potentially mitigate risks across the northern permafrost zone. Full article

Seismically induced soil liquefaction has been documented after the M7.8, 2016 Pedernales earthquake. In the city of Jama, the acceleration recorded by soil amplification yielded 1.05 g with an intensity of VIII to IX_ESI-07. The current study combines geological, geophysical, and geotechnical data in order to establish a geological characterization of the subsoil of the city of Jama in the Manabi province of Ecuador. Then, the liquefaction potential index (LPI) has been evaluated considering the PGA-rock values calculated from deterministic methods applied to nearby geological faults, as well as the soil acceleration records for the city of Jama since the Pedernales megathrust earthquake. The importance of conducting geotechnical evaluations of particular colluvial, alluvial, and floodplain deposits, for which the liquefaction probability profiles have been additionally obtained, may serve as a useful tool for edifices foundations or earthquake resistant designs. Finally, the site response analysis is presented using a linear equivalent analysis, where previously seismic records compatible with the target spectrum have been selected. Hereby, the results of ground surface effects have been compared with the spectra of the Ecuadorian Regulation of Construction (NEC) in the context of local seismic amplification. Full article

Meltwater from glaciers makes significant contributions to general streamflow and provides water for flora and fauna. Continuous glacier monitoring programs enhance our understanding of the impacts of global warming on glaciers and their topographical features. The objective of this study is to measure spatial and temporal changes in Canada’s Columbia Icefield glaciers. This study uses Landsat (TM 5 and OLI) images to delineate glacier extents in the Columbia Icefield between 1985 and 2018. The study also analyzes the retreat of the Athabasca, Castleguard, Columbia, Dome, Saskatchewan, and Stutfield Glaciers. The total area covered by the Icefield in 1985 was 227 km². By 2018, the Icefield had lost approximately 42 km² of its area coverage, representing 18% of its previous coverage. All glaciers in the study region retreated and decreased in area over the study period. The pattern observed in this study is one of general ice loss in the Columbia Icefield, which mirrors patterns observed in other mountain glaciers in Western Canada. Full article

Intense meteorological events are the primary cause of geohazard phenomena in mountain areas. In this paper, we present a study of the intense rainfall event that occurred in the provinces of Lecco and Sondrio from 11 to 12 June 2019. The aim of our work is to understand the effect of local topography on the spatial distribution of rainfall and to attempt the reconstruction of a realistic rainfall field relative to that extreme event. This task represents a challenge in the context of complex orography. Classical rain-gauge interpolation techniques, such as Kriging, may be too approximate, while meteorological models can be complex and often unable to accurately predict rainfall extremes. For these reasons, we tested the linear upslope model (LUM) designed for estimating rainfall records in orographic precipitation. This model explicitly addresses the dependence of rainfall intensification caused by the terrain elevation. In our case study, the available radio sounding data identified the convective nature of the event with a sustained and moist southern flow directed northward across the Pre-Alps, resulting in an orographic uplift. The simulation was conducted along a smoothed elevation profile of the local orography. The result was a reliable reconstruction of the rainfall field, validated with the ground-based rain gauge data. The error analysis revealed a good performance of the LUM with a realistic description of the interaction between the airflow and local orography. The areas subjected to rainfall extremes were correctly identified, confirming the determinant role of complex terrain in precipitation intensification. Full article

Generally seen as natural peripheries, upland landscapes present particular challenges both in terms of living, and of recording past human activity within. LiDAR (light detection and ranging) technology has now brought considerable improvement in our ability to record and map surviving archaeological features, but not necessarily increased our appreciation of local agency. Though the iconic landscape around the Iron Age Dacian capital of Sarmizegetusa Regia (Romania) and its Roman conquest have long caught the attention of specialists and the wider public, both previous research and more recent results from an airborne LiDAR survey leave considerable gaps in our understanding of networking potential across this challenging landscape. Based on LiDAR and satellite-generated high- and mid-resolution topographic data, our paper employs an innovative combination of GIS (geographic information system) spatial analysis tools to examine the spatial relationships between Roman military bases, Dacian targets, and the wider landscape as an integral part of a wider interdisciplinary archaeological research. This helped us formulate and test spatial and historical hypotheses, according to which all known and potential Roman military bases in the study area functioned as part of a system where each contributed individual advantages in securing their domination across the landscape. Our research highlighted the advantages and challenges for Comărnicelu as one of the key Roman logistical nodes, and for the attackers at Şesului and Muncelu working in tandem to besiege and subdue Sarmizegetusa Regia. Our study raises doubts with respect of the fall and destruction of the hillfort at Vârfu lui Hulpe as a result of a Roman siege, making space for alternative political narratives. Ultimately, our findings help build a better understanding of this iconic world heritage landscape and its Roman conquest. Full article

The purpose of the present study is to investigate simultaneously pre-earthquake ionospheric and atmospheric disturbances by the application of different methodologies, with the ultimate aim to detect their possible link with the impending seismic event. Three large earthquakes in Mexico are selected (8.2 Mw, 7.1 Mw and 6.6 Mw during 8 and 19 September 2017 and 21 January 2016 respectively), while ionospheric variations during the entire year 2017 prior to 37 earthquakes are also examined. In particular, Total Electron Content (TEC) retrieved from Global Navigation Satellite System (GNSS) networks and Atmospheric Chemical Potential (ACP) variations extracted from an atmospheric model are analyzed by performing statistical and spectral analysis on TEC measurements with the aid of Global Ionospheric Maps (GIMs), Ionospheric Precursor Mask (IPM) methodology and time series and regional maps of ACP. It is found that both large and short scale ionospheric anomalies occurring from few hours to a few days prior to the seismic events may be linked to the forthcoming events and most of them are nearly concurrent with atmospheric anomalies happening during the same day. This analysis also highlights that even in low-latitude areas it is possible to discern pre-earthquake ionospheric disturbances possibly linked with the imminent seismic events. Full article

Near-source effects can amplify seismic ground motion, causing large demand to structures and thus their identification and characterization is fundamental for engineering applications. Among the most relevant features, forward-directivity effects may generate near-fault records characterized by a large velocity pulse and unusual response spectral shape amplified in a narrow frequency-band. In this paper, we explore the main statistical features of acceleration and displacement response spectra of a suite of 230 pulse-like signals (impulsive waveforms) contained in the NESS1 (NEar Source Strong-motion) flat-file. These collected pulse-like signals are analyzed in terms of pulse period and pulse azimuthal orientation. We highlight the most relevant differences of the pulse-like spectra compared to the ordinary (i.e., no-pulse) ones, and quantify the contribution of the pulse through a corrective factor of the spectral ordinates. Results show that the proposed empirical factors are able to capture the amplification effect induced by near-fault directivity, and thus they could be usefully included in the framework of probabilistic seismic hazard analysis to adjust ground-motion model (GMM) predictions. Full article

Earthquake environmental effects (EEEs) were compiled for the earthquakes of 1626, 1759, 1819, and 1904 in the Fennoscandian Peninsula, northern Europe. The principal source of information was the contemporary newspaper press. Macroseismic questionnaires collected in 1759 and 1904 were also consulted. We prepared maps showing newly discovered EEEs together with previously known EEEs and analyzed their spatial distribution. We assigned intensities based on the 2007 Environmental Seismic Intensity (ESI) scale to 27 selected localities and compared them to intensities assigned based on the 1998 European Macroseismic Scale. While the overall agreement between the scales is good, intensities may remain uncertain due to the sparsity of written documentation. The collected data sets are most probably incomplete but still show that EEEs are not unprecedented cases in the target region. The findings include landslides and rockfalls as well as cascade effects with a risk potential and widespread water movements up to long distances. The winter earthquake of 1759 cracked ice over a large area. This investigation demonstrates that the ESI scale also has practical importance for regions with infrequent EEEs. Full article

The presence of the nonstationarity in flow datasets has challenged the flood hazard assessment. Nonstationary tools and evaluation metrics have been proposed to deal with the nonstationarity and guide the infrastructure design and mitigation measures. To date, the examination of how the flood hazards are affected by the nonstationarity is still very limited. This paper thus examined the association between the flood hazards and the nonstationary patterns and degrees of the underlying datasets. The Particle Filter, which allows for assessing the uncertainty of the point estimates, was adopted to conduct the nonstationary flood frequency analysis (NS-FFA) for subsequently estimating the flood hazards in three real study cases. The results suggested that the optimal and top NS-FFA models selected according to the fitting efficiency in general align with the pattern of nonstationarity, although they might not always be superior in terms of uncertainty. Moreover, the results demonstrated the association and the sensitivity of the flood hazards to the perceived patterns and degrees of nonstationarity. In particular, the variations of the flood hazards intensified with the increase in the degree of nonstationarity, which should be assessed in a more elaborate manner, i.e., considering multiple statistical moments. These advocate the potential of using the nonstationarity characteristics as a proxy for evaluating the evolutions of the flood hazards. Full article

In this paper, a study concerning the soil liquefaction potential in the city of Catania is presented. The stress-based liquefaction analysis framework for cohesionless soil includes a function that describes fundamental aspects of dynamic site response, i.e., the shear stress reduction coefficient, r_d, which depends on several factors (depth; earthquake and ground motion characteristics; dynamic soil properties). Various relationships of r_d are reported in literature because of the importance of assessment of CSR. Herein, new variations of r_d with depth have been obtained using different deterministic earthquake scenarios as input motion. The relationships are based on large numbers of site response analyses for different site conditions. The new relationships obtained have been used for the evaluation of the liquefaction potential in the area of the Catania Harbour. The liquefaction resistance has been evaluated by the horizontal stress index (K_D) from seismic dilatometer Marchetti tests (SDMTs). Various correlations were developed to estimate the CRR from K_D, expressed in form of CRR-K_D curves to differentiate between liquefiable and non-liquefiable zones. In this study three different CRR-K_D curves have been used. Full article

The thermal state of the lithosphere and related geothermal heat flow (GHF) is a crucial parameter to understand a variety of processes related to cryospheric, geospheric, and/or biospheric interactions. Indirect estimates of GHF in polar regions from magnetic, seismological, or petrological data often show large discrepancies when compared to thermal in situ observations. Here, the lack of in situ data represents a fundamental limitation for both investigating thermal processes of the lithosphere and validating indirect heat flow estimates. During RV Polarstern expeditions PS86 and PS118, we obtained in situ thermal measurements and present the derived GHF in key regions, such as the Antarctic Peninsula and the Gakkel Ridge in the Arctic. By comparison with indirect models, our results indicate (1) elevated geothermal heat flow (75 ± 5 mW m⁻² to 139 ± 26 mW m⁻²) to the west of the Antarctic Peninsula, which should be considered for future investigations of ice-sheet dynamics and the visco-elastic behavior of the crust. (2) The thermal signature of the Powell Basin characteristic for oceanic crust of an age between 32 and 18 Ma. Further, we propose (3) that at different heat sources at the slow-spreading Gakkel Ridge in the Aurora Vent Field region might explain the geothermal heat flow distribution. We conclude that in situ observations are urgently required to ground-truth and fine-tune existing models and that a multidisciplinary approach is of high importance for the scientific community’s understanding of this parameter. Full article

This work presents new apatite fission track LA–ICP–MS (Laser Ablation Inductively Coupled Plasma Mass Spectrometry) data from Mid–Late Paleozoic rocks, which form the substratum of the Swiss Jura mountains (the Tabular Jura and the Jura fold-and-thrust belt) and the northern margin of the Swiss Molasse Basin. Samples were collected from cores of deep boreholes drilled in North Switzerland in the 1980s, which reached the crystalline basement. Our thermochronological data show that the region experienced a multi-cycle history of heating and cooling that we ascribe to burial and exhumation, respectively. Sedimentation in the Swiss Jura Mountains occurred continuously from Early Triassic to Early Cretaceous, leading to the deposition of maximum 2 km of sediments. Subsequently, less than 1 km of Lower Cretaceous and Upper Jurassic sediments were slowly eroded during the Late Cretaceous, plausibly as a consequence of the northward migration of the forebulge of the neo-forming North Alpine Foreland Basin. Following this event, the whole region remained relatively stable throughout the Paleogene. Our data show that the Tabular Jura region resumed exhumation at low rates in early–middle Miocene times (≈20–15 Ma), whereas exhumation in the Jura fold-and-thrust belt probably re-started later, in the late Miocene (≈10–5 Ma). Erosional exhumation likely continues to the present day. Despite sampling limitations, our thermochronological data record discrete periods of slow cooling (rates of about 1°C/My), which might preclude models of elevated cooling (due to intense erosion) in the Jura Mountains during the Miocene. The denudation (≈1 km) of the Tabular Jura region and the Jura fold-and-thrust belt (≈500 m) has provided sediments to the Swiss Molasse Basin since at least 20 Ma. The southward migration of deformation in the Jura mountains suggests that the molasse basin started to uplift and exhume only after 5 Ma, as suggested also by previous authors. The data presented here show that the deformation of the whole region is occurring in an out-of-sequence trend, which is more likely associated with the reactivation of thrust faults beneath the foreland basin. This deformation trend suggests that tectonics is the most determinant factor controlling denudation and exhumation of the region, whereas the recently proposed “climate-induced exhumation” mechanism might play a secondary role. Full article

Cold-water coral (CWC) habitats are considered important centers of biodiversity in the deep sea, acting as spawning grounds and feeding area for many fish and invertebrates. Given their occurrence in remote parts of the planet, research on CWC habitats has largely been derived from remotely-sensed marine spatial data. However, with ever-developing marine data acquisition and processing methods and non-ubiquitous nature of infrastructure, many studies are completed in isolation resulting in large inconsistencies. Here, we present a concise review of marine remotely-sensed spatial raster data acquisition and processing methods in CWC habitats to highlight trends and knowledge gaps. Sixty-three studies that acquire and process marine spatial raster data since the year 2000 were reviewed, noting regional geographic location, data types (‘acquired data’) and how the data were analyzed (‘processing methods’). Results show that global efforts are not uniform with most studies concentrating in the NE Atlantic. Although side scan sonar was a popular mapping method between 2002 and 2012, since then, research has focused on the use of multibeam echosounder and photogrammetric methods. Despite advances in terrestrial mapping with machine learning, it is clear that manual processing methods are largely favored in marine mapping. On a broader scale, with large-scale mapping programs (INFOMAR, Mareano, Seabed2030), results from this review can help identify where more urgent research efforts can be concentrated for CWC habitats and other vulnerable marine ecosystems. Full article

In the northern South China Sea, pockmarks are widely distributed on the seabed offshore on the southwestern Xisha Uplift. The mineralogy and geochemistry of the clay minerals and surface sediments from the pockmark field were identified using X-ray diffraction (XRD) analysis and X-ray fluorescence (XRF) analysis to trace the provenance, weathering, and sediment transportation system in the area. The clay minerals are primarily comprised of illite, smectite, kaolinite, and chlorite, showing a distribution of average weight percentages of 35%, 35%, 18%, and 13%, respectively. Based on the surrounding fluvial drainage basins and various transport mechanisms (current or monsoon), illite and chlorite primarily originate from rivers in Taiwan and the Mekong and Red Rivers. Kaolinite primarily originates from the Pearl River, and smectite derived from the Luzon arc system is primarily transported by surface currents with significant influence from the Kuroshio intrusion. Full article