Through analysis of the constituent components of the observation data from 12 geoe-lectric field stations located approximately 600 km form the epicenter of the Jishishan Ms6.2 earthquake, this study reveals that during earthquake preparation and occurrence, the geoelec-tric field at each station exhibit non-uniform variations. These variations are characterized by significant decreases or increases in correlation coefficients between different azimuths (δxy) and between long and short polar distances within the same azimuth (δxx). Additionally, there are discrete or extruding changes observed in the dominant azimuth angles of the telluric field at the stations, ranging approximately from 40° to 90°. These anomalies indicate local site effects and exhibit quasi-synchronous temporal characteristics. The study found that stations located in regions with intense tectonic activities, such as Menyuan and Jingtai, experience more significant anomalous changes influenced by the tectonic activities of the Qaidam-Qilian Mountain block. When multiple stations with anomalous changes are used to determine the earthquake epicenter, earthquakes often occur within or near the intersection area of the predicted ranges of these sta-tions. The specific time of earthquake occurrence shows no clear relationship with the locations or tectonic activities of the stations.

Earthquake-triggered landslides represent a significant seismic-related disaster, posing threats to both the lives and properties of individuals in affected areas. Furthermore, they can result in road and river blockages, as well as other secondary disasters, significantly impacting post-earthquake rescue efforts. Efficient and accurate rapid assessment of high-risk landslide zones carries important implications for decision-making in disaster response and for mitigating potential secondary disasters. This study uses the high-intensity Zone VII to IX of the 5 September 2022 Luding Ms6.8 earthquake as a case study. The "difference method" and the "cumulative displacement method," based on the simple Newmark model, are employed to assess earthquake-triggered landslides. The evaluation results from both methods demonstrate that the areas posing an ex-tremely high risk of earthquake-triggered landslides are predominantly situated on the western side of the Xianshuihe Fault. Actual landslide data verification shows that both models have high predictive accuracy, with the difference method model slightly outper-forming the cumulative displacement method model. Moreover, the study recommends determining threshold values for each landslide hazard interval with significant physical meaning using past strong earthquake data as soon as possible when utilizing the difference method to assess the risk of earthquake-triggered landslides.

Relevant observations suggest that in geothermal power resources and volcanic regions in Earth the mean heat flow is either comparable or greater than the current global mean heat flow of this planet. We have inferred that critical or minimum heat flow required to sustain geothermal provinces (T >150⁰C) is at least 90 mWm^-2 and the same required to sustain active volcanoes in Earth is at least 100 mWm^-2.The volcanism in Venus is found have declined significantly during the past 300 Ma. It is suggested that Earth may lose significant part of its current active volcanoes and geothermal power resources in the next 300-500 Ma if the global mean heat flow in Earth is decreasing at a rate of 5.76 mW per 100 Ma. At the current inferred maximum rate of mantle cooling (20^o C per 100 Myrs) the plate tectonics in Earth is likely to cease within 750 Myrs from now. The critical mean surface heat flux for sustaining plate tectonics in Earth is inferred to be around 50-55 mW/m².

Earthquake forecasting is arguably one of the most challenging tasks in Earth sciences owing to the high complexity of the earthquake process. Over the past 40 years, there has been a plethora of work on finding credible, consistent and accurate earthquake precursors. This paper is a cumulative survey on earthquake precursor research, arranged into two broad categories: electromagnetic precursors and radon precursors. In the first category, methods related to measuring electromagnetic radiation in a wide frequency range, i.e. from a few hz to several MHz, are presented. Precursors based on optical and radar imaging acquired by space borne sensors are also considered, in the broad sense, as electromagnetic. In the second category, concentration measurements of radon gas found in soil and air, or even in ground water after being dissolved, form the basis of radon activity precursors. Well-established mathematical techniques for analysing data derived from electromagnetic radiation and radon concentration measurements are also described with an emphasis on fractal methods. Finally, physical models of earthquake generation and propagation aiming at interpreting the foundation of the aforementioned seismic precursors, are investigated

This work discusses the feasibility of multiple linear regression in predicting water/CO₂ relative permeability using training and testing datasets from two nearby wells, separately, of the Lower Cretaceous Lakota Sandstone, Jurassic Hulett Sandstone, and Pennsylvanian Minnelusa Formation at the Dry Fork Station site. The outcome is promising as the predicted and measured relative permeability data are decently comparable. Yet, whether this approach could be generally applicable needs more delicate models and larger training datasets to be determined.

Numerous investigations have suggested that seismic source parameters and properties of the surrounding medium harbor valuable information regarding alterations in stress fields and medium properties at the focal depth. Monitoring the spatial and temporal evolution of these parameters can yield insights into variations in stress fields or medium properties within the seismogenic zone, offering a critical avenue to overcome the Earth's inaccessible nature. In recent years, modern seismic parameters such as seismic moment, focal mechanism solution, source stress drop, corner frequency, radiated seismic energy, rupture radius, and apparent stress have been increasingly used to characterize source characteristics. The Sichuan region is situated on the southeastern edge of the Qinghai-Tibet Plateau and is known for its strong eastward extrusion and structural deformation. The main fault zones in the area include the Xianshuihe Fault, Anninghe Zemuhe Fault, and Longmenshan Fault. This paper estimates source parameters of small to medium-sized earthquakes (ML 1.5-5.2) in the main fault zones and adjacent areas in Sichuan. A total of 4,310 earthquake stress drop measurements from 2019 to 2023 were analyzed to create a stress distribution image along the fault. This earthquake parameters training catalog provides insight into how coseismic stresses change and their relationship with other geophysical factors, as well as their spatial and temporal evolution.

The relentless advancement and integration of mechanization and intelligence within the coal mining industry, combined with the increasing depths of mining operations, necessitate significantly higher safety standards. This necessity arises due to the increasingly complex geological and coal seam conditions encountered as mining operations delve deeper into the earth. Traditionally, 3D seismic exploration technology has been fundamental in ensuring both operational efficiency and safety in the coal mining sector. However, despite its critical role, this technology still has limitations in detecting and characterizing small-scale geological structures which can pose substantial risks if not accurately identified and mitigated. To address these challenges, this paper introduces a sparsity-constrained adaptive hard thresholding algorithm designed to tackle the spectral inversion optimization problem with an norm constraint for the enhanced retrieval of reflection coefficients, thereby establishing a comprehensive workflow for enhancing seismic data resolution through compressed sensing spectral inversion. Real data validation has demonstrated that this method significantly enhances the resolution of seismic data while maintaining the amplitude integrity and fidelity of the original signals. This improvement allows for a more detailed and accurate characterization of small-scale geological features, which are critical in identifying and mitigating potential hazards within coal mining operations, thus offering robust support for the implementation of more effective safety protocols and disaster prevention strategies in coal mines.

Weathering and erosion transport minerals and organics toward the sea or lake bottoms over geologic time. The finest solids are deposited in lower, less turbulent areas, such as lake bottoms and continental shelves. They sometimes stack up to thicknesses of kilometers, and begin compacting. These sediment sections are called shales, and as initially deposited in water, shales can have porosities up to 50-80{\%} water, As they are buried, many alteration products from oil to slate are produced due to overburden and temperature increasess, making them important to study. A preceeding paper \cite{1} showed that pressure solution is the primary mechanism for porosity reduction, with possible mechanical compaction at shallow depths. Without naming the mineral(s) involved, it postulated that the greater the product of the water and pore interfaces, the faster the reaction would proceed. This term is $ \varphi^{4m/3}(1-\varphi)^{4n/3}$ , where $\varphi$ is porosity and m and n are numbers close to unity. The large exponents, {\em 4/3}, recognize that the reaction occurs at the molecular scale at which the surfaces are rough. A second term, $\exp^{(-E/RT)}$ , indicates that the reaction is impeded by a quantum energy barrier, E, with diminished impeding power as increased available thermal energy, represented by the absolute temperature, T, becomes available at greater depths in the Earth. These two factors combine to allow porosity $\varphi$ to reduce with time, or equivalently for the fraction of solids, $(1-\varphi)$, to increase with time,$\left. \frac{\partial(ln(1-\varphi) ) }{\partial t } \right |_{\sigma} = (\varphi )^{4m/3}( 1-\varphi )^{4n/3}Ae^{-E/(RT)}$. { \em Here it is shown that this equation can cover quantitatively any actual deposition rates which may have been experienced by the six sections studied, the actual deposition rates being unknown for these cases. Hence a time-depth depositional history for any new shale section, known in detail, would allow determination of the parameters m, n, E and a lumped proportionality constant A. This was accomplished by showing that, for a wide range of depositional rates, r, the range of E for any of the studied sections is small compared to laboratory measurements of quartz solution\cite{2}, 24+-15kJ/mol. Results were obtained over this range of r's using the previously determined m and n, and porosity and temperature profiles. The presently existing porosity profiles necessarily incorporate any overpressure or underpressure conditions that may have existed in the past or currently, as the net difference between overburden and pore pressure is a primary driving force for pressure solution. $SiO_2$ usually comprises 20-50{\%} of shales. In conclusion, pressure solution of quartz can account for vertical compaction of shales quantitatively in the studied examples. Separately, a quasi-eqivalent method for discussing reduced porosities in Macron 2 versus Macron 1 as increased vertical stress, rather than additional horizontal stress, is illustrated. Experiments are suggested to clarify the pressure solution mechanism. The role of horizontal forces is discussed. From a global point of view, shale compaction by pressure solution is a quantified step in the cycling of continental rocks, beginning with rifting.

The Qiushuwan Cu-Mo deposit, located in the East Qinling molybdenum belt, which is a major mining site recognized for its enormous reserves of medium-sized molybdenum. In order to enhance our comprehension and ascertain supplementary mineral resources, it is critical to conduct an in-depth investigation of both the deep and peripheral regions. The primary goal of this study is to investigate the Qiushuwan Copper Molybdenum Mine, which is located in Zhenping County, Henan Province. The spread spectrum-induced polarization (SSIP) method is utilized as a means of conducting mining exploratory analysis. Significant findings have been obtained through the collection of SSIP data, geological analysis, and electrical parameter study of rock and ore samples along profile-80. By conducting a meticulous examination of the SSIP data at suitable frequencies, two prominent conductive ore bodies, specifically C2 and C4, have been successfully discovered. C2 exhibits a correlation with granite porphyry and molybdenum veins, while C4 is linked to a skarn deposit that contains a concentration of copper ore. The spatial differences in conductivity within the research area have been precisely characterized using resistivity models obtained from the SSIP data. The precision and dependability of our conclusions have been confirmed through the integration of borehole data and the matching of geological sections with inversion models. The recovered core samples show a diverse range of minerals, containing both valuable ores and minerals that have undergone alteration. Pyrite constantly appears with molybdenite in the majority of samples. The aforementioned results provide a strong foundation for future exploration initiatives and sustainable resource exploitation, enabling informed decision-making in the effective use and administration of natural resources. This study adds to the existing body of knowledge on the Qiushuwan Cu-Mo deposit and improves exploration methods in the region by utilizing the SSIP approach. The findings are a useful asset for geologists, mining experts, and stakeholders, offering essential knowledge for future exploration and sustainable utilization of resources.

Models of subsurface structures are important for successful deposit exploration but are challenged by the need to integrate data from different geophysical methods. In the present study, we evaluated a method of joint inversion in which audio-magnetotelluric (AMT) and dual frequency induced polarization (DFIP) data sets are inverted simultaneously to produce a consistent 2D resistivity model to show a clear image of subsurface structures. To achieve the objectives, we conducted AMT and DFIP surveys first on the same survey line within the Dongjun lead-zinc deposit in inner Mongolia by developing 31 AMT survey sites with an interstation separation of 40 m on a 1440 m survey track and operated in fifty-three frequencies in the range 1–10,400 Hz to record the resistivity distribution of subsurface to depths exceeding 800 m. The same survey setup up was applied to the DFIP method by utilizing the SQ-3C model with a pole-dipole array configuration and conducted measurements at frequencies of 4 Hz and 4/13 Hz. We first separately inverted AMT and DFIP measurements. The two-dimensional (2D) model obtained from AMT data revealed distinct low resistivity anomalies in the middle of the 2D inversion model. In contrast, the DFIP inversion model showed a high resistive body in the same region. In response to the discrepancies observed in the separate 2D inversion models, we implemented a joint inversion for both the AMT and DFIP datasets. The joint inversion resistivity model shows surficial conducting bodies and a high conductive body along the profile with relatively high Percent Frequency Effect (PEF) indicating high chargeability. The final joint inversion resistivity model clearly images the large silica alteration zone and the Pb-Zn mineralization. This study demonstrates the feasibility of a joint inversion methodology and highlights the value of integrating geophysical methods through joint inversion for enhanced characterization and exploration of lead-zinc ores.

The Parecis Basin is one of the largest intracratonic basins in Brazil. Despite its considerable size, the basin has yet to be extensively explored, with only five wildcat wells drilled for hydrocarbons. So far, no commercial discoveries have been announced. Regional studies have suggested Paleozoic sedimentation, while recent studies indicate a Neoproterozoic fill. However, no tectonic model accurately describes the basin’s development, and no detailed structural map for the entire basin has been published. The present work presents a new detailed structural map of the Parecis Basin based on a four-step interpretation workflow integrating seismic and gravimetric data. The first step includes converting the public 2D seismic lines to the depth domain. The second step is estimating the residual Bouguer anomaly, where the computed residual anomalies should relate to the basin’s tectonic features. The third step comprises the 2D forward modeling of the gravimetric anomalies using the 2D seismic interpretation as a constraint. The final step compiled all the interpreted features into our new structural map. This map reveals the top of the basement, forming a complex framework of horsts and grabens. Normal faults define the main structural style in the basin. Further, we could recognize thick, high-density bodies embedded in the crystalline basement. These bodies consist of Orosian-Calimian (1.8-1.6 Ga) mafic and ultramafic rocks, which may be a potential source for hydrogen exploration in the basin. Subsequent geophysical and geochemical surveys will assess the hydrogen potential in the area.

Successful geophysical exploration works in Gucheng-Yaxi area located in Gaochun District, Jiangsu province, China, have been limited partly due to complex geological condition of the area and high artificial noise in data acquired using electrical and electromagnetic methods. In this work, we deployed the new anti-interference spread spectrum induced polarization method (SSIP) and audio magnetotellurics (AMT) method to detect a copper-polymetallic deposit in the area. Two-dimensional (2D) inversion results at Gucheng-Yaxi Area revealed a high chargeability anomalous zone on SSIP profile that coincided with a zone of moderate resistivity located between two resistor bodies on AMT profile. A follow up 1200 meters drill hole was established at this high chargeability-moderate resistivity zone which encountered polymetallic (copper, lead, zinc, gold and silver) mineralization at the depth ≥ 400 m. Drill hole data analysis showed that mineralization occurs interspaced in marble rock mass at varies depth. Furthermore, several low resistivity-weak chargeability sections were revealed and inferred attributed to Cretaceous sediments and faults. These faults are thought to have played critical role in the polymetallic mineralization genesis. In summary, this work demonstrated the successful of application of SSIP and AMT at detecting a metallic deposit in an area with high artificial noise. Hence the geophysical prospection potential of Gucheng-Yaxi area is great.

Modern climatic changes have an impact on the bearing capacity of permafrost soils at the base of the foundations of buildings and structures in the urbanized territories of the Arctic and Subarctic. The activation of cryogenic processes leads to the destruction of the infrastructure of human settlements and, as a result, to social, economic, environmental consequences for the population. Based on the results of geothermy of frozen and thawing soil, GPR profiling of the territory of the city of Anadyr, it was concluded that the main risks of permafrost degradation are associated with the spread of hydrogenic melt zones. Maps of soil temperature in imaginary cross-sectional with a depth of 3, 5 and 10 m were compiled; maps of the capacity of thawing soils and the permafrost aquifer; a map of zones of dangerous spread of exogenous cryogenic processes. A permafrost monitoring system in the territory of Anadyr has been substantiated, which will consist of an automated network of observations of soil temperature in 35 wells at the border and at the epicenter of 20 zones of dangerous development of exogenous cryogenic processes and 12 control GPR profiles at the intersection of linear hydrogenic taliks.

The discovery of water ice on the moon has revolutionized our understanding of lunar science and has profound implications for future space exploration and habitation. Through a comprehensive analysis of remote sensing data, isotopic studies, and neutron detection experiments, this research elucidates the presence, distribution, and origins of lunar water ice. Radar observations from missions like Clementine and Chandrayaan-1 reveal distinct polar regions with elevated circular polarization ratios, indicative of volatile ice deposits. Spectral detections by instruments such as M3 and Diviner further confirm the presence of water ice in lunar-shadowed regions, with estimates suggesting reserves of up to 300 million tons at the poles. Neutron detection experiments on lunar orbiters like LRO provide additional evidence, detecting hydrogen signals consistent with water ice below the lunar surface. Isotopic analyses offer insights into the origin and evolution of lunar water ice, with comparisons to cometary and asteroidal sources. Proposed missions, such as Artemis-3 and Artemis-4, aim to collect samples from the lunar South Pole for detailed isotopic analysis, furthering our understanding of lunar water ice's composition and history. Scientific data, including the lunar bulk water content and estimated water ice reserves, underscore the significance of lunar water ice as a critical resource for future space exploration endeavors. Calculations estimating the total water-ice area and water abundance near specific lunar features provide valuable insights into the potential utilization of lunar water ice for sustained human presence on the Moon and beyond.

The aquifer delineation using a hydrogeophysical survey was carried out in De Aar, a small town in the eastern part of the Karoo supergroup, Northern Cape, South Africa. In this re-gion, the current approach for estimating recharge, water allocation plans, and water monitoring does not consider the existence of different types of aquifers. Instead, it assumes there is one aquifer system and therefore the area might have one value of recharge. This study aimed to delineate aquifers into unconfined and confined aquifers for estimating recharge and demonstrate the use of geophysics techniques in delineating aquifer systems. Six VES stations and seven profile lines were measured using a Wenner array configuration and ground magnetotelluric method, respectively. VES data was processed using IPI2win software while electric potential difference curve was in-terpreted together with lithological data. Four significant layers of low (0.9-8.1 Ωm), moderate (22.4-125 Ωm), and high 68-177 Ωm) resistivity values were delineated suggestive of the different lithological formations. The subsurface profile mapping using ground magentotelluric methods did not suggest groundwater accumulations below 50 m, instead supported shallow groundwater potentials boreholes. These results have provided a double layer aquifer system thereby suggesting that recharge estimation need to be carried out in two different aquifer systems to enhance water allocation plans, and water monitoring in groundwater-dependent areas similar to prevailing conditions in the current study area.

Low-frequency electromagnetic induction (EMI) is a non-invasive geophysical method that is based on the induction of electromagnetic (EM) waves into the subsurface to quantify changes in electrical conductivity. In this study we present an open (design details and software are accessible) and modular system for the collection of EMI data. The instrument proposed allows to adjust the separations between the transmitter and up to four receiving antennas as well as the acquisition frequency (in the range between 3 and 50 kHz) to permit measurements with variable depth of investigation. The sensor provides access to raw data and the software described in this study allows control of the signal processing chain. The design specifications permit apparent conductivity measurements in the range of between 1 mS/m and 1000 mS/m, with a resolution of 1.0 mS/m and with a sampling rate of up to 10 samples per second. The sensor allows for a synchronous acquisition of a time stamp and a location stamp for each data sample. The sensor has a mass of less than 5 kg, is portable and suitable for 1-person operation, provides 4 hours of operation time on one battery charge, and provides sufficient rigidity for practical field operations.

Deep underground laboratories offer advantages for conducting high-precision observations of weak geophysical signals, benefitting from a low background noise level. It is both valuable and feasible to enhance strong, noisy ground electromagnetic (EM) field data using synchronously recorded underground EM signals, which typically exhibit a high signal-to-noise ratio. In this study, we propose an EM field reconstruction method employing a Long Short-Term Memory (LSTM) recurrent neural network with referenced deep underground EM observations. Initially, a deep learning model was developed to capture the time-varying features of underground multi-component EM fields using the LSTM recurrent neural network. Subsequently, this model was applied to process synchronously observed strong, noisy data from other conventional observation systems, such as those at surface, to achieve noise suppression through signal reconstructions. Both theoretical analysis and practical observational data suggest that the proposed method effectively suppresses noise and reconstructs clean EM signals. This method is efficient and time-saving, representing an effective approach to fully utilizing the advantages of deep underground observation data. Furthermore, this method could be extended to the processing and analysis of other geophysical data.

The integration of magnetic observatories of many countries into the Intermagnet network has provided researchers with new opportunities to investigate many geophysical processes. The joint processing of measurement evidence from magnetic observatories located in various geolocations makes it possible to better understand the dynamics of the interaction of solar plasma fluxes with the Earth's magnetosphere. For the correct interpretation of the dynamics of physical processes, the accuracy of the time reference of measurement results is of great importance. The measurement results of local observatories were transformed to a unified ecliptic coordinate system oriented towards the Sun in order to extract small variations in the magnetic field occurring before the magnetic storm commencement. Before transforming the coordinates, a constant component equal to the average daily value of the magnetic field components was removed from the results of minute-by-minute measurements. Owing to the averaging evidence from 84 magnetic observatories of the Intermagnet network, variations in the magnetic field of cosmic origin, preceding the magnetic storm commencement were detected. The correlation analysis of these variations recorded by different observatories revealed the presence of time shifts between the data series of different observatories. Only 62% of the observatories had a relative data shift within the range of ± 1 minute. Before the joint data processing in order to efficient use of data recorded with a time shift a method for correcting time shifts has proposed.

The escalation of mining activities in the karst regions of Guizhou Province has heightened the occurrence of water-inrush incidents in deep coal mines. This study focused on water-inrush phenomena within the Xinhua mining area of Jinsha County, Guizhou Province, aiming to investigate the sources of these incidents. The findings indicated that the overlying limestone of the Changxing Formation in the coal seam served as a vulnerable aquifer under certain conditions, leading to water inrushes. The analysis of the spatiotemporal distribution patterns of water-inrush incidents at the working face indicated that previous mining operations damaged the shallow Changxing Formation limestone, resulting in the accumulation of goaf water and the formation of numerous mining-induced fractures. These fractures served as rapid conduits for water inrushes from both atmospheric precipitation and underground sources at the deep working face. The examination of surface water and mine water quality demonstrated that both exhibited similar characteristics, predominantly featuring bicarbonate, sulfate, and sodium compositions. Investigation into the relationship between mine water inflow and atmospheric precipitation established that atmospheric precipitation influenced the mine water supply cycle, with a replenishment period of ~10 months during the operational phase of the Jinyuan Coal Mine and about one month post-closure. The fractures induced by mining activities within the Changxing Formation limestone facilitated water flow, with atmospheric precipitation serving as the primary water source for the mine. This study offered a valuable scientific foundation for addressing water-related damage resulting from atmospheric precipitation in mines susceptible to water inrushes under analogous hydrogeological conditions.

Deep learning has increasingly been employed to identify minerals. However, deep learning can only be used to identify minerals in the distribution of the training set, while any mineral outside the spectrum of the training set is inevitably categorized erroneously within a predetermined class from the training set. To solve this problem, the study introduces the approach that amalgamates One-Class Support Vector Machines (OCSVM) with the ResNet architecture for the out-of-distribution mineral detection. Initially, ResNet undergoes training using a training set comprising well-defined minerals. Subsequently, the earlier layers of the trained ResNet are employed to extract the discriminative features of the mineral under consideration. These extracted mineral features then become the input for OCSVM. When OCSVM discerns the mineral in the training set's distribution, it triggers the subsequent layers within the trained ResNet, facilitating the accurate classification of the mineral into one of the predefined categories encompassing the known minerals. In the event OCSVM identifies the mineral out of the training set's distribution, it is unequivocally categorized as an unclassified or 'unknown' mineral. Empirical results substantiate the method's capability to identify out-of-distribution minerals while concurrently maintaining a commendably high accuracy rate for the classification of the 36 in-distribution minerals.

The Ayers Rock of Australia which is known for its charming colors at different time is actually the Celestial Body Ring of the planet Earth that implies the great Celestial Body Love happened on Earth long long ago. The Tasmania island of Australia which is known for the World Heart is actually in the state of face to face with the main island of Australia, which implies the image spokesperson for the inner core (the heart) of the planet Earth. In other words, with the outline of the main island of Australia, the basic key information of the inner core of the planet Earth could be easily obtained by deducing logically, just like the similarity between the fists and heart of human body. Similar to the mechanism of the blood circulation dominated by the heart of human body, it is dominated and motivated by the inner core of the Earth which forms the circulation of mantle magma between the crust and the inner core, and what leads to the formation of expanding and active area on the equatorial and middle-lower latitude area and the silent and convergent area on the North and South Pole central area. Therefore, the mantle magma whirlpools under the crust on the North and South Pole are supposed to exist logically. The walking sideways of the crabs just implies the celestial body language of the motive trend of the main island of Australia based on their similar outline.

The Huichang fault Basin in Jiangxi Province exhibits high terrestrial heat flow values and a dense concentration of hot springs. However, existing research is limited to individual geothermal fields without a comprehensive basin-wide analysis of geothermal water origins. This study investigates the source of geothermal water across five representative geothermal fields using hydrochemical components and isotopic analyses. The findings indicate: (1) The hydrochemical type transitions from bicarbonate sulphate-sodium-calcium in the north to bicarbonate sodium-calcium in the south; (2) Major chemical components of the geothermal water derive from evaporative salt rocks (Cretaceous sand conglomerate cover) and siliceous formations (granite/metamorphic rock basement); (3) The recharge elevation for geothermal water is approximately 1 km, sourced from adjacent hills. The central basin's geothermal water has an apparent age of about 26,000 years, with δ¹³C indicating participation in deep circulation; (4) Thermal reservoir temperatures range from 74 to 113℃, with the lowest temperatures at the basin's edges (Anzishan and Shangjin geothermal fields) and higher temperatures towards the center, suggesting increased warming and a circulation depth of 1.5 to 2.5 km; (5) Simulation results of reverse water-rock interaction indicate that mineral precipitation-dissolution processes are governed by the basin's main thermal control structures. The geothermal fields are predominantly influenced by reduction alterations associated with granite and are closely linked to paleo-geothermal waters in the enclosed basin, with the formation of these fields heavily reliant on the heat production or conduction from radioactive elements in granite.

A joint analysis of solar wind, geomagnetic field, and earthquake catalog data showed that before the catastrophic M7.8 and M7.5 Kahramanmaras earthquake sequence on February 6, 2023, the closed strong magnetic storm occurred on November 7, 2022, SYM/H=-117 nT. Storm started at 08:04 UT. At this time, the high latitudinal part of Turkey's longitudinal region of future epicenters was located under the polar cusp, where the solar wind plasma would directly access the Earth's environment. The time delay between storm onset and earthquake occurrence was ~91 days. We analyzed all seven strong (M7+) earthquakes from 1967 to 2020 to verify the initial findings. A similar pattern has been revealed for all events. The time delay between magnetic storm onset and earthquake occurrence varies from some days to some months. To continue these investigations, a retrospective analysis of seismic and other geophysical parameters just after preceded geomagnetic storms in the epicenter areas is desirable.

Carbon dioxide capture and storage (CCS) is a necessary requirement for high-emitting CO2 industries to significantly reduce volumes of greenhouse gases released to the atmosphere and mitigate the climate change. Geological CO2 storage into depleted oil and gas fields is the fastest and most accessible option for CCS deployment allowing for re-purposing of existing infrastructure and utilizing the significant knowledge about the subsurface acquired during the field production operations. Location of such depleted fields in neighbourhood of high-emitting CO2 industries is additional advantage available for matured on-shore Europe fields. Considering these advantages, the oil and gas operators are now evaluating different possibilities for CO2 sequestration projects for the fields approaching end of production. This article describes an integrated approach to reservoir simulations focused on evaluating CO2 injection pilot at one of such matured fields operated by MND and located in the Czech Republic. The potential CO2 injection site in focus is a naturally fractured carbonate reservoir. This oil-bearing formation has a gas cap at the top and a limited aquifer at the bottom of the reservoir and was produced mainly by pressure depletion with quite limited pressure support from water injection. The article summarizes results of efforts taken by the multi-disciplinary team to develop and apply an integrated approach to reservoir simulation started from geological modelling of the naturally fractured reservoir and integrating results of laboratory studies and their interpretation (geomechanical and geochemistry laboratory experiments) and dynamic field data analysis (pressure transient analysis including time-lapse) into the history matching of the reservoir model and simulations of pilot CO2 injection. Attention is given to evaluating and modelling stress-sensitive fracture properties and safe injection envelope preventing induced fracturing as well as impact of potential salt precipitation in the near wellbore area. These effects are considered in the context of pilot CO2 injection and addressed in reservoir simulations of the injection scenarios. Single porosity and permeability reservoir simulations with dominating fracture flow and black-oil formulation with CO2 simulated as a solvent were used in this study and arguments for such choice of the simulation approach for the site in focus are shortly discussed. The reservoir simulations have indicated larger site injection capacity than required for the pilot injection and allowed to evaluate CO2 migration pathway within the reservoir. The application of the approach to the site in focus has also revealed large uncertainties, related to fracture description and geomechanical evaluations, resulted in an uncertain safe injection envelope. These uncertainties should be addressed in further studies in preparation to the pilot. The article concludes with an overview of the outcomes of the integrated approach and its application to the field in focus including discussion of issues and uncertainties revealed.

Mexico is a well-known seismically active country, which is primarily affected by several tectonic plate interactions along the southern Pacific coastline and by active structures in the Gulf of Baja California. In this paper we investigate this seismicity in the context of a non-extensive statistical approach based on Tsallis entropy. The analysis is performed using data from the corrected Mexican seismic catalogue provided by the National Seismic Service, spanning a period from January 2000 to October 2023. The Gutenberg-Richter law fitting to the earthquake sub-catalogues for all six regions studied indicates magnitudes of completeness between 3.30 and 3.76. All six regions display values of the entropic index in the range 1.52≲q≲1.61, which are lower than the previously estimated range 1.54≲q≲1.70 using catalogue data from 1988 to 2010. The cause of this discrepancy is certainly due to the use of modern recording devices, which are sensitive enough to detect a larger number of low-magnitude events compared to older instrumentation.

Sea level rise due to global warming is a continuing and, disappointingly, accelerating process, which has already affected and will further impact coastal lowlands and the social and economic activities in these areas. Delos Island, situated in the middle of the Cyclades in the Aegean Sea, was considered the most sacred of all islands in ancient Greek culture and a trading hub for the entire eastern Mediterranean. Uninhabited since the 7th century AD, and consistently the focus of research and touristic attention, the island is designated as an archaeological site and inscribed on the Unesco World Heritage List. Previous studies on the relative sea level (rsl) changes suggest a steadily rsl rise during the last 6,300 years, starting from a sea level at -4.80±0.20m in the Late Neolithic. The seafront of the ancient city of Delos is subject to the effects of the rsl rise, which have caused significant coastline retreat and exposure to the northerly winds and waves, whereas parts of the coastal lowland where the remains of the ancient city lie are inundated, forming extended wetlands. The future impacts of rsl rise on the seafront of ancient Delos are illustrated on very high-resolution digital surface models, evaluating both the flooding risk under different climatic projections, as provided by the IPCC AR6 report, and the ongoing land subsidence, as recorded by GNSS data. An rsl rise ranging from 87cm (SSP1-2.6 scenario) to 148cm (SSP5-8.5 scenario) is anticipated by 2150, requiring both resilience strategies and adaptation solutions as well as mitigation policies to cope with the effects of climate change.

Earthquake-triggered landslides are one of the most significant hazards worldwide. These landslides, involving small or large volumes, can develop in debris flow or avalanche with high mobility and long-runout distance. The motion processes of this kind of landslide also quite complicated. This study reproduced the motion processes of two landslides using the LS-RAPID program. The Aso Bridge landslide occurred on April 16, 2016, with rapid and long-runout distance characteristics. The Kataragai landslide occurred on April 9, 2018, travelling with flow-like behavior. Through field investigations and laboratory experiments, features and mechanical parameters of these two landslides were obtained. In the LS-RAPID program, the pre-failure models and motion processes of these two landslides were established and reproduced. Results indicated that sliding speed and travelling morphology at different intervals were considered as two indicators to describe their motion process. However, while sliding speed of these two landslides have significant differences with sliding speeds of actual landslides, the final morphologies after failure were consistent with those observed in the field investigation.

On the afternoon of January 31, 2022, the northwest and central part of the city of Quito (Ecuador) was hit by a mud-flow that originated in one of the ravines of the Pichincha volcanic complex. According to the official reports, there were 28 fatalities, one missing person, 52 injured, and 53 families affected. The mudflow was triggered by unusual and heavy rainfall in the upper part of the drainage. However, the analysis of the deposits allowed us to establish the existence of anthropogenic elements as crucial components of a local disaster that affected daily life in the city for more than a week. Based on the preliminary data, we were able to determine that the mud and water flowed through different streets of the city, reaching a maximum distance of 3.2 km from the overflow point (collector, construction works to dam the drainage under the city). The analysis of the multimedia material published on social networks and in the press also made it possible to detect the change in the behavior of the flow during its evolution, which went from a hyper-concentrated flow to muddy water. A similar event occurred in 1975, and although it was a larger event (52000 m3, seven times larger than in 2022), it caused two deaths and little damage in the same area. The main causes of the disaster included urban development that did not meet the needs of the area and a lack of awareness or ignorance about the management and intervention of drainage in Quito City.

Climate change is increasing air temperatures and altering the precipitation and hydrological re-gime on a global scale. Challenges arise when assessing the impacts of climate change on the local scale for water resources management proposes, especially for low mountain headwater catch-ments that not only serve as important water towers for local communities, but also have distinct hydrological characteristics. Until now, no low flow or hydrological drought studies have been carried out on the Lauter River. This study is unique in that it compares the Lauter River, a transboundary Rhine tributary with a nearby station on the River Rhine just below its confluence at the French-German border. The Lauter catchment is a mostly natural, forested catchment, however, its water course has been influenced by past and present cultural activities. Climate change disturbances cascade through the hydrologic regime right down to the local scale. As we are expecting more low flow events, the decrease in water availability could cause conflicts be-tween different water user groups in the Lauter catchment. However, the choice of different methods for identifying low flow periods may cause confusion for local water resources manag-ers. Using flow rate time series of the Lauter River between 1956 and 2022 we compare for the first time three low flow identification methods: the variable threshold method (VT), the fixed threshold method (FT), and the Standardized Streamflow Index (SSI). Similar analyses are applied and compared to the adjacent Maxau station on the Rhine River for the same time period. This study aims at 1) interpreting the differences amongst the various low flow identification methods and 2) revealing the differences in low flow characteristics of the Lauter catchment compared to the Rhine River. It appears that the FT reacts faster to direct climate or anthropogenic impacts whereas VT is more sensitive to indirect factors such as decreasing subsurface flow which is typical for small headwater catchments, such as the Lauter where flow dynamics react faster to flow disturbances. Abnormally low flow during the early Spring in tributaries such as the Lauter can help predict low flow conditions of the Rhine River during the following half year and espe-cially for the summer. The results could facilitate Early Warning of hydrological droughts and drought management for water users in the Lauter catchment and further downstream along some of the Rhine.

The Double U-pipe ground heat exchanger, recognized for its straightforward process, cost-effectiveness, high heat exchange efficiency, and low thermal resistance, stands as the predominant ground heat exchanger type in today's shallow geothermal energy development and utilization. Recent years have seen significant research into the heat transfer influencing factors and heat exchange performance of Double U-pipe ground heat ex-changers via experimental testing methods. Yet, studies integrating numerical simulation with in-situ testing have been less frequent. Drawing on the cylindrical heat source model theory and the outcomes of regional in-situ thermal response tests, this paper develops a Double U-pipe ground heat transfer model by establishing physical, mathematical, and heat transfer geometric models. It assesses the impact of varying inlet temperatures, flow rates, and initial ground temperatures on heat exchange efficiency under heating conditions. The results validate the precision of the double U-pipe ground heat exchanger model based on in-situ testing. They show that enhancing the temperature differential between the inlet and initial temperatures, elevating the initial ground temperature, and moderately increasing the flow rate can boost the system's heat exchange efficiency.

Existing research indicates that mineralization of Carlin-type gold deposit is related to geological structures, but current studies mostly focus on the qualitative association, and lack quantitative analysis of the relationship between mineralization and structural features, which limits the potential for predictive mineral exploration. In this study, the Nibao Carlin-type gold deposit was selected as the research area to conduct three-dimensional geological modeling, spatial analysis, and three-dimensional prospectivity modeling. Spatial analysis revealed that the majority of gold ore bodies are situated near the main fault, with gold mineralization primarily occurring within a spatial range of approximately -80 meters to 20 meters from the main fault plane. This spatial distribution predominantly aligns in a SW-NE direction. Spatial characteristic extraction of fault slope indicators suggests that the slope angles are conducive to the migration and accumulation of ore-forming hydrothermal fluids, with ore body metal content rapidly increasing within an angle range of approximately 20° to 60°. The shape analysis of the mainfault indicates that the principal spaces of mineralization enrichment are within -10m to 20m of the fault plane undulations, where gold resources are concentrated in a relatively flat area. However, the degree of mineralization is more pronounced in the protruding sections than in the depressed areas. The areas with significant dip changes in faults display higher variations in orebody attitudes, with mineralization noticeably weaker in areas where the fault steepness changes markedly. Ore-bearing units exhibit significantly higher grades within the -17° to 3° interval of fault steepness transition. This interval also sees a rapid accumulation of gold resources. The Nibao gold mine's mineralization is derived from deep, ore-bearing hydrothermal fluids ascending along major fractures, precipitating, and accumulating in favorable ore-hosting spaces (faults, anticline cores) and lithological combinations (sedimentary tuffs, tuffs). Three deep potential gold targets were identified, highlighting the deep western edge of the F1 fault and the deep southeast-trending lower plate as prospective areas for future exploration. This study underscores the effectiveness of 3D prospective modeling as an indispensable tool for elucidating the spatial distribution of fault-controlled Carlin-type gold deposits and their tectonic determinants.

The research is devoted to the study of geotectonic and geodynamics of the Shu-Ile belt and petrographic composition and intrusive complexes geochemical features analysis. The Shu-Ile suture is an important geological structure of Kazakhstan. It is characterized by the presence of various intrusive formations and deep faults. The material composition of these intrusions is of interest for understanding the geological and metallogenic processes in this region. The research used modern scientific methods, including petrographic and geochemical studies of various compositions intrusion samples. Based on the results of material composition studies of intrusive complexes, their possible connections with other geological structures of Kazakhstan and the metallogeny of the region have been established. Based on the data obtained, a conclusion was made about the intrusive complexes origin and evolution of the Shu-Ile belt, which may be important for further geological research and forecasting of mineral resources in the region.

Identification and segmentation of rock particles is very important in geology, mining, engineering construction, and environmental monitoring. Traditional methods for recognizing the rock particle characteristics are time-consuming, labor-intensive, and limited by the geographical scope of the sampling sites, which makes it more difficult to represent the spatial distribution of the rock particles accurately. Recent advancements in deep learning techniques offer promising solutions to automate this process. This study aims at utilizing state-of-the-art (SOTA) Mask R-CNN deep learning algorithm from the Detectron2 framework, to capture global features of rocks that will enable the recognition of rock particles more efficiently, just from having the rock images. These images are subject to image processing techniques such as a Gaussian filter, to denoise and smoothen the image, and an Illumination Adaptive Transformer (IAT) framework, to adjust the lighting exposure; all of which add to improve the quality of the image dataset. Afterward, the Mask R-CNN and Detectron2 models are utilized in training the images using transfer learning. Experimental results from evaluating the performance of the proposed algorithms showcase the effectiveness of the approach, thereby highlighting its potential to revolutionize rock particle segmentation across various domains.

Surface–subsurface soil-rock modeling is crucial for infrastructure design and groundwater yield optimization, especially in the complex crystalline basement terrains of Penang Island, Malaysia. This study conducted large-scale characterization of soil-rock profiles using ERT, SRT, rock quality designation, and soil penetration test (SPT N-values) data, optimized through regression modeling. The approach effectively identified soil-rock structures and their suitability to prevent structural failures and enhance groundwater productivity. Correlations between borehole lithologic logs and velocity–resistivity models revealed distinct soil compositions across the study area. The results highlighted thick, saturated, and loose silty to sandy bodies in the eastern to northern sections, contrasting with sandy compositions and penetrative fractures in the southern part. The study also found good correlations between rock mass quality and N-values for different soil types. Depths of intra-bedrock weathered/fractured units varied between 12 and >35 m. This has significant implications for determining optimal foundation depths in the study area. Groundwater productivity was associated with intra-bedrock planes of weakness at depths exceeding 40 m. Overall, this study developed empirical relations between geophysical and geotechnical parameters for wet tropical granitic terrains, filling critical gaps in understanding subsurface characteristics.

Two strongly curved areas of the Carpathians in Romania are earthquake-prone: Vrancea - in the South-Eastern Carpathians, and Danubian – in the western edge of the Southern Carpathians. By far the highest seismicity episode in 2023 was recorded in the Southern Carpathians, near Targu Jiu city, of more than 4000 earthquakes. A first ML 5.2 mainshock was recorded on 13 February, followed by a second one, ML 5.7, on 14 February. Our purpose is to provide a first analysis of this sequence in order to explain its extent and unusual development, in time. The possible role of fluids in the middle crust, is discussed in connection with the rather unusual character of the space-time distribution of aftershocks. A large and high-quality dataset is provided by the dense array of seismic stations which has been gradually improved by adding new stations in the epicentral area. The thorough revision of the phase picking and location reveals a few interesting and intriguing features of seismicity that still remain worthy to be further investigated. The database built within this paper provides excellent information not only for Targu Jiu area, but for any other investigation in intra-continental regions characterized by diffuse and sporadic seismic activity.

Bulk composition, mineralogy and microchemistry of various fractions of the filling of the Crater No. 4 at Emmerting were investigated. Here we present mainly the fine fractions, which, in contrast to pebbles larger than ca. 4 cm, are relatively poor in material affected by melting. A 0.1 mm long meteorite fragment without ablation features was found, and remelted fragments of the same body have been identified. The projectile was dominated by enstatite, also containing forsterite, basic plagioclase, Ni-rich pyrrhotite and troilite, kamacite, silica, and probably merrillite. The molten fragments contain porous Mg-rich glass (indicating temperature > 1500 °C), minor Fe oxides, and occasionally uncommon Ni-rich minerals. The bulk composition and electron microprobe show that the amount of meteoritic material is low (<< 1 %), and it is not unequivocally manifested even in platinum-group elements abundances and osmium isotope ratios (the latter, however, were only measured in two samples). This can be explained by dilution of the filling by particles brought after the crater formation, and by probably achondritic composition of the projectile. Relatively large hollow spheroid Fe oxide particles were observed in thin sections from the magnetic fraction. A Mn oxide hollow spherule has also been found. In magnetically separated filling from the nearby Crater No. 5, bulk chemical composition points to contamination with several metals and metalloids probably due to influence of recent alluvial sediments. Therefore, the extent of meteoritic contamination in this crater may be even smaller.

In recent decades, as eutrophication in inland and coastal waters (ICWs) has increased due to anthropogenic activities and global warming, so has the need for timely monitoring. Compared with traditional sampling and laboratory analysis methods, satellite remote sensing technology can provide macro-scale, low-cost, and near real-time water quality monitoring services. The Geostationary Ocean Color Imager (GOCI), integrated onboard the Communication Ocean and Meteorological Satellite (COMS) from the Republic of Korea, was the first geostationary ocean color observation satellite and was operational from April 1, 2011 to March 31, 2021. Over ten years, GOCI has observed oceans, coastal waters, and inland waters within its 2,500 km×2,500 km target area centered on the Korean Peninsula. The most attractive feature of GOCI, compared with other commonly used watercolor sensors, was its high temporal resolution (1h, eight times daily from 0 UTC to 7 UTC), providing the opportunity to monitor the quality of ICWs, where optical properties can change rapidly throughout the day. This systematic review aims to comprehensively review GOCI features and applications in ICWs, analyzing progress in atmospheric correction algorithms and water quality monitoring. Analyzing 123 articles from the Web of Science and China National Knowledge Infrastructure (CNKI) through a bibliometric quantitative approach, we examined GOCI’s strength and performance with different processing methods. These articles reveal that GOCI played an essential role in monitoring the ecological health of ICWs in its observation area in East Asia. GOCI has led the way to a new era of geostationary ocean satellites, providing new technical means for monitoring water quality in oceans, coastal zones, and inland lakes. We also discuss the challenges encountered by geostationary satellites in monitoring water quality and provide suggestions for improvements.

Depressions containing partially molten rocks were investigated in Quaternary pebbly sedi-ments in Alpine foreland in Bavaria, Germany. Ruins of limekilns from 18-19th centuries are abundant around Seeshaupt. They are usually associated with earthworks for fuel feeding and protection from wind. Depressions contain partially molten stones from furnace walls, with lime and burnt loam remnants. On the surface of pebbles (typically gneisses), green glass formed with chemical composition possibly affected by moderately volatile components of wood biomass, mainly K and Cu. Thermoplastic deformation of pebbles is limited and the sur-face melt deformation points to slow flow. At Emmerting, on the other hand, two walled circu-lar depressions which lack additional earthworks have been identified as possible impact cra-ters. Thermally affected pebbles show much more intense deformation. Colorless to greenish glass, also K and Cu enriched, usually forms only thin surface coatings, in places covered by younger melt extruded from pebble interior. Burnt loam is missing. At Grabenstätt-Kaltenbach, a depression of unclear origin was investigated. Morphologically, it matches craters rather than limekilns, but anthropogenic contamination (brick-like material) is probable. The rare molten samples do not display extreme deformation features. The thermally affected stones, mainly those from Emmerting, have elevated magnetic susceptibility.

The catastrophic failure of two dams in Libya in the night of 10th and 11th September 2023, resulted in a significant amount of water and debris destroying a large part of the city of Derna, which is located downstream of the two dams. The flood caused more than 6000 fatalities and the dis-placement of thousands of people from the city. The dams failed due to heavy rainfall caused by the Medicane Daniel, resulting in intense precipitation and the dams overspilling after reaching full capacity. This paper analyses the consequences of the dam breach by modelling the flow discharge and the flooding of the city of Derna. The analysis compares the computed quantities with the findings of the post-event satellite image analysis from UNOSAT and Copernicus, the field reports, and information gathered from the social media through AI detection. The results allows to better understand the dynamics of the dam breach, its initial parameters and an assessment of the accu-racy of the results.

The gold deposits of the eastern region of Cuba are genetically related to the lithological formations of the island arc and the ophiolitic complex, which have been studied and exploited since the time of the Spanish colonization in the mid-sixteenth century. These deposits are part of the so-called Aguas Claras-Guajabales mineral field, in the Holguín Province (Cuba) and lie in an extensive elongated strip of approximately 15 km in length. The object of this work is to make an orderly, detailed, and chronological review of the geological and mining work carried out in this region as well as to highlight the degree of study achieved in each of the stages of research. To facilitate this work, an extensive bibliographic review of all available materials was planned and carried out, such as reports, manuals, monographs, and articles. The results obtained from this study highlight that gold mineralization is closely linked to metasomatic processes produced by the circulation of hydrothermal fluids that affected different lithologies of both the volcanic island arc and the ophiolitic complex, it is noteworthy that the highest gold contents are mainly where the lithologies contact each other. This highlights the fact that the presence of alterations such as serpentinization, listvenitization, rodingitization and propylitization have played a primary role in the deposition of gold during mineralization processes. This work could be a very useful guide for future research in this region, as it is an orderly, concrete, and practical compilation of the characteristics of the mineralization and alterations, as well as a precise indication of the spatial position, thicknesses, and contents of the gold horizons.

Analyses of the three predominant types of sedimentation – clastic, carbonate and black shale – revealed their discontinuity/continuity in different facies environments of the Middle and Late Devonian of the central part of the Volga-Urals petroleum province, namely (1) in depressions, both condensed and non-condensed, (2) on the slopes of depressions and (3) in shallow water environments represented by tectonic arches or localised biohermal uplifts. The regional zonation of conodonts, correlated with the global chronostratigraphic scheme, substantiates the stratigraphic completeness of the general Devonian sequence in this area. In particular, the sedimentation of organic-rich black shales, traditionally considered as source rocks, lasted more than 29 million years, starting from the Late Eifel and ending at the Devonian-Carboniferous boundary. There was only one hiatus in the deposition of both black shales and carbonate sediments, lasting about 2.5 million years at the end of the Frasnian. That interval is dominated by clastic sediments, which are considered to be the most productive reservoirs of the Volga-Ural province. In the time interval from the beginning of the Eifel to the end of the Frasnian (about 17 Ma), (a) clastic sediments show the highest accumulation rates; (b) the most complete sequences occur in depressions; and (c) arches and uplifts show the most significant gaps in the geological record. During Famennian, (a) clastic sedimentation has completely ceased; (b) differences between the slopes of troughs and shallow-water environments have diminished, resulting in stable accumulation of carbonate sediments in both settings; (c) the decrease in sequence thickness (up to 50 m) in depressions suggests reduced accumulation rates of deep-water organic-rich black shales. The results help to understand the complex spatial relationship between different coexisting sedimentation types and focus attention to the incompleteness of the geological record in different environments.

Variance-driven decomposition based on the singular spectrum analysis of the European Project for Ice Coring in Antarctica (EPICA) dD, CO2 and CH4 records allowed a novel quantitative structural interpretation of all glacial/interglacial cycles and glacial terminations of the last 800 kyr. This bottom-up approach used the response components of EPICA stacked records to reconstruct the envelope of the thermal response through a physical interference model. The aim was to improve the understanding regarding the features of intensity, amplitude, and asymmetry of 73 marine isotope stages (MIS) and seven terminations. The Antarctic stack record can be described by a variance-weighted interference model of ten thermal waves of different origins (mid-term oscillation, orbitals and suborbitals) that stochastically interfere at a given time according to their relative differences in frequency, intensity, and polarity. Interglacial/glacial stages resulted from constructive interference and bipolar amplification of warming/cooling responses, respectively. The low intensity MISs (including 90% of substages) and the unbiased-dated terminations fell in the low interference regions where dominant destructive patterns minimise the thermal envelope. The positive skewness of the EPICA stack resulted from the constructive interference with a strong bias in the warming direction, especially after the Mid-Brunhes Event. Duration analysis of short eccentricity hemicycles exhibited an intrinsic unexpectedly prolonged mean cooling in the nominal solution (5.8 kyr) and its EPICA response as well (8.6 kyr), along with an interference-induced asymmetry (21.1 kyr). The overall effect resulted in the saw-tooth shape of glacial cycles which is heavily induced by interference and related to the intermediate feedback mechanisms controlled by short eccentricity.

Under the same seismic load, different local topography will have different effects, which is an important research field in seismic engineering. Hills are connected with valleys is a very typical terrain. Taking this terrain as an example, its response to the most destructive SH waves in the earthquake is studied. By simplifying and dividing the model, the wave functions satisfying the wave equation and boundary conditions in each region are constructed by using the complex function and moving coordinate method. The equations satisfying the boundary continuity conditions are established. With the help of calculation software, the numerical results are obtained. The influence of different parameters on ground motion is discussed.

In this study, we define screening maps for Italy that classify sites based on their potential for triggering landslides. To this end, we analyze seismic hazard maps and hazard disaggregation results on a national scale considering four spectral periods (0.01s, 0.2s, 0.5s, and 1.0s) and three return periods (475, 975, and 2475 years). First, joint distributions of magnitude (M) and distance (R) from hazard disaggregation are analyzed by means of an innovative approach based on image processing techniques to find all modal scenarios contributing to the hazard. In order to obtain the M-R scenarios controlling the triggering of earthquake-induced landslides at any computation node, mean and modal M-R pairs are compared to empirical curves defining the M-R bounds associated with landslide triggering. Three types of landslides are considered (i.e., disrupted slides and falls, coherent slides, and lateral spreads and flows). As a result, screening maps for all of Italy showing the potential for triggering landslides based on the level of seismic hazard are obtained. The maps and the related data are freely accessible at the following web address: https://distav.unige.it/rsni/slideq.php.

The Ballynoe barite deposit is a conformable, mineralized horizon of Lower Carboniferous age overlying a diastem and mass faunal extinction demarking the transition from a quiet water environment to one of dynamic sedimentation. The geometry of the barite orebody correlates with the palaeo-topography of the footwall, which acted as an important control over the lateral extent, thickness and nature of the mineralization. Sedimentary features within the barite horizon suggest that it was precipitated in the form of a cryptocrystalline mud which underwent major diagenetic modification resulting in extensive stylolitisation, recrystallisation and remobilization. There is abundant and compelling geological and isotopic evidence for early local exhalation from the presence of a hydrothermal vent fauna consisting of a delicately pyritized worm tubes hosted and haematized filaments of apparent microbial origin. The worm tubes are remarkably similar to examples from modern and ancient volcanic-hosted massive sulphide deposits, and the filamentous microfossils have similarities to modern Fe-oxidizing bacteria. Strontium in the barite has a 87Sr/86Sr ratio indistinguishable from sea-water between 350 and 344Ma whilst oxygen isotopes from barite and chert suggest a diagenetic origin in equilibrium such sea-water around 60-70°C. Fluid inclusion studies have shown that, in general, low temperature inclusions are very saline (20-25%) whilst at higher homogenization temperatures they are more dilute (9-12%).

Antimony (Sb) has gained significance as a critical raw material (CRM) within the European Union (EU) due to its strategic importance in various industrial sectors, particularly in the textile industry for flame retardants and as a component of Sb-based semiconductor materials. Moreover, Sb is emerging as a potential alternative for anodes used in lithium-ion batteries, a key element in the Energy transition. This study focused on exploring the feasibility of identifying and quantifying Sb mineralizations through the spectral signature of soils using reflectance spectroscopy, a non-invasive remote sensing technique, and by employing deep learning algorithms such as Convolutional Neural Networks (CNNs). Common signal preprocessing techniques were applied to the spectral data, and the soils were analyzed by Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Despite achieving high R-squared values, the study faces a significant challenge of generalization of the model to new data. Despite the limitations, this study provides valuable insights into potential strategies for future research in this field.

The Sanaga River Valley is a deep-seated structure of about 220 km in width, bounded on its northern and southern sides by 025-0550 fault zones. Its internal structure has been very little described before this study. It is characterized by the presence of surface features where a thick stratum of shale and argillaceous sandstones have been documented. The use of electrical resistivity tomography for near-surface imaging and hydrogeological prospecting within shale succession is widespread. However, several methods of spatial modeling of resistivity contours are generally proposed for the interpretation of inverse resistivity model sections. To determine the most appropriate method for interpreting resistivity tomographies in the context of subsoil marked by significant heterogeneity, logarithmic contour intervals, and linear contour interval processing were implemented on data collected over five areas. The data consists of outcrops data, cuttings samples from groundwater boreholes, and resistivity sections ranging in length from 240 to 1000 m. Well-control of the inverse resistivity model sections derived from both processes shows that the linear contour intervals modeling gives a more realistic picture of the reservoir size of the wa-ter-bearing sand in both shale and argillaceous sandstone succession. The sedimentary architecture over Sanaga fault breccia is marked by vertical intrusion anomalies up to 80 m high. The shape of these anomalies, as well as the very low resistivity values that characterize them, suggests that they could be related to mud intrusion.

Due to the long duration and high cost of traditional pumping tests, the response of well water level to seismic wave, earth tide, and barometric pressure provides a feasible method for determining continuous changes in aquifer hydraulic parameters. Aiming at the problem that the response phase shift of well BB water level to tidal tide M2 is greater than that of tide O1, this paper preferentially calculated the time series changes of S, Be, n and BKu based on the response mechanism of well water level to barometric pressure and earth tide and with the help of smooth moving method. Then by using the mixed flow model, the variation of transmissivity and leakage coefficient over time were simultaneously obtained, and the evolution process and characteristics of aquifer parameters near well BB caused by the Wenchuan earthquake were ultimately discussed. The calculation results are basically consistent with the previous pumping test, and have been verified and compared by using other scholars' methods. This solving process avoids the problems such as excessive dependence on initial values, multiple solutions, and unstable tide O1, which has a promoting effect on the study of the impact of seismic activity on aquifer systems.

GNSS is nowadays a standard tool for monitoring and studying the Earth’s dynamic environment. However, the development of dense GNSS measurements remains limited in many experiments by the cost of geodetic-class equipments allowing to achieve the high precision required by many applications. Recently, multi-constellation, multi-frequency, low-power and, above all, less expensive GNSS electronic chips have become available. We present a prototype of a low-cost, open-source multi-GNSS station. Our prototype comprises a dual-frequency GNSS chip, a calibrated antenna, a Raspberry Pi card and a 4G key for data transmission. The system is easy to deploy in the field and allows precise positioning in real time and post-processing. We assess the performance of our prototype in terms of raw data quality, and quality of the obtained high rate and daily position one-year-long time series. Our results demonstrate a quality equivalent to geodetic-class equipment and better quality than other low-cost systems proposed so far.

The nephrite belt in Altun Mountain-Western Kunlun Mountain region spanning approximately 1300 km in Xinjiang, NW China, constituting the world's largest preeminent nephrite deposit. The Qiemo region in the Altun Mountain is a crucial nephrite-producing area in China with demonstrating substantial prospects for future explorationl. While extant research has extensively study on secondary deposits in the Karakash River and native black nephrite deposits in Guangxi Dahua, a lacuna exists in the comprehensive investigation of black nephrite from original deposits in Xinjiang. Margou nephrite represents a recent discovery of black-toned nephrite deposits in Qiemo County, Xinjiang, warranting a comprehensive investigation into this primary source of black nephrite within the region. This scholarly inquiry deploys a suite of sophisticated analytical techniques, encompassing Polarizing microscope, electric microprobe, Backscattered Electron Image, X-ray Fluorescence, Inductively Coupled Plasma Mass Spectrometry. The rigorous experimental test is dedicated to elucidating the chemical and mineralogical composition, and further clarify the type of its genetic types of black-black cyan nephrite from Margou deposit in Qiemo, Xinjiang. Results delineate the mainly composed of tremolite-actinolite, characterized by Mg/(Mg + Fe2+) ratios ranging from 0.86 to 1.0. Minor minerals include diopside, epidote, pargasite, apatite, zircon, pyrite and iron hydroxide. Bulk-rock rare earth element patterns exhibit distinctive features such as negative Eu anomalies (δEu=0.00~0.17), decreasing light rare earth elements, a relatively flat distribution of heavy rare earth elements, and low with total REE concentrations (1.64~38.9 μg/g), concurrently, the Cr (6~21 μg/g) and Ni (2.5~4.5 μg/g) contents are conspicuously low. The magmatic influence of granite emerges as pivotal factor in the genesis of magnesia skarn hosting the Margou nephrite. The distinctive black-black cyan colour is attributed to heightened iron content, mainly associated with FeO (0.08~6.29%). The discerned chemical composition characteristics posit the Margou nephrite as an exemplar of the magnesia skarn-type nephrite deposit.

With the increase of mining range from shallow coal seam outcrops to deep With the increase of mining range from shallow coal seam outcrops to deep mining and the improvement of mining intensity in the mining area, there has been a phenomenon of no water inrush in shallow coal mines and abnormal water inrush in deep coal mines in karst landform areas of Guizhou.The project conducted a study on the water source of water inrush in the Xinhua mining area of Jinsha County, Guizhou Province. The research results indicate that the overlying Changxing Formation limestone of the Longtan Formation is a weak (conductive) aquifer, and there is no possibility of water inrush in normal stratigraphic blocks. The spatiotemporal distribution characteristics of the water inrush working face have verified that the Changxing Formation limestone in the shallow mined area has been damaged by mining activities, resulting in a large number of cracks in the rock mass and the formation of a large-scale water accumulation area, becoming a water source for the water inrush in the working face below and a fast supply channel for atmospheric precipitation. The analysis of surface water and mine water quality data shows that the water quality types of mine water and surface water are basically the same, and they are HCO3•SO4-Na type water. The correlation curve between mine water inflow and atmospheric precipitation verifies that the cycle of atmospheric precipitation supplying the mine in Jinyuan Coal Mine is about 10 months before pit closure, and 1 month after pit closure; The Changxing Formation limestone serves as a water-conducting channel, and atmospheric precipitation serves as a source of replenishment for mine water. The research results provide scientific basis for the treatment of atmospheric precipitation water damage in sudden water inrush mines (areas) under similar hydrogeological conditions.

Due to the ongoing global warming on the Earth, permafrost degradation has been extensively taking place, which poses a substantial threat to civil and industrial facilities and infrastructure elements, as well as to the utilization of natural resources in the Arctic and high-latitude regions. In order to prevent the negative consequences of permafrost thawing under the foundations of constructions, various geophysical techniques for monitoring permafrost have been proposed and applied so far: temperature, electrical, seismic and many others. At the same time, non-stationary electromagnetic methods seem to have significant potential and a number of important practical advantages. We propose a cross-borehole exploration system for the highest localization of target objects in the cryolithozone. A novel mathematical apparatus for three-dimensional modeling of transient electromagnetic (TEM) signals by the vector finite element method has been developed. The original combination of the latter, the Sumudu integral transform and artificial neural networks makes it possible to examine spatially heterogeneous objects of the cryolithozone with a high contrast of geoelectric parameters, significantly reducing computational costs. We have studied the sensitivity of the TEM signals to the boundary between thawed and frozen rocks depending on the inter-borehole distance and spatial orientations of the transmitters and receivers. Numerical simulation results of the TEM monitoring of industrial facilities located on permafrost are considered. The formation of a talik has been shown to significantly manifest itself in the measured electromagnetic responses, which enables timely prevention of industrial disasters and environmental catastrophes.

Аннотация . Периодически развиваются самоподобные движения структурно-вещественных комплексов (СМК) из состава Солнечной системы в течение периода около 4,5 миллиардов лет земной коры, развития Земли в условиях самоорганизации. Синхронные или близкие во времени ряды генетических или парагенетических групп магматических, рудно-минеральных, метасоматических, осадочных образований, их структур и других сопутствующих им процессов, включают в себя каждый комплекс (ГМК). Каждый комплекс в его пределах ограничен во времени и пространстве, по причине признаков самопроизвольного динамического процесса; распространяется на территорию Земли в текущий исторический период; формирует состав и структуру земной коры на этом этапе, определяя ее полициклическое развитие. Количество циклов около 14,

The Badenian/Sarmatian boundary in the Central Paratethyan basins is characterized by a change from open marine conditions during the late Badenian to the assumed brackish conditions during the early Sarmatian. Foraminiferal and palynological results of the Badenian/Sarmatian boundary interval in the Babczyn 2 borehole (SE Poland) showed that the studied interval accumulated under variable, unstable sedimentary conditions. The Badenian/Sarmatian boundary as correlated with a sudden extinction of stenohaline foraminifera is interpreted as due to the shallowing of the basin. The lack of foraminifera and palynomorphs just above the Badenian/Sarmatian boundary can reflect short-term anoxia. The composition of euryhaline assemblages characteristic for the lower Sarmatian part of the studied succession indicates marine to hypersaline conditions.

Abstract: Recent years demonstrate the increased attention to radon from two scientific directions. After neglecting radon as earthquake precursors in 1990-th it becomes again the subject of earthquake-forecast papers discussions due to growing networks of the radon monitoring in different countries, especially, the technologies of real-time radon measurements where the gamma-spectrometers become the leader of interest as the sources of 222Rn identification. The second fast developing direction is including the radon in the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) models as a source of the boundary layer ionization. And here we encounter with second direction which is not connected with the earthquake forecast problems. It is the role of air ionization by radon as a source of the Global Electric Circuit (GEC) modification. In this publication we try to unite all these problems to present more complex view on radon as important element of our environment.

Well-exposed fluvio-deltaic Surma Group successions of the northeastern Surma Basin provide a direct analogue to its equivalent subsurface successions—the main hydrocarbon reservoirs in Bangladesh. The various levels of lateral and vertical heterogeneities lead to changes in the subsurface sandbody distributions and their reservoir qualities across fields. The present study aims to identify the heterogeneities in the fluvio-deltaic Surma Group successions by detailed sedimentological analysis of the Mio-Pliocene-aged Surma Group successions. Eleven lithofacies have been identified and based on the recurring pattern of these facies within the Surma Group, five facies associations have been identified. The facies associations indicate that these deposits originated from diverse depositional settings, ranging from the prodelta to the alluvial to delta plain. As a whole, they represent a progressive shallowing upward succession with two transgressive events in the upper and lower parts of the Surma Group. Five levels of heterogeneities such as prodelta mud (HL-I), tidal heterolithic deposits (HL-II & III), interbedded mud units in delta front (HL-IV), and intercalated mud laminae, drapes (HL-V) have been identified in the fluvio-deltaic Surma Group. The study suggests that the various levels of reservoir heterogeneities are controlled by the occurrence and distribution of various facies associations.

The main goal of the research is to assess the effectiveness of aerial ground-penetrating radar (GPR) in delineating boundaries between disturbed clayey ground and rock formations along the old Wadi. The study also tackles technical challenges associated with drone-based GPR applications in urban open areas, considering geological heterogeneity, topographic variations, and environmental conditions. The potential of using drone-based GPR equipped with an unshielded dipole antenna as a transformative tool in geospatial analysis was demonstrated. The research encompasses academic insights and field borehole drilling experiments to validate the accuracy and reliability of drone-based GPR data in real-world scenarios. Anticipated contributions include advancements in understanding drone-based GPR technology for mapping disturbed soil boundaries in foundation engineering applications and other relevant areas and recommendations for optimizing its performance in challenging terrains.

Abstract: Earth’s Phanerozoic history is marked by about 100 major events (“golden spikes”) indicating the beginnings of new stages in the Geologic Time Scale (GTS). Stage boundaries signify major deterministic events that are relatively well-known and can be correlated worldwide, mostly on the basis of appearances or disappearances of fossil species. The latest stage (Anthropocene) is human-made and currently involves rapid (approximately linear) increase in average atmospheric temperature. The main cause of Earth’s past and current temperature increases is addition of CO2 to the atmosphere. There is continuous exchange of greenhouse gases between atmosphere and oceans, which (per unit of volume) contain 50 to 60 times as much CO2 as the atmosphere. Abundance of new very precise worldwide observations is allowing multifractal modeling of weather and climate during the latest Quaternary stages during which some variables like temperature of the atmosphere and frequency of forest fires can be described by using the Pareto-lognormal frequency distribution model with parameters that are subject to continuous (deterministic) changes, similar to those in the earlier Phanerozoic.

Numerous experimental and theoretical studies undertaken to determine the effective-stress coefficient for seismic velocities in rocks stem from the importance of this geomechanical parameter both for monitoring changes in rock saturation and pore pressure distribution in connection with reservoir production, and for overpressure prediction in reservoirs and formations from seismic data. The present work pursues a task to determine, in the framework of a low-frequency laboratory study, the dependence of the elastic moduli of n-decane saturated sandstone on the relationship between pore and confining pressures. The study was conducted on a sandstone sample with high quartz and notable clay content in the quasi-static regime, when the 100 ml tank filled with n-decane was directly connected to the pore space of the sample. The measurements were carried out at a seismic frequency of 2 Hz and strains, controlled by semiconductor strain gauges, not exceeding 10-6. The study was performed using a forced-oscillation laboratory apparatus utilizing stress-strain relationship. The dynamic elastic moduli have been measured in two sets of experiments: at constant pore pressures of 0, 1 and 5 MPa and differential pressure (defined as a difference between confining and pore pressures) varied from 3 to 19 MPa; and at а constant confining pressure of 20 MPa and pore pressure varied from 1 to 17 MP. It was shown that the elastic moduli obtained in the measurements are in good agreement with the Gassmann moduli calculated for the range of differential pressures used in our experiments, which corresponds to the effective-stress coefficient equal to unity.

Translational rockslides caused by toe excavation are one of the commonly seen geohazards in mountainous regions due to line traffic construction. Quantifying their failure extension length (FEL) and travel distance are of significant interest as well as huge challenges in landslide hazard assessment. In this paper, a simple criterion is proposed for predicting these two factors based on the principles of rigid body limit equilibrium and kinetic energy theorem. Further, the proposed criterion is validated against the field observations and numerical results with a practical case of the Xinjianan landslide, a medium-sized translational rockslide that occurred in 2013 in Nanchuan District, Chongqing. The findings indicate that a tiny discrepancy can be found in FEL between the field observations and the proposed criteria, while this discrepancy could be considerable in travel distance between different methods. However, the relative discrepancies all fall within 20%, deemed acceptable.

The Gulf of Suez is a sedimentary rift-basin with thick evaporites and complex structural settings that presents challenges to achieving quality imaging based on seismic reflection data. Such subsurface settings and the lack of 3D seismic data poses challenges in building reservoir models to guide hydrocarbon exploration and development, especially in southern parts the Gulf of Suez. This study attempts to overcome the challenges of both 2D seismic interpretation and building a static reservoir model of two Miocene reservoir formations, Ras Ghara Field, Southern Gulf of Suez. Integration of well log data analysis of nine wells and regional-tectonics cognizant 2D seismic interpretation of nineteen 2D seismic lines enabled a reliable 3D static reservoir model. The Miocene reservoir characterization reveals heterogeneity in lithofacies distribution, with thick evaporite sequences causing seismic reflection distortions. Significant hydrocarbon accumulations are associated with horst fault blocks. The Miocene reservoirs demonstrate effective porosity of 14-33%, water saturation of 16-65%, and hydrocarbon saturation of 35-84%. Two new prospects are identified based on structure, facies distribution, and petrophysical results, with estimated original oil in place (OOIP) of 2,002,404 stock tank barrels (STB) and 1,175,580 STB for the Kareem Formation, and 6,302,916 STB and 2,841,894 STB for the Rudeis Formation.

Voluntary carbon markets open horizons for private companies, public institutions, and individuals developing CO2 removal projects in peatlands to reduce overall carbon footprint. These steps, however, should be in line with appropriate rewetting targets. To do so, the base line information about the status of area must be assessed. Here we follow the methodology set by the carbon certification program standards which define the necessary steps for reference conditions assessment. In this study we practically test the fulfilment of necessary drained peatland baseline evaluation for paludiculture and carbon certification activities. Estimates on the greenhouse gas emission mitigation potential were summarized to define priorities and propose quantifiable measures with suggested paludiculture implementation. Our hands-on approach shows that it is possible for companies to implement rewetting strategies without large EU level project fundings when the conservation and economic aspects are met, thus boosting climate mitigation actions.

The Lac Fallère area, located in the upper Clusellaz Valley (tributary of the middle Aosta Valley), is shaped in micaschist and gneiss referred to as the Mont Fort Unit (Middle Penninic) near the tectonic contact with calcschist and marble of the Aouilletta Unit (Combin Zone). Lac Fallère exibits an elongated shape and is hosted in a WSW–ENE trending depression, according to slope direction. This lake also shows a semi-submerged WSW–ENE rocky ridge, that longitudinally divides the lake. This evidence, in addition to the extremely fractured rocks, indicates the presence of a wide deep-seated gravitational slope deformation (DSGSD), even if this area is not yet included within the regional landslide inventory of the Aosta Valley Region. The Lac Fallère area also shows reliefs involved by glacial erosion (roche moutonnée), a wide cover of subglacial sediments and many moraines essentially referred to as Lateglacial. The DSGSD evolution in a glacial environment produced, as observed in other areas, significant effects both on the facies of Quaternary sediments and the formation of numerous and wide moraines. Glacial slope sectors and lateral moraines displaced by minor scarps and counterscarps and glaciers using trenches forming several arched moraines suggest an interplay between glacial and gravitational processes, which shared part of their evolution history.

Earthquakes often entail alterations in the groundwater flow regime, in the phreatic level, in the surge and loss of springs, and discharge in brooks and springs. As for the latter, there are diverse feasible mechanisms, so there is a variety of theoretical approaches to elucidate the post-earthquake effects on spring discharge. This study adopts a conceptual approach, primarily presenting various methods for estimating water released by earthquakes by calculating discharge surpluses in springs. This study delves into refined techniques rooted in groundwater hydraulics, displaying applications of analytical and simulation methodologies to quantify seismic-induced groundwater discharge in springs. Investigating springs as natural indicators, this research applies mathematical precipitation-runoff models, specifically the CREC model, for simulating hydrographs in post-earthquake scenarios. We apply analytical procedures or mathematical simulation techniques employed in groundwater hydraulics for natural aquifer recharge calculations. Firstly, we briefly describe the methods based on the analysis of depletion curves of hydrographs in spring discharge. Additionally, certain mathematical rainfall-runoff models utilized in simulating hydrographs of karstic springs, along with derived analytical approximations, are adapted for this scenario. These hydraulic calculations involve the depletion coefficient and hydrodynamic volumes of aquifers, parameters that reveal certain aspects of the relation between groundwater and earthquakes, such as: a) Acknowledging faults as the primary geological structures in transmitting pore pressures due to earthquakes. b) Clarifying mechanisms of permeability enhancement during earthquakes. c) Highlighting the substantial volumes in motion within the Earth's interior, which, for instance, could constitute a significant source for the origin of mineral deposits. The findings highlight diverse features of the complex interplay between hydrogeological and tectonic processes at challenging spatial and temporal scales.

A Global Historical Megatsunami Catalog (GHMCat) is presented, including all the events documented to date. The oldest event is from 1674 the Santorini event, ~3600 BP, is considered a particular case, and the most recent is from 2020. Although there is no formal definition of the term megatsunami, it is usually associated to the maximum height or runup reached by the waves when hitting the coast. Here we adopt this criterion and define a megatsunami when the wave reaches a height of 40 m or more. The catalog is initially based on tsunamis that meet this condition from the two existing Global Historical Tsunami Databases (GHTDs). Then, each of the events have been researched for its available references including earliest in order to verify the data, mainly on the maximum height of the waves and their causes. After a comprehensive review of published or on-line literature on historical tsunamis, including regional and local catalogs, papers, books, studies, reports and other available references, some of the oldest events have been ruled out. On the other hand, the literature review has also led to include in the catalog some new megatsunamis not included in the GHTDs. The global catalog contains 37 events from 1674 to 2023, 11 of which are either not included in the GHTDs or included with (erroneous) wave heights <40 m. All of them are described and characterized as far as can be gathered within the limitations of the scope of a global catalog and the available information on each historical case. The analysis of recorded historical megatsunamis shows that virtually all have been originated by large landslides, most of them subaerial landslides, with an unusual frequency of events in certain glaciated areas of Alaska and Norway, mainly in the first decades of the 20th century. In contrast, it is worth noting the small number of megatsunamis caused by submarine landslides or associated with volcanic explosions.

To study the oil and gas distribution law in the slope area of the hydrocarbon-bearing basin, the distribution of the effective dominant path of sand bodies transporting oil and gas is a fundamental subject. In this paper, we introduced an effective way to predict such transport paths and applied the method to study the sand bodies of the lower first member of the Shahejie Formation (E3s1L) in the southern area of the Qikou Depression in the Bohai Bay Basin (Qinan). We first define the dominant hydrocarbon transport paths of the sand body as the superimposed area of the contiguously distributed area of the sand body and the paleotectonic ridges at the top of the formation where the sand body is located. Then we define the oil and gas supply area of the oil source fault as the overlapped area of the favorable transport path of oil source faults and the width of the associated fracture zone. Stacking the two areas, we can delineate the distribution of the effective dominant transport paths of the sand bodies. The results from our case study of the Qinan area show that the effective dominant paths of sand bodies transporting oil and gas in the E3s1L Formation in the Qinan area are distributed throughout the region with better development in the west than in the east. The paths in the eastern part extend from north to south, while the paths in the western part extend from northeast to southwest, which is favorable for the convergence of oil and gas generated from the source rocks of the third member of Shahejie Formation (Es3) in the Northeastern main depression of Qikou and the formation of oil and gas accumulations in the l E3s1L Formation. The results are consistent with the fact that oil and gas discovered in the E3s1L Formation in the Qinan region are mainly distributed in the central and western regions, with a small amount in the southeast. Therefore, the suggested method is feasible to predict the effective dominant path of sand bodies transporting oil and gas.

The Sidi Bou Othmane (SBO) pegmatite district is in the Central Jebilet massif, Occidental Meseta domain, Morocco. The SBO district is hosted essentially in a volcano-sedimentary series composed of the Late-Devonian Sarhlef shales. Pegmatite bodies crop out as dykes, which are oriented from N-S to E-W and are generally variably deformed with ductile and/or brittle structures with ante, syn- or post kinematic criteria. Petrographic observations of pegmatite dykes show that feldspars (i.e., albite, microcline) are the most abundant mineral phases, followed by quartz, micas, with tourmaline and accessory mineral as garnet, and zircon and secondary minerals as clinochlore, sericite, and illite. The geochemical study indicates that, according to the major pegmatite elements dykes of SBO have mainly a S-type granitic composition and are peraluminous granite with calc-alkalic affinity. The study of trace elements pointed out that SBO pegmatite formed under post-orogenic syn-collision context during the Variscan orogeny, by the partial melting of argilliferous sediment and can be ascribed to the muscovite-bearing pegmatite; moreover, they have a good potential on ceramic and contain minerals, such as feldspar, recently assessed as Critical raw materials by European Union.

Sediment cores were collected from a mudflow lobe (80 m water depth) offshore of the Mississippi River’s Southwest Pass in 2017 to better understand the sedimentology near the lobe entraining the SS Virginia shipwreck (sunk by a German U-Boat in 1942) and surrounding Mississippi River Delta Front. Core analyses included 210Pb/137Cs geochronology, granulometry, and X-radiography. Sediment accumulation rates (SAR) calculated from excess 210Pb activity in multicores are 0.22–0.29 cm/yr at seabed depths less than 20 cm and 0.29–0.51 cm/year at depths greater than 20 cm. 137Cs accumulation rates are ~0.15 to ~ 0.37 cm/year since 1954 and 1963 respectively. Sediment accumulation rates from 210Pb and 137Cs geochronology and indicators of relative sedimentation and bioturbation from X-radiographs suggest that rates of sediment accumulation near Virginia have declined since ca. the mid-20th century. This may be explained by the multi-decade downslope mass transport of the mudflow lobe in which the shipwreck is embedded, and decreases in sediment supply delivered offshore from the Mississippi River. Mass-transport calculations of the Virginia lobe derived from core properties and published lobe advection rates suggest downslope mass transport is far higher than sediment resupply from the Mississippi River, consistent with recent studies of delta retreat.

Based on the theory that certain geological fault zones surpass the influence of other abiotic and biotic site factors for mound-building red wood ants (Formica rufa group), an area northwest of Tuttlingen with extremely high ant density was systematically surveyed, taking into account official geological maps. Indicators of active fault zones were recorded. These were additionally examined in a model area using a high-resolution digital terrain model (DTM 1-meter). The study area is located in the Swabian Jura at the western edge of the Albstadt shear zone, where recent seismic activity is documented by frequent earthquakes. The earthquakes are considered to be the cause of gas permeability in the strike-slip faults, around which the clustered nesting sites were mapped. More than 2900 nests of the Formica rufa group were mapped and their site conditions analysed in different geomorphological positions. In addition to recording forestry site conditions, a geological classification was carried out, as well as an analysis of soil gases/mofettes in areas with the highest nest densities. The correlation between forestry data and nest sites did not reveal a consistent pattern. In contrast, the relationships with gas-permeable strike-slip faults, though predominantly covered with Quaternary sediments, could be identified through CO₂- and CH₄-rich gas emissions, mofettes, and doline features. The interpretations were complicated due to different mapping ages and signatures of the geological maps. Fortunately, this issue could be analysed and evaluated in a sub-area with mappings from 1936/71 and 1995 regarding presumed or confirmed fault trajectories. New insights were gained regarding the relationships between Formica nest sites, gas-permeable faults, and groundwater levels on narrow mountain ridges with steep relief. Here, the lateral emission of CO₂ depends on the groundwater level and does not occur vertically directly above the fault. C- and H-isotope data of CO₂ and CH₄ do not show a consistent source. While methane is mainly derived from biological production by Archaea in the upper crust and in swamp areas, CO₂ mostly originates from the mantle. The results of the investigations revealed that a multitude of geological and ecological factors must be considered for the interpretation of the occurrences of Formica sites northwest of Tuttlingen. These factors confer a distinct originality to this region, contributing not only to an enhanced understanding of site conditions but also furnishing information regarding environmentally relevant gas emissions.

The past is an important lesson for facing the future with greater awareness. In this context, school plays a key role to spread knowledge on natural phenomena and to promote behavior change. Together with researchers, teachers can be strong allies to build more resilient future citizens. The Istituto Nazionale di Geofisica e Vulcanologia (INGV) school training activities provide tools to get ready for the next earthquake and/or tsunamis. About 5000 students from both Middle schools (ISCED 2) and High schools (ISCED 3) were involved in active learning activities based on a flipped-up approach during special online scientific events in the Pandemic. Online lab activities were conducted respectively during the European Researchers’ Night (“Earthquakes: history teaches us the future: researchers for a day with experimentation in didactics for ESL”) and during the World Water Day 2021 and the World Earth Day 2021 (“Tsunamis: history teaches us the future researchers for a day with experimentation in didactics for ESL”). These two Episodes of Situated Learning (ELS) experiences triggered students’ interest, favoring remote learning, developing life skills, focusing on historical seismic studies of both past earthquakes and tsunamis.

The ocean of the planet Earth, which is as thin as a cicada's wings, play a role of catalyst to the evolving of the ancient plates. The ocean and the ancient plates of the planet Earth formed a magic Whirlpool of the Cross “十” shape with double whirl arms. The central core of the Cross “十” shape Ocean-Plate Whirlpool is the Hint Black Hole (HBH) located at the Northwestern Pacific Ocean floor which is decorated with invisible geological age distribution of clam shell. The double whirl arms of ancient plates are distinguished to be left and right at Bering Strait and just formed the ring of fire of the Pacific Ocean.The Cross “十” shape Ocean-Plate Whirlpool of the planet Earth was born at the equatorial area about 66 Ma ago and is rotating anticlockwise and drifting northward by a mini footprint of spiral at s speed of as slow as a snail moves. Obviously, it would come to an end at the North Pole about 100 Ma in the future. So, the lifespan of the Cross “十” shape Ocean-Plate Whirlpool would be about 166 Ma, from birth to death on the planet Earth.What is hidden behind of the evolving of the oceans and the ancient plates is the golden law of geometrical progression of ladder pagoda of attributes evolving of them, which is “1➨2➨4➨2➨1”, including upgrading and downgrading process. It seems to hint that what is happening on the planet Earth, including the top-level species of human race, is actually unique in the Milky Way galaxy in essence.The property roles of expanding and shrinking are existing in the Cross “十” shape Ocean-Plate Whirlpool by face to face and the exchange of property roles took place tens of millions of years ago. And that situation would happen again tens of millions of years in the future. It seems a giant circle of evolving of the exchange of property roles would be drawn by the Cross “十” shape Ocean-Plate Whirlpool in her lifetime.

This paper presents the results of experiments on the study of a hydraulic fracture interaction with a preliminary created fracture that simulates a natural fracture. A distinctive feature of the conducted experiments is the ability to use an ultrasonic sounding to measure the fracture propagation and opening simultaneously with the fluid pore pressure measurements at the several points of the porous saturated sample. It allows to obtain the pressure distributions at various experiment stages and to establish a relation between the pore pressure distribution and the hydraulic fracture propagation and its interaction with macroscopic natural fracture. The possibilities of active ultrasonic monitoring have been expanded due to preliminary calibration experiments, which make it possible to estimate the fracture opening value by attenuation of ultrasonic pulses. The experiments show that the fracture propagation is limited by the natural fracture. This is caused by the hydraulic fracturing fluid leaks into the natural fracture, so both hydraulic fracture and natural fracture compose united hydraulic system.

A high-resolution catalog for the 2014 Ms6.5 Ludian aftershocks was constructed based on the deep learning phase-picking model (CERP) and the seismic-phase association technology (PALM). A specific training strategy, which combines the advantages of the conventional short-long window average energy ratio algorithm (STA/LTA) and AI algorithm, is employed to retrain the CERP model. The P & S phases were accurately detected and picked on continuous seismic waveforms by the retained AI model. Hypoinverse and HypoDD were utilized for precise location of 3286 events. Compared to the previous results, our new catalog exhibits superior performances in terms of location accuracy and the number of aftershock events, thereby enabling a more de-tailed depiction of the deep-seated tectonic features. According to the distribution of aftershocks, it can be inferred that, (1) the seismogenic fault of Ludian earthquake is the NW-trending Baogunao-Xiaohe Fault. (2) the Ludian aftershocks interconnected the discontinuous NW-trending Baogunao-Xiaohe Fault, and intersected with the Zhaotong-Ludian Fault. (3) It suggests that the NE-trending Zhaotong-Ludian Fault may have been intersected by the NW-trending Baogunao -Xiaohe Fault, indicating that the Baogunao-Xiaohe Fault is likely a rela-tively young Neogene fault.

Correlations between the Neoproterozoic belts surrounding the West African Craton and northern Brazilian cratons have long been a subject of interest and controversies. Due to the splitting of African and South American continents by the Atlantic oceanic domains, no direct links are preserved, requiring relying on various geological or geophysical characteristics to propose such correlations. In addition to the opening of the Atlantic oceanic domains, another difficulty arises from the covering of northern Brazilian belts by upper Paleozoic and Mesozoic basins, making these correlations speculative. Recently, four orogens have been evidenced in the Neoproterozoic belts of the western margin of the West Africa Craton, while the belts on the eastern side underwent only one orogeny. Similarities with the Pan-African I (900-650 Ma) and with the Pan-African II (650-480 Ma) orogenic events have been evidence in the western Brazilian belts (Araguay and Paraguay). The two first orogens on the western margin (Pan-African I and Pan-African II) can thus be extended to the western Brazilian belts and can be considered as parts of a single geodynamic system running from the Mauritania to the Paraguay including the “Gurupi rift” as an aulacogen connected to the NNW-SSE Panafrican I and II oceanic domains. Consequently, the eastern Brazilian belt should rather be linked the Eastern Trans-Saharan belts.

Increasing groundwater levels (GWL) may become one of the most serious issues for Odessa City, Ukraine. This study investigated the spatial distribution characteristics and multifractal scaling behaviour of groundwater level/depth fluctuation for a Quaternary aquifer in Odessa City using a geostatistical approach and a Multifractal Detrended Fluctuation Analysis (MF-DFA). These two methods were applied to monthly GWL fluctuation time series from 1970 to 2020 to monitor 72 hydrogeological wells situated in different parts of Odessa City. The spatial distribution of GWL revealed an overall trend of decline and recovery from 1970 to 2020 in the study area, except for most of the southern region, where a persistent recovery of groundwater depth was observed. The MF-DFA results suggest that the dynamics of GWL fluctuations have multifractal characteristics in the Odessa City area. In addition, both long-range correlation and fat-tail probability distribution contribute to multifractality. However, long-range correlations among fluctuations made a major contribution to the observed multifractality of the GWL fluctuations time series. The generalized Hurst exponent shows a wide range of change (0.20 < h(q) < 2.85), indicating the sensitivity of GWL fluctuations to changes in small scale factors and large-scale factors. Regarding the long range correlations of GWL depth, the Hurst exponents (q = 2) demonstrated the positive persistence of groundwater depth recovery in the southern region and the persistence of groundwater depth variation in the other regions of the study area. The dynamic changes in GWL depth in the Odessa City area may be affected by both natural structural and anthropogenic factors.

Permeability is an important transport property of porous materials, and its accurate evaluation is relevant in studying applied tasks such as CO2 injection into reservoirs and investigating groundwater quality issues. This study examines the dependence of permeability on total and connected porosity, hydraulic tortuosity, specific surface area, and mean pore radius based on the data of 408 cubic sub-volumes extracted from heterogeneous and naturally fractured cylindrical carbonate samples before and after injection of HCl solutions into them. These parameters were computed using pore-scale modeling of fluid flow. It has been shown that permeability correlates well with porosity and mean pore radius with correlation coefficients of R2≈0.65-0.79 for heterogeneous samples. It was found that the presence of natural fractures significantly influenced the relationship between permeability and other parameters. The relationship between permeability k, tortuosity τ, and specific surface area S is described by power laws k~τ-α and k~S-β, with coefficients α and β substantially exceeding those in the Kozeny-Carman equation. It was also found that there is a parabolic relationship between connected and total porosities both before and after rock dissolution with R2≈0.96-0.99. This allowed for estimation of percolation threshold porosity that are in good agreement with the literature data.

Landward displacement of coastal boulders is a morphodynamic signature of inundation and erosion processes. The study of these movements can provide useful elements for coastal hazards assessment. In this note, using open-access remote sensing resources, the 2017 to 2021 annual mobility of medium, coarse, and very coarse boulders spread over about 100 km of the eastern coast of the Gulf of Taranto (Italy) is detected. A large interannual variability was found. As established by multi-temporal object-based image analysis, one hundred and ten boulders changed position between July 2019 and June 2020, no boulder was displaced between July 2017 and July 2018, and few boulders were displaced between July 2018 and July 2019, and June 2020 and September 2021. A hydrodynamic analysis of the investigated boulder displacement led to the estimation of the energy of nearshore waves that hit the coast. The storm Detlef, that crossed over the Mediterranean Sea during 11-13 November 2019, is presumed to be the main cause of the massive displacement phenomenon.

The Kuznetsk Basin (Kuzbass) is one of the largest coal basins in Siberia and a reference area for the ancient Angaraland continent. The proximity of the Kuzbass and Siberian platform caused their biotic similarities in the Late Palaeozoic. However, due to biota endemism, the Kuzbass Upper Palaeozoic does not correlate directly with the International Chronostratigraphic Chart (ICC). This paper discusses radioisotopic (CA-ID-TIMS) dating of zircons from a volcanic tuff located in the Starokuznetsk Formation. This level matches the interval of the Balakhonka/Kolchugino (B/K) floral change in Kuzbass, i.e. the gradual replacement of cordaitoid-dominated wet forests (Balakhonka flora) by more arid fern-pteridosperm-cordaitoid assemblages (Kolchugino flora). New age (276.9 ± 0.4 Ma) directly correlates the Starokuznetsk Fm with the Upper Kungurian of the ICS. We compared the Kuzbass data with data of the Western Verkhoyanie, where Middle Permian ammonoids (Sverdrupites assemblage) occur in strata recording the B/K floral change. Available (ICC) and new datings indicate the lag between the B/K floral change in low (Kuzbass) and high (Verkhoyanie) latitudes of Angaraland. The B/K floral change in the Kuzbass began in the early Late Kungurian and was completed by the end of this age. In contrast, the B/K floral change in Verkhoyanie began at the end of the Late Kungurian and was completed in the Late Wordian. The delay in the floral changes at different latitudes of Angaraland suggests that existing interregional correlations need further improvement.

Seismic velocity inversion is one of the most critical issues in the field of seismic exploration and has long been the focus of numerous experts and scholars. In recent years, the advancement of machine learning technologies has infused new vitality into the research of seismic velocity inversion and yielded a wealth of research outcomes. Typically, seismic velocity inversion based on machine learning lacks control over physical processes and interpretability. Starting from wave theory and the physical processes of seismic data acquisition, this paper proposes a method for seismic velocity model inversion based on Physical Embedding Recurrent Neural Networks. Firstly, the wave equation is a mathematical representation of the physical process of acoustic waves propagating through a medium, and the finite difference method is an effective approach to solving the wave equation. With this in mind, we introduce the architecture of recurrent neural networks to describe the finite difference solution of the wave equation, realizing the embedding of physical processes into machine learning. Secondly, in seismic data acquisition, the propagation of acoustic waves from multiple sources through the medium represents a high-dimensional causal time series (wavefield snapshots), where the influential variable is the velocity model, and the received signals are the observations of the wavefield. This forms a forward modeling process as the forward simulation of the wavefield equation, and the use of error back-propagation between observations and calculations as the velocity inversion process. Through time-lapse inversion and incorporating the causal information of wavefield propagation, the non-uniqueness issue in velocity inversion is mitigated. Through mathematical derivations and theoretical model analyses, the effectiveness and rationality of the method are demonstrated. In conjunction with simulation results for complex models, the method proposed in this paper can achieve velocity inversion in complex geological structures.

Abstract: We present the Hellenic orogenic belt`s main geological structure and architecture of deformation, trying to reach a better interpretation of its geotectonic evolution during the Alpine orogeny. We based on our recent studies and experience about the deformational history of the Hellenides but also on the more modern views, published by other colleagues, concerning the Alpine geotectonic reconstruction of the Hellenides. The structural evolution started during the Permo-Triassic with the continental breaking of the Pangea super-continent and the birth of the Neotethyan ocean realm. Bimodal magmatism and A-type granitoid intrusions accompanied the initial stages of the continental rifting, followed by Triassic-Jurassic multiphase shallow- and deep-water sediment deposition on both formed continental margins. These margins were the Apulian margin, containing the Pelagonia, at the western part of the Neotethyan ocean and the European margin, containing continental parts of the Serbo-Macedonian and Rhodope massifs (upper Vertiskos-Kimi and middle Sidironero-Kerdylia units), at the eastern part of the Neotethyan ocean. Deformation and metamorphism are recorded in six main deformational stages from the Early-Middle Jurassic to the present day, getting started with an Early-Middle Jurassic Neotethyan intra-oceanic subduction and ensimatic island arc magmatism. Compression, nappe stacking, calc-alkaline magmatism and high-pressure metamorphic events related to subduction processes alternated successively through time with extension, orogenic collapse, medium- to high-temperature metamorphism, adakitic and calc-alkaline magmatism and partly migmatization, related to uplift and exhumation of deep crustal levels as tectonic windows or metamorphic core complexes. An S- to SW-ward migration of the dynamic peer compression vs. extension is recognized during the Tertiary Alpine orogenic stages in the Hellenides. In any case extension and crustal uplift follow compression and nappe stacking or compression and exten-sion act simultaneously at different parts of the orogen. We suggest that all ophiolite belts in the Hellenides originated from a single source and this was the Neotethyan Meliata/Maliac-Axios/Vardar ocean basin, parts of which obducted during the Mid-Late Jurassic on both continental margins, Apulian`s (containing Pelagonia) and European`s (containing units of Serbo-Macedonian/Rhodope nappes stack), W-SW-ward and E-NE-ward, respectively. In this case, the ophiolite nappes should be considered as far-traveled nappes on the Hellenides continental parts associated with the deposition of Middle-Late Jurassic ophiolitic mélanges in basins at the front of the adjacent ophiolite thrust sheets. The upper limit of the ophiolite emplacement is the Mid-Upper Jurassic, Callovian-Oxfordian, as shown by the deposition of the Kimmeridgian-Tithonian Upper Jurassic sedimentary carbonate series on the top of the obducted ophiolite nappes. The suture zone between the Pelagonian nappe and the External Hellenides as part of the Apulia plate has been dated about coeval in time with the Paleocene-Eocene to Oligocene SW-ward over thrusting (between 56 Ma and 32 Ma), along the Nestos thrust zone, of the Serbo-Macedonian/Rhodope Sidironero-Kerdylia unit on the tectonically lower-most Rhodope Pangaion unit. In this sense, the Nestos shear zone in the Rhodope province continues until the suturing between Pelagonian and External Hellenides and is the same with this. Therefore, the lower-most Rhodope Pangaion unit could be regarded, as the continuation of the marginal part of the Apulian plate (External Hellenides) that was underthrusted during the Paleocene-Eocene below the unified Sidironero-kerdylia unit and the Pelagonian nappe, following the Paleo-cene-Eocene subduction and closure of a small ocean basin at the western of the Pelagonian (Pindos-Cyclades ocean basin). Was preceded by the Late Cretaceous subduction of the Axios/Vardar ocean remnants below the European continental margin and the final closure of the Axios/Vardar ocean during the Paleocene-Eocene, which were associated with the overthrusting of the European origins Vertiskos-Kimi nappe on the Sidironero-Kerdylia nappe and subsequently the final collision of the European margin and Pelagonian fragment. In this scenario should be assumed, that the true Vardar/Axios suture zone should originate in the deep further to the East between the Rhodope units in the northern Rhodope province, while the Axios/Vardar zone in its today outcrops position should be considered as allochthonous. Subsequently, during synorogenic Oligocene-Miocene extension, the lower-most Rhodope Pangaion unit exhumed as a metamorphic core complex below the Serbo-Macedonian/Rhodope nappes, which were successively exhumed, from the higher to the middle, during Paleocene-Eocene and Eocene-Oligocene extension, respectively. A retreating subduction zone associated with rollback of the subducted lithosphere under the Pelagonian and the other Internal Hellenides nappes stack and/or slap break off and mantel delamination could well explain the Tertiary synorogenic extensional tectonics in the Internal Hellenides, taking place simultaneously with compression in the External Hellenides and the Hellenic foreland or accretionary wedge.

This paper provides a data set collected with an underground shaft mobile mapping system equipped with rotated Velodyne VLP-16, Velodyne Ultra Puck VLP-32c, Livox Tele-15, and IMU Xsens MTi-30. The ground truth data was acquired with a geodetic survey including 15 ground control points and 6 Faro Focus 3D Terrestrial Laser Scanner stations of a total 273,784,932 of 3D measurement points. The aim of the paper is to describe the data set of a mining shaft. Obtained mobile mapping data can be the basis for detailed analyzes of the technical condition of the shaft. This data set provides an end-user case study of realistic applications in mobile mapping technology. It provides ground truth data, mobile mapping data and software tools for manipulation and visualization. The project is released and maintained at https://michalpelka.github.io/mine-mapping-dataset/.

The overall seismicity of Pakistan from 24 B.C to 2020 is analysed in terms of its multifractal behaviour. Seismic events of magnitude ML=3.0 and above are spatially clustered into four distinct groups, each one corresponding to a different region of high seismic activity. The Multifractal Detrended Fluctuation Analysis (MFDA) method applied on each cluster reveals pronounced inter-cluster heterogeneity in terms of the resulting generalised Hurst exponent and fractality spectrum, possibly due to the particular tectonic characteristics of the regions under investigation. Additional results on the variability of the Gutenberg–Richter b‐value across the defined clusters further corroborate the uniqueness of the seismic profile of each region.

The paper presents the results of geological interpretation of the set of geophysical data with substantiation of geological-geophysical and structural geological criteria of oil and gas potential of the sedimentary cover deep levels and the upper part of the consolidated basement in the central West Siberian oil and gas basin. Gravity, magnetic and seismic data are used. The character of manifestation of the regional Pai-Khoi–Altai strike-slip zone covering an area up to 200 km wide is shown. The study provides a basis for assuming a spatial genetic relationship between the conditions for localizing hydrocarbon traps and the extensional and strike-slip and thrust structures of the zone that form so-called "flower structures". The substantiated structural geological criteria were applied to solve the problem of pattern recognition with training, highlighting the environments favourable for localization of hydrocarbon fields.

The search for a continuous continental record of interglacial periods in semi-arid regions is problematic due to the absence of stable and continuous sedimentary systems over time in this type of climate. In this work, a relatively stable basin is described and analyzed during the last interglacial in a semi-arid region of the western Mediterranean. For this purpose a geomorphological, stratigraphic and sedimentological study has been carried out, with dating through 230Th. A semi-endorrheic Pleistocene basin has been identified, with two units that correspond to a fluvial-lacustrine system (unit P1) and an alluvial system (unit P2). Unit P1 has been dated in the MIS 5e interglacial episode. A framework is described in which the Aljezares Quaternary basin can be considered as a possible source of paleoenvironmental and paleoclimatic information in semi-arid regions from the last interglacial period in future studies.

The development of mineral resources in China has gradually entered the stage of deep mining. Rock burst, collapse and other ground pressure disasters caused by deep high stress have become major hidden dangers restricting mine safety production. Microseismic monitoring technology is an important means of ground pressure risk prediction. In this paper, the wavelet coefficient threshold denoising method in time-frequency domain, STA/LTA method, AIC method, Skew and Kurtosis method are studied, and the automatic P-phase onset time picking model based on noise reduction and multiple detection indexes is established. Through the effect analysis of microseismic signals collected by microseismic monitoring system of coral Tungsten Mine in Guangxi, automatic P-phase onset time picking is realized. the reliability of the P-phase onset time picking method proposed in this paper based on noise reduction and multiple detection indexes is verified. The picking accuracy can still be guaranteed under the severe signal interference of background noise, power frequency interference and manual activity underground mine, Which is of great significance to the data processing and analysis of microseismic monitoring.

New detailed geological/structural mapping as well as field-based structural analysis were carried out to investigate the deformation pattern of well-preserved high-pressure rocks of the Blueschist Unit exposed on SE Syros (Cyclades, Greece). Our new geological mapping revealed the occurrence of a metasedimentary sequence underlain by a meta-igneous sequence. The contact between these two sequences displays typically interfingering patterns in map-scale due to folding. The earlier ductile deformation phase recognized in the mapped area, is associated with the development of a penetrative foliation, which was formed under eclogite/blueschist facies conditions at peak conditions. The subsequent main deformation phase occurred at blueschist facies conditions synchronous with the early stages of exhumation of the high-pressure rocks. This phase is mainly associated with the formation of map-scale WNW-trending folds and a pervasive axial planar foliation linked with ESE-directed shearing. The main deformation ceased within blueschist facies conditions and exhumation of the rocks to greenschist facies conditions took place under very weak deformation. Greenschist retrogression observed in the southwestern part of the mapped area seems to be controlled by fluids rather than by intense deformation. Our results indicate that the high-pressure rocks of the Blueschist Unit exposed on Syros Island represent a large-scale pod of low deformation under blueschist to greenschist facies conditions, likely occupying the core of an extrusion wedge.

Geoethics is a field of knowledge currently in full development. Researchers in geoethics are pri-marily concerned with the anthropogenic interaction with the Earth system. Due to its nature, ge-oethics holds particular importance in sustainable development due to its nature as it aims to promote ethical human behaviour that does not negatively impact the Earth system. In the present research, we implemented an intervention program addressing various issues related to the sus-tainability of the Earth system, such as the exploitation of geological resources, the management of geological risks, and the conservation and promotion of geopatrimony. The intervention program was applied to higher education students in the geosciences field. A sample of 90 students from various geosciences courses completed an initial questionnaire, revealing limited knowledge about geoethics. This study resorted to mixed-method research involving interviews with some students who volunteered (n=52). The results showed that after applying the intervention programme most students had developed a deeper understanding of the topics addressed and recognised the con-tributions this scientific area can make to sustainable development. Additional research in geoethics education is essential to foster the integration of geoethics into the curricula of higher education institutions.

Agricultural activities profoundly affect groundwater levels in Saudi Arabia. The Al Kharj and Wadi Sahba areas, central Saudi Arabia, produce groundwater mainly from the Jurassic Arab/Jubaila aquifer and the overlying Cretaceous Biyadh aquifer. A geographic information system analysis and linear regression analyses show that from 1978 to 2016, significant agricultural overexploitation of the area’s fossil groundwater was associated with groundwater level drops of ~104 m in Al Kharj and ~48 m in Wadi Sahba. The estimated total volume of groundwater withdrawal across both aquifers during this period of wheat and alfalfa irrigation was ~27 billion m3. Landsat image analysis shows that agricultural land use decreased by ~66% over the same period due to overexploitation. The 2016 phase-out of alfalfa farming was associated with increases in groundwater levels: ~26 m in the Arab/Jubaila aquifer and ~3 m in the Biyadh aquifer. This recovery was likely due to lateral flow from surrounding fossil water within the aquifer given the absence of significant surface recharge. Model-based evaluation of surface recharge of the aquifer reveals very low potential for surface infiltration, with almost no recharge into the deep aquifer. Assuming that recovery continues at the present rate, the average groundwater level for the Arab/Jubaila aquifer is expected to reach 98 m below ground level in 2025 and 74 m below ground level in 2030.

A recent research has identified an inverse amplitude link between obliquity damping and short eccentricity amplification during Mid-Late Pleistocene based on LR04 δ18O and equatorial Pacific Site 846 sea surface temperature records, which is associated with the Earth's long-term cooling. In the present study, new evidence of this anticorrelation is presented from Antarctic D-CO2-CH4 records, global benthic-planktic δ18O, and regional (Atlantic, Pacific, Mediterranean, and Indian) climate-related proxies. Based on a critical review of theoretical constrains (Earth's oblateness changes and ice-volume phase lag in the obliquity band <5.0 kyr), this widespread and symmetric (bipolar) obliquity response damping has been interpreted as an effect of the obliquity-oblateness feedback, which could be the latent physical mechanism at the origin of the Mid-Pleistocene Transition (MPT). Indeed, results and considerations of the present work suggest that fast and positive/negative net variation of the Earth's oblateness in the obliquity band was controlled by dominant glacio-eustatic water mass component and, assuming a rapid response of the ice-volume to surface temperature changes, the mean obliquity lag response is estimated to be <5.0 kyr over the past 800 kyr. These elements may explain the interglacial/glacial damping observed in the obliquity response. The consolidation of the Earth's long-term icy-state in the subtrend IV, culminated with the post-MPT obliquity damping, might have contributed to the strengthening of the short eccentricity response by mitigating the obliquity ‘ice-killing’ during obliquity maxima (interglacials), favouring the obliquity-cycle skipping and a feedback amplified ice-growth in the short eccentricity band (obliquity damping hypothesis). This could suggest a different impact of the climate friction than what is generally believed, which is presumably the latent physical mechanism that triggers transient ‘competitive’ interaction between obliquity and short eccentricity started early during the Piacenzian.

Landslide is one of the most common geological disasters in China, which is characterized by suddenness and uncertainty, and it is difficult to realize accurate identification, early warning and forecasting of landslide disaster by conventional means. With the development of high-resolution remote sensing satellites and InSAR surface deformation monitoring technology, the traditional means of landslide monitoring data sources are limited, and there is a lack of effective methods to excavate the characteristics of the spatial distribution of landslide hazards and their triggering factors and other issues. In this paper, the area extending 10 km outside the VII isobar of the Gengma earthquake is taken as the study area, and 13 evaluation factors are screened out by integrating the factors of InSAR surface deformation, topography and geological environment, and the Bayesian Optimized Convolutional Neural Network (BO-CNN) is used for the evaluation of landslide susceptibility, and the BO-RF and PSO-SVM models are selected for the comparative analysis. The model accuracy evaluation was carried out by three indexes: ROC curve, AUC value and FR value, in which the ROC curves of PSO-SVM, BO-RF and BO-CNN were all close to the upper-left corner of the corner, and the AUC values were 0.9388, 0.9529, and 0.9535, respectively, and the FR value of landslide in the high susceptibility area of BO-CNN was as high as 14.9, and was higher than that of PSO-SVM and BO-RF, respectively. SVM and BO-RF model is 4.55 and 3.69 higher, the experimental results show that the BO-CNN model used in this paper has a better effect in landslide susceptibility evaluation, and the research results of the local government's disaster prevention and mitigation measures are of great significance.

Based on the plate tectonic theory, using seismic, geology, and public literature data, the proto-type basins and lithofacies paleogeography during the main geological periods of the passive continental margin of East Africa are restored, and the architecture characteristics of the basins and the differences in sedimentary fillings are analyzed. Combined with the dissection of discovered oil and gas reservoirs, four types of hydrocarbon accumulation models are established, the favorable hydrocarbon accumulation combinations (plays), and the next exploration directions in this region are explored. The passive continental margin basins in East Africa experienced three proto-type stages: intracontinental aborted rift during the Late Carboniferous-Triassic Karoo period, intracontinental-intercontinental rift during the Jurassic period, and passive continental margin basin since the Cretaceous period. Affected by the difference in sedimentary filling thickness between the rift period and depression period, four types of passive continental margin basins are formed: rifted type, depressed type, faulted depression type, and reformed delta type. The maximum sedimentary thickness of the rifted basins during the rift period is greater than 3500 m, and that during the depression period is less than 4000 m, forming a "single-source structure type" hydrocarbon accumulation model. The exploration of rifted basins should focus on giant structural traps at the top of the rifted sequences. The maximum sedimentary thickness of the depressed basins during the depression period exceeds 4000 m, forming a "single-source fan group type" hydrocarbon accumulation model. The dominant exploration directions of depressed basins are sedimentary sand bodies and structural-lithologic traps during the depression and drift periods. The maximum sedimentary thickness of the faulted depression basins during the rift period is greater than 3500 m, and that during the depression period is greater than 4000 m, forming a "double-source and double-combination type" hydrocarbon accumulation model. The exploration of rifted basins should focus on giant slip-collapse-debris flow sedimentary sand bodies on the upper and middle slopes. In addition to the large sedimentary thickness of the early rift and the central depression, the delta thickness of the reformed delta basins during the late depression period since the Miocene is more than 4000 m, and four structural belts (growth fault, mud diapir, thrust fold, and foredeep gentle slope) were formed from onshore to offshore, forming a "three-source and multi-combination type" hydrocarbon accumulation model. All four structural belts can form giant oil and gas fields.

A Mw7.8 earthquake struck Turkey on February 6, 2023, causing severe casualties and economic losses. This paper investigates the characteristics of strong ground motion and seismogenic fault of the earthquake. We collected and processed the strong ground motion records of 379 stations using Matlab, SeismoSignal and Surfer software, and obtained the peak ground acceleration (PGA) contour map. We analyzed the near-fault effect, the fault locking segment effect and the trampoline effect of the earthquake based on the spatial distribution of PGA, the fault geometry and slip distribution. We found that the earthquake generated a very strong ground motion concentration effect in the near-fault area, with the maximum PGA exceeding 2000cm/s2. However, the presence of fault locking segments influenced the spatial distribution of ground motion, resulting in four significant PGA high-value concentration areas at a local dislocation, a turning point and the end of the East Anatolian Fault. We also revealed for the first time the typical manifestation of the trampoline effect in this earthquake, which was characterized by a large vertical acceleration with the positive direction significantly larger than the negative direction. This paper provides an important reference for understanding the seismogenic mechanism, damage mode, characteristics and strong earthquake law of Turkey earthquake.

Recent scholarship highlights the significant environmental impact of the expanding hydraulic fracturing industry, necessitating up-to-date, comprehensive, national-scale assessments. Addressing this gap, this literature review provides a thorough analysis of the contamination and induced seismicity associated with shale-gas hydraulic fracturing in the United States. Employing systematic reviews of literature and federal reports up until July 2023, this study reveals multiple areas of concern. Notably, water and soil contamination, seismic activity, and air pollution are identified as significant challenges linked to shale gas extraction. However, ongoing studies are actively exploring remediation strategies, including advancements in monitoring and treatment technologies, aimed at mitigating the impact of contamination. These findings underscore the need for stringent regulations, continuous research, and effective management practices. This work emphasizes the importance of a collaborative approach that leverages field studies, experimental investigations, and computational advancements to ensure the responsible development of shale gas resources.

The study was performed to get the first insight of distribution of the low field volume magnetic susceptibility (MS) throughout the Zagreb city area, based on in-situ field measurements. Most interesting locations were selected for soil sampling and geochemical content was determined using ICP-OES. A geostatistical approach was applied on MS and geochemical results. Median value of 0.245 x 10-3 SI units is proposed to be used as relevant average MS value in Zagreb area. Mean concentrations of heavy metals in Zagreb soils (in µg/g) are: Pb (36.82), Zn (87.77), Cu (30.84), Cd (0.66), Cr (29.04), Co (11.89) and Ni (28.40), being relatively low in comparison with Croatian legislation and with European and world average. Boxplot analyses showed that 45% of studied elements do not show any anomaly, while the rest of elements mostly show only one weak anomaly and they are located at the same sites as MS anomalies. Statistical analysis found significant correlations between MS and following elements: Cd, Co, Fe, Mn, Na, Pb, Sb and Zn. In-situ MS measurements proved to be an efficient tool for initial screening of large areas in terms of elevated heavy elements concentrations, enabling cheap and fast assessment of state of environment.

Forests play a pivotal role within the ecological milieu, with alterations to forested areas not only impacting the local environment but also serving as a significant contributing factor to the occurrence of landslide disasters. However, the significance of forests in connection with deep seated landslides is diminished. In instances of extreme events, heavy rainfall, highly weak slopes, or seismic events, it is improbable for forest cover to exert any influence. Due to extensive agricultural irrigation, the Heifangtai terrace of China has experienced a significant rise in groundwater levels, leading to the occurrence of numerous landslides. Through field investigations and remote sensing interpretation, this study has revealed dramatical geomorphological changes in the Heifangtai terrace due to the occurrence of widespread landslides. These changes are likely to disrupt the pre-existing hydrological balance within the study area, thereby altering the conditions for landslides. Therefore, understanding how the evolution patterns of loess landslides can be assessed is a crucial scientific issue for predicting the conditions of landslide occurrence. Taking the evolution of landslides along the edge of the Heifangtai terraces as an example, this study first classified different types of loess landslides. Subsequently, the temporal and spatial distribution characteristics of these landslides were analyzed. Finally, utilizing the FLAC3D software, the stability of the loess slopes at the edge of the Heifangtai terrace was examined under the lateral evolution mode. The research findings indicate that the current loess landslides in the Heifangtai terrace exhibits longitudinal multi-periodicity and lateral group-occurring behavior, which is significant for assessing loess landslides along the edge of Heifangtai terraces.

The nine plateau lake watersheds in Yunnan are important ecological security barriers for southwest China, and the prevention and control of landslides are important considerations in the management of these watersheds. Taking Dianchi Lake watershed as a typical study area, this paper puts forward an improved comprehensive process of landslide susceptibility evalua-tion, and discusses the model’s sensitivity regarding susceptibility to landslides. The compre-hensive process is based on the weight of evidence method (WoE) and integrates many analyti-cal techniques, such as cross-validation, multi-quantile cumulative Student’s comprehensive weight analysis, independence tests, step-by-step modeling, ROC analysis and ROC-based sus-ceptibility zoning. This paper established fourteen models, with high accuracy and validity AUC, reaching 0.83-0.87 and 0.85-0.88, respectively. Moreover, according to the susceptibility zoning map compiled using the optimal model, 80% of landslides can be predicted in the very high and high susceptibility zones, which only account for 19.58% of the study area. Finally, the strategies for geological disaster prevention and ecological restoration deployment were put forward, providing great guiding significance for deeply understanding the landslide susceptibility of lake watersheds (Dianchi Lake) in Yunnan Plateau and guiding the planning and deployment of landslide prevention and mitigation activities and ecological restoration in the watershed.

The Earth’s synthetic gravitational and density models can be used to validate numerical procedures applied for a global (or large-scale regional) gravimetric forward and inverse modelling. Since the Earth’s lithospheric structure is better constrained by tomographic surveys than a deep mantle, most of existing 3-D density models describe only a lithospheric density structure, while 1-D density models are typically used to describe a deep mantle density structure below the lithosphere-asthenosphere boundary. The accuracy of currently available lithospheric density models is examined in this study. The error analysis is established to assess the accuracy of modelling the sub-lithospheric mantle geoid, while focusing on the largest errors (according to our estimates) that are attributed to lithospheric thickness and lithospheric mantle density uncertainties. Since a forward modelling of the sub-lithospheric mantle geoid comprises also numerical procedures of adding and subtracting gravitational contributions of similar density structures, the error propagation is derived for actual rather than random errors (that are described by the Gauss’ error propagation law). Possible systematic errors then either lessen or sum up after applying particular corrections to a geoidal geometry that are attributed to individual lithospheric density structures (such as sediments) or density interfaces (such as a Moho density contrast). The analysis indicates that errors in modelling of the sub-lithospheric mantle geoid attributed to lithospheric thickness and lithospheric mantle density uncertainties could reach several hundreds of meters, particularly at locations with the largest lithospheric thickness under cratonic formations. This numerical finding is important for a calibration and further development of synthetic density models of which mass equals the Earth’s total mass (excluding the atmosphere). Consequently, the (long-to-medium wavelength) gravitational field generated by a synthetic density model should closely agree with the Earth’s gravitational field.

Knowledge of permafrost’s ice and unfrozen water content is critical for predicting the permafrost behavior during ice–water–ice transition. This is especially relevant when ice and permafrost are melting in many regions under the influence of global warming. It is well-known that only a part of the formation's pore water turns into ice at 0 ^oC. After the further lowering the temperature, the water phase transition continues, but at gradually decreasing rates. Thus, the porous space is filled with ice and unfrozen water. The laboratory data show that frozen formations' mechanical, thermal, and rheological properties strongly depend on the moisture content. Hence, porosity and temperature are essential parameters of permafrost. In this paper, it is shown how, and by combining research in three fields: (1) geophysical exploration, (2) numerical modeling, and (3) temperature logging, it is possible to estimate in-situ the porosity of permafrost. To demonstrate the procedure, five examples of numerical modeling (where all input parameters are specified) are given. This investigation is the first attempt to analyze the permafrost’s porosity in situ quantitively.

The diffusion coefficient (D) is a key parameter characterizing gas transport in coal seams. Usually, the D is calculated from the desorption curve of particle coal, which cannot accurately reflect the diffusion characteristics under the stress constraint conditions of in-situ coal seams. In this paper, based on Fick's law of counter diffusion, different metamorphic deformed coals of medium and high coal rank are taken as the research object. The change laws of D under different confining pressure, gas pressure, and temperature conditions are tested and analyzed, and the influencing mechanisms on the D are discussed. The results show that the D of different metamorphic deformed coals decreases exponentially with the increase of confining pressure, and increases exponentially with the increase of gas pressure and temperature. There is a limit diffusion coefficient. The influence of confining pressure on the D is essentially determined by the change of effective stress, and the D has a negative effect of effective stress similar to permeability. The effect of gas pressure on the D involves two mechanisms, namely, mechanical and adsorption effects, which are jointly restricted by effective stress and coal particle shrinkage and expansion deformation. The effect of temperature on the D is mainly achieved by changing the root mean square speed and average free path of gas molecules. Under the same temperature and pressure conditions, the D increases first and then decreases with the increase of deformation degree, and the D of fragmented coal is the largest. Under similar deformation conditions, the D of high-rank anthracite is greater than that of medium-rank fat coal. It is considered that porosity is the key factor affecting the change of D in different metamorphic deformed coals.

The unreliable prediction of the low frequency components from inverted acoustic impedance causes uncertainty in quantitative seismic interpretation. To address this issue, we first calculate various seismic and geological attributes that contain low frequency information, such as relative geological age, interval velocity, and integrated instantaneous amplitude. Then, we develop a method to predict the low frequency content of seismic data using these attributes, their high frequency components, and recurrent neural networks. Next, we test how to predict the low frequency components using stacking velocity obtained from velocity analysis. Using all the attributes and seismic data, we propose a supervised deep learning method to predict the low frequency components of the inverted acoustic impedance. The results obtained in both synthetic and real data cases show that the proposed method can improve the prediction accuracy of the low frequency components of the inverted acoustic impedance, with the best improvement in a real data example of 57.7% compared with the impedance predicted using well log interpolation.